Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamics 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Sorption, Structure and Dynamics of CO2 and Ethane in Silicalite at High Pressure: A Combined Monte Carlo and Molecular Dynamics Simulation Study

Silicalite is an important nanoporous material that finds applications in several industries, including gas separation and catalysis. While the sorption, structure, and dynamics of several molecules confined in the pores of silicalite have been reported, most of these studies have been restricted to low pressures. Here we report a comparative study of sorption, structure, and dynamics of CO2 and ethane in silicalite at high pressures (up to 100 bar) using a combination of Monte Carlo (MC) and molecular dynamics (MD) simulations. The behavior of the two fluids is studied in terms of the simulated sorption isotherms, the positional and orientational distribution of sorbed molecules in silicalite, and their translational diffusion, vibrational spectra, and rotational motion. Both CO2 and ethane are found to exhibit orientational ordering in silicalite pores; however, at high pressures, while CO2 prefers to reside in the channel intersections, ethane molecules reside mostly in the sinusoidal channels. While CO2 exhibits a higher self-diffusion coefficient than ethane at low pressures, at high pressures, it becomes slower than ethane. Both CO2 and ethane exhibit rotational motion at two time scales. At both time scales, the rotational motion of ethane is faster. The differences observed here in the behavior of CO2 and ethane in silicalite pores can be seen as a consequence of an interplay of the kinetic diameter of the two molecules and the quadrupole moment of CO2.

Gautam, Siddharth (ORCID:0000000314435382)↗
NASA NTRS1997

Comparison of Chain Conformation of Poly(vinyl alcohol) in Solutions and Melts from Quantum Chemistry Based Molecular Dynamics Simulations

Confirmations of 2,4-dihydroxypentane (DHP), a model molecule for poly(vinyl alcohol), have been studied by quantum chemistry (QC) calculations and molecular dynamics (MD) simulations. QC calculations at the 6-311G MP2 level show the meso tt conformer to be lowest in energy followed by the racemic tg, due to intramolecular hydrogen bond between the hydroxy groups. The Dreiding force field has been modified to reproduce the QC conformer energies for DHP. MD simulations using this force field have been carried out for DHP molecules in the gas phase, melt, and CHCl3 and water solutions. Extensive intramolecular hydrogen bonding is observed for the gas phase and CHCl3 solution, but not for the melt or aqueous solution, Such a condensed phase effect due to intermolecular interactions results in a drastic change in chain conformations, in agreement with experiments.

Jaffe, Richard↗
DOE OSTI2020

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

All-atom molecular dynamics simulations of weak polyionic brushes: influence of charge density on the properties of polyelectrolyte chains, brush-supported counterions, and water molecules

All atom molecular dynamics (MD) simulations of planar Na+-counterion-neutralized polyacrylic acid (PAA) brushes are performed for varying degrees of ionization (and thereby varying charge density) and varying grafting density. Variation in the PE charge density (or degree of ionization) and grafting density leads to massive changes of the properties of the PE molecules (quantified by the changes in the height and the mobility of the PE brushes) as well as the local arrangement and distribution of the brush-supported counterions and water molecules within the brushes. The effect on the counterions is manifested by the corresponding variation of the counterion mobility, counterion concentration, extent of counterion binding to the charged site of the PE brushes, water-in-salt-like structure formation, and counterion-water–oxygen radial distribution function within the PE brushes. On the other hand, the effect on water molecules is manifested by the corresponding variation of water–oxygen–water–oxygen RDF, local water density, water–water and water–PE functional group hydrogen bond networks, static dielectric constant of water molecules, orientational tetrahedral order parameter, and water mobility. Enforcing such varying degree of ionization of weak polyelectrolytes is possible by changing the pH of the surrounding medium. Thus, our results provide insights into the changes in microstructure (at the atomistic level) of weak polyionic brushes at varying pH. Finally, we anticipate that this knowledge will prove to be vital for the efficient design of several nano-scale systems employing PE brushes such as nanomechanical gates, current rectifiers, etc.

42 ENGINEERING↗
DOE OSTI2015

A wrench in the works of human acetylcholinesterase: Soman induced conformational changes revealed by molecular dynamics simulations

