Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion barriers”

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.

194 records · Page 11

Hydrogen permeation in iron-chromium-aluminum (FeCrAl) alloys and the effects of microstructure and surface oxide

Iron–chromium–aluminum (FeCrAl) class alloys are candidates for use as cladding for accident-tolerant fuels and moderators. In this context, hydrogen isotope permeation in FeCrAl alloys is an important material property. Here, in the present work, the apparent permeability, effective diffusivity, and apparent solubility of hydrogen in the FeCrAl alloys C26M and Kanthal D (KD) were measured with gas-driven hydrogen permeation. Permeation measurements were conducted at temperatures of 400 to 700 °C and at gas-driven pressures from 1 to 100 kPa. In particular, the effect of grain size on hydrogen transport was studied with KD samples with three different microstructures: nanocrystalline (NC), ultra-fine grained (UFG), and coarse-grained (CG). The UFG and NC specimens had higher apparent activation energies (73.4 kJ mol -1 and 65.2 kJ mol -1 , respectively) for hydrogen permeability than the CG sample (46.9 kJ mol -1 ). An aluminum oxide layer formed on the primary- and secondary-side surfaces of all samples subjected to permeation experiments which demonstrated the propensity of FeCrAl alloys to form these innate oxide permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Breaking Barriers in Chalcogenide Perovskite Synthesis: A Generalized Framework for Fabrication of BaMS 3 (M═Ti, Zr, Hf) Materials

Abstract Chalcogenide perovskites have garnered increasing attention as stable, non‐toxic alternatives to lead halide perovskites. However, their conventional synthesis at high temperatures (>1000 °C) has hindered widespread adoption. Recent studies have developed low‐to‐moderate temperature synthesis methods (<600 °C) using reactive precursors, yet a comprehensive understanding of the pivotal factors affecting reproducibility and repeatability remains elusive. This study delineates the critical factors in the low‐temperature synthesis of BaMS 3 (M═Zr, Hf, Ti) compounds and presents a generalized framework. Innovative approaches are developed for synthesizing BaMS 3 compounds using this framework involving organometallics for solution deposition. The molecular precursor routes, employing metal acetylacetonates to generate soluble metal–sulfur bonded complexes and metal–organic compounds to produce soluble metal‐thiolate, metal‐isothiocyanate, and metal‐trithiocarbonate species, are demonstrated to yield carbon‐free BaMS 3 . These methods have achieved the most contiguous films of BaZrS 3 and BaHfS 3 using solution deposition to date. Furthermore, a hybrid solution processing method involving stacking sputter‐deposited Zr and solution‐deposited BaS layers is employed to synthesize a contiguous, oxygen‐free BaZrS 3 film. The diffuse reflectance measurements indicate a direct bandgap of ≈ 1.85 eV for the BaZrS 3 films and ≈ 2.1 eV for the BaHfS 3 film under investigation.

25 ENERGY STORAGE↗

Multilayer Architecture Barriers for High Temperature Nuclear Applications

Innovative thin film multilayer metal/ceramic barrier coatings have been developed by Argonne National Laboratory to protect fuel cladding against different phenomena that can be encountered under the harsh operating conditions of advanced nuclear reactors. One phenomenon is Fuel-Cladding Chemical Interaction (FCCI), which can represent a challenge to the performance of steel-clad metallic fuels in sodium-cooled fast reactors (SFRs), especially when operating to high burnups at elevated cladding temperatures, and long fuel residence time. Another phenomenon is the corrosion of zirconium-based cladding under prototypic boiling water reactor (BWR) conditions and associated cladding during off normal conditions, including loss of coolant accidents (LOCA). This work provides details of the development and testing of a thin film (few microns) multilayer metal/ceramic (CrY/Y 2 O 3 ) coatings on steel cladding and its performance against FCCI. The coating architecture evolved through systematic mechanical testing, radiation tolerance assessments, and high-temperature diffusion studies from a preliminary multilayer concept to an optimized high-performance design. The optimized structure leverages stoichiometric Y 2 O 3 for its exceptional radiation stability, chemical inertness, and ultra-low fission product diffusivity (e.g., Ce diffusivity ~10 -20 m 2 /s in irradiated Y 2 O 3 at 650 °C). This is paired with ductile chromium-rich interlayers (Cr 95 Y 5 ), which significantly enhance fracture toughness and mechanical resilience. Design refinements, including increased ceramic layer thickness and optimized Cr–Y interlayer composition, eliminated fracture propagation and preserved coating integrity under extreme mechanical stress (>10× typical reactor conditions) and prolonged radiation exposure (~300 dpa). Those conditions are well beyond expected conditions during normal operations in SFR, and cover some of the transient conditions as well. In addition, the coating deposition was demonstrated to be applied uniformly to the inner surface of steel tubes representing prototypic SFR cladding. Another phenomenon under consideration in this work is the corrosion of zirconium-based cladding under prototypic boiling water reactor (BWR) conditions and associated off normal conditions including loss of coolant accident (LOCA). Another multilayer thin film coating is considered for mitigation of this phenomenon, where a coating of yttrium aluminum and yttrium aluminum ceramic/metal thin film coating. Preliminary results on the coating performance during high temperature thermal cycling (up to 1000 C) have shown promising performance of the coating. On-going optimization of the coating will continue in FY26.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development Fiber Optic Distributed System for Direct Detection of Subsurface Gases Leakages

