Search NASA⌕ Search

SEARCH · Search NASA

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations

We introduce an extension of the divide-and-conquer (DC) algorithmic paradigm called divide-conquer-recombine (DCR) to perform large quantum molecular dynamics (QMD) simulations on massively parallel supercomputers, in which interatomic forces are computed quantum mechanically in the framework of density functional theory (DFT). In DCR, the DC phase constructs globally informed, overlapping local-domain solutions, which in the recombine phase are synthesized into a global solution encompassing large spatiotemporal scales. For the DC phase, we design a lean divide-and-conquer (LDC) DFT algorithm, which significantly reduces the prefactor of the O(N) computational cost for N electrons by applying a density-adaptive boundary condition at the peripheries of the DC domains. Our globally scalable and locally efficient solver is based on a hybrid real-reciprocal space approach that combines: (1) a highly scalable real-space multigrid to represent the global charge density; and (2) a numerically efficient plane-wave basis for local electronic wave functions and charge density within each domain. Hybrid space-band decomposition is used to implement the LDC-DFT algorithm on parallel computers. A benchmark test on an IBM Blue Gene/Q computer exhibits an isogranular parallel efficiency of 0.984 on 786 432 cores for a 50.3 × 10{sup 6}-atom SiC system. As a test of production runs, LDC-DFT-based QMD simulation involving 16 661 atoms is performed on the Blue Gene/Q to study on-demand production of hydrogen gas from water using LiAl alloy particles. As an example of the recombine phase, LDC-DFT electronic structures are used as a basis set to describe global photoexcitation dynamics with nonadiabatic QMD (NAQMD) and kinetic Monte Carlo (KMC) methods. The NAQMD simulations are based on the linear response time-dependent density functional theory to describe electronic excited states and a surface-hopping approach to describe transitions between the excited states. A series of techniques are employed for efficiently calculating the long-range exact exchange correction and excited-state forces. The NAQMD trajectories are analyzed to extract the rates of various excitonic processes, which are then used in KMC simulation to study the dynamics of the global exciton flow network. This has allowed the study of large-scale photoexcitation dynamics in 6400-atom amorphous molecular solid, reaching the experimental time scales.

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

Crack propagation in silica from reactive classical molecular dynamics simulations

Mechanistic insight into the process of crack growth can be obtained through molecular dynamics (MD) simulations. In this investigation of fracture propagation, a slit crack was introduced into an atomistic amorphous silica model and mode I stress was applied through far–field loading until the crack propagates. Atomic displacements and forces and an Irving–Kirkwood method with a Lagrangian kernel estimator were used to calculate the J –integral of classical fracture mechanics around the crack tip. The resulting fracture toughness (K IC ), 0.76 ± 0.16 MPa√m, agrees with experimental values. In addition, the stress fields and dissipation energies around the slit crack indicate the development of an inelastic region ~30Å in diameter. This is one of the first reports of KIC values obtained from up–scaled atomic–level energies and stresses through the J –integral. In conclusion, the application of the ReaxFF classical MD force field in this study provides the basis for future research into crack growth in multicomponent oxides in a variety of environmental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗
DOE OSTI2014

Molecular-dynamics simulations of thin polyisoprene films confined between amorphous silica substrates

Constant temperature–constant pressure (NpT) molecular-dynamics computer simulations have been carried out for the united-atom model of a non-crosslinked (1,4) cis-polyisoprene (PI) melt confined between two amorphous, fully coordinated silica surfaces. The Lennard-Jones 12-6 potential was implemented to describe the polymer–silica interactions. The thickness H of the produced PI–silica film has been varied in a wide range, 1 < H/R{sub g} < 8, where R{sub g} is the individual PI chain radius of gyration measured under the imposed confinement. After a thorough equilibration, the PI film stratified structure and polymer segmental dynamics have been studied. The chain structure in the middle of the films resembles that in a corresponding bulk, but the polymer-density profile shows a pronounced ordering of the polymer segments in the vicinity of silica surfaces; this ordering disappears toward the film middles. Tremendous slowing down of the polymer segmental dynamics has been observed in the film surface layers, with the segmental relaxation more than 150 times slower as compared to that in a PI bulk. This effect increases with decreasing the polymer-film thickness. The segmental relaxation in the PI film middles shows additional relaxation process which is absent in a PI bulk. Even though there are fast relaxation processes in the film middle, its overall relaxation is slower as compared to that in a bulk sample. The interpretation of the results in terms of polymer glassy bridges has been discussed.

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

