Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MD 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

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

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 NTRS2016

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗
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 NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗
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 NTRS1991

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of evolution of life. Using computational and modeling techniques and guided by analogies from present day systems we, seek to understand the emergence of genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. In one possible scenario, the RNA enzymatic reaction plays a key role in the emergence of the self-replicating and offers a clue to the onset of enzymatic catalysis prior to the existence of the protein biosynthetic machinery. Our ultimate goal is to propose a simple RNA segment which contains the specificity and catalytic activity of the contemporary RNA enzyme and which could emerge in a primordial chemical environment. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation have been studied. Another important aspect of this reaction is the identity of the functional groups which are essential for the acid base catalysis. The structural characteristics of the target helices, particularly a possible role of G center dot T pair, is under examination by molecular dynamics (MD) simulation technique. Modeling of the ancestral aminoacyl-tRNA synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. Assuming that the catalytic function evolved before the elements of specific recognition of a particular amino acid, we are exploring the minimal structural requirements for the catalysis of tRNA aminoacylation. The molecular modeling system SYBYL was used for this study based on the high resolution crystallographic structures of the present day tyrosyl-adenylate:tyrRS and tRNA(Gln): ATP:glnRS complexes. The trinucleotide CCA of the 3'-end tRNA is placed into the active site pocket of tyrRS, based on the scheme of interaction between tRNA(Gln) and glnRS, and upon the stereochemistry of the tyrRS:tRNA:Tyr-AMP transition state. This provides a model of the non-specific recognition of a tRNA's 3'-end by an aRS, which might be similar to that of the ancestral aRS's. In the next step, modeling of the rest of the acceptor stem of tRNA (Tyr) with tyrRS is carried out.

Rein, Robert↗
NASA NTRS1991

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of the evolution of life. Using computational and modeling techniques and guided by analogies from present day systems, we seek to understand the emergence of the genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. Modeling of the ancestral aminoacyl-tRNA-synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. The minimal structural requirements for the catalysis of tRNA aminoacylation are being explored. A formation of an aminoacyl adenylate was studied in the framework of ab initio molecular orbital theory. The role of individual residues in the vicinity of the TyrRS active site was examined, and the effect of all possible amino acids substitutions near the active site was examined. A formation of aminoacyl tRNA was studied by the molecular modeling system SYBYL with the high resolution crystallographic structures of the present day tRNA, aRS's complexes. The ultimate goal is to propose a simple RNA segment that is small enough to be build in the primordial chemical environment but maintains the specificity and catalytic activity of the contemporary RNA enzyme. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation were studied. The structural characteristics of the target helices, particularly a possible role for the G-T pair, is also studied by a molecular dynamics (MD) simulation technique.

Rein, Robert↗
NASA NTRS1991

Molecular Dynamics study of Pb overlayer on Cu(100)

Isothermal-isobaric Molecular Dynamics (MD) simulation of a submonolayer Pb film in c(2x2) ordered structure adsorbed on a Cu(100) substrate showed retention of order to high T. The Embedded Atom Method (EAM) calculated the energy of atoms of overlayer and substrate. The time-averaged squared modulus of the two dimensional structure factor for the Pb overlayer measured the order of the overlayer. The results are for increasing T only, and require verification by simulated cooling.

Karimi, M.↗
NASA NTRS2001

Thermodynamic Behavior of Nano-sized Gold Clusters on the (001) Surface

We have studied thermal expansion of the surface layers of the hexagonally reconstructed Au (001) surface using a classical Molecular Dynamics (MD) simulation technique with an Embedded Atomic Method (EAM) type many-body potential. We find that the top-most hexagonal layer contracts as temperature increases, whereas the second layer expands or contracts depending on the system size. The magnitude of expansion coefficient of the top layer is much larger than that of the other layers. The calculated thermal expansion coefficients of the top-most layer are about -4.93 x 10(exp -5)Angstroms/Kelvin for the (262 x 227)Angstrom cluster and -3.05 x 10(exp -5)Angstroms/Kelvin for (101 x 87)Angstrom cluster. The Fast Fourier Transform (FFT) image of the atomic density shows that there exists a rotated domain of the top-most hexagonal cluster with rotation angle close to 1 degree at temperature T less than 1000Kelvin. As the temperature increases this domain undergoes a surface orientational phase transition. These predictions are in good agreement with previous phenomenological theories and experimental studies.

Paik, Sun M.↗
NASA NTRS2001

Thermal Expansion and Diffusion Coefficients of Carbon Nanotube-Polymer Composites

Classical molecular dynamics (MD) simulations employing Brenner potential for intra-nanotube interactions and van der Waals forces for polymer-nanotube interface have been used to investigate thermal expansion and diffusion characteristics of carbon nanotube-polyethylene composites. Addition of carbon nanotubes to polymer matrix is found to significantly increase the glass transition temperature Tg, and thermal expansion and diffusion coefficients in the composite above Tg. The increase has been attributed to the temperature dependent increase of the excluded volume for the polymer chains, and the findings could have implications in the composite processing, coating and painting applications.

Wei, Chengyu↗
NASA NTRS1998

Observation and Modeling of Single Wall Carbon Nanotube Bend Junctions

Single wall carbon nanotube (SWNT) bends, with diameters from approx. 1.0 to 2.5 nm and bend angles from 18 deg. to 34 deg., are observed in catalytic decomposition of hydrocarbons at 600 - 1200 C. An algorithm using molecular dynamics simulation (MD) techniques is developed to model these structures that are considered to be SWNT junctions formed by topological defects (i.e. pentagon-heptagon pairs). The algorithm is used to predict the tube helicities and defect configurations for bend junctions using the observed tube diameters and bend angles. The number and arrangement of the defects at the junction interfaces are found to depend on the tube helicities and bend angle. The structural and energetic calculations using the Brenner potential show a number of stable junction configurations for each bend angle with the 34 deg. bends being more stable than the others. Tight binding calculations for local density of state (LDOS) and transmission coefficients are carried out to investigate electrical properties of the bend junctions.

