Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD 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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Horne et al. (2026) supporting files - WRF-LES model outputs for a summer heatwave event on June 2025 in Baltimore, MD

Brief Description Shown is the supporting information for Horne et al. (2026). These files include all outputs from the WRF model simulations and the observational datasets used for comparison in the study. Scripts are provided so users can recreate the manuscript's figures using the provided observational and modeling data. For more information regarding the study, please contact the primary author of the associated manuscript, Jason Horne. Horne, J. P., Pan, Y., Davis, K. J., Waugh, D., Ahlswede, B. J., Prince, N. E. (2026). Simulating near-surface environments in urban neighborhoods using WRF-LES: A case study of classic atmospheric boundary layer (ABL) during a heatwave event JAMES. (to be submitted)

atmosphere↗
DOE OSTI2026

Horne et al. (2026) supporting files - WRF-LES model outputs for a summer heatwave event on June 2025 in Baltimore, MD

Brief Description Shown is the supporting information for Horne et al. (2026). These files include all outputs from the WRF model simulations and the observational datasets used for comparison in the study. Scripts are provided so users can recreate the manuscript's figures using the provided observational and modeling data. For more information regarding the study, please contact the primary author of the associated manuscript, Jason Horne. Horne, J. P., Pan, Y., Davis, K. J., Waugh, D., Ahlswede, B. J., Prince, N. E. (2026). Simulating near-surface environments in urban neighborhoods using WRF-LES: A case study of classic atmospheric boundary layer (ABL) during a heatwave event JAMES.

atmosphere↗
DOE OSTI2025

Radiative processes on a quantum computer

Radiative processes, where a photon/neutrino is emitted as a result of a collision or decay of a particle, play a central role in atomic, nuclear, and particle physics. Their rate is determined by certain off-diagonal matrix elements between different initial and final states. We propose a method to compute them using quantum computers. It relies on a single extra qubit that, in a certain sense, represents the photon/neutrino. The generic formula relating this matrix element to the amplitude and frequency of oscillations of the extra qubit is derived for the near-resonance case. Furthermore, we demonstrate the feasibility of the method by using it in actual quantum computations and simulations of simple systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Practical low-temperature gasoline combustion for very high efficiency off-road, medium- and heavy-duty engines

Low-temperature gasoline combustion (LTGC) with additive-mixing fuel injection (AMFI) is a new combustion strategy that has been demonstrated to deliver 9%–25% better brake thermal efficiency than similar-sized market-leading diesel engines over the operating map. Moreover, the LTGC-AMFI engine shows near-zero smoke, and NOx emissions are 4–100 times lower than those of a diesel, sufficiently low that no aftertreatment, or only passive NOx aftertreatment, would be sufficient (diesel exhaust fluid is not required). LTGC-AMFI combustion is based on kinetically controlled compression ignition of a dilute charge with a variable amount of low-to-moderate fuel stratification. Fast combustion control is provided by adding minute amounts of an ignition-enhancing additive into the fuel each engine cycle to control its reactivity. This strategy was used to operate a medium-duty (MD) LTGC-AMFI engine at loads from idle to 16.3 bar BMEP and speeds from 600 to 2400 rpm with regular E10 gasoline, which covers nearly the entire operating map of a typical MD engine. Turbine-out temperatures were sufficient for an oxidation catalyst to control hydrocarbon and CO emissions. Autonomie simulations over the GEM ARB Transient and the GEM 55 mph Cruise driving cycles for class-6 trucks using this technology showed fuel economies of 8.1 and 11.4 mpg-gasoline-equivalent, respectively, corresponding to 18.6% and 13.4% improvements over a similar-size diesel engine. Engine-out NOx emissions were 0.024 and 0.01 g/bhp-h, respectively, well below current U.S. emission standards. These results show that switching from diesel to LTGC-AMFI engines would greatly reduce greenhouse gas (GHG) emissions for off-road, MD and HD applications, which will continue to rely on combustion engines because electrification is not practical in the foreseeable future. Finally, with their reduced fuel consumption, the lower cost of gasoline compared to diesel fuel, and much lower aftertreatment costs, LTGC-AMFI engines also offer a significantly lower total cost of ownership.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Fundamental insights on ignition and combustion of natural gas in an active fueled pre-chamber spark-ignition system

Pre-chamber spark-ignition (PCSI), either fueled or non-fueled, is a leading concept with the potential to enable diesel-like efficiency in medium-duty (MD) and heavy-duty (HD) natural gas (NG) engines. However, the inadequate scientific base and simulation tools to describe/predict the underlying processes governing PCSI systems is one of the key barriers to market penetration of PCSI for MD/HD NG engines. To this end, experiments were performed in a heavy-duty, optical, single-cylinder engine fitted with an active fueled PCSI module. Furthermore, the spatial and temporal progress of ignition and subsequent combustion of lean-burn natural gas using PCSI system were studied using optical diagnostic imaging and heat release analysis based on main-chamber and pre-chamber pressure measurements.

