Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Understanding Binary Systems—a Comparison between COSMIC and MESA

Abstract We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, as they are promising candidates for the progenitor systems behind energetic, longer lasting, radio bright GRB jets. These systems also may be key contributors to the rates of binary black hole mergers throughout our Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Bootstrap current modeling in M3D-C1

Bootstrap current plays a crucial role in the equilibrium of magnetically confined plasmas, particularly in quasi-symmetric stellarators and in tokamaks, where it can represent bulk of the electric current density. Accurate modeling of this current is essential for understanding the magnetohydrodynamic (MHD) equilibrium and stability of these configurations. This study expands the modeling capabilities of M3D-C1, an extended-MHD code, by implementing self-consistent physics models for bootstrap current. It employs two analytical frameworks: a generalized Sauter model (Sauter et al. 1999 Phys. Plasmas vol. 6, no. 7, pp. 2834–2839), and a revised Sauter-like model (Redl et al. 2021 Phys. Plasmas vol. 28, no. 2, pp. 022502). The isomorphism described by Landreman et al. (2022 Phys. Rev. Lett. vol. 128, pp. 035001) is employed to apply these models to quasi-symmetric stellarators. The implementation in M3D-C1 is benchmarked against neoclassical codes, including NEO, XGCa and SFINCS, showing excellent agreement. These improvements allow M3D-C1 to self-consistently calculate the neoclassical contributions to plasma current in axisymmetric and quasi-symmetric configurations, providing a more accurate representation of the plasma behavior in these configurations. A workflow for evaluating the neoclassical transport using SFINCS with arbitrary toroidal equilibria calculated using M3D-C1 is also presented. This workflow enables a quantitative evaluation of the error in the Sauter-like model in cases that deviate from axi- or quasi-symmetry (e.g. through the development of an MHD instability).

fusion plasma↗
DOE OSTI2021

Separate-Effects Tests for Studying Temperature-Gradient-Driven Cracking in UO 2 Pellets

We report a variety of normal operation and accident scenarios can generate thermal stresses large enough to cause cracking in light-water reactor (LWR) fuel pellets. Cracking of fuel pellets can lead to reduced heat removal, higher centerline temperatures, and localized stress in cladding, all of which impact fuel performance. It is important to experimentally characterize the thermal and mechanical behavior in the pellet before and after cracking to improve cracking models in fuel performance codes. However, in-reactor observation and measurement of cracking is very challenging due to the harsh environment and logistics. Recently, an experimental pellet cracking test stand was developed for separate effects testing of normal operations and accident temperature conditions, using thermal imaging to capture the pellet surface temperature for evaluation of thermal stresses and optical imaging to capture the evolution of cracking in real time. Experiments were performed using depleted uranium dioxide (UO 2 ) pellets, which are useful for collecting data valuable for development and validation of cracking models. A combination of induction and resistance heating was used to create an average temperature gradient of 236°C/cm and 193°C/cm before and after cracking respectively. Characterization of the pellets were done before as well as after cracking. The cracking patterns are somewhat different than those expected in a typical reactor because of the differences in thermal conditions and pellet microstructure. However, if the actual conditions of these experiments are reproduced in computational models, these out-of-pile tests on UO 2 pellets provide relevant data for modeling purposes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Design and modeling of indirectly driven magnetized implosions on the NIF

The use of magnetic fields to improve the performance of hohlraum-driven implosions on the National Ignition Facility (NIF) is discussed. The focus is on magnetically insulated inertial confinement fusion, where the primary field effect is to reduce electron-thermal and alpha-particle loss from the compressed hotspot (magnetic pressure is of secondary importance). We summarize the requirements to achieve this state. The design of recent NIF magnetized hohlraum experiments is presented. These are close to earlier shots in the three-shock, high-adiabat (BigFoot) campaign, subject to the constraints that magnetized NIF targets must be fielded at room-temperature, and use ≲1 MJ of laser energy to avoid the risk of optics damage from stimulated Brillouin scattering. We present results from the original magnetized hohlraum platform, as well as a later variant that gives a higher hotspot temperature. In both platforms, imposed fields (at the capsule center) of up to 28 T increase the fusion yield and hotspot temperature. Integrated radiation-magneto-hydrodynamic modeling with the Lasnex code of these shots is shown, where laser power multipliers and a saturation clamp on cross-beam energy transfer are developed to match the time of peak capsule emission and the P2 Legendre moment of the hotspot x-ray image. The resulting fusion yield and ion temperature agree decently with the measured relative effects of the field, although the absolute simulated yields are higher than the data by 2.0−2.7×. The tuned parameters and yield discrepancy are comparable for experiments with and without an imposed field, indicating the model adequately captures the field effects. Self-generated and imposed fields are added sequentially to simulations of one BigFoot NIF shot to understand how they alter target dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗
DOE OSTI2021

