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 811 records · Page 45

Non-equilibrium critical scaling and universality in a quantum simulator

Universality and scaling laws are hallmarks of equilibrium phase transitions and critical phenomena. However, extending these concepts to non-equilibrium systems is an outstanding challenge. Despite recent progress in the study of dynamical phases, the universality classes and scaling laws for non-equilibrium phenomena are far less understood than those in equilibrium. In this work, using a trapped-ion quantum simulator with single-spin resolution, we investigate the non-equilibrium nature of critical fluctuations following a quantum quench to the critical point. We probe the scaling of spin fluctuations after a series of quenches to the critical Hamiltonian of a long-range Ising model. With systems of up to 50 spins, we show that the amplitude and timescale of the post-quench fluctuations scale with system size with distinct universal critical exponents, depending on the quench protocol. While a generic quench can lead to thermal critical behavior, we find that a second quench from one critical state to another (i.e. a double quench) results in a new universal non-equilibrium behavior, identified by a set of critical exponents distinct from their equilibrium counterparts. Our results demonstrate the ability of quantum simulators to explore universal scaling beyond equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Geometric Delocalization in Two Dimensions

We demonstrate the existence of transient two-dimensional surfaces where a random-walking particle escapes to infinity in contrast to localization in standard flat two-dimensional space. We first prove that any rotationally symmetric two-dimensional membrane embedded in flat three-dimensional space cannot be transient. Then we formulate a criterion for the transience of a general asymmetric two-dimensional membrane. We use it to explicitly construct a class of transient two-dimensional manifolds with a nontrivial metric and height function but “zero average curvature,” which we dub “tablecloth manifolds.” The absence of the logarithmic infrared divergence of the Laplace-Beltrami operator in turn implies the absence of weak localization, nonexistence of bound states in shallow potentials, and breakdown of the Mermin-Wagner theorem and Kosterlitz-Thouless transition on the tablecloth manifolds, which may be realizable in both quantum simulators and corrugated two-dimensional materials.

Anderson localization↗

Open Circuit Switch Fault Management Method of a Multi-Phase Synchronous Buck Converter for EV Charging Application

Multi-phase synchronous buck converter (MSBC) is becoming popular for the electric vehicle (EV) charging application because of its higher efficiency, smaller passive component sizes, bidirectional power flow, and output current ripple reduction through interleaving. This converter also provides higher reliability than the traditional single phase converter since its modular structure enables fault tolerant operation. Fault management in the converter requires the detection of switch fault, shut down of the faulty phase, and reconfiguration of the PWM signals for the healthy phases to avoid circulating current and high output current ripple. In addition, it is desirable that the fault management method can be implemented using the micro-controller unit (MCU) or digital signal processor (DSP), which are widely used to control power converters. This paper presents a fault management method for detecting switch faults and reconfiguring healthy phases using a DSP. The method is validated through simulations and experiments using a commercial off-the-shelf power stack.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Spatio-Temporal Clustering of Seismicity Enabled by Off-Fault Plasticity

While significant progress has been made in understanding earthquake source processes in linear elastic domains, the effect of more realistic rheologies including plasticity is poorly understood. Here, in this work, we simulate the sequence of earthquake and aseismic slip of a 2D antiplane rate-and-state fault embedded in a full-space elastic-plastic bulk. We show that off-fault plasticity may lead to partial ruptures as well as temporal clustering of seismic events. Furthermore, the interaction of fault slip and off-fault plasticity results in pockets of slip deficit. While the energy dissipated through plastic deformation remains a small fraction of the total energy budget, its impact on the source characteristics is disproportionally large through the redistribution of stresses and viscous relaxation. Our results suggest a new mechanism of dynamic heterogeneity in earthquake physics that may have important implications on earthquake size distribution and energy budget.

58 GEOSCIENCES↗

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗

Gigaton commercial-scale carbon storage and mineralization potential in stacked Columbia River basalt reservoirs

