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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structural origin of disorder-induced ion conduction in NaFePO 4 cathode materials

Diffusion in NaFePO 4 can be enhanced through amorphization. Based on computations using DFT and machine learning potentials, we ascribe this phenomenon to the formation of less constrained Na-ion environments upon disordering. Most modern battery technologies depend on solid-state crystalline cathode materials. However, some of these materials are constrained by the low ionic conductivity of their most stable phases. An example of this is maricite (NaFePO 4 ). Interestingly, experiments have shown that maricite can improve its rate capability through disordering (amorphization). However, experimental characterization of amorphous cathode materials remains a major challenge, hindering a clear understanding of the structural origin of the disorder-induced improvement in sodium-ion mobility. To address this, we here employ molecular dynamics simulations by first training a machine learning potential for NaFePO 4 based on the atomic cluster expansion approach and a batch active learning potential parameterization scheme. This potential is then applied to explore the structural and dynamical properties of NaFePO 4 glasses as cathode materials. Specifically, we investigate the effect of glass structure on sodium-ion diffusion, revealing the relative influences of short-range and medium-range order features. We find significant heterogeneity in sodium-ion diffusivity in the glass, with fast-conducting ions residing in less constrained atomic environments with fewer P and Fe neighbors. These more mobile ions are also surrounded by larger ring-type structures. Overall, the results and developed approach present promising avenues for developing high-performance glassy cathodes for next-generation batteries.

Christensen, Rasmus↗
DOE OSTI2026

Rotational memory function of SPC/E water

Memory effects are essential for the dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for water have never been calculated directly from molecular dynamics simulations. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The normalized memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, ≲1 fs, justifies the use of the mathematics of rotational diffusion to describe the dynamics of a single molecular dipole moment in bulk water. Here, an analytical equation for the rotational memory time is derived.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
DOE OSTI2020

Dislocation dynamics in polycrystalline materials

An extension of the dislocation dynamics method to polycrystalline materials is presented with emphasis on handling interactions between dislocations and grain boundaries. The key new features of the method include: (1) a grain boundary detection algorithm, (2) new dislocation emission criteria based on power dissipation and area growth, (3) mobility of dislocations at grain boundaries, and (4) extension of topological operations e.g. mesh adaption, collision, and dissociation to handle dislocations interacting with grain boundaries. Furthermore, verification and validation examples are executed to compare our calculations with existing analytic methods, molecular dynamics simulations, and experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A mesoscopic digital twin that bridges length and time scales for control of additively manufactured metal microstructures

We present our recent development of an integrated mesoscale digital twin (DT) framework for relating processing conditions, microstructures, and mechanical responses of additively manufactured (AM) metals. In particular, focusing on the laser powder bed fusion technique, we describe how individual modeling and simulation capabilities are coupled to investigate and control AM microstructural features at multiple length and time scales. We review our prior case studies that demonstrate the integrated modeling schemes, in which high-fidelity melt pool dynamics simulations provide accurate local thermal profiles and histories to subsequent AM microstructure simulations. We also report our new mechanical response modeling results for predicted AM microstructures. In addition, we illustrate how our DT framework has been validated through modeling–experiment integration, as well as how it has been practically utilized to guide and analyze AM experiments. Finally, we share our perspectives on future directions of further development of the DT framework for more efficient, accurate predictions and wider ranges of applications.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Internet Data Delivery for Future Space Missions

Ongoing work at National Aeronautics and Space Administration Goddard Space Flight Center (NASA/GSFC), seeks to apply standard Internet applications and protocols to meet the technology challenge of future satellite missions. Internet protocols and technologies are under study as a future means to provide seamless dynamic communication among heterogeneous instruments, spacecraft, ground stations, constellations of spacecraft, and science investigators. The primary objective is to design and demonstrate in the laboratory the automated end-to-end transport of files in a simulated dynamic space environment using off-the-shelf, low-cost, commodity-level standard applications and protocols. The demonstrated functions and capabilities will become increasingly significant in the years to come as both earth and space science missions fly more sensors and as the need increases for more network-oriented mission operations. Another element of increasing significance will be the increased cost effectiveness of designing, building, integrating, and operating instruments and spacecraft that will come to the fore as more missions take up the approach of using commodity-level standard communications technologies. This paper describes how an IP (Internet Protocol)-based communication architecture can support all existing operations concepts and how it will enable some new and complex communication and science concepts. The authors identify specific end-to-end data flows from the instruments to the control centers and scientists, and then describe how each data flow can be supported using standard Internet protocols and applications. The scenarios include normal data downlink and command uplink as well as recovery scenarios for both onboard and ground failures. The scenarios are based on an Earth orbiting spacecraft with downlink data rates from 300 Kbps to 4 Mbps. Included examples are based on designs currently being investigated for potential use by the Global Precipitation Measurement (GPM) mission.