Periodic boundary conditions for arbitrary deformations in molecular dynamics simulations

A generalization of the Lees-Edwards periodic boundary conditions (gLE-PBC) for molecular dynamics (MD) simulations is developed to allow for arbitrary deformations to be applied to the domain. The gLE-PBC domain remains a rectangular cuboid regardless of the applied deformation in contrast with the Lagrangian-rhomboid periodic boundary conditions (LR-PBC) where the domain deforms according to the applied deformation. Furthermore, the kinematics of gLE-PBC are validated against pure shear. The gLE-PBC method for interacting systems is then validated against the LR-PBC method and analytical solutions for a solid under isotropic compression, one-dimensional shearing, three-dimensional extension and shearing and for a liquid under Couette flow. Bulk physical properties extracted from the gLE-PBC simulations agree well with values calculated from equilibrium MD simulations. Three dimensional shearing and a deformation with a full velocity gradient matrix are also simulated, showing the range of problems gLE-PBC can explore.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Accelerating discrete dislocation dynamics simulations with graph neural networks

Discrete dislocation dynamics (DDD) is a widely employed computational method to study plasticity at the mesoscale that connects the motion of dislocation lines to the macroscopic response of crystalline materials. However, the computational cost of DDD simulations remains a bottleneck that limits its range of applicability. Here, we introduce a new DDD-GNN framework in which the expensive time-integration of dislocation motion is entirely substituted by a graph neural network (GNN) model trained on DDD trajectories. As a first application, we demonstrate the feasibility and potential of our method on a simple yet relevant model of a dislocation line gliding through an array of obstacles. We show that the DDD-GNN model is stable and reproduces very well unseen ground-truth DDD simulation responses for a range of straining rates and obstacle densities, without the need to explicitly compute nodal forces or dislocation mobilities during time-integration. Our approach opens new promising avenues to accelerate DDD simulations and to incorporate more complex dislocation motion behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Chemical effects on subcritical fracture in silica from molecular dynamics simulations

Fracture toughness of silicates is reduced in aqueous environments due to water-silica interactions at the crack tip. In this paper, to investigate this effect, classical molecular dynamics (MD) simulations using the bond-order based reactive force field (ReaxFF) were used to simulate silica fracture. The chemical and mechanical aspects were separated by simulating fracture in: (a) vacuum with dynamic loading, (b) an aqueous environment with dynamic loading, and (c) an aqueous environment with static subcritical mechanical loading to track silica dissolution. The addition of water to silica fracture reduced the silica fracture toughness by ~25%, a trend consistent with experimentally reported results. Analysis of Si-O bonds in the process zone and calculations of dissipation energy associated with fracture indicated that water relaxes the entire process zone, and not just the surface. Additionally, the crack tip sharpens during fracture in water and an increased number of microscopic propagation events occur. This results in earlier fracture in systems with increasing mechanical loading in aqueous conditions, despite the lack of significant silica dissolution. Finally and therefore, the threshold for Si-O bond breakage has been lowered in the presence of water and the reduction in fracture toughness is due to structural and energetic changes in the silica, rather than specific dissolution events.

58 GEOSCIENCES↗
DOE OSTI2015

The interplay between dynamic heterogeneities and structure of bulk liquid water: A molecular dynamics simulation study

