Search NASA⌕ Search

SEARCH · Search NASA

Results for “systems 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 OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Determining Levels of Detail for Simulators of Parallel and Distributed Computing Systems via Automated Calibration

There are two sources of inaccuracy when simulating parallel and distributed computing systems: (i) a simulator implemented at an insufficient level of detail; and (ii) incorrectly calibrated simulation parameter values. Increasing the simulator’s level of detail can improve accuracy, but at the cost of higher space, time, and/or software complexity. Furthermore, evaluating the intrinsic accuracy of a simulator requires that its parameters be well-calibrated. Making decisions regarding the level of detail is thus challenging. We propose a methodology for instantiating the simulation calibration process and a framework for automating this process, which makes it possible to pick appropriate levels of detail for any simulator. We demonstrate the usefulness of our approach via two case studies for two different domains.

McDonald, Jessie [University of Hawaii at Manoa, H↗
DOE OSTI2025

Quantum simulation of charge and exciton transfer in multi-mode models using engineered reservoirs

Quantum simulation enables studies of open-system dynamics in non-perturbative regimes by programming electronic, vibrational, and environmental interactions on comparable energy scales. Trapped ions offer this capability, combining spins, phonons, and tunable dissipation on one platform. We demonstrate an open-system quantum simulation of charge and exciton transfer in a multi-mode linear vibronic coupling model. Using tailored spin-phonon interactions with reservoir engineering, we emulate a system with two dissipative vibrational modes coupled to donor and acceptor sites and track its non-equilibrium dynamics. We continuously tune the system from the charge transfer regime to the vibrationally assisted exciton transfer regime and find that degenerate modes enhance transfer rates at large energy gaps, while non-degenerate modes activate pathways that reduce the energy-gap dependence. Thus, the presence of one additional vibration introduces interfering pathways and reshapes non-perturbative excitation transfer. Our results establish a scalable, hardware-efficient route to simulate vibronic processes with engineered environments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Solar Thermochemical Carbon Dioxide Splitting Using Ceria and Iron Aluminate Foam Devices and Simulation of a Plant System for Demonstration

An international research project has been undertaken to integrate a unique solar thermal processing reactor system with ceria and iron aluminate as active redox materials for CO2 splitting. Experimental investigations for CO2 splitting were conducted using a solar simulator and tube furnace at Niigata University, followed by demonstrations using a high-flux solar furnace (HFSF) at the National Renewable Energy Laboratory (NREL) in Golden, CO. Each experimental setup consisted of foam devices composed of reticulated porous ceramic (RPC). The RPC has a full ceria or iron aluminate body. It fabricated using the replica method and subjected to a two-step redox reaction, which iteratively separated a stream of CO2 into O2 and CO. Reactivity was evaluated using CO production per mass of the reactive material. The tubular furnace yielded a CO production of 6.41 mL/g at a reduction temperature of 1600degrees C, showcasing a higher CO production rate and total amount than those obtained from experiments conducted with solar simulators and solar furnace setups. For iron aluminate RPC, the productivity was measured as 3.57 mL/g using HFSF at a reduction temperature of 1450degrees C. These results are somewhat higher than those of the previous experiment at lower reduction temperatures of 1400degrees C-1500degrees C. Additionally, the production of CO in the case of ceria RPC was compared with the steady flow model simulation, which assumed chemical equilibrium at various levels of oxygen partial pressure during the reduction process. On the basis of these results, this study proposes a solar fuel system with an open receiver that uses a high-temperature heat transfer fluid.

carbon dioxide thermochemical splitting↗
DOE OSTI2025

Surrogate-constructed scalable-circuits adaptive variational quantum eigensolver in the Schwinger model

Inspired by recent advancements in simulating periodic systems on quantum computers, we develop an approach to further advance the simulation of these systems, named (SC) 2 -ADAPT-VQE. Our approach extends the scalable-circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable “surrogate constructed” method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC) 2 methodology. Our approach allows for a wider set of classical computations on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute the properties of the Schwinger model—quantum electrodynamics for a single, massive fermion in 1 +1 dimensions—and show that our method can be used to accurately extrapolate to the continuum limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dependencies of Simulated Convective Cell and System Growth Biases on Atmospheric Instability and Model Resolution