Rash, James↗
NASA NTRS2002

Internet Technology for Future Space Missions

Ongoing work at National Aeronautics and Space Administration Goddard Space Flight Center (NASA/GSFC), seeks to apply standard Internet applications and protocols to meet the technology challenge of future satellite missions. Internet protocols and technologies are under study as a future means to provide seamless dynamic communication among heterogeneous instruments, spacecraft, ground stations, constellations of spacecraft, and science investigators. The primary objective is to design and demonstrate in the laboratory the automated end-to-end transport of files in a simulated dynamic space environment using off-the-shelf, low-cost, commodity-level standard applications and protocols. The demonstrated functions and capabilities will become increasingly significant in the years to come as both earth and space science missions fly more sensors and the present labor-intensive, mission-specific techniques for processing and routing data become prohibitively. This paper describes how an IP-based communication architecture can support all existing operations concepts and how it will enable some new and complex communication and science concepts. The authors identify specific end-to-end data flows from the instruments to the control centers and scientists, and then describe how each data flow can be supported using standard Internet protocols and applications. The scenarios include normal data downlink and command uplink as well as recovery scenarios for both onboard and ground failures. The scenarios are based on an Earth orbiting spacecraft with downlink data rates from 300 Kbps to 4 Mbps. Included examples are based on designs currently being investigated for potential use by the Global Precipitation Measurement (GPM) mission.

Hennessy, Joseph F.↗
DOE OSTI2023

Modeling of the amplified optical stochastic cooling experiment at IOTA

Optical Stochastic Cooling (OSC), a beam cooling technique based on Stochastic Cooling, is in the early stages of experimental development. It uses radiation produced by the beam in an undulator magnet (the pickup) to correct the momentum deviation of particles downstream in another undulator (the kicker). The OSC mechanism was recently demonstrated at Fermilab’s IOTA ring using a passive configuration. However, the cooling rate of OSC can be dramatically increased by first amplifying the undulator radiation before applying the corrective kick. In collaboration with the IOTA experiment, we developed a computational model of the OSC mechanism. This paper presents beam-dynamics simulations of the amplified-OSC configuration. We implement a model of intrabeam scattering and study the effects on beam equilibrium and diffusion rate as a function of bunch charge. Finally, we investigate the phase-space dynamics with various coupling configurations between the transverse and longitudinal planes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Inter-Subunit Dynamics Controls Tunnel Formation During the Oxygenation Process in Hemocyanin Hexamers

Hemocyanin from horseshoe crab in its active form is a homo-hexameric protein. It exists in open and closed conformations when transitioning between deoxygenated and oxygenated states. Here, we present a detailed dynamic atomistic investigation of the oxygenated and deoxygenated states of the hexameric hemocyanin using explicit solvent molecular dynamics simulations. We focus on the variation in solvent cavities and the formation of tunnels in the two conformational states. By employing principal component analysis and CVAE-based deep learning, we are able to differentiate between the dynamics of the deoxy- and oxygenated states of hemocyanin. Finally, our results identify the deoxygenated open conformation, which adopts a stable, closed conformation after the oxygenation process.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1974

Specification and correlation of the sine vibration environment for Viking '75

Two Viking spacecraft will be individually launched on a new Titan IIIE/Centaur D-1T launch vehicle in August 1975. The method for the establishment of spacecraft sine vibration test levels prior to availability of any Titan IIIE/Centaur D-1T flight data by use of both computer simulations and data from previous Titan and Atlas Centaur vehicles is described. The specification level is compared with actual flight data obtained from a proof flight launch of the Titan IIIE/Centaur D-1T and a Viking dynamic simulator in January 1974. An objective of the proof flight launch was to obtain estimates of the flight loads and environments. The criteria used to minimize the structural weight that would result from an unmodified application of a sine test environment are described.

Snyder, R. E.↗
DOE OSTI2023

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics

In this article, strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2023

