Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic 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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Phase Transition and Liquid-like Superionic Conduction in Ag 2 S

Superionic conductivity in certain polymorphs of Ag 2 S has inspired numerous concepts for materials applications, but the relationship between the structure and the mobility of silver ions remains poorly explored. In this paper, we report ab initio molecular dynamics simulations for low- (acanthite) and high-temperature (argentite) Ag 2 S polymorphs that reveal the dynamical processes, giving rise to the superionic behavior in the latter. Similarities between their sulfur sublattices enable simulations of silver ion diffusivities and pathways on essentially an equal footing. For the higher temperature polymorph, calculated temperature-dependent mean square displacements and activation energies by the nudged elastic band method show good correspondence with expectations from the experiment. In the superionic state, silver atoms diffuse in a liquid-like behavior with no preferred diffusion pathways, within the relatively stable body-centered cubic sulfur framework. Conversely, conduction in acanthite appears to depend more on the mobilities of electronic charge carriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2023

The critical role of ultra-low-energy vibrations in the relaxation dynamics of molecular qubits

Improving the performance of molecular qubits is a fundamental milestone towards unleashing the power of molecular magnetism in the second quantum revolution. Taming spin relaxation and decoherence due to vibrations is crucial to reach this milestone, but this is hindered by our lack of understanding on the nature of vibrations and their coupling to spins. Here we propose a synergistic approach to study a prototypical molecular qubit. It combines inelastic X-ray scattering to measure phonon dispersions along the main symmetry directions of the crystal and spin dynamics simulations based on DFT. We show that the canonical Debye picture of lattice dynamics breaks down and that intra-molecular vibrations with very-low energies of 1-2 meV are largely responsible for spin relaxation up to ambient temperature. We identify the origin of these modes, thus providing a rationale for improving spin coherence. The power and flexibility of our approach open new avenues for the investigation of magnetic molecules with the potential of removing roadblocks toward their use in quantum devices.

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

Optimization of Solid Oxide Electrolysis Cell Systems Accounting for Long-Term Performance and Health Degradation

This study focuses on optimizing solid oxide electrolysis cell (SOEC) systems for efficient and durable long-term hydrogen (H2) production. While the elevated operating temperatures of SOECs offer advantages in terms of efficiency, they also lead to chemical degradation, which shortens cell lifespan. To address this challenge, dynamic degradation models are coupled with a steady-state, two-dimensional, non-isothermal SOEC model and steady-state auxiliary balance of plant equipment models, within the IDAES modeling and optimization framework. A quasi-steady state approach is presented to reduce model size and computational complexity. Long-term dynamic simulations at constant H2 production rate illustrate the thermal effects of chemical degradation. Dynamic optimization is used to minimize the lifetime cost of H2 production, accounting for SOEC replacement, operating, and energy expenses. Several optimized operating profiles are compared by calculating the Levelized Cost of Hydrogen (LCOH).

Giridhar, Nishant↗
DOE OSTI2024

Optimization of Solid Oxide Electrolysis Cell Systems Accounting for Long-Term Performance and Health Degradation

This study focuses on optimizing solid oxide electrolysis cell (SOEC) systems for efficient and durable long-term hydrogen (H2) production. While the elevated operating temperatures of SOECs offer advantages in terms of efficiency, they also lead to chemical degradation, which shortens cell lifespan. To address this challenge, dynamic degradation models are coupled with a steady-state, two-dimensional, non-isothermal SOEC model and steady-state auxiliary balance of plant equipment models, within the IDAES modeling and optimization framework. A quasi-steady state approach is presented to reduce model size and computational complexity. Long-term dynamic simulations at constant H2 production rate illustrate the thermal effects of chemical degradation. Dynamic optimization is used to minimize the lifetime cost of H2 production, accounting for SOEC replacement, operating, and energy expenses. Several optimized operating profiles are compared by calculating the Levelized Cost of Hydrogen (LCOH).

