Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive simulation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Axion Mass Prediction from Adaptive Mesh Refinement Cosmological Lattice Simulations

The quantum chromodynamics (QCD) axion arises as the pseudo-Goldstone mode of a spontaneously broken Abelian Peccei-Quinn (PQ) symmetry. If the scale of PQ symmetry breaking occurs below the inflationary reheat temperature and the domain wall number is unity, then there is a unique axion mass that gives the observed dark matter (DM) abundance. Computing this mass has been the subject of intensive numerical simulations for decades since the mass prediction informs laboratory experiments. Axion strings develop below the PQ symmetry-breaking temperature, and as the string network evolves, it emits axions that go on to become the DM. A key ingredient in the axion mass prediction is the spectral index of axion radiation emitted by the axion strings. We compute this index in this Letter using the most precise and accurate large-scale simulations to date of the axion-string network leveraging adaptive mesh refinement to achieve the precision that would, otherwise, require a static lattice with 262,144 3 lattice sites. We find a scale-invariant axion radiation spectrum to within 1% precision and find no evidence that the spectral index of radiation evolves with time. Accounting for axion production from strings prior to the QCD phase transition leads us to predict that the axion mass should be approximately 𝑚 𝑎 ∈ (45, 65) μ ⁢eV. However, we provide preliminary evidence that axions are produced in greater quantities from the string-domain-wall network collapse during the QCD phase transition, potentially increasing the mass prediction to as much as 300 μ ⁢eV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

When more data hurts: Optimizing data coverage while mitigating diversity-induced underfitting in an ultrafast machine-learned potential

Machine-learned interatomic potentials (MLIPs) are becoming an essential tool in materials modeling. However, optimizing the generation of training data used to parametrize the MLIPs remains a significant challenge. This is because MLIPs can fail when encountering local environments too different from those present in the training data. The difficulty of determining a priori the environments that will be encountered during molecular dynamics simulation necessitates diverse, high-quality training data. Here, this study investigates how training data diversity affects the performance of MLIPs using the Ultra-Fast force field (UF 3 ) to model amorphous silicon nitride. We employ expert and autonomously generated data to create the training data and fit four force field variants to subsets of the data. Our findings reveal a critical balance in training data diversity: insufficient diversity hinders generalization, while excessive diversity can exceed the MLIP's learning capacity, reducing simulation accuracy. Specifically, we found that the UF 3 variant trained on a subset of the training data, in which nitrogen-rich structures were removed, offered vastly better prediction and simulation accuracy than any other variant. By comparing these UF 3 variants, we highlight the nuanced requirements for creating accurate MLIPs, emphasizing the importance of application-specific training data to achieve optimal performance in modeling complex material behaviors.

ab initio molecular dynamics↗
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↗
DOE OSTI2024

Neural network kinetics for exploring diffusion multiplicity and chemical ordering in compositionally complex materials

Diffusion involving atom transport from one location to another governs many important processes and behaviors such as precipitation and phase nucleation. The inherent chemical complexity in compositionally complex materials poses challenges for modeling atomic diffusion and the resulting formation of chemically ordered structures. Here, we introduce a neural network kinetics (NNK) scheme that predicts and simulates diffusion-induced chemical and structural evolution in complex concentrated chemical environments. The framework is grounded on efficient on-lattice structure and chemistry representation combined with artificial neural networks, enabling precise prediction of all path-dependent migration barriers and individual atom jumps. To demonstrate the method, we study the temperature-dependent local chemical ordering in a refractory NbMoTa alloy and reveal a critical temperature at which the B2 order reaches a maximum. The atomic jump randomness map exhibits the highest diffusion heterogeneity (multiplicity) in the vicinity of this characteristic temperature, which is closely related to chemical ordering and B2 structure formation. The scalable NNK framework provides a promising new avenue to exploring diffusion-related properties in the vast compositional space within which extraordinary properties are hidden.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Origin of Heating-Induced Softening and Enthalpic Reinforcement in Elastomeric Nanocomposites

Molecular simulations demonstrate that the enthalpic softening of elastomeric nanocomposites upon heating can arise naturally from a Poisson’s ratio mismatch between elastomer and nanoparticle networks, providing a more parsimonious explanation for this phenomenon than the widely accepted interpretation based on glassy interparticle bridging. Despite a century of use, the mechanism of nanoparticle-driven mechanical reinforcement of elastomers is unresolved. Here, a major hypothesis attributes it to glassy interparticle bridges, supported by an observed inversion of the variation of the modulus E(T) on heating – from entropic stiffening in elastomers to enthalpic softening in nanocomposites. Here, molecular simulations reveal that elastomer enthalpic softening can instead emerge from a competition over the preferred volumes between elastomer and nanoparticulate networks. A theory for this competition accounting for softening of the bulk modulus on heating predicts the simulated E(T) inversion, suggesting that reinforcement is driven by a volume-competition mechanism unique to cocontinuous systems of soft and rigid networks.