Carbon, natural gas, and hydrogen gas storage is an emerging solution to safeguard us against pollution, support goals of negative carbon emission, and protect sources of renewable energy. Properly constructed storage wells provide a virtually impervious barrier to any unintended subsurface transmission. The ability to ensure the long-term integrity of such wells is vital to the success of any storage operation and be successful in the public eyes. Therefore, robust monitoring of any gas migration into the subsurface is highly sought. A fiber-optic distributed chemical sensor (DCS) enables monitoring of long-term well integrity along its depth, ensuring the success of any storage operation and bolsters public acceptance of the safety of the reservoir via leak early detection. The same technique can be applied to gas monitoring in pipeline networks and nuclear stockpile monitoring applications. Fiber based Raman spectroscopy enables DCS, as optical fibers can be deployed in virtually any environment and relay spectroscopic information over long distances back to the user. Hollow core fibers (HCF) make excellent DCSs as the air core of the fiber allows gas from the environment to diffuse into the core, which interacts with the laser signal that is carried in the air core. This work builds upon the previous LDRD project, Fiber Optic System for Direct Detection of Carbon Dioxide Leakage in Carbon Storage Wells (21-FS-003), in which the feasibility of Raman spectroscopy detection of Carbon Dioxide (CO2) in HCF detection was demonstrated. We mitigated the risk of this DCS technology by establishing and completing five objectives. The first objective was to model and optically characterize HCF uptake of CO2, establishing the relationship between HCF length, gas diffusion time, detectable gas concentration, and measured Raman intensity. In objective two, we developed a fiber core drilling recipe to enable additional diffusion ports in the fiber core and established a method for maintaining fiber strength and integrity post drilling. Objective three characterized the drilled fibers against the undrilled fibers, establishing the differences in the gas mechanics and optical properties and provided parameters to iterate the drilling process. In objective four, a fusion splicing technique was developed to join the HCF to conventional single-mode fibers, localizing the gas detection point at the drilled HCF hole, emulating a DCS. Lastly, objective five was the testing of the sensor in Edgar Mines at Colorado School of Mines on a CO2 pipeline with a simulated leak, to showcase the ability to detect CO2 leaks. This capstone result showed CO2 leak detection in < 10 minutes, raising the technology readiness level of HCF segments as deployable DCS.

organic↗

Promotional Effect of ZnO and ZrO 2 in K-Doped Fe Catalysts for CO 2 Hydrogenation to Light Olefins

Promoters play a critical role in tuning the activity and selectivity of Fe catalysts in CO 2 hydrogenation to produce light olefins, which are key building blocks in the petrochemical industry. Herein, by a combined experimental and theoretical approach, we show that high and stable performance of Fe catalysts could be achieved by taking advantage of the promotional effect of both Zn and Zr. Structural characterization indicates that ZnO could improve the dispersion and reducibility of Fe oxides and facilitate the formation of active Fe carbide species, whereas ZrO 2 could stabilize the structure and catalytic performance, especially the selectivity of hydrocarbon products. In situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) experiments suggest that carbonate, bicarbonate, formate, and methoxy are essential intermediates in the CO 2 hydrogenation to hydrocarbon products, including paraffins and olefins. The conversion kinetics of each intermediate species are dependent on the type of promoters as well as the phase structure of the active Fe species. DFT calculations revealed a strong correlation between the formation energy of surface oxygen vacancies and that of Fe carbide species in promoted Fe oxides, in accordance with experimental results. Moreover, the calculated energy profiles of CO 2 hydrogenation over different catalysts indicate that Zn could promote the activation of CO 2 and its transformation to the oxygenate intermediates, while Zr could facilitate the conversion of oxygenates to hydrocarbon precursors. Finally, the discrepancies in the evolution trend of various intermediate species on promoted Fe catalysts in the transient DRIFTS experiments can be rationalized by the differences in energy barriers of elementary or rate limiting steps.

