SEARCH · Search NASA

Results for “Molecular Dynamics”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Multi-Scale Modeling of a Graphite-Epoxy-Nanotube System

A multi-scale method is utilized to determine some of the constitutive properties of a three component graphite-epoxy-nanotube system. This system is of interest because carbon nanotubes have been proposed as stiffening and toughening agents in the interlaminar regions of carbon fiber/epoxy laminates. The multi-scale method uses molecular dynamics simulation and equivalent-continuum modeling to compute three of the elastic constants of the graphite-epoxy-nanotube system: C11, C22, and C33. The 1-direction is along the nanotube axis, and the graphene sheets lie in the 1-2 plane. It was found that the C11 is only 4% larger than the C22. The nanotube therefore does have a small, but positive effect on the constitutive properties in the interlaminar region.

Frankland, S. J. V.↗

NASA NTRS2005

An Application-Based Performance Characterization of the Columbia Supercluster

Columbia is a 10,240-processor supercluster consisting of 20 Altix nodes with 512 processors each, and currently ranked as the second-fastest computer in the world. In this paper, we present the performance characteristics of Columbia obtained on up to four computing nodes interconnected via the InfiniBand and/or NUMAlink4 communication fabrics. We evaluate floating-point performance, memory bandwidth, message passing communication speeds, and compilers using a subset of the HPC Challenge benchmarks, and some of the NAS Parallel Benchmarks including the multi-zone versions. We present detailed performance results for three scientific applications of interest to NASA, one from molecular dynamics, and two from computational fluid dynamics. Our results show that both the NUMAlink4 and the InfiniBand hold promise for application scaling to a large number of processors.

Biswas, Rupak↗

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 NTRS2006

Micromechanics Modeling of Functionally Graded Interphase Regions in Carbon Nanotube-Polymer Composites

The effective elastic properties of a unidirectional carbon fiber/epoxy lamina in which the carbon fibers are coated with single-walled carbon nanotubes are modeled herein through the use of a multi-scale method involving the molecular dynamics/equivalent continuum and micromechanics methods. The specific lamina representative volume element studied consists of a carbon fiber surrounded by a region of epoxy containing a radially varying concentration of carbon nanotubes which is then embedded in the pure epoxy matrix. The variable concentration of carbon nanotubes surrounding the carbon fiber results in a functionally graded interphase region as the properties of the interphase region vary according to the carbon nanotube volume fraction. Molecular dynamics and equivalent continuum methods are used to assess the local effective properties of the carbon nanotube/epoxy comprising the interphase region. Micromechanics in the form of the Mori-Tanaka method are then applied to obtain the global effective properties of the graded interphase region wherein the carbon nanotubes are randomly oriented. Finally, the multi-layer composite cylinders micromechanics approach is used to obtain the effective lamina properties from the lamina representative volume element. It was found that even very small quantities of carbon nanotubes (0.36% of lamina by volume) coating the surface of the carbon fibers in the lamina can have a significant effect (8% increase) on the transverse properties of the lamina (E22, k23, G23 and G12) with almost no affect on the lamina properties in the fiber direction (E11 and v12).

Seidel, Gary D.↗

NASA NTRS2004

Antigenic Protein In Microgravity-Grown Human Mixed Mullerian Tumor (LN1) Cells Preserved In RNA Stabilizing Agent

Cells treated with RNAlater(TradeMark) have previously been shown to contain antigenic proteins that can be visualized using Western blot analysis. These proteins seem to be stable for several months when stored in RNA stabilizer at 4 C. Antigenic protein can be recovered from cells that have been processed using an Ambion RNAqueous(Registered TradeMark) kit to remove RNA. In this set of experiments, human mixed Mullerian tumor (LNI) cells grown on the International Space Station during Expedition 3 were examined for antigenic stability after removal of RNA. The cells were stored for three months in RNAlater(TradeMark) and RNA was extracted. The RNA filtrate containing the protein was precipitated, washed, and suspended in buffer containing sodium dodecyl sulfate (SDS). Samples containing equal concentrations of protein were loaded onto SDS-polyacrylamide gels. Proteins were separated by electrophoresis and transferred by Western blot to polyvinylidene fluoride (PVDF) membrane. The Western blots were stained with an enhanced chemiluminescent ECL(Registered Trademark) Plus detection kit (Amersham) and scanned using a Storm 840 gel image analyzer (Amersham, Molecular Dynamics). ImageQuant(Registered TradeMark) software was used to quantify the densities of the protein bands. The ground control and flight LN1 cell samples showed a similar staining pattern over time with antibodies to vimentin, glyceraldehyde-3-phosphate dehydrogenase, and epithelial membrane antigens.

