Neutron Multiplication k-effective Calculator

Compute k-effective for any bare homogeneous reactor geometry using the one-group diffusion approximation with migration area.

🔗 Neutron Multiplication k-effective Calculator
k_eff (Effective Multiplication Factor)
Geometric Buckling B²
Non-Leakage Probability P_NL
Reactivity ρ (Δk/k)
Reactivity in pcm
Criticality Status
k_eff
Reactivity ρ
Reactivity in pcm
Criticality Status
Deviation from Critical

🔗 What is k-effective?

The effective neutron multiplication factor k_eff is the most fundamental quantity in nuclear reactor physics. It describes the ratio of neutrons produced by fission in one generation to the total number of neutrons lost in the previous generation, including both absorption within the core and leakage from the reactor boundary. When k_eff equals exactly 1.000, the reactor is critical and sustains a steady, self-propagating chain reaction. When k_eff exceeds 1, the reactor is supercritical and its power output rises. When k_eff falls below 1, the reactor is subcritical and the chain reaction is not self-sustaining.

The relationship between k_eff and the infinite multiplication factor k∞ is given by k_eff = k∞ x P_NL, where P_NL is the non-leakage probability. k∞ describes how many neutrons each generation would produce in a hypothetical reactor of infinite size (no boundary, no leakage), while k_eff corrects for the reality that neutrons near the reactor surface can escape. In the one-group diffusion approximation, P_NL = 1/(1 + M²B²), where M² is the migration area of the fuel-moderator lattice and B² is the geometric buckling of the reactor. The full formula becomes k_eff = k∞ / (1 + M²B²).

Geometric buckling B² depends only on the shape and size of the reactor, not on the fuel or moderator composition. For a sphere of radius R, B² = (π/R)². For a finite cylinder of radius R and height H, B² = (2.405/R)² + (π/H)². For a rectangular box with sides a, b, c, B² = (π/a)² + (π/b)² + (π/c)². The migration area M² is a material property of the moderator-fuel lattice, defined as M² = L² + τ, where L² is the thermal diffusion area and τ is the Fermi age. For light water reactor lattices, M² is typically 50 to 70 cm², while heavy water systems have M² of 5,000 to 7,000 cm².

Engineers use k_eff calculations throughout the reactor design and operation lifecycle. During the design phase, lattice calculations determine the minimum critical size for a given fuel composition. During operation, k_eff monitoring (via neutron flux detectors) confirms the reactor remains at its target power level. During refueling and maintenance, verified deep subcriticality (k_eff well below 1) provides the safety basis for working near the core. The reactivity ρ = (k_eff - 1)/k_eff, expressed in units of pcm (10^-5), quantifies how far the reactor is from the critical point and drives all reactor control system responses.

📐 Formula

k_eff = k∞ ÷ (1 + M² × B²)
k_eff = effective multiplication factor (dimensionless)
k∞ = infinite multiplication factor (from the four-factor formula: k∞ = ηεpf)
= migration area of the fuel-moderator lattice (cm²) = L² + τ
= geometric buckling (cm−2) - depends on geometry only
Geometric Buckling by Shape:
Sphere (radius R): B² = (π/R)²
Finite Cylinder (radius R, height H): B² = (2.405/R)² + (π/H)²
Box (sides a, b, c): B² = (π/a)² + (π/b)² + (π/c)²
Non-Leakage Probability: P_NL = 1 / (1 + M²B²)
Reactivity: ρ = (k_eff − 1) / k_eff
ρ in pcm = ρ × 105
Critical Radius (sphere): R_crit = π / √((k∞ − 1) / M²)

📖 How to Use This Calculator

Steps

1
Select calculation mode - use "Calculate k_eff" to find the effective multiplication factor from reactor geometry and material properties, or "Reactivity from k_eff" to convert a known k_eff directly to reactivity in percent and pcm.
2
Enter reactor parameters - input k∞ (from the four-factor formula or reactor code), migration area M² in cm² (50-70 for LWR, 5,000-7,000 for heavy water), and select the reactor geometry (sphere, cylinder, or box) then enter the dimensions in centimeters.
3
Click Calculate - the calculator returns k_eff, geometric buckling B² in scientific notation, non-leakage probability P_NL as a percentage, and reactivity ρ in both percent and pcm, plus a criticality status label.
4
Use reactivity mode for shutdown analysis - switch to "Reactivity from k_eff", enter the estimated or measured k_eff value, and read off the reactivity, pcm deviation from criticality, and criticality status directly.

