Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic models”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
NASA NTRSDate not supplied

Twin nucleation and growth mechanism in Ni-based superalloys

While micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 C, many aspects of twin nucleation and growth remain unexplored. Kolbe mechanism for micro-twinning, based on thermally activated reordering, is probably the only concept currently widely accepted by the scientific community. We propose a qualitatively different mechanism for nucleation and growth of twins. The mechanism can be briefly described as follows. Penetration of a   precipitate by two 1/2(110) edge dislocations travelling on adjacent {111} glide planes triggers nucleation (at the interface of the precipitate and the matrix) and emission of Shockley partial of screw character on the glide plane of the edge dislocation, which entered the precipitate first (generating trailing high-energy anti-phase boundary (APB)). Propagation of this Shockley partial into the precipitate converts the APB into super intrinsic stacking fault (SISF). The recurring arrival of additional edge dislocations on the glide planes adjacent to the configuration described above leads to formation of super extrinsic stacking fault SESF, subsequent micro-twin formation, and growth of the twinned region. We demonstrate the proposed mechanism via molecular dynamics simulations.

Ni-based superalloys↗
NASA NTRSDate not supplied

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗
NASA NTRSDate not supplied

Electrochemical Recycling of CO2

Carbon dioxide (CO2) is a molecule that plays multiple roles in our universe. Currently, its impact on global warming has rightfully given CO2 a label as a harmful waste product that must be amply handled. Nonetheless, under other circumstances, and from other perspectives, CO2 could be regarded as a valuable resource. As example, in nature CO2 is an essential building block in photosynthesis ultimately yielding molecules and materials that we harvest for numerous applications ranging from energy storage and construction to medicines and food. In a similar manner, there is a growing interest to utilize atmospheric CO2 to generate renewable fuels and platform chemicals through artificial processes. These processes can replace fossil based industrial processes on earth alleviating climate impact, but they can also find uses far away from our planet. In space, CO2 is an important source of precious oxygen, however, it can also serve as feedstock for carbon-based molecules and materials that are needed to sustain societies at distant planets such as Mars. In this presentation, I will discuss the contemporary understanding of electrochemical reduction of CO2 for selective production of desired molecules. We will zoom in to the very atomic scale to gain mechanistic understanding of how the structure and composition of the cathode electrode affects product distributions. This information is used as guide in the design of new materials that convert CO2 more efficiently and more selectively. On this basis, opportunities for recycling CO2 under different conditions will be discussed, as will remaining challenges for the ultimate application towards the sustainable human living on earth and in space.

electrochemistry↗
NASA NTRSDate not supplied

Effect of Alloying Additions on Twinning in Ni-Based Superalloys

Micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 C. We use atomistic simulations to study two mechanisms of twin nucleation and growth that are characterized by qualitatively different rate limiting processes. In case of the mechanism proposed by Kolbe, the rate limiting process is diffusion-mediated atomic reshuffling. In case of the other mechanism, we proposed recently, the rate limiting process is nucleation of Shockley partial dislocation. We demonstrate the effects of alloying additions on functionality of these mechanisms.

Ni-based superalloys↗
NASA NTRS2016

Multiscale Modeling of Carbon/Phenolic Composite Thermal Protection Materials: Atomistic to Effective Properties

Next generation ablative thermal protection systems are expected to consist of 3D woven composite architectures. It is well known that composites can be tailored to achieve desired mechanical and thermal properties in various directions and thus can be made fit-for-purpose if the proper combination of constituent materials and microstructures can be realized. In the present work, the first, multiscale, atomistically-informed, computational analysis of mechanical and thermal properties of a present day - Carbon/Phenolic composite Thermal Protection System (TPS) material is conducted. Model results are compared to measured in-plane and out-of-plane mechanical and thermal properties to validate the computational approach. Results indicate that given sufficient microstructural fidelity, along with lowerscale, constituent properties derived from molecular dynamics simulations, accurate composite level (effective) thermo-elastic properties can be obtained. This suggests that next generation TPS properties can be accurately estimated via atomistically informed multiscale analysis.

Materials Science↗
NASA NTRS2007

