Chemistry / Nuclear & Radiochemistry Decay Kinetics, Fission, Fusion & Radiometrics 100% Free Open Access
Chapter 5 β€’ Theory & Derivations

Unit 5: Nuclear Fission Mechanics & Reactor Physics

Comprehensive physical and technological exposition of nuclear fission: the historical discovery by Hahn, Strassmann, Meitner, and Frisch; the Bohr-Wheeler liquid drop fission barrier; the asymmetric double-hump fragment mass yield curve; prompt versus delayed neutron kinetics; neutron moderation dynamics; the four-factor and six-factor criticality equations; and modern nuclear reactor systems.

Β§5.1 Historical Discovery & Thermodynamics of Nuclear Fission: Hahn, Meitner & Frisch

In December 1938, German radiochemists Otto Hahn and Fritz Strassmann at the Kaiser Wilhelm Institute in Berlin irradiated natural uranium with thermal neutrons. Expecting to identify transuranic elements ($Z > 92$), they performed fractional crystallization with barium carriers. Astonishingly, the radioactivity precipitated identically with barium ($Z = 56$), an element whose mass is barely half that of uranium.

The Meitner-Frisch Theoretical Interpretation (January 1939)

In exile in Sweden, Lise Meitner and her nephew Otto Frisch interpreted the result using George Gamow and Niels Bohr's Liquid Drop Model of the nucleus:

$$^{235}_{92}\text{U} + ^{1}_{0}n \longrightarrow [^{236}_{92}\text{U}^*] \longrightarrow ^{141}_{56}\text{Ba} + ^{92}_{36}\text{Kr} + 3 \, ^{1}_{0}n + Q$$

Frisch coined the term fission by analogy with binary fission in cellular biology.

Energetic Partition of Nuclear Fission

The total energy released per fission of uranium-235 is approximately $200\text{ MeV}$ ($\approx 3.204 \times 10^{-11}\text{ J}$). This immense thermodynamic yield arises from the drop in binding energy per nucleon: from $\approx 7.6\text{ MeV/nucleon}$ for $^{235}\text{U}$ to $\approx 8.5\text{ MeV/nucleon}$ for mid-mass fission fragments:

$$Q \approx 235 \times (8.5 - 7.6)\text{ MeV} \approx 200\text{ MeV}$$

``` ENERGY COMPONENT TYPICAL VALUE FRACTION ───────────────────────────────────────────────────────────────────────────── Kinetic Energy of Fission Fragments (Prompt) ~ 168 MeV 84.0% Kinetic Energy of Prompt Neutrons (~2.4 per fission) ~ 5 MeV 2.5% Prompt Gamma-Ray Photons ~ 7 MeV 3.5% Beta- Particles from Radioactive Fission Fragments ~ 8 MeV 4.0% Antineutrinos (Ξ½Μ„_e, escape reactor core completely) ~ 12 MeV 6.0% Delayed Gamma Photons from Fragment Decay Chains ~ 7 MeV 3.5% ───────────────────────────────────────────────────────────────────────────── TOTAL FISSION ENERGY RELEASE: ~ 207 MeV 100.0% RECOVERABLE THERMAL CORE ENERGY: ~ 195 MeV 94.0% ```

Over $80\%$ of the energy appears as kinetic energy of the two massive, highly charged fission fragments. They travel a mere $\sim 10\,\mu\text{m}$ in uranium metal before stopping, converting their kinetic energy into intense localized thermal heat via Coulomb ionization collisions.

Β§5.2 The Bohr-Wheeler Liquid Drop Fission Barrier & Fissility Parameter

In their 1939 paper, Niels Bohr and John Archibald Wheeler formulated the definitive classical theory of nuclear fission. They modeled the nucleus as an incompressible, charged liquid drop undergoing ellipsoidal quadrupole deformation.

Deformation Energy of a Liquid Drop

Consider a spherical nucleus of radius $R_0$ distorted into an axisymmetric prolate spheroid with semi-major axis $a = R_0(1 + \epsilon)$ and semi-minor axis $b = R_0(1 - \epsilon/2)$, where $\epsilon$ is the eccentricity deformation parameter:

  • Volume: Conserved ($\frac{4}{3}\pi a b^2 = \frac{4}{3}\pi R_0^3$).
  • Surface Area: Increases with deformation:
$$S(\epsilon) = S_0 \left( 1 + \frac{2}{5} \epsilon^2 + \dots \right)$$
  • Coulomb Self-Energy: Decreases with deformation (charge is pushed further apart):
$$E_C(\epsilon) = E_C^0 \left( 1 - \frac{1}{5} \epsilon^2 + \dots \right)$$

The net change in nuclear potential energy $\Delta E$ relative to the spherical ground state is:

$$\Delta E = \Delta E_S + \Delta E_C = E_S^0 \left( \frac{2}{5} \epsilon^2 \right) - E_C^0 \left( \frac{1}{5} \epsilon^2 \right) = \frac{1}{5} \epsilon^2 \left( 2 E_S^0 - E_C^0 \right)$$

``` Potential Energy Ξ”E(Ξ΅) β–² β”‚ /\ Fission Barrier E_f β”‚ / \ β”‚ / \ β”‚ / \____ Spontaneous Fission (Scission) β”‚ ___ / 0 ┼─( O )/────────────────► Deformation Ξ΅ Spherical Ground State ```

The Fissility Parameter ($x$)

The condition for spontaneous instability against infinitesimal deformation ($\Delta E < 0$) is:

$$2 E_S^0 - E_C^0 < 0 \implies \frac{E_C^0}{2 E_S^0} > 1$$

We define the dimensionless Fissility Parameter $x$:

$$x \equiv \frac{E_C^0}{2 E_S^0}$$

Substituting the SEMF expressions $E_S^0 = a_s A^{2/3}$ and $E_C^0 = a_c Z^2 / A^{1/3}$:

$$x = \frac{a_c \frac{Z^2}{A^{1/3}}}{2 a_s A^{2/3}} = \frac{a_c}{2 a_s} \left(\frac{Z^2}{A}\right)$$

Using $a_c \approx 0.711\text{ MeV}$ and $a_s \approx 17.80\text{ MeV}$:

$$\left(\frac{Z^2}{A}\right)_{\text{crit}} = \frac{2 a_s}{a_c} \approx \frac{2(17.80)}{0.711} \approx 50.1$$
$$x = \frac{Z^2 / A}{(Z^2 / A)_{\text{crit}}} \approx \frac{Z^2 / A}{50.1}$$
  • If $x \ge 1.0$ ($Z^2/A \ge 50$): The spherical nucleus is unstable to immediate spontaneous fission ($\tau \sim 10^{-22}\text{ s}$).
  • If $x < 1.0$: A finite Fission Barrier $E_f$ opposes deformation.

Fissile Versus Fertile Radionuclides

  • Fissile Radionuclides: Capable of undergoing fission with zero-energy thermal neutrons ($0.025\text{ eV}$). Examples: $^{233}\text{U}, ^{235}\text{U}, ^{239}\text{Pu}, ^{241}\text{Pu}$.
  • Fertile Radionuclides: Do not fission with thermal neutrons because neutron binding energy $S_n < E_f$. They require fast neutrons ($E_n > 1.0\text{ MeV}$) to fission, but can capture a thermal neutron to breed a fissile isotope. Examples: $^{238}\text{U} \to ^{239}\text{Pu}$, $^{232}\text{Th} \to ^{233}\text{U}$.

