Reactor Criticality Simulator: k-Effective & Neutron Multiplication

simulator advanced ~12 min
Loading simulation...
k-eff ≈ 1.000 — reactor is exactly critical with balanced control rod insertion

At 3.5% enrichment with moderator ratio 4 and 50% control rod insertion, the reactor achieves k-eff ≈ 1.000 — steady-state criticality with stable neutron population and constant power output.

Formula

k_eff = η · f · p · ε · P_NL · T_NL (six-factor formula)
ρ = (k_eff - 1) / k_eff (reactivity definition)
T = ℓ / (k_eff - 1 - β) for period (prompt jump approximation)

The Chain Reaction

Nuclear fission releases 2–3 neutrons per split uranium-235 nucleus. If exactly one of those neutrons causes another fission, the reaction sustains itself at constant power — the reactor is critical. This deceptively simple balance is the central problem of nuclear engineering: maintaining k-effective at precisely 1.000 while the fuel depletes, fission products accumulate, and temperatures fluctuate.

The Six-Factor Formula

The effective multiplication factor k_eff is the product of six probabilities tracing a neutron's life from birth to the next fission. Each factor — fast fission, resonance escape, thermal utilization, reproduction, and two non-leakage probabilities — depends on geometry, materials, temperature, and isotopic composition. Changing any one parameter ripples through the entire neutron economy.

Delayed Neutrons: The Gift of Control

If all fission neutrons were emitted instantly (prompt neutrons), a reactor with k=1.001 would double in power in milliseconds — impossibly fast to control mechanically. Nature provides a saving grace: about 0.65% of neutrons are delayed, emitted seconds to minutes later from fission product decay. This tiny fraction stretches the effective neutron generation time from microseconds to seconds, making control rod adjustment feasible.

Temperature Feedback and Safety

This simulation models k-effective as a function of enrichment, moderation, control rod absorption, and temperature feedback. Adjust parameters to explore how the reactor transitions between subcritical, critical, and supercritical states. Watch how the Doppler effect at high temperatures provides inherent negative feedback — the single most important passive safety mechanism in thermal reactor design.

FAQ

What does k-effective mean in nuclear engineering?

k-effective (k_eff) is the ratio of neutrons in one generation to the previous generation. At k=1.0 (critical), the chain reaction is self-sustaining at constant power. At k<1.0 (subcritical), the reaction dies out. At k>1.0 (supercritical), power increases. Normal reactor operation adjusts k between 0.999 and 1.001.

What is the delayed neutron fraction?

About 0.65% of fission neutrons are emitted seconds to minutes after fission from decay of fission products (delayed neutrons). This tiny fraction is what makes reactor control possible — without delayed neutrons, the reactor period would be milliseconds instead of seconds, making mechanical control rod movement far too slow.

What is the Doppler effect in reactors?

As fuel heats up, thermal vibration broadens the resonance absorption peaks of U-238 (Doppler broadening), capturing more neutrons and reducing k-effective. This provides inherent negative temperature feedback — if power increases, fuel heats up, reactivity drops, and power self-corrects. This is the most important passive safety mechanism in thermal reactors.

How do control rods work?

Control rods contain neutron-absorbing materials (boron, hafnium, silver-indium-cadmium) that capture neutrons and remove them from the chain reaction. Inserting rods reduces k-effective; withdrawing them increases it. Operators adjust rod position to maintain k=1.0 at the desired power level, compensating for fuel burnup, xenon poisoning, and temperature changes.

Sources

Embed

<iframe src="https://homo-deus.com/lab/nuclear-engineering/reactor-criticality/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub