Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tritium Transport”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗
DOE OSTI2023

SAM Code Development for Source Term Modeling in Fluoride-salt-cooled High-temperature Reactors

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. Many concepts leverage higher operating temperatures for improved efficiency with a variety of coolants and structural materials to support those needs. Such design choices may also pose unique radiological source term risks. There are continuous efforts to incorporate new physics and phenomena relevant to advanced reactor concepts, including capabilities to enable system-level source term modeling. The Pebble-Bed Fluoride-salt-cooled High-temperature Reactor (PB-FHR) is a promising candidate among advanced nuclear reactor concepts with its improved passive safety characteristics and high thermal efficiency. In addition to past efforts to support the development and utilization of SAM for PB-FHR safety analysis, the species transport modeling capabilities in SAM have been extended to simulate source term phenomena in concepts like the PB-FHR. The PB-FHR concept utilizes pebble-form TRISO fuel and FLiBe salt coolant to provide robust barriers to the release of almost all fission products and radiological source terms. However, tritium poses a unique radiological risk due to its significant production from neutron interactions with 6 Li and 9 Be in the FLiBe as well as its high mobility at elevated temperatures where it can permeate through structural metals. Secondly, the use of graphite at elevated temperatures poses the risk of oxidation damage to any structural components or fuel pebbles that are exposed to accidental air ingress. This report summarizes progress made in modeling these phenomena in SAM, which leverages the system-level multi-physics thermal hydraulic simulation to support an effective engineering-scale source term modeling capability. A tritium transport model is developed to simulate the various transport pathways in FHR and MSR concepts and is presented through verification and validation examples as well as demonstrations simulating experimental test loops and reference plant FHR models. A graphite oxidation model is also introduced to calculate local oxidation rates and is presented with initial validation results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Tritium Pressure Enhancement on the TPBAR Cladding by Physical Processes at the Getter

Tritium for the U.S. nuclear weapon stockpile is produced in tritium producing burnable absorber rods (TPBARs) inserted into Tennessee Valley Authoritys (TVA) light-water nuclear reactors. The rods are stainless steel tubes with a permeation barrier coating and internal components that generate and contain the tritium. The TPBAR incorporates a Ni-plated Zircoloy getter tube to capture tritium and prevent it from reaching the rod cladding and permeating into the environment. Under the conventional view of getter behavior, the tritium pressure outside the getter tube is expected to be limited to the equilibrium vapor pressure of Zr hydride at the temperature of the rod as long as the total hydrogen concentration remains below the capacity of the hydride. Since the tritium pressure is higher within the rod core, this behavior relies on the thin getters ability to hold off a differential tritium pressure. The effective tritium pressure on the cladding can also be enhanced by isotope exchange. Hydrogen ingress through the cladding from the reactor coolant creates a hydrogen pressure on the outer surface of the getter that can exchange with tritium, allowing the tritium partial pressure to increase toward this hydrogen gettering pressure. The goal of this work was to use laboratory-scale experiments to examine these mechanisms and create a model of getter behavior that describes tritium transport within the TPBAR. A third mechanism wherein the concentration at the outer surface of the getter is increased by the temperature gradient within the getter tube wall (the Soret effect) is not experimentally tested but is captured in the model. While not conclusively demonstrated by the experiments due to low pressure, high temperature, and small gap volume conditions, the model shows that when combined, the three mechanisms can explain both the magnitude and time dependence of the tritium release observed for reactor fuel assemblies with TPBARs. The model also shows how various modifications of the TPBAR design can reduce this tritium release into the environment.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

First-Principles Studies of Tritium Species Diffusivity Across the Interfaces of Ni-NiZr Alloy-Zircaloy-4

A tritium producing assembly, also known as tritium producing burnable absorber rods (TPBARs), consists of a metal getter tube located between the cladding and γ-LiAlO2 pellets. The metal getter tube is composed of a nickel (Ni) layer coated on Zircaloy-4. The getter assembly is used to capture tritium (3H) species (mainly 3H2 and 3H2O) generated from γ-LiAlO2 pellets during irradiation. Exploring 3H species (3H2, 3H2O) dissociation on the Ni surface and diffusion in the Ni layer and across the interface of Ni-plated Zircaloy-4 getters can provide insights on tritium transport and retention in the pellets and the getter materials. . In FY25, based on the ideal interface model of Ni-Zircaloy-4 generated in FY24, we will further explore the diffusion pathways of 3H across such Ni-ZrNi-alloy layer-Zircalory-4 interfaces under different conditions, including with oxide/hydroxide clusters on the Ni layer, and impurities located in ZrNi-alloy and Zircalory-4 layers. Due to the size limitation of DFT simulation, we will separate the 3-layer system into 3 sub-systems: Ni-ZrNi interface, ZrNi alloy layer, ZrNi-alloy-Zircaloy-4 interface, and simulate the diffusion barriers for tritium.