Biopolymers↗
DOE OSTI2025

Monte Carlo N-Particle Transport Performance of Predicting Digital Radiographic IQI Inspection

The identification of porosity, geometric noncompliance, and other defect types are critical to the qualification of materials and components. X-ray radiographic nondestructive testing is a common industrial inspection method for process quality control and component qualification and certification. Digital radiography provides a quick and efficient alternative when compared to traditional film-based inspection. The quality of radiographic inspection is dependent on equipment specifications, such as the source spot size and detector pixel size, and the specific parameters selected for use for the radiographic technique. To evaluate if an x-ray system and technique is sufficient for a given requirement, a radiographic image quality indicator (IQI) can be used. Radiographic IQIs in hard to machine materials or hard to manufacture defects can be time consuming and expensive to manufacture. This study was conducted to evaluate current Savannah River National Laboratory (SRNL) x-ray imaging systems with a custom tantalum IQI and using Monte Carlo simulations to predict the performance of future systems. The tantalum IQI was tested using a Siefert Isovolt 420 keV x-ray tube with a Perkin Elmer XRD 1611 flat panel with 100-micron pixels. Using the Monte Carlo N-Particle transport software, the radiographic tally was used to simulate the photon flux through an identical tantalum IQI. These simulations provided a benchmark as to the best theoretical identification on a given system using our tantalum IQI. The simulations were refined to match SRNL’s current systems’ noise levels, leading to confidence in their ability to predict the performance of other systems that may be purchased and deployed in the future at the Savannah River Site. Future studies will be conducted to prove this research can be extended to artificially evaluate the ability for systems to identify critical defect sizes through x-ray radiographic inspection, drastically reducing the cost and time burdens of producing high-fidelity radiographic test articles.

digital X-ray radiography↗
DOE OSTI2025

High-Fidelity CFD Modeling of Cryogenic Hydrogen Isotope Extrusion for Fusion Reactor Pellet Fueling

This study investigates the extrusion processes of deuterium and protium using ANSYS-Polyflow. The geometries and computational fluid dynamics (CFD) settings closely replicate the experimental setups and data acquired from the extruder experiments at Oak Ridge National Laboratory (ORNL) for validation purposes. We explore the impacts of (1) slip versus non-slip boundary conditions and (2) the use of constant, temperature-, and shear rate–dependent viscosities, concluding that the implementation of non-slip wall boundary conditions combined with shear rate–dependent viscosity produced more accurate predictions. The simulations achieved excellent agreement with the experimental data, with relative differences of only 5% for deuterium, and 3% to 6% for protium. This is the first time that experimental extrusion data at ORNL have been accurately predicted through high-fidelity CFD modeling. In conclusion, the advancements offer valuable insights and a foundational modeling tool for optimizing pellet injectors for ITER and other future reactor-scale devices.

ANSYS-Polyflow↗
DOE OSTI2024

Optimization of ray-tracing simulations to confirm performance of the GP-SANS instrument at the High-Flux Isotope Reactor

The CG-2 beamline at the High Flux Isotope Reactor (HFIR) exhibits a notable discrepancy between observed count rates and the count rates we would expect based on a Monte-Carlo neutron ray-trace simulation. These simulations consistently predict count rates approximately five times greater than those observed in four separate experimental runs involving different instrument configurations. This discrepancy suggests that certain factors are causing losses in measurements that are not adequately accounted for in the simulation, in particular guide reflectivity or misalignment. To investigate these discrepancies, a high-dimensional simulation parameter approach is applied in order to understand the losses. Region of Interest (ROI) groups along the instrument are assigned to different surfaces of the guide components within the simulation. This allows the parameters of those guide components to be varied as a group to minimize the complexity of the search space. The result is an optimization of simulation parameters using an iterative scheme that aims to minimize the difference between experimentally measured count rates and simulated count rates across all tested collimator combinations. This proposed methodology holds the potential to reveal previously unrecognized sources of intensity loss in the CG-2 beamline at HFIR and improve the accuracy of simulations, leading to enhanced understanding and performance of the beamline for various scientific applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Computational and experimental analysis of structural and Thermophysical properties of LiX-KX (X = chloride, iodide, or bromide) molten salts