42 ENGINEERING↗
DOE OSTI2026

A new self-adaptive reconstruction method to identify defects through Wigner–Seitz approach

A new self-adaptive reconstruction method based on local atomic structure at any given molecular dynamics (MD) step has been developed in this article. The method can be used in Wigner–Seitz defect analysis approach to correctly and efficiently explore the information of both point defects and complex defect clusters (e.g. dislocation loops and voids) formed after a displacement cascade where the cascade interacts with grain boundaries and/or dislocations. The algorithm and validation are provided in detail. Results for identification of radiation defects during and after cascades interacting with a dislocation network show that the new method can well recognize all simple and complex defects and defect clusters. Thus, this new method provides a totally new way to explore the density and size of radiation defects at atomic scale after complex MD evolution processes, providing correct information to understand and predict radiation damage in materials through atomic simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Hydrogen Bonds and H 3 O + Formation at the Water Interface with Formic Acid Covered Anatase TiO 2

Carboxylic acid-modified TiO 2 surfaces in aqueous environment are of widespread interest, yet atomic-scale understanding of their structure is limited. In this work, we investigate formic acid (FA) on anatase TiO 2 (101) (A-101) in contact with water using density functional theory (DFT) and ab-initio molecular dynamics (AIMD). Isolated FA molecules adsorbed in a deprotonated bridging bidentate (BD) form on A-101 are found to remain stable at the interface with water, with the acid proton transferred to a surface oxygen to form a surface bridging hydroxyl (O br H). With increasing FA coverage, adsorbed monolayers of only BD and successively of alternating monodentate (MD) and BD species give rise to a higher concentration of surface O br H’s. Simulations of these adsorbed monolayers in water environment show that some protons are released from the surface O br H’s to water resulting in a negatively charged surface with nearby solvated H 3 O + ions. These results provide insight into the complex acid-base equilibrium between an oxide surface, adsorbates and water and can also help obtain a better understanding of the wetting properties of chemically modified TiO 2 surfaces.

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

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗
DOE OSTI2025

A comparison of probabilistic generative frameworks for molecular simulations

Generative artificial intelligence is now a widely used tool in molecular science. Despite the popularity of probabilistic generative models, numerical experiments benchmarking their performance on molecular data are lacking. Here, in this work, we introduce and explain several classes of generative models, broadly sorted into two categories: flow-based models and diffusion models. We select three representative models: neural spline flows, conditional flow matching, and denoising diffusion probabilistic models, and examine their accuracy, computational cost, and generation speed across datasets with tunable dimensionality, complexity, and modal asymmetry. Our findings are varied, with no one framework being the best for all purposes. In a nutshell, (i) neural spline flows do best at capturing mode asymmetry present in low-dimensional data, (ii) conditional flow matching outperforms other models for high-dimensional data with low complexity, and (iii) denoising diffusion probabilistic models appear the best for low-dimensional data with high complexity. Our datasets include a Gaussian mixture model and the dihedral torsion angle distribution of the Aib9 peptide, generated via a molecular dynamics simulation. We hope our taxonomy of probabilistic generative frameworks and numerical results may guide model selection for a wide range of molecular tasks.

Artificial intelligence↗
DOE OSTI2025

Operation and Control of Electric Vehicle Charger with Enhanced Dynamic Performance Under Non-Ideal Grid Voltage Condition

This paper presents a three-phase electric vehicle charger connected to the grid, featuring multiple boost converters on the DC side, specifically designed to ensure smooth, oscillation-free power transfer during unsymmetrical voltage sags. Precise control mechanisms are implemented on the boost converter side to regulate both voltage and current on the electric vehicle side, thereby maintaining optimal charging conditions. The control architecture for both the grid-connected and boost converter components is based on the Lyapunov energy function, which is employed to achieve superior dynamic performance and stability. The system's robustness and reliability are demonstrated through its ability to maintain stable operation and efficient power transfer despite fluctuations in grid conditions. Furthermore, the implementation of Lyapunov-based control ensures rapid response and minimal energy loss, enhancing the overall efficiency of the system. To validate the effectiveness of this approach, a comprehensive model of the system was developed and tested using MATLAB/Simulink, with detailed computer simulations conducted across various significant case studies.

DC-DC boost converter↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

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

