Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Ab Initio Simulations of Temperature Dependent Phase Stability and Martensitic Transitions in NiTi

For NiTi based alloys, the shape memory effect is governed by a transition from a low-temperature martensite phase to a high-temperature austenite phase. Despite considerable experimental and computational work, basic questions regarding the stability of the phases and the martensitic phase transition remain unclear even for the simple case of binary, equiatomic NiTi. We perform ab initio molecular dynamics simulations to describe the temperature-dependent behavior of NiTi and resolve several of these outstanding issues. Structural correlation functions and finite temperature phonon spectra are evaluated to determine phase stability. In particular, we show that finite temperature, entropic effects stabilize the experimentally observed martensite (B19') and austenite (B2) phases while destabilizing the theoretically predicted (B33) phase. Free energy computations based on ab initio thermodynamic integration confirm these results and permit estimates of the transition temperature between the phases. In addition to the martensitic phase transition, we predict a new transition between the B33 and B19' phases. The role of defects in suppressing these phase transformations is discussed.

nickel titanium↗
NASA NTRS2016

Yield and Failure Behavior Investigated for Cross-Linked Phenolic Resins Using Molecular Dynamics

Molecular dynamics simulations were conducted to fundamentally evaluate the yield and failure behavior of cross-linked phenolic resins at temperatures below the glass transition. Yield stress was investigated at various temperatures, strain rates, and degrees of cross-linking. The onset of non-linear behavior in the cross-linked phenolic structures was caused by localized irreversible molecular rearrangements through the rotation of methylene linkers followed by the formation or annihilation of neighboring hydrogen bonds. The yield stress results, with respect to temperature and strain rate, could be fit by existing models used to describe yield behavior of amorphous glasses. The degree of cross-linking only indirectly influences the maximum yield stress through its influence on glass transition temperature (Tg), however there is a strong relationship between the degree of cross-linking and the failure mechanism. Low cross-linked samples were able to separate through void formation, whereas the highly cross-linked structures exhibited bond scission.

Molecular dynamics↗
NASA NTRS2016

Yield and Failure Behavior Investigated for Cross-linked Phenolic Resins Using Molecular Dynamics

Molecular dynamics simulations were conducted to fundamentally evaluate the yield and failure behavior of cross-linked phenolic resins at temperatures below the glass transition. Yield stress was investigated at various temperatures, strain rates, and degrees of cross-linking. The onset of non-linear behavior in the cross-linked phenolic structures was caused by localized irreversible molecular rearrangements through the rotation of methylene linkers followed by the formation or annihilation of neighboring hydrogen bonds. The yield stress results, with respect to temperature and strain rate, could be fit by existing models used to describe yield behavior of amorphous glasses. The degree of cross-linking only indirectly influences the maximum yield stress through its influence on the glass transition temperature (Tg), however there is a strong relationship between the degree of cross-linking and the failure mechanism. Low cross-linked samples were able to separate through void formation, whereas the highly cross-linked structures exhibited bond scission.

Stress-strain curves↗
NASA NTRS2009

Prediction of Material Properties of Nanostructured Polymer Composites Using Atomistic Simulations

Atomistic models of epoxy polymers were built in order to assess the effect of structure at the nanometer scale on the resulting bulk properties such as elastic modulus and thermal conductivity. Atomistic models of both bulk polymer and carbon nanotube polymer composites were built. For the bulk models, the effect of moisture content and temperature on the resulting elastic constants was calculated. A relatively consistent decrease in modulus was seen with increasing temperature. The dependence of modulus on moisture content was less consistent. This behavior was seen for two different epoxy systems, one containing a difunctional epoxy molecule and the other a tetrafunctional epoxy molecule. Both epoxy structures were crosslinked with diamine curing agents. Multifunctional properties were calculated with the nanocomposite models. Molecular dynamics simulation was used to estimate the interfacial thermal (Kapitza) resistance between the carbon nanotube and the surrounding epoxy matrix. These estimated values were used in a multiscale model in order to predict the thermal conductivity of a nanocomposite as a function of the nanometer scaled molecular structure.

Hinkley, J.A.↗
NASA NTRS2004

Unfolding of Proteins: Thermal and Mechanical Unfolding

We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Hur, Joe S.↗
NASA NTRS2009

The Development of Directional Decohesion Finite Elements for Multiscale Failure Analysis of Metallic Polycrystals

Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.

Saether, Erik↗
NASA NTRS1998

Electronic Properties of Carbon Nanotubes and Junctions

Metallic and semiconducting Single Wall Carbon Nanotubes (CNT) have recently been characterized using scanning tunneling microscopy (STM) and the manipulation of individual CNT has been demonstrated. These developments make the prospect of using CNT as molecular wires and possibly as electronic devices an even more interesting one. We have been modeling various electronic properties such as the density of states and the transmission coefficient of CNT wires and junctions. These studies involve first calculating the stability of junctions using molecular dynamics simulations and then calculating the electronic properties using a pi-electron tight binding Hamiltonian. We have developed the expertise to calculate the electronic properties of both finite-sized CNT and CNT systems with semi-infinite boundary conditions. In this poster, we will present an overview of some of our results. The electronic application of CNT that is most promising at this time is their use as molecular wires. The conductance can however be greatly reduced because of reflection due to defects and contacts. We have modeled the transmission through CNT in the presence of two types of defects: weak uniform disorder and strong isolated scatterers. We find that the conductance is affected in significantly different manners due to these defects Junctions of CNT have also been imaged using STM. This makes it essential to derive rules for the formation of junctions between tubes of different chirality, study their relative energies and electronic properties. We have generalized the rules for connecting two different CNT and have calculated the transmission and density of states through CNT junctions. Metallic and semiconducting CNT can be joined to form a stable junction and their current versus voltage characteristics are asymmetric. CNT are deformed by the application of external forces including interactions with a substrate or other CNT. In many experiments, these deformation are expected to occur naturally. We will present some preliminary results of our calculations of the modification of CNT electronic properties as a result of deformations.

Anantram, M. P.↗
NASA NTRS2006

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized From Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, intersticies, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Glaessgen, Edward H.↗
NASA NTRS2011

Atomic Simulation of Complex DNA DSBs and the Interactions with the Ku70/80 Heterodimer

DNA double strand breaks (DSBs) induced by ionizing radiation (IR) usually contain modified bases such as 8-oxo-7,8-dihydroguanine (8-oxoG) and thymine glycol, apurinic/apyrimidinic (AP) sites, 2-deoxyribonolactone, or single-strand breaks (SSBs). The presence of such lesions in close proximity to the DSB terminus makes the DNA nicks more difficult to repair and rejoin than endogenously induced simple DSBs, and as such a major determinant of the biological effects of high linear energy transfer (LET) radiation as encountered in space travel. In this study we conducted molecular dynamics simulations on a series of DNA duplexes with various complex lesions of 8-oxoG and AP sites, in an effort to investigate the effects of such lesions to the structural integrity and stability of DNA after insulted by IR. We also simulated the interaction of such complex DSBs with the Ku70/80 heterodimer, the first protein in mammalian cells to embark the non-homologous end joining (NHEJ) DNA repair pathway. The results indicate, compared to DNA with simple DSBs, the complex lesions can enhance the hydrogen bonds opening rate at the DNA terminus, and increase the mobility of the whole duplex, thus they present more deleterious effects to the genome integrity if not captured and repaired promptly in cells. Simulations also demonstrate the binding of Ku drastically reduces structural disruption and flexibility caused by the complex lesions, and the interactions of Ku with complex DSBs have a different potential energy landscape from the bound structure with simple DSB. In all complex DSBs systems, the binding of DSB terminus with Ku70 is softened while the binding of the middle duplex with Ku80 is tightened. This energy shift may help the Ku protein to secure at the DSB terminus for a longer time, so that other end processing factors or repair pathways can proceed at the lesions before NHEJ repair process starts. These atomic simulations may provide valuable new insight into the selective action of repair proteins on damaged DNA.

Hu, Shaowen↗
NASA NTRS1992

Terraced spreading of simple liquids on solid surfaces

We have studied the spreading of liquid drops on a solid surface by molecular-dynamics simulations of coexisting three-phase Lennard-Jones systems of liquid, vapor, and solid. We consider both spherically symmetric atoms and diatomic molecules, and a range of interaction strengths. As the attraction between liquid and solid increases we observe a smooth transition in spreading regimes, from partial to complete to terraced wetting. In the terraced case, where distinct monomolecular layers spread with different velocities, the layers are ordered but not solid, with substantial molecular diffusion both within and between layers. The quantitative behavior resembles recent experimental findings, but the detailed dynamics differ. In particular, the layers exhibit an unusual spreading law, where their radii vary in time as R-squared approximately equal to log10t, which disagrees with experiments on polymeric liquids as well as recent calculations.

Yang, Ju-Xing↗
NASA NTRS1988

A molecular dynamics study of chloride binding by the cryptand SC24

The capture of chloride from water by the tetraprotonated form of the spherical macrotricyclic molecule SC24 was studied using molecular dynamics simulation methods. This model ionophore represents a broad class of molecules which remove ions from water. Two binding sites for the chloride were found, one inside and one outside the ligand. These sites are separated by a potential energy barrier of approximately 20 kcal mol-1. The major contribution to this barrier comes from dehydration of the chloride. The large, unfavorable dehydration effect is compensated for by an increase in electrostatic attraction between the oppositely charged chloride and cryptand, and by energetically favorable rearrangements of water structure. Additional assistance in crossing the barrier and completing the dehydration of the ion is provided by the shift of three positively charged hydrogen atoms of the cryptand towards the chloride. This structural rigidity is partially responsible for its selectivity.

NASA Discipline Exobiology↗
NASA NTRS2018

Simple Ion Channels: From Structure to Electrophysiology and Back

A reliable way to establish whether our understanding of a channel is satisfactory is to reproduce its measured ionic conductance over a broad range of applied voltages in computer simulations. In molecular dynamics (MD), this can be done by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. Since this approach is computationally very expensive, we have developed a markedly more efficient alternative in which MD is combined with the electrodiffusion (ED) equation. In this approach, the assumptions of the ED equation can be rigorously tested, and the precision and consistency of the calculated conductance can be determined. We have demonstrated that the full current/voltage dependence and the underlying free energy profile for a simple channel can be reliably calculated from equilibrium or non-equilibrium MD simulations at a single voltage. To carry out MD simulations, a structural model of a channel has to be assumed, which is an important constraint, considering that high-resolution structures are available for only very few simple channels. If the comparison of calculated ionic conductance with electrophysiological data is satisfactory, it greatly increases our confidence that the structure and the function are described sufficiently accurately. We examined the validity of the ED for several channels embedded in phospholipid membranes - four naturally occurring channels: trichotoxin, alamethicin, p7 from hepatitis C virus (HCV) and Vpu from the HIV-1 virus, and a synthetic, hexameric channel, formed by a 21-residue peptide that contains only leucine and serine. All these channels mediate transport of potassium and chloride ions. It was found that the ED equation is satisfactory for these systems. In some of them experimental and calculated electrophysiological properties are in good agreement, whereas in others there are strong indications that the structural models are incorrect.

molecular dynamics (MD↗
NASA NTRS2011

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian↗
NASA NTRS1997

Low-Density Aerodynamics of the Stardust Sample Return Capsule

The aerodynamics of the Stardust Sample Return Capsule are analyzed in the low- density, transitional flow regime using free-molecular, Direct Simulation Monte Carlo, Navier-Stokes, and Newtonian methods to provide inputs for constructing a transitional flow bridging relation. The accuracy of this bridging relation in reconstructing the aero- dynamic coefficients given by the more exact methods is presented for a range of flight conditions and vehicle attitudes. There is good agreement between the various prediction methods, and a simple sine-squared bridging relation is shown to provide a reasonably good description of the axial force, normal force, and pitching moment over a range of Knudsen numbers from 0.001 to 10. The predictions show a static instability of the Star- dust capsule in the free-molecular regime that persists well into the transitional flow. The addition of a thin disk to the base of the capsule is shown to remove this static instability. However, the extremely high entry velocity of 12.6 km/s for the proposed trajectory introduces difficult design issues for incorporating this disk caused by the high aerothermal loads that occur even under relatively rarefied conditions.

Wilmoth, Richard G.↗
NASA NTRS2020

Ion Channels – From Structure to Electrophysiology and Back

A reliable way to establish whether our understanding of a channel is satisfactory is to reproduce its measured ionic conductance over a broad range of applied voltages in computer simulations. In molecular dynamics (MD), this can be done by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. Since this approach is computationally very expensive, we have developed a markedly more efficient alternative in which MD is combined with the electrodiffusion (ED) equation. In this approach, the assumptions of the ED equation can be rigorously tested, and the precision and consistency of the calculated conductance can be determined. We have demonstrated that the full current/voltage dependence and the underlying free energy profile for a simple channel can be reliably calculated from equilibrium or non-equilibrium MD simulations at a single voltage. Free energy profiles can be obtained from non-equilibrium simulations without a loss of accuracy even without the knowledge of diffusion coefficient. To carry out MD simulations, a structural model of a channel has to be assumed, which is an important constraint, considering that high-resolution structures are available for only very few simple channels. If the comparison of calculated ionic conductance with electrophysiological data is satisfactory, it greatly increases our confidence that the structure and the function are described sufficiently accurately. We examined the validity of the ED for several channels embedded in phospholipid membranes - four naturally occurring channels: trichotoxin, alamethicin, p7 from hepatitis C virus (HCV) and Vpu from the HIV-1 virus, a synthetic, hexameric channel, formed by a 21-residue peptide that contains only leucine and serine and a bacterial pentameric ligand-gated ion channel, GLIC. All these channels mediate transport of potassium and chloride ions. It was found that the ED equation is satisfactory for these systems. In some of them, experimental and calculated electrophysiological properties are in good agreement, whereas in others there are strong indications that the structural models are incorrect.

Pohorille, Andrew↗
NASA NTRS1994

A random distribution reacting mixing layer model

A methodology for simulation of molecular mixing and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and non-reacting shear layer present in the facility given basic assumptions about turbulence properties.

Jones, Richard A.↗
NASA NTRS1994

A random distribution reacting mixing layer model

A methodology for simulation of molecular mixing, and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Research Center Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and nonreacting shear layer present in the facility given basic assumptions about turbulence properties.

Jones, Richard A.↗
NASA NTRS2012

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