Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular 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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Linking Friction Scales from Nano to Macro via Avalanches

Abstract Steady-state fluctuations in the friction force of molybdenum disulfide (MoS 2 ), a prototypical lamellar solid, were analyzed experimentally for newton-scale forces and computationally via molecular dynamics simulations for nanonewton-scale forces. A mean field model links the statics and the dynamics of the friction behavior across these eight orders of magnitude in friction force and six orders of magnitude in friction force fluctuations (i.e., avalanches). Both the statistics and dynamics of the avalanches match model predictions, indicating that friction can be characterized as a series of avalanches with properties that are predictable over a wide range of scales.

Salners, Tyler (ORCID:0000000332647944)↗
DOE OSTI2025

Support of Adhesion Mechanisms in Al 2 O 3 Aerosol Deposition Through Laser-Induced Particle Impact Testing

Aerosol deposition (AD) is a kinetic spray process capable of depositing ceramic coatings at room temperature, but AD process development is generally a laborious exploration of a large process parameter space. Here, this paper presents a case study investigating whether laser-induced particle impact testing (LIPIT) could be applied to expedite development of an alumina (Al 2 O 3 ) coating on nickel (Ni): Specifically, whether LIPIT measurements could predict critical velocities of adhesion on Ni and Al 2 O 3 , and the effect of ball milling the Al 2 O 3 powder. Because LIPIT has a diffraction-limited lower bound on imageable particle size, the usefulness of Al 2 O 3 powder agglomerates as a proxy for single particles was additionally studied. Overall, LIPIT measurements and AD sprays agreed that ball milling dramatically improves adhesion. Additionally, LIPIT measurements of critical velocity of adhesion of Al 2 O 3 powder agglomerates on Ni and Al 2 O 3 substrates (150 meters per second [m/s] and 250 m/s, respectively) quantitatively agreed with predictions from a previously published model based on picoindentation and molecular dynamics simulations. Together, these findings support the established hypothesis that Al 2 O 3 adheres via a dislocation-mediated mechanism in AD, that Al 2 O 3 powder agglomerates adhere as individual constituent particles rather than collectively, and that, for this case study, LIPIT measurements were predictive of AD process parameters.

Al2O3↗
DOE OSTI2021

Observation of Fundamental Mechanisms in Compression-Induced Phase Transformations Using Ultrafast X-ray Diffraction

As theoretically hypothesized for several decades in group IV transition metals, we have discovered a dynamically stabilized body-centered cubic (bcc) intermediate state in Zr under uniaxial loading at sub-nanosecond timescales. Under ultrafast shock wave compression, rather than the transformation from α-Zr to the more disordered hex-3 equilibrium ω-Zr phase, in its place we find the formation of a previously unobserved nonequilibrium bcc metastable intermediate. We probe the compression-induced phase transition pathway in zirconium using time-resolved sub-picosecond x-ray diffraction analysis at the Linac Coherent Light Source. We also present molecular dynamics simulations using a potential derived from first-principles methods which independently predict this intermediate phase under ultrafast shock conditions. In contrast with experiments on longer timescale (> 10 ns) where the phase diagram alone is an adequate predictor of the crystalline structure of a material, our recent study highlights the importance of metastability and time dependence in the kinetics of phase transformations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-Driven Insights into the Structural Essence of Plasticity in High-Entropy Alloys

The heterogeneous mechanical response of a crystalline alloy with multiple principal elements was investigated using molecular dynamics simulations. The local configuration of the alloy in its quiescent state was characterized by the variables derived from the gyration tensor and the atomic electronegativity. A multivariate analysis identified the geometric and chemical factors that influenced the atomic packing variations. Further, upon straining, the non-affine displacement exhibited spatial heterogeneity. A statistical correlation was established between the local yield events and the specific features of the local configuration. Our findings, validated by the performance metrics analysis, provided a structural criterion for the instability mechanisms in high-entropy alloys (HEAs) and enhanced the understanding of their plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards an Understanding of Oxidative Damage in an a-L-Arabinofuranosidase of Trichoderma reesei : a Molecular Dynamics Approach

