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

Molecular Dynamics Simulation and Analysis of Interfacial Water at Selected Sulfide Mineral Surfaces under Anaerobic Conditions

In this paper, we report on a molecular dynamics simulation (MDS) study of the behavior of interfacial water at selected sulfide mineral surfaces under anaerobic conditions. The study revealed the interfacial water structure and wetting characteristics of the pyrite (100) surface, galena (100) surface, chalcopyrite (012) surface, sphalerite (110) surface, and molybdenite surfaces (i.e., the face, armchair-edge, and zigzag-edge surfaces), including simulated contact angles, relative number density profiles, water dipole orientations, hydrogen-bonding, and residence times. For force fields of the metal and sulfur atoms in selected sulfide minerals used in the MDS, we used the universal force field (UFF) and another set of force fields optimized by quantum chemical calculations for interactions with interfacial water molecules at selected sulfide mineral surfaces. Simulation results for the structural and dynamic properties of interfacial water molecules indicate the natural hydrophobic character for the selected sulfide mineral surfaces under anaerobic conditions as well as the relatively weak hydrophobicity for the sphalerite (110) surface and two molybdenite edge surfaces. Part of the financial support for this study was provided by the U.S. Department of Energy (DOE) under Basic Science Grant No. DE-FG-03-93ER14315. The Division of Chemical Sciences, Geosciences, and Biosciences, Office of Basic Energy Sciences (BES), of the DOE, funded work performed by Liem X. Dang. Battelle operates Pacific Northwest National Laboratory for DOE. The calculations were carried out using computer resources provided by BES. The authors are grateful to Professor Tsun-Mei Chang for valuable discussions.

Jin, Jiaqi↗
DOE OSTI2014

PuReMD-GPU: A reactive molecular dynamics simulation package for GPUs

We present an efficient and highly accurate GP-GPU implementation of our community code, PuReMD, for reactive molecular dynamics simulations using the ReaxFF force field. PuReMD and its incorporation into LAMMPS (Reax/C) is used by a large number of research groups worldwide for simulating diverse systems ranging from biomembranes to explosives (RDX) at atomistic level of detail. The sub-femtosecond time-steps associated with ReaxFF strongly motivate significant improvements to per-timestep simulation time through effective use of GPUs. This paper presents, in detail, the design and implementation of PuReMD-GPU, which enables ReaxFF simulations on GPUs, as well as various performance optimization techniques we developed to obtain high performance on state-of-the-art hardware. Comprehensive experiments on model systems (bulk water and amorphous silica) are presented to quantify the performance improvements achieved by PuReMD-GPU and to verify its accuracy. In particular, our experiments show up to 16× improvement in runtime compared to our highly optimized CPU-only single-core ReaxFF implementation. PuReMD-GPU is a unique production code, and is currently available on request from the authors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The structure of PX{sub 3} (X = Cl, Br, I) molecular liquids from X-ray diffraction, molecular dynamics simulations, and reverse Monte Carlo modeling

Synchrotron X-ray diffraction measurements have been conducted on liquid phosphorus trichloride, tribromide, and triiodide. Molecular Dynamics simulations for these molecular liquids were performed with a dual purpose: (1) to establish whether existing intermolecular potential functions can provide a picture that is consistent with diffraction data and (2) to generate reliable starting configurations for subsequent Reverse Monte Carlo modelling. Structural models (i.e., sets of coordinates of thousands of atoms) that were fully consistent with experimental diffraction information, within errors, have been prepared by means of the Reverse Monte Carlo method. Comparison with reference systems, generated by hard sphere-like Monte Carlo simulations, was also carried out to demonstrate the extent to which simple space filling effects determine the structure of the liquids (and thus, also estimating the information content of measured data). Total scattering structure factors, partial radial distribution functions and orientational correlations as a function of distances between the molecular centres have been calculated from the models. In general, more or less antiparallel arrangements of the primary molecular axes that are found to be the most favourable orientation of two neighbouring molecules. In liquid PBr{sub 3} electrostatic interactions seem to play a more important role in determining intermolecular correlations than in the other two liquids; molecular arrangements in both PCl{sub 3} and PI{sub 3} are largely driven by steric effects.

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

