Multiplication Factor Calculator

When k = 1, the reactor is exactly critical: each generation of neutrons produces exactly the same number of neutrons as the previous generation, maintaining a steady-state fission chain reaction. This is the condition for stable power operation. If k > 1, the reactor is supercritical and the neutron population (and power) increases exponentially. If k < 1, the reactor is subcritical and the neutron population decreases, eventually shutting down. Achieving k = 1 precisely is the goal of reactor control systems.

k∞ is the multiplication factor for an infinite reactor (no neutron leakage). It depends only on the material composition and is given by the four‑factor formula: k∞ = η ε p f. k_eff is the multiplication factor for a finite reactor, which accounts for neutron leakage: k_eff = k∞ × (probability of non‑leakage). For a critical reactor, k_eff = 1, while k∞ is always greater than 1 because some neutrons leak out. The difference between k∞ and k_eff represents the fraction of neutrons lost to leakage. In a large power reactor, the leakage is typically 3-5%, so k∞ is about 1.03-1.05.

Control rods contain strong neutron absorbers (e.g., B₄C, Ag-In-Cd) that increase the absorption term in the denominator of k = Production / (Absorption + Leakage). Inserting control rods increases absorption, reducing k. Reactivity (ρ) is defined as ρ = (k - 1) / k and is often expressed in units of 'dollars' or 'pcm' (percent mille, 1 pcm = 10⁻⁵). When control rods are fully inserted, k drops below 1, shutting down the reactor. The amount of reactivity worth of a control rod is the change in k when it is moved from fully withdrawn to fully inserted, typically a few thousand pcm.

As the reactor operates, fissile isotopes (U‑235, Pu‑239) are depleted, reducing the production term in k = Production / (Absorption + Leakage). Additionally, fission products (like Xe‑135) and transuranics (like Pu‑240) build up, increasing the absorption term. This gradually reduces k. To maintain k = 1, the reactor uses: (1) burnable absorbers (e.g., gadolinia) that initially suppress excess reactivity and burn out over time, (2) control rods that are gradually withdrawn, and (3) soluble boron (in PWRs) that is diluted to maintain criticality. At the end of the cycle, the reactor is shut down for refueling.

Delayed neutrons are emitted by fission products (precursors) with half-lives ranging from milliseconds to seconds. They account for only about 0.65% of all neutrons in a thermal reactor (for U‑235 fission). While they don't change the steady-state value of k, they determine the reactor's time response. The prompt multiplication factor (k_p) is based only on prompt neutrons; if k_p > 1, the reactor is prompt critical, leading to a very rapid power increase. Delayed neutrons allow the reactor to be controlled by adjusting k within a narrow range (from k = 1 to k = 1 + β, where β is the delayed neutron fraction, about 0.0065 for U‑235). This makes control possible by mechanical means.

The neutron diffusion equation is a form of the neutron balance: D∇²φ - Σ_a φ + νΣ_f φ = (1/v) ∂φ/∂t. In steady state, the criticality condition is that the neutron source (νΣ_f φ) equals the sum of absorption and leakage. The multiplication factor is defined as k = (νΣ_f φ) / (Σ_a φ - D∇²φ) for a finite system. In the one‑group diffusion approximation, k_eff = (νΣ_f) / (Σ_a + D B²), where B² is the buckling (a measure of leakage). This shows that k depends on the material cross‑sections and the core geometry via the buckling.

In most thermal reactors, k decreases with increasing temperature due to two main effects: (1) the Doppler broadening of U‑238 resonances, which increases resonance absorption (reducing the resonance escape probability, p), and (2) the moderator temperature coefficient, which in LWRs is negative because as water heats up, its density decreases, reducing moderation and increasing leakage. These effects combine to give a negative temperature coefficient of reactivity: when the reactor heats up, k decreases, reducing the fission rate and providing inherent stability. This is a key safety feature: if the reactor overheats, it naturally shuts down. The magnitude of this feedback is typically a few pcm per °C.

The six‑factor formula is an extension of the four‑factor formula that accounts for non‑leakage probabilities of fast and thermal neutrons. It is written as k_eff = η ε p f P_FNL P_TNL, where P_FNL is the fast non‑leakage probability (fraction of fast neutrons that do not leak), and P_TNL is the thermal non‑leakage probability. These factors depend on the core dimensions and the diffusion lengths. The formula helps in understanding how k_eff is influenced by the material properties and the geometry. It is used in nuclear design codes to calculate the multiplication factor for different core configurations.

