Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulator”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Coarse-grained discrete particle simulations of particle segregation in rotating fluidized beds in vortex chambers [Discrete particle simulations of particle segregation in rotating fluidized beds in vortex chambers]

Vortex chambers allow the generation of rotating fluidized beds, offering high-G intensified gas-solid contact, gas-solids separation and solids-solids segregation. Focusing on binary particle mixtures and fixing the density and diameter of the heavy/large particles, transient batch CFD-coarse-grained DPM simulations were carried out with varying densities or sizes of the light/small particles to evaluate to what extent combining these three functionalities is possible within a vortex chamber of given design. Both the rate and quality of segregation were analyzed. Within a relatively wide density and size range, fast and efficient segregation takes place, with an inner and slower rotating bed of the lighter/small particles forming within the outer and faster rotating bed of the heavier/large particles. Simulations show that the contamination of the outer bed with lighter particles occurs more easily than contamination of the inner bed with heavier particles and increases with decreasing difference in size or density of the particles. Bubbling in the inner bed is observed with an inner bed of very low density or small particles. Porosity plots show that vortex chambers with a sufficient number of gas inlet slots have to be used to guarantee a uniform gas distribution and particle bed. Lastly, the flexibility of particle segregation in vortex chambers with respect to the gas flow rate is demonstrated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Equilibria and Rate Phenomena from Atomistic to Mesoscale: Simulation Studies of Magnetite

Batteries are complex devices composed of multiple components that operate far from equilibrium and may operate under extreme stress and varying loads. Research of isolated battery components are valuable to the fundamental understanding of the physical processes occurring within each constituent element. When the components are integrated into a full device and operated under realistic conditions, it can be difficult to decouple the physical processes that occur across multiple interfaces and multiple length scales. Thus, the physical processes studied in isolated components may change in a full battery setup or may be irrelevant to performance. Simulation studies on many length scales play a key role in the analysis of experiments and in the elucidation of the relevant physical processes impacting performance. In this Account, we aim to highlight the use of modeling on multiple length scales to identify rate limiting phenomena in lithium-ion batteries. To illustrate the utility of modeling, we examine lithium-ion batteries with nanostructured magnetite, Fe 3 O 4 , as the positive electrode active material against a solid Li 0 negative electrode. Due to continuous operation away from equilibrium, batteries exhibit highly nonideal behavior, and a model that aims to reproduce behavior under realistic operating conditions must be able to capture the physics occurring on the length scales relevant to the performance of the system. It is our experience that limiting behavior in lithium-ion batteries can be observed on the atomic scale and up through the electrode scale and thus, predictive models must be capable of integrating and communicating physics across multiple length scales. Magnetite is studied as an electrode material for lithium-ion batteries, but it is found to suffer from slow solid-state transport of lithium, slow reaction kinetics, and poor cycling. Magnetite (Fe 3 O 4 ) is a material capable of undergoing multiple electron transfers (MET), and can accept up to eight lithium per formula unit (Li 8 Fe 3 O 4 ). Magnetite, (Fe 8a 3+ )[Fe 3+ Fe 2+ ] 16d O 4,32e 2– , has a close-packed inverse spinel structure and undergoes intercalation, structural rearrangement, and conversion reactions upon full lithiation. To overcome solid-state transport resistances, magnetite can be nanostructured to decrease Li + diffusion lengths, and this has been shown to increase capacity. Additionally, unique architectures incorporating both carbon and Fe 3 O 4 have shown to alleviate transport and cycling issues in the material. Here, we solely address traditional composite electrodes, in which Fe 3 O 4 is synthesized as nanoparticles and combined with additives to fabricate the electrode. In the case of nanoparticulate magnetite, it has been found that the electrode fabrication process results in the formation of micrometer-sized agglomerates of the Fe 3 O 4 nanoparticles, introducing a secondary structural motif. The agglomerates may form in one or more fabrication processes, and their elimination may not be straightforward or warranted. Here, we highlight the impact of these secondary formations on the performance of the Fe 3 O 4 lithium-ion battery. We illustrate how simulations can be used to design experiments, prioritize research efforts, and predict performance.

25 ENERGY STORAGE↗
DOE OSTI2023

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

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

Combined Force-Frequency Sampling for Simulation of Systems Having Rugged Free Energy Landscapes.

An adaptive, machine learning-based sampling method is presented for simulation of systems having rugged, multidimensional free energy landscapes. The method's main strength resides in its ability to learn both from the frequency of visits to distinct states and the generalized force estimates that arise in a system as it evolves in phase space. This is accomplished by introducing a self-integrating artificial neural network, which generates an estimate of the free energy directly from its derivatives. The usefulness of the proposed combined approach is examined in the context of two concrete examples, namely, an alanine dipeptide molecule in water and a polymer diffusing through a narrow pore. This new method is found to be robust, faster, and more accurate than approaches that rely only on frequency-based or generalized force-based estimations. After combining the proposed approach with overfill protection and support for sparse data storage and training, the method is shown to be more effective than comparable, previously available techniques and capable of scaling efficiently to larger numbers of collective variables.

