Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility

Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.

buoyant mixing↗
DOE OSTI2024

Complete Development of Critical Capabilities for TRISO Fission Product Source Term Calculations and Quantify Mechanisms for Pd Penetration of SiC

Overall fission product (FP) release will be an important consideration for the licensing and deployment of advanced reactors utilizing tristructural isotropic (TRISO) fuels. This work focuses on enhancing and applying the BISON models needed to predict FP transport within TRISO particles and particle failure probability, both of which factor directly into release predictions. Specifically, this report details (1) the development of the models needed to predict palladium (Pd) conservation at the engineering scale and the application of those models to characterize Pd fluxes for input into a mechanistic multiscale model for Pd penetration; (2) the refinement of sorption mass transfer models and the development of models for trapping in porous layers, which were applied and compared to particle scans from AGR-2 to provide proof of concept for a method of particle-scale validation that may reduce uncertainties compared to compact-scale validation using data from integral effects tests; (3) the development of a failure-statistics-informed, mesh-independent methodology for applying smeared cracking, enabling further study of the localized multiphysics behaviors associated with cascading particle failure mechanisms; and (4) the preliminary characterization of those coupled multiphysics particle failure behaviors using smeared, nonretentive diffusivities to provide a baseline for future study and to guide ongoing engineering applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modeling the Charge and Exciton Distributions in Phosphorescent White Organic Light-Emitting Diodes

A comprehensive numerical model is developed to quantitatively analyze the charge and exciton distributions in phosphorescent white organic light emitting diodes (WOLEDs). The model encompasses bipolar charge drift-diffusion transport with doping-dependent charge mobilities, field-dependent charge hopping across multiple interfaces, exciton generation, diffusion, energy transfer, radiative decay, and non-radiative exciton annihilation. The model is validated by direct measurement of the exciton distribution in monochromatic phosphorescent OLEDs. It is then applied to complex, multilayer WOLEDs where it is used to accurately predict their intensity and current-dependent spectral shifts. The magnitudes of different exciton loss channels in WOLEDs are calculated and quantitatively compared. This model provides an alternative for the many experimental iterations commonly used to optimize the color balance of complex OLEDs. In this context, the model inspired the addition of a hole blocking layer in the WOLED to balance the charge conduction and reduce spectral shifts – a design that is verified by experiment.

42 ENGINEERING↗
DOE OSTI2021

Deterioration mechanism of alkali-activated materials in sulfuric acid and the influence of Cu: A micro-to-nano structural, elemental and stable isotopic multi-proxy study

In this study, a multi-proxy approach combining {sup 29}Si, {sup 27}Al and {sup 1}H MAS-NMR, FEG-EPMA, XANES at the Cu K-edge and XRD analytics with hydrochemical tools such as ICP-OES analyses, oxygen-isotope signatures, and thermodynamic modelling was applied to K-silicate-activated metakaolin specimens - with and without CuSO{sub 4}·5H{sub 2}O addition - exposed to sulfuric acid at pH = 2 for 35 days. The results revealed a multistage deterioration mechanism governed by (i) acid diffusion, (ii) leaching of K-A-S-H, (iii) microstructural damage related to precipitation of expansive (K,Ca,Al)-sulfate-hydrate phases (iv) complete dissolution of the K-A-S-H framework, (v) and formation of silica gel in the outermost corroded regions. Copper ions were mainly located in layered spertiniite-chrysocolla-like phases in the as-cured materials. The results demonstrate an overall negative effect of Cu addition on chemical material durability, implying that the reported higher durability of Cu-doped AAM in biocorrosion environments can be best explained by bacteriostatic effects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Terrestrial photosynthesis inferred from plant carbonyl sulfide uptake