Burnable absorbers (like Gd₂O₃ or B₄C) are added to the fuel to absorb excess reactivity at the beginning of the cycle (BOC). They increase the absorption term in the denominator, reducing k initially. As the reactor operates, these absorbers burn out (via neutron capture), gradually reducing their absorption, which allows k to increase and compensate for the depletion of fissile material. The net effect is a relatively flat reactivity curve over the cycle. At the end of the cycle (EOC), most burnable absorbers are gone, and k is maintained solely by the remaining fissile material and control rod movements.

Yes, an external neutron source (like a startup source) adds neutrons to the system, increasing the neutron population. However, the multiplication factor k (defined as the ratio of neutrons in one generation to the next) is a property of the reactor itself, independent of the source. If k < 1, the neutron population will eventually die out once the source is removed; the source only sustains a steady flux by continuously adding neutrons. The effective multiplication factor remains k = 0.99 (for example), but the neutron flux is maintained by the source. This is used in subcritical assemblies for research or for neutron multiplication measurements.

The reactor period is the time required for the neutron population (or power) to change by a factor of e (2.718). For a small reactivity insertion (ρ small), the reactor period is approximately given by the inhour equation, which relates the reactivity to the delayed neutron parameters and the prompt neutron lifetime. In a supercritical reactor (k > 1), the period is determined by the excess reactivity (k - 1). For large reactivities, the period is short, leading to rapid power excursions. The multiplication factor is the steady-state ratio, while the period governs the transient behavior. They are linked by the reactor kinetics equations.

In fast reactors, the neutron spectrum is hard (energies > 100 keV), and the fission cross‑sections are lower than in thermal reactors. The multiplication factor in a fast reactor depends heavily on the number of neutrons emitted per fission (ν), which is slightly higher for Pu‑239 than for U‑235, and on the capture‑to‑fission ratio (α). Fast reactors have a higher η (neutrons produced per absorption) than thermal reactors, but they also have higher leakage due to the long mean free paths. The multiplication factor in fast reactors is often expressed in terms of the breeding ratio (BR), because fast reactors can breed more fissile material. The criticality condition is achieved by a compact core with high fissile concentration, making k_eff sensitive to the core geometry and the presence of structural materials.

The shutdown margin is the reactivity worth (in pcm or dollars) of all control rods when they are fully inserted, which ensures that the reactor can be safely shut down from any operating condition. It is typically defined as the negative reactivity required to make the reactor subcritical (k_eff < 1) even with the most reactive control rod stuck out. In terms of k, the shutdown margin is the difference between the k_eff with all rods in and the k_eff with the stuck rod. A typical shutdown margin for a PWR is about 1,000-2,000 pcm (1-2% Δk/k), ensuring that if the reactor needs to be scrammed, it will remain subcritical with a high degree of confidence.

Xe‑135 and Sm‑149 are strong neutron absorbers (with thermal absorption cross‑sections of ~2.6×10⁶ barns and ~41,000 barns, respectively). As they build up from fission product decay, they increase the absorption term in the multiplication factor equation, reducing k. This is called reactivity poisoning. Xe‑135 is particularly significant because of its large cross‑section and its decay chain (I‑135 → Xe‑135 → Cs‑135). After reactor shutdown, the Xe‑135 concentration can increase (due to decay of I‑135), leading to 'reactor dead time' where the reactor cannot be restarted until the Xe‑135 decays (about 10-20 hours). Sm‑149 is less transient but builds up over longer times. Reactor designers must include these poisons in the k calculation and provide sufficient excess reactivity to overcome them.

The multiplication factor cannot be measured directly; it is inferred from the neutron flux response to a known reactivity change. During startup, the reactor is brought to critical by slowly withdrawing control rods while monitoring the neutron count rate. The inverse count rate (1/CR) is plotted against control rod position. When the inverse count rate extrapolates to zero, the reactor is critical (k = 1). This method is called the 'approach to criticality' test. Alternatively, the reactivity can be measured using the reactor period method or using a reactivity meter that compares the neutron flux to the delayed neutron precursors. Once critical, the control rod position is recorded, and the multiplication factor is determined by the known reactivity worth of the rods.