Skip to content

Radioactive Waste

In fusion reactors, high-energy neutrons generated by the DT reaction activate the reactor structural materials. However, the radioactive waste produced from fusion reactors has fundamentally different characteristics from that of fission reactors. The ability to significantly reduce environmental impact through appropriate material selection and waste management is one of the key advantages of fusion energy.

This chapter provides a quantitative and systematic explanation of the physical and chemical characteristics unique to fusion reactor radioactive waste, activation mechanisms, development status of reduced-activation materials, waste management strategies, regulatory frameworks, and international initiatives.

Differences Between Fusion and Fission Waste

Section titled “Differences Between Fusion and Fission Waste”

The nature of radioactive waste generated from fusion and fission reactors is fundamentally different. This difference originates from the essential differences in their respective nuclear reaction mechanisms.

In fission reactors, the following radioactive materials are produced by the fission reaction of uranium or plutonium:

  1. Fission products (cesium-137, strontium-90, etc.)
  2. Transuranic elements (plutonium, americium, curium, etc.)

Some of these transuranic elements (actinides) have half-lives ranging from tens of thousands to hundreds of thousands of years, requiring long-term isolation through geological disposal. High-level radioactive waste from spent fuel reprocessing requires management for nearly one million years.

The radioactivity of fission products AFP(t)A_{\text{FP}}(t) is composed of contributions from numerous nuclides:

AFP(t)=iAi,0exp(ln2T1/2,it)A_{\text{FP}}(t) = \sum_i A_{i,0} \exp\left(-\frac{\ln 2}{T_{1/2,i}} t\right)

where Ai,0A_{i,0} is the initial radioactivity of nuclide ii, and T1/2,iT_{1/2,i} is its half-life. Due to the presence of long-lived nuclides (e.g., 239^{239}Pu with a half-life of 24,100 years), the overall radioactivity decay is extremely slow.

The production of transuranic elements proceeds through a chain of neutron capture reactions on uranium fuel:

238U(n,γ)239Uβ239Npβ239Pu^{238}\text{U} \xrightarrow{(n,\gamma)} ^{239}\text{U} \xrightarrow{\beta^-} ^{239}\text{Np} \xrightarrow{\beta^-} ^{239}\text{Pu}

Furthermore, additional neutron capture on 239^{239}Pu leads to transitions to 240^{240}Pu, 241^{241}Pu, and then to americium and curium. These reactions are unique to fission reactors and do not occur in fusion reactors.

In fusion reactors, fission reactions do not occur, so high-level radioactive waste is not generated. All radioactive waste produced originates from the activation of structural materials by neutron irradiation and is classified as low-level radioactive waste.

A characteristic feature of fusion reactor activation products is that the half-lives of the generated radioactive nuclides are relatively short. This can be controlled through material selection, and by using reduced-activation materials, the radioactivity decays to natural uranium ore levels within approximately 100 years after the end of operation.

The fusion reaction itself is:

2D+3T4He(3.5MeV)+n(14.1MeV)^2\text{D} + ^3\text{T} \rightarrow ^4\text{He} (3.5 \, \text{MeV}) + n (14.1 \, \text{MeV})

and all products are stable nuclides or non-radioactive particles. Radioactive waste is produced secondarily when the 14.1 MeV neutrons interact with structural materials.

It is important to understand from a physics perspective why transuranic elements are not produced in fusion reactors.

In fission reactors, elements heavier than uranium (atomic number 92) are produced through successive neutron capture. This is because uranium fuel, a “heavy seed,” is present.

On the other hand, fusion reactor fuels are deuterium (atomic number 1) and tritium (atomic number 1), and it is physically impossible to produce transuranic elements from these light elements. Fusion reactor structural materials also consist of elements lighter than uranium, such as iron (atomic number 26) and tungsten (atomic number 74).

The increase in mass number due to neutron capture is:

ZAX+nZA+1X+γ^A_Z\text{X} + n \rightarrow ^{A+1}_Z\text{X} + \gamma

and the atomic number Z does not change. Therefore, elements heavier than uranium cannot be produced from structural materials.

The main differences are shown below:

ItemFission ReactorFusion Reactor
High-level wasteYes (spent fuel)None
Long-lived nuclidesActinides (half-life tens of thousands of years or more)None
Management periodTens of thousands to hundreds of thousands of yearsApproximately 100 years
Disposal methodGeological disposal requiredNear-surface disposal possible
Toxicity decayVery slowDecays to one millionth in 100 years

Fusion reactor radioactive waste is more than two orders of magnitude less toxic than fission reactor waste even immediately after shutdown, and after 100 years it decreases to one millionth.

A quantitative comparison of radioactivity over time for a 1 GWe power plant:

Time ElapsedFission Reactor (relative value)Fusion Reactor (relative value)
Immediately after shutdown1.00.01
1 year later0.10.001
10 years later0.0510410^{-4}
100 years later0.0110610^{-6}
1000 years later0.00510810^{-8}

This difference is due to the fact that transuranic elements are not produced in fusion reactors.

Mathematical Description of Waste Characteristics

Section titled “Mathematical Description of Waste Characteristics”

The radioactivity decay characteristics of both are determined by the half-life distribution of dominant nuclides. The long-term radioactivity of fission reactor waste is described by the actinide series:

Afission(t)APuexp(tτPu)+AAmexp(tτAm)+A_{\text{fission}}(t) \approx A_{\text{Pu}} \exp\left(-\frac{t}{\tau_{\text{Pu}}}\right) + A_{\text{Am}} \exp\left(-\frac{t}{\tau_{\text{Am}}}\right) + \cdots

where τ=T1/2/ln2\tau = T_{1/2} / \ln 2 is the mean lifetime. The τ\tau of 239^{239}Pu is approximately 35,000 years, and that of 241^{241}Am is approximately 623 years, which dominate long-term radioactivity.

For fusion reactor waste:

Afusion(t)AMn-54et/τMn+AFe-55et/τFe+ACo-60et/τCo+A_{\text{fusion}}(t) \approx A_{\text{Mn-54}} e^{-t/\tau_{\text{Mn}}} + A_{\text{Fe-55}} e^{-t/\tau_{\text{Fe}}} + A_{\text{Co-60}} e^{-t/\tau_{\text{Co}}} + \cdots

Even Co-60, the longest-lived, has τ7.6\tau \approx 7.6 years, and essentially all will have decayed after 100 years.

In DT fusion reactions, 14.1 MeV high-energy neutrons are generated. These neutrons interact with reactor structural materials, causing material activation.

The neutron generation rate in DT fusion reactions is directly related to fusion power:

Sn=PfusionEfusion=Pfusion17.6MeVS_n = \frac{P_{\text{fusion}}}{E_{\text{fusion}}} = \frac{P_{\text{fusion}}}{17.6 \, \text{MeV}}

For 1 GW of fusion power:

Sn=109W17.6×1.602×1013J3.5×1020n/sS_n = \frac{10^9 \, \text{W}}{17.6 \times 1.602 \times 10^{-13} \, \text{J}} \approx 3.5 \times 10^{20} \, \text{n/s}

This enormous neutron flux continuously irradiates the structural materials.

The neutron flux ϕ(E)\phi(E) in fusion reactors is expressed as a function of energy EE. Neutrons from the DT fusion neutron source have an initial energy of 14.1 MeV, but through moderation processes in materials, they form a broad energy spectrum.

Typical values of neutron flux at the first wall are:

ϕwall10141015n/cm2/s\phi_{\text{wall}} \approx 10^{14} - 10^{15} \, \text{n/cm}^2\text{/s}

This is comparable to the flux near the core fuel in fission reactors, but the average energy is significantly higher.

The neutron spectrum ϕ(E)\phi(E) is obtained as a solution of the transport equation with the fusion source term S(E)S(E):

Ωϕ(r,E,Ω)+Σt(E)ϕ=0dE4πdΩΣs(EE,ΩΩ)ϕ+S(r,E)\Omega \cdot \nabla \phi(\mathbf{r}, E, \Omega) + \Sigma_t(E) \phi = \int_0^\infty dE' \int_{4\pi} d\Omega' \, \Sigma_s(E' \to E, \Omega' \to \Omega) \phi + S(\mathbf{r}, E)

where Σt\Sigma_t is the total cross section and Σs\Sigma_s is the scattering cross section.

The spectrum with a 14 MeV peak unique to fusion reactors can be approximately described by group constants considering energy dependence:

ϕg=EgEg1ϕ(E)dE\phi_g = \int_{E_g}^{E_{g-1}} \phi(E) dE

Generally, calculations are performed by dividing into tens to hundreds of energy groups. Libraries such as VITAMIN-J (175 groups) and FENDL (211 groups) are used.

The processes by which neutrons interact with atomic nuclei in materials are mainly classified into the following three types:

  1. Elastic scattering: Neutron and nucleus collide and exchange kinetic energy
  2. Inelastic scattering: Nucleus becomes excited and emits gamma rays
  3. Nuclear reactions: Neutron is absorbed by nucleus and converted to another nuclide

In elastic scattering, the energy loss of neutrons in collision with a nucleus of mass number AA is:

ΔEmax=En4A(1+A)2\Delta E_{\text{max}} = E_n \cdot \frac{4A}{(1+A)^2}

This shows that lighter nuclides can moderate neutrons more efficiently.

The range of neutron energy EE' after scattering is:

αEEE\alpha E \leq E' \leq E

where α=(A1A+1)2\alpha = \left(\frac{A-1}{A+1}\right)^2.