CO2 conversion↗

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗

Understanding and Tailoring Diffusion and Co-Adsorption Inside the Confined Pores of Metal-Organic Frameworks (Final Scientific/Technical Report for Award DE-SC0019902)

The aim of this program was to gain a fundamental understanding of the behavior of various guest molecules in nano-confined environments, such as metal organic frameworks (MOFs), using a combination of novel synthesis, ab initio modeling, and in situ characterization. Through this project, we developed a concise understanding of the mechanisms that control adsorption/desorption of gaseous molecules and their mixtures, leading to design/synthesis guidelines for MOFs with desired functionality. We further developed methods to disentangle kinetic from thermodynamic effects during adsorption, as well as to characterize the interactions at play. In the first funding cycle, the focus was on the unambiguously characterization of co-adsorption and diffusion of gasses/vapors and their mixtures. In the second funding cycle, the focus was on characterizing the effects of the nano-confinement on the kinetics and thermodynamics of adsorption processes inside MOFs, again with an emphasis on mixtures of gasses and vapors. The nano-confinement can tip the thermodynamic vs. kinetic balance, and current understanding and theory based on single-component analysis can lead to incorrect predictions for mixtures. This is of particular interest in real-world applications, where gasses/vapors are typically mixed, contain impurities, or are often exposed to humid conditions. Our main findings were: (i) within confined environments the adsorption behavior of mixed gasses/vapors can be drastically different from the “sum” of the corresponding single phases; (ii) co-adsorption is often competitive and detrimental to performance, but it can also be cooperative and beneficial; (iii) in some co-adsorbed gasses/vapors, molecules that are strongly bound in the single-component phase can be replaced by molecules that are nominally weaker bound (molecular exchange) due to guest-guest interactions that lower the kinetic barriers and favor the final adsorption state; (iv) kinetic and thermodynamic effects can be precisely controlled through pore-size engineering and synthesis; and, (v) kinetic effects can be identified and disentangled from thermodynamic effects during adsorption through a series of sequential and simultaneous gas loading measurements. The short-term goal of this program was the controlling and understanding of common MOF systems in real-world situations where gasses/vapors are mixed, which will have an important impact on industrial processes and applications from gas storage and sequestration to catalysis and sensors. The long-term goals include the development of theoretical and experimental methods for gaining a fundamental understanding of adsorption/reaction processes within MOFs, as well as new guidelines for synthesizing MOFs with tailored physical and chemical properties.

36 MATERIALS SCIENCE↗

Optical properties of vacancies in aluminum oxide (𝛼−Al 2 ⁢O 3 ) from first principles

We employ first-principles calculations based on hybrid density functional theory to investigate the structural and optical properties of the oxygen vacancy (𝑉 O ) and aluminum vacancy (𝑉 Al ) in 𝛼−Al 2 ⁢O 3 , the most stable (corundum) phase of alumina. Our calculations facilitate the identification of experimental excitation and luminescence spectra with specific electronic transitions at the vacancy sites. The absorption line shape for excitation of an electron at 𝑉$^0_O$ to the conduction-band minimum (CBM) is in excellent agreement with the 6.1 eV band detected by optical absorption spectroscopy, and we find that the 5.9 eV absorption band is generated by an internal electron transition at 𝑉$^0_O$. We confirm that the slowly decaying 3.0 eV emission band of the 𝐹 center is due to a triplet-singlet transition at 𝑉$^0_O$. Our calculations also reveal that the 4.8 eV/5.4 eV absorption and 3.8 eV emission bands assigned to the 𝐹 + center are generated by internal transitions at 𝑉$^{+1}_O$. The line shape for the excitation of an electron from 𝑉$^{+1}_O$ to the CBM agrees well with an absorption band at 6.4 eV, while the recombination of an electron with 𝑉$^{+2}_O$ also produces luminescence around 3.8 eV, but with a considerably broader line shape. For Al vacancies, we confirm that the most prevalent configuration in 𝛼−Al 2 ⁢O 3 is a split-vacancy configuration 𝑉 Al,s , and we calculate migration barriers for different directions in the corundum crystal. We predict absorption and emission spectra for the excitation and recombination of an electron localized at 𝑉$^{−3}_{Al}$ and 𝑉$^{−3}_{Al,s}$ sites with the CBM. The line shapes of the two 𝑉 Al configurations overlap and are considerably broader than the spectra corresponding to 𝑉 O .

