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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ab initio Geochemistry of Hydrous Phases (Final Report)

Hydrous phases are among the most important Earth components. They are important for a broad suite of Earth processes, including the origin of life. From both thermodynamic and structural perspectives, however, they represent some of the most complex naturally occurring materials: their bonding often includes a combination of covalent solid, weak ionic, van der Waals, and hydrogen bonding, all within large unit cells. Most are solid solutions, and many are prone to variations in layer packing. Thermodynamic modeling of these materials is fundamental for understanding present and past natural processes, including those that shaped—and continue to shape—the structure and evolution of our planet. Yet, the thermodynamic properties of these materials at appropriate conditions are challenging to measure. To make significant progress and attain a deep understanding of these materials requires an atomistic theoretical approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unraveling the Ion Adsorption Kinetics in Microporous Carbon Electrodes: A Multiscale Quantum-Continuum Simulation and Experimental Approach

With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies

36 MATERIALS SCIENCE↗
NASA NTRS2004

Atomistic Modeling of Pd Site Preference in NiTi

An analysis of the site subsitution behavior of Pd in NiTi was performed using the BFS method for alloys. Through a combination of Monte Carlo simulations and detailed atom-by-atom energetic analyses of various computational cells, representing compositions of NiTi with up to 10 at% Pd, a detailed understanding of site occupancy of Pd in NiTi was revealed. Pd subsituted at the expense of Ni in a NiTi alloy will prefer the Ni-sites. Pd subsituted at the expense of Ti shows a very weak preference for Ti-sites that diminishes as the amount of Pd in the alloy increases and as the temperature increases.

Bozzolo, Guillermo↗
DOE OSTI2021

Computational study of the effect of core–skin structure on the mechanical properties of carbon nanofibers

The effect of the core–skin structure on the mechanical properties of carbon nanofibers is investigated in large-scale molecular dynamics simulations of tensile deformation of carbon nanofibers with the core–skin and homogeneous structures. Contrary to an established notion of the deleterious effect of the skin layer on the strength of carbon fibers, the presence of a high-quality skin layer is found to increase both the Young’s modulus and tensile strength of the nanofiber. A detailed analysis of the fracture process indicates that the nanofiber strengthening is related to the ability of skin layer to suppress crack nucleation at the core–skin interface. Furthermore, the computational predictions suggest that the design of new approaches to carbon fiber manufacturing and processing leading to the generation of a high-quality skin layer while avoiding the introduction of structural defects at the core–skin interface may yield a significant enhancement of the mechanical properties of carbon fibers.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Multiscale Modeling, Simulation and Visualization and Their Potential for Future Aerospace Systems

This document contains the proceedings of the Training Workshop on Multiscale Modeling, Simulation and Visualization and Their Potential for Future Aerospace Systems held at NASA Langley Research Center, Hampton, Virginia, March 5 - 6, 2002. The workshop was jointly sponsored by Old Dominion University's Center for Advanced Engineering Environments and NASA. Workshop attendees were from NASA, other government agencies, industry, and universities. The objectives of the workshop were to give overviews of the diverse activities in hierarchical approach to material modeling from continuum to atomistics; applications of multiscale modeling to advanced and improved material synthesis; defects, dislocations, and material deformation; fracture and friction; thin-film growth; characterization at nano and micro scales; and, verification and validation of numerical simulations, and to identify their potential for future aerospace systems.

Noor, Ahmed K.↗
DOE OSTI2024

Exploring dielectric properties in atomistic models of amorphous boron nitride

Abstract We report a theoretical study of dielectric properties of models of amorphous Boron Nitride, using interatomic potentials generated by machine learning. We first perform first-principles simulations on small (about 100 atoms in the periodic cell) sample sizes to explore the emergence of mid-gap states and its correlation with structural features. Next, by using a simplified tight-binding electronic model, we analyse the dielectric functions for complex three dimensional models (containing about 10.000 atoms) embedding varying concentrations of sp 1 , sp 2 and sp 3 bonds between B and N atoms. Within the limits of these methodologies, the resulting value of the zero-frequency dielectric constant is shown to be influenced by the population density of such mid-gap states and their localization characteristics. We observe nontrivial correlations between the structure-induced electronic fluctuations and the resulting dielectric constant values. Our findings are however just a first step in the quest of accessing fully accurate dielectric properties of as-grown amorphous BN of relevance for interconnect technologies and beyond.

