SEARCH · Search NASA

Results for “simulation framework”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

DOE OSTI2025

CaptureFOAM

CaptureFOAM provides utilities for fast and accurate direct air capture (DAC) contactor simulations, based on the OpenFOAM framework. At the core, CaptureFOAM is a suite of mass transfer and chemistry models built upon the OpenFOAM thin film model, which allows coarser mesh resolution and less constrictive time step requirements when compared to a traditional volume of fluid approach. This allows CaptureFOAM to conduct device-scale simulations for complex geometries including commercial packing materials and novel contactor shapes, matching experimental results in a wide variety of applications.

Kincaid, Kellis [Oak Ridge National Laboratory (OR↗

DOE OSTI2014

Formation mechanism of the secondary building unit in a chromium terephthalate metal-organic framework

Based on density functional theory calculations and simulation, a detailed mechanism is presented on the formation of the secondary building unit (SBU) of MIL-101, a chromium terephthalate metal-organic framework (MOF). SBU formation is key to MOF nucleation, the rate-limiting step in the formation process of many MOFs. A series of reactions that lead to the formation of the SBU of MIL-101 is proposed in this work. Initial rate-limiting reactions form the metal cluster with three chromium (III) atoms linked to a central bridging oxygen. Terephthalate linkers play a key role as chromium (III) atoms are joined to linker carboxylate groups prior to the placement of the central bridging oxygen. Multiple linker addition reactions, which follow in different paths due to structural isomers, are limited by the removal of water molecules in the first chromium coordination shell. The least energy path is identified were all linkers on one face of the metal center plane are added first. A simple kinetic model based on transition state theory shows the rate of secondary building unit formation similar to the rate metal-organic framework nucleation. The authors are thankful to Dr. R. Rousseau for a critical reading of the manuscript. This research would not have been possible without the support of the Office of Fossil Energy, U.S. Department of Energy. This research was performed using EMSL, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and the PNNL Institutional Computing (PIC) program located at Pacific Northwest National Laboratory.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Agentic traffic intelligence: Augmented human-in-the-loop scenario generation for microscopic traffic simulation

Traditional microscopic traffic simulation generation often relies on static datasets and manual design, limiting its ability to simulate complex conditions easily. This paper presents a novel framework, Agentic Traffic Intelligence, which combines human approval large language models (LLMs), the Real-Twin tool, and multi-agent systems to perform realistic microscopic traffic simulation scenario generation. The proposed framework incorporates human-in-the-loop (HIL) control, retrieval-augmented generation (RAG), and multi-agent control mechanisms. HIL mechanisms are used to guide multiple LLMs focused on attributes for microscopic simulation generation and to improve the interpretability and transparency of LLM execution for users. RAG enhances context extraction by dynamically integrating external knowledge sources for traffic scenario generation foundations. A multi-agent architecture with supervisory control coordinates the interaction of simulation components, including traffic simulators, control logic, and calibration tools. This enables the synthesis of simulation-ready scenarios that reflect dynamic demand profiles and behavior controls. Furthermore, the framework fuses multisource traffic data with unstructured context and supports iterative refinement through interactive user feedback. Validated through microscopic simulation using Simulation of Urban Mobility, the generated scenarios demonstrate high-fidelity network generation with inflow and turn movement and behavioral calibration, offering a robust and efficient tool for stress-testing and optimizing urban mobility systems.

Hierarchical multi-agent control↗

DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Hydrodynamic characterization of the coastal pioneer array ocean observing system

Ocean observation buoys require relatively small amounts of power, yet traditionally necessitate costly resupply trips for battery replacement. With the offshore location of the buoys and small power requirements, wave energy may be an effective solution for providing consistent and reliable power to support the buoy instrumentation. The US National Science Foundation Ocean Observatories Initiative (OOI) includes arrays of point absorber-like buoy systems used for ocean observation that have been deployed at multiple locations including the Southern Mid-Atlantic Bight. A study is currently underway to design a pitch resonator wave energy converter to supplement existing renewable energy generation for powering observation instrumentation. This paper details field measurements from surface moorings of the OOI Coastal Pioneer Array, which informs the subsequent development of a numerical model for the moored observation system. The model is developed in Wave Energy Converter Simulator (WEC-Sim), which leverages the Simscape multibody solver within the MATLAB/Simulink framework and linear potential flow theory to simulate the hydrodynamic interactions and multibody dynamics in 6 degrees of freedom. Multiple tuning variables are considered to produce a model for the system that matches well with empirical data (about 8% error). In conclusion, the WEC-Sim model will serve as a platform for integrating the pitch resonator wave energy converter concept and deployment preparation (detailed design including power take-off and control systems, response evaluation, etc.).

hydrodynamic modeling↗

DOE OSTI2024

Surrogate Model Guided Optimization of Expensive Black-Box Multi-Objective Problems: A Posteriori Methods

Many engineering applications require the simultaneous optimization of multiple conflicting objective functions. Often, these objective functions are evaluated using highly accurate computer simulations that are computationally too expensive to be evaluated hundreds or thousands of times during optimization. Thus, the goal is to find good approximations of the Pareto front using as few of these expensive simulations as possible. Here, we describe an optimization approach based on surrogate models and diverse sampling strategies to accelerate the search for the Pareto solutions. We use a separate surrogate model for approximating each objective function and then we use the surrogate models to inform where additional expensive simulations should be run. The surrogate models are updated in an active learning framework whenever new information from the expensive simulations becomes available. The sampling strategies aim at balancing local improvements of the approximate Pareto front and global exploration to identify the extrema and fill in large gaps of the approximate Pareto front. We demonstrate on a large set of benchmark problems the effectiveness of the method for finding good approximations of the Pareto front.

MATHEMATICS AND COMPUTING↗

DOE OSTI2013

Simulation of background from low-level tritium and radon emanation in the KATRIN spectrometers

The KArlsruhe TRItium Neutrino (KATRIN) experiment is a large-scale experiment for the model independent determination of the mass of electron anti-neutrinos with a sensitivity of 200 meV/c{sup 2}. It investigates the kinematics of electrons from tritium beta decay close to the endpoint of the energy spectrum at 18.6 keV. To achieve a good signal to background ratio at the endpoint, a low background rate below 10{sup −2} counts per second is required. The KATRIN setup thus consists of a high luminosity windowless gaseous tritium source (WGTS), a magnetic electron transport system with differential and cryogenic pumping for tritium retention, and electro-static retarding spectrometers (pre-spectrometer and main spectrometer) for energy analysis, followed by a segmented detector system for counting transmitted beta-electrons. A major source of background comes from magnetically trapped electrons in the main spectrometer (vacuum vessel: 1240 m{sup 3}, 10{sup −11} mbar) produced by nuclear decays in the magnetic flux tube of the spectrometer. Major contributions are expected from short-lived radon isotopes and tritium. Primary electrons, originating from these decays, can be trapped for hours, until having lost almost all their energy through inelastic scattering on residual gas particles. Depending on the initial energy of the primary electron, up to hundreds of low energetic secondary electrons can be produced. Leaving the spectrometer, these electrons will contribute to the background rate. This contribution describes results from simulations for the various background sources. Decays of {sup 219}Rn, emanating from the main vacuum pump, and tritium from the WGTS that reaches the spectrometers are expected to account for most of the background. As a result of the radon alpha decay, electrons are emitted through various processes, such as shake-off, internal conversion and the Auger deexcitations. The corresponding simulations were done using the KASSIOPEIA framework, which has been developed for the KATRIN experiment for low-energy electron tracking, field calculation and detector simulation. The results of the simulations have been used to optimize the design parameters of the vacuum system with regard to radon emanation and tritium pumping, in order to reach the stringent requirements of the neutrino mass measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗

DOE OSTI2013

Risk-Based Comparison of Carbon Capture Technologies

In this paper, we describe an integrated probabilistic risk assessment methodological framework and a decision-support tool suite for implementing systematic comparisons of competing carbon capture technologies. Culminating from a collaborative effort among national laboratories under the Carbon Capture Simulation Initiative (CCSI), the risk assessment framework and the decision-support tool suite encapsulate three interconnected probabilistic modeling and simulation components. The technology readiness level (TRL) assessment component identifies specific scientific and engineering targets required by each readiness level and applies probabilistic estimation techniques to calculate the likelihood of graded as well as nonlinear advancement in technology maturity. The technical risk assessment component focuses on identifying and quantifying risk contributors, especially stochastic distributions for significant risk contributors, performing scenario-based risk analysis, and integrating with carbon capture process model simulations and optimization. The financial risk component estimates the long-term return on investment based on energy retail pricing, production cost, operating and power replacement cost, plan construction and retrofit expenses, and potential tax relief, expressed probabilistically as the net present value distributions over various forecast horizons.

Engel, David W.↗

DOE OSTI2017

An ice sheet model validation framework for the Greenland ice sheet

We propose a new ice sheet model validation framework the Cryospheric Model Comparison Tool (CMCT) that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013 using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quanti- tative metrics for use in evaluating the different model simulations against the observations. We find 10 that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin- and whole-ice-sheet scale metrics, the model initial condition as well as output from idealized and dynamic models all provide an equally reasonable representation of the ice sheet surface (mean elevation differences of <1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CMCT, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate predictive skill with respect to observed dynamic changes occurring on Greenland over the past few decades. An extensible design will allow for continued use of the CMCT as future altimetry, gravimetry, and other remotely sensed data become available for use in ice sheet model validation.

58 GEOSCIENCES↗

DOE OSTI2025

Paradigm for global gyrokinetic turbulence

Using a novel wave-number-advection algorithm, we describe a method to add global profile curvature terms to the local gyrokinetic equations. This approach enables a high-precision global simulation capability without sacrificing the efficiency and spectral accuracy of local simulations. Preliminary numerical results suggest that there are two distinct types of global correction to the local result: (1) a trivial shift of the local result in radius, (2) a non-trivial decorrelation that is roughly equivalent to the well-known E×B shear effect. This interpretation provides a new, intuitive framework for understanding global gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗

DOE OSTI2018

The Effect of Topology in Lewis Pair Functionalized Metal Organic Frameworks on CO 2 Adsorption and Hydrogenation

We have used density functional theory and classical grand canonical Monte Carlo simulations to identify two functionalized metal organic frameworks (MOFs) that have the potential to be used for both CO 2 capture from flue gas and catalytic conversion of CO 2 to valuable chemicals. These new materials based on MIL-140B and MIL-140C functionalized with Lewis pair (acid and base) moieties, which are integrated into the framework linkers. In this work, we show that the Lewis pair functional groups are capable of catalyzing heterolytic dissociation of H 2 and subsequent hydrogenation of CO 2 through concerted 2-H addition. We have examined the effect of pore size and framework topology on the competitive binding of H 2 and CO 2 . We show that the small pore size of functionalized MIL-140B stabilizes the formation of a pre-activated CO 2 species and that this pre-activated CO 2 has a lower overall reaction barrier for hydrogenation of CO 2 to formic acid than a competing pathway in the same material that does not go through a pre-activated complex. We demonstrate that steric hindrance can potentially break energy scaling relationships, which limit the ability to optimize traditional heterogeneous catalysts, by independently changing one part of the CO 2 hydrogenation pathway, without negatively impacting other parts of the pathway. Specifically, we show that steric effects can reduce the CO 2 hydrogenation barrier without impacting the H 2 dissociation barrier or binding energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Sensitivity and Effective Parameterization of a Multi-Scale Model of Proton-Exchange-Membrane Water Electrolysis

The development of proton-exchange-membrane water electrolysis (PEMWE) is crucial to establishing a green-hydrogen infrastructure. Active research in hydrogen technologies is focused on minimizing costs, improving efficiency, and mitigating safety risks. In this study, sensitivity analyses were performed using a physics-based model and used to identify critical parameters capable of reducing overpotentials and safety risks. These parameters include operational (i.e. temperature), structural (i.e. thickness of the components), and material properties (i.e. exchange current density). Simulated results suggest that better engineering of the membrane electrode assemblies (MEA) can offset intrinsic limitations. Membrane thickness and conductivity yielded the largest percent change on performance. Gas crossover was observed to be sensitive to diffusion, solubility, and bubble formation. Current density regimes were identified, where each crossover parameter was most influential. MEA parameterization serves as a framework for theoretical PEMWE optimization simulations and input for design criteria for future research targets.

Dizon, A↗

DOE OSTI2018

A Non-hydrostatic Variable Resolution Atmospheric Model in ACME

The Accelerated Climate Modeling for Energy (ACME) Project is designed to perform coupled climate simulations at high resolution on current and next generation DOE leadership class supercomputers. Prior to this project, the ACME atmospheric model employed the primitive equations of motion, which constrained its maximum resolution and limited the physical phenomena that it could resolve. Our team removed this constraint by implementing a non-hydrostatic (NH) atmospheric model in the ACME Framework. Collaboration with SciDAC FASTMath institutes facilitated improvements to key technologies needed for the non-hydrostatic model. This model was tested in a variable resolution framework, which enables high resolution simulation regions at reasonable computational costs. The correct performance of the NH model was verified with standard dynamical-core tests, and by demonstrating that it reproduced expected climate statistics in the hydrostatic regime. The completed model was applied to evaluate the limitations of existing turbulent mountain stress and gravity wave drag parameterizations, in order to demonstrate the power of the new nonhydrostatic model to model new phenomena and enhance the quality of present day 1/4 degree simulations.

Hall, David Matthew↗

DOE OSTI2023

System and method for simulation-assisted additive manufacturing involving estimating anisotropic deformation

A simulation-assisted additive manufacturing system and method for generating an anisotropic compensation to account for non-uniform deformation due to additive manufacturing and service loading. The predicted deformation may not be fully defined. The present disclosure provides a system and method for estimating the missing deformation data through regression analysis. The present disclosure also provides an integrated framework where the various simulated-assisted design modules are configured for two-way communication and sharing access to changes to the model.

Kim, Seokpum↗

DOE OSTI2024

Simulation of triggering and evolution of ELM by pellet injection in EAST under BOUT++ framework

A BOUT ++ three-field magnetohydrodynamic model is employed to study the triggering and evolution of edge localized mode (ELM) by Li pellets injected along the outer mid-plane in the EAST configuration. The linear simulation shows that compared with a large deposition on the pedestal top (scenario I), a smaller deposition within the steep-gradient pedestal region (scenario II) can stimulate much larger linear growth rates of all-n peeling-ballooning modes (PBMs). The nonlinear simulation shows that there exists a pellet size threshold for ELM triggering for two deposition locations; the threshold for scenario I predicted in the present study matches the EAST observation well. Comparison of the two scenarios reveals that a smaller deposition is sufficient to trigger an ELM in a much shorter time in scenario II, whose ELM size is comparable to that in scenario I. This conclusion confirms previous DIII-D and ASDEX-Upgrade observations, suggesting that the steep-gradient pedestal region is a favorable deposition location for ELM triggering with minimum pellet size. Simulation analyses also find that the positive radial gradient of the hump-like pressure profile in the outer mid-plane induced by the pellet deposition plays a different role in the two scenarios. In scenario I, the force resulting from the gradient hinders the outflow of core plasmas and in return, the perturbation is suppressed from spreading inwards after ELM crashes. In scenario II, with a sizable deposition, the gradient results in another competitive perturbation growth region during the linear phase, thus dispersing the free energy and reducing the efficiency of destabilizing PBMs by pellet injection. The suppressing effect of saturated zonal flow on other modes, the short ELM fast crash phase, and the restricting transport effect of the positive radial pressure gradient work together to constrain the pedestal energy loss, especially when the pellet deposition amount is high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