Irreversible inactivation of human acetylcholinesterase (hAChE) by organophosphorous pesticides (OPs) and chemical weapon agents (CWA) has severe morbidity and mortality consequences. We present data from quantum mechanics/molecular mechanics (QM/MM) and 80 classical molecular dynamics (MD) simulations of the apo and soman-adducted forms of hAChE to investigate the effects on the dynamics and protein structure when the catalytic Serine 203 is phosphonylated. We find that the soman phosphonylation of the active site Ser203 follows a water assisted addition-elimination mechanism with the elimination of the fluoride ion being the highest energy barrier at 6.5 kcal/mole. We observe soman-dependent changes in backbone and sidechain motions compared to the apo form of the protein. These alterations restrict the soman-adducted hAChE to a structural state that is primed for the soman adduct to be cleaved and removed from the active site. The altered motions and resulting structures provide alternative pathways into and out of the hAChE active site. In the soman-adducted protein both side and back door pathways are viable for soman adduct access. Correlation analysis of the apo and soman adducted MD trajectories shows that the correlation of gorge entrance and back door motion is disrupted when hAChE is adducted. This supports the hypothesis that substrate and product can use two different pathways as entry and exit sites in the apo form of the protein. These alternative pathways have important implications for the rational design of medical countermeasures.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2018

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗
DOE OSTI2020

Application of First-Principles-Based Artificial Neural Network Potentials to Multiscale-Shock Dynamics Simulations on Solid Materials

The use of artificial neural network (ANN) potentials trained with first-principles calculations has emerged as a promising approach for molecular dynamics (MD) simulations encompassing large space and time scales while retaining first-principles accuracy. To date, however, the application of ANN-MD has been limited to near-equilibrium processes. Here we combine first-principles-trained ANN-MD with multiscale shock theory (MSST) to successfully describe far-from-equilibrium shock phenomena. Overall, our ANN-MSST-MD approach describes shock-wave propagation in solids with first-principles accuracy but a 5000 times shorter computing time. Accordingly, ANN-MD-MSST was able to resolve fine, long-time elastic deformation at low shock speed, which was impossible with first-principles MD because of the high computational cost. This work thus lays a foundation of ANN-MD simulation to study a wide range of far-from-equilibrium processes.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Time Step Considerations when Simulating Dynamic Behavior of High Performance Homes

Building energy simulations, especially those concerning pre-cooling strategies and cooling/heating peak demand management, require careful analysis and detailed understanding of building characteristics. Accurate modeling of the building thermal response and material properties for thermally massive walls or advanced materials like phase change materials (PCMs) are critically important.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Low-mass molecular dynamics simulation: A simple and generic technique to enhance configurational sampling

Highlights: • Reducing atomic masses by 10-fold vastly improves sampling in MD simulations. • CLN025 folded in 4 of 10 × 0.5-μs MD simulations when masses were reduced by 10-fold. • CLN025 folded as early as 96.2 ns in 1 of the 4 simulations that captured folding. • CLN025 did not fold in 10 × 0.5-μs MD simulations when standard masses were used. • Low-mass MD simulation is a simple and generic sampling enhancement technique. - Abstract: CLN025 is one of the smallest fast-folding proteins. Until now it has not been reported that CLN025 can autonomously fold to its native conformation in a classical, all-atom, and isothermal–isobaric molecular dynamics (MD) simulation. This article reports the autonomous and repeated folding of CLN025 from a fully extended backbone conformation to its native conformation in explicit solvent in multiple 500-ns MD simulations at 277 K and 1 atm with the first folding event occurring as early as 66.1 ns. These simulations were accomplished by using AMBER forcefield derivatives with atomic masses reduced by 10-fold on Apple Mac Pros. By contrast, no folding event was observed when the simulations were repeated using the original AMBER forcefields of FF12SB and FF14SB. The results demonstrate that low-mass MD simulation is a simple and generic technique to enhance configurational sampling. This technique may propel autonomous folding of a wide range of miniature proteins in classical, all-atom, and isothermal–isobaric MD simulations performed on commodity computers—an important step forward in quantitative biology.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Structure, Dynamics, and Hydrogen Transport in Amorphous Polymers: An Analysis of the Interplay between Free Volume Element Distribution and Local Segmental Dynamics from Molecular Dynamics Simulations

