Monte Carlo Fundamentals

❮ Probability TheoryGeometry ❯

Transport Theory

The Boltzmann Equation: The Master Equation

Boltzmann Transport Equation
Ω^⋅∇ψ+Σt ψ=∫ ⁣ ⁣∫Σs(E′ ⁣→E, Ω^′ ⁣→Ω^) ψ(E′,Ω^′) dE′ dΩ^′+S\hat{\Omega}\cdot\nabla\psi + \Sigma_t\,\psi = \int\!\!\int \Sigma_s(E'\!\to E,\,\hat{\Omega}'\!\to\hat{\Omega})\,\psi(E',\hat{\Omega}')\,dE'\,d\hat{\Omega}' + S

Neutron balance: streaming + removal = scattering source + external source

  • Ω⃗·∇ψ: Neutrons streaming into or out of a point due to their motion
  • Σₜψ: Neutrons removed from this energy and direction by collisions
  • ∫∫ Σₛψ dE'dΩ⃗': Neutrons added by scattering from other energies and directions
  • S: Neutrons added by external sources (fission, fixed sources)

Monte Carlo codes solve this equation by tracking individual neutrons as they stream, scatter, and cause fission. The statistical average of many histories gives the solution — no spatial, angular, or energy discretization needed.

Deeper explanation: why Monte Carlo solves this naturally

Every Monte Carlo neutron transport calculation is solving the Boltzmann equation. Instead of discretizing the equation mathematically, codes track individual neutrons as they stream through space, collide with nuclei, scatter to new directions and energies, and create new neutrons through fission. Each particle event directly represents a term in the equation.

This is why Monte Carlo methods are so powerful for complex geometries — they follow the physics directly. No mesh, no energy groups, no angular quadrature sets.

Why this matters for code users

Understanding the Boltzmann equation helps you recognize when approximations might be valid. In a beam port, streaming dominates (left side of the equation). In a thermal reactor core, the scattering integral is crucial. This insight helps you choose appropriate physics options.

Cross-Sections: The Language of Nuclear Interactions

Macroscopic Cross-Section
Σ=N⋅σ\Sigma = N \cdot \sigma

N = number density, σ = microscopic cross-section

Mean Free Path
λ=1Σt\lambda = \frac{1}{\Sigma_t}

Average distance a neutron travels between collisions

Microscopic σ is an intrinsic property of the nucleus. Macroscopic Σ depends on how many nuclei are packed into a volume. Typical thermal mean free paths: water ~0.4 cm, uranium ~1.5 cm, graphite ~2.6 cm. Monte Carlo codes sample collision distances from exp(-s/λ).

Enrichment, resonances, and energy dependence

Natural uranium absorbs thermal neutrons readily, but most absorptions occur parasitically in U-238 rather than causing fission in U-235. Enriching uranium improves the fissile-to-parasitic absorption ratio, which is why enrichment is so important in reactor design.

The mean free path helps you understand reactor design intuitively: it tells you how much moderator you need between fuel pins. If pins are closer than several mean free paths, fast neutrons hit fuel before being fully moderated.

Energy dependence

At thermal energies (~0.025 eV), many nuclei follow the "1/v law" — absorption cross-sections inversely proportional to velocity. At higher energies, cross-sections exhibit resonances. U-238 has major resonances at 6.7, 20.9, and 36.7 eV; 15-30% of slowing-down neutrons in PWRs are captured in these peaks. This is why Monte Carlo tracks energy continuously rather than using groups.

Criticality: The Heart of Reactor Physics

k-effective Definition
keff=Neutron production in generation n+1Neutron production in generation nk_{\text{eff}} = \frac{\text{Neutron production in generation } n+1}{\text{Neutron production in generation } n}

k > 1 → supercritical; k < 1 → subcritical; k = 1 → critical

Monte Carlo codes calculate keff by tracking neutrons through many fission generations and measuring how many new neutrons each generation produces. This automatically accounts for fast fission, resonance absorption, leakage, and spectrum effects.

Four-Factor Formula (Thermal Reactors)
k∞=η⋅f⋅p⋅εk_{\infty} = \eta \cdot f \cdot p \cdot \varepsilon

Reproduction factor × thermal utilization × resonance escape × fast fission factor

  • η (eta): Neutrons produced per thermal absorption in fuel (U-235: ~2.07, Pu-239: ~2.11)
  • f: Fraction of thermal absorptions in fuel vs. moderator/structure
  • p: Fraction of fast neutrons escaping resonance capture (0.7–0.85 in PWRs)
  • ε: Fast fission bonus (typically 1.02–1.08)

For a finite reactor to be critical, keff must equal 1. Because leakage removes neutrons, k∞ must be greater than 1 — the excess above unity is what leakage consumes.

Design implications and criticality safety

Enriching uranium increases η. Optimizing fuel-to-moderator ratio increases f. Reducing parasitic absorbers increases f. Heterogeneous fuel design increases p by reducing resonance self-shielding.

For fresh fuel at beginning of life, k > 1 is needed to compensate for control rods, fission product buildup, and burnup. For spent fuel storage, k < 1 with adequate margin.

Criticality safety

The regulatory requirement is often k + 2σ < 0.95, where σ is the Monte Carlo statistical uncertainty. You're not just calculating a number — you're calculating a number with an uncertainty, and that uncertainty directly affects safety margins.

The Monte Carlo Algorithm

  1. Sample collision distance: Exponential distribution with mean free path λ = 1/Σₜ
  2. Move the particle: Transport in current direction
  3. Determine collision type: Scatter, absorb, or fission — sampled from relative cross-sections
  4. Process the collision: Scattering → new energy/direction; fission → new neutrons
  5. Repeat: Continue until absorbed or escaped
When transport theory matters most
  • Shielding: Attenuation and buildup factors for shield design
  • Reactor control: Control rod worth depends on spectrum and spatial distribution
  • Fuel management: Burnup changes spectrum and reactivity
  • Safety analysis: Reactivity changes with temperature, voiding, or geometry

Quick Reference: Physics Options

Thermal reactors

Enable thermal scattering with appropriate S(α,β) libraries. Water molecule dynamics significantly affect thermal neutron behavior.

Fast reactors

Inelastic scattering models and cross-section data above 1 MeV are critical.

Shielding

Include photon transport for dose — capture gammas often dominate behind thick shields.

Interpreting results and building intuition
  • Large k-eff changes from small geometry changes → cliff edge in parameter space
  • Unreasonable flux spectrum → check materials and cross-section libraries
  • Large statistical uncertainties → insufficient sampling in important regions
  • Monte Carlo gives the exact solution to your input — if results look wrong, check the model first

Related Topics

❮ Probability TheoryGeometry ❯