Criticality K Effective Interactive Calculator

← Back to Engineering Library

Determining if a fissile system is subcritical, critical, or supercritical is a fundamental part of nuclear engineering. Mistakes here can lead to anything from non-startup to a criticality accident, so these calculations are not taken lightly. This K-Effective (keff) calculator helps you compute the effective neutron multiplication factor using actual values for the six-factor formula, cross-sections, buckling, reactivity, and period. You'll see keff calculations used in reactor design, spent fuel storage safety checks, and fuel loading steps where maintaining enough distance from criticality allows work to proceed without remote handling. On this page you get the main equations, a PWR startup worked example, some basics on neutron kinetics and temperature effects, and a Q&A for edge cases and field situations.

What is k-effective (keff)?

K-effective indicates if a chain reaction will fizzle, sustain, or grow. At keff = 1.0, the neutron chain reaction holds steady. If it drops below 1.0, the chain reaction dies off. Above 1.0, neutron population increases each generation.

Simple Explanation

Keff is just a measure of how many neutrons from one generation will survive to produce fissions in the next. 1.0 means the numbers balance out—a steady reaction. In any real reactor, keff is nudged just above 1.0 for power generation and kept well below 1.0 during storage or shutdown by control rods, boron, and physical layout—whatever gets the job done with as little excess as practical.

📐 Browse all 1000+ Interactive Calculators

Visual Diagram: Neutron Multiplication Cycle

Criticality K Effective Interactive Calculator Technical Diagram

K-Effective Interactive Calculator

How to Use This Calculator

Engineering calculation notice

This calculator is intended for education, concept evaluation, and preliminary design. Results are based on the equations and assumptions described on this page, but cannot account for every real-world load case, tolerance, material property, environmental condition, installation detail, safety factor, code, or regulatory requirement. Verify all inputs, assumptions, units, and results independently before selecting components or using the result in a real application. Safety-critical, structural, medical, lifting, transportation, or regulated applications must be reviewed by a qualified engineer.

Found a calculation error? Message us

  1. Pick the kind of calculation you want—six-factor k-effective, η, critical mass, buckling, reactivity, or period.
  2. Enter the numbers actually relevant for your scenario—things like eta, fast fission factor, cross sections, or reactivity in dollars.
  3. The "Try Example" button is there if you want to see typical numbers or check your method against a worked example.
  4. Press Calculate. Your answer comes out with the inputs you've set.

Criticality K-Effective Interactive Visualizer

Visualize neutron multiplication through reactor generations to understand how the six-factor formula determines whether your nuclear system goes critical. Watch neutron populations grow, decay, or stabilize as you adjust multiplication factors.

η (Eta) 2.07
ε (Fast Fission) 1.03
p (Resonance Escape) 0.875
f (Thermal Utilization) 0.72
Lf (Fast Non-Leakage) 0.95
Lth (Thermal Non-Leakage) 0.98

K-EFFECTIVE

1.249

REACTIVITY

+0.199

STATUS

SUPER

FIRGELLI Automations — Interactive Engineering Calculators

Core Equations for Criticality Analysis

Here’s how you combine the main parameters to get k-effective with the six-factor formula.

Six-Factor Formula for k-effective

keff = η × ε × p × f × Lf × Lth

Where:

  • η (eta) = Reproduction factor (neutrons produced per neutron absorbed in fuel)
  • ε (epsilon) = Fast fission factor (≥1, typically 1.02-1.10)
  • p = Resonance escape probability (fraction avoiding capture during slowing down, 0-1)
  • f = Thermal utilization factor (fraction absorbed in fuel vs. all materials, 0-1)
  • Lf = Fast non-leakage probability (fraction not leaking during slowing down, 0-1)
  • Lth = Thermal non-leakage probability (fraction not leaking as thermal neutrons, 0-1)

Here’s how you get the reproduction factor η from actual cross-section numbers.

Reproduction Factor (η)

η = ν × (σf / σa)

Where:

  • ν = Average number of neutrons released per fission event
  • σf = Microscopic fission cross section (barns)
  • σa = Microscopic absorption cross section (barns)

If you know k-effective, here’s the quickest way to estimate reactivity.

Reactivity

ρ = (keff - 1) / keff

Where:

  • ρ (rho) = Reactivity (dimensionless, or expressed in pcm, dollars, cents)
  • keff = Effective multiplication factor

Note: 1 pcm = 10-5 Δk/k; 1 dollar = βeff (typically 0.0065 for thermal reactors)

If you need material buckling for a geometry problem, use this expression.

Material Buckling

B² = (k - keff) / (M² × keff)

Where:

  • = Material buckling (cm-2)
  • k = Infinite multiplication factor
  • = Migration area (cm²)

This gets the reactor period if you know your reactivity and generation time.

Reactor Period

T = Λ / (ρ - β)