Trichoderma reesei is a ‘workhorse’ fungus that produces glycosyl hydrolases (e.g., cellulases) at high titers for use in industrial bioprocessing. In this study, we focused on two side-chain hemicellulases, a-L-arabinofuranosidases and a-D-galactosidases, which are important for the treatment of lignocellulosic biomass, but susceptible to oxidative damage that can occur during industrial processing. The molecular details that render these enzymes inactive have not yet been identified. To approach this issue, we used proteomics to identify amino acid residues that were oxidized after a relevant oxidative treatment (Fenton reaction). These oxidative modifications were included in the 3D-protein structures, and using molecular dynamics simulations, we then studied the behaviors of non-modified and oxidized enzymes. These simulations showed significant alterations of the conformational stability of proteins when oxidized, as evidenced by changes in root mean square deviation (RMSD) and principal component analyses (PCA) trajectories. Likewise, enzyme-ligand interactions such as hydrogen bonds were greatly reduced in quantity and quality in the oxidized proteins. Finally, free energy landscape plots showed that there was a more rugged energy surface in the oxidized proteins, implying a less favorable reaction pathway. These results reveal the basis for loss of function in these carbohydrate active enzymes (CAZYs) in T. reesei.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Differences in Sb 2 Te 3 growth by pulsed laser and sputter deposition

High quality van der Waals chalcogenides are important for phase change data storage, thermoelectrics, and spintronics. Using a combination of statistical design of experiments and density functional theory, we clarify how out-of-equilibrium van der Waals epitaxial deposition methods can improve the crystal quality of Sb 2 Te 3 films. We compare films grown by radio frequency sputtering and pulsed laser deposition (PLD). The growth factors that influence the crystal quality for each method are different. For PLD grown films a thin amorphous Sb 2 Te 3 seed layer most significantly influences the crystal quality. In contrast, the crystalline quality of films grown by sputtering is rather sensitive to the deposition temperature and less affected by the presence of a seed layer. This difference is somewhat surprising as both methods are out-of-thermal-equilibrium plasma-based methods. Non-adiabatic quantum molecular dynamics simulations show that this difference originates from the density of excited atoms in the plasma. The PLD plasma is more intense and with higher energy than that used in sputtering, and this increases the electronic temperature of the deposited atoms, which concomitantly increases the adatom diffusion lengths in PLD. In contrast, the adatom diffusivity is dominated by the thermal temperature for sputter grown films. These results explain the wide range of Sb 2 Te 3 and superlattice crystal qualities observed in the literature. Finally, these results indicate that, contrary to popular belief, plasma-based deposition methods are suitable for growing high quality crystalline chalcogenides.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Short- and medium-range orders in Al90Tb10 glass and their relation to the structures of competing crystalline phases

Molecular dynamics simulations using an interatomic potential developed by artificial neural network deep machine learning are performed to study the local structural order in Al90Tb10 metallic glass. We show that more than 80% of the Tb-centered clusters in Al90Tb10 glass have short-range order (SRO) with their 17 first coordination shell atoms stacked in a ‘3661’ or ‘15551’ sequence. Medium-range order (MRO) in Bergman-type packing extended out to the second and third coordination shells is also clearly observed. Analysis of the network formed by the ‘3661’ and ‘15551’ clusters show that ~82% of such SRO units share their faces or vertexes, while only ~6% of neighboring SRO pairs are interpenetrating. Such a network topology is consistent with the Bergman-type MRO around the Tb-centers. Moreover, crystal structure searches using genetic algorithm and the neural network interatomic potential reveal several low-energy metastable crystalline structures in the composition range close to Al 90 Tb 10 . Some of these crystalline structures have the ‘3661’ SRO while others have the ‘15551’ SRO. While the crystalline structures with the ‘3661’ SRO also exhibit the MRO very similar to that observed in the glass, the ones with the ‘15551’ SRO have very different atomic packing in the second and third shells around the Tb centers from that of the Bergman-type MRO observed in the glassy phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nanotwin assisted reversible formation of low angle grain boundary upon reciprocating shear load