This work presents a detailed supercritical CO 2 storage resource estimation for the stacked basalt reservoirs in the Grande Ronde Basalt of the Columbia River Basalt Group in eastern Washington and Oregon. The assessment aims to derisk the commercialization potential of geologic carbon storage in basalt by leveraging both structural and mineralization trapping of CO 2 in basalt. The structural closures formed by anticlinal ridges and synclinal valleys in Yakima Fold Belt are excellent physical traps to accommodate injected supercritical CO 2 . Rigorous hydraulic testing, well logs and simulation results from the Wallula Basalt Pilot #1 well showed the occurrence of 17 suitable permeable injection zones (up to 2,496 mD) intercalated with dense seals (~2.6E-10 mD) in the Grand Ronde Basalt. In addition, geochemical studies showed fast reactions between supercritical CO 2 and dissolved basalt minerals to form stable carbonates. In conclusion, our calculation indicates up to 40 gigatons (P90) of mineralization storage resources exist in the Grande Ronde Basalt reservoirs.

58 GEOSCIENCES↗

First principles molecular dynamics maps out complete mineral surface acidity landscape

In “Interfacial structures and acidity constants (pKa) of goethite from first principles molecular dynamics simulations,” authors Y. Zhang, X. Lui, J. Cheng, and X. Lu apply First Principles molecular dynamics (FPMD, also called Density Functional Theory MD, DFT/MD, or ab initio MD, AIMD), to evaluable the complete set of acidity constants (pK a ) of the hydroxyl groups on the most prominent goethite crystal facets. The pK a of these OH and OH$^+_2$ groups are compared with available data from the multisite complexation (MUSIC) model traditionally used to estimate pK a on mineral surfaces. The authors have presented eloquent rationale for the importance and implications of understanding goethite acidity constants in room temperature geochemistry settings. In this paper, I focus on the computational aspects, the strengths of FPMD, and its possibilities.

58 GEOSCIENCES↗

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗

Comprehensive parametric study of CO 2 sequestration in deep saline aquifers

Carbon dioxide injection in deep saline aquifers is a key method for permanently sequestering anthropogenic CO 2 . Here this study employs a reactive transport model to explore mineral precipitation/dissolution and its impact on reservoir properties in deep saline aquifers. We also assess capillary pressure and relative permeability hysteresis on various CO 2 trapping mechanisms. Results from this study reveal the significant influence of initial brine composition on mineral precipitation/dissolution. The dissolution and precipitation of minerals have different effects around the wellbores compared to the overall reservoir. Additionally, salt concentration (Ca ++ and Mg ++ ) and quartz surface area affect CO 2 mineralization, while Na + impacts halite precipitation, altering flow properties. The effect of capillary pressure is significant, as including the capillary pressure in the simulation case resulted in significantly improved CO 2 trapping, achieving almost total dissolution of the injected CO 2 in around 300 years. This study offers novel insights into the interactions of reservoir minerals, brine properties, and the injected CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Medium Voltage Solid State Transformer for Extreme Fast Charging Applications

A modular and scalable solid state transformer (SST) with direct medium voltage (MV) AC connectivity is proposed to enable DC extreme fast charging (XFC) of electric vehicles. Single-phase-modules (SPMs), each consisting of an active-front-end (AFE) stage and an isolated DC-DC stage, are connected in input-series-output-parallel (ISOP) configuration. The modular hardware is co-designed with decentralized control of the DC-DC stages where voltage and power balancing are achieved by each SPM using only its local sensor feedback; a centralized controller (CC) regulates the low voltage (LV) DC bus through the AFE stages without any sensor feedback form the SPMs. The controller architecture contrasts sharply with the prior art for MV AC to LV DC SSTs where high-speed bidirectional communication among SPMs and a CC are required for module-level voltage and power balancing, which severely limits the scalability and practical realization of higher voltage and higher power units. Detailed small-signal analysis and controller design guidelines are developed. Furthermore, a soft start-up strategy is presented. The proposed converter and control structure are validated through simulation and experimental results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Hardware-based Advanced Electromagnetic Transient Simulation for A Large-Scale PV Plant in Real Time Digital Simulator

