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:
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:
``` 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:
- Coulomb Self-Energy: Decreases with deformation (charge is pushed further apart):
The net change in nuclear potential energy $\Delta E$ relative to the spherical ground state is:
``` 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:
We define the dimensionless Fissility Parameter $x$:
Substituting the SEMF expressions $E_S^0 = a_s A^{2/3}$ and $E_C^0 = a_c Z^2 / A^{1/3}$:
Using $a_c \approx 0.711\text{ MeV}$ and $a_s \approx 17.80\text{ MeV}$:
- 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:
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):
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$:
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:
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:
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:
- Minimum Energy (Head-On Collision, $\theta = \pi$):
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:
For isotropic s-wave scattering in the CM frame, $P(E'/E) = \frac{1}{1 - \alpha}$. Evaluating the integral:
For $A > 10$, this is approximated by:
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:
| 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:
- 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:
``` 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:
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:
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:
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:
The Six-Factor Formula for Finite Reactors ($k_{\text{eff}}$)
In a real finite reactor core, neutrons can leak out across the boundary:
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:
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 |
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$
- Calculate the infinite multiplication factor $k_\infty$.
- 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}}$.
- 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:
Multiply step-by-step:
Step 2: Effective Multiplication Factor $k_{\text{eff}}$
Using the Six-Factor Formula:
Step 3: Criticality Evaluation and Reactivity
Because $k_{\text{eff}} = 0.99605 < 1.0000$, the finite reactor is subcritical. Reactivity $\rho$:
In percent:
In parts per hundred thousand ($\text{pcm} = 10^{-5}$):
The reactor has a subcritical shutdown margin of $397\text{ pcm}$.
Using the SEMF surface coefficient $a_s = 17.80\text{ MeV}$ and Coulomb coefficient $a_c = 0.711\text{ MeV}$:
- Calculate the critical fissility ratio $(Z^2/A)_{\text{crit}}$.
- Compute the fissility parameter $x$ for uranium-235 ($^{235}_{92}\text{U}$) and californium-252 ($^{252}_{98}\text{Cf}$).
- 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:
Step 2: Compute Fissility Parameter $x$
- For Uranium-235 ($Z = 92, A = 235$):
- For Californium-252 ($Z = 98, A = 252$):
Step 3: Physical Explanation
The fission barrier height $E_f$ decreases with fissility parameter $x$:
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.
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}$.
- Calculate the core thermal power $P_{\text{th}}$ in megawatts (MWth).
- Determine the required number of fission events per second.
- 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
Step 2: Fissions per Second
Energy released per fission:
Fission rate $\dot{N}_f$:
Step 3: Fissioned Mass per Day and Year
Number of fissions per day ($86,400\text{ s}$):
Moles of $^{235}\text{U}$:
Mass of $^{235}\text{U}$ fissioned per day:
Annual burnup:
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!
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}$.
- Formulate the simplified one-delayed-group inhour equation for asymptotic reactor period $T$.
- Calculate the reactor period $T$ in seconds.
- 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}$):
In the single-delayed-group approximation for small reactivity ($\rho \ll \beta$ and $T \gg \Lambda$):
Because $\bar{\lambda} T$ is typically $\gg 1$:
More precisely, solving $\rho (1 + \bar{\lambda} T) \approx \beta$:
Step 2: Calculate Reactor Period $T$
Substitute values ($\rho = 0.0010$, $\beta = 0.0065$, $\bar{\lambda} = 0.080\text{ s}^{-1}$):
The reactor power rises with an asymptotic period of $T \approx 68.8\text{ seconds}$.
Step 3: Power Multiplication After 60 Seconds
Reactor power increases by a manageable factor of $2.39$ over one minute, allowing operator and automated rod adjustments.
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}$)
- Formulate the steady-state concentrations of $^{135}\text{I}$ and $^{135}\text{Xe}$ in terms of macroscopic fission cross-section $\Sigma_f$.
- 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:
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$):
Substituting $\lambda_I N_I = \gamma_I \Sigma_f \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:
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.
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$
- Calculate the reproduction factor $\eta = \nu / (1 + \alpha_c)$ for $^{239}\text{Pu}$ in this spectrum.
- Formulate and compute the breeding ratio $BR$ (excess fissile nuclei produced per fissile nucleus destroyed).
- 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:
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}$:
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:
At $1\text{ MWe}$ ($3\text{ MWth}$), annual fission burnup is $\approx 1.15\text{ kg/yr}$. Annual fissile surplus:
Doubling time for $3.0\text{ kg}$ inventory:
The reactor doubles its initial fuel load in $\approx 12\text{ years}$.
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$
- Calculate the reproduction factor $\eta$ for natural uranium.
- 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:
The effective fission cross-section is:
Reproduction factor:
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:
- 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:
Heavy water easily sustains criticality with unenriched natural uranium!
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)
- Formulate the thermal utilization factor $f$ incorporating heterogeneous flux depression:
- Calculate the thermal utilization factor $f$.
- 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}}$:
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:
Step 3: Calculation for $V_{\text{mod}} = 65.0\text{ L}$
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!
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$.
- 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}$.
- 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$):
Given $k_\infty = 1.080$ and $M^2 = 32.0\text{ cm}^2$:
Step 2: Critical Radius and Volume
For a sphere:
Taking the square root:
Critical volume:
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.