SEARCH · Engineering Papers

Results for “molecular dynamic 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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗

NASA NTRSDate not supplied

Molecular dynamics simulation of effects of solutes on dislocation propagation in Ni-based superalloys

Ni-based superalloys are used in the hot sections of jet turbine engines because of their high strength, stability and resistance to oxidation. The properties of these alloys can be further optimized by adding alloying elements. Therefore, a fundamental knowledge on the effect of different elements on properties of Ni-based superalloys is required. Molecular dynamics simulation could shed light here but its application is hindered by the absence of reliable and computationally cheap semi-empirical potential of the interatomic interaction for 4 and more element alloys. We will present a new Ni-Al-Cr-Nb Finnis-Sinclair (FS) potential specially designed to simulate the dislocation propagation from to  phase. In order to construct this potential, we designed a special algorithm to incorporate the data on element partitioning in the potential development procedure. For example, it is known from experiment, that Cr is mostly present in the gamma phase. Figure 1 shows a snapshot obtained after equilibration of the model of the Ni68Al17Cr15 alloy at T=1000 K using the hybrid Monte-Carlo (MC)/molecular dynamics (MD) simulation with the developed semi-empirical potential. One can clearly see that the Cr partitioning is in agreement with the experimental data. We will discuss the developed algorithm to incorporate the solute partition data in details. Using the developed semi-empirical potential, we first investigated the effect of anti-site defects in the  phase on the single dislocation propagation. It was found that the dislocation velocity increases with the increasing of the anti-site defect concentration. This effect was attributed to smaller number of Al-Al pairs forming during the dislocation migration in the presence of the anti-site defects. Next, we investigated the effect of Nb on the dislocation pair propagation in the Kolbe mechanism. It was found that the addition of Nb leads to considerable decrease in the dislocation propagation rate. This is in agreement with the experimental data on the effect of Nb on the creep resistance of the Ni-based superalloys. We will discuss the origin of this effect.

Mikhail I. Mendelev↗

NASA NTRS1997

Comparison of Chain Conformation of Poly(vinyl alcohol) in Solutions and Melts from Quantum Chemistry Based Molecular Dynamics Simulations

Confirmations of 2,4-dihydroxypentane (DHP), a model molecule for poly(vinyl alcohol), have been studied by quantum chemistry (QC) calculations and molecular dynamics (MD) simulations. QC calculations at the 6-311G MP2 level show the meso tt conformer to be lowest in energy followed by the racemic tg, due to intramolecular hydrogen bond between the hydroxy groups. The Dreiding force field has been modified to reproduce the QC conformer energies for DHP. MD simulations using this force field have been carried out for DHP molecules in the gas phase, melt, and CHCl3 and water solutions. Extensive intramolecular hydrogen bonding is observed for the gas phase and CHCl3 solution, but not for the melt or aqueous solution, Such a condensed phase effect due to intermolecular interactions results in a drastic change in chain conformations, in agreement with experiments.

Jaffe, Richard↗

NASA NTRS1997

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗

NASA NTRS2018

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗

NASA NTRS2023

Molecular Dynamics Simulations of Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

The unique properties of shape memory alloys (SMAs) arise from a reversible martensitic transformation. The nucleation and migration of austenite-martensite interfaces are the key to understanding the SMA properties. Molecular dynamics (MD) simulations can provide important atomic-scale information about these aspects of the transformation, but their time scales prevent the interface formation under near-equilibrium conditions relevant to experiment. We present a new MD methodology which allows for the natural formation of energetically preferred austenite-martensite interfaces under near-equilibrium conditions. Our simulation demonstrates that the interfaces in NiTi are semi-coherent, composed of a series of terrace planes and structural disconnections, and they migrate rapidly through single crystals with only a small thermodynamic driving force. In bi-crystals and polycrystals, the migration of these same interfaces is significantly impeded by grain boundaries and stored elastic energy. This behavior can result in SMA hysteresis via several mechanisms associated with nucleation and non-elastic strain accommodation.

Gabriel Plummer↗

NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗

NASA NTRS2019

Molecular Dynamics Simulations of Silicon Carbide, Boron Nitride and Silicon for Ceramic Matrix Composite Applications

A comprehensive computational molecular dynamics study is presented for crystalline α-SiC (6H, 4H, and 2H SiC), β-SiC (3C SiC), layered boron nitride, amorphous boron nitride and silicon, the constituent materials for high-temperature SiC/SiC compositions. Large-scale Atomic/Molecular Parallel Simulator software package was used. The Tersoff Potential force field was utilized to evaluate their mechanical characteristics of most of the materials, and the Reax force field was used to model silicon when the Tersoff Potential did not provide accurate results. Their mechanical behaviors were evaluated at a strain rate of 10(exp 7)/s and the results agree with the experimental data in the literature. The results are foundational for linking constituent behavior to composite performance, particularly when test data is unavailable or suspect.

Aluko, Olanrewaju↗

NASA NTRSDate not supplied

Scalable Variable Charge Molecular Dynamics Simulations of Metal-Oxide Systems