Non-Empirical and Self-Interaction Corrections for DFTB: Towards Accurate Quantum Simulations for Large Mesoscale Systems (Final Report)

This project was comprised of two complementary (but parallel) thrusts: (1) implementing massively-parallelized computing hardware (with new computational hardware that may replace GPUs) in the density functional tight binding (DFTB) approach and (2) developing new capabilities in DFTB to calculate the electronic structure and dynamics of large chemical systems. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of chemical systems at the quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum mechanical nature of chemical systems; however, these methods cannot tackle the large sizes and length scales relevant to dynamics simulations of realistic systems. The DFTB formalism utilized in this project provides a viable approach for probing these large systems at a quantum mechanical level of detail. However, to utilize the DFTB approach for accurate calculations of electronic properties, it is crucial to incorporate quantum-based corrections in DFTB since exchange-correlation effects can still remain very strong in these large systems. At the same time, enhancing the computational efficiency of DFTB is also essential since optimal computational performance is required for addressing the large size scales associated with realistic chemical systems. As such, the new non-empirical corrections and computing hardware enhancements implemented in this project will enable accurate and computationally efficient approaches to directly probe electronic properties in these large, complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From 2D to 4D: a containerized workflow and browser to explore dynamic chromatin architecture

Background Characterizing the physical organization of the genome is essential for understanding long-range gene regulation, chromatin compartmentalization, and epigenetic accessibility. Hi-C experiments generate two-dimensional (2D) genome-wide contact maps of chromatin interactions by capturing the spatial proximity between genomic loci, which reveal interaction frequencies but lack the spatial resolution needed to interpret the three-dimensional (3D) genome structure(s). Emerging evidence suggests that epigenetic regulation is closely linked to 3D genome architecture, and that structural changes over time (4D) drive key biological processes in development, disease, and environmental response. Thus, integrating 3D structure with functional data is critical for a more complete understanding of genome regulation. Previous work, most notably the 4DHiC chromosome modeling framework, has shown that physical multi-dimensional modeling approaches rooted in polymer physics and molecular dynamics can resolve these structures at biologically meaningful resolutions by integrating temporal Hi-C data with physical constraints to uncover dynamic chromosome reorganization. Thus, molecular dynamics simulations, constrained by Hi-C contact matrices, can resolve fine-scale structural changes and reveal functionally significant transitions in chromatin conformation. Results Herein, we present the 4D Genome Browser Workflow (4DGBWorkflow) and the 4D Genome Browser (4DGB). The algorithm is based on the 4DHiC method, and the containerized tool is an end-to-end workflow that can transform, filter, and view 4D epigenomics and chromatin datasets, allowing non-specialists to apply three-dimensional modeling principles to diverse datasets and experimental conditions. The software executes on a laptop running macOS, Linux or Windows. From input Hi-C files (.hic), the 4DGBWorkflow produces 3D reconstructions of chromosomes, integrates the reconstruction with track data (e.g., epigenetic marks, transcriptome profiles), and provides comparative visualization of the results in a single workflow. Conclusions The 4DGBWorkflow and 4D Genome Browser are open-source tools for comparative analysis and visualization of 4D chromosome datasets, including chromatin architecture and epigenomic signals. Automatic integration of Hi-C data with molecular dynamics democratizes the construction of time resolved 3D genome structures, simplifying complex simulations and data integration schemes.

3D Genome Browser↗
DOE OSTI2024

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗
DOE OSTI2024

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Learning stochastic dynamics and predicting emergent behavior using transformers

We show that a neural network originally designed for language processing can learn the dynamical rules of a stochastic system by observation of a single dynamical trajectory of the system, and can accurately predict its emergent behavior under conditions not observed during training. We consider a lattice model of active matter undergoing continuous-time Monte Carlo dynamics, simulated at a density at which its steady state comprises small, dispersed clusters. We train a neural network called a transformer on a single trajectory of the model. The transformer, which we show has the capacity to represent dynamical rules that are numerous and nonlocal, learns that the dynamics of this model consists of a small number of processes. Forward-propagated trajectories of the trained transformer, at densities not encountered during training, exhibit motility-induced phase separation and so predict the existence of a nonequilibrium phase transition. Transformers have the flexibility to learn dynamical rules from observation without explicit enumeration of rates or coarse-graining of configuration space, and so the procedure used here can be applied to a wide range of physical systems, including those with large and complex dynamical generators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2024

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α–CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Chemistry↗