Biogeochemistry simulations for the Salt Marsh Accretion Response to Temperature eXperiment (SMARTX)

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.Included in this data package are PFLOTRAN input (PFLOTRAN input files and chemical database) files for simulating single column biogeochemistry, root, and tide interactions at the Global Change Research Wetland (Kirkpatrick Marsh; Edgewater, MD). The biogeochemical network includes soil organic matter decomposition, nitrogen, iron, and sulfur cycling, and methanogenesis. Reduced species can be oxidized and plant processes include oxygen and nutrient priming, methane release, and nutrient uptake.Inputs:TAI_database.dat - geochemical database for reactions, more information on database structure and variables can be found here https://www.pflotran.org/documentation/user_guide/cards/pages/geochemical_database.htmlswamp.in - input file for biogeochemical network in PFLOTRANswamp_eCO2.in - input file for biogeochemical network in PFLOTRAN with input gas partial pressures/concentrations adjusted for elevated CO2 treatmentsOutputs:swamp_obs_0.tec - hourly porewater concentrations from from multiple depths in the soil columnswamp_eCO2_obs_0.tec - hourly porewater concentrations from multiple depths in the soil column for elevated CO2 treatmentsPFLOTRAN code access: https://github.com/fmyuan/pflotran-elm-interface.git

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of U10Mo Fuel Plate Performance Modeling Over Hot Isostatic Press and Hydraulic Bending for MURR DDE Plates

The United States High Performance Research Reactor Program’s objective is to reduce the amount of highly enriched uranium currently implemented in research reactors. The conversion of these research reactors requires designing a monolithic U10Mo plate fuel, with the fuel plate geometry being dependent on each research reactor. The process of forming the plates includes a hot isostatic pressing (HIP) to manufacture a prototypic plate. In the case of the Missouri University Research Reactor (MURR) design demonstration element (DDE) plate manufacture, plates that have been through HIP are then curved using dies and a hydraulic press to impart the desired curvature. Both fabrication processes impart residual stresses into each fuel plate region, with the curvature of the plates taking some regions of the fuel plate up to their material yield stresses, accompanied by plastic strain. The amount of plastic strain and stress imparted onto each MURR DDE plate is determined by the radius of curvature, thickness of each region, and overall width of the fuel plates. Furthermore, this work aims to predict the yield stresses and strain using ABAQUS to simulate the proposed fabrication process of the MURR DDE plates, accompanied by discussion over the stresses and strains as to their relation to nuclear fuel performance and the impact they will have during early irradiation.

ABAQUS↗
DOE OSTI2025

Sensitivity of the Shallow‐To‐Deep Convective Transition to Moisture and Wind Shear in the Amazon

Deep convection is the primary influence on weather and climate in tropical regions. However, understanding and simulating the shallow‐to‐deep (STD) convective transition has long been challenging. Here, we conduct high‐resolution numerical simulations to assess the environmental controls on the evolution of isolated convection in the Amazon during the wet season. The large‐scale forcing derived through a constrained variational analysis approach for the GoAmazon2014/5 Experiment is used in the simulations. Through sensitivity experiments, we examine the relative importance of moisture and wind shear in controlling the shallow‐to‐deep convective transition for isolated convective events. Convection exhibits the greatest sensitivity to humidity within the lowest 1.5 km, where a 4 mm reduction in column water vapor nearly suppresses ice water formation on deep convective days. In contrast, a reduction in column water vapor in the free troposphere by a factor of two or more is necessary to produce a comparable impact on convection. Increasing low‐level wind speed from 6 to 9 m s -1 enhances afternoon deep convection, raising the cloud ice mixing ratio by approximately 25%. Conversely, upper‐level wind shear reveals the weakest correlation with daytime convection in our simulations. Our results help characterize the role of moisture and wind shear on the STD transition and our understanding of the underlying mechanisms.

environmental controls↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2022

Transient nature of fast relaxation in metallic glass

Metallic glasses exhibit fast mechanical relaxations at temperatures well below the glass transition, one of which shows little variation with temperature known as nearly constant loss (NCL). Despite the important implications of this phenomenon to deformation, the origin of the relaxation is unclear. Through molecular dynamics simulations of a model metallic glass, Cu 64.5 Zr 35.5 , we implement molecular dynamics dynamical mechanical spectroscopy (MD-DMS) with system stress decomposed into atomic-level stresses to identify the group of atoms responsible for NCL. This work demonstrates that NCL relaxation is due to transient groups of atoms that revert to the typical atomic-level viscoelastic behavior over picosecond timescales. They are homogenously distributed throughout the glass and have no outstanding features, rather than having defect-like local structure as previously reported.

36 MATERIALS SCIENCE↗