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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Ionic Polymers Under Dynamic Conditions

Combining neutron scattering tools and exascale molecular dynamics simulations, research was carried out to resolve the structure and motion of ionic polymers under dynamic conditions. Ionizable polymers constitute a key component in current and potential transformative lightweight energy generation and storage devices. However, the inherent segregation of these polymers to ionic clusters and non-ionic domains, or clusters, that on one hand enables their function, confines the dynamics of the polymer chains, leading to the formation of far-from equilibrium states. These non-equilibrium states morph under in situ conditions, affecting the performance-longevity-safety-cost foundation that underlines the use of these polymers. On route to controlling the assemblies of ionizable polymers into stable functional materials, the research focused on a) understanding the inherent confinement of the polymer chains by the ionic clusters and b) resolving of the response of these polymers to flow fields that are integral to the processing and function of these macromolecules.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

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

Self-Healing Ceramics and Mechanical Metamaterials

We performed large scale molecular dynamics (MD) simulations and reactive & quantum molecular dynamics simulations to investigate: (1) faceting, grain growth, and crack healing in ceramics; (2) stress-release mechanism and healing in fibers; and (3) mechanical properties of two-dimensional metamaterials. We also designed highly accurate interatomic potentials and forces using artificial neural networks (ANN) and first-principles molecular dynamics simulation data for use in our MD and reactive MD simulations. We demonstrated the efficacy of modeling quantum materials atomistic simulation data using unsupervised deep learning approaches on a quantum computer. These achievements are elaborated in the following sections.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation of polymerization induced phase separation in model thermosets

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T g ) and chain length and thus have the potential for phase separation. Here, one crosslinker has a high T g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T g . A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

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

Acclimation of Phytoplankton Fe:C Ratios Dampens the Biogeochemical Response to Varying Atmospheric Deposition of Soluble Iron

Abstract Dissolved iron (dFe) plays an important role in regulating marine productivity. In high nutrient, low chlorophyll regions (>33% of the global ocean), iron is the primary growth limiting nutrient, and elsewhere iron can regulate nitrogen fixation by diazotrophs. The link between iron availability and carbon export is strongly dependent on the phytoplankton iron quotas or cellular Fe:C ratios. This ratio varies by more than an order of magnitude in the open ocean and is positively correlated with ambient dFe concentrations in field observations. Representing Fe:C ratios within models is necessary to investigate how ocean carbon cycling will interact with perturbations to iron cycling in a changing climate. The Community Earth System Model ocean component was modified to simulate dynamic, group‐specific, phytoplankton Fe:C that varies as a function of ambient iron concentration. The simulated Fe:C ratios improve the representation of the spatial trends in the observed Fe:C ratios. The acclimation of phytoplankton Fe:C ratios dampens the biogeochemical response to varying atmospheric deposition of soluble iron, compared to a fixed Fe:C ratio. However, varying atmospheric soluble iron supply has first order impacts on global carbon and nitrogen fluxes and on nutrient limitation spatial patterns. Our results suggest that pyrogenic Fe is a significant dFe source that rivals mineral dust inputs in some regions. Changes in dust flux and iron combustion sources (anthropogenic and wildfires) will modify atmospheric Fe inputs in the future. Accounting for dynamic phytoplankton iron quotas is critical for understanding ocean biogeochemistry and projecting its response to variations in atmospheric deposition.

Environmental Sciences & Ecology↗
DOE OSTI2021

X-composer: enabling cross-environments in-situ workflows between HPC and cloud

As large-scale scientific simulations and big data analyses become more popular, it is increasingly more expensive to store huge amounts of raw simulation results to perform post-analysis. To minimize the expensive data I/O, "in-situ" analysis is a promising approach, where data analysis applications analyze the simulation generated data on the fly without storing it first. However, it is challenging to organize, transform, and transport data at scales between two semantically different ecosystems due to the distinct software and hardware difference. To tackle these challenges, we design and implement the X-Composer framework. X-Composer connects cross-ecosystem applications to form an "in-situ" scientific workflow, and provides a unified approach and recipe for supporting such hybrid in-situ workflows on distributed heterogeneous resources. X-Composer reorganizes simulation data as continuous data streams and feeds them seamlessly into the Cloud-based stream processing services to minimize I/O overheads. For evaluation, we use X-Composer to set up and execute a cross-ecosystem workflow, which consists of a parallel Computational Fluid Dynamics simulation running on HPC, and a distributed Dynamic Mode Decomposition analysis application running on Cloud. Our experimental results show that X-Composer can seamlessly couple HPC and Big Data jobs in their own native environments, achieve good scalability, and provide high-fidelity analytics for ongoing simulations in real-time.

Wang, Dali↗
DOE OSTI2021

Berni Julian Alder, theoretical physicist and inventor of molecular dynamics, 1925–2020

Berni Julian Alder, one of the leading figures in the invention of molecular dynamics simulations used for a wide array of problems in physics and chemistry, died on September 7th, 2020. His career, spanning more than 65 years, transformed statistical mechanics, many body physics, the study of chemistry and the microscopic dynamics of fluids, by making atomistic computational simulation (in parallel with traditional theory and experiment) a new pathway to unexpected discoveries. Among his many honors, the CECAM prize, recognizing exceptional contributions to the simulation of the microscopic properties of matter is named for him. He was awarded the National Medal of Science by President Obama in 2008.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Amplifying Nanoparticle Reinforcement through Low Volume Topologically Controlled Chemical Coupling