Hammond, Dianne K.↗

NASA NTRS2006

Computer Modeling of Protocellular Functions: Peptide Insertion in Membranes

Lipid vesicles became the precursors to protocells by acquiring the capabilities needed to survive and reproduce. These include transport of ions, nutrients and waste products across cell walls and capture of energy and its conversion into a chemically usable form. In modem organisms these functions are carried out by membrane-bound proteins (about 30% of the genome codes for this kind of proteins). A number of properties of alpha-helical peptides suggest that their associations are excellent candidates for protobiological precursors of proteins. In particular, some simple a-helical peptides can aggregate spontaneously and form functional channels. This process can be described conceptually by a three-step thermodynamic cycle: 1 - folding of helices at the water-membrane interface, 2 - helix insertion into the lipid bilayer and 3 - specific interactions of these helices that result in functional tertiary structures. Although a crucial step, helix insertion has not been adequately studied because of the insolubility and aggregation of hydrophobic peptides. In this work, we use computer simulation methods (Molecular Dynamics) to characterize the energetics of helix insertion and we discuss its importance in an evolutionary context. Specifically, helices could self-assemble only if their interactions were sufficiently strong to compensate the unfavorable Free Energy of insertion of individual helices into membranes, providing a selection mechanism for protobiological evolution.

Rodriquez-Gomez, D.↗

NASA NTRS2006

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗

NASA NTRS2006

Experimental Partitioning of Cr(3+) and Sc(3+) into Olivine: Mechanisms and Implications

Olivine (Mg, Fe)Si2O4 does not, by stoichiometry, accept cations such as Sc(3+) or Cr(3+). However, the partition coefficients of Sc and Cr between olivine and liquid are significant 0.2-1.0. We have measured Cr(3+) partition coefficients of near unity and have grown olivines with nearly 3 wt.% Sc2O3. Therefore, there must be a simple means of charge balancing 3+ ions in a crystal structure that was obviously not designed to receive other than 2+ ions on the olivine M sites. The simplest explanation is that two 3+ ions enter the olivine structure by displacing three 2+ ions and creating an M site vacancy. Even this explanation has difficulties. For minor elements in our experiments (~1 wt.%) the odds of a minor element 3+ ion finding a second 3+ for charge balance are of the order of 100:1 against. Because of the reducing conditions of our experiments, Fe(3+) will not suffice; and Al(3+) is not in sufficient quantity in olivine for charge balance. Therefore, Cr or Sc must, in effect, charge balance itself. For true trace elements, the problem is compounded many times. For an ion at the 10 ppm level the chances of finding a second (for example) Sc ion is approx.10(exp 5):1 against. Of course, any other 3+ ion would suffice but comparisons between percent level doping experiments and trace level partitioning indicate that Henry s law is obeyed. This implies that the same substitution mechanism occurs at both the percent and tens of ppm levels. There are two simple solutions to this problem: (i) The electrical conductivity of olivine is such that charge balance need not be local. This requires substantial domains within the olivine crystal in electrical contact by migration of vacancies or electronic defects. (ii) The 3+ cation brings along its own charge-balancing ion because it existed as a dimer in the silicate liquid. Olivine is not a true insulator but is actually a p-type semiconductor. Even so, electrical communication by this means is unlikely over the tens or hundreds of unit cells that would be required for charge balance to be local. Therefore, we cautiously favor the idea that melt speciation is the means by which 3+ ions enter the olivine structure. Possibly this model might be tested by in situ XAFS measurements or by molecular dynamical calculations.

Jones, John↗

NASA NTRS2004