Terrestrial photosynthesis, or gross primary production (GPP), is the largest carbon flux in the biosphere, but its global magnitude and spatiotemporal dynamics remain uncertain. The global annual mean GPP is historically thought to be around 120 PgC yr −1 , which is about 30–50 PgC yr −1 lower than GPP inferred from the oxygen-18 ( 18 O) isotope and soil respiration. This disparity is a source of uncertainty in predicting climate–carbon cycle feedbacks. Here, in this study, we infer GPP from carbonyl sulfide, an innovative tracer for CO 2 diffusion from ambient air to leaf chloroplasts through stomata and mesophyll layers. We demonstrate that explicitly representing mesophyll diffusion is important for accurately quantifying the spatiotemporal dynamics of carbonyl sulfide uptake by plants. From the estimate of carbonyl sulfide uptake by plants, we infer a global contemporary GPP of 157 (±8.5) PgC yr −1 , which is consistent with estimates from 18 O (150–175 PgC yr −1 ) and soil respiration ($149^{+29}_{-23}$) PgC yr −1 ), but with an improved confidence level. Our global GPP is higher than satellite optical observation-driven estimates (120–140 PgC yr –1 ) that are used for Earth system model benchmarking. This difference predominantly occurs in the pan-tropical rainforests and is corroborated by ground measurements, suggesting a more productive tropics than satellite-based GPP products indicated. As GPP is a primary determinant of terrestrial carbon sinks and may shape climate trajectories, our findings lay a physiological foundation on which the understanding and prediction of carbon–climate feedbacks can be advanced.

Lai, Jiameng↗
DOE OSTI2026

Hydrogen Infrastructure Modeling and Optimization

HyET, established in 2008, is a small business developing electro-chemical hydrogen compressors. HyET’s goals are for a low-cost, efficient, high-pressure compressor ready for market and installation into hydrogen infrastructure stations servicing mobile fuel cell applications. NLR has been working on hydrogen infrastructure research for over a decade and includes electrolyzer production via electrolysis, storage, compression, dispensing, safety, and performance validation. Accelerating technologies into the marketplace is one goal of NLR’s hydrogen research. Therefore, HyET and NLR have teamed up to provide one year of support for HyET’s modeling and small-scale experiments optimization efforts for HyET’s next-generation electro-chemical compressor.

08 HYDROGEN↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗
DOE OSTI2021

Thermal properties of U-Mo alloys irradiated under high fission power density

A variety of physical and thermal property measurements made on irradiated U-Mo alloy monolithic fuel samples with a Zr diffusion barrier and clad in Al alloy 6061 are reported. These properties form the basis for additional investigations in thermal behavior of U-Mo fuel under high performance research reactor irradiation conditions. Measurements were performed on samples harvested from fuel plates irradiated under high fission power density and surface heat flux and complement previous measurements on samples subjected to moderate fission power density and surface heat flux. Specific heat capacity of irradiated monolithic U-Mo increases with increasing temperature and fission density, but carries large uncertainties associated with thickness measurements used to determine the mass fraction of the layer constituents for extraction of the U-Mo specific heat capacity, specifically the very thin and nonuniform Zr diffusion barrier. Density of irradiated monolithic U-Mo decreases with increasing fission density resulting from fission gas swelling and porosity formation. Thermal diffusivity of irradiated monolithic U-Mo increased with increasing temperature and decreased with increasing fission density. Thermal conductivity of irradiated monolithic U-Mo increased with increasing temperature, and in general, appeared to decrease and become less sensitive to temperature with increasing fission density. Thermal conductivity measurements compared to a semi-empirical model agreed within ±30%. It is suspected that the main driver in deviation between model calculations and measurements is the fission gas swelling determination, which does not account for fission power density or irradiation temperature.

uranium (U), molybdenum, nuclear fuel, thermal con↗
DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
DOE OSTI2022

A Phase-Field Study on Internal to External Oxidation Transition in High-Temperature Structural Alloys

Structural alloys applied at high temperatures rely on an external dense layer of oxide scale for protection. As some reactive alloy components are selectively oxidized internally, understanding how these dispersed metal oxide particles reach the surface to form a protective scale—the so-called internal to external oxidation transition—is crucial for designing these alloys. While the literature is replete with experimental studies on oxidation of alloys, there is a lack of computational studies in this realm due to the complex nature of coupled reaction and diffusion processes in multicomponent multi-phase alloy systems. In this work, we apply a recently developed phase-field model to simulate the oxidation processes under different compositions and nucleation scenarios to gain insights into how a continuous oxide scale can be established. The results show that while alloy composition is critical for internal to external oxidation transition, the oxide nuclei size, shape and distribution also have significant impact on the transition kinetics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

BISON analyses of TRISO fuel performance, its dependence on time-at-temperature, and possible implications for fuel design and qualification