Abstract This study evaluates convective cell properties and their relationships with convective and stratiform rainfall within a season‐long convection‐permitting weather research and forecasting simulation over central Argentina using radar, satellite, and radiosonde measurements from the RELAMPAGO‐CACTI field campaign. The simulation slightly underestimates radar‐estimated rainfall over the ∼3.5‐month evaluation period but underestimates stratiform rainfall by 46% and overestimates convective rainfall by 43%. As convective available potential energy (CAPE) increases, the convective rainfall overestimation decreases, but the stratiform rainfall underestimation increases such that the contribution of convective to total rainfall remains constantly high biased by ∼26%. Overestimated convective rainfall arises from the simulation generating 2.6 times more precipitating convective cells (14,299) than observed by radar (5,662) despite similar observed and simulated cell growth processes, with relatively wide cells contributing mostly to excessive convective rainfall. Relatively shallow cells, typically reaching heights of 4–7 km, contribute most to the cell number bias. This cell number bias increases as CAPE decreases, potentially because cells and their updrafts become narrower and more under‐resolved as CAPE decreases. The gross overproduction of precipitating shallow cells leads to overly efficient precipitation and inadequate detrainment of ice aloft, thereby diminishing the formation of robust stratiform rainfall regions. Decreasing model horizontal grid spacing from 3 to 1 or 0.333 km for low (<300 J kg −1 ) and high CAPE (>1,000 J kg −1 ) cases results in minimal change to cell number, depth, and convective‐to‐stratiform partitioning biases. This suggests that improving prediction of these convective properties depends on factors beyond solely increasing model resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multiscale Molecular Dynamics Simulations: Accelerating Conformational Sampling of Biomolecular Systems by Iterating All-Atom and Coarse-Grained Simulations

We developed the atomistic-coarse-grained multiscale MD simulation method in the OpenMM simulation package by iterating between the all-atom (AA) and coarse-grained (CG) MD simulations to enhance the sampling of biomolecular conformations. As the free energy surfaces are flattened during CG MD simulations, we can accelerate the transitions between different low-energy conformations. The AA-CG-AA cycles are repeated, facilitating the accelerated sampling of biomolecular conformations at a CG level, while the finer atomistic interactions are refined with AA simulators.

Do, Hung Nguyen↗
DOE OSTI2025

Stability Analysis of Power Systems with High Penetration of State-of-the-Art Inverter Technologies

With the increasing level of inverter-based resources (IBRs) in modern power systems, this paper presents a small-signal stability analysis for power systems comprising synchronous generators (SGs) and IBRs. Four types of inverter controls are considered: two grid-following (GFL) controls, with or without grid support functions; droop-based grid-forming (GFM) controls; and virtual oscillator control-based GFM. We also analyze the impact of STATCOM and synchronous condensers on system stability to assess their role in the energy mix transition. With the small-signal dynamic behavior of the major technologies modeled, this paper provides stringent stability assessments using the IEEE 39-bus benchmark system modified to simulate future power systems. The exhaustive test cases allow for (a) assessing the impacts of different types and controls of generation and supplementary grid assets, as well as system inertia and line impedance on grid stability, and (b) elucidating pathways for the stabilization of IBR-dominated power systems. The analysis also indicates that future power systems can be stabilized with only a fraction of the total generation as voltage sources without SGs or significant system inertia if they are well distributed. This study provides insights into future power system operations with a high level of IBRs that can also be used for planning and operation studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Impact of Hydrological Data on Power System Operational Studies: Preprint