Materials Science↗
DOE OSTI2024

Enhanced mobility of dislocation network nodes and its effect on dislocation multiplication and strain hardening

Understanding plastic deformation of crystals in terms of the fundamental physics of dislocations has remained a grand challenge in materials science for decades. To overcome this, the Discrete Dislocation Dynamics (DDD) method has been developed, but its lack of atomistic resolution leaves open the possibility that certain key mechanisms may be overlooked. Here, by comparing large-scale Molecular Dynamics (MD) with DDD simulations performed under identical conditions we uncover significant discrepancies in the predicted strength and microstructure evolution in BCC crystals under high-strain rate conditions. These are traced to unexpected behaviors of dislocation network nodes forming at dislocation intersections, that can move in ways not previously anticipated as revealed by MD. Once these newfound freedoms of nodal motion are incorporated, DDD simulations begin to closely match plastic evolution observed in MD. This additional mechanism of motion whereby non-screw dislocations can change their glide plane profoundly affects fundamental processes of dislocation multiplication, recovery and storage that define strength of metals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A molecular dynamics study of domain switching in BiFeO 3 nanofilm under DC electric field

The performance of ferroelectric materials is deeply rooted in the domain switching process triggered by an applied electric field. To understand the effect of controlling factors on the domain switching characteristics, nanoscale studies are essential. In this work, an atomistic-scale computational framework is developed to study the domain switching process of bismuth ferrite (BiFeO 3 ) nanofilms. Through a series of molecular dynamics simulations, the effects of duration and magnitude of an electric field, domain size and nanofilm thickness on the domain switching behavior are investigated. The evolution of domains under different magnitudes and durations of the applied electric field is studied, and the results of unswitched and switched cases are summarized in a switching map. Analysis of various nanofilms reveals that by increasing the film thickness, the switching time is increased. Additionally, the nanofilms with different domain sizes and configurations are examined, and the minimum domain size for complete reversal switching is obtained. Our results are helpful to predict and tune the switching process in ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Formation and coarsening of epitaxially-supported metal nanoclusters

This mini-review describes developments over the last ~30 years in characterizing the nucleation & growth of epitaxially-supported metal nanoclusters (NCs) or islands during vapor deposition, as well as their post-deposition coarsening. A beyond-mean-field treatment for homogeneous nucleation & growth corrects the deficiencies of traditional treatments in describing, e.g., the island size distribution, but also necessitates consideration of the spatial distribution of islands and their capture zones. We discuss advances in modeling capabilities, including those based upon on an ab-initio level treatment of periphery diffusion kinetics, for description of the non-equilibrium growth shapes of these NCs, focusing on 2D NCs. For post-deposition coarsening of arrays of NCs, there is generally a competition between Ostwald Ripening (OR) and Smoluchowski Ripening (SR). SR is also known as Particle Migration & Coalescence. For 2D NCs in homoepitaxial systems, conventional OR is observed on pristine fcc(111) surfaces, dramatically enhanced OR in the presence of even trace amounts of chalcogens for Cu(111) and Ag(111), and anomalous OR on anisotropic fcc(110) surfaces. The unexpected discovery of SR for fcc(100) homoepitaxial systems prompted extensive analysis of the underlying diffusivities of 2D NCs as a function of size, as well as of NC coalescence dynamics. A comprehensive understanding of these processes is now available. Self-assembly of 3D NCs during deposition, issues related to heterogeneous nucleation, directed assembly, NC growth structure selection, and coarsening are addressed. Finally, for SR of 3D epitaxial NCs, recent insights into the size-dependence of diffusivity are described.

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

Fundamental Investigations into the Properties and Performance of Advanced Materials

