Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamics simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Molecular Dynamics Simulation of Fission Fragment Damage in Nuclear Fuel and Surrogate Material

We have performed classical molecular dynamics simulations of swift heavy ion damage, typical of fission fragments, in nuclear fuel (UO2) for energy deposition per unit length of 3.9 keV/nm. We did not observe amorphization. The damage mainly consisted of isolated point defects. Only about 1% of the displacements occur on the uranium sublattice. Oxygen Frenkel pairs are an order of magnitude more numerous than uranium Frenkel pairs in the primary damage state. In contrast, previous results show that the ratio of Frenkel pairs on the two sublattices is close to the stoichiometric ratio in ceria. These differences in the primary damage state may lead to differences in radiation response of UO2 and CeO2.

Molecular dynamics, radiation damage in materials,↗
DOE OSTI2025

Differentiable hybrid neural network approach for enhancing reactor dynamics simulations

Reactor dynamics simulations provide essential insights into the time-dependent behavior of nuclear reactors under various operating conditions. However, high-fidelity simulations can be computationally intensive, requiring significant computational resources. Here, to address this challenge, this study employs a differentiable hybrid model that utilizes neural networks as a corrector to enhance the performance of a low-fidelity simulation, aligning its predictions with those of a high-fidelity simulation. Low-fidelity and high-fidelity simulations were obtained by adjusting the mesh size in the System Dynamics Analysis Tool. The differentiable hybrid model was trained in two approaches: time-step-wise and sequence-wise. It was then applied to simulate various transients in a molten salt reactor. Its performance was evaluated by comparing its responses to transients against those of the high-fidelity simulation. An additional approach was performed using a data-driven model to correct the low-fidelity simulation. In comparison, the differentiable hybrid model showed significant improvements in transient prediction, effectively addressing the limitations of the low-fidelity simulations. The results highlighted the robustness of the differentiable hybrid model in both training approaches. It delivered simulations that were at least 3.8 times faster than high-fidelity models. In the time-step-wise approach, it achieved at least a 39% improvement in accuracy. In the sequence-wise approach, it showed at least an 81% accuracy improvement over the full transient. This approach offers a promising path for improving computational efficiency without compromising accuracy in nuclear reactor simulations, making it suitable for real-time digital twin applications.

42 - ENGINEERING↗
DOE OSTI2015

Testing thermal gradient driving force for grain boundary migration using molecular dynamics simulations

Strong thermal gradients in low-thermal-conductivity ceramics may drive extended defects, such as grain boundaries and voids, to migrate in preferential directions. In this work, molecular dynamics simulations are conducted to study thermal gradient driven grain boundary migration and to verify a previously proposed thermal gradient driving force equation, using uranium dioxide as a model system. It is found that a thermal gradient drives grain boundaries to migrate up the gradient and the migration velocity increases under a constant gradient owing to the increase in mobility with temperature. Different grain boundaries migrate at very different rates due to their different intrinsic mobilities. The extracted mobilities from the thermal gradient driven simulations are compared with those calculated from two other well-established methods and good agreement between the three different methods is found, demonstrating that the theoretical equation of the thermal gradient driving force is valid, although a correction of one input parameter should be made. The discrepancy in the grain boundary mobilities between modeling and experiments is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Classical molecular dynamics simulation of electronically non-adiabatic processes

Both classical and quantum mechanics (as well as hybrids thereof, i.e., semiclassical approaches) find widespread use in simulating dynamical processes in molecular systems. For large chemical systems, however, which involve potential energy surfaces (PES) of general/arbitrary form, it is usually the case that only classical molecular dynamics (MD) approaches are feasible, and their use is thus ubiquitous nowadays, at least for chemical processes involving dynamics on a single PES (i.e., within a single Born-Oppenheimer electronic state). Our paper reviews recent developments in an approach which extends standard classical MD methods to the treatment of electronically non-adiabatic processes, i.e., those that involve transitions between different electronic states. This approach treats nuclear and electronic degrees of freedom (DOF) equivalently (i.e., by classical mechanics, thereby retaining the simplicity of standard MD), and provides "quantization" of the electronic states through a symmetrical quasi-classical (SQC) windowing model. The approach is seen to be capable of treating extreme regimes of strong and weak coupling between the electronic states, as well as accurately describing coherence effects in the electronic DOF (including the de-coherence of such effects caused by coupling to the nuclear DOF). A survey of recent applications is presented to illustrate the performance of the approach. We also described is a newly developed variation on the original SQC model (found universally superior to the original) and a general extension of the SQC model to obtain the full electronic density matrix (at no additional cost/complexity).

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

Protein–Mineral Interactions: Molecular Dynamics Simulations Capture Importance of Variations in Mineral Surface Composition and Structure

