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

Nanoparticle generation and transport resulting from femtosecond laser ablation of ultrathin metal films: Time-resolved measurements and molecular dynamics simulations

The synthesis of metal nanoparticles by ultrafast laser ablation of nanometers-thick metal films has been studied experimentally and computationally. Near-threshold backside laser ablation of 2–20 nm-thick Pt films deposited on fused silica substrates was found to produce nanoparticles with size distributions that were bimodal for the thicker films, but collapsed into a single mode distribution for the thinnest film. Time-resolved imaging of blackbody emission from the Pt nanoparticles was used to reveal the nanoparticle propagation dynamics and estimate their temperatures. The observed nanoparticle plume was compact and highly forward-directed with a well-defined collective velocity that permitted multiple rebounds with substrates to be revealed. Large-scale molecular dynamics simulations were used to understand the evolution of compressive and tensile stresses in the thicker melted liquid films that lead to their breakup and ejection of two groups of nanoparticles with different velocity and size distributions. Ultrafast laser irradiation of ultrathin (few nm) metal films avoids the splitting of the film and appears to be a method well-suited to cleanly synthesize and deposit nanoparticles from semitransparent thin film targets in highly directed beams.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

System-size corrections for self-diffusion coefficients calculated from molecular dynamics simulations: The case of CO{sub 2}, n-alkanes, and poly(ethylene glycol) dimethyl ethers

Molecular dynamics simulations were carried out to study the self-diffusion coefficients of CO{sub 2}, methane, propane, n-hexane, n-hexadecane, and various poly(ethylene glycol) dimethyl ethers (glymes in short, CH{sub 3}O–(CH{sub 2}CH{sub 2}O){sub n}–CH{sub 3} with n = 1, 2, 3, and 4, labeled as G1, G2, G3, and G4, respectively) at different conditions. Various system sizes were examined. The widely used Yeh and Hummer [J. Phys. Chem. B 108, 15873 (2004)] correction for the prediction of diffusion coefficient at the thermodynamic limit was applied and shown to be accurate in all cases compared to extrapolated values at infinite system size. The magnitude of correction, in all cases examined, is significant, with the smallest systems examined giving for some cases a self-diffusion coefficient approximately 15% lower than the infinite system-size extrapolated value. The results suggest that finite size corrections to computed self-diffusivities must be used in order to obtain accurate results.

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

Understanding nucleic acid structural changes by comparing wide-angle x-ray scattering (WAXS) experiments to molecular dynamics simulations

Wide-angle x-ray scattering (WAXS) is emerging as a powerful tool for increasing the resolution of solution structure measurements of biomolecules. Compared to its better known complement, small angle x-ray scattering (SAXS), WAXS targets higher scattering angles and can enhance structural studies of molecules by accessing finer details of solution structures. Although the extension from SAXS to WAXS is easy to implement experimentally, the computational tools required to fully harness the power of WAXS are still under development. Currently, WAXS is employed to study structural changes and ligand binding in proteins; however, the methods are not as fully developed for nucleic acids. Here, we show how WAXS can qualitatively characterize nucleic acid structures as well as the small but significant structural changes driven by the addition of multivalent ions. We show the potential of WAXS to test all-atom molecular dynamics (MD) simulations and to provide insight into understanding how the trivalent ion cobalt(III) hexammine (CoHex) affects the structure of RNA and DNA helices. We find that MD simulations capture the RNA structural change that occurs due to addition of CoHex.

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

Molecular dynamics simulations of zinc oxide solubility: From bulk down to nanoparticles

The solubility of metal oxides is one of the key descriptors for the evaluation of their potential toxic effects, both in the bulk form and in nanoparticulated aggregates. Current work presents a new methodology for the in silico assessment of the solubility of metal oxides, which is demonstrated using a well-studied system, ZnO. The calculation of the solubility is based on statistical thermodynamics tools combined with Density Functional Tight Binding theory for the evaluation of the free energy exchange during the dissolution process. Models of small ZnO clusters are used for describing the final dissolved material, since the complete ionic dissolution of ZnO is hindered by the formation of O2- anions in solution, which are highly unstable. Results show very good agreement between the computed solubility values and experimental data for ZnO bulk, up to 0.5 mg·L-1 and equivalents of 50 g·L-1 for the free Zn2+ cation in solution. However, the reference model for solid nanoparticles formed by free space nanoparticles can only give a limited quantitative solubility evaluation for ZnO nanoparticles.

