Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic simulations”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Atomic structure evolution during solidification of liquid niobium from ab initio molecular dynamics simulations

Atomic structure transitions of liquid niobium during solidification, at different temperatures from 3200 to 1500 K, were studied by using ab initio molecular dynamics simulations. The local atomic structure variations with temperature are investigated by using the pair-correlation function, the structure factor, the bond-angle distribution function, the Honeycutt–Anderson index, Voronoi tessellation and the cluster alignment methods. Our results clearly show that, upon quenching, the icosahedral short-range order dominates in the stable liquid and supercooled liquid states before the system transforms to crystalline body-center cubic phase at a temperature of about 1830 K.

Condensed matter↗
DOE OSTI2023

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗
DOE OSTI2015

Molecular dynamics simulation of annealed ZnO surfaces

The effect of thermally annealing a slab of wurtzite ZnO, terminated by two surfaces, (0001) (which is oxygen-terminated) and (0001{sup ¯}) (which is Zn-terminated), is investigated via molecular dynamics simulation by using reactive force field (ReaxFF). We found that upon heating beyond a threshold temperature of ∼700 K, surface oxygen atoms begin to sublimate from the (0001) surface. The ratio of oxygen leaving the surface at a given temperature increases as the heating temperature increases. A range of phenomena occurring at the atomic level on the (0001) surface has also been explored, such as formation of oxygen dimers on the surface and evolution of partial charge distribution in the slab during the annealing process. It was found that the partial charge distribution as a function of the depth from the surface undergoes a qualitative change when the annealing temperature is above the threshold temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Three-dimensional ordering of cold ion beams in a storage ring: A molecular-dynamics simulation study

Three-dimensional (3D) ordering of a charged-particle beams circulating in a storage ring is systematically studied with a molecular-dynamics simulation code. An ion beam can exhibit a 3D ordered configuration at ultralow temperature as a result of powerful 3D laser cooling. Various unique characteristics of the ordered beams, different from those of crystalline beams, are revealed in detail, such as the single-particle motion in the transverse and longitudinal directions, and the dependence of the tune depression and the Coulomb coupling constant on the operating points.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Field Programmable Gate Arrays for Enhancing the Speed and Energy Efficiency of Quantum Dynamics Simulations

We present the first application of field programmable gate arrays (FPGAs) as new, customizable hardware architectures for carrying out fast and energy-efficient quantum dynamics simulations of large chemical/material systems. Instead of tailoring the software to fixed hardware, which is the typical case for writing quantum chemistry code for central processing units (CPUs) and graphics processing units (GPUs), FPGAs allow us to directly customize the underlying hardware (even at the level of specific electrical signals in the circuit) to give a truly optimized computational performance for quantum dynamics calculations. By offloading the most intensive and repetitive calculations onto an FPGA, we show that the computational performance of our real-time electron dynamics calculations can even exceed that of optimized commercial mathematical libraries running on high-performance GPUs. In addition to this impressive computational speedup, we show that FPGAs are immensely energy-efficient and consume 4 times less energy than modern GPU or CPU architectures. These energy savings are a practical and important metric for supercomputing centers (many of which exceed over $1 million in power costs alone), as exascale computing capabilities become more widespread and commonplace. Taken together, the implementation techniques and performance metrics of our study demonstrate that FPGAs could play a promising role in upcoming quantum chemistry and materials science applications, particularly for the acceleration and energy-efficient execution of quantum dynamics calculations.

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

Beam Dynamics simulations for ERDC project -- SRF linac for industrial use

Compact conductively cooled SRF industrial linacs can provide unique parameters of the electron beam for industrial applications. (up to 10MeV, 1MW). For ERDC project we designed normal conducting RF injector with thermal RF gridded gun integrated in first cell of multi-cell cavities. For design of the RF gun we used MICHELLE software to simulate and optimize parameters of the beam. Output file was converted to ASTRA format and most beam dynamic simulations in multi-cell normal conducting cavity and cryomodule were performed by using ASTRA software. For cross-checking we compare results of MICHELLE and AS-TRA in first few cells. At the end of injector beam reach ~250keV energy which allow to trap bunch in acceleration regime without losses in TESLA like 1.3 GHz cavity. Short solenoid at the end of injector allow to regulate transverse beam size in cryomodule to match beam to extraction system and also reduce charge losses in accelerator.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