Severe plastic deformation of metals is known to lead to superior properties that cannot be achieved by any traditional metallurgic process. Origin of the superior properties is perceived to be closely associated with grain refinement, a fundamental process during the severe plastic deformation, which is essentially the formation of new grain boundaries. However, the atomistic mechanism of grain boundary formation remains largely obscure. Here, by using in-situ transmission electron microscopy and molecular dynamic simulation, we reveal, for the first time at atomic level, the shear-induced low-angle grain boundary (LAGB) formation processes in Au nanocrystal. In this work, we discover the LAGB formation is accomplished through inward propagation of nanotwins accompanied by dislocations gliding on the twin boundaries, a nanotwin-mediated dislocation slip mechanism, which shows reversible characteristic under reciprocating shear load and is affected by the nanocrystal microstructure and orientation. Our result unveils unprecedented atomistic insights on shear driven grain refinement towards nanostructure of superior properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic modeling of surface and grain boundary dislocation nucleation in FCC metals

Dislocation nucleation plays a critical role in the plastic deformation of crystalline materials. However, it is challenging to predict the active mode and associated rate of dislocation nucleation under typical experimental loading conditions through molecular dynamics simulation due to timescale limitations. In this report we use the free-end nudged elastic band method to determine the activation energies and activation volumes of dislocation nucleation in four typical face-centered cubic metals of Au, Al, Cu and Ni. We focus on the representative processes of surface and grain boundary dislocation nucleation. The atomistically determined activation volumes of these dislocation nucleation processes are larger than 10b3 (with b being the Burgers vector length) under typical experimental loading conditions. These results are compared with experimentally measured activation volumes in ultrafine-grained and nanocrystalline metals, thereby providing mechanistic insight into their rate-controlling deformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanical properties, failure mechanisms, and scaling laws of bicontinuous nanoporous metallic glasses

Molecular dynamics simulations are employed to study the mechanical properties of nanoporous Cu x Zr 1-x metallic glasses (MGs) with five different compositions, x = 0.28, 0.36, 0.50, 0.64, and 0.72, and porosity in the range 0.1 < Φ < 0.7. Results from tensile loading simulations indicate a strong dependence of Young's modulus, E, and Ultimate Tensile Strength (UTS) on porosity and composition. By increasing the porosity from Φ = 0.1 to Φ = 0.7, the topology of the nanoporous MG shifts from closed cell to open-cell bicontinuous. The change in nanoporous topology enables a brittle-to-ductile transition in deformation and failure mechanisms from a single critical shear band to necking and rupture of ligaments. Genetic Programming (GP) is employed to find scaling laws for E and UTS as a function of porosity and composition. Here, a comparison of the GP-derived scaling laws against existing relationships shows that the GP method is able to uncover expressions that can predict accurately both the values of E and UTS in the whole range of porosity and compositions considered.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Detecting Thermally-Induced Spinodal Decomposition with Picosecond Ultrasonics in Cast Austenitic Stainless Steels

Destructive techniques to monitor nuclear reactor component health may not always be available during service, as they are time-consuming and often require pre-installed inspection coupons. Non-destructive evaluation (NDE) techniques can bridge this gap by rapidly identifying the state of mission-critical reactor components, via inference between NDE-measurable material properties and those of ultimate interest, such as ductility and toughness. Here, we demonstrate one such inference about the health of thermally aged cast austenitic stainless steels. Observations of surface acoustic wave peak (SAW) splitting correlate with spinodal decomposition-induced embrittlement as destructively measured by Charpy impact energy. Elastodynamic calculations and molecular dynamics simulations of the effects of spinodal decomposition on elastic moduli support that the new acoustic modes present are due to stiffening in the δ-ferrite domains. Finally, this discovery enables one to probe structure-property relationships in materials in a greatly accelerated manner, suggesting that similar inference methods can be used to determine material fitness-for-service, or to quickly uncover new structure-property relationships.

304-type stainless steel↗
DOE OSTI2023

High-temperature active oxidation of nanocrystalline silicon-carbide: A reactive force-field molecular dynamics study