Antigenic Protein In Microgravity-Grown Human Mixed Mullerian Tumor (LN1) Cells Preserved In RNA Stabilizing Agent

Cells treated with RNAlater(TradeMark) have previously been shown to contain antigenic proteins that can be visualized using Western blot analysis. These proteins seem to be stable for several months when stored in RNA stabilizer at 4 C. Antigenic protein can be recovered from cells that have been processed using an Ambion RNAqueous(Registered TradeMark) kit to remove RNA. In this set of experiments, human mixed Mullerian tumor (LN1) cells grown on the International Space Station during Expedition 3 were examined for antigenic stability after removal of RNA. The cells were stored for three months in RNAlater(TradeMark) and RNA was extracted. The RNA filtrate Containing the protein was precipitated, washed, and suspended in buffer containing sodium dodecyl sulfate (SDS). Samples containing equal concentrations of protein were loaded onto SDS-polyacrylamide gels. Proteins were separated by electrophoresis and transferred by Western blot to polyvinylidene fluoride (PVDF) membrane. The Western blots were stained with an enhanced chemiluminescent ECL(Registered TradeMark)Plus detection kit (Amersham) and scanned using a Storm 840 gel image analyzer (Amersham, Molecular Dynamics). ImageQuant(Registered TradeMark)a software was used to quantify the densities of the protein bands. The ground control and flight LN1 cell samples showed a similar staining pattern over time with antibodies to vimentin, glyceraldehyde-3-phosphate dehydrogenase, and epithelial membrane antigens.

Hammond, Dianne K.↗

NASA NTRS2008

Multiscale Modeling of Intergranular Fracture in Aluminum: Constitutive Relation For Interface Debonding

Intergranular fracture is a dominant mode of failure in ultrafine grained materials. In the present study, the atomistic mechanisms of grain-boundary debonding during intergranular fracture in aluminum are modeled using a coupled molecular dynamics finite element simulation. Using a statistical mechanics approach, a cohesive-zone law in the form of a traction-displacement constitutive relationship, characterizing the load transfer across the plane of a growing edge crack, is extracted from atomistic simulations and then recast in a form suitable for inclusion within a continuum finite element model. The cohesive-zone law derived by the presented technique is free of finite size effects and is statistically representative for describing the interfacial debonding of a grain boundary (GB) interface examined at atomic length scales. By incorporating the cohesive-zone law in cohesive-zone finite elements, the debonding of a GB interface can be simulated in a coupled continuum-atomistic model, in which a crack starts in the continuum environment, smoothly penetrates the continuum-atomistic interface, and continues its propagation in the atomistic environment. This study is a step towards relating atomistically derived decohesion laws to macroscopic predictions of fracture and constructing multiscale models for nanocrystalline and ultrafine grained materials.

Yamakov, V.↗

NASA NTRS2006

Modeling of Interfacial Modification Effects on Thermal Conductivity of Carbon Nanotube Composites

The effect of functionalization of carbon nanotubes on the thermal conductivity of nanocomposites has been studied using a multi-scale modeling approach. These results predict that grafting linear hydrocarbon chains to the surface of a single wall carbon nanotube with covalent chemical bonds should result in a significant increase in the thermal conductivity of these nanocomposites. This is due to the decrease in the interfacial thermal (Kapitza) resistance between the single wall carbon nanotube and the surrounding polymer matrix upon chemical functionalization. The nanocomposites studied here consist of single wall carbon nanotubes in a bulk poly(ethylene vinyl acetate) matrix. The nanotubes are functionalized by end-grafting linear hydrocarbon chains of varying length to the surface of the nanotube. The effect which this functionalization has on the interfacial thermal resistance is studied by molecular dynamics simulation. Interfacial thermal resistance values are calculated for a range of chemical grafting densities and with several chain lengths. These results are subsequently used in an analytical model to predict the resulting effect on the bulk thermal conductivity of the nanocomposite.

Clancy, Thomas C.↗

NASA NTRS2008

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗

NASA NTRS2008

State of the Art Assessment of Simulation in Advanced Materials Development