The average logarithmic energy decrement (lethargy) is:

ξ=1+αlnα1α2A+2/3\xi = 1 + \frac{\alpha \ln \alpha}{1 - \alpha} \approx \frac{2}{A + 2/3}

For iron (A=56A = 56), ξ0.035\xi \approx 0.035, and for hydrogen (A=1A = 1), ξ=1\xi = 1.

When fast neutrons collide with atomic nuclei in materials, atoms are knocked out of their normal lattice positions. This phenomenon is called “displacement damage.”

The degree of displacement damage is expressed in dpa (Displacement per Atom). This indicates how many times an atom is displaced. dpa is calculated by the following equation:

dpa=0tdt0dEϕ(E,t)σd(E)\text{dpa} = \int_0^t dt' \int_0^\infty dE \, \phi(E, t') \sigma_d(E)

where σd(E)\sigma_d(E) is the displacement cross section for neutrons of energy EE.

The displacement cross section is evaluated by the Norgett-Robinson-Torrens (NRT) model:

σd(E)=EdTmaxσ(E,T)ν(T)dT\sigma_d(E) = \int_{E_d}^{T_{\text{max}}} \sigma(E, T) \nu(T) dT

where EdE_d is the threshold displacement energy of lattice atoms (approximately 40 eV for iron), TT is the recoil atom energy, and ν(T)\nu(T) is the number of displacements produced by one recoil atom:

ν(T)=0.8Tdam2Ed\nu(T) = \frac{0.8 T_{\text{dam}}}{2 E_d}

TdamT_{\text{dam}} is the damage energy after subtracting electronic excitation losses.

The separation of electronic stopping power and displacement stopping power is described by the Lindhard model:

Tdam=T1+kLg(ϵ)T_{\text{dam}} = \frac{T}{1 + k_L g(\epsilon)}

where kLk_L is the Lindhard constant and g(ϵ)g(\epsilon) is a function of reduced energy ϵ\epsilon.

Fusion reactor structural materials are expected to receive damage of about 150-200 dpa during operation. This is about 10 times that of ITER, at prototype reactor levels.

Displacement damage causes the following phenomena:

  • Irradiation embrittlement: Material becomes brittle
  • Irradiation creep: Deformation progresses under constant stress
  • Swelling: Volume expansion due to vacancy accumulation

The volume change ΔV/V\Delta V/V due to swelling is expressed as a function of vacancy concentration CvC_v and void radius rr:

ΔVV=4π3Nvr3\frac{\Delta V}{V} = \frac{4\pi}{3} N_v r^3

where NvN_v is the void number density per unit volume.

The time evolution of swelling is described by the diffusion equations for vacancies and interstitials:

Cvt=Dv2Cv+GvRivCiCvkv2DvCv\frac{\partial C_v}{\partial t} = D_v \nabla^2 C_v + G_v - R_{iv} C_i C_v - k_v^2 D_v C_v Cit=Di2Ci+GiRivCiCvki2DiCi\frac{\partial C_i}{\partial t} = D_i \nabla^2 C_i + G_i - R_{iv} C_i C_v - k_i^2 D_i C_i

where GG is the production rate, RivR_{iv} is the recombination coefficient, and k2k^2 is the sink strength.

When neutrons are absorbed by atomic nuclei, transmutation occurs, changing them to different elements. Main reactions include:

  • (n, p) reaction: Absorbs neutron and emits proton
  • (n, α) reaction: Absorbs neutron and emits alpha particle (helium)
  • (n, 2n) reaction: Absorbs neutron and emits two neutrons
  • (n, γ) reaction: Gamma ray emission by neutron capture

Gas production rate by transmutation is an important parameter unique to fusion reactors. The helium production rate in iron by 14 MeV neutrons:

n˙He=ϕσ(n,α)N\dot{n}_{\text{He}} = \phi \cdot \sigma_{(n,\alpha)} \cdot N

where σ(n,α)\sigma_{(n,\alpha)} is the (n, α) reaction cross section and NN is the atomic number density.

Typical gas production rates in reduced-activation ferritic steel:

GasProduction Rate (appm/dpa)
Helium10 - 15
Hydrogen40 - 50

Here, appm is atomic parts per million (number of produced atoms per million atoms).

Helium and hydrogen produced by these reactions accumulate in materials and cause degradation of mechanical properties. Helium in particular accumulates at grain boundaries, forming voids and causing embrittlement at high temperatures (helium embrittlement).

The mechanism of helium embrittlement is described by bubble growth at grain boundaries and grain boundary embrittlement:

rbubble=(3nHekBT4πPeq)1/3r_{\text{bubble}} = \left( \frac{3 n_{\text{He}} k_B T}{4 \pi P_{\text{eq}}} \right)^{1/3}

where nHen_{\text{He}} is the number of helium atoms in the bubble and PeqP_{\text{eq}} is the equilibrium pressure.

Among the nuclides produced by transmutation, some are radioactive. When these radioactive nuclides emit beta and gamma rays, the material becomes radioactive. This phenomenon is called “activation,” and the resulting radioactivity is called “induced radioactivity.”

The time evolution of the atomic number density NiN_i of a nuclide ii is described by the Bateman equations:

dNidt=jϕσjiNj+kλkfkiNk(ϕσi+λi)Ni\frac{dN_i}{dt} = \sum_j \phi \sigma_{j \to i} N_j + \sum_k \lambda_k f_{k \to i} N_k - \left( \phi \sigma_i + \lambda_i \right) N_i

where:

  • σji\sigma_{j \to i}: Cross section for conversion from nuclide jj to ii
  • λk\lambda_k: Decay constant of nuclide kk (λ=ln2/T1/2\lambda = \ln 2 / T_{1/2})
  • fkif_{k \to i}: Decay branching ratio from nuclide kk to ii
  • σi\sigma_i: Absorption cross section of nuclide ii

By solving this system of coupled differential equations, the concentration of each nuclide at any time can be obtained.

In actual calculations, thousands to tens of thousands of nuclides must be handled simultaneously. Written in matrix form:

dNdt=AN\frac{d\mathbf{N}}{dt} = \mathbf{A} \mathbf{N}

where A\mathbf{A} is the transition matrix. The solution is:

N(t)=exp(At)N(0)\mathbf{N}(t) = \exp(\mathbf{A} t) \mathbf{N}(0)

Various numerical methods are used to calculate the matrix exponential.

In actual activation, there are nuclides produced through multiple pathways from parent nuclides. For example, the production pathway of Co-60:

59Co(n,γ)60Co^{59}\text{Co} \xrightarrow{(n,\gamma)} ^{60}\text{Co}

Also, the pathway from Ni-58:

58Ni(n,p)58Coβ+58Fe^{58}\text{Ni} \xrightarrow{(n,p)} ^{58}\text{Co} \xrightarrow{\beta^+} ^{58}\text{Fe}

Analysis considering all these competing pathways is necessary.

Under constant neutron flux with long-term irradiation, radioactivity approaches a saturation value. Considering simple one-stage production and decay:

Asat=ϕσNtargetA_{\text{sat}} = \phi \sigma N_{\text{target}}

The radioactivity at irradiation time tt is:

A(t)=Asat(1eλt)A(t) = A_{\text{sat}} \left( 1 - e^{-\lambda t} \right)

When irradiation time exceeds 5 times the half-life, radioactivity reaches more than 97% of the saturation value.

More generally, considering production and destruction:

dNdt=ϕσprodNparent(λ+ϕσabs)N\frac{dN}{dt} = \phi \sigma_{\text{prod}} N_{\text{parent}} - (\lambda + \phi \sigma_{\text{abs}}) N

The steady-state solution is:

Nsat=ϕσprodNparentλ+ϕσabsN_{\text{sat}} = \frac{\phi \sigma_{\text{prod}} N_{\text{parent}}}{\lambda + \phi \sigma_{\text{abs}}}

In the high-flux environment of fusion reactors, the absorption term ϕσabs\phi \sigma_{\text{abs}} may not be negligible.

After reactor shutdown (end of irradiation), radioactivity decays. For a single nuclide:

A(tc)=A0exp(λtc)=A0exp(ln2tcT1/2)A(t_c) = A_0 \exp\left( -\lambda t_c \right) = A_0 \exp\left( -\frac{\ln 2 \cdot t_c}{T_{1/2}} \right)

where tct_c is the cooling time and A0A_0 is the radioactivity at the end of irradiation.

In actual activated materials, numerous nuclides coexist, so the total radioactivity is:

Atotal(tc)=iAi,0exp(λitc)A_{\text{total}}(t_c) = \sum_i A_{i,0} \exp\left( -\lambda_i t_c \right)

This shows a two-stage behavior where decay is rapid in the initial stage dominated by short half-life nuclides, and slow in the later stage when long half-life nuclides remain.

This behavior can be understood as a superposition of multiple components with different half-lives:

A(t)=Ashortet/τshort+Amediumet/τmedium+Alonget/τlongA(t) = A_{\text{short}} e^{-t/\tau_{\text{short}}} + A_{\text{medium}} e^{-t/\tau_{\text{medium}}} + A_{\text{long}} e^{-t/\tau_{\text{long}}}

The relative contribution of each component changes with time, with short-lived components dominating initially and long-lived components dominating later.

Radioactive decay generates thermal energy. The decay heat P(t)P(t) is:

P(t)=iλiNi(t)EˉiP(t) = \sum_i \lambda_i N_i(t) \bar{E}_i