Molten salts containing lithium (Li) and potassium (K) halides, including chlorides (Cl), bromides (Br), and iodides (I), are pivotal in advanced technological applications including nuclear reactor technologies, thermal batteries, and industrial pyrochemical processes. Despite their significance, the availability of the temperature-dependent structural and thermophysical properties of these molten salts remain scarce. Here, this study addresses this knowledge gap by predicting and measuring the key properties such as coordination number, radial distribution function, density, heat capacity, and volumetric thermal expansion of molten LiI-KI, LiBr-KBr and LiCl-KCl eutectic mixtures. Utilizing ab initio molecular dynamics (AIMD) simulations, we predict these properties of LiX-KX (X=Br,I, Cl) across various temperatures, with chlorides as a benchmark. Experimental measurements using Archimedes methods and Differential Scanning Calorimetry validate the densities and heat capacities of these molten salts. Our comprehensive analysis of the structural and thermophysical properties provides critical insights into the behavior of these molten salts at various temperatures, enhancing the understanding necessary for their application in advanced technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Stellarator Design Exploration Using Symbolic-Regression Neutronics Surrogates

Systems codes require fast, simplified models to rapidly evaluate fusion power plant concepts, but neutronics analyses are often a computational bottleneck. Here, to address this, surrogate models for key neutronics responses have been developed using 3-D neutronics-ready models built with the open-source code ParaStell from a database of stellarator equilibria. Neutronics responses such as tritium breeding ratio (TBR), nuclear heating, and neutron-induced radiation damage displacements per atom (dpa) were simulated using OpenMC. Through sensitivity analysis and symbolic regression (SR), simple power-law formulas were derived connecting these neutronics responses to global stellarator parameters, including fusion power, plasma surface area, and plasma elongation. Validation shows these formulas can predict the simulation results with low error, enabling quick and accurate assessment of neutronics requirements in stellarator design exploration activities with systems codes.

Modeling↗
DOE OSTI2025

Perturbative model for the saturation of energetic-particle-driven modes limited by self-generated zonal modes

We present a simplified energy-conserving approach to incorporate wave–wave nonlinear effects within the framework commonly used to describe wave–particle nonlinearities. In particular, the effects of zonal mode (ZM) generation on the determination of the saturation amplitude of energetic particle (EP)-driven Alfvénic instabilities is studied. The model assumes that the zonal perturbations grow at a rate twice that of the original (pump) wave, consistent with a beat-driven (or force-driven) generation mechanism. The evolution and saturation of the mode amplitude are investigated both analytically and numerically within our reduced model assumptions, in both the collisionless and scattering-dominated regimes. These studies underscore the crucial role of sources and sinks in capturing the impact and the role of beat-driven zonal perturbations on mode evolution. In the realistic case of saturation set by sources and sinks, we discuss the role of a finite amplitude ZM in reducing microturbulent particle scattering, thus limiting the energy source for the resonant mode. We then discuss comparisons between the model’s predictions and simulation results. The model reproduces key features observed in gyrokinetic simulations as the reduction in saturated mode amplitude and the onset of wave–wave nonlinear effects as functions of mode growth rate and amplitude. Thanks to its simplicity, it can be readily implemented into codes based on reduced models, thereby improving their predictive capability for strongly driven instabilities.

Alfvén eigenmodes↗
DOE OSTI2025

Structural uncertainty assessment for low-Mach wall-resolved large-eddy simulations: Plane Channel and Periodic Hill Use Cases

A structural uncertainty assessment is provided for a set of wall-resolved large-eddy simulations (WRLES) using a turbulent channel flow (𝑅𝑒 𝜏 395) and a periodic hill in the turbulent regime (𝑅𝑒 𝑏 2800) validation suite. The standard one-equation subgrid scale turbulent kinetic energy model (𝑘 𝑠𝑔𝑠 ) using both a static and dynamic coefficient approach along with Wall Adapting Local Eddy Viscosity (WALE) are exercised. The 𝑘 𝑠𝑔𝑠 model activated in this study also allows for a correction that supports theoretical turbulent viscosity and total kinetic energy scaling (cubit and quadratic, respectively) in the near-wall regime. Channel and hill results showcase very good agreement between prediction and simulation for the WALE and the dynamic coefficient 𝑘 𝑠𝑔𝑠 model. The study also demonstrates that near-wall damping improves predictivity for the static coefficient 𝑘 𝑠𝑔𝑠 model, while proving an unnecessary additional contribution for the dynamic coefficient 𝑘 𝑠𝑔𝑠 model as near-wall scaling is naturally captured. Overall conclusions are that this suite of WRLES models are well suited for production usage in the Sierra-based low-Mach codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Magnon-Magnon Interaction Induced by Dynamic Coupling in a Hybrid Magnonic Crystal