CO 2 Capture Characteristics of Hyperbranched Poly(alkylene imine): A Molecular Dynamics Simulation Approach

This study explores the CO 2 capture characteristics of hyperbranched poly(ethylenimine) (HB-PEI) and poly- (propyleneimine) (HB-PPI) through molecular dynamics simulations using density functional theory-calibrated force fields. Key features such as density, free volume, glass transition temperature, CO 2 /H 2 O distribution, and molecular diffusion are systematically investigated to elucidate structure−function relationships under dry and hydrated conditions. HB-PEI demonstrates a slightly higher density and lower free volume compared to HB-PPI yet shows superior CO 2 capture due to the high amine concentration. Glass transition analysis indicates a higher thermal mobility in HBPEI, enhancing the CO 2 diffusivity. Pair correlation and coordination analyses confirm a stronger affinity of CO 2 with primary and secondary amines, particularly in hydrated environments where water competes with CO 2 for binding sites. Despite its more compact structure, HB-PEI outperformed HB-PPI in CO 2 and H 2 O transport, as confirmed by higher diffusion coefficients across all hydration levels. These findings highlight a critical balance among polymer architecture, amine accessibility, and hydration in designing next-generation solid amine sorbents for efficient direct air capture applications.

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

Cavity molecular dynamics simulations of liquid water under vibrational ultrastrong coupling

We simulate vibrational strong coupling (VSC) and vibrational ultrastrong coupling (V-USC) for liquid water with classical molecular dynamics simulations. When the cavity modes are resonantly coupled to the O–H stretch mode of liquid water, the infrared spectrum shows asymmetric Rabi splitting. The lower polariton (LP) may be suppressed or enhanced relative to the upper polariton (UP) depending on the frequency of the cavity mode. Furthermore, although the static properties and the translational diffusion of water are not changed under VSC or V-USC, we do find the modification of the orientational autocorrelation function of H 2 O molecules especially under V-USC, which could play a role in ground-state chemistry.

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

Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory

Iterative energy minimization with the aim of achieving self-consistency is a common feature of Born-Oppenheimer molecular dynamics (BOMD) and classical molecular dynamics with polarizable force fields. In the former, the electronic degrees of freedom are optimized, while the latter often involves an iterative determination of induced point dipoles. The computational effort of the self-consistency procedure can be reduced by re-using converged solutions from previous time steps. However, this must be done carefully, as not to break time-reversal symmetry, which negatively impacts energy conservation. Self-consistent schemes based on the extended Lagrangian formalism, where the initial guesses for the optimized quantities are treated as auxiliary degrees of freedom, constitute one elegant solution. We report on the performance of two integration schemes with the same underlying extended Lagrangian structure, which we both employ in two radically distinct regimes—in classical molecular dynamics simulations with the AMOEBA polarizable force field and in BOMD simulations with the Onetep linear-scaling density functional theory (LS-DFT) approach. Furthermore, both integration schemes are found to offer significant improvements over the standard (unpropagated) molecular dynamics formulation in both the classical and LS-DFT regimes.

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

Membrane Sculpting by F-BAR Domains Studied by Molecular Dynamics Simulations