Polymers are attractive membrane materials owing to their mechanical robustness and relatively inexpensive fabrication. An important indicator of membrane performance are free volume elements (FVE): microporous void spaces created by the inefficient packing of bulky groups along the polymer chain. FVEs tend to degrade over time, as polymer chains reorganize irreversibly. While it is widely accepted that polymer flexibility has an impact on membrane transport properties, the molecular nature of this impact is still not well understood. By the establishment of a correlation between local chain dynamics and the distribution of free volume elements (FVEs), penetrant transport can be regulated more efficiently in amorphous polymer membranes. In this work, we implement all-atom molecular dynamics (MD) simulations to explore the relationship between chain dynamics and free volume in three polymers with different levels of backbone flexibility: polymethylpentene (PMP), polystyrene (PS), and HAB-6FDA thermally rearranged polymer (TRP). We construct these polymers at different temperatures and examine how temperature impacts the FVE distribution and segmental mobility. Our analysis shows that chain segments near FVEs have higher mobility compared with the atoms in the bulk; the extent of this difference increases with chain flexibility. Increasing the chain flexibility by increasing the temperature results in a broader FVE distribution. Rigid polymers such as TRP show the most robust FVE distribution and are not significantly affected by the temperature change. To capture penetrant diffusion through the polymer matrix, hydrogen is inserted and the diffusion is measured at different temperatures; hydrogen mobility is influenced by the FVE structure and overall mobility of polymer chains. Lastly, at low temperatures, hydrogen mobility is influenced by void distribution, while at high temperatures, polymer dynamics dictate hydrogen transport.

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

Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine

The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.

Artificial neural networks↗
DOE OSTI2015

Effect of Tapering on Morphology and Interfacial Behavior of Diblock Copolymers from Molecular Dynamics Simulations

Tapered diblock copolymers are AB diblock copolymers modified by adding a gradient region between the blocks in which composition varies smoothly from one species to the other. This gives additional control parameters, the length of the taper and its direction, to control the microphase separation behavior. Using tapers can also increase the accessibility of the bicontinuous gyroid phase, potentially of interest in transport applications. Recently, the phase diagram of these systems was predicted using self-consistent field theory (SCFT). Here, we perform coarse-grained molecular dynamics (MD) simulations to investigate both the structural and dynamic properties of tapered diblock copolymers. The MD results are consistent with the SCFT phase diagram, and the density profiles are also similar as a function of taper length and direction. We additionally compute mean-squared displacements and end-to-end relaxation times of lamellar systems and observe individual polymer conformations. Increasing taper length effectively lowers segregation strength and leads to a smaller domain spacing; the changes in dynamics are also discussed. The short (30%) inverse taper has a significantly shorter domain spacing than the diblock and qualitatively different polymer conformations (chains can fold back and forth across the interface).

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

Operator learning for energy-efficient building ventilation control with computational fluid dynamics simulation of a real-world classroom

Energy-efficient ventilation control plays an important role in reducing building energy consumption while ensuring occupant health and comfort. While Computational Fluid Dynamics (CFD) simulations provide detailed and physically accurate representations of indoor airflow, their high computational cost limits their use in real-time building control. In this work, we present a neural operator learning framework that combines the physical accuracy of CFD with the computational efficiency of machine learning to enable building ventilation control with the high-fidelity fluid dynamics models. Our method jointly optimizes the airflow supply rates and vent angles to reduce energy use and adhere to air quality constraints. We train an ensemble of neural operator transformer models to learn the mapping from building control actions to airflow fields using high-resolution CFD data. This learned neural operator is then embedded in an optimization-based control framework for building ventilation control. Experimental results show that our approach achieves significant energy savings compared to maximum airflow rate control, rule-based control, as well as data-driven control methods using spatially averaged CO 2 prediction and deep learning–based reduced-order models, while consistently maintaining safe indoor air quality. These results highlight the practicality and scalability of our method in maintaining energy efficiency and indoor air quality in real-world buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2018

CuCl Complexation in the Vapor Phase: Insights from Ab Initio Molecular Dynamics Simulations