Hydration and reaction mechanisms on sodium silicate glass surfaces from molecular dynamics simulations with reactive force fields

Recent development of reactive force fields have enabled molecular dynamics simulations of interactions between silicate glasses and water at the atomistic scale. While multicomponent silicate glasses encompass a wide variety of compositions and properties, one common structural feature in these glasses is the combination of the network structure that is made up of silica tetrahedra linked through corner sharing interspersed with network modifiers like alkali and alkaline-earth ions that break up the Si–O–Si linkages by forming nonbridging oxygen. In reactions with water, ion exchange between alkali ions in the glass and proton or hydronium in the solution, as well as hydrolysis reaction of the Si–O–Si linkages and subsequent silanol formation, is observed and well documented. Furthermore, we have used a set of recently developed reactive force field to investigate the reactions between water and the surfaces of silica and sodium silicate glasses of different compositions for reactions up to 8 nanoseconds. Our results indicate sodium leaching into water and diffusion of water molecules up to 25 Å into the glass surface. We examined the structural and compositional changes inside the glass and around the diffused ions and use these to explain the rates of silanol formation at the surface. We also observed proton transport in the glass which has an indirect influence on the silanol formation rates. While the surface of the glass was rough to start with, it undergoes further modification into a hydrated gel-like structure in the glass for up to 5 Å in the higher alkali containing glasses. It was found that the leached sodium ions remain close to the interface and that fragments of silicate network from the surface is capable of dislodging from the bulk glass and enter the aqueous solution. These simulations thus provide insights into the formation and structure of an alteration layers commonly observed in multicomponent silicate glasses corroded in aqueous solutions.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Data Systems Dynamic Simulation - A total system for data system design assessments and trade studies

Data Systems Dynamic Simulation is a simulation system designed to reduce cost and time and increase the confidence and comprehensiveness of Data Systems Simulation. It is designed to simulate large data processing and communications systems from end-to-end or by subsystem. Those features relevant to system timing, control, sizing, personnel support activities, cost and external influences are modeled. Emphasis is placed on ease of use, comprehensive system performance measures, and extensive post simulation analysis capability. The system has been used to support trade studies of the NASA data system needs in the 1985 to 1990 time frame.

Hooper, J. W.↗
DOE OSTI2015

Molecular dynamical simulations of melting behaviors of metal clusters

The melting behaviors of metal clusters are studied in a wide range by molecular dynamics simulations. The calculated results show that there are fluctuations in the heat capacity curves of some metal clusters due to the strong structural competition; For the 13-, 55- and 147-atom clusters, variations of the melting points with atomic number are almost the same; It is found that for different metal clusters the dynamical stabilities of the octahedral structures can be inferred in general by a criterion proposed earlier by F. Baletto et al. [J. Chem. Phys. 116 3856 (2002)] for the statically stable structures.

36 MATERIALS SCIENCE↗
NASA NTRS1974

The effect of spurious angular accelerations on tracking in dynamic simulation

A laboratory study was conducted to investigate the effect of spurious simulator yaw motions on a pilot's control performance. A second objective was to compare the efficiency of static and dynamic simulator tracking in previously unexamined vehicle dynamics. Twelve airline pilots served as subjects in a moving-base flight simulator under congruent-motion, spurious-motion, and no-motion conditions. The results indicated a significant increase in the amount of error with increasing levels of spurious motion during the initially administered series of trials. The influence of spurious motion, however, was absent in a second series of trials. The data suggest that the pilots learned to compensate in their performance for the spurious inputs. It was also found that congruent visual and rotational cueing produced superior performance to that of tracking with visual information alone.