Where:

  • T = Reactor period (seconds)
  • Λ = Neutron generation time (seconds)
  • ρ = Reactivity (absolute value)
  • β = Delayed neutron fraction

Note: Valid for ρ greater than β (subcritical or slightly supercritical)

Simple Example

With the six-factor model, try these values:

  • η = 2.07, ε = 1.03, p = 0.875, f = 0.720
  • Lf = 0.950, Lth = 0.980

k = 2.07 × 1.03 × 0.875 × 0.720 = 1.3432

keff = 1.3432 × 0.950 × 0.980 = 1.2494 — so this system is supercritical with these assumptions.

Theory & Engineering Applications of K-Effective

Fundamental Neutron Physics

Keff represents how the neutron population changes generation to generation in real systems. It pulls together the fundamental physics of neutron creation, moderation, absorption, and system geometry all in one number. At keff = 1, power is steady because each generation exactly replaces itself; below 1, power decays; above 1, it grows exponentially.

The six-factor formula breaks keff into key parameters. k (the product η × ε × p × f) assumes an infinite system—so no leakage losses. The non-leakage terms Lf and Lth pull in the effect of system size and shape. In practice, k is set by fuel and moderator chemistry, while geometric leakage is a factor of your core shape and surface-to-volume ratio. For example, a fresh fuel assembly can be calculated at k ≈ 1.35, but installed in a small core with lots of surface area, actual keff could be nearer 1.05, sometimes less.

The Reproduction Factor and Fuel Selection

The realistic limit for criticality begins with η. For U-235 at thermal energies, the typical numbers (ν ≈ 2.43, σfa ≈ 0.85) give η ≈ 2.07, so roughly two neutrons survive per absorption. Natural uranium only has 0.711% U-235, so η for bulk natural uranium is much lower—too low for criticality unless you're using heavy water or pure graphite as a moderator. U-238 absorbs neutrons but doesn't fission at thermal energies, so it just 'steals' them, and plutonium breeding from those absorptions plays a role in later cycles or in fast reactors.

η isn't a fixed value—it shifts with neutron energy. For U-235, η is highest a bit above thermal, then drops; for Pu-239, η stays higher over a wider energy band. Fast breeder and MOX fuel designs use these differences. Breeding of Pu-239 from U-238 (via neutron capture followed by beta decay through Np-239) is how you shift a thermal spectrum core toward a plutonium-producing one if left long enough.

Criticality Safety and Subcritical Margin

Maintaining a comfortable subcritical margin is a non-negotiable in nuclear criticality safety. For most storage and handling, keff ≤ 0.95 is required, and margin tightens for certain accident conditions. That 5% margin handles everything from uncertainties in cross-sections (often up to 1%), computational model errors, assembly tolerances, and real-world chaos. Monte Carlo methods (like MCNP or KENO) can nail down keff to a few tenths of a percent on known configurations, but the cumulative error sources add up. The sum is a built-in safety buffer—the difference between just barely subcritical in calculation and subcritical even if a few things line up against you in the field.

For spent fuel pools: A burned PWR fuel assembly might still have k over 1.2. Stack these up close in a pool, and you could create a critical system. The fix? Spreading the fuel out, and putting in neutron absorbers like Boral panels to lower keff. Changing rack spacing from 10 to 25 cm with absorbers can make the difference between a safe and totally unsafe facility—well worth a geometry check at the early design stage, before you spend on concrete and steel.

Reactor Period and Kinetics

The link between reactivity and reactor period controls how fast a reactor responds to changes. For positive reactivity below prompt critical, periods are set by delayed neutrons—so even with a few dollars of excess, the system responds in fractions of a second. In a typical thermal reactor (Λ ≈ 0.0005 s, ρ ≈ 0.003), you get a period around 0.14 s—not instant, but quick. The reality is tempered further by delayed neutron groups with different decay times, so power changes aren't as fast as the math alone would suggest.

If you cross prompt critical (ρ > β), only prompt neutrons matter, and the response becomes extremely fast. In a fast reactor (Λ ~ 10-7 s), adding 0.01 reactivity gives a doubling period on the order of microseconds. No mechanical or electronic control will keep up; at that point, it's all about passive physics and negative coefficients, not clever feedbacks or operators.

Worked Example: Complete PWR Startup Analysis

Scenario: Bringing a PWR from shutdown to critical, a startup engineer estimates required critical boron concentration.

Given Data:

  • Beginning of cycle, fresh fuel, η = 2.09
  • Fast fission factor ε = 1.028 (from U-238)
  • Resonance escape p = 0.883 (based on fuel-moderator ratio)
  • Thermal utilization f = 0.745 (measured before boron’s added)
  • Non-leakage Lf = 0.972, Lth = 0.990 (large core)
  • Boron worth: -8.5 pcm per ppm

Step 1: k without boron

k = 2.09 × 1.028 × 0.883 × 0.745 = 1.4154