We report a combined experimental and numerical investigation of spin-wave dynamics in a hybrid magnonic crystal consisting of a CoFeB artificial spin ice (ASI) of stadium-shaped nanoelements patterned atop a continuous NiFe film separated by a 5 nm Al 2 O 3 spacer. Using Brillouin light scattering spectroscopy, we probe the frequency dependence of thermal spin waves as functions of applied magnetic field and wavevector, revealing the decisive role of interlayer dipolar coupling in the magnetization dynamics. Micromagnetic simulations complement the experiments, showing a strong interplay between ASI edge modes and backward volume modes in the NiFe film. The contrast in saturation magnetization between CoFeB and NiFe enhances this coupling, leading to a pronounced hybridization manifested as a triplet of peaks in the BLS spectra predicted by simulations and observed experimentally. This magnon−magnon coupling persists over a wide magnetic field range, shaping both the spin-wave dispersion at fixed fields and the full frequency-field response throughout the magnetic hysteresis loop. Our findings establish how ASI geometry can selectively enhance specific spin-wave wavelengths in the underlying film, thereby boosting their amplitude and identifying them as preferential channels for spin wave transmission and manipulation.

Brillouin light scattering spectroscopy↗
DOE OSTI2025

Machine learning interatomic potential for predicting the thermal properties of uranium nitride

We present a combined computational and experimental investigation of the thermal properties of uranium nitride (UN), focusing on the development of a machine learning interatomic potential (MLIP) using the moment tensor potential framework. The MLIP was trained on density functional theory (DFT) data and validated against various quantities including energies, forces, elastic constants, phonon dispersion, and defect formation energies, achieving excellent agreement with DFT calculations, prior experimental results, and our thermal conductivity measurement. The potential was then employed in molecular dynamics simulations to predict key thermal properties such as melting point, thermal expansion, specific heat, and lattice thermal conductivity. To further assess model accuracy, we fabricated a UN sample and performed new thermal conductivity measurements representative of single-crystal properties, which showed strong agreement with the MLIP predictions. This work confirms the reliability and predictive capability of the developed potential for determining the thermal properties of UN.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Establishment of a Vertically Integrated Domestic Manufacturing Process for Production of Substrates Needed for Manufacture of Gas Diffusion Layers

In this project, AvCarb, LLC evaluated the baseline performance metrics of commercial carbon veils and their corresponding Gas Diffusion Layers (GDLs) with the goal of establishing an optimized, vertically integrated production system for wet-laid nonwoven substrates used in gas diffusion media for electrochemical energy storage and conversion devices. Mechanical testing and microstructural characterization were conducted and used to develop a multiscale computational model capable of simulating and predicting the performance of GDLs in fuel cells. Although the project successfully generated foundational transport and modeling data, it was terminated prior to identifying the critical GDL design parameters necessary for full optimization. The program aimed to improve carbon veil fabrication through enhanced fiber dispersion, fiber-fiber adhesion control, and improved web formation, enabling the production of high-quality, uniform substrates. Simulations were intended to guide mixing and solution delivery system design and process conditions, followed by production-scale trials to evaluate fiber dispersion, web uniformity, and mechanical robustness. At full deployment, the proposed production line would have been capable of producing approximately 650,000 m² of carbon veil annually. This capability remains strategically important, as the United States currently lacks a domestic source of wet-laid nonwoven carbon substrates that satisfy the stringent quality requirements for fuel-cell GDLs and electrolyzers representing an ongoing supply-chain vulnerability. Beyond supply-chain benefits, the project established a robust benchmarking dataset for existing commercial carbon veils while advancing next-generation material concepts targeting improved performance and manufacturing consistency.

Olson, Cynthia Lemay↗
DOE OSTI2024

An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations

We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the Illustris-TNG suite. We create an HOD mock with Illustris-like correlations using our method, constraining the associated IA model parameters, with the $\mathcal{X}$$^{2}_{dof}$ between our model’s correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation–position and orientation–orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to 0.1 $\mathcal{h}$ –1 Mpc. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Data assimilation for combustion ignition delay time simulation using schlieren image velocimetry

This study sought to improve the accuracy of simulating spray penetration and combustion ignition delay by means of data assimilation (DA). The simulations were conducted using the Reynolds-averaged Navier–Stokes (RANS) equations and assimilating the schlieren image data. In DA, an ensemble square root filter (EnSRF) was used to build the statistical model, making the simulation results more accurate without any change in the governing equations. Recognizing that the spray-cone injection angle has a large effect on penetration, we created ensemble members with different injection angles. And we applied the two-component velocity distribution calculated via SIV and updated both velocity and temperature by using a DA statistical model derived from RANS ensemble simulations. The ignition delay time is generally known to vary even under the same experimental conditions because it is influenced by many factors. In this study, we attempted the transient DA-assisted RANS simulation to predict the ignition delay time even when the temporal resolution and accuracy of the observation data ware insufficient. Our trials offer an example of how a combination of techniques can be effectively used to assimilate experimental data obtained under restricted conditions.

Combustion simulation↗