Elastic Response and Failure Studies of Multi-Wall Carbon Nanotube Twisted Yarns

Experimental data on the stress-strain behavior of a polymer multiwall carbon nanotube (MWCNT) yarn composite are used to motivate an initial study in multi-scale modeling of strength and stiffness. Atomistic and continuum length scale modeling methods are outlined to illustrate the range of parameters required to accurately model behavior. The carbon nanotubes yarns are four-ply, twisted, and combined with an elastomer to form a single-layer, unidirectional composite. Due to this textile structure, the yarn is a complicated system of unique geometric relationships subjected to combined loads. Experimental data illustrate the local failure modes induced by static, tensile tests. Key structure-property relationships are highlighted at each length scale indicating opportunities for parametric studies to assist the selection of advantageous material development and manufacturing methods.

Gates, Thomas S.↗
NASA NTRS2003

Alloy Design Workbench-Surface Modeling Package Developed

NASA Glenn Research Center's Computational Materials Group has integrated a graphical user interface with in-house-developed surface modeling capabilities, with the goal of using computationally efficient atomistic simulations to aid the development of advanced aerospace materials, through the modeling of alloy surfaces, surface alloys, and segregation. The software is also ideal for modeling nanomaterials, since surface and interfacial effects can dominate material behavior and properties at this level. Through the combination of an accurate atomistic surface modeling methodology and an efficient computational engine, it is now possible to directly model these types of surface phenomenon and metallic nanostructures without a supercomputer. Fulfilling a High Operating Temperature Propulsion Components (HOTPC) project level-I milestone, a graphical user interface was created for a suite of quantum approximate atomistic materials modeling Fortran programs developed at Glenn. The resulting "Alloy Design Workbench-Surface Modeling Package" (ADW-SMP) is the combination of proven quantum approximate Bozzolo-Ferrante-Smith (BFS) algorithms (refs. 1 and 2) with a productivity-enhancing graphical front end. Written in the portable, platform independent Java programming language, the graphical user interface calls on extensively tested Fortran programs running in the background for the detailed computational tasks. Designed to run on desktop computers, the package has been deployed on PC, Mac, and SGI computer systems. The graphical user interface integrates two modes of computational materials exploration. One mode uses Monte Carlo simulations to determine lowest energy equilibrium configurations. The second approach is an interactive "what if" comparison of atomic configuration energies, designed to provide real-time insight into the underlying drivers of alloying processes.

Abel, Phillip B.↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Systematic improvement of redox potential calculation of Fe(III)/Fe(II) complexes using a three-layer micro-solvation model

Electrochemical transformations of metal ions in aqueous media are challenging to model accurately due to the dynamic solvation structure surrounding ions at different charge states. Predictive modeling at the atomistic scale is essential for understanding these solvation architectures but is often computationally prohibitive. In this contribution, we present a simple, fast, and accurate three-layer micro-solvation model to evaluate the redox potential of metal ions in aqueous solutions. Our model, developed and validated for Fe 3+ /Fe 2+ redox potentials, combines the DFT-based geometry optimizations of the octahedral Fe complex with two layers of explicit water molecules to capture solute–solvent interactions and an implicit solvation model to account for bulk solvent effects. This approach yields accurate predictions for Fe 3+ /Fe 2+ redox potentials in water, achieving errors of 0.02 V with ωB97X-V, 0.01 V with ωB97X-D3, 0.04 V with ωB97M-V, and 0.02 V with B3LYP-D3 functionals. We further demonstrate the generality of our model by applying it to additional metal complexes, including the challenging Fe(CN) 6 3−/4− system, where our model successfully achieves close agreement with experimental values, with an error of 0.07 V and an average error of 0.21 V for all five systems. In summary, the presented simple solvation model has broad applicability and potential for enhancing computational efficiency in redox potential predictions across various chemical and industrial processes of metal ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗
DOE OSTI2025

A Universal Model of Cation Effects in Electrocatalysis