Giridhar, Nishant↗
DOE OSTI2021

Direct Observation of Ultrafast Hydrogen Bond Strengthening in Liquid Water

Water is one of the most important, yet least understood, liquids in nature. Many anomalous properties of liquid water originate from its well-connected hydrogen bond network, including unusually efficient vibrational energy redistribution and relaxation. An accurate description of the ultrafast vibrational motion of water molecules is essential for understanding the nature of hydrogen bonds and many solution-phase chemical reactions. Most existing knowledge of vibrational relaxation in water is built upon ultrafast spectroscopy experiments. However, these experiments cannot directly resolve the motion of the atomic positions and require difficult translation of spectral dynamics into hydrogen bond dynamics. Here we measured the ultrafast structural response to the excitation of the OH stretching vibration in liquid water with femtosecond temporal and atomic spatial resolution using liquid ultrafast electron scattering. We observed a transient hydrogen bond contraction of roughly 0.04 Å on a timescale of 80 femtoseconds, followed by a thermalization on a timescale of ~1 picosecond. Molecular dynamics simulations reveal the need to treat the distribution of the shared proton in the hydrogen bond quantum mechanically to capture the structural dynamics on the femtosecond timescales. Our experiment and simulations unveiled the intermolecular character of the water vibration preceding the relaxation of the OH stretch.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SynopFrame: Multiscale time-dependent visual abstraction framework for analyzing DNA nanotechnology simulations

We present an open-source framework, SynopFrame, that allows DNA nanotechnology (DNA-nano) experts to analyze and understand molecular dynamics simulation trajectories of their designs. We use a multiscale multi-dimensional abstraction space, connect the representations to a projected conformational space plot of the structure’s temporal sequence, and thus enable experts to analyze the dynamics of their structural designs and, specifically, failure cases of the assembly. In addition, our time-dependent abstraction representation allows the biologists, for the first time in a smooth and structurally clear way, to identify and observe temporal transitions of a DNA-nano design from one configuration to another, and to highlight important periods of the simulation for further analysis. We realize SynopFrame as a dashboard of the different synchronized 3D spatial and 2D schematic visual representations, with a color overlay to show essential properties such as the status of hydrogen bonds. The linking of the spatial, schematic, and abstract views ensures that users can effectively analyze the high-frequency motion. We also categorize the status of the hydrogen bonds into a new format to allow us to color-encode it and overlay it on the representations. To demonstrate the utility of SynopFrame, we describe example usage scenarios and report user feedback.

Abstraction space↗
DOE OSTI2023

Composition-dependent fracture energy in metallic glasses

The interplay between metallic glasses (MGs) mechanical properties, fracture energy (G), and glass-forming ability (GFA) and their dependence on alloy composition remain poorly understood. Here, we perform molecular dynamics simulations to investigate the intrinsic composition-dependence of G in Cu x Zr 100-x MGs (x = 20, 30, 40, 50, 64). The results indicate that the value of G increases with Cu content. In addition, it is revealed that MGs with higher G values display higher Poisson’s ratio (ν) and GFA, suggesting a close correlation between fracture toughness, mechanical properties, and GFA. This correlation between G, ν and GFA can be understood based on the fragility (m) of supercooled liquids, which is directly related to the structural heterogeneity in MGs. Larger m values are related to dynamic slowdown and supercooled liquid stabilization, which enhance GFA and the formation of pronounced structural heterogeneity, comprised of loosely packed regions that favor β-relaxation and the activation of shear transformation zones. Those concurrently promote the expansion of the plastic zone at the crack tip, enhancing the observed value of G. Furthermore, these simulation results shed light on the intrinsic relationship between fracture toughness, mechanical properties, and alloy composition in MGs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Ab initio molecular dynamics simulations combined with enhanced sampling techniques are becoming widespread methods to investigate chemical phenomena in catalytic systems. These techniques automatically include finite temperature effects, anharmonicity, and collective dynamics in their robust description of enthalpic and entropic contributions, which can have significant impact on reaction free energy landscapes. This contrasts with standard ab initio static approaches that are based on assessing reaction free energies from various coarse-grained descriptions of the reaction potential energy surface. Enhanced sampling ab initio molecular dynamics opens the way to first principles simulations of systems of increasing complexity like solid/liquid catalytic interfaces. Here, we aim at guiding the reader through the basis of these techniques, summarizing their fundamental theoretical and practical aspects, and reviewing the relevant literature in the field. After a brief introduction to the problem, we will illustrate the advantage of using molecular simulations to include finite temperature effects, examine the most common ab initio techniques currently in use, describe their application to solid state heterogeneous catalysts, and finally critically review the most popular enhanced sampling techniques used in computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2024