Sevgen, Emre↗
DOE OSTI2025

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

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

Crystal Growth Simulations To Establish Physically Relevant Kinetic Parameters from the Empirical Kolmogorov-Johnson-Mehl-Avrami Model

A series of simulations was performed to enable interpretation of the material and physical significance of the parameters defined in the Kolmogorov, Johnson and Mehl, and Avrami (KJMA) rate expression commonly used to describe phase boundary controlled reactions of condensed matter. The parameters k, n, and t 0 are shown to be highly correlated, which if unaccounted for seriously challenge mechanistic interpretation. It is demonstrated that rate measurements exhibit an intrinsic uncertainty without precise knowledge of the location and orientation of nucleation with respect to the free volume into which it grows. More significantly, it is demonstrated that the KJMA rate constant k is highly dependent on sample size. However, under the simulated conditions of slow nucleation relative to crystal growth, sample volume and sample anisotropy correction affords a means to eliminate the experimental condition dependence of the KJMA rate constant, k, producing the material-specific parameter, the velocity of the phase boundary, v pb .

Dill, Eric D.↗
DOE OSTI2012

First Principles Simulations of the Electrochemical Lithiation and Delithiation of Faceted Crystalline Silicon

Silicon is of significant interest as a next-generation anode material for lithium-ion batteries due to its extremely high capacity. The reaction of lithium with crystalline silicon is known to present a rich range of phenomena, including electrochemical solid state amorphization, crystallization at full lithiation of a Li15Si4 phase, hysteresis in the first lithiation–delithiation cycle, and highly anisotropic lithiation in crystalline samples. Very little is known about these processes at an atomistic level, however. To provide fundamental insights into these issues, we develop and apply a first principles, history-dependent, lithium insertion and removal algorithm to model the process of lithiation and subsequent delithiation of crystalline Si. The simulations give a realistic atomistic picture of lithiation demonstrating, for the first time, the amorphization process and hinting at the formation of the Li15Si4 phase. Voltages obtained from the simulations show that lithiation of the (110) surface is thermodynamically more favorable than lithiation of the (100) or (111) surfaces, providing an explanation for the drastic lithiation anisotropy seen in experiments on Si micro- and nanostructures. Analysis of the delithiation and relithiation processes also provides insights into the underlying physics of the lithiation–delithiation hysteresis, thus providing firm conceptual foundations for future design of improved Si-based anodes for Li ion battery applications.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2024

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. Here, in this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

Hamiltonians↗
DOE OSTI2012

An investigation of the effect of pore scale flow on average geochemical reaction rates using direct numerical simulation

The scale-dependence of geochemical reaction rates hinders their use in continuum scale models intended for the interpretation and prediction of chemical fate and transport in subsurface environments such as those considered for geologic sequestration of CO 2 . Processes that take place at the pore scale, especially those involving mass transport limitations to reactive surfaces, may contribute to the discrepancy commonly observed between laboratory-determined and continuum-scale or field rates. In this study we investigate the dependence of mineral dissolution rates on the pore structure of the porous media by means of pore scale modeling of flow and multicomponent reactive transport. The pore scale model is composed of high-performance simulation tools and algorithms for incompressible flow and conservative transport combined with a general-purpose multicomponent geochemical reaction code. The model performs direct numerical simulation of reactive transport based on an operator-splitting approach to coupling transport and reactions. The approach is validated with a Poiseuille flow single-pore experiment and verified with an equivalent 1-D continuum-scale model of a capillary tube packed with calcite spheres. Using the case of calcite dissolution as an example, the high-resolution model is used to demonstrate that nonuniformity in the flow field at the pore scale has the effect of decreasing the overall reactivity of the system, even when systems with identical reactive surface area are considered. In conclusion, the effect becomes more pronounced as the heterogeneity of the reactive grain packing increases, particularly where the flow slows sufficiently such that the solution approaches equilibrium locally and the average rate becomes transport-limited.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Interactions Between the Amazonian Rainforest and Cumuli Clouds: A Large‐Eddy Simulation, High‐Resolution ECMWF, and Observational Intercomparison Study

Abstract The explicit coupling at meter and second scales of vegetation's responses to the atmospheric‐boundary layer dynamics drives a dynamic heterogeneity that influences canopy‐top fluxes and cloud formation. Focusing on a representative day during the Amazonian dry season, we investigate the diurnal cycle of energy, moisture and carbon dioxide at the canopy top, and the transition from clear to cloudy conditions. To this end, we compare results from a large‐eddy simulation technique, a high‐resolution global weather model, and a complete observational data set collected during the GoAmazon14/15 campaign. The overall model‐observation comparisons of radiation and canopy‐top fluxes, turbulence, and cloud dynamics are very satisfactory, with all the modeled variables lying within the standard deviation of the monthly aggregated observations. Our analysis indicates that the timing of the change in the daylight carbon exchange, from a sink to a source, remains uncertain and is probably related to the stomata closure caused by the increase in vapor pressure deficit during the afternoon. We demonstrate quantitatively that heat and moisture transport from the subcloud layer into the cloud layer are misrepresented by the global model, yielding low values of specific humidity and thermal instability above the cloud base. Finally, the numerical simulations and observational data are adequate settings for benchmarking more comprehensive studies of plant responses, microphysics, and radiation.