Intelligent engineering designs require an accurate understanding of material behavior, since any uncertainties or gaps in knowledge must be counterbalanced with heightened factors of safety, leading to overdesign. Therefore, building better structures and pushing the performance of new components requires an improved understanding of the thermomechanical response of advanced materials under service conditions. This dissertation provides fundamental investigations of several advanced materials: thermoset polymers, a common matrix material for fiber-reinforced composites and nanocomposites; aluminum alloy 7075-T6 (AA7075-T6), a high-performance aerospace material; and ceramic matrix composites (CMCs), an advanced composite for extreme-temperature applications. To understand matrix interactions with various interfaces and nanoinclusions at their fundamental scale, the properties of thermoset polymers are studied at the atomistic scale. An improved proximity-based molecular dynamics (MD) technique for modeling the crosslinking of thermoset polymers is carefully established, enabling realistic curing simulations through its ability to dynamically and probabilistically perform complex topology transformations. The proximity-based MD curing methodology is then used to explore damage initiation and the local anisotropic evolution of mechanical properties in thermoset polymers under uniaxial tension with an emphasis on changes in stiffness through a series of tensile loading, unloading, and reloading experiments. Aluminum alloys in aerospace applications often require a fatigue life of over 109 cycles, which is well over ii the number of cycles that can be practically tested using conventional fatigue testing equipment. In order to study these high-life regimes, a detailed ultrasonic cycle fatigue study is presented for AA7075-T6 under fully reversed tension-compression loading. The geometric sensitivity, frequency effects, size effects, surface roughness effects, and the corresponding failure mechanisms for ultrasonic fatigue across different fatigue regimes are investigated. Finally, because CMCs are utilized in extreme environments, oxidation plays an important role in their degradation. A multiphysics modeling methodology is thus developed to address the complex coupling between oxidation, mechanical stress, and oxygen diffusion in heterogeneous carbon fiber-reinforced CMC microstructures.

Schichtel, Jacob J.↗
DOE OSTI2023

Voltage-Dependent First-Principles Simulation of Insertion of Chloride Ions into Al/Al 2 O 3 Interfaces Using the Quantum Continuum Approximation

Experiments have shown that pitting corrosion can develop in aluminum surfaces at potentials > − 0.5 V relative to the standard hydrogen electrode (SHE). Until recently, the onset of pitting corrosion in aluminum has not been rigorously explored at an atomistic scale because of the difficulty of incorporating a voltage into density functional theory (DFT) calculations. We introduce the Quantum Continuum Approximation (QCA) which self-consistently couples explicit DFT calculations of the metal-insulator and insulator-solution interfaces to continuum Poisson-Boltzmann electrostatic distributions describing the bulk of the insulating region. By decreasing the number of atoms necessary to explicitly simulate with DFT by an order of magnitude, QCA makes the first-principles prediction of the voltage of realistic electrochemical interfaces feasible. After developing this technique, we apply QCA to predict the formation energy of chloride atoms inserting into oxygen vacancies in Al(111)/α-Al 2 O 3 (0001) interfaces as a function of applied voltage. We predict that chloride insertion is only favorable in systems with a grain boundary in the Al 2 O 3 for voltages > − 0.2 V (SHE). Here our results roughly agree with the experimentally demonstrated onset of corrosion, demonstrating QCA's utility in modeling realistic electrochemical systems at reasonable computational cost.

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

Accurate and scalable graph neural network force field and molecular dynamics with direct force architecture

Abstract Recently, machine learning (ML) has been used to address the computational cost that has been limiting ab initio molecular dynamics (AIMD). Here, we present GNNFF, a graph neural network framework to directly predict atomic forces from automatically extracted features of the local atomic environment that are translationally-invariant, but rotationally-covariant to the coordinate of the atoms. We demonstrate that GNNFF not only achieves high performance in terms of force prediction accuracy and computational speed on various materials systems, but also accurately predicts the forces of a large MD system after being trained on forces obtained from a smaller system. Finally, we use our framework to perform an MD simulation of Li 7 P 3 S 11 , a superionic conductor, and show that resulting Li diffusion coefficient is within 14% of that obtained directly from AIMD. The high performance exhibited by GNNFF can be easily generalized to study atomistic level dynamics of other material systems.

Chemistry↗
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Chemical order transitions within extended interfacial segregation zones in NbMoTaW

Interfacial segregation and chemical short-range ordering influence the behavior of grain boundaries in complex concentrated alloys. In this study, we use atomistic modeling of a NbMoTaW refractory complex concentrated alloy to provide insight into the interplay between these two phenomena. Hybrid Monte Carlo and molecular dynamics simulations are performed on columnar grain models to identify equilibrium grain boundary structures. Here our results reveal extended near-boundary segregation zones that are much larger than traditional segregation regions, which also exhibit chemical patterning that bridges the interfacial and grain interior regions. Furthermore, structural transitions pertaining to an A2-to-B2 transformation are observed within these extended segregation zones. Both grain size and temperature are found to significantly alter the widths of these regions. An analysis of chemical short-range order indicates that not all pairwise elemental interactions are affected by the presence of a grain boundary equally, as only a subset of elemental clustering types are more likely to reside near certain boundaries. The results emphasize the increased chemical complexity that is associated with near-boundary segregation zones and demonstrate the unique nature of interfacial segregation in complex concentrated alloys.

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

