Search NASA⌕ Search

SEARCH · Search NASA

Results for “complex 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗
DOE OSTI2025

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials

The historic December 5, 2022 experiment at Lawrence Livermore National Lab’s (LLNL) National Ignition Facility (NIF) reached fusion energy ignition for the first time. This is the most important scientific breakthrough of the 21st century paves the way to future clean inertial fusion energy (IFE). The diamond (high density carbon (HDC)) ablator material used in this experiment displays detrimental effects due to the development of hydrodynamic instabilities at the diamond/fuel interface under shock compression. New alternatives to diamond ablators are required to step up the energy yield in ICF experiments. The unique combination of mechanical strength (approaching that of diamond), the ability to accommodate high-Z dopants (in contrast to diamond), and the tunability of the properties (through synthesis material with varying sp 3 content) make amorphous carbon (a-C) a promising material for next-generation IFE ablative capsules. However, despite its critical importance to the IFE program, the behavior of a-C carbon at extreme temperatures and pressures remains largely unexplored. The primary goals of this project were to perform groundbreaking dynamic compression experiments and predictive simulations to uncover the fundamental high-energy-density physics of amorphous carbon. Our goals were (1) to uncover the metastability range of amorphous carbon and probe phase transitions to diamond or metastable supercooled liquid carbon; (2) to acquire high-quality equation of state (EOS) data and develop an experimentally validated EOS from machine-learning MD simulations of the complex states of carbon; and (3) to uncover the complex behavior of carbon liquid in both thermodynamically stable and metastable supercooled states by accessing large areas of carbon phase diagram with amorphous samples with variable sp 3 content. Our proposed experimental program included measurements of equation of state and diffraction measurements using the Omega EP laser at the Laboratory of Laser Energetics at the University of Rochester. The theoretical/simulation program involved the development of machine-learning models of the complex response of amorphous carbon under dynamic compression by performing molecular dynamics simulations at experimental time and length scales using leadership class DOE supercomputers. Simulations guided experiments to observe predicted phenomena and acquire critical experimental data in specific pressure-temperature domains to validate theoretical models. This research delivered fundamental properties of novel amorphous carbon IFE ablator material, including phase diagram and EOS. These results will aid in IFE target design and implosion experiments. A unique combination of predictive simulations and dynamic and static experiments provided a highly inspirational intellectual environment for graduate students and postdocs involved in this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

AutonomieAI: An efficient and deployable vehicle energy consumption estimation toolkit

Here, this paper presents AutonomieAI, a novel toolkit designed for efficient energy estimation of vehicles across diverse trip scenarios, routes, and drive cycles, applicable to a broad range of vehicle powertrain technologies. It leverages state-of-the-art Machine Learning techniques to deliver real-time energy prediction of vehicles, enabling co-simulation with transportation level system tools and opening doors for large-scale optimization at city, network or national level. Benchmark results show that AutonomieAI achieves high accuracy, with an average percentage error below 2% for most powertrain types, and computational efficiency capable of processing over 10,000 trips per second. Applications of AutonomieAI have potential to offer the flexibility to assist in solving eco-routing problems, optimize for vehicle and powertrain selection, study charging decision behavior, and optimize for charging station placement. AutonomieAI is the result of large neural network based model architectures, trained on very large and unique high fidelity vehicle simulation data. It is lightweight, deployable, efficient and has accuracy comparable to specialized and complex physics based simulation softwares.

Autonomie↗
DOE OSTI2025

Organic Acid-Assisted Thermal Dehalogenation of Halide Salt Nuclear Wastes: From Waste Salts to Borosilicate Glass

Only a handful of high-halide salt waste forms have been demonstrated for vitrification-based immobilization strategies for halide-salt nuclear waste streams (e.g., pyroprocessing wastes, molten salt reactor wastes) and they all have low waste loading potential and most have low chemical durabilities for high-alkali streams. An alternative approach to direct salt immobilization is salt partitioning prior to waste form fabrication and one option for partitioning is halide removal (called dehalogenation). Removing the halogen fraction through dehalogenation can significantly reduce the waste volume required for disposal in the primary waste form. Furthermore, when dehalogenation is performed using organic acids, the dehalogenation reagent can decompose during high-temperature vitrification, reducing waste loading limitations in the waste form. In the current work, different organic acids (i.e., oxalic, formic, acetic, oxamic, and citric) were evaluated for dehalogenation efficiency of a simple chloride salt simulant (7.19% LaCl 3 , 53.77% LiCl, and 39.04% KCl, by mole) and a more complex chloride salt simulant called ERV3 (electrorefiner version 3) at 150 °C–300 °C and using H + /Cl – molar ratios of 1:1, 2:1, and 3:1. Additionally, a borosilicate glass waste form called TARS (or the average of refined specifications) was formulated, produced, and characterized for dehalogenated ERV3.

Amorphous materials↗
DOE OSTI2024

Developing Data-Driven Synthetic Infrastructure Models for Resilience Analysis