diffusion barrier↗
DOE OSTI2025

Effect of Mg and Ni impurities on tritium diffusion in lithium ceramics through cluster dynamics simulations

This study investigates the impact of Mg and Ni doping on tritium diffusion in LiAlO 2 and LiAl 5 O 8 ceramics, that are used in tritium-producing burnable absorber rods (TPBARs). Utilizing Centipede simulations across a broad temperature range (500 K to 1250 K), we explore the interplay between defect dynamics, cluster formation, and tritium mobility. In LiAlO 2 , Mg doping significantly enhances tritium diffusivity by increasing tritium interstitial concentrations and diffusion coefficients of key species, thereby doubling the overall tritium diffusivity. Ni doping, while shifting the dominant defect to Li vacancies, maintains high tritium mobility due to the low binding energy of Li vacancy-tritium complexes, which ensures effective tritium migration. In LiAl 5 O 8 , Mg and Ni doping results in a slight reduction in the diffusion coefficients of key species, yet the dramatic increase in tritium interstitial concentrations compensates, leading to a net small increase in tritium diffusivity. In conclusion, the findings highlight the critical role of defects in tritium transport and the effect of Mg and Ni defects on the performance of these ceramics in demanding nuclear environments.

Cluster dynamics↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Deuterium Extraction from Helium with a Vanadium Vacuum Permeator: Commissioning of the Tritium Extraction eXperiment (TEX)

Here, the Tritium Extraction eXperiment (TEX) is a forced-convection lead-lithium (PbLi) loop in the Safety and Tritium Applied Research (STAR) facility at Idaho National Laboratory (INL) with the purpose of providing validation data for the vacuum permeator tritium extraction concept. A vanadium tube of 1000 mm length, 12.7 mm outside diameter, and 0.50 mm wall thickness is installed in the test section of TEX. The installed vanadium tube is characterized to quantify impurity concentrations, surface chemistry, and microstructure to elucidate permeation phenomena observed in experimentation. Herein, the permeation properties of the vanadium tube are characterized by measuring deuterium permeation at 300 °C, 325 °C, and 350 °C at 100 kPa, 125 kPa, and 150 kPa total pressures with 5000 ppm deuterium in helium gas mixture in a once-through flow configuration. The hydrogen isotope permeation through the vanadium tube in the test section is measured with quadrupole mass spectrometers and the hydrogen isotope concentration in the feed and retentate gas stream is measured with gas chromatography. The transient permeation results are modeled with MELCOR-TMAP, a thermal-hydraulic tritium transport code, and compared well with literature data.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

MOOSE ProbML: Parallelized probabilistic machine learning and uncertainty quantification for computational energy applications

Here, this paper presents the development and demonstration of massively parallel probabilistic machine learning (ML) and uncertainty quantification (UQ) capabilities within the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source computational platform for parallel finite element and finite volume analyses. In addressing the computational expense and uncertainties inherent in complex multiphysics simulations, this paper integrates Gaussian process (GP) variants, active learning, Bayesian inverse UQ, adaptive forward UQ, Bayesian optimization, evolutionary optimization, and Markov chain Monte Carlo (MCMC) within MOOSE. It also elaborates on the interaction among key MOOSE systems — Sampler, MultiApp, Reporter, and Surrogate — in enabling these capabilities. The modularity offered by these systems enables development of a multitude of probabilistic ML and UQ algorithms in MOOSE. Example code demonstrations include parallel active learning and parallel Bayesian inference via active learning. The impact of these developments is illustrated through five applications relevant to computational energy applications: UQ of nuclear fuel fission product release, using parallel active learning Bayesian inference; very rare events analysis in nuclear microreactors using active learning; advanced manufacturing process modeling using multi-output GPs (MOGPs) and dimensionality reduction; fluid flow using deep GPs (DGPs); and tritium transport model parameter optimization for fusion energy, using batch Bayesian optimization. These capabilities are part of the MOOSE framework.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

SAM Code Development for Transient Safety Analyses of Fluoride-salt-cooled High-temperature Reactors