Flexible woven SiC ceramics are prone to accelerated fiber embrittlement under high temperature oxidation in dynamic oxygen environments. The nanocrystalline structure of the constituent fibers impacts the reaction kinetics and phase transformations during active oxidation. However, fundamental understanding and quantification of grain boundary effects on oxidation behavior in nanocrystalline SiC remain elusive when temperatures exceed 1500 K. This study deploys large-scale molecular dynamics simulations with a reactive force-field to elucidate the complex roles of atomic oxygen reservoir conditions and grain size on oxidation kinetics and the nature of oxides produced in both monocrystalline and nanocrystalline 3C-SiC between 1100 K and 2000 K. The simulations with dynamically replenished oxygen provide good agreement with oxidation kinetics and activation energies for the monocrystalline Si(100) and C(100) orientations published in the available literature. This study reveals that, by contrast, nanocrystalline SiC samples exhibit two distinct oxidation kinetics with a transition point at 1500 K due to surface melting, which is supported by experimental evidence. The introduction of a grain-boundary network produces a two-fold decrease in oxidation activation energies compared to monocrystalline SiC below 1500 K. Above 1500 K, however, the activation energies rise substantially due to the formation of a liquid Si phase at the SiC/Si oxide interface. Further, it is shown that the stability of the interfacial liquid phase is promoted by incoherent grain boundaries in the crystalline SiC. These findings are important for the deployment of nanocrystalline SiC fibers in advanced thermal protection systems for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Athermal annealing of pre-existing defects in crystalline silicon

Systematic investigations of electronic energy loss (S e ) effects on pre-existing defects in crystalline silicon (Si) are crucial to provide reliance on the use of ionizing irradiation to anneal pre-existing defects, leading to successful implementation of this technology in the fabrication of Si-based devices. In this regard, the S e effects on nonequilibrium defect evolution in pre-damaged Si single crystals at 300 K has been investigated using intermediate-energy ions (12 MeV O and Si ions) that interact with the pre-damaged surface layers of Si mainly by ionization, except at the end of their range where the nuclear energy loss (S n ) is no longer negligible. Furthermore, under these irradiation conditions, experimental results and molecular dynamics simulations have revealed that pre-existing disorder in Si can be almost fully annealed by subsequent irradiation with intermediate-energy incident ions with S e values as low as 1.5–3.0 keV/nm. Selective annealing of pre-existing defect levels in Si at room temperature can be considered as an effective strategy to mediate the transient enhanced diffusion of dopants in Si. This approach is more desirable than the regular thermal annealing, which is not compatible with the processing requirements that fall below the typical thermal budget.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mobility of twinning dislocations in copper up to supersonic speeds

Understanding the mobility of twinning dislocations is important for multiscale modeling of crystal plasticity, especially at high strain rates, where such dislocations may reach transonic or supersonic speeds. Here, we used molecular dynamics simulations to investigate the relationship between dislocation velocity and the applied resolved shear stress of an edge twinning dislocation in copper up to supersonic speeds. The twinning dislocation mobility relation is composed of two branches separated by a band of forbidden velocities. The lower velocity branch is limited by the first transverse sound speed ~2000 m/s while the upper branch stretches from ~3500 m/s in the transonic regime to supersonic velocities. Twinning dislocations cannot undergo uniform steady-state motion at velocities within the forbidden band. Our simulation results also reveal that edge twinning dislocation motion in copper is kink-mediated. We discuss the implications of our findings for the motion of twins, twinning dislocations, and twinning dislocation kinks in copper.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring NaCl-PuCl 3 molten salts with machine learning interatomic potentials and graph theory