Research on infrastructure resilience has produced promising methods to simulate and optimize complex networks to improve performance. However, restrictions on sharing infrastructure models and the steep cost of developing and maintaining infrastructure models presents a roadblock to adoption. To overcome this limitation, this research focuses on methods to create data-driven infrastructure models that will help improve infrastructure resilience and security. The analysis couples incomplete utility data, geospatial data, machine learning, and synthetic network generation methods to rapidly develop and update infrastructure models. The methods are validated using realistic utility models and site-specific data, with a focus on Puerto Rico due to its unique infrastructure challenges and available data. This research highlights promising opportunities for the use of synthetic network generation and machine learning to create infrastructure models when very little data is available. Results demonstrate that hybrid methods, which combine sparse utility data with synthetic models, can enhance model accuracy, and machine learning can predict model attributes using training data from other models. However, the complexity of infrastructure systems means that even minor changes in network connectivity can significantly impact simulation results. Resilience analysis using synthetic infrastructure models shows that while some system behaviors are preserved, the magnitude of disruptions may not be accurately represented, indicating the need for more research and validation before using synthetic models for critical infrastructure investment decisions. The framework outlined in this report represents a significant advance to infrastructure model development and could be applied to additional domains and sites. Future research will continue to streamline and validate methods to help reduce roadblocks to resilience analysis.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

The effect of composition on helium bubble nucleation in complex concentrated alloys

Molecular dynamics simulations of helium implantation have been performed in the MoNbTaTi complex concentrated alloy (CCA) to determine the impact of composition on helium clustering and bubble nucleation. It was found that compositions with lower interstitial atomic volume resulted in smaller cluster sizes, higher helium-to-vacancy ratios, and larger migration barriers. As the atomic volume decreases, it is harder to accommodate additional helium in the lattice and for helium to diffuse between interstitial sites, resulting in smaller cluster sizes. In particular, increasing molybdenum directly correlated with a decreased interstitial atomic volume and increased migration barrier. Niobium exhibited the opposite trend, where more niobium in the material resulted in larger interstitial atomic volumes and cluster sizes but lower migration barriers. The atomic volume was determined to be an indicator of the susceptibility of the material to helium bubble growth and could potentially be used as a metric to screen compositions for resistance to helium damage. In this way, properties such as available interstitial volumes could be used to screen materials more rapidly given the vast compositional space of CCAs.