The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. The molten-salt-cooled pebble-bed reactor, or pebble-bed FHR (PB-FHR) is a promising candidate among advanced nuclear reactor concepts with its improved passive safety characteristics and high thermal efficiency. To support the development and utilization of the SAM code for PB-FHR safety analysis, activities on SAM code enhancements, reference plant model developments, and code validations have been performed in the past a few years to support near-term industry and NRC needs. This report summarizes recent progress under DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program in SAM code development and demonstration for transient safety analysis of Fluoride-salt-cooled High-temperature Reactors. SAM capabilities has been significantly enhanced over the years to add FHR specific modeling features, including salt freezing and thawing, spherical core channel and pebble bed core modeling, solid-fluid thermal radiation, tritium transport and general species transport in fluids and solids, and the general code enhancements on solver schemes of point kinetics module and reactivity feedback models. A reference PB-FHR model is developed, based on publicly available information from Kairos Power’s generic FHR design and the University of California, Berkeley (UCB) Mk1 design. A reference reactor model is foundational to the methodologies employed by NRC to verify the adequacy of computer codes and evaluation models. The reference FHR model was utilized for a number of selected FHR design basis accidents, including station blackout, loss of heat sink, loss of flow, transient overpower, and overcooling events.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Estimation of China Fusion Engineering Test Reactor performance and burning fraction in different pellet fueling scenarios by a multi-species radial transport model

Tritium self-sufficiency in future deuterium–tritium fusion reactors is a crucial challenge. As an engineering test reactor, the China Fusion Engineering Test Reactor requires a burning fraction of 3% for the goal to test the accessibility to the future fusion plant. To self-consistently simulate burning plasmas with profile changes in pellet injection scenarios and to estimate the corresponding burning fraction, a one-dimensional multi-species radial transport model is developed in the BOUT++ framework. Here several pellet-fueling scenarios are then tested in the model. The results show that the increased fueling depth improves the burning fraction by particle confinement improvement and fusion power increase. Nevertheless, by increasing the depth, the pellet cooling-down may significantly lower the temperature in the core region. Taking the density perturbation into consideration, the reasonable parameters of the fueling scenario in these simulations are estimated as pellet radius ${r}_{{\rm{p}}}=3\,{\rm{mm}},$ injection rate $=\,4\,\mathrm{Hz},$ and pellet injection velocity $=\,1000\mbox{-}2000\,{\rm{m}}\,{{\rm{s}}}^{-1}$ without drift or $450\,{\rm{m}}\,{{\rm{s}}}^{-1}$ with high-field-side drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of Transport Processes in Liquid-Metal Fusion Blankets: Past, Present, and Future

The successful development of robust breeding blanket systems will strongly rely on computational tools for predicting the complex behavior of the electrically conducting liquid-metal (LM) breeder flowing in the complex-shaped blanket ducts in the presence of a strong plasma-confining magnetic field, volumetric heating, and tritium generation. Associated transport processes involve magnetohydrodynamic (MHD) flows, heat transfer, corrosion, and tritium transport. This paper is an overview of past and present efforts in the development, application, and verification and validation (V&V) of such computational tools. As a result of the ongoing campaign on V&V of computer codes for LM blankets, the international fusion community has identified several candidates that promise to become real blanket design and analysis tools in the near future. Among them are HIMAG, MHD-UCAS, COMSOL Multiphysics, ANSYS FLUENT, ANSYS CFX, and OpenFOAM. The progress, over the last decade, in the application of such codes in blanket studies is tremendous. This is illustrated with two examples for a dual-coolant lead-lithium (DCLL) blanket: (1) integrated computer modeling for the recently designed DCLL blanket in the United States and (2) application of the code MHD-UCAS to the analysis of PbLi flows and heat transfer in a generic DCLL blanket prototype at high Hartmann (Ha ~ 10 4 ) and Grashof numbers (Gr ~ 10 12 ). Here, this paper also presents an approach to the development of a new integrated computational tool called the virtual dual-coolant lead-lithium (VDCLL) blanket, which elaborates the existing U.S. MHD code HIMAG.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Thin Alumina Film as a Tritium Adsorption Inhibitor for Stainless Steel 316