Step 2: keff without boron

keff = 1.4154 × 0.972 × 0.990 = 1.3619

Step 3: Excess reactivity

ρexcess = (1.3619 - 1) / 1.3619 = 0.2657

Which is 0.2657 × 100,000 = 26,570 pcm

Step 4: Boron concentration needed

Boron = 26,570 pcm ÷ 8.5 pcm/ppm = 3,126 ppm

Step 5: Safety check

With all rods out and 3,126 ppm boron, keff should be around 1.000. As burnup progresses, boron is reduced to keep things at critical. Actual startup always cross-checks these predictions as criticality is approached with careful monitoring—never blindly by calculation.

Practical note: Expect the prediction to be off ±50 ppm for real startup due to actual moderator temperature, xenon buildup, and other secondary effects. This method still gives a reliable ballpark for planning and safety margins.

To try other calculations like these—or to address anything from fluid flow to bolted joint stress—see the full calculator library here.

Temperature and Void Coefficients

The effect of temperature on keff is a main driver for any reactor safety case. For light water reactors, hotter moderator means lower density, so fewer hydrogen atoms to slow neutrons, shifting the spectrum higher and dropping η for U-235. At the same time, thermal expansion opens up the fuel geometry, increasing leakage. A typical PWR design sees a moderator coefficient of -30 pcm/°C, so a 10°C swing will knock out ~300 pcm—enough to help keep power under control just by the feedback from heating itself.

The void coefficient is even more important in boiling water reactors. If you get too much steam in the core (loss of moderation), the reactivity drops—unless your fuel is heavily enriched or your neutron spectrum is already hard, in which case things can swing positive. The RBMK (Chernobyl) disaster came from a positive void coefficient—mainly an effect of the core’s enrichment and the graphite moderator, not just the operational mistake. In most commercial BWRs, the void coefficient is negative by design to avoid this scenario.

Practical Applications

Scenario: Research Reactor Startup Procedure

After a weekend shutdown, Dr. Chen checks her TRIGA reactor will go critical with the expected rod position. Her numbers show keff = 1.035 for the fully withdrawn state, but needs to be keff = 1.000 for startup. Rod worth is about 900 pcm each (four rods). She does the math: that's 0.0337 (3,370 pcm) excess, so expects to reach criticality with rods about 94% out. She withdraws slowly, checks neutron count rates, and expects any deviation to be within ±2%, flagging unexpected configuration shifts if it’s off more than that.

Scenario: Spent Fuel Storage Criticality Analysis

Marcus is reviewing a redesign for spent fuel storage racks. Each old PWR assembly has k ≈ 1.18 (even after a lot of burnup). MCNP says keff = 0.983 fully loaded—too close to the 0.95 limit. Using the calculator’s buckling mode, he sees he needs to reduce keff by 0.033, options being either adding more spacing (takes up construction cost and floor space) or adding more absorber. He checks that a 0.02" Boral panel between units drops keff far enough (by up to -4,000 pcm) for the needed margin. This directs his next round of layout before plowing time into massive MCNP reruns.

Scenario: Fast Reactor Core Design Optimization

Elena is tuning a sodium fast reactor core (Pu-239/U-238, 18% fissile). No moderator means some terms drop out: p ≈ 1, but ε is much higher (1.15–1.25 is normal for U-238-rich fuel in a fast spectrum). She works with η = 2.45, ε = 1.20, "fast utilization" f = 0.83, and non-leakage Lf = 0.88. The tool predicts keff = 1.045—enough to handle fuel depletion and still stay safely subcritical during refueling. For transients, Elena looks at βeff = 0.0035, Λ = 0.5 μs, and finds even 50 cents of extra reactivity gives a very fast period (0.35 s). This tells her, in practice, she can't rely only on rods or external equipment for control—she needs enough negative feedback and maybe passive systems for any transient the core could see.

Frequently Asked Questions

What is the difference between k-effective and k-infinity? +

Why is k-eff slightly greater than 1.0 in operating reactors? +

How does fuel enrichment affect k-effective? +

What is prompt criticality and why is it dangerous? +

How do neutron reflectors affect k-effective? +

What uncertainties affect k-effective calculations? +

Free Engineering Calculators

Explore our complete library of free engineering and physics calculators.

Browse All Calculators →

About the Author

Robbie Dickson — Chief Engineer & Founder, FIRGELLI Automations

Robbie Dickson brings over two decades of engineering expertise to FIRGELLI Automations. With a distinguished career at Rolls-Royce, BMW, and Ford, he has deep expertise in mechanical systems, actuator technology, and precision engineering.

Wikipedia · Full Bio

📹 Video Walkthrough — How to Use This Calculator

📹 Video Walkthrough — How to Use This Calculator

Criticality K Effective Interactive Calculator

Need to implement these calculations?

Explore the precision-engineered motion control solutions used by top engineers.

Share This Article
Tags