Han, Jie↗
NASA NTRS1998

Nanoscale Etching and Indentation of Silicon Surfaces with Carbon Nanotubes

The possibility of nanolithography of silicon and germanium surfaces with bare carbon nanotube tips of scanning probe microscopy devices is considered with large scale classical molecular dynamics (MD) simulations employing Tersoff's reactive many-body potential for heteroatomic C/Si/Ge system. Lithography plays a key role in semiconductor manufacturing, and it is expected that future molecular and quantum electronic devices will be fabricated with nanolithographic and nanodeposition techniques. Carbon nanotubes, rolled up sheets of graphene made of carbon, are excellent candidates for use in nanolithography because they are extremely strong along axial direction and yet extremely elastic along radial direction. In the simulations, the interaction of a carbon nanotube tip with silicon surfaces is explored in two regimes. In the first scenario, the nanotubes barely touch the surface, while in the second they are pushed into the surface to make "nano holes". The first - gentle scenario mimics the nanotube-surface chemical reaction induced by the vertical mechanical manipulation of the nanotube. The second -digging - scenario intends to study the indentation profiles. The following results are reported in the two cases. In the first regime, depending on the surface impact site, two major outcomes outcomes are the selective removal of either a single surface atom or a surface dimer off the silicon surface. In the second regime, the indentation of a silicon substrate by the nanotube is observed. Upon the nanotube withdrawal, several surface silicon atoms are adsorbed at the tip of the nanotube causing significant rearrangements of atoms comprising the surface layer of the silicon substrate. The results are explained in terms of relative strength of C-C, C-Si, and Si-Si bonds. The proposed method is very robust and does not require applied voltage between the nanotube tips and the surface. The implications of the reported controllable etching and hole-creating for nanolithography on silicon are discussed in detail.

Dzegilenko, Fedor N.↗
NASA NTRS2002

Analysis of Carbon Nanotube Pull-out from a Polymer Matrix

Molecular dynamics (MD) simulations of carbon nanotube (NT) pull-out from a polymer matrix are carried out. As the NT pull-out develops in the simulation, variations in the displacement and velocities of the NT are monitored. The existence of a carbon-ring-based period in NT sliding during pull-out is identified. Linear trends in the NT velocity-force relation are observed and used to estimate an effective viscosity coefficient for interfacial sliding at the NT/polymer interface. As a result, the entire process of NT pull-out is characterized by an interfacial friction model that is based on a critical pull-out force, and an analog of Newton's friction law used to describe the NT/polymer interfacial sliding.

Frankland, S. J. V.↗
NASA NTRS2003

Theory of Transport of Long Polymer Molecules through Carbon Nanotube Channels

A theory of transport of long chain polymer molecules through carbon nanotube (CNT) channels is developed using Fokker-Planck equation and direct molecular dynamics (MD) simulations. The mean transport or translocation time tau is found to depend on the chemical potential energy, entropy and diffusion coefficient. A power law dependence tau approx. N(sup 2)is found where N is number of monomers in a molecule. For 10(exp 5)-unit long polyethylene molecules, tau is estimated to be approx. 1micro-s. The diffusion coefficient of long polymer molecules inside CNTs, like that of short ones, are found to be few orders of magnitude larger than in ordinary silicate based zeolite systems.

Wei, Chenyu↗
NASA NTRS2003

A Dynamic Pathway for Stone-Wales Bond Rotation on Carbon Nanotubes through Diamond-Like Bonds

A new lower energy barrier with a two-step pathway of Stone-Wales (SW) ,ond rotation on carbon nanotubes (CNTs) is found through molecular dynamics (MD) simulations of CNTs under tension. The first step involves going over to a stable sp3-like metastable configuration with half rotated and partially tilted C-C bond. The second step involves going over to the fully rotated C-C bond with the formation of a SW defect in the nanotube. The energy barrier for this two-step dynamic pathway is significantly lower than the previously known static barrier for in-plane rotation of the C-C bond on a tensile strained (> 4%) CNT.

Wei, Chen-Yu↗
NASA NTRS2007

Interactions of the SAP Domain of Human Ku70 with DNA Substrate: A Molecular Dynamics Study

NASA is developing a systems biology approach to improve the assessment of health risks associated with space radiation. The primary toxic and mutagenic lesion following radiation exposure is the DNA double strand break (DSB), thus a model incorporating proteins and pathways important in response and repair of this lesion is critical. One key protein heterodimer for systems models of radiation effects is the Ku70/80 complex. The Ku70/80 complex is important in the initial binding of DSB ends following DNA damage, and is a component of nonhomologous end joining repair, the primary pathway for DSB repair in mammalian cells. The SAP domain of Ku70 (residues 556-609), contains an a helix-extended strand-helix motif and similar motifs have been found in other nucleic acid-binding proteins critical for DNA repair. However, the exact mechanism of damage recognition and substrate specificity for the Ku heterodimer remains unclear in part due to the absence of a high-resolution structure of the SAP/DNA complex. We performed a series of molecular dynamics (MD) simulations on a system with the SAP domain of Ku70 and a 10 base pairs DNA duplex. Large-scale conformational changes were observed and some putative binding modes were suggested based on energetic analysis. These modes are consistent with previous experimental investigations. In addition, the results indicate that cooperation of SAP with other domains of Ku70/80 is necessary to explain the high affinity of binding as observed in experiments.

Hu, Shaowen↗
NASA NTRS2009

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [111] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.

Yamakov, V.↗