Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Short- to Intermediate-Range Structure, Transport, and Thermophysical Properties of LiF–NaF–ZrF4 Molten Salts

LiF–NaF–ZrF 4 multicomponent molten salts are identified as promising candidates for coolant salts in molten salt reactors and advanced high-temperature reactors. This study focused on low-melting point salt compositions of interest: 38LiF–51NaF–11ZrF 4 , 42LiF–29NaF–29ZrF 4 , and 26LiF–37NaF–37ZrF 4 . Ab-initio molecular dynamics (AIMD) calculations were performed and compared with available experimental data to assess the ability of rigid ion models (RIM) to reproduce short- to intermediate-range structure, transport, and thermophysical properties of the LiF–NaF–ZrF 4 salt mixtures. It is found that as ZrF 4 mol% increases, the average cation–anion coordination number (CN) of monovalent cations (Li + , Na + ) obtained from RIM calculations decreases, while multivalent Zr 4+ CN varied from 15% to 19% in comparison to corresponding AIMD values. In addition, RIM is found to predict the existence of 7, 8, and 9 coordinated fluorozirconate complexes, while AIMD and the available experimental data showed an occurrence of 6, 7, and 8 coordinated complexes in the melt. The intermediate-range structure analysis revealed that while the RIM parameters are able to reproduce a local structure for lower ZrF 4 mol% salts such as in 38LiF–51NaF–11ZrF 4 , an extensive fluorozirconate network formation is observed in RIM simulations for higher ZrF 4 mol% compositions. The network generated by RIM parameters is found to be mainly connected by “corner-sharing” fluorozirconate complexes as opposed to both “edge-sharing” and “corner-sharing” connectively portrayed by AIMD. It is found that a close agreement between AIMD and the RIM salt structure for the 11-mol% ZrF 4 salt resulted in good agreement in the calculated Zr diffusivities and the viscosity values. However, due to the inaccurate short- to intermediate-range structure prediction by RIM for higher ZrF 4 mol% compositions, thermophysical properties such as densities and heat capacity differ by up to 26% and 27%, respectively, upon comparison with AIMD and experimental values. Also, the network-dominated properties such as diffusion coefficients and viscosities differed by up to two and three orders of magnitude, respectively. This study signifies the importance of accurate salt structure generation for an accurate prediction of transport and thermophysical properties of multicomponent molten salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Fundamental microscopic properties as predictors of large-scale quantities of interest: Validation through grain boundary energy trends

Correlations between fundamental microscopic properties computable from first principles, which we term canonical properties, and complex large-scale quantities of interest (QoIs) provide an avenue to predictive materials discovery. Here, we propose that such correlations can be efficiently discovered through simulations utilizing approximate interatomic potentials (IPs), which serve as an ensemble of “synthetic materials”. As a proof of principle we build a regression model relating canonical properties to the symmetric tilt grain boundary (GB) energy curves in face-centered cubic crystals, characterized by the scaling factor in the universal lattice matching model of Runnels et al. (2016), which we take to be our QoI. Our analysis recovers known correlations of GB energy to other properties and discovers new ones. We also demonstrate, using available density functional theory (DFT) GB energy data, that the regression model constructed from IP data is consistent with DFT results, confirming the assumption that the IPs and DFT belong to same statistical pool and thereby validating the approach. Regression models constructed in this fashion can be used to predict large-scale QoIs based on first-principles data and provide a general method for training IPs for QoIs beyond the scope of first-principles calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Decoupling Proton and Cation Contributions to Capacitive Charge Storage in Birnessite in Aqueous Electrolytes

Abstract Nanostructured birnessite is of interest as an electrode material for aqueous high power electrochemical energy storage as well as desalination devices. In neutral pH aqueous electrolytes, birnessite exhibits a capacitive response attributed to the adsorption of cations and protons at the outer surface and within the hydrated interlayer. Here, we utilize the understanding of proton‐coupled electron transfer (PCET) in buffered electrolytes to decouple the role of protons and cations in the capacitive charge storage mechanism of birnessite at neutral pH. We find that without buffer, birnessite exhibits primarily potential‐independent (capacitive) behavior with excellent cycling stability. Upon the addition of buffer, the capacity initially increases and the cyclic voltammograms become more potential‐dependent, features attributed to the presence of PCET with the birnessite. However, long‐term cycling in the buffered electrolyte leads to significant capacity fade and dissolution, which is corroborated through ex situ characterization. ReaxFF atomistic scale simulations support the observations that proton adsorption leads to birnessite degradation and that capacitive charge storage in birnessite is primarily attributed to cation adsorption at the outer surface and within the interlayer.