In order to study the interplay between dynamical heterogeneities and structural properties of bulk liquid water in the temperature range 130–350 K, thus including the supercooled regime, we use the explicit trend of the distribution functions of some molecular properties, namely, the rotational relaxation constants, the atomic mean-square displacements, the relaxation of the cross correlation functions between the linear and squared displacements of H and O atoms of each molecule, the tetrahedral order parameter q and, finally, the number of nearest neighbors (NNs) and of hydrogen bonds (HBs) per molecule. Two different potentials are considered: TIP4P-Ew and a model developed in this laboratory for the study of nanoconfined water. The results are similar for the dynamical properties, but are markedly different for the structural characteristics. In particular, for temperatures higher than that of the dynamic crossover between “fragile” (at higher temperatures) and “strong” (at lower temperatures) liquid behaviors detected around 207 K, the rotational relaxation of supercooled water appears to be remarkably homogeneous. However, the structural parameters (number of NNs and of HBs, as well as q) do not show homogeneous distributions, and these distributions are different for the two water models. Another dynamic crossover between “fragile” (at lower temperatures) and “strong” (at higher temperatures) liquid behaviors, corresponding to the one found experimentally at T{sup ∗} ∼ 315 ± 5 K, was spotted at T{sup ∗} ∼ 283 K and T{sup ∗} ∼ 276 K for the TIP4P-Ew and the model developed in this laboratory, respectively. It was detected from the trend of Arrhenius plots of dynamic quantities and from the onset of a further heterogeneity in the rotational relaxation. To our best knowledge, it is the first time that this dynamical crossover is detected in computer simulations of bulk water. On the basis of the simulation results, the possible mechanisms of the two crossovers at molecular level are discussed.

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

Structural and chemical orders in N i 64.5 Z r 35.5 metallic glass by molecular dynamics simulation

The atomic structure of Ni 64.5 Zr 35.5 metallic glass has been investigated by molecular dynamics (MD) simulations. The calculated structure factors from the MD glassy sample at room temperature agree well with the X-ray diffraction (XRD) and neutron diffraction (ND) experimental data. Using the pairwise cluster alignment and clique analysis methods, we show that there are three types dominant short-range order (SRO) motifs around Ni atoms in the glass sample of Ni 64.5 Zr 35.5 , i.e., Mixed- Icosahedron(ICO)-Cube, Twined-Cube and icosahedron-like clusters. Furthermore, chemical order and medium-range order (MRO) analysis show that the Mixed-ICOCube and Twined-Cube clusters exhibit the characteristics of the crystalline B2 phase. In conclusion, our simulation results suggest that the weak glass-forming ability (GFA) of Ni 64.5 Zr 35.5 can be attributed to the competition between the glass forming ICO SRO and the crystalline Mixed-ICO-Cube and Twined-Cube motifs.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Executive Summary: Special Section on Credible Computational Fluid Dynamics Simulations

This summary presents the motivation for the Special Section on the credibility of computational fluid dynamics (CFD) simulations, its objective, its background and context, its content, and its major conclusions. Verification and validation (V&V) are the processes for establishing the credibility of CFD simulations. Validation assesses whether correct things are performed and verification assesses whether they are performed correctly. Various aspects of V&V are discussed. Progress is made in verification of simulation models. Considerable effort is still needed for developing a systematic validation method that can assess the credibility of simulated reality.

Mehta, Unmeel B.↗
DOE OSTI2015

Melting behaviour of gold-platinum nanoalloy clusters by molecular dynamics simulations

The melting behavior of bimetallic gold-platinum nanoclusters is studied by applying Brownian-type isothermal molecular dynamics (MD) simulation, a program modified from the cubic coupling scheme (CCS). The process begins with the ground-state structures obtained from global minimum search algorithm and proceeds with the investigation of the effect of temperature on the thermal properties of gold-platinum nanoalloy clusters. N-body Gupta potential has been employed in order to account for the interactions between gold and platinum atoms. The ground states of the nanoalloy clusters, which are core-shell segregated, are heated until they become thermally segregated. The detailed melting mechanism of the nanoalloy clusters is studied via this approach to provide insight into the thermal stability of the nanoalloy clusters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Molecular Dynamics Simulation of Carbon Nanotube Based Gears

We used molecular dynamics to investigate the properties and design space of molecular gears fashioned from carbon nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. One gear was powered by forcing the atoms near the end of the buckytube to rotate, and a second gear was allowed.to rotate by keeping the atoms near the end of its buckytube on a cylinder. The meshing aromatic gear teeth transfer angular momentum from the powered gear to the driven gear. A number of gear and gear/shaft configurations were simulated. Cases in vacuum and with an inert atmosphere were examined. In an extension to molecular dynamics technology, some simulations used a thermostat on the atmosphere while the hydrocarbon gear's temperature was allowed to fluctuate. This models cooling the gears with an atmosphere. Results suggest that these gears can operate at up to 50-100 gigahertz in a vacuum or inert atmosphere at room temperature. The failure mode involves tooth slip, not bond breaking, so failed gears can be returned to operation by lowering temperature and/or rotation rate. Videos and atomic trajectory files in xyz format are presented.

