Part 6, Chapter 2

Inertial Confinement

Ablation-driven implosion, Rayleigh-Taylor instability, and laser fusion

2.1 Ablation-Driven Implosion

In inertial confinement fusion (ICF), a small spherical capsule containing deuterium-tritium fuel is compressed to extreme densities using intense laser or X-ray radiation. The outer surface of the capsule ablates outward, and by Newton's third law the remaining shell is driven inward (the "rocket effect"). The implosion velocity typically reaches 300-400 km/s, compressing the fuel to densities of 300-1000 g/cm^3.

The key ignition condition for ICF is expressed in terms of the areal density rho R of the compressed fuel:

$$\rho R > 0.3\;\text{g/cm}^2 \quad \text{(ignition threshold)}$$

This condition ensures that the 3.5 MeV alpha particles produced in D-T fusion reactions deposit their energy within the fuel before escaping. The alpha particle range in compressed DT is approximately lambda_alpha = 0.3 g/cm^2 / rho, so rho R greater than 0.3 g/cm^2 means the alphas are self-absorbed. The ablation pressure driving the implosion scales as:

$$P_{abl} \approx 12\;\text{Mbar}\left(\frac{I}{10^{14}\;\text{W/cm}^2}\right)^{2/3}\left(\frac{\lambda}{0.35\;\mu\text{m}}\right)^{-2/3}$$

where I is the laser intensity and lambda is the laser wavelength. Shorter wavelengths couple energy more efficiently because the critical density n_c = epsilon_0 m_e omega^2 / e^2 increases as lambda^-2, allowing the laser to penetrate closer to the ablation surface.

2.2 Rayleigh-Taylor Instability in ICF

The Rayleigh-Taylor (RT) instability is the primary threat to symmetric implosion. It occurs whenever a heavy fluid is accelerated by a lighter fluid -- in ICF, the dense shell is pushed by the low-density ablation plasma. The classical RT growth rate is:

$$\gamma_{RT} = \sqrt{k g \mathcal{A}}$$

where k is the perturbation wavenumber, g is the acceleration, and the Atwood number is:

$$\mathcal{A} = \frac{\rho_h - \rho_l}{\rho_h + \rho_l}$$

In ICF, the Atwood number is close to 1 because the density contrast is extreme. However, mass ablation stabilizes short-wavelength modes. The ablative RT growth rate is modified to:

$$\gamma_{abl} = \alpha\sqrt{kg} - \beta k v_a$$

where v_a is the ablation velocity and alpha approximately 0.9, beta approximately 1.4-3.1 are constants. The second term shows that ablation flow carries away the perturbation peaks, providing a stabilizing "fire-polishing" effect. Modes with wavelength shorter than a critical value lambda_c = 4 pi^2 beta^2 v_a^2 / (alpha^2 g) are completely stabilized.

During deceleration of the compressed shell, a second RT instability phase occurs at the inner surface (the hot-spot boundary), where the decelerating dense shell pushes against the lighter hot-spot plasma. This deceleration-phase RT determines the integrity of the hot spot and sets stringent requirements on implosion symmetry.

2.3 Hot-Spot Self-Heating and Ignition

The conventional ICF ignition scheme creates a central hot spot surrounded by a dense cold fuel layer. The hot spot must satisfy a more stringent rho R condition for self-heating to propagate a burn wave into the surrounding cold fuel:

$$(\rho R)_{hs} T_{hs} > 1\;\text{g/cm}^2 \cdot \text{keV}$$

This is the ICF analog of the Lawson criterion. For a hot spot at T approximately 5 keV and rho R approximately 0.3 g/cm^2, the product is 1.5 g keV/cm^2, which is marginally above ignition threshold. The total fuel gain G relates the fusion energy output to the laser energy input:

$$G = \eta_{abs}\,\eta_{hydro}\,\frac{E_{fus}}{E_{laser}} = \eta_{abs}\,\eta_{hydro}\,\Phi\left(\frac{\rho R}{6\;\text{g/cm}^2}\right)$$

where eta_abs is the absorption efficiency, eta_hydro is the hydrodynamic efficiency, and Phi is the burn fraction Phi = rho R / (rho R + 6 g/cm^2) for D-T. The National Ignition Facility (NIF) achieved scientific energy breakeven in December 2022, producing 3.15 MJ of fusion energy from 2.05 MJ of laser input -- a historic milestone for ICF.

Two drive geometries exist: direct drive (laser beams directly illuminate the capsule) and indirect drive (lasers heat a hohlraum that produces thermal X-rays to compress the capsule). NIF uses indirect drive for better symmetry, while direct drive offers higher coupling efficiency but requires more uniform beam illumination.

Interactive Simulations

ICF Gain Curves and RT Growth Rates

Python

script.py113 lines

import numpy as np

# ICF Target Gain Curves and Rayleigh-Taylor Analysis
print("Inertial Confinement Fusion: Target Gain and RT Instability")
print("=" * 62)

# 1) Target gain vs laser energy
# Empirical scaling: G ~ eta * phi * E_fus/E_laser
# For indirect drive: eta_coupling ~ 0.08-0.15
# For direct drive: eta_coupling ~ 0.05-0.10
# Burn fraction: phi = rhoR / (rhoR + 6 g/cm^2)
# Implosion velocity scales as: v_imp ~ (I * E_laser)^0.3