The difference is rooted in the pairing energy: When an odd-neutron nucleus ($^{235}_{92}\text{U}$) captures a neutron, it forms an even-even compound state ($[^{236}_{92}\text{U}^*]$), releasing an extra $+1.2\text{ MeV}$ of pairing binding energy ($S_n = 6.55\text{ MeV}$), which exceeds the $5.7\text{ MeV}$ fission barrier! Capturing a neutron on even-even $^{238}_{92}\text{U}$ forms an odd-neutron state ($[^{239}_{92}\text{U}^*]$) with zero pairing bonus ($S_n = 4.80\text{ MeV}$), falling short of the $6.2\text{ MeV}$ barrier.

Β§5.3 Energetics & Mass Distribution of Fission Fragments: Asymmetric Double-Hump Yields

When a heavy nucleus such as uranium-235 undergoes low-energy thermal fission, it does not split into two equal halves ($A_1 = A_2 \approx 117$). Instead, the fragment mass yield distribution exhibits a pronounced asymmetric double-humped morphology.

``` Fragment Yield (%) 10 β–² Peak 1 (Light) Peak 2 (Heavy) β”‚ A β‰ˆ 95 A β‰ˆ 138 8 β”‚ (Kr, Sr, Zr) (Xe, Cs, Ba) β”‚ __ __ 6 β”‚ / \ / \ β”‚ / \ / \ 4 β”‚ / \ / \ β”‚ / \ / \ 2 β”‚ / \ Valley / \ β”‚____/ \__(A β‰ˆ 117)/ \____ 0 └────┴─────────────┴──────────┴─────────────┴────► Mass Number A 70 95 117 140 160 ```

The Double-Hump Peak Characteristics

For thermal neutron fission of $^{235}\text{U}$:

1. Light Fragment Group: Centered around $A_L \approx 95$ ($Z \approx 36 - 40$, elements $\text{Kr}, \text{Rb}, \text{Sr}, \text{Y}, \text{Zr}$).

2. Heavy Fragment Group: Centered around $A_H \approx 138$ ($Z \approx 53 - 57$, elements $\text{I}, \text{Xe}, \text{Cs}, \text{Ba}, \text{La}$).

3. Symmetric Valley: Symmetric fission ($A \approx 117$) occurs with a probability of less than $0.01\%$ ($\sim 600$ times less frequent than asymmetric fission).

Physical Origin: Quantum Shell Effects

The dominance of asymmetric fission is governed by nuclear shell closures: The heavy fragment peak is anchored by the proximity of the doubly magic spherical shell closure:

$$Z = 50 \text{ (Protons)}, \quad N = 82 \text{ (Neutrons)} \implies ^{132}_{50}\text{Sn}_{82}$$

The extraordinary shell stabilization of nascent fragments near $Z = 50$ and $N = 82$ energetically favors asymmetric neck scission.

At high incident excitation energies ($E_n > 40\text{ MeV}$), individual single-particle shell structures wash out, and the fragment yield curve transitions into a single symmetric Gaussian peak centered at $A/2$.

Β§5.4 Prompt and Delayed Neutrons: Precursor Kinetics & Reactor Control Stability

In every fission event, an average of $\bar{\nu} \approx 2.4 - 3.0$ neutrons are released. Crucially for nuclear engineering, these neutrons are emitted via two distinct physical mechanisms:

1. Prompt Neutrons ($>99\%$)

Emitted directly from the scission neck and evaporating fragments within $\tau_{\text{prompt}} \sim 10^{-14}\text{ seconds}$ of scission. They exhibit a continuous Maxwellian fission spectrum (Watt spectrum):

$$\chi(E) = c \cdot \sqrt{E} \exp(-E / E_0) \quad \text{with average energy } \bar{E} \approx 2.0\text{ MeV}$$

2. Delayed Neutrons ($<1\%$)

Fission fragments are neutron-rich and undergo cascades of negative beta decays. In certain daughter nuclei (termed delayed neutron precursors), the beta decay Q-value exceeds the neutron separation energy ($Q_\beta > S_n$). The daughter is formed in an excited state that promptly ($<10^{-14}\text{ s}$) de-excites by boiling off a neutron:

``` FISSION FRAGMENT PRECURSOR (e.g., ⁸⁷₃₅Br, T₁/β‚‚ = 55.6 s) β”‚ β–Ό (β⁻ decay, slow: governed by T₁/β‚‚) EXCITED DAUGHTER [⁸⁷₃₆Kr] (Excitation E > S_n) β”‚ β–Ό (Prompt neutron emission, fast: <10⁻¹⁴ s) STABLE RESIDUAL ⁸⁢₃₆Kr + ΒΉβ‚€n (DELAYED NEUTRON!) ```

The appearance of the delayed neutron is rate-limited not by nuclear forces, but by the half-life of the preceding beta decay!

Precursor Groups and the Delayed Neutron Fraction ($\beta$)

The fraction of all fission neutrons that are delayed is defined as $\beta$:

$$\beta \equiv \frac{\text{Delayed Neutrons}}{\text{Total Neutrons}} = \begin{cases} 0.0065 \quad (0.65\%) & \text{for } ^{235}\text{U} \\ 0.0021 \quad (0.21\%) & \text{for } ^{239}\text{Pu} \end{cases}$$

Delayed neutrons are categorized into six empirical precursor groups:

| Group $i$ | Representative Precursor | Half-Life $T_{1/2,i}$ | Decay Constant $\lambda_i$ ($\text{s}^{-1}$) | Yield Fraction $\beta_i$ ($^{235}\text{U}$) | | :--- | :--- | :--- | :--- | :--- | | 1 | $^{87}\text{Br}$ | $55.6\text{ s}$ | $0.0124$ | $0.00021$ | | 2 | $^{137}\text{I}$ | $24.5\text{ s}$ | $0.0283$ | $0.00142$ | | 3 | $^{138}\text{I}, ^{89}\text{Br}$ | $16.3\text{ s}$ | $0.0425$ | $0.00127$ | | 4 | $^{139}\text{I}, ^{93}\text{Kr}$ | $5.21\text{ s}$ | $0.1330$ | $0.00257$ | | 5 | $^{140}\text{I}, ^{91}\text{Br}$ | $2.37\text{ s}$ | $0.2920$ | $0.00075$ | | 6 | $^{97}\text{Rb}$ | $0.23\text{ s}$ | $3.0100$ | $0.00028$ | | Total | β€” | β€” | $\bar{\tau}_d \approx 12.7\text{ s}$ | $\beta = 0.00650$ |

Criticality and Reactor Control Physics

Without delayed neutrons, the average neutron generation lifetime in a thermal reactor would be $\Lambda \approx l_p \approx 10^{-4}\text{ seconds}$. If reactivity increased by just $\Delta k = +0.001$, reactor power would escalate as:

$$P(t) = P_0 \exp\left(\frac{\Delta k}{\Lambda} t\right) = P_0 \exp\left(\frac{0.001}{10^{-4}} t\right) = P_0 e^{10 t} \approx P_0 (22,000)^t$$

Power would multiply by 22,000 every single second, rendering mechanical control impossible.