NRAP-Open-IAM Analytical Reservoir Model: Development and Testing

Geological carbon sequestration (GCS) is a key technology for reducing global carbon dioxide (CO 2 ) emissions. Over the last decade, the U.S. Department of Energy has invested in understanding the science base, developing practical implementation methods, and demonstrating secure GCS technologies to mitigate the environmental impacts associated with the atmospheric release of CO 2 . As part of the National Risk Assessment Partnership, a systems-level risk assessment tool, called the NRAP-Open-IAM, has been developed to conduct risk assessment and enable safe operations at a GCS site. The current NRAP-Open-IAM contains a simple reservoir model component that calculates the evolution of CO 2 saturation and fluid pressure in a storage reservoir during CO 2 injection operations. This report presents the development and testing of a new analytical reservoir reduced-order model (ROM), which is extended from an existing semi-analytical model for estimation of CO 2 and brine leakage along legacy wells, and enhances the capability of the NRAP-Open-IAM to simulate more types of reservoir conditions. The developed model is validated against three reference studies, and the results indicate that the new ROM predicts the behavior of the two-phase fluids (brine and injected CO 2 ) well and is applicable to different reservoir simulation boundary conditions (i.e., constant pressure boundary and infinite-acting boundary) without a priori user specification of the boundary type. Sensitivity analysis for a set of model parameters is performed using 4,000 synthetic cases prepared via a fully automated process and using machine-learning-based feature selection. The stochastic analysis identifies gravitational number (i.e., ratio of gravitational forces to viscous force) and distance between the injection well and observation location as the most impactful parameters for matching the pressure and CO 2 saturation, respectively, between the numerical simulations and the ROM. This report details the possible ROM uncertainties and serves as a guide for users to understand the use and limitations of this ROM. The code implementation of the model will be released as a module within the NRAP-Open-IAM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Recent MCNP® Code Developments and Improvements for Nuclear Energy Applications [Slides]

The MCNP® code is a general-purpose radiation transport code developed at LANL over the past 46+ years. This code is capable of modeling the fundamental physics of particles as they move through and interact with materials. It is used in many LANL and worldwide radiation transport applications (see next slide) for example in the Fundamental nuclear physics and data experiments at LANSCE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Characterization of sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements with 3D hybrid particle-in-cell code

Hard x-rays produced by intense laser-produced fast electrons interacting with solids are a vital source for producing radiographs of high-density objects and implosion cores for inertial confinement fusion. Accurate calculation of hard x-ray sources requires a three-dimensional (3D) simulation geometry that fully models the electron transport dynamics, including electron recirculation and the generation of absolute photon yields. To date, 3D simulations of laser-produced bremsstrahlung photons over tens of picoseconds and code benchmarking have not been performed definitively. In this study, we characterize sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements for refluxing and non-refluxing targets using the 3D hybrid particle-in-cell (PIC), Large Scale Plasma code. Bremsstrahlung radiation and escaped electron data were obtained by focusing a 50-TW Leopard laser (15 J, 0.35 ps, 2 × 10 19 W/cm 2 ) on a 100-μm-thick Cu foil and a Cu with a large plastic backing (Cu–CH target). Data for both the Cu and Cu–CH targets were reproduced for simulations with a given set of electron parameters. Comparison of the simulations revealed that the hard x-ray emission from the Cu target was significantly longer in duration than that from the Cu–CH target. Here, the benchmarked hybrid PIC code could prove to be a powerful tool in the design and optimization of time- and angular-dependent bremsstrahlung sources for flash x-ray and gamma-ray radiography.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evaluation of ceria as a surrogate material for UO 2 in experiments on fuel cracking driven by resistive heating