Advances in both the underlying theory and in the practical implementation of molecular modeling techniques have increased their value in the advanced materials development process. The objective is to accelerate the maturation of emerging materials by tightly integrating modeling with the other critical processes: synthesis, processing, and characterization. The aims of this report are to summarize the state of the art of existing modeling tools and to highlight a number of areas in which additional development is required. In an effort to maintain focus and limit length, this survey is restricted to classical simulation techniques including molecular dynamics and Monte Carlo simulations.

Wise, Kristopher E.↗

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 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.↗

NASA NTRS2009

Filtering Water by Use of Ultrasonically Vibrated Nanotubes

Devices that could be characterized as acoustically driven molecular sieves have been proposed for filtering water to remove all biological contaminants and all molecules larger than water molecules. Originally intended for purifying wastewater for reuse aboard spacecraft, these devices could also be attractive for use on Earth in numerous settings in which there are requirements to obtain potable, medical-grade, or otherwise pure water from contaminated water supplies. These devices could also serve as efficient means of removing some or all water from chemical products . for example, they might be useful as adjuncts or substitutes for stills in the removal of water from alcohols and alcoholic beverages. These devices may be constructed using various materials, such as ceramics, metallics, or polymers, depending on end-use requirements. A representative device of this type (see figure) would include a polymeric disk, about 1 mm in diameter and between 1 and 40 microns thick, within which would be embedded single-wall carbon nanotubes aligned along the thickness axis. The polymeric disk would be part of a unitary polymeric ring assembly. An acoustic transducer in the form of a piezoelectric-film-and-electrode subassembly - typically 9 microns thick and made of poly(vinylidene fluoride) coated with copper 150 nm thick -. would be affixed to the outside of the outer polymeric ring by means of an electrically nonconductive epoxy. The nanotubes would be chosen to have diameters between about 8 and about 13.5 A because water molecules could fit into the nanotubes, but larger molecules could not. Water to be purified would be placed in contact with one face (typically, the upper face) of the filter disk. The surface tension of water is low enough that water molecules should enter and travel along the nanotubes, and computational simulations of molecular dynamics and experimental measurements have shown that the water molecules inside the nanotubes in this size range can be expected to become aligned into helical columns that exhibit properties of both hexagonal ice crystals and liquid water

Gavalas, Lillian Susan↗

NASA NTRS2010

Multiscale Modeling of Thermal Conductivity of Polymer/Carbon Nanocomposites

Molecular dynamics simulation was used to estimate the interfacial thermal (Kapitza) resistance between nanoparticles and amorphous and crystalline polymer matrices. Bulk thermal conductivities of the nanocomposites were then estimated using an established effective medium approach. To study functionalization, oligomeric ethylene-vinyl alcohol copolymers were chemically bonded to a single wall carbon nanotube. The results, in a poly(ethylene-vinyl acetate) matrix, are similar to those obtained previously for grafted linear hydrocarbon chains. To study the effect of noncovalent functionalization, two types of polyethylene matrices. -- aligned (extended-chain crystalline) vs. amorphous (random coils) were modeled. Both matrices produced the same interfacial thermal resistance values. Finally, functionalization of edges and faces of plate-like graphite nanoparticles was found to be only modestly effective in reducing the interfacial thermal resistance and improving the composite thermal conductivity

Clancy, Thomas C.↗

NASA NTRS2010

Modeling Initial Stage of Ablation Material Pyrolysis: Graphitic Precursor Formation and Interfacial Effects

Reactive molecular dynamics simulations are used to study initial stage of pyrolysis of ablation materials and their composites with carbon nanotubes and carbon fibers. The products formed during pyrolysis are characterized and water is found as the primary product in all cases. The water formation mechanisms are analyzed and the value of the activation energy for water formation is estimated. A detailed study on graphitic precursor formation reveals the presence of two temperature zones. In the lower temperature zone (less than 2000 K) polymerization occurs resulting in formation of large, stable graphitic precursors, and in the high temperature zone (greater than 2000 K) polymer scission results in formation of short polymer chains/molecules. Simulations performed in the high temperature zone on the phenolic resin composites (with carbon nanotubes and carbon fibers) shows that the presence of interfaces had no substantial effect on the chain scission rate or the activation energy value for water formation.

Desai, Tapan G.↗