Why a reactor is controllable — and why it stops itself.
A nuclear reactor lives on a knife-edge: the neutron population can double in a fraction of a second, yet operators steer it by hand. Two things make that possible — the delayed neutrons that slow the response to human time scales, and the temperature feedback that makes the core fight its own power rise. This program models both, from six-group point kinetics up through Doppler self-limiting, the xenon poison that can lock a reactor out of restart, load-follow control, and emergency shutdown.

Delayed neutrons buy the time.
If a reactor ran on prompt neutrons alone, a small reactivity insertion would blow the power up in milliseconds — uncontrollable. But a fraction of a percent of neutrons emerge seconds to minutes later, from decaying fission products, and those delayed neutrons set the pace: insert a little reactivity and the power climbs on a period of tens of seconds, slow enough to steer. The model captures the prompt jump (a +0.5 dollar step leaps power to 2.07×, matching theory), the inhour relation between reactivity and period, and the cliff at one dollar where the delayed neutrons stop mattering and the reactor goes prompt-critical.


The core that fights its own power rise.
Reactivity isn't fixed — it depends on temperature. As fuel heats, the Doppler broadening of absorption resonances captures more neutrons, pushing reactivity down. That negative feedback is what makes a power excursion self-limit instead of running away. The model shows it starkly: a +0.6 dollar insertion that explodes to nearly thirty-thousand times power with no feedback instead peaks at 2.4× and settles, once Doppler is in the loop. It's the same physics that makes a well-designed reactor inherently stable.



Every number is one you can re-run.
The sign-off notebook re-derives each requirement from the same point-kinetics model the program builds.
| Result | Requirement | |
|---|---|---|
| Prompt jump (+0.5 $ step) | 2.07× | theory 2.0 |
| Doppler-limited excursion | 2.4× peak | bounded, not runaway |
| Xenon iodine-pit | −4459 pcm @ 8 h | post-shutdown |
| Load-follow tracking error | 2.4 % | 100→60→100 % |
| Requirements verified | 6 / 6 | PASS |
0-D neutronics, a teaching & controls model.
The reactor is 0-D point kinetics — the neutron population as a single number, with no spatial flux shape, so there are no spatial xenon oscillation modes (a noted follow-up). Thermal feedback is lumped, the xenon is single-node, and the cross-sections are illustrative. This is a teaching and control-design model, not a licensing-grade neutron-transport code. It is exactly the fidelity that builds intuition and supports control logic — the prompt jump, the inhour relation, the safety feedbacks, the xenon transient, load-follow and SCRAM — on reproducible numbers, before a spatial diffusion or Monte-Carlo transport solver.
Model your own reactor kinetics.
Book a walkthrough and we'll set up your delayed-neutron data, feedback coefficients and control scheme, and run the kinetics, feedback and xenon studies live.