where Eˉi\bar{E}_i is the average energy release per decay of nuclide ii.

Decay heat is calculated as the sum of beta ray, gamma ray, and alpha ray energies:

Eˉi=Eˉβ,i+Eˉγ,i+Eˉα,i\bar{E}_i = \bar{E}_{\beta,i} + \bar{E}_{\gamma,i} + \bar{E}_{\alpha,i}

In beta decay, since some energy is carried away by neutrinos:

Eˉβ=13Eβ,max\bar{E}_{\beta} = \frac{1}{3} E_{\beta,\text{max}}

is approximated (assuming Fermi distribution).

For ITER, the decay heat immediately after shutdown is about 11 MW, but it decreases rapidly with time:

Cooling TimeDecay Heat
Immediately after shutdown11 MW
1 hour later2 MW
1 day later0.6 MW
1 week later0.2 MW
1 month later0.1 MW

Fusion reactor decay heat is very small compared to fission reactors:

  • Does not require emergency cooling systems
  • Natural circulation cooling is sufficient
  • Temperature rise is limited due to large heat capacity of large structures such as vacuum vessel

The decay heat density q(t)q(t) is:

q(t)=P(t)Vq(t) = \frac{P(t)}{V}

The temperature rise in the vacuum vessel (volume V1000V \approx 1000 m³, mass M5000M \approx 5000 tons) is:

ΔT=PtMcp\Delta T = \frac{P \cdot t}{M \cdot c_p}

Using the specific heat of iron cp500c_p \approx 500 J/(kg·K), for 1 hour in adiabatic conditions:

ΔT2×106×36005×106×5002.9K\Delta T \approx \frac{2 \times 10^6 \times 3600}{5 \times 10^6 \times 500} \approx 2.9 \, \text{K}

This level of temperature rise does not affect structural integrity.

The dose rate D˙\dot{D} from activated material is, in point source approximation at distance rr:

D˙(r)=AΓr2B(E,r)eμr\dot{D}(r) = \frac{A \cdot \Gamma}{r^2} \cdot B(E, r) \cdot e^{-\mu r}

where Γ\Gamma is the gamma-ray constant, BB is the buildup factor, and μ\mu is the attenuation coefficient.

The buildup factor represents the contribution of scattered radiation, and in Taylor approximation:

B(μr)=A1eα1μr+(1A1)eα2μrB(\mu r) = A_1 e^{-\alpha_1 \mu r} + (1 - A_1) e^{-\alpha_2 \mu r}

where A1A_1, α1\alpha_1, α2\alpha_2 are parameters that depend on material and energy.

In actual evaluations, Monte Carlo calculation codes (MCNP, Serpent, etc.) are used to calculate dose distributions in complex geometries.

We explain in detail the major radioactive nuclides produced by activation of fusion reactor materials.

Nuclides produced from iron, chromium, and tungsten, which are the main constituent elements of reduced-activation ferritic steels (such as F82H):

  • Half-life: 2.737 years
  • Decay mode: Electron capture (EC)
  • Production reaction: 56^{56}Fe(n, 2n)55^{55}Fe
  • Threshold energy: 11.2 MeV
  • Characteristics: Emits only low-energy X-rays (5.9 keV), small contribution to external exposure

The production cross section is about 460 mb (millibarns) for 14 MeV neutrons, making it one of the dominant activated nuclides in the fusion reactor environment.

The cross section for threshold reactions depends strongly on energy:

σ(E)={0(E<Eth)σ0(1EthE)n(EEth)\sigma(E) = \begin{cases} 0 & (E < E_{\text{th}}) \\ \sigma_0 \left(1 - \frac{E_{\text{th}}}{E}\right)^n & (E \geq E_{\text{th}}) \end{cases}

where nn is a reaction-dependent parameter, typically n12n \approx 1-2 for (n, 2n) reactions.

  • Half-life: 312.2 days
  • Decay mode: Electron capture
  • Gamma-ray energy: 835 keV
  • Production reactions: 54^{54}Fe(n, p)54^{54}Mn, 55^{55}Mn(n, 2n)54^{54}Mn
  • Characteristics: Emits relatively high-energy gamma rays, important for dose evaluation

Mn-54 is the nuclide that dominates dose rate for several years after shutdown and is important in maintenance work planning.

The contribution to dose is:

D˙Mn-54=AMn-54ΓMn-54=AMn-541.18×103(mSv/h)/(MBq at 1 m)\dot{D}_{\text{Mn-54}} = A_{\text{Mn-54}} \cdot \Gamma_{\text{Mn-54}} = A_{\text{Mn-54}} \cdot 1.18 \times 10^{-3} \, \text{(mSv/h)/(MBq at 1 m)}
  • Half-life: 5.27 years
  • Decay mode: Beta decay
  • Gamma-ray energies: 1.17 MeV, 1.33 MeV
  • Production reaction: 59^{59}Co(n, γ)60^{60}Co
  • Characteristics: High-energy gamma rays, large contribution to dose rate

Cobalt exists as an impurity in reduced-activation materials, and its concentration management is extremely important. By keeping impurity concentration below 10 ppm, Co-60 production can be minimized.

The relationship between Co-60 production and impurity concentration:

ACo-60=ϕσγNCo=ϕσγρMCocCoA_{\text{Co-60}} = \phi \cdot \sigma_{\gamma} \cdot N_{\text{Co}} = \phi \cdot \sigma_{\gamma} \cdot \frac{\rho}{M_{\text{Co}}} \cdot c_{\text{Co}}

where cCoc_{\text{Co}} is the mass concentration of cobalt (ppm).

With impurity concentrations of 10 ppm versus 100 ppm, the dose rate from Co-60 differs by a factor of 10.

  • Half-life: 27.7 days
  • Production reactions: 52^{52}Cr(n, 2n)51^{51}Cr, 50^{50}Cr(n, γ)51^{51}Cr
  • Characteristics: Short half-life, disappears during initial cooling stage
  • Half-life: 121.2 days
  • Production reaction: 182^{182}W(n, 2n)181^{181}W
  • Characteristics: Contributes to initial radioactivity of tungsten alloys
  • Half-life: 75.1 days
  • Production reaction: 184^{184}W(n, γ)185^{185}W
  • Decay mode: Beta decay
  • Half-life: 23.7 hours
  • Production reaction: 186^{186}W(n, γ)187^{187}W
  • Characteristics: Dominant nuclide immediately after shutdown, decays rapidly

The activation characteristics of tungsten vary greatly depending on isotope composition. Isotope composition of natural tungsten:

IsotopeNatural AbundanceMain Activation Reaction
180^{180}W0.12%(n, γ) → 181^{181}W
182^{182}W26.50%(n, 2n) → 181^{181}W
183^{183}W14.31%(n, γ) → 184^{184}W (stable)
184^{184}W30.64%(n, γ) → 185^{185}W
186^{186}W28.43%(n, γ) → 187^{187}W
  • Half-life: 20,300 years
  • Production reaction: 93^{93}Nb(n, γ)94^{94}Nb
  • Characteristics: Extremely long half-life, reason for excluding niobium from reduced-activation materials

The presence of niobium has a significant impact on long-term management of fusion reactor waste. In F82H, niobium concentration is controlled to below 1 ppm.

The reason why Nb-94 is problematic is shown quantitatively. The residual radioactivity after 100 years of cooling is:

ANb-94(100yr)=A0eln2100/203000.997A0A_{\text{Nb-94}}(100 \, \text{yr}) = A_0 \cdot e^{-\ln 2 \cdot 100/20300} \approx 0.997 \cdot A_0

That is, it hardly decays even after 100 years. On the other hand, Mn-54:

AMn-54(100yr)=A0eln2100365/3121035A00A_{\text{Mn-54}}(100 \, \text{yr}) = A_0 \cdot e^{-\ln 2 \cdot 100 \cdot 365/312} \approx 10^{-35} \cdot A_0 \approx 0

It completely disappears.

  • Half-life: 4,000 years
  • Production reaction: 92^{92}Mo(n, γ)93^{93}Mo
  • Characteristics: Problematic in molybdenum-containing materials, reduced-activation steels replace molybdenum with tungsten
  • Half-life: 76,000 years
  • Production reaction: 58^{58}Ni(n, γ)59^{59}Ni
  • Characteristics: Cause of long-term radioactivity in nickel-based alloys
  • Half-life: 100.1 years
  • Production reaction: 62^{62}Ni(n, γ)63^{63}Ni
  • Decay mode: Pure beta (66.9 keV)
  • Characteristics: Contributes to medium-term radioactivity

In nickel-containing materials (316SS, etc.), Ni-59 and Ni-63 dominate long-term radioactivity.

In aluminum alloys being considered as lightweight structural materials, there are special problems:

  • Half-life: 717,000 years
  • Production reaction: 27^{27}Al(n, 2n)26^{26}Al
  • Threshold energy: 13.5 MeV
  • Characteristics: Extremely long half-life, produced by 14 MeV neutrons

This reaction is a problem unique to fusion. Since the threshold is 13.5 MeV, it does not occur in fission reactors, but it becomes a problem in fusion reactors with 14 MeV neutrons.

Calculation of Radioactive Nuclide Production

Section titled “Calculation of Radioactive Nuclide Production”

The atomic number density of nuclide ii at the end of irradiation is, in a simplified model:

Ni=ϕσprodλiNparent(1eλitirr)N_i = \frac{\phi \sigma_{\text{prod}}}{\lambda_i} N_{\text{parent}} \left( 1 - e^{-\lambda_i t_{\text{irr}}} \right)

where σprod\sigma_{\text{prod}} is the production cross section and tirrt_{\text{irr}} is the irradiation time.

In actual calculations, activation calculation codes such as FISPACT-II, ACAB, and ORIGEN are used, considering complete decay chains and cross-section libraries.

Activation Calculation Codes and Evaluation Methods

Section titled “Activation Calculation Codes and Evaluation Methods”

Dedicated calculation codes and nuclear data libraries are used for activation evaluation of fusion reactors.

Activation calculation code developed by UK Atomic Energy Authority (UKAEA), a standard tool for fusion applications.

Features:

  • Matrix exponential solution of Bateman equations
  • Uncertainty propagation analysis capability
  • Sensitivity analysis capability
  • Support for multiple nuclear data libraries

The calculation flow is:

  1. Input neutron spectrum
  2. Calculate group-averaged cross sections
  3. Build transition matrix
  4. Calculate time evolution
  5. Output radioactivity, decay heat, dose rate

Group-averaged cross section:

σˉg=EgEg1σ(E)ϕ(E)dEEgEg1ϕ(E)dE\bar{\sigma}_g = \frac{\int_{E_g}^{E_{g-1}} \sigma(E) \phi(E) dE}{\int_{E_g}^{E_{g-1}} \phi(E) dE}

Code developed by Oak Ridge National Laboratory (ORNL), provided as part of the SCALE package.

Features:

  • Extensive decay data library
  • Consideration of spontaneous fission, (α, n) reactions
  • Neutron source term calculation
  • Processing of multi-stage irradiation history

Activation calculation code developed by Japan Atomic Energy Agency (JAEA).

Features:

  • Compatible with JENDL (Japanese Evaluated Nuclear Data Library)
  • Optimized for Japanese fusion research
  • Used in ITER TBM evaluation

Standard activation analysis system used by EUROfusion.

Components:

  • FISPACT-II (inventory calculation)
  • EAF (European Activation File)
  • TENDL (nuclear data library)

The accuracy of activation calculations depends greatly on the quality of nuclear data.

Activation-specific nuclear data library developed in Europe:

  • 816 target nuclides
  • Approximately 62,000 reactions
  • Optimized for 14 MeV neutrons

Evaluated nuclear data library based on the TALYS nuclear reaction code:

  • Systematic evaluation based on theoretical models
  • Includes uncertainty information
  • Updated annually

Japanese nuclear data library for activation:

  • Added activation reactions based on JENDL-4.0
  • Approximately 770 target nuclides
  • Standard use in Japanese fusion research

FENDL (Fusion Evaluated Nuclear Data Library)

Section titled “FENDL (Fusion Evaluated Nuclear Data Library)”

Library maintained by IAEA for fusion applications:

  • International standard
  • Applicable to both transport and activation calculations
  • Used in ITER design

Activation calculations have various uncertainty sources:

  1. Nuclear data uncertainty

    • Cross section measurement errors
    • Covariance data
  2. Neutron spectrum uncertainty

    • Transport calculation errors
    • Geometric model approximations
  3. Material composition uncertainty

    • Impurity concentration variations
    • Variability between manufacturing lots

Uncertainty propagation is:

σA2=i(Aσi)2σσi2+i(ANi)2σNi2\sigma_A^2 = \sum_i \left( \frac{\partial A}{\partial \sigma_i} \right)^2 \sigma_{\sigma_i}^2 + \sum_i \left( \frac{\partial A}{\partial N_i} \right)^2 \sigma_{N_i}^2

In FISPACT-II, uncertainty analysis by Monte Carlo sampling is possible:

  1. Set probability distributions for input parameters
  2. Multiple sampling calculations
  3. Statistical analysis of outputs

The validity of activation calculations is verified by comparison with experimental data.

Accelerator neutron source of Japan Atomic Energy Agency:

  • 14 MeV neutron beam
  • Sample irradiation and activation measurement
  • Comparison verification with calculations (C/E values)

Representative verification results:

NuclideC/E ValueUncertainty
Mn-540.95±10%
Fe-551.02±15%
Co-601.10±12%

The closer the C/E (Calculation/Experiment) value is to 1.0, the higher the calculation accuracy.

High-intensity neutron irradiation facility under construction:

  • Neutron spectrum equivalent to fusion
  • High irradiation dose of 40 dpa/year
  • Acquisition of material irradiation data

Activation calculation example of F82H steel (1 MW/m² wall load, 5-year irradiation):

Input conditions:

  • Material: F82H (Fe-8Cr-2W-0.2V-0.04Ta)
  • Neutron flux: 5×10145 \times 10^{14} n/cm²/s
  • Neutron spectrum: ITER first wall conditions
  • Irradiation time: 5 years
  • Cooling time: 0-100 years

Calculation results (specific radioactivity):

Cooling TimeSpecific Radioactivity (Bq/kg)Major Nuclides
0 (immediately after shutdown)3×10123 \times 10^{12}Mn-56, W-187
1 day1×10121 \times 10^{12}Mn-54, Fe-55
1 year5×10115 \times 10^{11}Mn-54, Fe-55
10 years1×10111 \times 10^{11}Fe-55, Co-60
100 years1×1081 \times 10^{8}Nb-94, Mo-93

Radioactive waste from fusion reactors is classified based on radioactivity concentration and half-life.

IAEA (International Atomic Energy Agency) classification system:

  1. Exempt Waste (EW)

    • Below clearance levels
    • Can be processed as non-regulated
  2. Very Short Lived Waste (VSLW)

    • Dominated by short half-life nuclides
    • Decays with several years of storage
  3. Very Low Level Waste (VLLW)

    • Near-surface disposal possible
    • Management period: Several decades
  4. Low Level Waste (LLW)

    • Near-surface disposal
    • Management period: Several hundred years
  5. Intermediate Level Waste (ILW)

    • Requires deeper disposal
    • Low heat generation
  6. High Level Waste (HLW)

    • Requires deep geological disposal
    • Not generated in fusion reactors

IAEA classification and disposal method correspondence:

CategoryRadioactivity LevelHeatDisposal Method
EWBelow clearanceNoneGeneral waste
VSLWVery lowNoneDecay storage
VLLWVery lowNoneNear-surface disposal
LLWLowNegligibleNear-surface disposal
ILWIntermediateLowIntermediate depth disposal
HLWHighHighDeep geological disposal

Disposal categories for low-level radioactive waste defined by the Atomic Energy Commission:

CategoryRadioactivity ConcentrationDisposal MethodManagement Period
L3Very low levelTrench disposalSeveral decades
L2Low levelNear-surface pit disposal300 years
L1Relatively highIntermediate depth disposalSeveral thousand years

Activated materials from fusion reactors are classified as L2 or L3 after appropriate cooling periods.

Concentration upper limits for each category (representative nuclides):

NuclideL3 Upper Limit (Bq/t)L2 Upper Limit (Bq/t)L1 Upper Limit (Bq/t)
Co-60101010^{10}101210^{12}101410^{14}
Ni-63101110^{11}101310^{13}101510^{15}
Nb-9410710^{7}10910^{9}101110^{11}

Classification by NRC (Nuclear Regulatory Commission) 10 CFR 61:

  • Lowest level
  • No stabilization required
  • Decays in 100 years
  • Structural stabilization required
  • Decays in 300 years
  • Intrusion prevention design
  • Highest level of LLW
  • Burial at depth of 5 m or more
  • Intrusion prevention barriers
  • Level exceeding Class C
  • Special disposal required
  • Fusion reactor waste does not fall into this category

Class C concentration upper limits:

NuclideConcentration Upper Limit
Ni-63700 Ci/m³
Nb-940.2 Ci/m³
Co-60Unlimited
C-148 Ci/m³

Fusion Reactor-Specific Waste Classification

Section titled “Fusion Reactor-Specific Waste Classification”

Practical classification considering characteristics of fusion reactor waste:

The highest level for near-surface disposal under US NRC classification. Highly activated components of fusion reactors (first wall, divertor, etc.) fall within this level after cooling.

Nuclide concentration upper limits for determination:

NuclideConcentration Upper Limit (Ci/m³)
Ni-63700
Nb-940.2
Co-60Unlimited

Radioactivity level that allows remote-operated recycling. Defined by dose rate:

D˙<10mSv/h at 1 m\dot{D} < 10 \, \text{mSv/h at 1 m}

Low radioactivity materials allowing direct work:

D˙<10μSv/h at 1 m\dot{D} < 10 \, \mu\text{Sv/h at 1 m}

Estimated waste generation from a 1 GWe fusion power plant:

ComponentReplacement CycleWeight (tons/year)Major Nuclides
First wall2-3 years50-100Mn-54, Fe-55
Blanket4-5 years200-400W isotopes, Cr-51
Divertor2 years20-50W-181, W-185
Vacuum vesselPlant lifetime1000-2000Co-60, Ni-63

After 40 years of operation, total waste volume is estimated at about 10,000 tons. This is comparable to a fission reactor of the same output (including spent fuel), but the radioactivity and toxicity are much lower.

Time evolution of waste volume:

Vwaste(t)=componentsMcomponenttcycletV_{\text{waste}}(t) = \sum_{\text{components}} \frac{M_{\text{component}}}{t_{\text{cycle}}} \cdot t

Estimated waste amounts by category after 40 years of operation:

CategoryMass (tons)Volume (m³)Proportion
L1 equivalent2,00040020%
L2 equivalent5,0001,00050%
L3 equivalent2,00040020%
Clearance1,00020010%

Transition between categories is possible with cooling periods.

To reduce radioactive waste from fusion reactors, development of “reduced-activation materials” that do not readily produce long-lived radioactive nuclides under neutron irradiation is being advanced.

The basic principle of reduced-activation material design is to exclude elements that produce long half-life nuclides.

The activation characteristics of elements are determined by the following factors:

  1. Abundance ratios of stable isotopes
  2. Neutron reaction cross sections
  3. Half-lives of produced nuclides
  4. Decay modes (gamma-ray energies)

Definition of Activation Potential (AP):

AP=iσifiAiMPCiT1/2,i\text{AP} = \sum_i \sigma_i \cdot f_i \cdot \frac{A_i}{\text{MPC}_i} \cdot T_{1/2,i}

where fif_i is the production branching ratio and MPC is the Maximum Permissible Concentration.

Elements suitable as reduced-activation materials, when irradiated with 14 MeV neutrons:

  • Do not produce long half-life radioactive nuclides
  • Produced radioactive nuclides decay relatively quickly

As a quantitative evaluation index, Contact Dose Rate (CDR) is used:

CDR(t)=iAi(t)Γi\text{CDR}(t) = \sum_i A_i(t) \cdot \Gamma_i

where Γi\Gamma_i is the gamma-ray constant of nuclide ii.

The target for reduced-activation materials is to achieve after 100 years of cooling:

CDR(100yr)<10μSv/h\text{CDR}(100 \, \text{yr}) < 10 \, \mu\text{Sv/h}

This corresponds to a level allowing direct hands-on work.

Elements can be classified into 3 groups by activation characteristics:

Decay to near clearance levels after 100 years of cooling:

  • C, Si, Ti, V, Cr, Mn, Fe, Ta, W

LLW level with several hundred years of cooling:

  • Cu, Y, Zr

Produce long-lived nuclides, elements that should be excluded:

  • Nb, Mo, Ni, Co, Al, Ag

The following elements are excluded from reduced-activation materials because they produce long-lived radioactive nuclides:

ElementProduced NuclideHalf-lifeProduction Reaction
NbNb-9420,300 years(n, γ)
MoMo-934,000 years(n, γ)
NiNi-5976,000 years(n, γ)
CoCo-605.27 years(n, γ)
AlAl-26717,000 years(n, 2n)

These elements also require management when present as impurities.

Relationship between impurity concentration and long-term radioactivity:

Along(t)cimpurityϕσ(1eλt)A_{\text{long}}(t) \propto c_{\text{impurity}} \cdot \phi \cdot \sigma \cdot (1 - e^{-\lambda t})

Reduced-Activation Ferritic/Martensitic Steel (RAFM Steel)

Section titled “Reduced-Activation Ferritic/Martensitic Steel (RAFM Steel)”

RAFM (Reduced Activation Ferritic/Martensitic) steel is the primary candidate for fusion reactor structural materials.

Development of RAFM steel began in the 1980s:

  • 1985: Start of conceptual design
  • 1990s: Laboratory-scale manufacturing
  • 2000s: Establishment of industrial-scale manufacturing
  • 2010s: Adoption for ITER TBM

F82H, developed by Japan Atomic Energy Research Institute (now National Institutes for Quantum Science and Technology), is a representative reduced-activation ferritic steel.

Basic composition:

ElementContentRole
FeBalanceBase material
Cr7.5-8.5%Corrosion resistance
W1.5-2.5%Solid solution strengthening
V0.15-0.25%Precipitation strengthening
Ta0.02-0.08%Precipitation strengthening
C0.08-0.12%Carbide formation
Mn0.1-0.5%Deoxidation

Impurity management upper limits:

ImpurityUpper LimitReason
Nb< 1 ppmSuppress Nb-94 production
Mo< 50 ppmSuppress Mo-93 production
Ni< 200 ppmSuppress Ni-59/63 production
Co< 10 ppmSuppress Co-60 production
Cu< 100 ppmSuppress irradiation embrittlement

Mechanical properties of F82H:

PropertyValueCondition
Tensile strength550 MPaRoom temperature
Yield strength450 MPaRoom temperature
Elongation20%Room temperature
Impact transition temperature-70°CDBTT
Creep rupture time10⁵ h550°C, 120 MPa

Features of F82H:

  • Replaced molybdenum with tungsten from conventional chrome-molybdenum steel
  • Manufacturing experience of 5-ton scale through IEA (International Energy Agency) international joint research
  • Adopted as structural material for ITER Test Blanket Module
  • Primary candidate material for prototype reactor

The upper limit of operating temperature is about 550°C, used in combination with water-cooled blankets.

RAFM steel developed in Europe with a composition similar to F82H.

Composition: Fe-9Cr-1.1W-0.2V-0.12Ta-0.11C

Main differences:

  • Slightly higher tungsten content (1.1%)
  • Higher vanadium content (0.2%)
  • Reference material for DEMO design

Industrial-scale manufacturing of EUROFER97 was established through cooperation of 15 European countries.

Groups of reduced-activation steels developed in the United States. Various variations are being studied.

ODS Steel (Oxide Dispersion Strengthened Steel)

Section titled “ODS Steel (Oxide Dispersion Strengthened Steel)”

To overcome the operating temperature limit of conventional RAFM steel (550°C), development of oxide dispersion strengthened (ODS) steel is being advanced.

By dispersing fine oxide particles (Y₂O₃, etc.) in steel, high-temperature strength is improved:

σy=σ0+σOrowan\sigma_y = \sigma_0 + \sigma_{\text{Orowan}}

Orowan strengthening:

σOrowan=0.8MGb2π1νln(d/b)λ\sigma_{\text{Orowan}} = \frac{0.8 M G b}{2\pi\sqrt{1-\nu}} \cdot \frac{\ln(d/b)}{\lambda}

where MM is the Taylor factor, GG is the shear modulus, bb is the Burgers vector, dd is the particle diameter, and λ\lambda is the particle spacing.

Composition: Fe-9Cr-2W-0.2V-0.35Y₂O₃

Properties:

  • Operating temperature: Up to 650°C
  • Excellent high-temperature creep strength
  • Stable oxide dispersion under irradiation

Manufacturing process:

  1. Mechanical alloying
  2. Sintering and hot working
  3. Heat treatment

Challenges:

  • High cost
  • Anisotropy
  • Joining technology

Features of V-4Cr-4Ti (vanadium-4% chromium-4% titanium) alloy:

  • Extremely low induced radioactivity
  • Excellent compatibility with liquid lithium
  • High-temperature strength
  • Challenges: Irradiation embrittlement, workability

Comparison of radioactivity decay (100 years after irradiation):

MaterialRelative Radioactivity
316SS (conventional material)1.0
F82H0.01
V-4Cr-4Ti0.001

Main produced nuclides in vanadium alloys:

NuclideHalf-lifeProduction Reaction
Ti-453.1 hTi-46(n, 2n)
V-523.8 minV-51(n, γ)
Cr-5127.7 dCr-52(n, 2n)
Ca-45163 dTi-48(n, α)

Short-lived nuclides are dominant, and radioactivity after 100 years of cooling is very low.

Silicon carbide fiber-reinforced silicon carbide matrix composites:

  • Inherently reduced activation (Si, C are difficult to activate)
  • Ultra-high temperature capability (above 1400°C)
  • High thermal efficiency (> 50%)
  • Challenges: Property changes under irradiation, joining technology

Structure of SiC/SiC composite materials:

  1. SiC fibers (diameter 10-15 μm)
  2. Interface layer (PyC or BN, thickness 0.1-1 μm)
  3. SiC matrix

Manufacturing methods:

  • CVI (Chemical Vapor Infiltration)
  • PIP (Polymer Infiltration and Pyrolysis)
  • MI (Melt Infiltration)

Main nuclides produced:

NuclideHalf-lifeProduction Reaction
C-145,730 yearsN-14(n, p) impurity origin
Si-32172 yearsSi-30(n, γ)Si-31(β)Si-32
H-312.3 yearsLi-6(n, α) impurity origin

C-14 production depends on nitrogen impurities:

AC-14=ϕσ(n,p)NNA_{\text{C-14}} = \phi \cdot \sigma_{(n,p)} \cdot N_{\text{N}}

By controlling nitrogen impurities below 10 ppm, C-14 production can be minimized.

Activation characteristics of tungsten and beryllium used as plasma-facing materials:

  • High melting point, low sputtering rate, ideal for divertor
  • Activation: W-181, W-185, W-187, etc.
  • Long-lived nuclides are not readily produced
  • Low-level waste after 100 years of cooling

Tungsten activation calculation example (divertor conditions):

Cooling TimeSpecific Radioactivity (Bq/kg)Major Nuclides
0101310^{13}W-187
1 day101210^{12}W-185, W-181
1 year101010^{10}W-181, Ta-182
100 years10710^{7}Ta-182, Hf-178m
  • Neutron multiplier, first wall coating material
  • Activation: Be-10 (half-life 1.6 million years)
  • Production amount is small but long half-life
  • Also issues with tritium accumulation

Be-10 production reaction:

9Be+n10Be+γ^9\text{Be} + n \rightarrow ^{10}\text{Be} + \gamma 9Be(n,2n)8Be2α^9\text{Be}(n, 2n)^8\text{Be} \rightarrow 2\alpha

Helium is produced by (n, 2n) reactions, affecting material properties.

Clearance systems and recycling are being considered for volume reduction of radioactive waste.

Clearance is a system that excludes materials with radioactivity concentrations below a certain level (clearance level) from regulation as radioactive waste. Materials below clearance levels can be processed and disposed of as general waste.

Development of clearance systems:

  • 1988: IAEA proposes concept
  • 1996: IAEA Safety Guide RS-G-1.7
  • 2004: Institutionalized in Japan (amendment of Reactor Regulation Act)
  • 2005: EU Directive

Based on IAEA Safety Guide RS-G-1.7, determined under the following conditions:

  • Individual dose: 10 μSv/year or less
  • Collective dose: 1 man·Sv/year or less

Dose evaluation considers various exposure pathways:

  1. External exposure

    • Direct gamma rays
    • Exposure from surface contamination
  2. Internal exposure

    • Inhalation intake
    • Ingestion

The most restrictive pathway is identified through scenario analysis:

  • Metal recycling scenario
  • Landfill disposal scenario
  • Reuse scenario

The clearance level CLC_L is determined from the most restrictive pathway:

CL=DlimitDCFintakefactorC_L = \frac{D_{\text{limit}}}{\text{DCF} \cdot \text{intake} \cdot \text{factor}}

where DCF (Dose Conversion Factor) is the dose conversion factor.

NuclideClearance Level (Bq/g)
Fe-551000
Mn-540.1
Co-600.1
Ni-63100
Nb-940.1

Nuclide-specific clearance levels in Japanese nuclear regulation (Reactor Regulation Act):

NuclideRadioactivity Concentration (Bq/g)
H-3100
C-141
Fe-551000
Co-600.1
Ni-63100

Clearable rate of fusion reactor components after cooling:

Component10 Years Cooling50 Years Cooling100 Years Cooling
Outer blanket30%60%80%
Vacuum vessel40%70%90%
Shield60%85%95%
Support structure80%95%99%

By providing appropriate cooling periods, waste amounts can be significantly reduced.

Calculation of clearance achievement time:

For each nuclide, the clearance achievement time tct_c is obtained:

Ai(tc)=Ai,0exp(λitc)=CL,iA_i(t_c) = A_{i,0} \exp(-\lambda_i t_c) = C_{L,i} tc=1λiln(Ai,0CL,i)=T1/2,iln2ln(Ai,0CL,i)t_c = \frac{1}{\lambda_i} \ln\left(\frac{A_{i,0}}{C_{L,i}}\right) = \frac{T_{1/2,i}}{\ln 2} \ln\left(\frac{A_{i,0}}{C_{L,i}}\right)

Time when all nuclides satisfy the condition:

tclearance=maxi{tc,i}t_{\text{clearance}} = \max_i \{ t_{c,i} \}

Even materials processed as radioactive waste have recycling possibilities.

  1. Unrestricted reuse

    • Below clearance level
    • Reuse as general materials
  2. Restricted reuse

    • Reuse within nuclear facilities
    • Processing into new fusion reactor components
  3. Remote-operated recycling

    • High radioactivity materials
    • Processing in hot cells

Comparison of recycling cost CRC_R and disposal cost CDC_D:

CR=Cseparation+Cprocessing+CfabricationC_R = C_{\text{separation}} + C_{\text{processing}} + C_{\text{fabrication}} CD=Cconditioning+Ctransport+CdisposalC_D = C_{\text{conditioning}} + C_{\text{transport}} + C_{\text{disposal}}

Economic decision criterion:

Recycling advantageousCR+Vmaterial<CD\text{Recycling advantageous} \Leftrightarrow C_R + V_{\text{material}} < C_D

where VmaterialV_{\text{material}} is the market value of recycled materials.

Rare metals such as tungsten and vanadium have high economic value for recycling.

Recycling value by material:

MaterialPrice ($/kg)Annual Generation (tons)Potential Value ($/year)
Tungsten30501.5M
Vanadium20200.4M
Beryllium50052.5M
Steel15000.5M
  • Separation by material type
  • Processing technology under residual radioactivity
  • Quality assurance (management of impurity accumulation)
  • Development of regulatory framework

Recycling process for radioactive materials:

  1. Dismantling and cutting
  2. Decontamination
  3. Separation (material type, radioactivity level)
  4. Melting and refining
  5. Processing into new materials
  6. Quality inspection

Waste management of fusion reactors requires a consistent strategy from design stage through operation to decommissioning.

  1. Selection of reduced-activation materials
  2. Formulation of impurity management specifications
  3. Consideration of ease of dismantling
  4. Design for waste minimization

Waste impact assessment formula in design:

Iwaste=imifi(Ci,tcool)I_{\text{waste}} = \sum_i m_i \cdot f_i(C_i, t_{\text{cool}})

where mim_i is the component mass and fif_i is a waste coefficient depending on cooling time tcoolt_{\text{cool}} and composition CiC_i.

Objective function for design optimization:

mindesign[Cconstruction+Coperation+Cdecommissioning]\min_{design} \left[ C_{\text{construction}} + C_{\text{operation}} + C_{\text{decommissioning}} \right]

With constraints:

  • Safety requirements
  • Performance requirements
  • Waste limitations
  1. Verification of material purity
  2. Obtaining impurity certificates
  3. Ensuring traceability
  4. Storage of manufacturing records

Quality assurance program for impurity management:

  • Analysis certification of raw materials
  • Control of manufacturing processes
  • Inspection of final products
  • Long-term storage of records
  1. Activation monitoring
  2. Management of replacement components
  3. Operation of interim storage facilities
  4. Management of cooling periods

Tracking of activation inventory during operation:

dAdt=ϕ(t)σNλA\frac{dA}{dt} = \phi(t) \sigma N - \lambda A

Calculations considering time-varying neutron flux ϕ(t)\phi(t) are necessary.

  1. Formulation of dismantling plan
  2. Implementation of decontamination work
  3. Separation and classification
  4. Final disposal

Decommissioning cost model:

CD&D=Cprep+Cremoval+Cprocessing+CdisposalC_{\text{D\&D}} = C_{\text{prep}} + C_{\text{removal}} + C_{\text{processing}} + C_{\text{disposal}}

Replacement components and decommissioning waste require interim storage for cooling.

  • Shielding capability
  • Cooling capability (initial stage)
  • Environmental monitoring
  • Safety management system

Shielding thickness design:

D˙(d)=D˙0B(d)eμd<D˙limit\dot{D}(d) = \dot{D}_0 \cdot B(d) \cdot e^{-\mu d} < \dot{D}_{\text{limit}}

For concrete shielding, μ0.1\mu \approx 0.1 cm⁻¹ (1 MeV gamma rays).

The optimal storage period toptt_{\text{opt}} minimizes the sum of storage cost CS(t)C_S(t) and disposal cost CD(t)C_D(t):

topt=argmint[CS(t)+CD(t)]t_{\text{opt}} = \arg\min_t \left[ C_S(t) + C_D(t) \right]

Extending the storage period:

  • Lowers disposal category (reduces disposal cost)
  • Increases storage cost
  • Increases land use cost

Generally, a cooling period of 50-100 years is the economic optimum.

Detailed economic model:

Ctotal(t)=CSt+CD(A(t))C_{\text{total}}(t) = C_S \cdot t + C_D(A(t))

where CD(A)C_D(A) is the disposal cost depending on radioactivity level AA.

Measures to minimize waste generation:

  1. Optimization of material selection

    • Use of reduced-activation materials
    • Thorough impurity management
  2. Design optimization

    • Ingenious shielding arrangement
    • Localization of activated regions
  3. Operational management

    • Uniformization of neutron fluence distribution
    • Efficient use of materials
  4. Dismantling strategy

    • Effective separation
    • Utilization of clearance

The waste volume reduction factor RR is:

R=VdisposedVtotal=1fclearancefrecycleR = \frac{V_{\text{disposed}}}{V_{\text{total}}} = 1 - f_{\text{clearance}} - f_{\text{recycle}}

The target is R<0.3R < 0.3 (70% or more to clearance or recycling).

Dismantling and Decontamination Technology

Section titled “Dismantling and Decontamination Technology”

Decommissioning of fusion reactors requires specialized dismantling and decontamination technologies.

Remote operation technology is essential for dismantling highly activated components.

Main technologies:

  • Master-slave manipulators
  • Program-controlled robots
  • Remote cutting devices

Comparison of cutting methods:

MethodCutting SpeedSecondary WasteApplication Range
Plasma cuttingHighMediumThick plates
Laser cuttingMediumLowPrecision machining
Mechanical cuttingLowLowThin plates
Abrasive water jetMediumMediumComplex shapes

Physics of each cutting method:

Plasma cutting:

P=12ρv3A+ρvΔHP = \frac{1}{2} \rho v^3 A + \rho v \Delta H

where ρ\rho is material density, vv is cutting speed, AA is cutting area, and ΔH\Delta H is melting enthalpy.

Laser cutting:

vcut=αIρ(cpΔT+Lm)v_{\text{cut}} = \frac{\alpha I}{\rho (c_p \Delta T + L_m)}

where II is laser intensity, α\alpha is absorptivity, and LmL_m is latent heat of melting.

Dismantling strategies for large structures such as vacuum vessels:

  1. Segmentation approach

    • Divisible structure from design stage
    • No large on-site cutting required
  2. On-site dismantling approach

    • Use of large cutting equipment
    • Securing of removal routes
  3. One-piece removal

    • Remove entire structure
    • Dismantling at dedicated facility

Surface contamination removal aims to reduce waste volume and improve working environment.

Surface contamination removal using acids and chelating agents:

M+nH++mLnMLm(nmn)++n2H2M + nH^+ + mL^{n-} \rightarrow ML_m^{(n-mn)+} + \frac{n}{2}H_2

Main decontaminants:

  • Acids (nitric acid, sulfuric acid, phosphoric acid)
  • Chelating agents (EDTA, citric acid)
  • Redox agents (potassium permanganate)

Decontamination Factor (DF):

DF=AbeforeAafter\text{DF} = \frac{A_{\text{before}}}{A_{\text{after}}}

Typical DF for chemical decontamination is about 10-100.

Surface layer removal by electrolytic polishing:

MMn++neM \rightarrow M^{n+} + ne^-

According to Faraday’s law, the removal amount is:

m=MItnFm = \frac{M \cdot I \cdot t}{n \cdot F}

where MM is atomic weight, II is current, tt is time, and FF is Faraday constant.

Features:

  • Uniform surface removal
  • Applicable to complex shapes
  • Secondary effluent treatment required

Physical surface removal:

  • Blast treatment
  • Grinding
  • Ultra-high pressure water washing

Relationship between surface removal thickness dd and decontamination effect (when contamination is limited to surface layer):

A(d)=A0exp(dλ)A(d) = A_0 \exp\left(-\frac{d}{\lambda}\right)

where λ\lambda is the contamination penetration depth.

Pressure and removal effect of ultra-high pressure water washing:

m˙removalP0.5Q\dot{m}_{\text{removal}} \propto P^{0.5} \cdot Q

where PP is pressure and QQ is flow rate.

UseTarget Surface Contamination Density
Unrestricted release< 0.4 Bq/cm² (β, γ)
Restricted reuse< 4 Bq/cm² (β, γ)
Within controlled areaBelow regulatory limits

ITER decommissioning will be a model case for fusion reactor decommissioning.

  1. End of operation and transition period (5 years)

    • Tritium removal
    • System decontamination
    • Initial dismantling preparation
  2. Dismantling period (10 years)

    • Removal of highly activated components
    • Remote dismantling work
    • Waste processing and packaging
  3. Site restoration (5 years)

    • Building demolition
    • Environmental restoration
    • Final measurement and release
CategoryGenerated Amount (tons)
High activation (L1 equivalent)4,000
Medium activation (L2 equivalent)8,000
Low activation (L3 equivalent)12,000
Clearance target6,000

Main challenges for ITER decommissioning:

  1. Management of beryllium dust
  2. Processing of tritium-contaminated components
  3. Dismantling of large superconducting coils
  4. Remote removal of highly activated divertor

Final disposal of fusion reactor waste applies methods according to radioactivity levels.

The majority of fusion reactor waste can be disposed of by near-surface disposal.

Applied to very low level waste (L3):

  • Burial in trenches near the surface
  • Shielding by natural soil
  • Management period: About 50 years

Trench design criteria:

  • Depth: 3-5 m
  • Width: 10-20 m
  • Cover thickness: 1-2 m

Dose calculation in safety assessment:

D=CVfreleasefdilutionDCFD = C \cdot V \cdot f_{\text{release}} \cdot f_{\text{dilution}} \cdot \text{DCF}

where CC is the waste package concentration, VV is the volume, freleasef_{\text{release}} is the release factor, and fdilutionf_{\text{dilution}} is the dilution factor.

Applied to low level waste (L2):

  • Storage in concrete pits
  • Containment by artificial barriers
  • Management period: 300-400 years

Multi-barrier system:

  1. Waste form (solidified body)
  2. Container (metal drum)
  3. Filling material (mortar)
  4. Concrete pit
  5. Cover soil

Functions of each barrier:

BarrierMain FunctionDesign Life
Waste formFixation of radionuclides300 years
ContainerPhysical containment50 years
Filling materialVoid filling, pH buffering300 years
ConcreteStructure, water barrier300 years
Cover soilIntrusion prevention, water barrierPermanent

Applied to relatively high radioactivity waste (L1):

  • Underground 50-100 m
  • More robust barrier system
  • Management period: Several thousand years

Highly activated components of fusion reactors (first wall, divertor, etc.) can be disposed of in this category after appropriate cooling periods.

Design requirements for intermediate depth disposal:

  • Rock permeability coefficient: <107< 10^{-7} m/s
  • Groundwater flow velocity: <1< 1 m/year
  • Depth: 50-100 m

High-level waste (HLW) from fission reactors requires deep geological disposal.

ItemNear-Surface DisposalDeep Geological Disposal
DepthSeveral m to 100 m300 m or more
TargetLLW, ILWHLW
Management periodSeveral hundred yearsTens of thousands of years or more
CostLowExtremely high
Social acceptanceRelatively easyDifficult

Fusion reactors do not generate waste requiring deep geological disposal.

Reasons why deep geological disposal is not needed:

  1. Actinides are not produced
  2. Production of long-lived nuclides can be suppressed by material selection
  3. Radioactivity decays significantly with about 100 years of cooling

Estimated waste disposal costs for 1 GWe power plants:

Reactor TypeDisposal CostNotes
Fission (PWR)100-200 billion yenIncluding HLW disposal
Fusion20-40 billion yenNear-surface disposal only

Fusion reactors have a significant advantage in disposal costs.

Breakdown of disposal costs (fusion reactor):

ItemCost (billion yen)Proportion
Interim storage515%
Processing and packaging825%
Transport26%
Disposal facility construction1031%
Disposal operation515%
Post-closure management2.58%
  • Stable ground
  • Low groundwater flow velocity
  • Appropriate depth

Groundwater migration assessment (1D advection-dispersion model):

C(x,t)=C02[erfc(xvt2Dt)+exp(vxD)erfc(x+vt2Dt)]exp(λt)C(x, t) = \frac{C_0}{2} \left[ \text{erfc}\left( \frac{x - vt}{2\sqrt{Dt}} \right) + \exp\left(\frac{vx}{D}\right) \text{erfc}\left( \frac{x + vt}{2\sqrt{Dt}} \right) \right] \exp(-\lambda t)

where vv is groundwater flow velocity, DD is dispersion coefficient, and λ\lambda is decay constant.

Nuclide retardation factor (depends on distribution coefficient KdK_d):

R=1+ρbϵKdR = 1 + \frac{\rho_b}{\epsilon} K_d

where ρb\rho_b is soil bulk density and ϵ\epsilon is porosity.

  • Local community agreement
  • Accessibility
  • Long-term management system

Site selection process:

  1. Literature survey (national screening)
  2. Preliminary investigation (investigation from ground surface)
  3. Detailed investigation (boring survey)
  4. Safety review
  5. Construction and operation

Both international standards and national regulations apply to radioactive waste management of fusion reactors.

IAEA has established a safety standards framework for radioactive waste management:

  1. Safety Fundamentals (SF-1)

    • Basic principles of radiation protection
    • Intergenerational equity
  2. Safety Requirements (GSR Part 5)

    • Predisposal management of radioactive waste
    • Safety assessment requirements
  3. Safety Guides

    • Specific technical standards
    • Clearance levels (RS-G-1.7)
  1. Justification: Benefits must outweigh risks
  2. Optimization: ALARA (As Low As Reasonably Achievable)
  3. Dose limits: Limitation of individual doses

Mathematical expression (optimization):

min[S+αX+βY]\min \left[ S + \alpha X + \beta Y \right]

where SS is protection cost, XX is collective dose, YY is individual dose, and α\alpha, β\beta are cost coefficients.

Regulated by the Nuclear Regulation Authority:

  1. Act on the Regulation of Nuclear Source Material, Nuclear Fuel Material and Reactors (Reactor Regulation Act)
  2. Act on the Regulation of Radioisotopes, etc. (RI Act)
  3. Specified Radioactive Waste Final Disposal Act

Application to fusion reactors:

  • ITER: Regulated as nuclear material use facility
  • From prototype reactor onwards: Regulation equivalent to power reactor is expected

Clearance system (Article 61-2 of Reactor Regulation Act):

  • Radioactivity concentration confirmation system
  • Concentration standards by nuclide
  • Confirmation by third-party organization

Regulation under the EURATOM Treaty:

  1. Basic Safety Standards Directive (BSS Directive 2013/59/EURATOM)
  2. Radioactive Waste Directive (2011/70/EURATOM)

Member states implement through national law.

Regulated by NRC (Nuclear Regulatory Commission):

  1. 10 CFR 20: Radiation protection standards
  2. 10 CFR 61: Low-level waste disposal
  3. 10 CFR 30: Use of radioactive materials

Fusion reactor-specific regulations have not yet been developed; NRC is studying.

Current regulatory systems are based on fission reactors, and there are issues in applying them to fusion reactors:

  1. Waste classification criteria

    • Classification assuming fission products
    • Application to activation products
  2. Tritium regulation

    • Large-scale tritium use unique to fusion
    • Setting of release and management standards
  3. Clearance standards

    • Response to fusion-specific nuclides
    • Evaluation of W isotopes, etc.
  4. Licensing procedures

    • Applicability of existing procedures
    • Fusion-specific review requirements

Regulatory values are set based on risk assessment.

General public: 1 mSv/year Occupational exposure: 20 mSv/year (5-year average)

These values are derived from acceptable risk levels:

Risk=D×Risk coefficient<Acceptable risk\text{Risk} = D \times \text{Risk coefficient} < \text{Acceptable risk}

ICRP risk coefficient: 5×1025 \times 10^{-2} Sv⁻¹ (fatal cancer)

For 1 mSv/year:

Risk=103×5×102=5×105/year\text{Risk} = 10^{-3} \times 5 \times 10^{-2} = 5 \times 10^{-5} \, \text{/year}

This is comparable to other social risks (traffic accidents, etc.).

Derivation of clearance level (10 μSv/year):

  1. Set scenarios (metal recycling, landfill, etc.)
  2. Identify exposure pathways
  3. Build dose assessment model
  4. Derive levels by back-calculation
CL=10μSv/yearDCF×intake×utilization factorC_L = \frac{10 \, \mu\text{Sv/year}}{\text{DCF} \times \text{intake} \times \text{utilization factor}}

We quantitatively evaluate the characteristics of fusion reactor waste.

Radioactivity at end of operation for a 1 GWe fusion reactor (structural materials only):

Atotal10181019BqA_{\text{total}} \approx 10^{18} - 10^{19} \, \text{Bq}

This is about 1/100 of a fission reactor of the same output.

Contribution of major nuclides (immediately after shutdown):

NuclideRadioactivity Fraction
Mn-5420%
Fe-5535%
W isotopes15%
Cr-5110%
Others20%

Time evolution of radioactivity:

A(t)=A0ifieλitA(t) = A_0 \sum_i f_i e^{-\lambda_i t}

where fif_i is the initial fraction of each nuclide.

Radiotoxicity is evaluated by VDT (Volume Dilution for Toxicity):

VDT=iAiMPCi\text{VDT} = \sum_i \frac{A_i}{\text{MPC}_i}

where MPC (Maximum Permissible Concentration) is the maximum permissible concentration.

The VDT of fusion reactor waste decreases to 10⁻⁶ times that of fission reactors after 100 years.

Comparison (per 1 GWe·year):

TimeFusion Reactor VDT (m³)Fission Reactor VDT (m³)
0101310^{13}101510^{15}
100 years10910^{9}101510^{15}
1000 years10610^{6}101410^{14}

Long-term potential hazard (RI: Radiotoxicity Index):

RI(t)=iAi(t)eing,i\text{RI}(t) = \sum_i A_i(t) \cdot e_{\text{ing},i}

where eing,ie_{\text{ing},i} is the dose coefficient for ingestion.

Comparison (after 1 GWe operation):

TimeFusion Reactor RI (Sv)Fission Reactor RI (Sv)Ratio
1 year10710^710910^910210^{-2}
100 years10310^310810^810510^{-5}
1000 years10110^110710^710610^{-6}

Contact dose rate at 100 years after irradiation for each material (1 MW/m² wall load, 5-year irradiation):

MaterialCDR (μSv/h)Evaluation
316SS10⁵Remote operation required
F82H10²Limited access possible
V-4Cr-4Ti10Direct access possible
SiC/SiC1Near clearance

The effect of reduced-activation materials is clearly demonstrated.

Uncertainty in activation evaluation is mainly due to the following factors:

  1. Nuclear data uncertainty

    • Cross sections: 5-30%
    • Decay data: 1-5%
  2. Neutron spectrum uncertainty

    • Transport calculation: 10-20%
    • Geometry approximation: 5-10%
  3. Material composition uncertainty

    • Main components: 1-2%
    • Impurities: 50-100%

Propagation of total uncertainty:

σA2=i(Api)2σpi2+2i<jApiApjcov(pi,pj)\sigma_A^2 = \sum_i \left(\frac{\partial A}{\partial p_i}\right)^2 \sigma_{p_i}^2 + 2\sum_{i<j} \frac{\partial A}{\partial p_i} \frac{\partial A}{\partial p_j} \text{cov}(p_i, p_j)

Typical total uncertainty is about 20-50%.

International research and development activities on fusion reactor waste management are being advanced.

ITER will be a pioneering case for fusion reactor waste management:

  1. Activation evaluation

    • Detailed calculations with FISPACT-II
    • Implementation of uncertainty evaluation
  2. Waste classification

    • Compliance with French regulations
    • Harmonization with international standards
  3. Decommissioning planning

    • Lifecycle cost evaluation
    • Development of dismantling technology

ITER activation inventory (at end of operation):

ComponentRadioactivity (Bq)Major Nuclides
First wall5×10175 \times 10^{17}Fe-55, Mn-54
Blanket2×10182 \times 10^{18}W isotopes
Vacuum vessel3×10173 \times 10^{17}Co-60, Ni-63
Divertor1×10171 \times 10^{17}W isotopes

European fusion development is conducting waste research in coordination with DEMO design:

  1. Development of reduced-activation materials

    • Improvement of EUROFER97
    • Development of ODS steel
  2. Recycling technology

    • Remote recycling process
    • Quality assurance methods
  3. Development of regulatory framework

    • Formulation of EU unified standards
    • Harmonization of clearance systems

Activities centered on QST (National Institutes for Quantum Science and Technology):

  1. Reduced-activation materials

    • Advancement of F82H
    • Impurity management technology
  2. Activation evaluation

    • Development of ACT-4 code
    • Maintenance of JENDL/AD
  3. Neutron irradiation testing

    • Participation in IFMIF-DONES
    • Testing at domestic facilities

DOE (Department of Energy)-led programs:

  1. Material development

    • 9Cr RAFM steels
    • SiC/SiC composites
  2. Activation database

    • ENDF/B nuclear data
    • Activation calculation codes
  3. Regulatory research

    • Collaboration with NRC
    • Study of fusion-specific regulation

A problem unique to fusion reactors is the management of waste contaminated with tritium.

Tritium is:

  • Half-life: 12.3 years
  • Decay mode: Pure beta (maximum 18.6 keV)
  • Chemical behavior: Identical to hydrogen

Tritium permeation into materials:

C(x,t)=C0[1erf(x2Dt)]C(x, t) = C_0 \left[1 - \text{erf}\left(\frac{x}{2\sqrt{Dt}}\right)\right]

DD is the diffusion coefficient, and its temperature dependence is:

D=D0exp(EaRT)D = D_0 \exp\left(-\frac{E_a}{RT}\right)

Tritium diffusion coefficient in iron: D0104D_0 \approx 10^{-4} m²/s, Ea10E_a \approx 10 kJ/mol

Classification of Tritium-Contaminated Waste

Section titled “Classification of Tritium-Contaminated Waste”
Contamination LevelTritium ConcentrationManagement Method
Low< 10⁶ Bq/gManaged as normal LLW
Medium10⁶ - 10⁹ Bq/gSpecial containment
High> 10⁹ Bq/gConsider recovery and reuse

Tritium removal from tritium-contaminated waste:

  1. Thermal desorption

    • Degassing at high temperature (300-500°C)
    • In vacuum or inert gas
  2. Oxidative detritiation

    • Heating in presence of oxygen
    • Recovery as HTO
  3. Isotope exchange

    • Isotope exchange with hydrogen
    • Catalyst use

Detritiation efficiency:

η=1exp(DtL2)\eta = 1 - \exp\left(-\frac{D t}{L^2}\right)

where LL is the characteristic length of the material.

The following initiatives are important in radioactive waste management of fusion reactors.

  1. Advancement of reduced-activation materials

    • Stricter impurity management
    • Development of new materials
    • Establishment of manufacturing processes
  2. Activation evaluation technology

    • High-precision cross-section data
    • Verification of calculation codes
    • Uncertainty evaluation
  3. Dismantling and decontamination technology

    • Advancement of remote technology
    • Efficient decontamination methods
    • Waste processing technology
  4. Disposal technology

    • Optimization of disposal concepts
    • Establishment of safety assessment methods
    • Long-term stability evaluation

A regulatory framework specific to fusion reactor waste is needed:

  1. Establishment of clearance standards

    • Evaluation of fusion-specific nuclides
    • International harmonization
  2. Formulation of recycling standards

    • Reuse standards for radioactive materials
    • Quality assurance requirements
  3. Development of disposal standards

    • Classification of fusion waste
    • Requirements for disposal facilities
  4. Rationalization of licensing procedures

    • Fusion-specific review requirements
    • Phased approach
  1. Promotion of public understanding

    • Explanation of differences from fission reactors
    • Risk communication
  2. Securing disposal sites

    • Site selection process
    • Consensus building with local communities
  3. International cooperation

    • Formulation of common standards
    • Sharing of experience and knowledge

Toward practical use of fusion power generation after 2050:

  1. First generation (2050-2070)

    • Use of RAFM steel
    • Response by near-surface disposal
    • Clearance rate 50%
  2. Second generation (2070-2100)

    • Introduction of ODS steel, SiC/SiC
    • Improved recycling rate
    • Clearance rate 70%
  3. Third generation (after 2100)

    • Complete reduced-activation materials
    • Nearly complete recycling
    • Clearance rate 90%

Through these initiatives, fusion energy is expected to contribute to the realization of a sustainable society as an energy source with low environmental impact.

Advantages of fusion reactor waste management:

  1. No high-level waste
  2. Short management period (approximately 100 years)
  3. Can be handled by near-surface disposal
  4. Low disposal costs
  5. High social acceptability

These characteristics will be significant advantages in the social implementation of fusion energy.

Radioactive waste from fusion reactors can minimize the burden on the environment and future generations through appropriate material selection and management strategies. This demonstrates the important role that fusion will play in building a sustainable energy system.