E_laser = np.logspace(-1, 1, 50)  # MJ

# Indirect drive scaling (NIF-like)
# Approximate empirical scaling based on hydrodynamic simulations
def target_gain_indirect(E_MJ):
    # rhoR scales with implosion velocity and convergence
    eta_abs = 0.85   # hohlraum absorption
    eta_x = 0.15     # X-ray to kinetic
    eta_hydro = eta_abs * eta_x
    # v_imp ~ 350 km/s for ~1 MJ
    v_imp = 350e3 * (E_MJ / 1.0)**0.12  # weak scaling
    # rhoR ~ rho0 * (v_imp / v0)^2 * R0, simplified scaling
    rhoR = 1.2 * (E_MJ / 1.0)**0.33  # g/cm^2
    phi = rhoR / (rhoR + 6.0)
    E_fus = eta_hydro * E_MJ * 1e6 * phi * 338  # 338 MJ/kg for DT
    G = E_fus / (E_MJ * 1e6)
    # Clamp unrealistic values
    return min(G, 200), rhoR, phi

# Direct drive scaling
def target_gain_direct(E_MJ):
    eta_hydro = 0.08
    rhoR = 1.5 * (E_MJ / 1.0)**0.33
    phi = rhoR / (rhoR + 6.0)
    E_fus = eta_hydro * E_MJ * 1e6 * phi * 338
    G = E_fus / (E_MJ * 1e6)
    return min(G, 200), rhoR, phi

print()
print("Target Gain vs Laser Energy:")
print(f"{'E_laser (MJ)':>13} {'G (indirect)':>13} {'G (direct)':>12} {'rhoR (g/cm2)':>13} {'Burn frac':>10}")
print("-" * 63)

E_table = [0.1, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 10.0]
for E in E_table:
    Gi, rhoR_i, phi_i = target_gain_indirect(E)
    Gd, rhoR_d, phi_d = target_gain_direct(E)
    print(f"{E:13.1f} {Gi:13.1f} {Gd:12.1f} {rhoR_i:13.2f} {phi_i:10.3f}")

print()
print("NIF milestone (Dec 2022): E_laser = 2.05 MJ, E_fusion = 3.15 MJ")
print(f"  Measured gain G = {3.15/2.05:.2f}")

# 2) Rayleigh-Taylor growth rates
print()
print("=" * 62)
print("Rayleigh-Taylor Growth Rates During Implosion")
print("=" * 62)

# Parameters
g_accel = 1e14  # m/s^2 (typical shell acceleration)
v_abl = 3e3     # m/s (ablation velocity)
alpha_RT = 0.9
beta_RT = 1.4

# Mode numbers (spherical harmonics l)
R_shell = 1e-3  # 1 mm initial shell radius
l_modes = np.array([2, 4, 6, 10, 20, 50, 100, 200, 500])
k_modes = l_modes / R_shell  # wavenumber

# Classical RT
gamma_classical = np.sqrt(k_modes * g_accel)  # Atwood ~ 1

# Ablative RT
gamma_ablative = alpha_RT * np.sqrt(k_modes * g_accel) - beta_RT * k_modes * v_abl
gamma_ablative = np.maximum(gamma_ablative, 0)

# Implosion time ~ R/v_imp
v_imp = 3.5e5  # 350 km/s
t_imp = R_shell / v_imp  # ~3 ns

print(f"Shell radius: {R_shell*1e3:.1f} mm")
print(f"Acceleration: {g_accel:.1e} m/s^2")
print(f"Ablation velocity: {v_abl:.0f} m/s")
print(f"Implosion velocity: {v_imp/1e3:.0f} km/s")
print(f"Implosion time: {t_imp*1e9:.1f} ns")
print()

print(f"{'Mode l':>8} {'gamma_cl (1/ns)':>16} {'gamma_abl (1/ns)':>17} {'e-folds (cl)':>13} {'e-folds (abl)':>14}")
print("-" * 70)

for i, l in enumerate(l_modes):
    gc = gamma_classical[i] * 1e-9  # convert to 1/ns
    ga = gamma_ablative[i] * 1e-9
    ef_c = gamma_classical[i] * t_imp
    ef_a = gamma_ablative[i] * t_imp
    print(f"{l:8d} {gc:16.2f} {ga:17.2f} {ef_c:13.1f} {ef_a:14.1f}")

# Critical mode number (cutoff)
l_cutoff = alpha_RT**2 * g_accel * R_shell / (beta_RT**2 * v_abl**2)
print(f"\nCritical cutoff mode: l_c = {l_cutoff:.0f}")
print(f"Modes with l > {l_cutoff:.0f} are stabilized by ablation")

# Convergence ratio
CR = 30  # typical
print(f"\nConvergence ratio CR = {CR}")
print(f"Final hot-spot radius: {R_shell*1e3/CR*1e3:.1f} um")
print(f"If initial perturbation delta_0/R = 1%, final: {0.01 * np.exp(5):.1f}x (for 5 e-folds)")
print("=> Perturbation amplitudes must be kept below ~0.1% for modes l=10-50")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← Previous Next →