The Advanced Gas Reactor Fuel Development and Qualification (AGR) program has established a substantial technical foundation to support private entry into the U.S. high-temperature gas-cooled reactor market. However, emerging tristructural isotropic (TRISO)-fueled reactor applications include small modular reactors and microreactors with longer fuel residence times, which may expose fuels to higher time-at-temperature (TAT) values than were explored by the AGR program. Increased TAT could affect diffusive and thermomechanical behaviors such as Pd penetration, fission gas release, creep, and fission product transport. In this work, we applied multiscale best-estimate BISON fuel performance modeling to assess these effects within a representative design space based on the AGR-5/6/7 experiment and analyzed trends in predicted particle and compact fuel performance metrics with possible implications for near-term fuel design and qualification. BISON unambiguously predicted that TRISO fuel performance is sensitive to TAT. Increasing TAT was not predicted to increase the magnitude of failure-inducing tangential stresses in particle coating layers. Predictions obtained using a mechanistic model for Pd penetration indicated that penetration depth does not depend strongly on TAT. While these observations suggest that AGR testing provides a conservative upper bound for the steady-state operation of TRISO particles at lower powers and higher residence times, BISON also predicted that the release of poorly retained Ag would increase with TAT. Because these analyses applied models to extrapolate beyond the available experimental data, the authors recommend performing targeted experiments to confirm these predictions. Nevertheless, these predictions may provide reactor developers with enough confidence to make near-term design decisions associated with the potential fuel performance trade-offs of increasing TAT.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Importance of the Antarctic Slope Current in the Southern Ocean Response to Ice Sheet Melt and Wind Stress Change

Abstract We use two coupled climate models, GFDL‐CM4 and GFDL‐ESM4, to investigate the physical response of the Southern Ocean to changes in surface wind stress, Antarctic meltwater, and the combined forcing of the two in a pre‐industrial control simulation. The meltwater cools the ocean surface in all regions except the Weddell Sea, where the wind stress warms the near‐surface layer. The limited sensitivity of the Weddell Sea surface layer to the meltwater is due to the spatial distribution of the meltwater fluxes, regional bathymetry, and large‐scale circulation patterns. The meltwater forcing dominates the Antarctic shelf response and the models yield strikingly different responses along West Antarctica. The disagreement is attributable to the mean‐state representation and meltwater‐driven acceleration of the Antarctic Slope Current (ASC). In CM4, the meltwater is efficiently trapped on the shelf by a well resolved, strong, and accelerating ASC which isolates the West Antarctic shelf from warm offshore waters, leading to strong subsurface cooling. In ESM4, a weaker and diffuse ASC allows more meltwater to escape to the open ocean, the West Antarctic shelf does not become isolated, and instead strong subsurface warming occurs. The CM4 results suggest a possible negative feedback mechanism that acts to limit future melting, while the ESM4 results suggest a possible positive feedback mechanism that acts to accelerate melt. Our results demonstrate the strong influence the ASC has on governing changes along the shelf, highlighting the importance of coupling interactive ice sheet models to ocean models that can resolve these dynamical processes.

Beadling, R. L.↗
DOE OSTI2025

Observation of kinetic mix enhancement in thin-shell OMEGA implosions

Recent separated reactant experiments for thin-shell (6 µ⁢m) shock-driven implosions on OMEGA have demonstrated significant mix from a buried deuterated layer of the shell into the hot spot. Time resolved D 3 He-p reaction history data demonstrate a (50 ± 20)⁢ ps shift earlier in peak nuclear emission for separated reactant experiments relative to control, in contrast to past experimental data for thicker, 20 µ⁢m shells with no laser burn through that show a 75 ps delay due to the time required for hydrodynamic instabilities to develop. This contrast suggests that the mix mechanism was not hydrodynamic. Ion kinetic simulations utilizing fall line analyses show much closer agreement with mix yield and temperature than diffusion models, predicting a D 3 He-p mix yield of 1.7 × 10 9 as compared to the experimental value of 9.3⁢ (±2.1) × 10 8 . This is three orders of magnitude closer than the fall line analysis from a hydrodynamic simulation with an inline diffusive mix model, which suggests minimal mix and D 3 He-p yields of 5×10 5 . This makes kinetic mechanisms the only feasible explanation for the mix seen, demonstrating impact of a non-standard mix mechanism. An analytical model of this kinetic mix mechanism suggests that it can remain significant in situations when the shell expands significantly to low densities, and diffusive models predict negligible mix. Finally, kinetic mix will impact multiple types of high energy density, laser-driven fusion experiments including high-adiabat direct drive cryoexperiments, nuclear cross section experiments, and thin-shell polar direct drive experiments used to tune heat conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Time-lapse approach to correct deficiencies of 2D soil zymography