DFTB modelling↗
DOE OSTI2015

A coupled ordinates method for solution acceleration of rarefied gas dynamics simulations

Non-equilibrium rarefied flows are frequently encountered in a wide range of applications, including atmospheric re-entry vehicles, vacuum technology, and microscale devices. Rarefied flows at the microscale can be effectively modeled using the ellipsoidal statistical Bhatnagar–Gross–Krook (ESBGK) form of the Boltzmann kinetic equation. Numerical solutions of these equations are often based on the finite volume method (FVM) in physical space and the discrete ordinates method in velocity space. However, existing solvers use a sequential solution procedure wherein the velocity distribution functions are implicitly coupled in physical space, but are solved sequentially in velocity space. This leads to explicit coupling of the distribution function values in velocity space and slows down convergence in systems with low Knudsen numbers. Furthermore, this also makes it difficult to solve multiscale problems or problems in which there is a large range of Knudsen numbers. In this paper, we extend the coupled ordinates method (COMET), previously developed to study participating radiative heat transfer, to solve the ESBGK equations. In this method, at each cell in the physical domain, distribution function values for all velocity ordinates are solved simultaneously. This coupled solution is used as a relaxation sweep in a geometric multigrid method in the spatial domain. Enhancements to COMET to account for the non-linearity of the ESBGK equations, as well as the coupled implementation of boundary conditions, are presented. The methodology works well with arbitrary convex polyhedral meshes, and is shown to give significantly faster solutions than the conventional sequential solution procedure. Acceleration factors of 5–9 are obtained for low to moderate Knudsen numbers on single processor platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Molecular Dynamics Simulations of Water Structure and Diffusion in a 1 nm Diameter Silica Nanopore as a Function of Surface Charge and Alkali Metal Counterion Identity

Water confined in nanopores-particularly in pores narrower than 2 nm-displays distinct physicochemical properties that remain incompletely examined despite their importance in nanofluidics, molecular biology, geology, and materials sciences. In this study, we use molecular dynamics simulations to investigate the coordination structure and mobility of water and alkali metals (Li, Na, K, Cs) inside a 1 nm diameter cylindrical silica nanopore as a function of surface charge density, a model system particularly relevant to the alteration kinetics of silicate glasses and minerals in geologic formations. We find that the presence of a negative surface charge and adsorbed counterions within the pore strongly impacts water structure and dynamics. In particular, it significantly orients water O-H bonds toward the surface and slows water diffusion by almost 1 order of magnitude. Ion crowding in the charged nanopore enhances the tendency of counterions to coordinate closely with the silica surface, which moderates the impact of ions on water dynamics. Co-ions are strongly excluded from the nanopore at all surface charges, suggesting that the 1 nm diameter cylindrical silica nanopores likely exhibit nearly ideal semipermeable membrane transport properties.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Non-equilibrium dynamics in disordered materials: Ab initio molecular dynamics simulations

The dynamic properties of liquid B{sub 2}O{sub 3} under pressure and highly-charged bromophenol molecule are studied by using molecular dynamics (MD) simulations based on density functional theory (DFT). Diffusion properties of covalent liquids under high pressure are very interesting in the sense that they show unexpected pressure dependence. It is found from our simulation that the magnitude relation of diffusion coefficients for boron and oxygen in liquid B{sub 2}O{sub 3} shows the anomalous pressure dependence. The simulation clarified the microscopic origin of the anomalous diffusion properties. Our simulation also reveals the dissociation mechanism in the coulomb explosion of the highly-charged bromophenol molecule. When the charge state n is 6, hydrogen atom in the hydroxyl group dissociates at times shorter than 20 fs while all hydrogen atoms dissociate when n is 8. After the hydrogen dissociation, the carbon ring breaks at about 100 fs. There is also a difference on the mechanism of the ring breaking depending on charge states, in which the ring breaks with expanding (n = 6) or shrink (n = 8)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Visual Analytics for Deep Embeddings of Large Scale Molecular Dynamics Simulations