Saeed, Saeed↗
DOE OSTI2021

Grain size dependencies of intergranular solute segregation in nanocrystalline materials

Grain boundary segregation is of increasing importance for stabilizing nanocrystalline alloys, and while segregation is typically modeled as a scale-independent phenomenon, at the finest nanocrystalline grain sizes there are several sources of grain size dependence. Here we provide quantitative detail on two grain size dependencies not previously elaborated, both of which arise from the spectrality of intergranular segregation sites. The first of these arises because a change in grain size affects which portions of the segregation spectrum are occupied by the solutes, with no change in the shape or nature of the spectrum itself. The second of these arises from the presence of higher-order grain junctions, which have unique segregation behaviors and an increasing relative presence at finer grain sizes. This paper uses molecular statics to quantify these size dependencies in detail for the case of an Al(Mg) binary alloy. Together, these two effects combined have a significant effect on the effective McLean segregation energy of more than 6 kJ/mol (from approximately -18.5 kJ/mol at 5 nm to -12.5 kJ/mol at more than 40 nm), with the majority of grain size dependence arising from the spectrality of segregation sites.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calcium silicate hydrate colloid at different humidities: Microstructure, deformation mechanism, and mechanical properties

Calcium silicate hydrate (C-S-H) is the main hydration product of Portland cement, while fundamental questions persist regarding its colloidal structure and water content. In this work, we developed a full atomistic C-S-H colloid model and investigated the physical and mechanical properties at different gel water concentrations. Water confined in nanometer gel pores degrades the mechanical properties of C-S-H because of its hydrolysis effect. A “ductile to brittle” transition was observed in the C-S-H colloid model with increasing gel water contents, evolving from transgranular failure to intergranular fracture. The gel water confined in nanometer-sized pores displays glassy characteristics, and it acts as a lubricant and facilitates structure distortions during the deformation. We elucidated the topology effect of the hydrogen bond network on the evolution of the viscoelasticity of C-S-H at different humidities. In conclusion, this work sheds light on the role of water in C-S-H colloids from nanostructure perspectives.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Energetics of the nucleation and glide of disconnection modes in symmetric tilt grain boundaries

Grain boundaries (GBs) evolve by the nucleation and glide of disconnections, which are dislocations with a step character. In this work, motivated by recent success in predicting GB properties such as the shear coupling factor and mobility from the intrinsic properties of disconnections, we develop a systematic method to calculate the energy barriers for the nucleation and glide of individual disconnection modes under arbitrary driving forces and a quasi-2D setting. This method combines tools from bicrystallography to enumerate disconnection modes and the Nudged elastic band (NEB) method to calculate their energetics, yielding minimum energy paths and atomistic mechanisms for the nucleation and glide of each disconnection mode. We apply the method to accurately predict shear coupling factors of $[001]$ symmetric tilt grain boundaries in Cu. Particular attention is paid to the boundaries where the dislocation-based disconnection nucleation model produces incorrect nucleation barriers. We demonstrate that the method can accurately compute energy barriers and predict shear-coupling factors in the low-temperature regime. For certain disconnection modes in which the assumptions underlying our method do not hold, we report upper bounds on the energy barriers for disconnection nucleation and glide. In addition, the NEB trajectories reveal interesting phenomena such as the dissociation of a higher energy mode into lower energy modes, and in some cases, shear coupling being mediated by partial disconnections, wherein the GB structure temporarily changes to a metastable state before reverting back to its original structure. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low-temperature carbonization of polyacrylonitrile/graphene carbon fibers: A combined ReaxFF molecular dynamics and experimental study

Graphene inclusion in a polymer matrix is a promising route to significantly enhance the mechanical properties of low-grade carbon fibers (CFs). Using ReaxFF molecular dynamics simulation, the atomistic mechanism leading to this enhancement is investigated. We demonstrate that the graphene edges along with the nitrogen and oxygen functional groups play a catalytic role and act as seeds to expedite alignment of the all-carbon rings, which are starting sites for the growth of graphitic structures. To examine the role of this proposed mechanism that enhances the graphitic structure of PAN/graphene CFs, we discuss the experimental results wherein the PAN/graphene CFs carbonized at 1250 C demonstrate 91% (from 632 to 1207 MPa) increase in strength and 101% (from 88 to 177 GPa) enhancement in Young’s modulus compared to PAN-based CFs carbonized at 1500 C. In conclusion, these enhanced mechanical properties of low-grade carbon fibers achieved via graphene inclusion at decreased carbonization temperature provide a means to realize both energy savings and cost reduction.

42 ENGINEERING↗
DOE OSTI2023

Soliton motion induced along ferromagnetic skyrmion chains in chiral thin nanotracks