Vilà‐Guerau de Arellano, J.↗
DOE OSTI2020

The Role of the Mean State on MJO Simulation in CESM2 Ensemble Simulation

Here we examine the role of the mean state in the propagation of the Madden-Julian oscillation (MJO) over the Maritime Continent (MC). We use an ensemble of simulations made with a single model—the Community Earth System Model version 2—to assess the effect of the mean state that is unaffected by that of model components such as parameterization schemes. Results show that the background meridional moisture gradient is much steeper over the MC region in the periods with a stronger MJO propagation. Column water vapor budget of the MJO strongly suggests that the simulated mean state affects MJO via its impacts on moisture dynamics—a greater advection of mean moisture by MJO wind in the MC region is responsible for the anomalous MJO activity.

58 GEOSCIENCES↗
DOE OSTI2021

The Influence of Ocean Coupling on Simulated and Projected Tropical Cyclone Precipitation in the HighResMIP–PRIMAVERA Simulations

This study aims to quantify the impacts of atmosphere–ocean coupling on simulated and projected tropical cyclone (TC) precipitation globally. We used global climate model (GCM) simulations over 1950–2050 from the High Resolution Model Intercomparison Project (HighResMIP) and compared its fully coupled atmosphere–ocean GCMs (AOGCMs) with atmosphere-only GCMs (AGCMs). We find that ocean coupling generally leads to decreased TC precipitation over ocean and land. Large-scale sea surface temperature (SST) biases are critical drivers of the precipitation difference, with secondary contributions from local TC–ocean feedbacks via SST cold wakes. The two driving factors, attributed to ocean coupling in the AOGCMs, influence TC precipitation in association with decreased TC intensity and specific humidity. The AOGCMs and AGCMs consistently project TC precipitation increases in 2015–2050 relative to 1950–2014 over ocean for all basins, and for landfalling TCs in the North Atlantic and western North Pacific.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of E3SM Simulated Aerosols and Aerosol‐Cloud Interactions Across GCM and Convection‐Permitting Scales

This paper introduces an Earth system modeling testbed for predicting aerosols and aerosol‐cloud interactions (ACIs) at convection‐permitting scales. Using the Energy Exascale Earth System Model (E3SM) version 2 with a four‐mode Modal Aerosol Module, we conduct simulations at 3.25 km resolution on a regionally refined mesh (RRM) across four regions with distinct aerosol and cloud regimes. Results are compared with the standard 100 km E3SM configuration and evaluated against satellite, aircraft, and ground‐based observations. We find that increasing model resolution improves heavy precipitation simulation but amplifies positive bias in light drizzle at coarse resolution. These resolution‐induced changes affect cloud and aerosol properties to varying degrees across regions. Generally, cloud cover and liquid water path (LWP) show better agreement with satellite retrievals at 3.25 km, though surface‐based comparisons suggest otherwise. Aerosol composition remains poorly represented at both resolutions. The RRM increases Aitken mode aerosol number concentrations via enhanced new particle formation. However, accumulation mode aerosols are decreased at higher resolution as aerosol removals become more efficient. This partially contributes to fewer cloud condensation nuclei (CCN) and lower cloud droplet number concentrations (N d ), which produces larger model biases in some scenarios. These findings suggest that solely increasing horizontal resolution to kilometer scales is insufficient to broadly improve aerosol and cloud predictions without concurrent advancements in physical and chemical process representations. Nonetheless, the RRM moderately improves key ACI relationships such as CCN‐N d correlation, reflecting enhanced aerosol activation representation. The LWP‐N d relationship is also better captured by RRM, suggesting a better characterization of LWP adjustment.

Huang, Meng [Pacific Northwest National Laboratory↗
DOE OSTI2012

ELMy H-mode linear simulation with 3-field model on experimental advanced superconducting tokamak using BOUT++

H-mode plasmas with ELM (edge localized mode) have been realized on experimental advanced superconducting tokamak (EAST) with 2.45 GHz low hybrid wave at P{sub LHW}{approx}1 MW in 2010. Data from EAST experiments including magnetic geometry, measured pressure profiles, and calculated current profiles are used to investigate the physics of ELM utilizing the BOUT++ code. Results from linear simulations show that the ELMs in EAST are dominated by resistive ballooning modes. When the Lundquist number (dimensionless ratio of the resistive diffusion time to the Alfven time) is equal to or less than 10{sup 7}, the resistive ballooning modes are found to become unstable in the ELMy H-mode plasma. For a fixed pedestal pressure profile, increasing plasma current generates more activities of low-n ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