The effect of a thin alumina coating on stainless-steel 316 (SS316) samples on tritium adsorption and transport are discussed. Compact films of alumina were produced on the surfaces of pristine SS316 samples using an atomic layer deposition (ALD) technique. Subsequently, these samples were exposed for 24 hours to a deuterium-tritium gas mixture (P T =0.5 atm, 25°C). A combination of methods including selective-etching and programmed thermal desorption were employed to assess both, the depth profile of the tritium concentration in the sample and the total quantity of tritium absorbed respectively. Tritium was quantitatively determined through the measurement of beta radioactivity using liquid-scintillation counting techniques. Data suggest that SS316 with a thin film of alumina reduces the total tritium uptake by ~25% relative to uncoated samples. Significantly, such films appear to reduce, by a factor of 200, tritium diffusion into SS316 and therefore constitute an effective barrier against tritium transport. This observation is of practical importance for tritium and, generally, reactive gas handling.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Isotope mass scaling and transport comparison between JET Deuterium and Tritium L-mode plasmas

The dimensionless isotope mass scaling experiment between pure Deuterium and pure Tritium plasmas with matched $\rho^*$, $\nu^*$ , β n , q and $T_e/T_i$ has been achieved in JET L-mode with dominant electron heating (NBI+ohmic) conditions. 28% higher scaled energy confinement time $B_t\tau_{E,\mathrm{th}}/A$ is found in favour of the Tritium plasma. This can be cast in the form of the dimensionless energy confinement scaling law as $\Omega_i \tau_{E,\mathrm{th}} \sim A^{0.48 \pm 0.16}$. This significant isotope mass scaling is consequently seen in the scaled one-fluid heat diffusion coefficient $A \chi_{\mathrm{eff}}/B_t$ which is around 50% lower in the Tritium plasma throughout the whole plasma radius. The isotope mass dependence in the particle transport channel is negligible, supported also by the perturbative particle transport analysis with gas puff modulation. The comparison of the edge particle fuelling or ionisation profiles from the EDGE2D-EIRENE simulations show that the absolute density differences that are necessary for the dimensionless match in the confined plasma dominate over any isotope mass dependencies of particle fuelling and ionization profiles at the plasma edge. Local GENE simulation results indicate a mild anti-gyroBohm effect at $\rho_\mathrm{tor}$ = 0.6 and thereby a small isotope mass dependence in favour of Tritium on heat transport and a negligible effect on particle transport. A significant fraction of the isotope scaling and reduced heat transport observed in the Tritium plasma is not captured in the GENE and ASTRA-TGLF-SAT2 simulations by simply changing the isotope mass for the same input profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Species Transport Framework Development in SAM for System-Level Tritium Source Term Analysis

The SAM code is under development as a modern system-level modeling and simulation tool for advanced non–light water reactor safety analyses, with recent efforts to add capabilities to evaluate radiological source term risks in these novel reactor concepts. By leveraging the established system-level multiphysics thermal-hydraulic models in SAM, a framework for tightly coupled species transport modeling has been integrated into the code for engineering-scale source term evaluation. This species transport framework was first applied to the simulation of tritium, which is a well-known source term in conventional light water reactors. Tritium poses a unique risk in salt-cooled reactors, especially those with lithium-bearing salts such as the fluoride salt–cooled high-temperature reactor (FHR) concept, as tritium is generated in the salt coolant in significant quantities due to neutron interactions. A compounding factor is the increased mobility of tritium at high temperatures, which is able to permeate through metals while also potentially being retained in graphite pebbles and structures. Engineering-scale models for the tritium transport pathways in a FHR have been developed using the new species transport framework in SAM. The capabilities are assessed through analytical verification problems and validated with data from a graphite retention experiment. In conclusion, the system-level model is demonstrated by performing an initial estimate of baseline tritium generation and flows in a generic reference SAM FHR model, setting a foundation for future studies of source term transient analysis with the potential for further multiscale and multiphysics integration.

SAM↗
DOE OSTI2021

Asymmetry between deuterium and tritium turbulent particle flows

The asymmetry between deuterium (D) and tritium (T) turbulent particle fluxes in mixed D–T plasmas is studied with numerical simulations of nonlinear gyrokinetic turbulence in ion temperature gradient-dominated and trapped electron mode-dominated regimes. At 50–50 D–T concentration, the asymmetry, or flow separation, between D and T fluxes is such that the tritium is better confined than the deuterium in both regimes. To supplement the nonlinear simulations, an analytic quasilinear theory of the particle flux symmetry breaking is developed and is valid for general electron dynamics. Furthermore, this theory correctly predicts the numerically computed deviation in the ion density fraction (from 50–50), or ion density gradient (from the electron gradient), required to restore equal deuterium and tritium fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