Beck, L. J.↗
DOE OSTI2023

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Modern HPC systems are increasingly relying on greater core counts and wider vector registers. Thus, applications need to be adapted to fully utilize these hardware capabilities. One class of applications that can benefit from this increase in parallelism are molecular dynamics simulations. In this paper, we describe our efforts at modernizing the ESPResSo++ simulation package for molecular dynamics by restructuring its particle data layout for efficient memory accesses and applying vectorization techniques to benefit the calculation of short-range non-bonded forces, which results in an overall three times speedup and serves as a baseline for further optimizations. We also implement fine-grained parallelism for multi-core CPUs through HPX, a C++ runtime system which uses lightweight threads and an asynchronous many-task approach to maximize concurrency. Our goal is to evaluate the performance of an HPX-based approach compared to the bulk-synchronous MPI-based implementation. This requires the introduction of an additional layer to the domain decomposition scheme that defines the task granularity. On spatially inhomogeneous systems, which impose a corresponding load-imbalance in traditional MPI-based approaches, we demonstrate that by choosing an optimal task size, the efficient work-stealing mechanisms of HPX can overcome the overhead of communication resulting in an overall 1.4 times speedup compared to the baseline MPI version.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Calibration and validation of a spar-type floating offshore wind turbine model using the FAST dynamic simulation tool

In this study, high-quality computer simulations are required when designing floating wind turbines because of the complex dynamic responses that are inherent with a high number of degrees of freedom and variable metocean conditions. In 2007, the FAST wind turbine simulation tool, developed and maintained by the U.S. Department of Energy's (DOE's) National Renewable Energy Laboratory (NREL), was expanded to include capabilities that are suitable for modeling floating offshore wind turbines. In an effort to validate FAST and other offshore wind energy modeling tools, DOE funded the DeepCwind project that tested three prototype floating wind turbines at 1/50 th scale in a wave basin, including a semisubmersible, a tension-leg platform, and a spar buoy. This paper describes the use of the results of the spar wave basin tests to calibrate and validate the FAST offshore floating simulation tool, and presents some initial results of simulated dynamic responses of the spar to several combinations of wind and sea states. Wave basin tests with the spar attached to a scale model of the NREL 5-megawatt reference wind turbine were performed at the Maritime Research Institute Netherlands under the DeepCwind project. This project included free-decay tests, tests with steady or turbulent wind and still water (both periodic and irregular waves with no wind), and combined wind/wave tests. The resulting data from the 1/50th model was scaled using Froude scaling to full size and used to calibrate and validate a full-size simulated model in FAST. Results of the model calibration and validation include successes, subtleties, and limitations of both wave basin testing and FAST modeling capabilities.

17 WIND ENERGY↗
DOE OSTI2020

Dynamic Simulation Technoeconomic Model for Power Generation

Sandia National Laboratories has built and successfully tested a dynamic simulation technoeconomic model of the Palo Verde Generating Station that is now being updated to help other US power plants improve operations. Palo Verde, located west of Phoenix, Arizona, is the largest electricity generator in the US at 4 GW. Palo Verde uses — 60 million gallons per day of treated wastewater from Phoenix to cool reactors, and disposes of blowdown in evaporation ponds. The model built for Palo Verde numerically evaluates the economic impact of changing, for example, alternative cooling technologies, water usage and treatment, and influent water chemistry, and is based on detailed accounting of mass, energy, and cash flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2002

The Stress-strain Behavior of Polymer-Nanotube Composites from Molecular Dynamics Simulations

Stress-strain curves of polymer-carbon nanotube composites are derived from molecular dynamics simulations of a single-walled carbon nanotube embedded in polyethylene. A comparison is made between the response to mechanical loading of a composite with a long, continuous nanotube (replicated via periodic boundary conditions) and the response of a composite with a short, discontinuous nanotube. Both composites are mechanically loaded in the direction of and transverse to the NT axis. The long-nanotube composite shows an increase in the stiffness relative to the polymer and behaves anisotropically under the different loading conditions. The short-nanotube composite shows no enhancement relative to the polymer, most probably because of its low aspect ratio. The stress-strain curves are compared with rule-of-mixtures predictions.