Power electronics-based resources, such as high-voltage direct current (HVdc) substations, photovoltaic (PV) plants, wind plants, electric vehicle charging stations, and energy storage systems, are increasingly being integrated within the power grid. Recently, multiple reports have emphasized the necessity for high-fidelity electromagnetic transient (EMT) simulations of these large-scale power electronics-based resources to accurately understand their behavior in power grids. However, performing hardware-based EMT simulations with high-fidelity models for such large power electronics systems is challenging due to the small time-step requirements and the involvement of a large number of states. This paper presents the implementation of hardware-based high-fidelity EMT dynamic model of a large-scale PV plant, accomplished through custom model development using the specific-C language in real-time digital simulator hardware (RTDS) and software (RSCAD).

Choi, Jongchan↗

A Fast Framework for Generating Radioactive Mixture Spectra and Its Application to Remote High-Performance Mixture Identification

Remote detection of radioactive materials in mixtures using handheld or portal detectors remains a challenge because of factors such as low concentration, environmental interference, sensor noise, and other complications. This work introduces a fast framework for generating realistic mixture spectra. Moreover, we present mixture isotope identification using data generated by the fast framework. Researchers have examined a range of conventional and recent algorithms within the fields of machine learning and deep learning. An application to uranium enrichment-level prediction has been included. Extensive simulation experiments validated the efficacy of the proposed framework.

GADRAS↗

Quantum tensor network algorithms for evaluation of spectral functions on quantum computers

We investigate quantum algorithms derived from tensor networks to simulate the static and dynamic properties of quantum many-body systems. Using a sequentially prepared quantum circuit representation of a matrix product state (MPS) that we call a quantum tensor network (QTN), we demonstrate algorithms to prepare ground and excited states on a quantum computer and apply them to molecular nanomagnets (MNMs) as a paradigmatic example. In this setting, we develop two approaches for extracting the spectral correlation functions measured in neutron-scattering experiments: (a) a generalization of the SWAP test for computing wave function overlaps and, (b) a generalization of the notion of matrix product operators to the QTN setting which generates a linear combination of unitaries. The latter method is discussed in detail for translationally invariant spin-half systems, where it is shown to reduce the qubit resource requirements compared with the SWAP method and may be generalized to other systems. We demonstrate the versatility of our approaches by simulating spin-1/2 and spin-3/2 MNMs, with the latter being an experimentally relevant model of a Cr$^{3+}_{8}$ ring. Here, our approach has qubit requirements that are independent of the number of constituents of the many-body system and scale only logarithmically with the bond dimension of the MPS representation, making them appealing for implementation on near-term quantum hardware with mid-circuit measurement and reset.

Neutron scattering↗

Petrophysical and geomechanical properties of gas hydrate-bearing sediments recovered from Alaska North Slope 2018 Hydrate-01 Stratigraphic Test Well