💡 Example Calculations

Example 1 - Large LWR Sphere (k∞=1.30, M²=60 cm², R=400 cm)

Spherical LWR core: k∞ = 1.30, M² = 60 cm², radius R = 400 cm

1
B² = (π/400)² = (0.007854)² = 6.168 × 10-5 cm-2
2
M²B² = 60 × 6.168 × 10-5 = 0.003701
3
k_eff = 1.30 / (1 + 0.003701) = 1.30 / 1.003701 = 1.2952
4
P_NL = 1.2952 / 1.30 = 99.63%; ρ = (1.2952-1)/1.2952 = 22.79% = 22,792 pcm
k_eff = 1.2952 | P_NL = 99.63% | Reactivity = 22,792 pcm
Try this example →

Example 2 - Finite Cylinder (k∞=1.28, M²=60 cm², R=150 cm, H=300 cm)

Cylindrical reactor core: k∞ = 1.28, M² = 60 cm², R = 150 cm, H = 300 cm

1
B² = (2.405/150)² + (π/300)² = (0.01603)² + (0.01047)² = 2.570 × 10-4 + 1.097 × 10-4 = 3.667 × 10-4 cm-2
2
M²B² = 60 × 3.667 × 10-4 = 0.02200
3
k_eff = 1.28 / (1 + 0.02200) = 1.28 / 1.02200 = 1.2524
4
P_NL = 1.2524 / 1.28 = 97.84%; ρ = (1.2524-1)/1.2524 = 20.15% = 20,147 pcm
k_eff = 1.2524 | P_NL = 97.84% | Reactivity = 20,147 pcm
Try this example →

Example 3 - Shutdown Reactivity Analysis (k_eff = 0.95)

Control rods fully inserted, reactor shutdown: k_eff = 0.95

1
ρ = (0.95 - 1) / 0.95 = -0.05 / 0.95 = -0.05263
2
Reactivity in % = -5.263%; in pcm = -5,263 pcm
3
Deviation from criticality: |ρ| = 5,263 pcm subcritical (well within shutdown margin requirements)
k_eff = 0.950 | Reactivity = -5,263 pcm | Status: Subcritical
Try this example →

❓ Frequently Asked Questions