Frankland, S. J. V.↗
DOE OSTI2014

Molecular Dynamics Simulations of Tri-n-butyl-phosphate/n-Dodecane Mixture: Thermophysical Properties and Molecular Structure

Molecular dynamics simulations of tri-n-butyl-phosphate (TBP)/n-dodecane mixture in the liquid phase have been carried out using two recently developed TBP force field models (J. Phys. Chem. B 2012, 116, 305) in combination with the all-atom optimized potentials for liquid simulations (OPLS-AA) force field model for n-dodecane. Specifically, the electric dipole moment of TBP, mass density of the mixture, and the excess volume of mixing were computed with TBP mole fraction ranging from 0 to 1. It is found that the aforementioned force field models accurately predict the mass density of the mixture in the entire mole fraction range. Commensurate with experimental measurements, the electric dipole moment of the TBP was found to slightly increase with the mole fraction of TBP in the mixture. Also, in accord with experimental data, the excess volume of mixing is positive in the entire mole fraction range, peaking at TBP mole fraction range 0.3 0.5. Finally, a close examination of the spatial pair correlation functions between TBP molecules, and between TBP and n-dodecane molecules, revealed formation of TBP dimers through self-association at close distance, a phenomenon with ample experimental evidence.

de Almeida, Valmor F↗
NASA NTRS1993

Real-time dynamic simulation of the Cassini spacecraft using DARTS. Part 2: Parallel/vectorized real-time implementation

Part 1 of this paper presented the requirements for the real-time simulation of Cassini spacecraft along with some discussion of the DARTS algorithm. Here, in Part 2 we discuss the development and implementation of parallel/vectorized DARTS algorithm and architecture for real-time simulation. Development of the fast algorithms and architecture for real-time hardware-in-the-loop simulation of spacecraft dynamics is motivated by the fact that it represents a hard real-time problem, in the sense that the correctness of the simulation depends on both the numerical accuracy and the exact timing of the computation. For a given model fidelity, the computation should be computed within a predefined time period. Further reduction in computation time allows increasing the fidelity of the model (i.e., inclusion of more flexible modes) and the integration routine.

Fijany, A.↗
DOE OSTI2022

Molecular Dynamics Simulation and Cryo-Electron Microscopy Investigation of AOT Surfactant Structure at the Hydrated Mica Surface

Structural properties of the anionic surfactant dioctyl sodium sulfosuccinate (AOT or Aerosol-OT) adsorbed on the mica surface were investigated by molecular dynamics simulation, including the effect of surface loading in the presence of monovalent and divalent cations. The simulations confirmed recent neutron reflectivity experiments that revealed the binding of anionic surfactant to the negatively charged surface via adsorbed cations. At low loading, cylindrical micelles formed on the surface, with sulfate head groups bound to the surface by water molecules or adsorbed cations. Cation bridging was observed in the presence of weakly hydrating monovalent cations, while sulfate groups interacted with strongly hydrating divalent cations through water bridges. The adsorbed micelle structure was confirmed experimentally with cryogenic electronic microscopy, which revealed micelles approximately 2 nm in diameter at the basal surface. At higher AOT loading, the simulations reveal adsorbed bilayers with similar surface binding mechanisms. Adsorbed micelles were slightly thicker (2.2–3.0 nm) than the corresponding bilayers (2.0–2.4 nm). Upon heating the low loading systems from 300 K to 350 K, the adsorbed micelles transformed to a more planar configuration resembling bilayers. The driving force for this transition is an increase in the number of sulfate head groups interacting directly with adsorbed cations.

58 GEOSCIENCES↗
DOE OSTI2022

Propagating Parameter Uncertainty in Power System Nonlinear Dynamic Simulations Using a Koopman Operator-Based Surrogate Model