Hydropower is expected to play an important role in maintaining grid reliability and flexibility as the share of of variable renewable energy increases. While the current hydropower operational models have been studied and used widely, they haven't been updated for decades to meet new performance standards. For example, current steady state and dynamic models often neglect hydrological conditions, which may lead to unrealistic expectations when relying on hydropower for energy and ancillary services. To study this impact, a multi-timescale hydrological model was created by leveraging the National Renewable Energy Laboratory-developed Multi-timescale Integrated Dynamics and Scheduling (MIDAS) tool. Using MIDAS, we compare the impact of considering hydrological conditions in a day-ahead unit commitment (DAUC) schedule on the reduced 240-bus Western Interconnect (WI) test system under winter and summer case studies. We show that neglecting current hydrological conditions of hydropower plants in power system models can lead to an overestimation of hydropower capabilities, which could lead to power balancing issues. For example, power system DAUC simulation results of our reduced test system show that in the case where hydrological conditions are not considered in the model, an approximate 31% overestimation of hydropower capabilities occurs in the summer case and approximately 60% occurs in the winter case compared to what is available. Additionally, results show an underestimation of WI day-ahead power system generation costs by approximately $54M - $80M in the weekly summer scenario and $116M - $126M in the weekly winter scenario. This analysis helps to underscore the importance of considering hydrological data in power system operational studies.

ENERGY PLANNING, POLICY, AND ECONOMY,HYDRO ENERGY↗
DOE OSTI2024

Impacts of Sea‐Level Rise on Coastal Groundwater Table Simulated by an Earth System Model With a Land‐Ocean Coupling Scheme

Abstract Sea‐level rise (SLR) poses a severe threat to the coastal environment through seawater intrusion into freshwater aquifers. The rising groundwater table also exacerbates the risk of pluvial, fluvial, and groundwater flooding in coastal regions. However, current Earth system models (ESMs) commonly ignore the exchanges of water at the land‐ocean interface. To address this gap, we developed a novel land‐ocean hydrologic coupling scheme in a state‐of‐the‐science ESM, the Energy Exascale Earth System Model version 2 (E3SMv2). The new scheme includes the lateral exchange between seawater and groundwater and the vertical infiltration of seawater driven by the SLR‐induced inundation. Simulations were performed with the updated E3SMv2 for the global land‐ocean interface to assess the impacts of SLR on coastal groundwater under a high CO 2 emission scenario. By the middle of this century, seawater infiltration on the inundated areas will be the dominant component in the land‐ocean coupling process, while the lateral subsurface flow exchange will be much smaller. The SLR‐induced seawater infiltration will raise the groundwater levels, enhance evapotranspiration, and increase runoff with distinct spatial patterns globally in the future. Although the coupling process is induced by SLR, we found topography and warming temperature have more control on the coupling impacts, probably due to the relatively modest magnitude of SLR during the selected future period. Overall, our study suggests significant groundwater and seawater exchange at the land‐ocean interface, which needs to be considered in ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

SCALE HTR-PROTEUS Benchmark Model

This dataset contains input and result files of computational simulations of HTR-PROTEUS benchmark with the latest version of SCALE code system. The simulations cover criticality control rod worth calculations as well as sensitivity analysis and uncertainty quantification. Users wanting to reproduce results from this dataset are required to obtain a license to the SCALE code system for which details on the distribution can be found here: https://www.ornl.gov/scale/releases

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Ocean Wave Energy Harvester with Oak Ridge Converter

Oceans can provide great potential for the American energy dominance. There is significant potential to utilize marine energy resources. In the United States, the total amount of marine energy available is equivalent to about 57% of the country's total power generation in 2019. Even if a fraction of this technical potential is harnessed, marine energy technologies could play a crucial role in fulfilling the nation's energy requirements. Marine energy resources are spread out geographically, and because more than 50% of the United States' population resides within 50 miles of the coastline, they are well-positioned to power local communities. These resources are also very dependable, making them a viable option for contributing to a consistent, trustworthy energy grid. Due to their predictable daily and seasonal patterns, marine energy resources can be integrated into our energy generation portfolio. On the other hand, ocean environment presents many challenges for cost-effective renewable energy conversion, including optimal control of ocean wave energy. This report presents a novel cost-effective energy conversion technique for ocean wave energy harvesters. The proposed system is simulated by using the Oak Ridge Converter to directly interface ocean wave energy source with the utility grid. The system description and simulation results are presented in detail. The results show that the proposed system is a cost effective and promising technology to reduce the infrastructure cost for ocean wave energy harvesters.

Sutton, Elizabeth [ORNL] (ORCID:0009000078885935)↗
DOE OSTI2026

Probabilistic Resource Adequacy Suite (PRAS) v0.8 Model Documentation