Interfaces between metals and oxides are an important feature in many technologically relevant materials, e.g., oxidation of metal surfaces, oxide-dispersion strengthened (ODS) alloys, dielectric components, and thermal barrier coatings among others. Experimental studies of such interfaces are challenging since the majority are buried within the bulk, making computational modeling an attractive alternative. Molecular dynamics (MD) simulations operate at the length scales relevant to many interface-mediated mechanisms, but the requisite interatomic potentials for metal-oxide systems require computationally expensive variable charge schemes to account for the disparate bonding types, thus often limiting their effectiveness. Here we introduce several improvements to the charge transfer interatomic potential (CTIP) model which enable greater computational efficiency for large scale MD simulations. Then, using a new CTIP parametrization for the Ni-O system, we demonstrate its capabilities to capture critical atomic scale mechanisms associated with metal-oxide interfaces. Long time scale simulations (>10 ns) are used to investigate high temperature oxidation and oxide precipitation from the melt, and large length scale simulations (> 1 million atoms) are used to study the interaction of dislocations with oxide particles. We have implemented the new CTIP model in the widely used, open-source MD code LAMMPS.

Gabriel Plummer↗

NASA NTRS2024

Molecular Dynamics Simulations of Microstructural Effects on Austenite-Martensite Interfaces in NiTi

Formation and migration of austenite-martensite interfaces plays the key role in reversible martensitic transformations of shape memory alloys (SMAs). How these interfaces interact with the SMA microstructure is a primary determining factor in important functional properties such as hysteresis and transformation span. As such, successful microstructural engineering of SMAs requires in-depth knowledge of interface behavior. The rapid nature of martensitic transformations makes experimental observations of moving austenite-martensite interfaces challenging. Molecular dynamics (MD) simulation is a unique tool which can probe the atomic-scale details of austenite-martensite interfaces as they migrate through different microstructures. However, in focusing on the entire transformation process, including the nucleation of new phases, MD studies are usually performed so far from equilibrium that their relevance to experiment is questionable. Here, we demonstrate new MD techniques to generate energetically preferred austenite-martensite interfaces in NiTi under near-equilibrium conditions. The interfaces are semi-coherent, exhibiting a series of structural disconnections, and they can migrate rapidly through single crystals under only small thermodynamic driving forces. In contrast, when interfaces migrate in polycrystals, their motion is impeded by thermoelastic effects as well as changes in orientation relationships at grain boundaries. Microstructures which accumulate large amounts of elastic energy tend to release some fraction through irreversible, hysteresis-inducing mechanisms. We demonstrate that engineering microstructures with less constraints is a viable strategy to produce SMAs with reduced hysteresis and transformation span. Similar thermoelastic and hysteresis-inducing mechanisms also arise when austenite-martensite interfaces encounter precipitates and can be controlled by tuning characteristics of the precipitate distribution.

Gabriel Plummer↗

NASA NTRSDate not supplied

Molecular Dynamics Simulations of Microstructural Effects on Austenite-Martensite Interfaces in NiTi

Gabriel Plummer↗

NASA NTRSDate not supplied

Molecular Dynamics Simulations of Microstructural Effects on Austenite-Martensite Interfaces in NiTi

The formation and migration of austenite-martensite interfaces plays the key role in reversible martensitic transformations of shape memory alloys (SMAs). How these interfaces interact with the SMA microstructure is a primary determining factor in important functional properties such as hysteresis and transformation span. As such, successful microstructural engineering of SMAs requires in-depth knowledge of interface behavior. The rapid nature of martensitic transformations makes experimental observations of moving austenite-martensite interfaces challenging. Molecular dynamics (MD) simulation is a unique tool which can probe the atomic-scale details of austenite-martensite interfaces as they migrate through different microstructures. However, in focusing on the entire transformation process, including the nucleation of new phases, MD studies are usually performed so far from equilibrium that their relevance to experiment is questionable. Here, we demonstrate new MD techniques to generate energetically preferred austenite-martensite interfaces in NiTi under near-equilibrium conditions. The interfaces are semi-coherent, exhibiting a series of structural disconnections, and they can migrate rapidly through single crystals under only small thermodynamic driving forces. In contrast, when interfaces migrate in polycrystals, their motion is impeded by thermoelastic effects as well as changes in orientation relationships at grain boundaries. Microstructures which accumulate large amounts of elastic energy tend to release some fraction through irreversible, hysteresis-inducing mechanisms. We demonstrate that engineering microstructures with less constraints is a viable strategy to produce SMAs with reduced hysteresis and transformation span. Similar thermoelastic and hysteresis-inducing mechanisms also arise when austenite-martensite interfaces encounter precipitates and can be controlled by tuning characteristics of the precipitate distribution.

Gabriel Plummer↗

NASA NTRSDate not supplied

Characterizing Porous and Nonporous Phenolic Resins from Molecular Dynamics Simulations