36 MATERIALS SCIENCE↗

Structure–Activity Relationships for Ethanol Dehydrogenation to Acetaldehyde by Silica-Supported Zinc Oxide Catalysts

Silica-supported ZnO efficiently catalyzes the nonoxidative dehydrogenation of ethanol to acetaldehyde, which is relevant for production of 1,3-butadiene from bioethanol. Characterization with in situ spectroscopies under dehydrated conditions (high sensitivity-low energy ion scattering (HS-LEIS), diffuse reflectance (DR) UV–vis, X-ray absorption spectroscopy (XAS), diffuse reflectance Fourier transform infrared spectroscopy (DRIFTS), inelastic neutron scattering (INS), and UV Raman), and ammonia adsorption probed by temperature-programmed desorption followed by DRIFTS and mass spectrometry (DRIFTS-MS NH 3 -TPD), and DFT calculations revealed that the supported ZnO x phase was present as isolated surface ZnO x sites on SiO 2 , with the vast majority coordinated by two siloxane bonds and one silicon atom with two nonbridging oxygens ((≡SiO) 2 Zn 2+ O 2 Si=), anchored at 4-, 5-, and 6-membered siloxane rings. A minor fraction of surface ZnO x sites possessed Lewis acidity, and even fewer sites possessed a Bro̷nsted acidic Zn(OH) + Si moiety. Ethanol temperature-programmed surface reaction-mass spectrometry (TPSR-MS) with various oxidative or ethanol reaction pretreatments indicated that only sites with Lewis and Bro̷nsted acidic character (Zn(OH) + Si) were active for ethanol dehydrogenation, while the majority surface (≡SiO) 2 Zn 2+ O 2 Si= sites were inactive. Greater heterogeneity among all surface ZnO x sites, as assessed by in situ DR UV–vis spectroscopy, was associated with a greater number of ZnO x sites that were active for ethanol dehydrogenation as well as lower enthalpic barriers for acetaldehyde production among the most active surface ZnO x sites. Turnover frequencies and the apparent activation energy for ethanol dehydrogenation were determined from steady-state kinetics. Together, these findings suggested that anchoring inactive surface (≡SiO) 2 Zn 2+ O 2 Si= sites on the silica support caused a greater number of active surface ZnO x sites to adopt a more strained configuration, promoting ethanol dehydrogenation catalysis. Pretreatments and catalysts that promoted desorption of ethanol during TPSR, taken as a marker of surface dehydroxylation, were associated with an increased number of the most active surface (Zn(OH) + Si) sites. Such findings suggested that inactive surface ZnO x sites were activated for ethanol dehydrogenation by dehydroxylation of the support and/or decreased coordination to hemilabile siloxane ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature Evolution of the Activation Barriers Leads to Meyer−Neldel Rules for Structural Relaxation and Transport in Polymers

Understanding activation barriers controlling structural relaxation in glass-forming liquids, molecular transport, and ionic conductivity in amorphous polymers is a grand challenge of fundamental scientific and materials engineering interest across disciplines. Over decades, intriguing but puzzling empirical correlations between the elementary time scale of activated barrier crossing and the apparent Arrhenius activation energy, the so-called Meyer−Neldel (MN) rules, have been discovered in diverse liquids and glasses. Here, in this study, we formulate and successfully apply a new experimental analysis and an explicitly dynamical theoretical framework which provides an understanding of the origin, validity, and failure of such correlations, that bridge and unify the three fields of structural relaxation, molecular transport, and ionic conductivity in liquids and quenched glasses. Distinct quasi-universal laws are predicted in equilibrated liquids and nonequilibrium glasses, consistent with experiments. Our analysis reveals that even if the relaxation appears Arrhenius over a limited temperature range, the physical activation barrier is generally temperature-dependent in polymeric systems even below glass transition temperature. In addition, we show that the approximate validity of classical MN rules hinges on a linear temperature dependence of this barrier and the temperature range probed in experiments. Our findings are relevant for controlling the activation barrier in functional soft polymeric materials relevant to molecular separations, barrier coatings, and charge transport, and also provide new constraints on the theoretical understanding of the mechanism underlying slow activated dynamics in glass-forming condensed matter.