McCarthy, Megan Jeanne [Sandia National Laboratori↗
DOE OSTI2025

Prelude to Simulation Diagnostics in MOOSE: Diagnosing Meshes

As finite element simulations increase in complexity, so does the need for high quality meshes. Without proper diagnostics, issues pertaining to poor mesh quality or improper mesh construction may go undetected, causing simulations to fail unexpectedly. This paper details the new mesh diagnostic system available in MOOSE, which provides users with the tools to inspect mesh quality for both externally and internally generated meshes. This system can detect issues arising from element volume size, non-conformal elements, intersecting edges, missing sidesets, as well as other common mesh issues. As meshes increase in size and complexity, such issues become more difficult to diagnose through visual inspection alone.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

SANDWake3D

SAND2026-23958O SANDWake3D models atmospheric surface layers and turbine wakes using a 3D parabolic RANS solver. This Python-based tool uses an efficient method to simulate wakes, capturing complex wake behaviors such as wind shear, veer, and atmospheric stratification, naturally. It has been calibrated and validated against higher-fidelity large eddy simulations (LES), providing reliable results for wind energy and atmospheric research applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2024

The Effect of Air Separations on Fast Pyrolysis Products for Forest Residue Feedstocks

This study investigates the intricate relationship between biomass preprocessing and pyrolysis product yields, employing the air classification technique for the treatment of loblolly pine residues with varying moisture content. A comprehensive exploration of the physicochemical properties of air-classified loblolly pine informs a sophisticated pyrolysis simulation model. Given the complex and multifaceted nature of biomass pyrolysis, operating across diverse temporal and spatial scales, a pyrolysis kinetics-based CFD–DEM simulation method is employed to predict product yields. Results showed that the elevated moisture content amplifies particle adhesiveness, necessitating augmented air velocities for effective separation, thereby influencing the efficiency of the separation process. While carbon and hydrogen contents exhibit relative stability across diverse moisture contents and blower frequencies, the oxygen content undergoes noticeable changes. For example, the oxygen contents were measured as 29.2 and 38.6 wt% in the light fraction of 30% moisture content sample at blower frequencies of 10 and 20 Hz, respectively. An intriguing finding emerges from pyrolysis simulation, indicating that a lower blower frequency in air classification moderately enhances bio-oil yield and significantly improves its quality, particularly in terms of water content. For instance, the water content in the bio-oil was about 1.5% and 10% in the heavy and light fractions, respectively from 10% moisture sample under 15 Hz blower frequency. In summary, a detailed understanding and strategic manipulation of critical material attributes in biomass through efficient fractionation techniques are imperative for advancing fast pyrolysis as a sustainable avenue for renewable energy and chemical production.

09 BIOMASS FUELS↗
DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Lowering entry barriers to developing custom simulators of distributed applications and platforms with SimGrid

Researchers in parallel and distributed computing (PDC) often resort to simulation because experiments conducted using a simulator can be for arbitrary experimental scenarios, are less resource-, labor-, and time-consuming than their real-world counterparts, and are perfectly repeatable and observable. Many frameworks have been developed to ease the development of PDC simulators, and these frameworks provide different levels of accuracy, scalability, versatility, extensibility, and usability. Further, the SimGrid framework has been used by many PDC researchers to produce a wide range of simulators for over two decades. Its popularity is due to a large emphasis placed on accuracy, scalability, and versatility, and is in spite of shortcomings in terms of extensibility and usability. Although SimGrid provides sensible simulation models for the common case, it was difficult for users to extend these models to meet domain-specific needs. Furthermore, SimGrid only provided relatively low-level simulation abstractions, making the implementation of a simulator of a complex system a labor-intensive undertaking. In this work we describe developments in the last decade that have contributed to vastly improving extensibility and usability, thus lowering or removing entry barriers for users to develop custom SimGrid simulators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Revolutionizing Materials Design: The Intersection of Quantum Mechanics and Data Modeling

The field of materials design is currently experiencing a notable evolution, driven by the convergence of sophisticated computational methodologies based on first principles and data-driven modeling approaches. I will review our recent endeavors employing AI/ML to expedite first-principles simulations and mitigate traditional methods' temporal and spatial limitations. Central to our efforts is developing and utilizing ML interatomic potentials (MLPs) across a diverse spectrum of materials. We show that MLPs serve as invaluable tools for navigating the complexities of the simulations, such as understanding the behavior of MgO at extreme environments of ~1 terapascal and temperatures >10,000 Kelvin. Moreover, we show that MLPs can provide precise details of the intricate dynamics governing the oxidation processes of binary alloy systems due to the competition between surface segregation and reconstruction tendencies. In summation, advancements in MLPs open the door to fresh possibilities in material modeling and, ultimately, discovery.

Saidi, Wissam↗
DOE OSTI2025

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗
DOE OSTI2025

2024 Incubator Final Report: Active Hybrid Mooring

This Incubator project has defined and developed a representation of an active hybrid mooring (AHM) system in Simulink to assess whether the speed and accuracy of AHM is sufficient to warrant further investigation in a larger project. We first built a representative mooring system for the VolturnUS-S 15MW floating offshore wind turbine (FOWT) in 1,000-meter water depth to provide baseline loads and displacements needed in the AHM system, with the results showing off-the-shelf actuators would be practical. We then defined model configurations for the numerical and physical substructures, using standalone MoorDyn v2 with a novel input strategy and chained hydraulic actuators represented in Simulink Simscape, respectively. The remainder of the project has focused on verifying the existing models for speed and accuracy against OpenFAST simulations, increasing their complexity to enable six degree of freedom (6DOF) operation, and connecting the numerical and physical substructures into a single holistic model. The verification phase of the project has shown promising results, though further refinement is needed to resolve accuracy discrepancies and for robust performance across a wider range of metocean conditions. The numerical model we developed is computationally stable and operates faster than real time using a time step of 0.01 seconds or shorter, indicating the model is likely fast enough for our AHM design to operate successfully. Additionally, the Simscape model of a hydraulic actuator operates with less than 0.3% relative error when actuating mooring loads from an OpenFAST simulation of the VolturnUS-S 1,000-meter mooring system. The aforementined refinements will be performed in the next year as a Continuing Incubator, with key remaining tasks include programming the predictor-corrector loop to minimize movement error, testing against more complex nonlinear wave cases, and expanding to enable loads from connected shared mooring lines.

17 WIND ENERGY↗
DOE OSTI2025

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗
DOE OSTI2025

Building confidence in models for complex barrier systems for radionuclides

The modeling and simulation of the Cement-clay Interaction-Diffusion field (CI-D) experiment at the Mont Terri site in Switzerland presented here demonstrates that it is possible to capture the multiscale physical and chemical features of natural and engineered barrier systems for radionuclides. The simulations are successfully carried out with the newly developed CrunchODiTi high-performance computing software that accounts for multiple continua, including a continuum representing the electrical double layer (EDL) developed along negatively charged clay particles in clay rock. The simulation also accounts for both the complex three-dimensional (3D) geometry, expected as the norm in a geological waste repository, and the anisotropy of the geological formation. In addition, the high resolution of the model makes it possible to include "skin effects" developed at the interface between highly reactive materials, in this case between the high pH cement and the circumneutral but electrostatic Opalinus Clay. The successful history matching with the field experiment demonstrates that the distinct geochemical and physical properties of the cement and the Opalinus Clay in the CI-D experiment can be accounted for. Such analyses are essential for developing a defensible safety case for the underground storage of radioactive waste.

Sarsenbayev, Dauren↗