We investigated the hydration of the CuCl 0 complex in HCl-bearing water vapor at 350°C and a vapor-like fluid density between 0.02 and 0.09 g/cm 3 using ab initio molecular dynamics (MD) simulations. The simulations reveal that one water molecule is strongly bonded to Cu(I) (first coordination shell), forming a linear [H 2 O-Cu-Cl] 0 moiety. The second hydration shell is highly dynamic in nature, and individual configurations have short life-spans in such low-density vapors, resulting in large fluctuations in instantaneous hydration numbers over a timescale of picoseconds. The average hydration number in the second shell ( m ) increased from ~0.5 to ~3.5 and the calculated number of hydrogen bonds per water molecule increased from 0.09 to 0.25 when fluid density (which is correlated to water activity) increased from 0.02 to 0.09 g/cm 3 ( f H 2 O 1.72 to 2.05). These changes of hydration number are qualitatively consistent with previous solubility studies under similar conditions, although the absolute hydration numbers from MD were much lower than the values inferred by correlating experimental Cu fugacity with water fugacity. This could be due to the uncertainties in the MD simulations and uncertainty in the estimation of the fugacity coefficients for these highly nonideal “vapors” in the experiments. Finally, our study provides the first theoretical confirmation that beyond-first-shell hydrated metal complexes play an important role in metal transport in low-density hydrothermal fluids, even if it is highly disordered and dynamic in nature.

58 GEOSCIENCES↗
NASA NTRS1986

Development and application of dynamic simulations of a subsonic wind tunnel

Efforts are currently underway at NASA Lewis to improve and expand ground test facilities and to develop supporting technologies to meet anticipated aeropropulsion research needs. Many of these efforts have been focused on a proposed rehabilitation of the Altitude Wind Tunnel (AWT). In order to insure a technically sound design, an AWT modeling program (both analytical and physical) was initiated to provide input to the AWT final design process. This paper describes the approach taken to develop analytical, dynamic computer simulations of the AWT, and the use of these simulations as test-beds for: (1) predicting the dynamic response characteristics of the AWT, and (2) evaluating proposed AWT control concepts. Plans for developing a portable, real-time simulator for the AWT facility are also described.

Szuch, J. R.↗
NASA NTRS1986

Development and application of dynamic simulations of a subsonic wind tunnel

Efforts are currently underway at NASA Lewis to improve and expand ground test facilities and to develop supporting technologies to meet anticipated aeropropulsion research needs. Many of these efforts have been focused on a proposed rehabilitation of the Altitude Wind Tunnel (AWT). In order to insure a technically sound design, an AWT modeling program (both analytical and physical) was initiated to provide input to the AWT final design process. This paper describes the approach taken to develop analytical, dynamic computer simulations of the AWT, and the use of these simulations as test-beds for (1) predicting the dynamic response characteristics of the AWT and (2) evaluating proposed AWT control concepts. Plans for developing a portable, real-time simulator for the AWT facility are also described.

Szuch, J. R.↗
DOE OSTI2014

Molecular dynamics simulation of diffusion coefficients and structural properties of some alkylbenzenes in supercritical carbon dioxide at infinite dilution

The binary infinite dilute diffusion coefficients, D{sub 12}{sup ∞}, of some alkylbenzenes (Ph-C{sub n}, from Ph-H to Ph-C{sub 12}) from 313 K to 333 K at 15 MPa in supercritical carbon dioxide (scCO{sub 2}) have been studied by molecular dynamics (MD) simulation. The MD values agree well with the experimental ones, which indicate MD simulation technique is a powerful way to predict and obtain diffusion coefficients of solutes in supercritical fluids. Besides, the local structures of Ph-C{sub n}/CO{sub 2} fluids are further investigated by calculating radial distribution functions and coordination numbers. It qualitatively convinces that the first solvation shell of Ph-C{sub n} in scCO{sub 2} is significantly influenced by the structure of Ph-C{sub n} solute. Meanwhile, the mean end-to-end distance, the mean radius of gyration and dihedral angle distribution are calculated to gain an insight into the structural properties of Ph-C{sub n} in scCO{sub 2}. The abnormal trends of radial distribution functions and coordination numbers can be reasonably explained in term of molecular flexibility. Moreover, the computed results of dihedral angle clarify that flexibility of long-chain Ph-C{sub n} is the result of internal rotation of C-C single bond (σ{sub c-c}) in alkyl chain. It is interesting that compared with n-alkane, because of the existence of benzene ring, the flexibility of alkyl chain in Ph-C{sub n} with same carbon atom number is significantly reduced, as a result, the carbon chain dependence of diffusion behaviors for long-chain n-alkane (n ≥ 5) and long-chain Ph-C{sub n} (n ≥ 4) in scCO{sub 2} are different.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