Meyer-Neldel rules↗

Phase Field Modeling of Chemical Reaction Related Damage Evolution in Environmental Barrier Coatings

The advent of next-generation engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the thermally grown oxides (TGO) and phase transformation in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and porosity evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

fast-diffusion path↗

Flexibility of oxygen sublattice and hydrogen bond length predict proton mobility in ternary metal oxides

Discovery of fast proton conductors is important for advancing clean energy technologies. This requires a better understanding of proton migration mechanisms. While structural and chemical traits of ternary metal oxides have been related to proton migration barriers, lattice dynamical effects have not been resolved quantitatively. Here, in this work, we introduce a phonon-based dynamic descriptor, termed ‘‘thermal O…O fluctuation,’’ quantifying the flexibility of donor-acceptor oxide-ion pairs. This enables direct comparison of O-sublattice flexibility across diverse metal oxides. Using regression models, we ranked physical descriptors as predictors of proton mobility, finding that H-bond length and thermal O…O fluctuation were the strongest descriptors. Further analysis revealed a critical O…O spacing of 2.4 A˚ at the transition state, which is easier to reach by more flexible donor-acceptor pairs, enabling facile proton transfer. Our results demonstrate oxygen sublattice flexibility as a dynamic descriptor and provide guiding principles for enhancing proton mobility in ternary metal oxides.

diffusivity↗

Materials Science of the Interstitial Doping Process

Particle accelerators are an increasingly important tool for frontier science. Growing initial and operating costs are a significant barrier for upgrades and for new machines. While everything matters, the major cost contributor is the SRF cavities and their ancillary facilities (e.g., cryoplant). Accordingly, the accelerator science community devotes much R&D effort to improving their energy efficiency (increased Q o ) and gradient (E acc ). While improved gradient is at the forefront for certain machines (ILC), improved quality factor has broader impact, benefitting all SRF applications. An important opportunity for accelerator science and technology to move forward arose in the course of building the LCLS-II, the second generation Linac Coherent Light Source at SLAC. At more or less the same time, researchers at Fermilab discovered that introducing a small amount of nitrogen to the niobium surface could improve the mid-range quality factor as much as three-fold. A firm resolution of how nitrogen confers its benefit attracts much current research interest. The key elements of the “nitrogen doping” process were vacuum bake at 800 °C, brief exposure to mTorr of nitrogen at several hundred degrees followed by electropolish (EP) to remove several microns from the surface to eliminate unwanted nitrides: While the process as a whole was novel, the comprising unit operations are familiar to the accelerator community. It was judged reasonable to adopt it as a cost-reduction technology for LCLS-II. Researchers carried out a program of varying process parameters and measuring performance in single-cell cavities, leading to a consensus stable protocol for the project. The transition to vendor fabrication and multiple niobium sources has presented unforeseen challenges of performance variation evidently not connected to anything that could be incorporated in a purchase specification. Moving beyond the high temperature N-doping process above, researchers reported a simplified process consisting entirely of tens of hours anneal in a N atmosphere at low temperatures (~120°C – 160°C) after 800 °C UHV bake. These processes typically yield a few-nm doped layer while the high temperature process yields at least a few to many micron doped layer. Even more recently, oxygen has been used to dope instead of nitrogen which leaves a few-µm O-alloyed layer upon vacuum annealing for 300 °C for ~3 hours. The investigation of these materials is just beginning, but the process simplification they may offer is surely attractive. The very low quantity of material that appears to be significant in the “infusion” process indicates the need for very careful control of gas species available for diffusion into the surface during low temperature treatment, both for process control and research to characterize the underlying dynamics. Oxygen alloying offers the further opportunity to utilize the decomposition of the surface native oxide as the dopant source. It is necessary to understand and (thus) manage this process. We have been supported by the Department of Energy Offices of High Energy Physics and Nuclear Physics to pursue this goal.

36 MATERIALS SCIENCE↗

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