Membrane zymography is commonly used in rhizosphere ecology to map enzyme activities in intact soil samples and plant roots. The method consists of incubating a membrane saturated with an enzyme-specific fluorogenic substrate on the soil/root surface followed by measurements of fluorescence intensity of the product in the membrane. The traditional zymography is based on assumptions that fluorescence on membrane images is linearly increasing with time during zymography and an increase in product content is numerically equal to enzyme activity in soil below the membrane. These assumptions are unlikely to hold in experimental settings. In this work we introduce a new zymography technique, time-lapse zymography (TLZ); the approach that eliminates the need for assumptions of the traditional zymography and provides more realistic estimates of enzymatic activities. We assessed the new technique in a series of laboratory and modeling experiments, including quantification of the fluorescent product diffusion (e.g. MUF: 4-methylumbelliferone) from plant roots, enzyme activity measurements on the roots, and HYDRUS-2D & HP2 software calibration with obtained data. The calibrated model was used to analyze the processes governing spatial and temporal dynamics of MUF contents in the membrane. The results indicated that the enzyme diffusion within the membrane-soil system was negligible, and measured zymograms were adequately reproduced solely by accounting for substrate and product diffusions and for catalytic enzyme reaction described by the Michaelis-Menten equation. TLZ enabled identifying and using linear parts on MUF time series and accounting for MUF losses in each zymogram pixel, considerably improving the accuracy as compared to the traditional method. Results demonstrated that enzymatic activity from only a thin soil layer (~0.2 mm) is adequately represented in zymograms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Compositional zoning of the Otowi Member of the Bandelier Tuff, Valles caldera, New Mexico, USA

The Otowi Member of the Bandelier Tuff erupted at ca. 1.60 Ma from the Valles caldera (New Mexico, USA). It consists of as much as 400 km 3 (dense rock equivalent) of strongly differentiated high-silica rhyolite and shows systematic upward variations in crystallinity, mineral chemistry, and trace element concentrations through its thickness, but the major element composition is almost constant and is near the low-pressure granite minimum. Incompatible trace elements in whole pumice fragments and glasses show well-correlated linear covariations. Upward zoning to lower abundances of incompatible trace elements is accompanied by development of overgrowths on quartz and alkali feldspar, although earlier-formed interiors of quartz and feldspar have near-constant compositions throughout the tuff, modified by cation diffusion in the case of feldspar. Melt inclusions in remnant quartz cores show diverse Pb isotope ratios, pointing to a wide range of distinct protoliths that contributed rhyolitic melt to the Otowi magma. Mineral thermometers suggest a modest temperature gradient through the melt body, perhaps of 40 °C, at the time of eruption. Chemical, textural, and mineralogical variations and volume-composition relations through the tuff are consistent with an origin for zoning by melting of a high-crystallinity cumulate layer beneath cognate supernatant liquid to produce denser, remobilized liquid of accumulative composition (i.e., the “modified mush model”). Melting may have occurred in several episodes. The latest of these episodes, probably thousands of years prior to eruption, introduced new rhyolitic liquid into the system and was associated with a thermal excursion, recorded in core compositions of pyroxene, during which much of the earlier crystal mass was dissolved. This left inherited cores and interiors of accumulated quartz and feldspar mantled with new growth having less-evolved compositions (higher Ti, Sr, and Ba). Changing solubility of zircon during cumulate melting produced a reversal of Zr concentrations. There is no clear petrologic evidence of a recharge eruption trigger; nonetheless, compositional zoning resulted mainly from repeated recharge-induced remobilization of quartz-feldspar cumulate. The Otowi system was built, evolved, and modified by several events over the course of a few hundred thousand years.

Wolff, J. A.↗
DOE OSTI2025

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↗