Molecular Dynamics Study of the Effect of Grafting Density on Ion Diffusivity in a MARTINI Coarse-Grained Strong Polyelectrolyte Brush

Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separation technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact the PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and the counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI force field. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N = 150 monomers; M W = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm 2 . Using 5 μs simulations, we investigate the effects of grafting density on the PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. Finally, the interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sliceable, Moldable, and Highly Conductive Electrolytes for All-Solid-State Batteries

All-solid-state batteries (ASSBs) require solid electrolytes with high ionic conductivity, stability, and deformability for optimal energy and power density. Here, we developed lithium-deficient lithium yttrium bromide (LYB) solid electrolytes, Li 3–x YBr 6–x (0 ≤ x ≤ 0.50), using a comelting method with controlled lithium deficiency. These electrolytes exhibit favorable mechanical properties such as high moldability and sliceability. The Li 2.65 YBr 5.65 composition has an ionic conductivity of 4.49 mS cm –1 at 25 °C and an activation energy of 0.28 eV. Compared to Li 3 YBr 6 , Li 2.65 YBr 5.65 demonstrates improved rate performance and cycling stability in ASSBs. High-resolution X-ray diffraction confirms the formation of the LYB phase with a C2/m space group. Structural analysis reveals increased cation disorder and larger polyhedral volumes for x > 0 in Li 3–x YBr 6–x , contributing to reduced Li + migration energy barriers. Bond valence site energy calculations and molecular dynamics simulations reveal enhanced 3D lithium-ion transport. NMR spectroscopy further highlights increased Li + dynamics and impurity elimination.