C-C Coupling Mechanism on Cu(100) A Molecular Dynamics Study at 298K

The electrochemical reduction of carbon dioxide (CO2) into valuable fuels such as C1 (syngas, methane) and C2 (ethylene, ethanol) products is a key strategy for achieving a carbon-neutral economy. Computational studies of C-C coupling, a critical step in CO2 reduction, are essential for designing more efficient catalysts. However, simulating these processes under realistic electrochemical conditions, including temperature and solvent effects, is computationally demanding. In this work, we develop a machine learning-based atomistic potential to study CO2 reduction on Cu(100) surfaces, accounting for temperature and explicit water solvent effects. We compute thermodynamic free energies of the possible C-C coupling pathways, CO*+CO*→OCCO*, CO*+CHO*→OCCHO*, CO*+COH*→OCCOH*, CHO*+COH*→OHCCOH*, COH*-COH*→HOCCOH*, and CHO*-CHO*→OHCCHO*. Our results quantify the thermodynamic tendencies of these reactions and reveal that, in addition to the well-established CO* + CO* → OCCO* pathway, CHO* is a critical intermediate in the formation of C2 products on Cu(100). Furthermore, we demonstrate that the machine learning approach offers a cost-efficient framework for studying CO2 reduction on diverse catalysts under realistic electrochemical conditions.

machine learning↗
DOE OSTI2021

Atomistic Assessment of Solute-Solute Interactions during Grain Boundary Segregation

Grain boundary solute segregation is becoming increasingly common as a means of stabilizing nanocrystalline alloys. Thermodynamic models for grain boundary segregation have recently revealed the need for spectral information, i.e., the full distribution of environments available at the grain boundary during segregation, in order to capture the essential physics of the problem for complex systems like nanocrystalline materials. However, there has been only one proposed method of extending spectral segregation models beyond the dilute limit, and it is based on simple, fitted parameters that are not atomistically informed. In this work, we present a physically motived atomistic method to measure the full distribution of solute-solute interaction energies at the grain boundaries in a polycrystalline environment. We then cast the results into a simple thermodynamic model, analyze the Al(Mg) system as a case study, and demonstrate strong agreement with physically rigorous hybrid Monte Carlo/molecular statics simulations. This approach provides a means of rapidly measuring key interactions for non-dilute grain boundary segregation for any system with an interatomic potential.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring structural transitions at grain boundaries in Nb using a generalized embedded atom interatomic potential

The advancement in experimental techniques, like the atom probe tomography and high resolution electron microscopy, is fueling interest in studying structural transformations of grain boundaries in metal and alloys to uncover correlations between mechanical properties and solute or impurity segregation to grain boundaries. Atomistic modeling is an important tool that can pinpoint the intricate dynamics of grain boundary phase transitions, but the lack of accurate interatomic potentials needed to simulate the complex dynamics of grain boundary structural transitions and identify different metastable phases has been the bottleneck. To this end, we use niobium as a model body centered cubic (BCC) metal and develop an interatomic potential to study grain boundary phase transitions. The potential for Nb is based on a generalization of the embedded atomic method potential and has sufficient flexibility to learn complex energy landscapes using a small set of training structures. We systematically test and validate the using data from ab initio density functional theory calculations and experiments. Using this potential, we calculate energies of multiple symmetric-tilt grain boundaries spanning a wide range of misorientation angles. Additionally, we explore different metastable structures of the Σ 27(552) [$1\overline{1}0$] grain boundary and use molecular dynamic simulations to study the coexistence of metastable phases and grain boundary transitions at finite temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Automated calculation and convergence of defect transport tensors

Defect diffusion is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end-user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewide sampling completeness, to build atomistic kinetic Monte Carlo models on a state-space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for drift and diffusion coefficients, providing a convergence metric by calculating the Kullback–Leibler divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and Hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

36 MATERIALS SCIENCE↗