What is the effective multiplication factor k_eff in a nuclear reactor?+
The effective multiplication factor k_eff is the ratio of neutrons produced by fission in one generation to the neutrons lost (by absorption or leakage) in the previous generation. When k_eff = 1, the reactor is critical and maintains a steady chain reaction. When k_eff > 1, the reactor is supercritical and power rises. When k_eff < 1, the reactor is subcritical and the chain reaction dies out. k_eff differs from k∞ by including neutron leakage from the finite reactor boundary.
What is the one-group diffusion formula for k_eff?+
In the one-group diffusion approximation, k_eff = k∞ / (1 + M²B²), where k∞ is the infinite multiplication factor, M² is the migration area (cm²), and B² is the geometric buckling (cm⁻²). The factor 1/(1 + M²B²) equals the non-leakage probability P_NL. This formula applies to bare, homogeneous thermal reactors and is accurate when M²B² is small compared to 1.
What is geometric buckling and how is it calculated for different shapes?+
Geometric buckling B² measures how rapidly the neutron flux varies across the reactor and depends only on geometry. For a sphere of radius R: B² = (π/R)². For a finite cylinder of radius R and height H: B² = (2.405/R)² + (π/H)². For a rectangular box with sides a, b, c: B² = (π/a)² + (π/b)² + (π/c)². The unit is cm⁻² when dimensions are in cm. Larger reactors have smaller buckling and lower leakage.
What is the migration area M² and what values does it take for common reactor types?+
Migration area M² = L² + τ, where L² is the diffusion area and τ is the Fermi age. M² quantifies the mean squared distance a neutron travels from birth to absorption. For light water with UO2 fuel, M² is roughly 50-70 cm². For heavy water with UO2 fuel, M² is 5,000-7,000 cm². For graphite-moderated reactors, M² is 300-500 cm². Larger M² means neutrons travel farther, so leakage is more significant for a given core size.
What is the difference between k∞ and k_eff?+
k∞ is the multiplication factor in a hypothetical infinite reactor where no neutrons leak. k_eff accounts for the actual finite size: k_eff = k∞ × P_NL, where P_NL is the non-leakage probability. For a critical reactor, k_eff = 1.000. Large power reactors have P_NL of 0.95 to 0.99, so k_eff is only slightly below k∞. Smaller reactors lose proportionally more neutrons through leakage.
What is reactivity and how is it measured in pcm?+
Reactivity ρ = (k_eff - 1) / k_eff. It is zero at criticality, positive when supercritical, and negative when subcritical. The unit pcm (percent-milli) equals 10^-5. A typical power reactor at beginning-of-life has excess reactivity of 15,000-25,000 pcm, controlled by soluble boron, control rods, and burnable poisons. One dollar of reactivity equals the delayed neutron fraction β, approximately 650 pcm for U-235 fuel.
What is the optimum height-to-diameter ratio for a cylindrical reactor?+
The optimum height-to-diameter ratio that minimizes critical volume for a finite cylinder is H/D ≈ 0.924. At this ratio, the radial and axial components of buckling are equal, meaning neither direction contributes disproportionately to leakage. Real reactor designs deviate from this optimum for engineering reasons (containment dimensions, coolant flow), but it serves as the baseline for bare-reactor lattice calculations.
How does increasing reactor size affect k_eff?+
As reactor dimensions increase, geometric buckling B² decreases (proportional to 1/R² for a sphere), so the term M²B² decreases and k_eff approaches k∞. In the limit of an infinite reactor, B² approaches zero and k_eff = k∞. For a given fuel composition, there is a minimum critical size at which k_eff = 1. Reactors smaller than this critical size cannot sustain a chain reaction regardless of fuel enrichment.
What is the critical radius of a spherical reactor?+
For a bare spherical reactor, criticality requires B² = (k∞ - 1)/M². Since B² = (π/R)² for a sphere, the critical radius is R_crit = π / √((k∞ - 1)/M²). For a typical LWR with k∞ = 1.30 and M² = 60 cm², R_crit = π × √(60/0.30) = π × √200 ≈ 44.4 cm. This is the extrapolated radius; the physical radius is slightly smaller by the extrapolation distance (about 0.71 × transport mean free path).
What values of k_eff are used during reactor startup and shutdown?+
During normal operation, a critical reactor maintains k_eff = 1.000 exactly. During startup, k_eff is brought slightly above 1 to allow controlled power increase. For safe shutdown, control rods are inserted to achieve a deeply subcritical k_eff (often 0.95 or below), providing a shutdown margin of at least 5,000 pcm. Regulatory requirements typically demand a shutdown margin of at least 1,000-2,000 pcm with the most reactive control rod stuck out.
How do fission product poisons like xenon-135 affect k_eff?+
Xenon-135 is the most powerful neutron absorber in reactor operations, with a thermal absorption cross-section of 2.65 million barns. At equilibrium in a high-flux reactor, Xe-135 reduces k_eff by several thousand pcm, typically 2,000-3,000 pcm for LWRs. After shutdown or power reduction, the Xe-135 concentration peaks within 6-12 hours due to I-135 decay, causing a temporary additional negative reactivity that can prevent immediate restart. This xenon poisoning must be managed carefully during reactor operations.
What is the two-group model and how does it differ from the one-group model?+
The two-group model separates neutrons into fast and thermal energy groups. It gives k_eff = k∞ × P_FNL × P_TNL, where P_FNL = exp(-B²τ) is the fast non-leakage probability and P_TNL = 1/(1 + L²B²) is the thermal non-leakage probability. The one-group model combines these into 1/(1 + M²B²) = 1/(1 + (L² + τ)B²). The two-group model is more accurate for small reactors or heavy-water systems with large M².