Actinide molten salts are the basis of the liquid fuels used in molten salt reactors. Due to the inherent difficulties associated with high temperature and hazardous conditions, experimental investigations of fundamental properties of these materials are usually challenging. In this work, we describe the structure and transport of NaCl-PuCl 3 mixtures using computational techniques. Three compositions were considered (16, 25, and 36 mol% PuCl 3 ) over a temperature range (730 – 1257K) using ab initio molecular dynamics, which provided the necessary data sets for training machine learned interatomic potentials. Further, molecular dynamics simulations based on these potentials were then used to determine structure and transport properties. A substantial change was noted in the structure factor when increasing the PuCl 3 content from 25 to 36 mol%. This change is linked to the aggregation of larger Pu 3+ clusters. In addition, the similarity of the atomic environments of metal cations in molten salt systems to their solid states counterparts was investigated using an unsupervised learning technique. Finally, graph theory was employed to explore the structure and size of actinide networks. Consistent with the structure factor, a dense Pu 3+ intermolecular structure is observed within the 36 mol% PuCl 3 mixture. The structure of cation-cation inter-junctions is also discussed. In all cases, the diffusion of Pu 3+ is significantly lower than that of Na + and Cl - .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deformation behavior and phase transformation of nanotwinned Al/Ti multilayers

We report nanotwinned Al/Ti multilayers have exhibited size-dependent microstructure evolution and high strength. However, their deformation mechanisms are less well understood. In this work, we investigated the deformation mechanisms of nanotwinned Al/Ti multilayers with FCC/HCP layer interfaces by using in situ micropillar compression tests. Nanotwinned Al/Ti multilayers exhibit compressive strength up to 2.4 GPa and good work hardening capability. Post-compression TEM analyses reveal high-density stacking faults and the HCP-to-FCC phase transformations in Ti. Molecular dynamics simulations elucidate the mechanisms of deformation induced phase transformation in Ti and the influence of collective movement of partial dislocations on the deformability of Al/Ti multilayers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formation of pyrophosphates across grain boundaries induces the formation of mismatched but oriented interfaces in silver phosphate polypods

The interfaces and their misfit defects determine the materials properties for a wide range of applications, such as electronic devices, photocatalysis, and mechanical engineering, etc. However, current understanding of atomic interfacial structures is limited. Here we discover a special interfacial structure, mismatched but oriented interface via two distinct facets. Transmission and scanning electron microscopy results suggest that, in Ag 3 PO 4 polypods structures, interfaces of {1 0 0} and {1 1 0}, ({1 0 0}/(1 1 0}), {1 0 0}/{1 1 1}, {1 1 0}/{1 1 1}, and {1 0 0}/{1 0 0}, etc., have a certain orientation relationship, corresponding to the energy minima and coincident site lattice of interfacial atoms as demonstrated by molecular dynamics simulations. Density functional theory demonstrates that the formation of pyrophosphate and/or phosphates rotation to bond across the interface compensate the lattice mismatch at the interfaces, as well as deformations of Ag-O bonds. Furthermore, our work opens up a new avenue for a much wider range of interfacial structures, allow for a higher diversity of structures, and shine light on tailoring crystal structures, morphologies, and the resulting properties.

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

Thermal energy evolution and mechanical deformation of monocrystalline yttria-stabilized zirconia nanoparticles in aerosol deposition processes

Aerosol deposition (AD) is a coating process wherein aerosol particles are impacted on a target substrate. There are fundamental differences between the AD process (cold impact), where particle translational kinetic energy is high and thermal energy is low, and thermal spray deposition (thermal impact), where translational energy is lower but thermal energy is high. To better compare cold and thermal impact effects on particles, we carried out molecular dynamics simulations using yttria-stabilized zirconia (YSZ) nanoparticles on YSZ substrates as a model system. We performed cold impact simulations at 300 K with variable impact velocity in the 500 ms –1 –1500 ms –1 range to understand how increasing translational kinetic energy affects thermal energy and mechanical evolution. We then performed thermal impact simulations at variable temperature and impact velocity, but where the total kinetic energy of the nanoparticle was equivalent to that of a 300 K, 1000 impact. In cold impact, the temperature increases in YSZ nanoparticles at a rate of 10 13 – 10 14 K s –1 , and large temperature gradients result. Conversely, in thermal impact, nanoparticle temperatures remain uniform. Furthermore, the temperature gradients during cold impact coincide with plastic deformation in nanoparticles, while with larger thermal energies, plastic deformation is reduced.

42 ENGINEERING↗