Electrolyte cations are conventionally viewed as inert spectators in electrocatalysis. However, a wealth of observations show that catalytic rates are often highly sensitive to cation identity. Despite their prevalence, these cation effects have resisted a unified mechanistic explanation, with different physical phenomena implicated across reaction chemistries, catalyst compositions, and choice of solvent. In this perspective, we describe a general framework for understanding cation effects in electrocatalysis based on electrostatics. We argue that cations influence reaction rates by modifying the strength of the electric field present at the catalyst surface, which alters the energetics of adsorbed intermediates and transition states according to their dipole moments and polarizabilities. The magnitude of this field depends on how cations arrange at the electrode surface, controlled by their size, shape, solvation, and packing efficiency. Cations that can arrange more densely result in a steeper potential drop at the electrode surface and consequently a stronger electric field. Our model further identifies two criteria for observing cation effects: (1) the operating potential must be negative of the electrode’s potential of zero total charge, ensuring that cations accumulate at the interface, and (2) the energetics of the kinetically relevant elementary step must be field sensitive. This framework reconciles previously inconsistent trends, including why cation effects appear only for some catalysts, why reaction selectivity is sensitive to cation identity, and why activity can increase with cation size on certain metals but decrease on others. Supported by kinetic measurements, spectroscopy, and atomistic simulations, the model provides both conceptual value for building intuition about catalysis at charged interfaces and predictive value for anticipating trends for new reactions, catalysts, and electrolytes. We conclude by highlighting the importance of electric fields across electrochemical, thermochemical, and biological catalysis and propose that considering the electrostatic environment around active sites offers new opportunities for improving activity and selectivity.

cation effects↗
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2026

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗
DOE OSTI2024

Towards engineering the growth morphology of anisotropic bicrystals

Eutectic alloys, such as Pb-free Bi-Sn-Ag or Sn-Ag-Cu, are essential for their low melting points in solder joints but suffer from brittle intermetallic phases, reducing reliability. Traditional solder design relies on equilibrium phase diagrams, often neglecting the kinetic influence of crystalline anisotropy during solidification. This study hypothesizes that the energetics at the solid-solid-liquid trijunction dictate crystalline facet selection. Using the Al-Al 3 Ni alloy as a model, this project combines atomistic simulations and experimental methods to understand interfacial energies and develop predictive models. This approach aims to mitigate the negative effects of intermetallics, enhancing component reliability in electronic applications.

36 MATERIALS SCIENCE↗
NASA NTRS1969

Self-consistent many-electron theory of electron work functions and surface potential characteristics for selected metals

Electron work functions, surface potentials, and electron number density distributions and electric fields in the surface region of 26 metals were calculated from first principles within the free electron model. Calculation proceeded from an expression of the total energy as a functional of the electron number density, including exchange and correlation energies, as well as a first inhomogeneity term. The self-consistent solution was obtained via a variational procedure. Surface barriers were due principally to many-body effects; dipole barriers were small only for some alkali metals, becoming quite large for the transition metals. Surface energies were inadequately described by this model, which neglects atomistic effects. Reasonable results were obtained for electron work functions and surface potential characteristics, maximum electron densities varying by a factor of over 60.

Smith, J. R.↗
NASA NTRS2006

Computer Simulation of Fracture in Aerogels

Aerogels are of interest to the aerospace community primarily for their thermal properties, notably their low thermal conductivities. While the gels are typically fragile, recent advances in the application of conformal polymer layers to these gels has made them potentially useful as lightweight structural materials as well. In this work, we investigate the strength and fracture behavior of silica aerogels using a molecular statics-based computer simulation technique. The gels' structure is simulated via a Diffusion Limited Cluster Aggregation (DLCA) algorithm, which produces fractal structures representing experimentally observed aggregates of so-called secondary particles, themselves composed of amorphous silica primary particles an order of magnitude smaller. We have performed multi-length-scale simulations of fracture in silica aerogels, in which the interaction b e e n two secondary particles is assumed to be described by a Morse pair potential parameterized such that the potential range is much smaller than the secondary particle size. These Morse parameters are obtained by atomistic simulation of models of the experimentally-observed amorphous silica "bridges," with the fracture behavior of these bridges modeled via molecular statics using a Morse/Coulomb potential for silica. We consider the energetics of the fracture, and compare qualitative features of low-and high-density gel fracture.

Good, Brian S.↗