Molecular Dynamics (MD) simulation have been emerging as an excellent candidate for understanding complex atomic and molecular scale mechanism of bio-molecules that control essential bio-physical phenomenon in a living organism. But this MD technique produces large-size and long-timescale data that are inherently high-dimensional and occupies many terabytes of data. Processing this immense amount of data in a meaningful way is becoming increasingly difficult. Therefore, specific dimensionality reduction algorithm using deep learning technique has been employed here to embed the high-dimensional data in a lower-dimension latent space that still preserves the inherent molecular characteristics i.e. retains biologically meaningful information. Subsequently, the results of the embedding models are visualized for model evaluation and analysis of the extracted underlying features. However, most of the existing visualizations for embeddings have limitations in evaluating the embedding models and understanding the complex simulation data. We propose an interactive visual analytics system for embeddings of MD simulations to not only evaluate and explain an embedding model but also analyze various characteristics of the simulations. Our system enables exploration and discovery of meaningful and semantic embedding results and supports the understanding and evaluation of results by the quantitatively described features of the MD simulations (even without specific labels).

Chae, Junghoon↗
DOE OSTI2015

Look-ahead Dynamic Simulation

Look-ahead dynamic simulation software system incorporates the high performance parallel computing technologies, significantly reduces the solution time for each transient simulation case, and brings the dynamic simulation analysis into on-line applications to enable more transparency for better reliability and asset utilization. It takes the snapshot of the current power grid status, functions in parallel computing the system dynamic simulation, and outputs the transient response of the power system in real time.

Jin, Shuangshuang↗
DOE OSTI2021

Interfacial Reactivity and Speciation Emerging from Na-Montmorillonite Interactions with Water and Formic Acid at 200 °C: Insights from Reactive Molecular Dynamics Simulations, Infrared Spectroscopy, and X-ray Scattering Measurements

Reactive organic fluid - mineral interactions at elevated temperatures contribute to the evolution of planetary matter. One of the less studied but important transformations in this regard involves the reactions of formic acid with naturally occurring clays such as sodium montmorillonite. To advance a mechanistic understanding of these interactions, we use ReaxFF reactive molecular dynamics simulations in conjunction with infrared (IR) spectroscopy and X-ray scattering experiments to investigate the speciation behavior of water-formic acid mixtures on sodium montmorillonite interfaces at 473 K and 1 atm. Using a newly developed reactive forcefield, we show that the experimental IR spectra of unreacted and reacted mixture can be accurately reproduced by ReaxFF/MD. We further benchmark the simulation predictions of sodium carbonate and bicarbonate formation in the clay interlayers using Small and Wide-Angle X-ray Scattering measurements. Subsequently, leveraging the benchmarked forcefield, we interrogate the pathway of speciation reactions with emphasis on carbonate, formate, and hydroxide groups elucidating the energetics, transition states, intermediates, and preferred products. Further, we also delineate the differences in reactivities and catalytic effects of clay edges, facets, and interlayers owing to their local chemical environments, which have far reaching consequences in their speciation behavior. The experimental and simulation approaches described in this study and the transferable forcefields can be applied translationally to advance the science of clay-fluid interactions for several applications including subsurface fluid storage and recovery and clay-pollutant dynamics.

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

High-Performance Molecular Dynamics Simulation for Biological and Materials Sciences: Challenges of Performance Portability, Dallas, TX, November 2018.

Highly-optimized parallel molecular dynamics programs have allowed researchers to achieve ground-breaking results in biological and materials sciences. This type of performance has come at the expense of portability: a significant effort is required for performance optimization on each new architecture. Using a metric that emphasizes speedup, we assess key accelerating programming components of four different best-performing molecular dynamics programs - GROMACS, NAMD, LAMMPS and CP2K-- each having a particular scope of application, for contribution to performance and for portability. We use builds with and without these components, tested on HPC systems. We also analyze the code-bases to determine compliance with portability recommendations. We find that for all four programs, the contributions of the non-portable components to speed are essential to the programs' performances; without them we see a reduction in time-to-solution of a magnitude that is insufferable to domain scientists. This characterizes the performance efficiency that must be approached for good performance portability on a programmatic level, suggesting solutions to this difficult problem, which should come from developers, industry and funding institutions, and possibly new research in programming languages.

Computer Science↗