Phenolic resins are an important component of many ablative heat shield materials, which protect spacecrafts from the extreme temperatures reached during atmospheric entry. Examples include the high-density Heritage Carbon Phenolic (HCP) used in the Pioneer-Venus and Galileo missions, as well as the low-density Phenolic Impregnated Carbon Ablator (PICA) used in the Mars Science Laboratory and Mars 2020 missions. Additionally, recent developments within NASA have produced the mid-density Heatshield for Extreme Entry Environment Technology (HEEET) and its derivative 3D Woven Mid-Density Carbon Phenolic (3MDCP). Unlike the nonporous phenolic in HCP, PICA and HEEET/3MDCP are fabricated by infusing preforms with diluted phenolic formulations to obtain a lower density porous matrix. Despite the importance of the phenolic phase to the material response during entry, the variation in properties of porous and nonporous phenolic is not well understood. Here, we present an investigation of porous and nonporous phenolic resins using molecular dynamics (MD) simulations. Resin cure is mimicked in the simulations through the inclusion of representative reaction templates to generate accurate models of the complex crosslinked structures. To create porous models, explicit solvent molecules are included during the cure simulations. We observe nanoscale separation of the phenolic and solvent phases, which results in significant differences in the final structures of porous and nonporous models. In addition to a quantitative assessment of the network structure and porosity, we elucidate the effects of the phenolic formulation on the final material properties. These results are compared with experimental data as appropriate.

phenolic↗

NASA NTRS1998

Quantum Molecular Dynamics Simulations of Nanotube Tip Assisted Reactions

In this report we detail the development and application of an efficient quantum molecular dynamics computational algorithm and its application to the nanotube-tip assisted reactions on silicon and diamond surfaces. The calculations shed interesting insights into the microscopic picture of tip surface interactions.

Menon, Madhu↗

NASA NTRS2018

Aladyn - Adaptive Neural Network Molecular Dynamics Simulation Code: Computational Materials Mini-Application

This report provides an overview and commands description of the Computational Materials mini-application, Aladyn. Aladyn is a simple molecular dynamics code written in FORTRAN 2008, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of quantum mechanics-based energy calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. The Aladyn code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing many core architectures based on graphic processing units (GPUs). The effort is related to the High Performance Computing Incubator (HPCI) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗

NASA NTRS2019

AladynPi – Adaptive Neural Network Molecular Dynamics Simulation Code with Physically Informed Potential: Computational Materials Mini-Application

This report provides an overview and description of commands used in the Computational Materials mini-application, AladynPi. AladynPi is an extension of a previously released mini-application, Aladyn (https://github.com/nasa/aladyn; Yamakov, V.I., and Glaessgen, E.H., NASA/TM-2018-220104). Aladyn and AladynPi are basic molecular dynamics codes written in FORTRAN 2003, which are designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory method. An input for the ANN is a set of structure coefficients, characterizing the local atomic environment of each atom, for which the atomic energy is obtained in the ANN inference process. In Aladyn, the ANN gives directly the energy of interatomic interactions. In AladynPi, the ANN gives optimized parameters for a predefined empirical function, known as bond-order-potential (BOP). The parameterized BOP function is then used to calculate the energy. AladynPi code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing manycore architectures based on graphic processing units (GPUs). The effort is supported by the High Performance Computing incubator (HPCi) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗

NASA NTRS2000

Molecular Dynamics Simulations of the Thermal Conductivity of Single-Wall Carbon Nanotubes

Carbon nanotubes (CNT) have very attractive electronic, mechanical. and thermal properties. Recently, measurements of thermal conductivity in single wall CNT mats showed estimated thermal conductivity magnitudes ranging from 17.5 to 58 W/cm-K at room temperature. which are better than bulk graphite. The cylinderical symmetry of CNT leads to large thermal conductivity along the tube axis, additionally, unlike graphite. CNTs can be made into ropes that can be used as heat conducting pipes for nanoscale applications. The thermal conductivity of several single wall carbon nanotubes has been calculated over temperature range from l00 K to 600 K using non-equilibrium molecular dynamics using Tersoff-Brenner potential for C-C interactions. Thermal conductivity of single wall CNTs shows a peaking behavior as a function of temperature. Dependence of the peak position on the chirality and radius of the tube will be discussed and explained in this presentation.

Osman, M.↗

NASA NTRS2021

Molecular dynamics simulations of ultrafast radiation induced melting at metal–semiconductor interfaces

Metal–semiconductor contacts in silicon carbide (SiC) diodes endure damages at the interface when exposed to harsh radiation environments. Due to the rapid rise in temperature and ultrafast cooling that follows the radiation impact, the structural properties of the materials can be altered through melting, recrystallization, and amorphization. A detailed understanding of the material failure modes at the interface is lacking, specifically at the nanoscale. We use molecular simulations to investigate the ultrafast melting at tungsten (W)–SiC interfaces following radiation damage and apply deep learning techniques to track the transient evolution of the local molecular structures. We show that W near the radiation track undergoes melting and, eventually, most of it recrystallizes with a noticeable degree of undercooling, while SiC is rendered permanently amorphous. The observation of local undercooling in W films is important as it can affect the device performance even before the bulk melting temperature of the material is reached. We also show that at high temperatures, the interface undergoes a fracture-like failure. The results presented here are significant in understating the different failure modes of SiC diode materials.

Ashwin Ravichandran↗