A variety of normal operation and accident scenarios can generate thermal stresses large enough to cause cracking in light-water reactor (LWR) fuel pellets. Cracking of fuel pellets can lead to reduced heat removal, larger centerline temperatures, and localized stress in cladding all of which impact fuel performance. Furthermore, pellet cracking also contributes to a temperature reduction in the pellet since the pellet fragments tend to move towards the heat sink (cladding), and the heat flow remains predominantly radial despite the presence of cracks. It is important to understand the temperature profile on the pellet before and after cracking to improve cracking models in fuel performance codes However, in-reactor observation and measurement of cracking is very challenging owing to the harsh environment and design of fuel rods. Recently, an experimental pellet cracking test stand was developed for separate effects testing of normal operations and accident temperature conditions, using thermal imaging to capture the pellet surface temperature for evaluation of thermal stresses and optical imaging to capture the evolution of cracking in real time. Cracking experiments were initially performed using ceria (CeO 2 ) as a surrogate fuel material, which is useful for developing and demonstrating the experimental approaches but is also valuable in its own right for cracking model development and validation. A combination of induction and resistance heating was used for volumetric heat generation in the pellet creating a thermal gradient. The material properties of CeO 2 and UO 2 are reviewed and compared for use in model development. Simulations of the experiment were performed to evaluate the behavior of the surrogate (CeO 2 ) fuel in BISON. The measured temperature profiles from BISON models match reasonably well with the observed experiments for the ceria pellets before cracking. The findings from this work will help improve confidence in fracture models used for fuel pellets under similar in-reactor conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Algorithmic and GPU enhancements for molecular dynamics in Cabana and LAMMPS

The focus of this milestone was to implement and test new algorithms and optimizations for improved GPU performance by molecular dynamics (MD) and related particle codes. The work was done in the CoPA Cabana library, the CabanaMD mini-app which uses Cabana to implement and benchmark prototypical MD models, and the LAMMPS MD code which is used by the ECP EXAALT project as well as for materials modeling generally. Efforts included (a) implementation of new or optimized algorithms for fundamental MD kernels like neighbor-list building and FFTs which are also widely used in other particle methods, (b) performance testing of data structure layouts for both simple pairwise and computationally intense neural network potentials, (c) optimizations aimed at improving GPU performance for models with small atom counts, and (d) optimizations specific to the SNAP interatomic potential due to its importance for the EXAALT project.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The cis-regulatory codes of response to combined heat and drought stress in Arabidopsis thaliana

Abstract Plants respond to their environment by dynamically modulating gene expression. A powerful approach for understanding how these responses are regulated is to integrate information about cis-regulatory elements (CREs) into models called cis-regulatory codes. Transcriptional response to combined stress is typically not the sum of the responses to the individual stresses. However, cis-regulatory codes underlying combined stress response have not been established. Here we modeled transcriptional response to single and combined heat and drought stress in Arabidopsis thaliana. We grouped genes by their pattern of response (independent, antagonistic and synergistic) and trained machine learning models to predict their response using putative CREs (pCREs) as features (median F-measure = 0.64). We then developed a deep learning approach to integrate additional omics information (sequence conservation, chromatin accessibility and histone modification) into our models, improving performance by 6.2%. While pCREs important for predicting independent and antagonistic responses tended to resemble binding motifs of transcription factors associated with heat and/or drought stress, important synergistic pCREs resembled binding motifs of transcription factors not known to be associated with stress. These findings demonstrate how in silico approaches can improve our understanding of the complex codes regulating response to combined stress and help us identify prime targets for future characterization.

Azodi, Christina B.↗
DOE OSTI2023

Initial Fracture Propagation Modeling of Graphite Components with Grizzly

Graphite has historically been extensively used in power reactor cores and will be used in multiple types of advanced reactors currently under development. These graphite structural components can experience significant stresses due to nonuniform volumetric strains induced by irradiation and thermal expansion, which can lead to fracture. Robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors are important for evaluating component integrity, developing design standards, and interpreting experimental results to characterize graphite performance. The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program has been developing degradation models for other structural components in nuclear reactors within the Grizzly and BlackBear codes. This report documents an effort to develop initial capabilities for modeling graphite fracture within these codes, building on prior efforts to model fracture in other materials. Major elements of this effort include developing a new system for modeling fracture nucleation and growth in two dimensions using the extended finite element method and incorporating a damage and plasticity model. These capabilities are applied here to model a representative graphite component and a splitting disc experiment used to obtain tensile strength.