Poudel, Tej P. [Florida State Univ., Tallahassee, ↗
DOE OSTI2025

Ultrafast Events in Electrocyclic Ring-Opening Reactions

Electrocyclic reactions are characterized by the concerted formation and cleavage of multiple σ and π bonds in a molecular system and have been extensively studied since they were introduced by Robert Burns Woodward and Roald Hoffmann in 1965. Recent advances and the integration of time-resolved experiments and nonadiabatic quantum molecular dynamics simulations have transformed the traditional understanding of electrocyclic reactions beyond the Woodward–Hoffmann rules. In this review, we focus on recent studies of 1,3-cyclohexadiene and two of its derivatives, α-phellandrene and α-terpinene, to shed light on the underlying mechanisms of electrocyclic photochemical reactions. We highlight recent progress in ultrafast electron diffraction techniques and the simulation approach of ab initio multiple spawning. Together, these approaches can elucidate molecular structure dynamics from femtosecond to picosecond timescales as well as nuclear and electronic responses at conical intersections.

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

Glovebox Fires Phase 1: Simulations of Open Burner Experiments

Understanding the duration of a fire required to breach a glovebox via glove ports and/or windows within a facility is important in developing mitigation strategies for safety. To facilitate this understanding, both simulation and experiments are utilized. The approach is to first validate the Computation Fluid Dynamics (CFD) code, Fire Dynamics Simulator (FDS), with experimental data collected by New Mexico Tech (NMT) and then perform simulations of full-scale rooms containing gloveboxes to assess numerous scenarios which would otherwise be cost-prohibitive experimentally. This report provides comparison to first-phase experiments involving a fire without a glovebox. A subsequent report will provide comparison to experiments involving a glovebox. The results indicate good agreement with FDS tending to over predict the pre-mixed and diffusion-mode tests by 6% and 10%, respectively.

42 ENGINEERING↗
DOE OSTI2021

Biomembrane Structure and Material Properties Studied With Neutron Scattering

Cell membranes and their associated structures are dynamical supramolecular structures where different physiological processes take place. Detailed knowledge of their static and dynamic structures is therefore needed, to better understand membrane biology. The structure–function relationship is a basic tenet in biology and has been pursued using a range of different experimental approaches. In this review, we will discuss one approach, namely the use of neutron scattering techniques as applied, primarily, to model membrane systems composed of lipid bilayers. An advantage of neutron scattering, compared to other scattering techniques, is the differential sensitivity of neutrons to isotopes of hydrogen and, as a result, the relative ease of altering sample contrast by substituting protium for deuterium. This property makes neutrons an ideal probe for the study of hydrogen-rich materials, such as biomembranes. In this review article, we describe isotopic labeling studies of model and viable membranes, and discuss novel applications of neutron contrast variation in order to gain unique insights into the structure, dynamics, and molecular interactions of biological membranes. We specifically focus on how small-angle neutron scattering data is modeled using different contrast data and molecular dynamics simulations. We also briefly discuss neutron reflectometry and present a few recent advances that have taken place in neutron spin echo spectroscopy studies and the unique membrane mechanical data that can be derived from them, primarily due to new models used to fit the data.

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

3D structure of anisotropic flow in small collision systems at energies available at the BNL Relativistic Heavy Ion Collider

Here, we present (3 + 1)-dimensional [(3 + 1)D] dynamical simulations of asymmetric nuclear collisions at the BNL Relativistic Heavy Ion Collider (RHIC). Employing a dynamical initial state model coupled to (3 + 1)D viscous relativistic hydrodynamics, we explore the rapidity dependence of anisotropic flow in the RHIC small system scan at 200 GeV center-of-mass energy. We calibrate parameters to describe central 3 He + Au collisions and make extrapolations to d + Au and p + Au collisions. Our calculations demonstrate that approximately 50% of the v 3 (p T ) difference between the measurements by the STAR and PHENIX Collaborations can be explained by the use of reference flow vectors from different rapidity regions. This emphasizes the importance of longitudinal flow decorrelation for anisotropic flow measurements in asymmetric nuclear collisions, and the need for (3 + 1)D simulations. We also present results for the beam energy dependence of particle spectra and anisotropic flow in d + Au collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Simulating non-native cubic interactions on noisy quantum machines

As a milestone for general-purpose computing machines, we demonstrate that quantum processors can be programed to efficiently simulate dynamics that are not native to the hardware. Moreover, on noisy devices without error correction, we show that simulation results are significantly improved when the quantum program is compiled using modular gates instead of a restricted set of standard gates. We demonstrate the general methodology by solving a cubic interaction problem, which appears in nonlinear optics, gauge theories, as well as plasma and fluid dynamics. To encode the non-native Hamiltonian evolution, we decompose the Hilbert space into a direct sum of invariant subspaces in which the nonlinear problem is mapped to a finite-dimensional Hamiltonian simulation problem. Furthermore, in a three-states example, the resultant unitary evolution is realized by a product of approximately 20 standard gates, using which approximately ten simulation steps can be carried out on state-of-the-art quantum hardware before results are corrupted by decoherence. In comparison, the simulation depth is improved by more than an order of magnitude when the unitary evolution is realized as a single cubic gate, which is compiled directly using optimal control. Alternatively, parametric gates may also be compiled by interpolating control pulses. Modular gates thus obtained provide high-fidelity building blocks for quantum Hamiltonian simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