Molecular dynamics simulations, conventional and metadynamics, were performed to determine the interaction of model protein Gb1 over kaolinite (001), Na+-montmorillonite (001), Ca2+-montmorillonite (001), goethite (100), and Na+-birnessite (001) mineral surfaces. Gb1, a small (56 residue) protein with a well-characterized solution-state nuclear magnetic resonance (NMR) structure and having α-helix, four-fold β-sheet, and hydrophobic core features, is used as a model protein to study protein soil mineral interactions and gain insights on structural changes and potential degradation of protein. From our simulations, we observe little change to the hydrated Gb1 structure over the kaolinite, montmorillonite, and goethite surfaces relative to its solvated structure without these mineral surfaces present. Over the Na+-birnessite basal surface, however, the Gb1 structure is highly disturbed as a result of interaction with this birnessite surface. Unraveling of the Gb1 β-sheet at specific turns and a partial unraveling of the α-helix is observed over birnessite, which suggests specific vulnerable residue sites for oxidation or hydrolysis possibly leading to fragmentation.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2014

Development of EEM based silicon–water and silica–water wall potentials for non-reactive molecular dynamics simulations

Molecular dynamics simulations of heat transfer in gases are computationally expensive when the wall molecules are explicitly modeled. To save computational time, an implicit boundary function is often used. Steele's potential has been used in studies of fluid–solid interface for a long time. In this work, the conceptual idea of Steele's potential was extended in order to simulate water–silicon and water–silica interfaces. A new wall potential model is developed by using the electronegativity-equalization method (EEM), a ReaxFF empirical force field and a non-reactive molecular dynamics package PumMa. Contact angle simulations were performed in order to validate the wall potential model. Contact angle simulations with the resulting tabulated wall potentials gave a silicon–water contact angle of 129°, a quartz–water contact angle of 0°, and a cristobalite–water contact angle of 40°, which are in reasonable agreement with experimental values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Mean ionic activity coefficients in aqueous NaCl solutions from molecular dynamics simulations

The mean ionic activity coefficients of aqueous NaCl solutions of varying concentrations at 298.15 K and 1 bar have been obtained from molecular dynamics simulations by gradually turning on the interactions of an ion pair inserted into the solution. Several common non-polarizable water and ion models have been used in the simulations. Gibbs-Duhem equation calculations of the thermodynamic activity of water are used to confirm the thermodynamic consistency of the mean ionic activity coefficients. While the majority of model combinations predict the correct trends in mean ionic activity coefficients, they overestimate their values at high salt concentrations. The solubility predictions also suffer from inaccuracies, with all models underpredicting the experimental values, some by large factors. These results point to the need for further ion and water model development.

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

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗
DOE OSTI2015

Extended Lagrangian Born-Oppenheimer molecular dynamics simulations of the shock-induced chemistry of phenylacetylene

The initial chemical events that occur during the shock compression of liquid phenylacetylene have been investigated using self-consistent tight binding molecular dynamics simulations. The extended Lagrangian Born-Oppenheimer molecular dynamics formalism enabled us to compute microcanonical trajectories with precise conservation of the total energy. Our simulations revealed that the first density-increasing step under shock compression arises from the polymerization of phenylacetylene molecules at the acetylene moiety. The application of electronic structure-based molecular dynamics with long-term conservation of the total energy enabled us to identify electronic signatures of reactivity via monitoring changes in the HOMO-LUMO gap, and to capture directly adiabatic shock heating, transient non-equilibrium states, and changes in temperature arising from exothermic chemistry in classical molecular dynamics trajectories.

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

A Case Study on Pathogen Transport, Deposition, Evaporation and Transmission: Linking High-Fidelity Computational Fluid Dynamics Simulations to Probability of Infection

A high-fidelity, low-Mach computational fluid dynamics simulation tool that includes evaporating droplets and variable-density turbulent flow coupling is well-suited to ascertain transmission probability and supports risk mitigation methods development for airborne infectious diseases such as COVID-19. A multi-physics large-eddy simulation-based paradigm is used to explore droplet and aerosol pathogen transport from a synthetic cough emanating from a kneeling humanoid. For an outdoor configuration that mimics the recent open-space social distance strategy of San Francisco, maximum primary droplet deposition distances are shown to approach 8.1 m in a moderate wind configuration with the aerosol plume transported in excess of 15 m. In quiescent conditions, the aerosol plume extends to approximately 4 m before the emanating pulsed jet becomes neutrally buoyant. A dose–response model, which is based on previous SARS coronavirus (SARS-CoV) data, is exercised on the high-fidelity aerosol transport database to establish relative risk at eighteen virtual receptor probe locations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Understanding Viscoelasticity of an Entangled Silicone Copolymer via Coarse-Grained Molecular Dynamics Simulations

Entangled dynamics is important for understanding rheological properties of long-chain polymers. For entangled homopolymers, the classic tube-reptation model and its refinements have been successfully applied to quantify properties like diffusion coefficient and zero-rate viscosity. However, the application of such models to copolymers has been limited despite scientific and industrial importance. Here, we study the entangled melt dynamics of poly(dimethyl-co-diphenyl)siloxane random copolymer for a range of mean-composition-ratio ϕ of the diphenyl component via long-term molecular dynamics simulation with a recently developed coarse-grained model. We found that the segmental relaxation is heterogeneous at the monomeric level because of compositional fluctuations. However, at the chain-entanglement level and higher length scales, the viscoelastic response is homogeneous with compositional dependence only through the overall diphenyl fraction ϕ. The relaxation modulus of the entangled copolymer melt conforms to the Likhtman–McLeish model, and the viscosity predicted using our current coarse-grained parameters is in good quantitative agreement with experimental data.