With delayed neutrons, the effective neutron lifetime is dominated by precursor half-lives:

$$\Lambda_{\text{eff}} \approx (1 - \beta) l_p + \beta \bar{\tau}_d \approx (1 - 0.0065)(10^{-4}) + 0.0065(12.7\text{ s}) \approx 0.083\text{ seconds}$$

Power response slows by a factor of nearly $1,000$, enabling mechanical control rods to manage power levels safely.

  • Prompt Criticality ($\rho \ge \beta$ or $k \ge 1 + \beta$): The reactor is critical on prompt neutrons alone; power explodes uncontrollably (Chernobyl scenario).
  • Delayed Criticality ($1.000 \le k < 1 + \beta$): The normal, stable operational regime where criticality requires delayed neutrons.

Β§5.5 Neutron Moderation Kinetics: Elastic Collisions, Logarithmic Decrement & Moderator Selection

Prompt fission neutrons are born with high kinetic energies averaging $E_0 \approx 2\text{ MeV}$. However, the fission cross-section of $^{235}\text{U}$ is hundreds of times higher for thermal neutrons ($E_{\text{th}} \approx 0.025\text{ eV}$). To sustain a thermal chain reaction, fast neutrons must be slowed down via elastic collisions with moderator nuclei.

Elastic Collision Kinematics in the LAB Frame

Consider a neutron of mass $m = 1$ colliding elastically with a stationary moderator nucleus of mass number $A$. By conservation of momentum and energy in the center-of-mass frame, the ratio of final neutron energy $E'$ to initial energy $E$ after scattering through CM angle $\theta$ is:

$$\frac{E'}{E} = \frac{1 + A^2 + 2A \cos\theta}{(1 + A)^2}$$
  • Minimum Energy (Head-On Collision, $\theta = \pi$):
$$\left(\frac{E'}{E}\right)_{\min} = \left(\frac{A - 1}{A + 1}\right)^2 \equiv \alpha$$

where $\alpha = \left(\frac{A - 1}{A + 1}\right)^2$ is the collision parameter.

  • For hydrogen ($A = 1$): $\alpha = 0$. A neutron can transfer $100\%$ of its kinetic energy in a single collision!
  • For carbon ($A = 12$): $\alpha = (11/13)^2 \approx 0.716$. The neutron retains at least $71.6\%$ of its energy.

The Average Logarithmic Energy Decrement ($\xi$)

Because neutron energy loss is multiplicative rather than additive, the slowing-down power is characterized by the average decrease in the natural logarithm of energy per collision:

$$\xi \equiv \left\langle \ln\left(\frac{E}{E'}\right) \right\rangle = \int_\alpha^1 \ln\left(\frac{E}{E'}\right) P\left(\frac{E'}{E}\right) d\left(\frac{E'}{E}\right)$$

For isotropic s-wave scattering in the CM frame, $P(E'/E) = \frac{1}{1 - \alpha}$. Evaluating the integral:

$$\xi = 1 + \frac{(A - 1)^2}{2A} \ln\left(\frac{A - 1}{A + 1}\right) = 1 + \frac{\alpha \ln\alpha}{1 - \alpha}$$

For $A > 10$, this is approximated by:

$$\xi \approx \frac{2}{A + 2/3}$$

Number of Collisions to Thermalize ($N_{\text{coll}}$)

The average number of collisions required to slow a neutron from fission energy $E_0 = 2\text{ MeV}$ to thermal energy $E_{\text{th}} = 0.025\text{ eV}$ is:

$$N_{\text{coll}} = \frac{\ln(E_0 / E_{\text{th}})}{\xi} = \frac{\ln(2 \times 10^6 / 0.025)}{\xi} = \frac{\ln(8.0 \times 10^7)}{\xi} = \frac{18.2}{\xi}$$

| Moderator | Mass Number $A$ | $\alpha = \left(\frac{A-1}{A+1}\right)^2$ | Decrement $\xi$ | Collisions to Thermalize $N_{\text{coll}}$ | Moderating Ratio ($\xi \Sigma_s / \Sigma_a$) | | :--- | :--- | :--- | :--- | :--- | :--- | | Light Water ($\text{H}_2\text{O}$) | $1$ | $0.000$ | $1.000$ | $18$ | $71$ | | Heavy Water ($\text{D}_2\text{O}$) | $2$ | $0.111$ | $0.725$ | $25$ | $5,670$ | | Beryllium ($\text{Be}$) | $9$ | $0.640$ | $0.207$ | $86$ | $143$ | | Graphite ($\text{C}$) | $12$ | $0.716$ | $0.158$ | $115$ | $192$ |

Heavy water ($\text{D}_2\text{O}$) has the highest moderating ratio ($\xi \Sigma_s / \Sigma_a = 5,670$) because deuterium has an exceptionally small neutron capture cross-section ($0.5\text{ mb}$ vs $332\text{ mb}$ for $^1\text{H}$), enabling CANDU reactors to operate using natural, unenriched uranium.

Β§5.6 Nuclear Chain Reactions: The Four-Factor Formula, Six-Factor Formula & Criticality Metrics

The operational state of a nuclear reactor is governed by the effective neutron multiplication factor $k_{\text{eff}}$, defined as:

$$k_{\text{eff}} \equiv \frac{\text{Neutrons produced in generation } n+1}{\text{Neutrons absorbed or lost in generation } n}$$
  • Subcritical ($k_{\text{eff}} < 1$): Chain reaction dies out exponentially; power decreases.
  • Critical ($k_{\text{eff}} = 1.0000$): Steady-state chain reaction; power is exactly constant.
  • Supercritical ($k_{\text{eff}} > 1$): Chain reaction diverges; power increases exponentially.

The Four-Factor Formula for Infinite Media ($k_\infty$)

In an infinitely large reactor (where leakage is zero), the multiplication factor is governed by Fermi's Four-Factor Formula:

$$k_\infty = \eta \cdot \epsilon \cdot p \cdot f$$

``` Fast Neutrons from Thermal Fission: Ξ· * f β”‚ β–Ό Fast Fission Bonus (* Ξ΅) Total Fast Neutrons: Ξ΅ Ξ· f β”‚ β–Ό Resonance Escape (* p) Neutrons Reaching Thermal Energy: p Ξ΅ Ξ· * f β”‚ β–Ό Thermal Utilization (* f) Thermal Neutrons Absorbed in Fuel: f p Ξ΅ Ξ· f ──► Next Generation! ```

1. Reproduction Factor ($\eta$):

Average number of fission neutrons produced per thermal neutron absorbed in the fuel:

$$\eta = \nu \frac{\Sigma_f^F}{\Sigma_a^F} = \nu \frac{\sigma_f^{235} N_{235}}{\sigma_a^{235} N_{235} + \sigma_a^{238} N_{238}}$$

For natural uranium, $\eta \approx 1.34$; for $3.5\%$ enriched fuel, $\eta \approx 1.80$.

2. Fast Fission Factor ($\epsilon$):

Ratio of total fast neutrons (including those from fast fission of $^{238}\text{U}$) to neutrons from thermal fission alone:

$$\epsilon \approx 1.03 - 1.07$$

3. Resonance Escape Probability ($p$):

Probability that a fast neutron slows down through the broad $^{238}\text{U}$ resonance absorption capture peaks ($10\text{ eV} - 1\text{ keV}$) without being absorbed:

$$p = \exp\left( -\frac{N_{238}}{\xi \Sigma_s} I_{\text{eff}} \right) \approx 0.85 - 0.92$$

Using heterogeneous fuel rods separated by moderator increases $p$ because neutrons slow down in the moderator away from $^{238}\text{U}$.

4. Thermal Utilization Factor ($f$):

Fraction of thermal neutrons absorbed in the nuclear fuel compared to total thermal absorptions across fuel, moderator, cladding, and structure:

$$f = \frac{\Sigma_a^{\text{fuel}}}{\Sigma_a^{\text{fuel}} + \Sigma_a^{\text{mod}} + \Sigma_a^{\text{clad}} + \Sigma_a^{\text{poisons}}} \approx 0.88 - 0.95$$

The Six-Factor Formula for Finite Reactors ($k_{\text{eff}}$)

In a real finite reactor core, neutrons can leak out across the boundary:

$$k_{\text{eff}} = k_\infty \cdot P_{NL,f} \cdot P_{NL,\text{th}} = \eta \cdot \epsilon \cdot p \cdot f \cdot P_{NL,f} \cdot P_{NL,\text{th}}$$

where:

  • $P_{NL,f} = \frac{1}{1 + M_f^2 B_g^2} \approx \frac{1}{1 + \tau_F B_g^2}$ is the fast non-leakage probability (Fermi age $\tau_F$).
  • $P_{NL,\text{th}} = \frac{1}{1 + L_{\text{th}}^2 B_g^2}$ is the thermal non-leakage probability ($L_{\text{th}}$ = thermal diffusion length).
  • $B_g^2$ is the geometric buckling of the core geometry (e.g., $B_g^2 = (\pi/H)^2 + (2.405/R)^2$ for a finite cylinder).

The reactivity $\rho$ of a reactor core is defined as:

$$\rho \equiv \frac{k_{\text{eff}} - 1}{k_{\text{eff}}}$$

Reactivity is measured in percent ($\%\Delta k/k$), parts per hundred thousand ($\text{pcm} = 10^{-5}$), or dollars ($\$ = \rho / \beta$).

Β§5.7 Commercial Reactor Architectures: PWR, BWR, CANDU, Fast Breeder Reactors & Safety Systems

Commercial nuclear power generation utilizes distinct reactor designs engineered around specific neutron moderation, cooling, and isotopic fuel cycles.

1. Pressurized Water Reactor (PWR)

The dominant global reactor technology ($>65\%$ of worldwide capacity).

  • Coolant & Moderator: Light water ($\text{H}_2\text{O}$) under high pressure ($15.5\text{ MPa} \approx 155\text{ bar}$) to prevent bulk boiling at operating temperatures of $\sim 315^\circ\text{C}$.
  • Steam Cycle: Two distinct circuits (primary and secondary). The hot radioactive primary water passes through U-tube steam generators, boiling secondary water to drive the steam turbine.
  • Fuel: Low-enriched uranium ($\text{UO}_2$ pellets, $3.0 - 5.0\%$ $^{235}\text{U}$) in Zircaloy cladding.

2. Boiling Water Reactor (BWR)

  • Coolant & Moderator: Light water at lower operating pressure ($7.0\text{ MPa} \approx 70\text{ bar}$).
  • Steam Cycle: Direct single loop. Water boils directly in the reactor core ($12 - 15\%$ steam void fraction at core exit), and radioactive steam flows straight to the turbine.
  • Control Rods: Inserted from the bottom of the pressure vessel (because steam voids at the top reduce moderation, shifting flux to the bottom).

3. CANDU (Canada Deuterium Uranium)

  • Coolant & Moderator: Heavy water ($\text{D}_2\text{O}$) in separate circuits. The low-pressure moderator is contained in a horizontal calandria vessel traversed by hundreds of pressurized fuel channels.
  • Fuel: Natural uranium ($0.72\% \, ^{235}\text{U}$). Enabled by the low neutron capture of deuterium.
  • Refueling: On-line refueling while operating at full power.

4. Liquid Metal Fast Breeder Reactor (LMFBR)

  • Moderator: None (operates on fast neutron spectrum, $E_n > 100\text{ keV}$).
  • Coolant: Liquid sodium ($\text{Na}$) or lead-bismuth eutectic at atmospheric pressure, with thermal conductivity $\sim 100\times$ water.
  • Breeding Cycle: Breeds more fissile $^{239}\text{Pu}$ from $^{238}\text{U}$ blankets than it consumes ($BR > 1.0$), multiplying available nuclear fuel reserves by a factor of 60.

``` REACTOR TYPE MODERATOR COOLANT FUEL ENRICHMENT THERMAL EFFICIENCY ────────────────────────────────────────────────────────────────────────────── PWR Light Water Water (15 MPa) 3.2 - 4.9% ²³⁡U ~ 33 - 34% BWR Light Water Steam/Water 3.0 - 4.5% ²³⁡U ~ 33 - 34% CANDU Heavy Water Heavy Water Natural (0.72%) ~ 30 - 31% LMFBR None (Fast) Liquid Sodium 15 - 20% ²³⁹Pu/U ~ 39 - 41% HTGR Graphite Helium Gas 8 - 15% ²³⁡U ~ 45 - 50% ```

Inherent Passive Safety Principles

Modern Generation III+ and IV reactors incorporate passive safety mechanisms that operate without electrical power or operator intervention:

  • Negative Doppler Temperature Coefficient: As fuel temperature rises, thermal agitation broadens $^{238}\text{U}$ resonance capture peaks (Doppler broadening), absorbing more neutrons and automatically shutting down the chain reaction.
  • Negative Moderator Void Coefficient: Boiling of water removes moderator, decreasing reactivity.

Β§5.8 Advanced Fission Physics: Ternary Fission, Delayed Precursor Chemistry & Actinide Incineration

While binary fission splits a nucleus into two primary fragments, approximately 1 in every 500 thermal fissions of $^{235}\text{U}$ is a ternary fission event in which a third light charged particle (LCP) is emitted simultaneously from the scission neck.

Ternary Fission Mechanics

  • Over $90\%$ of ternary particles are energetic alpha particles ($^{4}\text{He}^{2+}$ with kinetic energy $\bar{E} \approx 16\text{ MeV}$).
  • Other light particles include tritons ($^3\text{H}$, $\sim 7\%$), deuterons ($^2\text{H}$), and trace lithium, beryllium, and carbon ions ($^8\text{Be}, ^{10}\text{Be}, ^{14}\text{C}$).
  • The ternary particle is ejected nearly perpendicular to the fission axis ($90^\circ$) by the intense mutual Coulomb repulsion of the two receding heavy fragments. Ternary fission is the primary source of radioactive tritium ($^3\text{H}$) generated inside nuclear fuel rods.

Precursor Chemistry of Delayed Neutron Emitters

Delayed neutron precursors reside in specific chemical groups within the fission fragment distribution:

1. Halogen Precursors: Bromine ($^{87}\text{Br}, ^{88}\text{Br}, ^{89}\text{Br}$) and Iodine ($^{137}\text{I}, ^{138}\text{I}, ^{139}\text{I}$). They possess high beta-decay energies ($Q_\beta \sim 6 - 8\text{ MeV}$) feeding states above the neutron separation energy $S_n$ of the noble gas daughters ($\text{Kr}$ and $\text{Xe}$).

2. Alkali Precursors: Rubidium ($^{92}\text{Rb}, ^{93}\text{Rb}$) and Cesium ($^{141}\text{Cs}, ^{142}\text{Cs}$).

Minor Actinide Partitioning & Transmutation (Actinide Incineration)

In spent nuclear fuel, the dominant long-term radiotoxicity ($>1,000\text{ years}$) arises from minor actinides: neptunium ($^{237}\text{Np}$), americium ($^{241}\text{Am}, ^{243}\text{Am}$), and curium ($^{244}\text{Cm}$). To eliminate multi-millennial geological repository hazards, advanced fuel cycles develop Partitioning and Transmutation (P&T):

1. Pyrochemical Pyroprocessing: Spent oxide fuel is reduced to metal in molten lithium chloride-potassium chloride ($\text{LiCl-KCl}$) eutectic salt at $500^\circ\text{C}$ and electrorefined, separating actinides from fission products.

2. Fast Reactor & Accelerator-Driven System (ADS) Incineration:

In a fast neutron spectrum, the fission-to-capture cross-section ratio $\sigma_f / \sigma_c$ increases dramatically. Minor actinides are loaded into subcritical fast cores driven by spallation proton accelerators, where fast neutrons fission them into short-lived fission products, reducing waste storage lifespans from 300,000 years to under 300 years!

Comprehensive Comparison of Global Commercial Nuclear Reactor Architectures

| Reactor System | Full Name | Moderator | Coolant | Core Pressure | Outlet Temp | Fuel Type & Enrichment | Worldwide Fleet Share | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | PWR | Pressurized Water Reactor | $\text{H}_2\text{O}$ | $\text{H}_2\text{O}$ (liquid) | $15.5\text{ MPa}$ | $325^\circ\text{C}$ | $\text{UO}_2$ ($3.2 - 4.95\%$) | $68\%$ | | BWR | Boiling Water Reactor | $\text{H}_2\text{O}$ | Steam/Water | $7.2\text{ MPa}$ | $285^\circ\text{C}$ | $\text{UO}_2$ ($3.0 - 4.5\%$) | $15\%$ | | PHWR / CANDU | Pressurized Heavy Water | $\text{D}_2\text{O}$ | $\text{D}_2\text{O}$ | $10.0\text{ MPa}$ | $310^\circ\text{C}$ | Natural Uranium ($0.72\%$) | $11\%$ | | AGR | Advanced Gas-Cooled | Graphite | $\text{CO}_2$ gas | $4.2\text{ MPa}$ | $650^\circ\text{C}$ | $\text{UO}_2$ ($2.5 - 3.5\%$) in Stainless | $2\%$ | | SFR / LMFBR | Sodium Fast Reactor | None | Liquid Sodium | $0.1\text{ MPa}$ | $550^\circ\text{C}$ | $\text{MOX}$ ($15 - 20\% \, \text{Pu}$) | $<1\%$ | | VHTR | Very High Temperature Gas | Graphite | Helium gas | $7.0\text{ MPa}$ | $950 - 1000^\circ\text{C}$ | TRISO fuel particles ($10 - 15\%$) | Emerging Gen-IV |

Easy Example 5.1: Four-Factor Formula and Infinite Multiplication Factor Calculation

A low-enriched thermal reactor fuel lattice has the following neutron parameters:

  • Reproduction factor $\eta = 1.340$
  • Fast fission factor $\epsilon = 1.035$
  • Resonance escape probability $p = 0.890$
  • Thermal utilization factor $f = 0.885$
  1. Calculate the infinite multiplication factor $k_\infty$.
  2. If non-leakage probabilities are $P_{NL,f} = 0.940$ and $P_{NL,\text{th}} = 0.970$, determine the effective multiplication factor $k_{\text{eff}}$.
  3. Determine whether the finite reactor is subcritical, critical, or supercritical, and compute core reactivity $\rho$ in pcm.

Step 1: Infinite Multiplication Factor $k_\infty$

Using the Four-Factor Formula:

$$k_\infty = \eta \cdot \epsilon \cdot p \cdot f$$
$$k_\infty = 1.340 \times 1.035 \times 0.890 \times 0.885$$

Multiply step-by-step:

$$\eta \cdot \epsilon = 1.340 \times 1.035 = 1.38690$$
$$1.38690 \times 0.890 = 1.23434$$
$$k_\infty = 1.23434 \times 0.885 \approx 1.0924$$

Step 2: Effective Multiplication Factor $k_{\text{eff}}$

Using the Six-Factor Formula:

$$k_{\text{eff}} = k_\infty \cdot P_{NL,f} \cdot P_{NL,\text{th}}$$
$$k_{\text{eff}} = 1.0924 \times 0.940 \times 0.970 = 1.0924 \times 0.9118 \approx 0.99605$$

Step 3: Criticality Evaluation and Reactivity

Because $k_{\text{eff}} = 0.99605 < 1.0000$, the finite reactor is subcritical. Reactivity $\rho$:

$$\rho = \frac{k_{\text{eff}} - 1}{k_{\text{eff}}} = \frac{0.99605 - 1.0000}{0.99605} = \frac{-0.00395}{0.99605} \approx -0.003966$$

In percent:

$$\rho = -0.397\% \, \Delta k/k$$

In parts per hundred thousand ($\text{pcm} = 10^{-5}$):

$$\rho = -0.003966 \times 10^5 \approx -397\text{ pcm}$$

The reactor has a subcritical shutdown margin of $397\text{ pcm}$.

Easy Example 5.2: Bohr-Wheeler Fissility Parameter and Critical Limit Comparison

Using the SEMF surface coefficient $a_s = 17.80\text{ MeV}$ and Coulomb coefficient $a_c = 0.711\text{ MeV}$:

  1. Calculate the critical fissility ratio $(Z^2/A)_{\text{crit}}$.
  2. Compute the fissility parameter $x$ for uranium-235 ($^{235}_{92}\text{U}$) and californium-252 ($^{252}_{98}\text{Cf}$).
  3. Explain why $^{252}\text{Cf}$ has a high spontaneous fission branch ($3.09\%$) compared to $^{235}\text{U}$ ($7 \times 10^{-9}\%$).

Step 1: Critical Fissility Parameter

The liquid drop stability threshold against spontaneous deformation is:

$$\left(\frac{Z^2}{A}\right)_{\text{crit}} = \frac{2 a_s}{a_c} = \frac{2(17.80\text{ MeV})}{0.711\text{ MeV}} \approx 50.07$$

Step 2: Compute Fissility Parameter $x$

  1. For Uranium-235 ($Z = 92, A = 235$):
$$\frac{Z^2}{A} = \frac{92^2}{235} = \frac{8464}{235} \approx 36.017$$
$$x(^{235}\text{U}) = \frac{36.017}{50.07} \approx 0.719$$
  1. For Californium-252 ($Z = 98, A = 252$):
$$\frac{Z^2}{A} = \frac{98^2}{252} = \frac{9604}{252} \approx 38.111$$
$$x(^{252}\text{Cf}) = \frac{38.111}{50.07} \approx 0.761$$

Step 3: Physical Explanation

The fission barrier height $E_f$ decreases with fissility parameter $x$:

$$E_f \approx 0.83(1 - x)^3 E_S^0$$

For $^{235}\text{U}$, $x \approx 0.72 \implies E_f \approx 5.7\text{ MeV}$. The quantum tunneling probability is extremely small, giving a spontaneous fission half-life of $10^{19}\text{ years}$. For $^{252}\text{Cf}$, $x \approx 0.76 \implies E_f \approx 3.7\text{ MeV}$. The barrier is $2\text{ MeV}$ lower and narrower, increasing quantum tunneling probability by 11 orders of magnitude! Californium-252 undergoes spontaneous fission with $T_{1/2,\text{SF}} = 85.5\text{ years}$, emitting $3.77$ neutrons per fission.

Intermediate Example 5.3: Nuclear Reactor Uranium-235 Burnup and Thermal Power Output

A commercial nuclear power plant operates at a steady electrical output of $P_e = 1,000\text{ MWe}$ with a net thermodynamic thermal efficiency $\eta_{\text{th}} = 33.3\%$. Each fission of $^{235}\text{U}$ releases an average recoverable thermal energy $E_{\text{fiss}} = 200\text{ MeV}$.

  1. Calculate the core thermal power $P_{\text{th}}$ in megawatts (MWth).
  2. Determine the required number of fission events per second.
  3. Calculate the mass of $^{235}\text{U}$ consumed (fissioned) per day in kilograms, and over a full operating year ($365\text{ days}$).

Step 1: Core Thermal Power

$$P_{\text{th}} = \frac{P_e}{\eta_{\text{th}}} = \frac{1000\text{ MWe}}{0.33333} \approx 3,000\text{ MWth} = 3.00 \times 10^9\text{ J/s}$$

Step 2: Fissions per Second

Energy released per fission:

$$E_{\text{fiss}} = 200\text{ MeV} \times 1.60218 \times 10^{-13}\text{ J/MeV} = 3.20436 \times 10^{-11}\text{ J}$$

Fission rate $\dot{N}_f$:

$$\dot{N}_f = \frac{P_{\text{th}}}{E_{\text{fiss}}} = \frac{3.00 \times 10^9\text{ J/s}}{3.20436 \times 10^{-11}\text{ J/fission}} \approx 9.362 \times 10^{19}\text{ fissions/second}$$

Step 3: Fissioned Mass per Day and Year

Number of fissions per day ($86,400\text{ s}$):

$$N_{\text{day}} = (9.362 \times 10^{19}\text{ s}^{-1}) \times 86400\text{ s} \approx 8.089 \times 10^{24}\text{ fissions/day}$$

Moles of $^{235}\text{U}$:

$$n = \frac{8.089 \times 10^{24}}{6.02214 \times 10^{23}\text{ mol}^{-1}} \approx 13.432\text{ moles/day}$$

Mass of $^{235}\text{U}$ fissioned per day:

$$m_{\text{day}} = 13.432\text{ mol} \times 235.04\text{ g/mol} \approx 3,157\text{ g/day} \approx 3.16\text{ kg/day}$$

Annual burnup:

$$m_{\text{year}} = 3.157\text{ kg/day} \times 365\text{ days} \approx 1,152\text{ kg/year} \approx 1.15\text{ metric tons/year}$$

A $1,000\text{ MWe}$ coal power plant burns over $3,000,000\text{ metric tons}$ of coal per year to produce the same energy as just $1.15\text{ tons}$ of fissioned uranium!

Advanced Example 5.4: Delayed Neutron Inhour Equation and Reactor Period

A critical thermal reactor operating at $k = 1.0000$ experiences a small step reactivity insertion $\rho = +0.0010$ ($100\text{ pcm}$). The delayed neutron fraction is $\beta = 0.0065$ and the average precursor decay constant is $\bar{\lambda} = 0.080\text{ s}^{-1}$ ($\bar{\tau} = 12.5\text{ s}$). The prompt neutron generation time is $\Lambda = 1.0 \times 10^{-4}\text{ s}$.

  1. Formulate the simplified one-delayed-group inhour equation for asymptotic reactor period $T$.
  2. Calculate the reactor period $T$ in seconds.
  3. Determine the reactor power multiplication factor after 60 seconds.

Step 1: One-Group Inhour Equation

The inhour equation relates reactivity $\rho$ to the asymptotic reactor period $T$ (where $P(t) = P_0 e^{t/T}$):

$$\rho = \frac{\Lambda}{T + \Lambda} + \sum_{i=1}^6 \frac{\beta_i}{1 + \lambda_i T}$$

In the single-delayed-group approximation for small reactivity ($\rho \ll \beta$ and $T \gg \Lambda$):

$$\rho \approx \frac{\Lambda}{T} + \frac{\beta}{1 + \bar{\lambda} T}$$

Because $\bar{\lambda} T$ is typically $\gg 1$:

$$\rho \approx \frac{\beta}{\bar{\lambda} T} \implies T \approx \frac{\beta - \rho}{\bar{\lambda} \rho} \approx \frac{\beta}{\bar{\lambda} \rho}$$

More precisely, solving $\rho (1 + \bar{\lambda} T) \approx \beta$:

$$\rho + \rho \bar{\lambda} T = \beta \implies T = \frac{\beta - \rho}{\rho \bar{\lambda}}$$

Step 2: Calculate Reactor Period $T$

Substitute values ($\rho = 0.0010$, $\beta = 0.0065$, $\bar{\lambda} = 0.080\text{ s}^{-1}$):

$$\beta - \rho = 0.0065 - 0.0010 = 0.0055$$
$$\rho \bar{\lambda} = 0.0010 \times 0.080\text{ s}^{-1} = 8.0 \times 10^{-5}\text{ s}^{-1}$$
$$T = \frac{0.0055}{8.0 \times 10^{-5}\text{ s}^{-1}} = 68.75\text{ seconds}$$

The reactor power rises with an asymptotic period of $T \approx 68.8\text{ seconds}$.

Step 3: Power Multiplication After 60 Seconds

$$\frac{P(60)}{P_0} = e^{t / T} = \exp\left(\frac{60\text{ s}}{68.75\text{ s}}\right) = e^{0.8727} \approx 2.39$$

Reactor power increases by a manageable factor of $2.39$ over one minute, allowing operator and automated rod adjustments.

Advanced Example 5.5: Fission Product Poisoning Kinetics: Iodine-135 and Xenon-135 Pit

Xenon-135 is the strongest thermal neutron poison known ($\sigma_a = 2.65 \times 10^6\text{ barns}$). In a reactor operating at thermal neutron flux $\Phi = 1.0 \times 10^{14}\text{ n/cm}^2\cdot\text{s}$:

  • Fission yield of $^{135}\text{I}$: $\gamma_I = 0.0639$ ($T_{1/2,I} = 6.57\text{ h} \implies \lambda_I = 2.93 \times 10^{-5}\text{ s}^{-1}$)
  • Direct fission yield of $^{135}\text{Xe}$: $\gamma_X = 0.0023$ ($T_{1/2,X} = 9.14\text{ h} \implies \lambda_X = 2.11 \times 10^{-5}\text{ s}^{-1}$)
  1. Formulate the steady-state concentrations of $^{135}\text{I}$ and $^{135}\text{Xe}$ in terms of macroscopic fission cross-section $\Sigma_f$.
  2. Explain the "xenon pit" (xenon poisoning peak) that occurs several hours after a sudden reactor shutdown.

Step 1: Steady-State Concentrations

1. Iodine-135 Balance:

At steady state, production by fission equals radioactive decay:

$$\gamma_I \Sigma_f \Phi = \lambda_I N_I \implies N_I^0 = \frac{\gamma_I \Sigma_f \Phi}{\lambda_I}$$

2. Xenon-135 Balance:

Production occurs via direct fission ($\gamma_X \Sigma_f \Phi$) and decay of $^{135}\text{I}$ ($\lambda_I N_I$). Disappearance occurs via radioactive decay ($\lambda_X N_X$) and neutron burnup absorption ($\sigma_a^X \Phi N_X$):

$$\gamma_X \Sigma_f \Phi + \lambda_I N_I = \lambda_X N_X + \sigma_a^X \Phi N_X$$

Substituting $\lambda_I N_I = \gamma_I \Sigma_f \Phi$:

$$(\gamma_X + \gamma_I) \Sigma_f \Phi = (\lambda_X + \sigma_a^X \Phi) N_X$$
$$N_X^0 = \frac{(\gamma_I + \gamma_X) \Sigma_f \Phi}{\lambda_X + \sigma_a^X \Phi}$$

Step 2: The Xenon Pit (Shutdown Transient)

Upon reactor shutdown, neutron flux drops to zero ($\Phi \to 0$):

1. Loss of Burnup: The primary destruction channel for xenonβ€”neutron burnup ($\sigma_a^X \Phi N_X$)β€”instantly vanishes.

2. Continued Production: The massive reservoir of accumulated iodine-135 ($N_I^0 \gg N_X^0$) continues to decay into $^{135}\text{Xe}$ with half-life $6.57\text{ hours}$.

3. Transient Peak: Xenon-135 concentration builds up rapidly to a maximum at:

$$t_{\max} = \frac{1}{\lambda_X - \lambda_I} \ln\left[ \frac{\lambda_X}{\lambda_I} \left(1 - \frac{\lambda_X - \lambda_I}{\lambda_X} \frac{N_X^0}{N_I^0}\right) \right] \approx 10 - 11\text{ hours}$$

At $t \approx 11\text{ hours}$, xenon negative reactivity reaches a deep maximum ("xenon pit"). If control margins are inadequate, the reactor cannot be restarted until the xenon decays away ($\sim 30 - 40\text{ hours}$ later). Attempting an improper restart during a xenon transient contributed to the 1986 Chernobyl disaster.

Intermediate Example 5.6: Breeding Ratio in Liquid Metal Fast Breeder Reactor

A sodium-cooled Fast Breeder Reactor operates with mixed oxide fuel ($^{239}\text{Pu}\text{O}_2 / ^{238}\text{U}\text{O}_2$). In the fast neutron spectrum:

  • Average neutrons per fission of $^{239}\text{Pu}$: $\nu = 2.92$
  • Ratio of capture to fission cross-sections for $^{239}\text{Pu}$: $\alpha_c = \sigma_c / \sigma_f = 0.15$
  • Fractional parasitic absorption in structure/sodium: $L_p = 0.18$
  • Core leakage fraction: $L = 0.08$
  1. Calculate the reproduction factor $\eta = \nu / (1 + \alpha_c)$ for $^{239}\text{Pu}$ in this spectrum.
  2. Formulate and compute the breeding ratio $BR$ (excess fissile nuclei produced per fissile nucleus destroyed).
  3. If $BR = 1.25$, calculate the fuel doubling time $T_D$ in years for a specific inventory of $3.0\text{ kg/MWe}$ and capacity factor $85\%$.

Step 1: Reproduction Factor $\eta$

Thermal neutrons yield $\eta \approx 2.11$ for Pu-239; in a fast spectrum, $\nu$ rises and capture $\alpha_c$ drops:

$$\eta = \frac{\nu}{1 + \alpha_c} = \frac{2.92}{1 + 0.15} = \frac{2.92}{1.15} \approx 2.539$$

Step 2: Breeding Ratio Formulation

Of the $\eta$ neutrons produced per fissile $^{239}\text{Pu}$ destroyed:

  • Exactly $1.00$ neutron must be absorbed in $^{239}\text{Pu}$ to sustain the chain reaction.
  • $L_p$ neutrons are lost to parasitic capture.
  • $L$ neutrons leak out of the blanket.

The remaining neutrons are captured by fertile $^{238}\text{U}$ to breed $^{239}\text{Pu}$:

$$BR = \eta - 1 - L_p - L$$
$$BR = 2.539 - 1.000 - 0.180 - 0.080 = 1.279$$

Because $BR = 1.28 > 1.00$, the reactor breeds $28\%$ more fuel than it consumes!

Step 3: Doubling Time $T_D$

For $BR = 1.25$, the net breeding gain is $G = BR - 1 = 0.25$. The fuel doubling time (simple compound formulation) is:

$$T_D = \frac{M_{\text{inv}}}{G \cdot \dot{M}_{\text{fiss}} \cdot CF}$$

At $1\text{ MWe}$ ($3\text{ MWth}$), annual fission burnup is $\approx 1.15\text{ kg/yr}$. Annual fissile surplus:

$$\Delta M = 0.25 \times 1.15\text{ kg/yr} \times 0.85 \approx 0.244\text{ kg/yr}$$

Doubling time for $3.0\text{ kg}$ inventory:

$$T_D = \frac{3.0\text{ kg}}{0.244\text{ kg/yr}} \approx 12.3\text{ years}$$

The reactor doubles its initial fuel load in $\approx 12\text{ years}$.

Easy Example 5.7: CANDU Natural Uranium Heavy Water Moderation Optimization

In a CANDU reactor utilizing natural uranium fuel ($0.720\% \, ^{235}\text{U}$, $99.280\% \, ^{238}\text{U}$):

  • Thermal capture cross-section of $^{235}\text{U}$: $\sigma_a^{235} = 680\text{ b}$ ($\sigma_f^{235} = 585\text{ b}$)
  • Thermal capture cross-section of $^{238}\text{U}$: $\sigma_a^{238} = 2.70\text{ b}$
  • Average neutrons per thermal fission of $^{235}\text{U}$: $\nu = 2.42$
  1. Calculate the reproduction factor $\eta$ for natural uranium.
  2. Explain why a light water moderated reactor cannot achieve criticality with natural uranium, whereas a heavy water moderated reactor can.

Step 1: Reproduction Factor $\eta$

The effective absorption cross-section per atom of natural uranium is:

$$\bar{\sigma}_a = 0.00720(680\text{ b}) + 0.99280(2.70\text{ b}) = 4.896 + 2.681 = 7.577\text{ barns}$$

The effective fission cross-section is:

$$\bar{\sigma}_f = 0.00720(585\text{ b}) = 4.212\text{ barns}$$

Reproduction factor:

$$\eta = \nu \frac{\bar{\sigma}_f}{\bar{\sigma}_a} = 2.42 \times \frac{4.212}{7.577} = 2.42 \times 0.5559 \approx 1.345$$

Step 2: Physical Explanation of Moderator Choice

With $\eta \approx 1.345$, the product $\epsilon \cdot p \cdot f$ in the Four-Factor Formula must satisfy:

$$\epsilon \cdot p \cdot f \ge \frac{1}{\eta} = \frac{1}{1.345} \approx 0.743$$
  • In Light Water ($\text{H}_2\text{O}$):

Hydrogen has a significant thermal neutron capture cross-section ($\sigma_a(H) = 0.332\text{ barns}$). The thermal utilization factor $f$ drops severely ($f \sim 0.70$), driving $k_\infty = \eta \epsilon p f \approx 1.345 \times 1.03 \times 0.85 \times 0.70 \approx 0.825 < 1.000$. Light water absorbs too many neutrons to achieve criticality with natural uranium.

  • In Heavy Water ($\text{D}_2\text{O}$):

Deuterium has a negligible capture cross-section ($\sigma_a(D) = 0.0005\text{ barns}$β€”nearly $700\times$ lower than hydrogen). Thermal utilization remains high ($f \approx 0.94$), and with optimal lattice pitch, $p \approx 0.90$, yielding:

$$k_\infty \approx 1.345 \times 1.03 \times 0.90 \times 0.94 \approx 1.17 > 1.00$$

Heavy water easily sustains criticality with unenriched natural uranium!

Intermediate Example 5.8: Thermal Utilization and Heterogeneous Fuel Rod Lattice Pitch Optimization

In a heterogeneous graphite-moderated reactor core: Fuel rods of natural uranium metal ($V_{\text{fuel}} = 1.00\text{ L}$) are arranged in a lattice with graphite moderator volume $V_{\text{mod}} = 45.0\text{ L}$. Thermal absorption parameters:

  • Uranium fuel: $\Sigma_a^{\text{fuel}} = 0.367\text{ cm}^{-1}$
  • Graphite moderator: $\Sigma_a^{\text{mod}} = 0.000385\text{ cm}^{-1}$
  • Average thermal neutron flux ratio in fuel to moderator: $\bar{\Phi}_{\text{fuel}} / \bar{\Phi}_{\text{mod}} = 0.720$ (flux depression factor)
  1. Formulate the thermal utilization factor $f$ incorporating heterogeneous flux depression:
$$f = \frac{\Sigma_a^{\text{fuel}} V_{\text{fuel}} \bar{\Phi}_{\text{fuel}}}{\Sigma_a^{\text{fuel}} V_{\text{fuel}} \bar{\Phi}_{\text{fuel}} + \Sigma_a^{\text{mod}} V_{\text{mod}} \bar{\Phi}_{\text{mod}}}$$
  1. Calculate the thermal utilization factor $f$.
  2. If the lattice spacing is increased such that $V_{\text{mod}} = 65.0\text{ L}$, compute the new thermal utilization factor and describe the physical tradeoff with resonance escape probability $p$.

Step 1: Thermal Utilization Formulation

Dividing numerator and denominator by $\Sigma_a^{\text{fuel}} V_{\text{fuel}} \bar{\Phi}_{\text{fuel}}$:

$$f = \frac{1}{1 + \left(\frac{\Sigma_a^{\text{mod}}}{\Sigma_a^{\text{fuel}}}\right) \left(\frac{V_{\text{mod}}}{V_{\text{fuel}}}\right) \left(\frac{\bar{\Phi}_{\text{mod}}}{\bar{\Phi}_{\text{fuel}}}\right)}$$

Step 2: Numerical Calculation for $V_{\text{mod}} = 45.0\text{ L}$

Given:

  • $\Sigma_a^{\text{mod}} / \Sigma_a^{\text{fuel}} = 0.000385 / 0.367 \approx 1.049 \times 10^{-3}$
  • $V_{\text{mod}} / V_{\text{fuel}} = 45.0 / 1.00 = 45.0$
  • $\bar{\Phi}_{\text{mod}} / \bar{\Phi}_{\text{fuel}} = 1 / 0.720 \approx 1.3889$

Product:

$$\text{Term} = (1.049 \times 10^{-3}) \times (45.0) \times (1.3889) \approx 0.06556$$
$$f = \frac{1}{1 + 0.06556} \approx 0.9385 \quad (93.85\%)$$

Step 3: Calculation for $V_{\text{mod}} = 65.0\text{ L}$

$$\text{Term} = (1.049 \times 10^{-3}) \times (65.0) \times (1.3889) \approx 0.09470$$
$$f' = \frac{1}{1 + 0.09470} \approx 0.9135 \quad (91.35\%)$$

The thermal utilization drops by $2.5\%$. Physical Tradeoff: Increasing moderator volume decreases $f$ (more parasitic absorption in graphite), but increases the resonance escape probability $p$ (neutrons slow down safely in the moderator without encountering $^{238}\text{U}$ resonance capture peaks). Reactor engineers optimize lattice pitch where the product $p \cdot f$ reaches its global maximum!

Intermediate Example 5.9: Thermal Reactor Critical Size and Buckling for Spherical Core

A bare, unreflected spherical nuclear reactor core has an infinite multiplication factor $k_\infty = 1.080$ and a neutron migration area $M^2 = 32.0\text{ cm}^2$.

  1. Using one-group diffusion theory ($k_{\text{eff}} = \frac{k_\infty}{1 + M^2 B_g^2} = 1.000$), calculate the required material buckling $B_m^2$ in $\text{cm}^{-2}$.
  2. For a spherical reactor of radius $R$, the geometric buckling is $B_g^2 = (\pi / R)^2$. Determine the critical radius $R_{\text{crit}}$ in centimeters and critical volume $V_{\text{crit}}$ in cubic meters.

Step 1: Material Buckling $B_m^2$

At exact criticality ($k_{\text{eff}} = 1.000$):

$$\frac{k_\infty}{1 + M^2 B_m^2} = 1.000 \implies 1 + M^2 B_m^2 = k_\infty$$
$$B_m^2 = \frac{k_\infty - 1}{M^2}$$

Given $k_\infty = 1.080$ and $M^2 = 32.0\text{ cm}^2$:

$$B_m^2 = \frac{1.080 - 1.000}{32.0\text{ cm}^2} = \frac{0.080}{32.0} = 0.00250\text{ cm}^{-2}$$

Step 2: Critical Radius and Volume

For a sphere:

$$B_g^2 = \left(\frac{\pi}{R}\right)^2 = B_m^2 = 0.00250\text{ cm}^{-2}$$

Taking the square root:

$$\frac{\pi}{R_{\text{crit}}} = \sqrt{0.00250} = 0.0500\text{ cm}^{-1}$$
$$R_{\text{crit}} = \frac{\pi}{0.0500\text{ cm}^{-1}} = \frac{3.14159}{0.0500} \approx 62.83\text{ cm}$$

Critical volume:

$$V_{\text{crit}} = \frac{4}{3}\pi R_{\text{crit}}^3 = \frac{4}{3}\pi (62.83\text{ cm})^3 = \frac{4}{3}\pi (248,091\text{ cm}^3) \approx 1,039,200\text{ cm}^3 \approx 1.039\text{ m}^3$$

The critical core has a radius of $62.8\text{ cm}$ and volume of $1.04\text{ m}^3$.

Solved Honors Problems & Derivations

Step-by-step rigorous solutions with full physical, thermodynamic, and nuclear kinematic validation.