Here, we present a streamlined method to covalently bond hydroxylated carbon nanotubes (CNOH) within a polyphenol matrix, all achieved through a direct, solvent-free process. Employing an extremely small concentration of CNOH (0.01% w/w) along with topologically contrasting linkers led to a maximum of 5-fold increase in modulus and a 25% enhancement in tensile strength compared to the unaltered matrix, an order of magnitude greater reinforcement (w/w) compared to state-of-the-art melt-processed nanocomposites. Through dynamic mechanical analysis, low field solid-state nuclear magnetic resonance spectroscopy, and molecular dynamics simulations, we uncovered the profound influence of linker’s conformational degrees of freedom on the segmental dynamics and therefore the material’s properties.

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

Ligand modulation of the conformational dynamics of the A 2A adenosine receptor revealed by single-molecule fluorescence

G protein-coupled receptors (GPCRs) are the largest class of transmembrane proteins, making them an important target for therapeutics. Activation of these receptors is modulated by orthosteric ligands, which stabilize one or several states within a complex conformational ensemble. The intra- and inter-state dynamics, however, is not well documented. Here, we used single-molecule fluorescence to measure ligand-modulated conformational dynamics of the adenosine A 2A receptor (A 2A R) on nanosecond to millisecond timescales. Experiments were performed on detergent-purified A 2R in either the ligand-free (apo) state, or when bound to an inverse, partial or full agonist ligand. Single-molecule Förster resonance energy transfer (smFRET) was performed on detergent-solubilized A 2A R to resolve active and inactive states via the separation between transmembrane (TM) helices 4 and 6. The ligand-dependent changes of the smFRET distributions are consistent with conformational selection and with inter-state exchange lifetimes ≥ 3 ms. Local conformational dynamics around residue 229 6.31 on TM6 was measured using fluorescence correlation spectroscopy (FCS), which captures dynamic quenching due to photoinduced electron transfer (PET) between a covalently-attached dye and proximal aromatic residues. Global analysis of PET-FCS data revealed fast (150–350 ns), intermediate (50–60 μs) and slow (200–300 μs) conformational dynamics in A 2A R, with lifetimes and amplitudes modulated by ligands and a G-protein mimetic (mini-G s ). Most notably, the agonist binding and the coupling to mini-G s accelerates and increases the relative contribution of the sub-microsecond phase. Molecular dynamics simulations identified three tyrosine residues (Y112, Y288 7.53 , and Y290 7.55 ) as being responsible for the dynamic quenching observed by PET-FCS and revealed associated helical motions around residue 229 6.31 on TM6. As a result, this study provides a quantitative description of conformational dynamics in A 2A R and supports the idea that ligands bias not only GPCR conformations but also the dynamics within and between distinct conformational states of the receptor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Fast Na diffusion and anharmonic phonon dynamics in superionic Na 3 PS 4

The design of new solid electrolytes (SEs) hinges on identifying and tuning relevant descriptors. Phonons describe the atomic dynamics in crystalline materials and provide a basis to encode possible minimum energy pathways for ion migration, but anharmonic effects can be large in SEs. Identifying and controlling the pertinent phonon modes coupled most strongly with ionic conductivity, and assessing the role of anharmonicity, could therefore pave the way for discovering and designing new SEs via phonon engineering. In this work, we investigate phonons in Na 3 PS 4 and their coupling to fast Na diffusion, using a combination of neutron scattering, ab initio molecular dynamics (AIMD), and extended molecular dynamics based on machine-learned potentials. We identify that anharmonic soft-modes at the Brillouin zone boundary of the anharmonically stabilized cubic phase constitute key phonon modes that control the Na diffusion process in Na 3 PS 4 . We demonstrate how these strongly anharmonic phonon modes enable Na-ions to hop along the minimum energy pathways. Further, the quasi-elastic neutron scattering (QENS) measurements, supplemented with ab initio and machine-learned molecular dynamics simulations, probe the Na diffusivity and diffusion mechanism. These results offer detailed microscopic insights into the dynamic mechanism of fast Na diffusion and provide an avenue to search for further Na solid electrolytes.

42 ENGINEERING↗
DOE OSTI2021

Effects of interaction strength of associating groups on linear and star polymer dynamics

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. In this study, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead–spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 k B T. We find that, for all melts and blends, interaction of a few k B T between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few k B T and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Improving the electrostatic design of the Jefferson Lab 300 kV DC photogun

The 300 kV DC high voltage photogun at Jefferson Lab was redesigned to deliver electron beams with a much higher bunch charge and improved beam properties. The original design provided only a modest longitudinal electric field ( E z ) at the photocathode, which limited the achievable extracted bunch charge. To reach the bunch charge goal of approximately few nC with 75 ps full-width at half-maximum Gaussian laser pulse width, the existing DC high voltage photogun electrodes and anode–cathode gap were modified to increase E z at the photocathode. In addition, the anode aperture was spatially shifted with respect to the beamline longitudinal axis to minimize the beam deflection introduced by the non-symmetric nature of the inverted insulator photogun design. We present the electrostatic design of the original photogun and the modified photogun and beam dynamics simulations that predict vastly improved performance. We also quantify the impact of the photocathode recess on beam quality, where recess describes the actual location of the photocathode inside the photogun cathode electrode relative to the intended location. A photocathode unintentionally recessed/misplaced by sub-millimeter distance can significantly impact the downstream beam size.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2021

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