The Probabilistic Resource Adequacy Suite, or PRAS, is a software package for studying power system resource adequacy. It allows the user to simulate power system operations under a wide range of operating conditions, in order to study the system's risk of failing to meet demand due to a resource shortfall, and identify the time periods and regions in which that risk occurs. This reports documents version 0.8 of the tool.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Evaluation of Global Fire Simulations in CMIP6 Earth System Models

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of the El Niño–Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in CMIP5: (1) the simulated global burned area being less than half of that of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past 2 decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere midlatitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet–dry conditions, indicating the need to improve fuel wet-ness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest a postprocessing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products.

Wildfire, Earth system models↗
DOE OSTI2025

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Spatiotemporal quenches for efficient critical ground state preparation in the two-dimensional transverse field Ising model

Quantum simulators have the potential to shed light on the study of quantum many-body systems and materials, offering unique insights into various quantum phenomena. Although adiabatic evolution has been conventionally employed for state preparation, it faces challenges when the system evolves too quickly or the coherence time is limited. In such cases, shortcuts to adiabaticity, such as spatiotemporal quenches, provide a promising alternative. This paper numerically investigates the application of spatiotemporal quenches in the two-dimensional transverse field Ising model with ferromagnetic interactions, focusing on the emergence of the ground state and its correlation properties at criticality when the gap vanishes. We demonstrate the effectiveness of these quenches in rapidly preparing ground states in critical systems. Our simulations reveal the existence of an optimal quench front velocity at the emergent speed of light, leading to minimal excitation energy density and correlation lengths of the order of finite system sizes we can simulate. These findings emphasize the potential of spatiotemporal quenches for efficient ground state preparation in quantum systems, with implications for the exploration of strongly correlated phases and programmable quantum computing.

2-dimensional systems↗
DOE OSTI2025

Surrogate Constructed Scalable Circuits ADAPT-VQE in the Schwinger model

Inspired by recent advancements of simulating periodic systems on quantum computers, we develop a new approach, (SC)$^2$-ADAPT-VQE, to further advance the simulation of these systems. Our approach extends the scalable circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable ``Surrogate Constructed'' method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC)$^2$ methodology. Our approach allows for a wider set of classical computations, on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute properties of the Schwinger model - quantum electrodynamics for a single, massive fermion in $1+1$ dimensions - and show that our method can be used to accurately extrapolate to the continuum limit.

Gustafson, Erik [RIACS, Mtn. View] (ORCID:00000001↗
DOE OSTI2026

Simulations of Interacting Binary Systems—Pathways to Radio-bright Gamma-Ray Burst Progenitors

Although the association of gamma-ray bursts (GRBs) with massive stellar death is on firm footing, the nature of the progenitor system and the key ingredients required for a massive star to produce a GRB remain open questions. Here, we investigate the evolution of a 15–25M ⊙ massive star with a 10–15M ⊙ black hole (BH) using Modules for Experiments in Stellar Astrophysics. We quantify companion-influenced angular-momentum evolution over a stellar lifetime for orbital periods where tides are significant, varying stellar and BH masses, initial stellar spin, and accretion and dynamo prescriptions while tracking mass loss and angular momentum. The final spin is set by tidal torques versus stellar winds. For binaries that initially avoid Roche-lobe overflow, tides can spin up the star, but a late-stage expansion can drive tidal stripping; associated mass and angular-momentum loss can suppress spin-up. We find that massive-star BH binaries at comparable mass ratios may be potential GRB progenitors for short orbital periods (∼20–5 × 10 2 days) and long orbital periods (∼2 × 10 3 –4 × 10 3 days), where our suite of lifetime simulations reveals a favored parameter space with negligible mass loss and enough spin angular momentum to power a GRB jet. For initially nonrotating stars, this provides a lower limit on the final spin above a threshold estimate consistent with forming a postcollapse BH mass of 5–10M ⊙ with a spin parameter ≥ 0.5. For initially rapidly rotating stars, tidal interactions may sustain high spin when mass loss is negligible because the binary is not tidally synchronized.

Hernandez, Angel [Univ. of Colorado, Boulder, CO (↗