Using atomistic magnetic simulations we investigate the soliton motion along a pinned skyrmion chain containing an interstitial skyrmion. We find that the soliton can exhibit stable motion along the chain without a skyrmion Hall effect for an extended range of drives. Under a constant drive the solitons have a constant velocity. We also measure the skyrmion velocity–current curves and identify the signatures of different phases including a pinned phase, stable soliton motion, and quasi-free motion at higher drives where all of the skyrmions depin from the pinning centers and move along the rigid wall. In the quasi-free motion regime, the velocity is oscillatory due to the motion of the skyrmions over the pinning sites. For increasing pinning strength, the onset of soliton motion shifts to higher values of current density. We also find that for stronger pinning, the characteristic velocity–current shape is affected by the annihilation of single or multiple skyrmions in the drive interval over which the soliton motion occurs. Our results indicate that stable skyrmion soliton motion is possible and that the solitons could be used as information carriers instead of the skyrmions themselves for technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Local damage in grain boundary stabilized nanocrystalline aluminum

The excess free volume in grain boundaries of pure nanocrystalline aluminum leads to the grain coarsening, which degrades their desirable mechanical properties. The preferential segregation of dopants to the grain boundaries has been proposed to enhance the stability of the nanostructure. Here, to shed light on the stabilization mechanisms, we present results from virtual mechanical tests on nanocrystalline aluminum samples with similar grain sizes but different magnesium dopant contents. The segregation of magnesium decreases the excess free volume within the grain boundaries and prohibits the nucleation and growth of intergranular cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nanotwin stability in alloyed copper under ambient and cryo-temperature dependent deformation states

Nanometer size grains, solid solution strengthening, low temperatures, and a high density of nanotwins can provide significant strengthening to a material. Nevertheless, under loading the microstructure may evolve altering the effect of grain size, twin density, etc. on this strengthening. In some cases, these effects can be coupled as in detwinning, which is a mechanism for cryogenic grain growth. Here, in this work, a series of copper nanotwinned columnar grains are solid solution alloyed with aluminum and subject to ambient and cryogenic (–196 °C) environments during mechanical indentation, with the influence of alloying on the detwinning evolution studied. Under indentation, at a lower aluminum solute content (2 at.%), the columnar grains bent corresponding to a loss of twins. At a higher aluminum solute content (8 at.%), the grains shear and retain the nanotwins. Atomistic-scale simulations reveal an increase in the incoherent twin boundary velocity with decreasing temperature, but this boundary velocity is reduced with increasing alloy solute content. The coupling of solid solution strengthening and nanotwin stabilization is then used to explain the microstructural stability as a function of indentation at different temperature conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-temperature lean Cu alloys with Cr-to-Nb atomic ratio of 2

Two Cu-Cr-Nb alloys, denoted as alloy 1 (comprising Cu-0.89 at% Cr-0.42 at% Nb) and alloy 2 (comprising Cu-1.84 at% Cr-0.99 at% Nb), were produced through a series of manufacturing processes including vacuum induction melting, melt spinning, consolidation, brazing, and baking, with both alloys aimed at achieving a nominal Cr-to-Nb atomic ratio of 2. Microstructural characterization using transmission electron microscopy and X-ray diffraction identified the cubic C15 Laves-phase Cr 2 Nb as the dominant precipitate in both alloys, cross-validated by thermodynamic calculations and atomistic simulation-based density functional theory (DFT). Besides cubic C15 Cr 2 Nb, hexagonal C14-phase Cr 2 Nb and α-BiF3 cubic structured Cr 3 Nb were also observed in the alloys, including a coherent interface formed between the Cr 3 Nb precipitate and the Cu matrix. The hardness of the alloys increases, and the electrical conductivity decreases with increasing alloying addition content; two practical equations described the trends. Further DFT simulations revealed that the electrical conductivity (conductance) of the Cu/Cr 2 Nb interface is an order of magnitude higher than the intrinsic Cu high-angle grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Insights into grain boundary energy structure-property relationships by examining computed [100] disorientation axis grain boundaries in Nickel

A study of 346 computed [100] disorientation axis grain boundaries provides insight into grain boundary energy structure-property relationships. The grain boundaries come from 6 coincidence site lattice misorientations at regularly spaced disorientation angles, and cover the full range of boundary plane orientations. The computed energy has clear trends in disorientation angle and boundary plane orientation. Here, the data are used to reexamine the 'special' attribute frequently applied to low Σ coincidence site lattice misorientations and low-angle boundaries and to assess the accuracy of a face-centered cubic grain boundary energy function.

36 MATERIALS SCIENCE↗