Copolymers↗
DOE OSTI2015

Growth of bi- and tri-layered graphene on silicon carbide substrate via molecular dynamics simulation

Molecular dynamics (MD) simulation with simulated annealing method is used to study the growth process of bi- and tri-layered graphene on a 6H-SiC (0001) substrate via molecular dynamics simulation. Tersoff-Albe-Erhart (TEA) potential is used to describe the inter-atomic interactions among the atoms in the system. The formation temperature, averaged carbon-carbon bond length, pair correlation function, binding energy and the distance between the graphene formed and the SiC substrate are quantified. The growth mechanism, graphitization of graphene on the SiC substrate and characteristics of the surface morphology of the graphene sheet obtained in our MD simulation compare well to that observed in epitaxially grown graphene experiments and other simulation works.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Characterization of radiation damage in TiO 2 using molecular dynamics simulations

Molecular dynamics simulations are carried out to characterize irradiation effects in TiO 2 rutile, for wide ranges of temperatures (300–900 K) and primary knock-on atom (PKA) energies (1–10 keV). The number of residual defects decreases with increased temperature and decreased PKA energy, but is independent of PKA type. In the ballistic phase, more oxygen than titanium defects are produced, however, the primary residual defects are titanium vacancies and interstitials. Defect clustering depends on the PKA energy, temperature, and defect production. For some 10 keV PKAs, the largest cluster of vacancies at the peak of the ballistic phase and after annealing has up to ≈1200 and 100 vacancies, respectively. For the 10 keV PKAs at 300 K, the energy storage, primarily in residual Ti vacancies and interstitials, is estimated at 140–310 eV. It decreases with increased temperature to as little as 5–180 eV at 900 K. In conclusion, selected area electron diffraction patterns and radial distribution functions confirm that although localized amorphous regions form during the ballistic phase, TiO 2 regains full crystallinity after annealing.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Bergman-type medium range order in amorphous Zr77Rh23 alloy studied by ab initio molecular dynamics simulations

In recent years, some arguments about the existence of medium-range order (MRO) in the Zr-Rh system have been put forward. However, research on the structural features of the Zr-Rh binary alloy at the atomic level is still lacking. This study uses ab initio molecular dynamics simulations to systematically study the local structures of Zr77Rh23 from the liquid to the glassy states. Pair correlation function (PCF), coordination number (CN), Honeycutt–Anderson(HA) index, bond-angle distribution functions, and the Voronoi tessellation method are used to reveal a clear icosahedral-like configuration in the amorphous Zr77Rh23 alloy. It is noteworthy that the splitting in the second peak of the partial PDF implies the existence of a medium range order (MRO) in the Zr77Rh23 system. We obtain the local order in three-dimensional atomic density distributions by using a new atomistic cluster alignment (ACA) method. Interestingly, a Bergman-type MRO is observed in the glassy Zr77Rh23. Furthermore, the spatial distribution and connections of icosahedral-like clusters are shown to further demonstrate the MRO network. Our findings shed light on the nature of atomic local structures of amorphous Zr77Rh23 alloy and have important implications to understanding the formation of various MROs in metallic glasses.

Guo, Y. R.↗
DOE OSTI2020

Ab initio molecular dynamics simulation of Nd3+ incorporation in calcite

Trivalent rare-earth elements (REEs) readily incorporate in Ca-bearing minerals such as calcite (CaCO3), a phenomenon with important implications for the use of REEs as biogeochemical tracers as well as for environmental remediation. Despite substantial efforts, determining the incorporation modes of trivalent REEs in calcite has remained challenging due to their presence in multiple coordination environments. This is particularly problematic for the interpretation of extended X-ray absorption fine structure (EXAFS) spectroscopy. In this work, ab initio molecular dynamics (AIMD) simulations of Nd3+ incorporated in calcite were performed to resolve its incorporation modes. Nd3+ and Ca2+ have very similar ionic radii but the additional positive charge requires a charge compensation scheme (CCS). Five CCSs were considered in AIMD simulations, and the resulting atomic trajectories were then used to compute EXAFS spectra for direct fitting to a published experimental spectrum. The results of the fits indicated that Nd3+ incorporates at the Ca2+ site in calcite in a mixture of six- and seven-fold coordination environments and through three main incorporation modes: associated with a protonated nearest-neighbor Ca2+ vacancy, incorporated as NdOH2+, and charge compensated by a remote species (on the spatial scale probed by EXAFS). The incorporation of Nd3+ in calcite via multiple modes was consistent with previous work, but the AIMD-EXAFS approach allowed for resolving the nature, proportion, and structure of the different coordination environments.

Kerisit, Sebastien N.↗
DOE OSTI2021

Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data

The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.

Prange, Micah P.↗