Interplay between cellular membranes and their peripheral proteins drives many processes in eukaryotic cells. Proteins of the Bin/Amphiphysin/Rvs (BAR) domain family, in particular, play a role in cellular morphogenesis, for example curving planar membranes into tubular membranes. However, it is still unclear how F-BAR domain proteins act on membranes. Electron microscopy revealed that, in vitro, F-BAR proteins form regular lattices on cylindrically deformed membrane surfaces. Using all-atom and coarse-grained (CG) molecular dynamics simulations, we show that such lattices, indeed, induce tubes of observed radii. A 250 ns all-atom simulation reveals that F-BAR domain curves membranes via the so-called scaffolding mechanism. Plasticity of the F-BAR domain permits conformational change in response to membrane interaction, via partial unwinding of the domains 3-helix bundle structure. A CG simulation covering more than 350 µs provides a dynamic picture of membrane tubulation by lattices of F-BAR domains. A series of CG simulations identified the optimal lattice type for membrane sculpting, which matches closely the lattices seen through cryo-electron microscopy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Molecular dynamics simulations of a dicationic ionic liquid for CO 2 capture

A dicationic ionic liquid ([DBU-PEG][Tf 2 N] 2 ) was studied using classical molecular dynamics simulations to examine its structural and gas separation properties. The dication was designed in an attempt to improve CO 2 solubility by means of tuning the cation-anion interactions of the ionic liquid (IL). The computational model was compared to experimentally obtained density, viscosity, and powder X-ray diffraction spectra. The structure of the IL was further investigated with radial distribution functions and free volume analysis through cavity distributions. It was found that the shape and charge distribution of the dication enhances CO 2 interaction: the CO 2 molecule is hugged by the dication along the PEG linker and close to one of the cationic ends. Here, the geminal design of the dication allows for strong interaction with CO 2 , showing promise as a means of carbon capture.

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

A Scalable O(N) Algorithm for Large-Scale Parallel First-Principles Molecular Dynamics Simulations

Traditional algorithms for first-principles molecular dynamics (FPMD) simulations only gain a modest capability increase from current petascale computers, due to their O(N 3 ) complexity and their heavy use of global communications. To address this issue, we are developing a truly scalable O(N) complexity FPMD algorithm, based on density functional theory (DFT), which avoids global communications. The computational model uses a general nonorthogonal orbital formulation for the DFT energy functional, which requires knowledge of selected elements of the inverse of the associated overlap matrix. We present a scalable algorithm for approximately computing selected entries of the inverse of the overlap matrix, based on an approximate inverse technique, by inverting local blocks corresponding to principal submatrices of the global overlap matrix. The new FPMD algorithm exploits sparsity and uses nearest neighbor communication to provide a computational scheme capable of extreme scalability. Accuracy is controlled by the mesh spacing of the finite difference discretization, the size of the localization regions in which the electronic orbitals are confined, and a cutoff beyond which the entries of the overlap matrix can be omitted when computing selected entries of its inverse. We demonstrate the algorithm's excellent parallel scaling for up to O(100K) atoms on O(100K) processors, with a wall-clock time of O(1) minute per molecular dynamics time step.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics of plasma atomic layer etching: Molecular dynamics simulations and optical emission spectroscopy

Atomic layer etching is intrinsically dynamic as it involves sequential and repeated exposures of a surface to be etched with different species at different energies. The composition and structure of the near surface region change in both time and depth. Full understanding of this process requires resolving both temporal and spatial variations. In this work, we consider silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ). Molecular dynamics (MD) simulations are compared to experimental measurements with the aim of better understanding the dynamics of ALE and to test the simulation procedure. The simulations help to more fully interpret the experimental measurements. Further, optical emission measured just above the surface being etched can be related to etch products and can, therefore, be directly compared to simulation predictions. The simulations capture the measured initial product distribution leaving the surface and match the measured etch per cycle reasonably well. While simulations demonstrate the importance of ion-induced surface damage and mixing into a layer below the surface, the depth of which depends mainly on ion energy, the experiments suggest there is more Cl mixed into the layer than the MD procedure predicts.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Elastic Properties of Plasma-Exposed Tungsten Predicted by Molecular-Dynamics Simulations