Knowledge of petrophysical and geomechanical properties of gas hydrate-bearing sediments are essential for predicting reservoir response due to gas production from gas hydrate reservoirs. That information will be critical parameters for designing production well completion such as specification of depressurization pump, water storage tank, and mesh size of sand screen. In December 2018, Stratigraphic Test Well Hydrate-01 was drilled in the western part of the Prudhoe Bay Unit, Alaska North Slope as part of the technical planning effort for a future long-term production test by collaborative team of DOE/NETL, USGS, and MH21-S (Boswell et al., 2020, Collett et al., 2020, Okinaka et al., 2020). Data set of logging-while-drilling (LWD) were acquired (Haines et al., 2020, Suzuki et al., 2019) and core sampling depth was determined on-site.Side-wall pressure coring was conducted to recover gas hydrate-bearing sediments from two reservoir sections named Unit-B and Unit-D. A total of 34 cores were successfully recovered by 5 runs of a wire-line deployed pressure corer (CoreVault® System - Halliburton). Core analysis scheme of this project are shown in Figure 1. All cores were quenched in liquid nitrogen while at high pressure in the laboratory of Stratum Reservoir, LLC. at Anchorage (Figure 1, a)). And the cores were removed from the pressure corer autoclave with temperature support of dry ice and stored under liquid nitrogen at atmospheric pressure. 19 damaged cores were processed to index property measurements directly such as grain size, grain density. 4 of another 17 cores were depressurized and trimmed for making a plug to analyze petrophysical properties of host sediments. Unsteady-state permeability test was conducted to obtain relative water permeability to gas and core scale NMR T2 distribution measurement was performed for evaluating pore size distribution at Houston (Figure 1, b)). Remained high quality 13 cores were preserved with gas hydrate for advanced laboratory analysis. National Institute of Advanced Industrial Science and Technology, as a part of the Japanese National Hydrate Research Program (MH21-S, funded by Ministry of Economy, Trade and Industry), received the samples at Sapporo, Japan for advanced core analysis. High-resolution X-ray CT was used to analyze the quality of the samples, which showed undisturbed lithological layers. Cores were lathed into cylindrical shape and distributed for multi property measurements (Figure 1, c)).<p>As a result, sediment from Unit-D is categorized as silty sand at ~37% porosity with ~80% gas hydrate saturation. An average hydration number n = 6.16 was measured by Raman spectroscopy. An average intrinsic permeability of ~400 mD and in situ effective permeability (with hydrate) on the order of ~10 mD. The Unit B recovered cores consisted of well sorted sand at ~40% porosity with ~95% gas hydrate saturation. An average intrinsic permeability of ~1 Darcy and in situ effective permeability on the order of ~30 mD was measured for the Unit B cores. There was a small permeability reduction due to porosity loss with increasing effective stress that simulated consolidation behavior along with depressurization in the highly permeable sandy sediment. The apparent minimum change in porosity and permeability may be caused by the low compressibility of quartz sand grains in the recovered cores. XRD and thermal conductivity analysis also suggested high quartz content. Triaxial compression tests established internal friction angles based on the Mohr-Coulomb's failure criterion, which are 40° for hydrate-bearing sediment and 29.8° for hydrate free sediment.</p>

Yoneda, Jun↗

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗

Para-particle oscillator simulations on a trapped-ion quantum computer

Deformed oscillators allow for a generalization of the standard fermions and bosons, namely, for the description of para-particles. Such particles remain hypothetical and unobserved in nature; yet, they can model physical phenomena, such as topological phases of matter. Here, we report the digital quantum simulation of para-particle oscillators by mapping para-particle states to the state of a qubit register, which allows us to identify the para-particle oscillator Hamiltonian as an XY model and further digitize the system onto a universal set of gates. In both instances, the gate depth grows polynomially with the number of qubits used. To establish the validity of our results, we experimentally simulate the dynamics of para-fermions and para-bosons, demonstrating full control of para-particle oscillators on a quantum computer. Furthermore, we compare the overall performance of the digital simulation of dynamics of the driven para-Fermi oscillator to a recent analog quantum simulation result.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗

Asymptotic errors in adiabatic evolution

The adiabatic theorem in quantum mechanics implies that if a system is in a discrete eigenstate of a Hamiltonian and the Hamiltonian evolves in time arbitrarily slowly, the system will remain in the corresponding eigenstate of the evolved Hamiltonian. Understanding corrections to the adiabatic result that arise when the evolution of the Hamiltonian is slow—but not arbitrarily slow—has become increasingly important, especially since adiabatic evolution has been proposed as a method of state preparation in quantum computing. Here, this paper identifies two regimes, an adiabatic regime in which corrections are generically small and can depend on details of the evolution throughout the path, and a hyperadiabatic regime in which the error is given by a form similar to an asymptotic expansion in the inverse of the evolution time with the coefficients depending principally on the behavior at the endpoints. However, the error in this hyperadiabatic regime is neither given by a true asymptotic series nor solely dependent on the endpoints: the coefficients combine the contributions from both endpoints, with relative phase factors that depend on the average spectral gaps along the trajectory, multiplied by the evolution time. The central result of this paper is to identify a quantity, referred to as the typical error, which is obtained by appropriately averaging the error over evolution times that are small compared to the evolution time itself. This typical error is characterized by an asymptotic series and depends solely on the endpoints of the evolution, remaining independent of the details of the intermediate evolution.

adiabatic approximation↗