36 MATERIALS SCIENCE↗
DOE OSTI2026

NGEE Arctic Integrated Modeling (IM3): Improved snow-vegetation interaction

This data product represents the integration of new code capability for arctic tundra snow-vegetation-terrain interactions into the Energy Exascale Earth System Model (E3SM), through the E3SM Land Model (ELM) component. This code integration is the result of collaborative effort between the NGEE Arctic project and the E3SM project. The NGEE Arctic project developed a total of six Integrated Modeling (IM) modules informed by observations and experiments. New ELM capability represented by this data product (IM3) falls into three categories: 1) Downscaling from gridcell to topographic unit level when working through the existing coupler bypass code. 2) Four new parameters (taper, stocking, bendresist, and vegshape) have been added to ELM to allow for flexible definition of snow-vegetation interactions. 3) Vegshape and bendresist parameters are used to calculate the fraction of leaf area and/or stem area buried by snow for a given snow depth. This data record consists of a single document (pdf format) that describes the theoretical basis for the snow-vegetation-terrain interactions added to ELM, and describes the modifications made to the ELM code. The Methods section of this metadata record includes a link to the public E3SM code repository where the exact code modifications as integrated in E3SM can be accessed. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

Thornton, Peter E [ORNL] (ORCID:0000000247595158)↗
DOE OSTI2025

A generative modeling approach to reconstructing 21 cm tomographic data

Abstract Analyses of the cosmic 21 cm signal are hampered by astrophysical foregrounds that are far stronger than the signal itself. These foregrounds, typically confined to a wedge-shaped region in Fourier space, often necessitate the removal of a vast majority of modes, thereby degrading the quality of the data anisotropically. To address this challenge, we introduce a novel deep generative model based on stochastic interpolants to reconstruct the 21 cm data lost to wedge filtering. Our method leverages the non-Gaussian nature of the 21 cm signal to effectively map wedge-filtered 3D lightcones to samples from the conditional distribution of wedge-recovered lightcones. We demonstrate how our method is able to restore spatial information effectively, considering both varying cosmological initial conditions and astrophysics. Furthermore, we discuss a number of future avenues where this approach could be applied in analyses of the 21 cm signal, potentially offering new opportunities to improve our understanding of the Universe during the epochs of cosmic dawn and reionization. Code, pre-trained models, and scripts for making plots in this paper can be found here .

Sabti, Nashwan (ORCID:000000027924546X)↗
DOE OSTI2026

NGEE Arctic Integrated Modeling (IM2): Improved subgrid hillslope hydrologic connectivity

This data product represents the integration of new code capability for arctic tundra hillslope hydrologic processes into the Energy Exascale Earth System Model (E3SM), through the E3SM Land Model (ELM) component. This code integration is the result of collaborative effort between the NGEE Arctic project and the E3SM project. The current ELM represents water movement primarily through vertical processes, such as precipitation, canopy interception, evaporation, infiltration, and soil water movement. Lateral water movement—such as surface runoff, subsurface flow, and river transport—plays a significant role in the hydrological cycle, especially in regions with varied topography. While E3SM includes a runoff routing component representing water transport in the river network, the lateral transport of water at the subgrid scale within the land model has previously not been taken into account. With the recent development of topographic units within the ELM subgrid data structure, there is an opportunity to simulate hillslope hydrologic connectivity by introducing water transport along topographic gradients. We expect that more realistic representation of hillslope hydrologic processes will lead to improved predictions of both soil water content and river network flows. Lateral transport of water at and near the surface is represented as a sub-grid process in this new code development. Water is tracked as it moves from higher to lower elevations within a gridcell. This capability uses the nested hierarchical sub-grid scheme within ELM to connect water fluxes from sub-grid elements with higher elevation to those with lower elevation. This data record consists of a single document (pdf format) that describes the theoretical basis for the hillslope hydrology processes added to ELM, and describes the modifications made to the ELM code. The Methods section of this metadata record includes a link to the public E3SM code repository where the exact code modifications as integrated in E3SM can be accessed. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

Thornton, Peter E [ORNL] (ORCID:0000000247595158)↗