Han, Jie↗
DOE OSTI2016

Counter-ion binding and mobility in the presence of hydrophobic polyions – combining molecular dynamics simulations and NMR

Counter-ion binding and mobility in aqueous solutions of partially hydrophobic ionene oligoions is studied here by a combination of all-atomic molecular dynamics (MD) simulations and NMR ({sup 19}F and {sup 81}Br nuclei) measurements. We present results for 12, 12–ionenes in the presence of different halide ions (F{sup −}, Cl{sup −}, Br{sup −} and I{sup −}), as well as their mixtures; the latter allowing us to probe counter-ion selectivity of these oligoions. We consolidate both structural and dynamic information, in particular simulated radial distribution functions and average residence times of counter-ions in the vicinity of ionenes and NMR data in the form of counter-ion chemical shift and self-diffusion coefficients. On one hand, previously reported enthalpy of dilution and mixing measurements show a reverse counter-ion sequence for 12, 12–ionenes with respect to their less hydrophobic 3, 3– and 6, 6– analogues. On the other hand, the current MD and NMR data, reflecting the counter-ion binding tendencies to the ionene chain, give evidence for the same ordering as that observed by MD for 3, 3–ionenes. This is not seen as a contradiction and can be rationalized on the basis of increasing chain hydrophobicity, which has different consequences for enthalpy and ion-binding. The latter is reflecting free energy changes and as such includes both enthalpic and entropic contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

The Molecular Structure of a Phosphatidylserine Bilayer Determined by Scattering and Molecular Dynamics Simulations

Phosphatidylserine (PS) lipids play essential roles in biological processes, including enzyme activation and apoptosis. We report on the molecular structure and atomic scale interactions of a fluid bilayer composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS). A scattering density profile model, aided by molecular dynamics (MD) simulations, was developed to jointly refine different contrast small-angle neutron and X-ray scattering data, which yielded a lipid area of 62.7 A2 at 25 C. MD simulations with POPS lipid area constrained at different values were also performed using all-atom and aliphatic united-atom models. The optimal simulated bilayer was obtained using a model-free comparison approach. Examination of the simulated bilayer, which agrees best with the experimental scattering data, reveals a preferential interaction between Na+ ions and the terminal serine and phosphate moieties. Long-range inter-lipid interactions were identified, primarily between the positively charged ammonium, and the negatively charged carboxylic and phosphate oxygens. The area compressibility modulus KA of the POPS bilayer was derived by quantifying lipid area as a function of surface tension from area-constrained MD simulations. It was found that POPS bilayers possess a much larger KA than that of neutral phosphatidylcholine lipid bilayers. We propose that the unique molecular features of POPS bilayers may play an important role in certain physiological functions.

biomembranes↗
DOE OSTI2013

Slow dynamics of a protein backbone in molecular dynamics simulation revealed by time-structure based independent component analysis

We recently proposed the method of time-structure based independent component analysis (tICA) to examine the slow dynamics involved in conformational fluctuations of a protein as estimated by molecular dynamics (MD) simulation [Y. Naritomi and S. Fuchigami, J. Chem. Phys. 134, 065101 (2011)]. Our previous study focused on domain motions of the protein and examined its dynamics by using rigid-body domain analysis and tICA. However, the protein changes its conformation not only through domain motions but also by various types of motions involving its backbone and side chains. Some of these motions might occur on a slow time scale: we hypothesize that if so, we could effectively detect and characterize them using tICA. In the present study, we investigated slow dynamics of the protein backbone using MD simulation and tICA. The selected target protein was lysine-, arginine-, ornithine-binding protein (LAO), which comprises two domains and undergoes large domain motions. MD simulation of LAO in explicit water was performed for 1 μs, and the obtained trajectory of C{sub α} atoms in the backbone was analyzed by tICA. This analysis successfully provided us with slow modes for LAO that represented either domain motions or local movements of the backbone. Further analysis elucidated the atomic details of the suggested local motions and confirmed that these motions truly occurred on the expected slow time scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