In this work, we propose a Koopman operator-based surrogate model for propagating parameter uncertainties in power system nonlinear dynamic simulations. First, we augment a priori known state-space model by reformulating parameters deemed uncertain as pseudo-state variables. Then, we apply the Koopman operator theory to the resulting state-space model and obtain a linear dynamical system model. This transformation allows us to analyze the evolution of the system dynamics through its Koopman eigenfunctions, eigenvalues, and modes. Of particular importance for this letter, the obtained linear dynamical system is a surrogate that enables the evaluation of parameter uncertainties by simply perturbing the initial conditions of the Koopman eigenfunctions associated with the pseudo-state variables. Simulations carried out on the New England test system reveal the excellent performance of the proposed method in terms of accuracy and computational efficiency.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Coarse-Grained Reactive Molecular Dynamics Simulations of Heterogeneities in Shocked Energetic Materials (LDRD Final Report)

This report summarizes the result of LDRD project 16-0161, titled "Coarse-Grained Re- active Molecular Dynamics Simulations of Heterogeneities in Shocked Energetic Materials." The purpose of the project was to develop a coarse-grained reactive molecular dynamics capability in LAMMPS enabling simulations of initiation in energetic materials comparable in accuracy to what is currently possible using large-scale reactive molecular dynamics, but with greatly reduced computational cost. The starting point for this work was the reactive dissipative particle dynamics (DPD) approach, which has been implemented as the new USER-DPD package in LAMMPS by researchers at Army Research Laboratory. Using modified versions of the examples provided with the new package, we examined the computational efficiency of the method, as well as its ability to model energy release in energetic materials. We observed that the Shardlow splitting method provides a great speed and accuracy advantage over conventional velocity Verlet time integration. We observed that the generic model of an energetic material provided with the USER-DPD package exhibited incomplete reaction under constant volume conditions. This was caused by quenching of the internal temperature of the molecules due to a rapid build-up of repulsive interactions between product gas components under constant volume conditions. Under constant pressure conditions, complete reaction was observed, as volume expansion prevented the buildup of strong repulsive interactions. Finally, a more realistic model calibrated to reproduce the equation of state of the RDX molecular crystal was examined. This model exhibited much less quenching of the internal temperature under constant volume conditions and reacted very rapidly under constant pressure conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Near-surface damage and mixing in Si-Cl 2 -Ar atomic layer etching processes: Insights from molecular dynamics simulations

Silicon-chlorine-argon (Si-Cl 2 -Ar) atomic layer etching (ALE) is simulated using classical molecular dynamics (MD). The simulations provide a detailed view into the near-surface region during ALE processing. Bombardment of Ar + ions creates a mixed amorphous region that significantly differs from the picture of ideal ALE. There is also a significant change in the Si etch yield and the etch product distribution as a function of Ar + ion fluence. The Si etch yield is the highest at the beginning of the bombardment step but eventually decays to the physical sputtering yield. Atomic Cl and silicon chlorides are major etch products at the start of an ion bombardment step, but quickly decay. Atomic Si yields remain relatively constant as a function of Ar + ion fluence. As a result, a new schematic of Si-Cl 2 -Ar ALE is presented in order to emphasize the complex behavior observed in MD simulations.

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

Rupture Dynamics Simulations for Shallow Crustal Earthquakes

In this study we perform simulations of rupture dynamics for shallow crustal earthquakes. The objective is to investigate the implication of multi-scale and depth-dependent stress drop variations on rupture dynamics, and creation of strong motion generation areas (SMGAs) and long period ground motion generation areas (LMGAs) for a M6.7crustal earthquake on a strike slip vertical fault. In our fault rupture simulation we used heterogenous stress drop models with correlated stochastic perturbations, following the work of Pitarka et al. (2009). The stress model with multiscale variability is designed to accommodate high-stress drop asperity areas and small-scale spatial stress variations. The stress drop and size of asperities are based on the IM2011 (Irikura and Miyake, 2011) recipe.

58 GEOSCIENCES↗