We report results of systematic molecular-dynamics computations of the elastic properties of single-crystalline tungsten containing structural defects, voids and overpressurized He nanobubbles, related to plasma exposure of tungsten serving as a plasma- facing component (PFC) in nuclear fusion devices. Our computations reveal that the empty voids are centers of dilatation resulting in development of tensile stress in the tungsten matrix, whereas He-filled voids (nanobubbles) introduce compressive stress in the plasma-exposed tungsten. We find that the dependence of the elastic moduli of plasma-exposed tungsten, namely, the bulk, Young, and shear modulus, on its void fraction follows a universal exponential scaling relation. We also find that the elastic moduli of plasma-exposed tungsten soften substantially as a function of He content in the tungsten matrix, following an exponential scaling relation; this He-induced exponential softening is in addition to the softening caused in the matrix with increasing temperature. A systematic characterization of the dependence of the elastic moduli on the He bubble size reveals that He bubble growth affects significantly both the bulk modulus and the Poisson ratio of plasma-exposed tungsten, while its effect on the Young and shear moduli of the plasma-exposed material is weak. Our findings contribute directly to the development of a structure-properties database that is required for the predictive modeling of the dynamical response of PFCs in nuclear fusion devices.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Restrained Proton Indicator in Combined Quantum-Mechanics/Molecular-Mechanics Dynamics Simulations of Proton Transfer through a Carbon Nanotube

Recently, a collective variable “proton indicator” was purposed for tracking an excess proton solvated in bulk water in molecular dynamics simulations. In this work, we demonstrate the feasibility of utilizing the position of this proton indicator as a reaction coordinate to model an excess proton migrating through a hydrophobic carbon nanotube in combined quantum-mechanics/molecular-mechanics simulations. Our results indicate that applying a harmonic restraint to the proton indicator in the bulk solvent near the nanotube pore entrance leads to the recruitment of water molecules into the pore. This is consistent with an earlier study that employed a multistate empirical valence bond potential and a different representation (center of excess charge) of the proton. We attribute this water recruitment to the delocalized nature of the solvated proton, which prefers to be in high-dielectric bulk solvent. While water recruitment into the pore is considered an artifact in the present simulations (because of the artificially imposed restraint on the proton), if the proton were naturally restrained, it could assist in building water wires prior to proton transfer through the pore. The potential of mean force for a proton translocation through the water-filled pore was computed by umbrella sampling, where the bias potentials were applied to the proton indicator. The free energy curve and barrier heights agree reasonably with those in the literature. Finally, the results suggest that the proton indicator can be used as a reaction coordinate in simulations of proton transport in confined environments.

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

Final Project Report: Data Locality Enhancement of Dynamic Simulations for Exascale Computing

The goal of this project is to develop a set of techniques and software tools to enhance the matching between memory accesses in dynamic simulations and the prominent features of modern and future manycore systems, alleviating the memory performance issues for exascale computing. In the first three years, the PI and his group have achieves some significant progress towards the goal, producing a set of novel techniques for improving the memory performance and data locality in manycore systems, yielding 18 conference and workshop papers and 4 journal papers and graduating 6 Ph.Ds. This report summarizes the research results of this project through that period.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

EFFECT OF STRAIN FIELD ON THRESHOLD DISPLACEMENT ENERGY OF TUNGSTEN STUDIED BY MOLECULAR DYNAMICS SIMULATION

The influence of hydrostatic strain on point defect formation energy and threshold displacement energy (Ed) in body-centered cubic (BCC) tungsten was studied with molecular dynamics simulations. Two different tungsten potentials (Fikar and Juslin) were used. The minimum Ed direction calculated with the Fikar-potential was <100>, but with the Juslin-potential it was <111>. The most stable self-interstitial (SIA) configuration was a <111>-crowdion for both potentials. The stable SIA configuration did not change with applied strain. Varying the strain from compression to tension increased the vacancy formation energy but decreased the SIA formation energy. The SIA formation energy changed more significantly than for a vacancy such that Ed decreased with applied strain from compression to tension.

Wang, D.↗