Search NASA⌕ Search

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.

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 NTRS2006

MODELING FUNCTIONALLY GRADED INTERPHASE REGIONS IN CARBON NANOTUBE REINFORCED COMPOSITES

A combination of micromechanics methods and molecular dynamics simulations are used to obtain the effective properties of the carbon nanotube reinforced composites with functionally graded interphase regions. The multilayer composite cylinders method accounts for the effects of non-perfect load transfer in carbon nanotube reinforced polymer matrix composites using a piecewise functionally graded interphase. The functional form of the properties in the interphase region, as well as the interphase thickness, is derived from molecular dynamics simulations of carbon nanotubes in a polymer matrix. Results indicate that the functional form of the interphase can have a significant effect on all the effective elastic constants except for the effective axial modulus for which no noticeable effects are evident.

Seidel, G. D.↗
NASA NTRS1999

Computer Simulation Studies of Ion Channel Gating: Characteristics of the M2 Channel of Influenza-A Virus in a Phospholipid Bilayer

The 25 amino acids long, transmembrane fragment of the Influenza virus M2 protein forms a homotetrameric channel that transports protons across lipid bilayers. It has been postulated that high efficiency and selectivity of this process is due to gating by four histidine residues that occlude the channel lumen in the closed state. Two mechanisms of gating have been postulated. In one mechanism, the proton is "shuttled" through the gate by attaching to the delta nitrogen atom on the extracellular side of the imidazole ring, followed by the release of the proton attached to the epsilon nitrogen atom on the opposite side. In the second mechanism, the four histidines move away from each other due to electrostatic repulsion upon protonation, thus opening the gate sufficiently that a wire of water molecules can penetrate the gate. Then, protons are transported by "hopping" along the wire. In this paper, both mechanisms are evaluated in a series of molecular dynamics simulations by investigating stability of different protonation states of the channel that are involved in these mechanisms. For the shuttle mechanism, these are states with all epsilon protonated histidines, one biprotonated residue or one histidine protonated in the delta position. For the gate opening mechanism, this is the state in which all four histidines are biprotonated. In addition, a state with two biprotonated histidines is considered. For each system, composed of the protein channel embedded in phospholipid bilayer located between two water lamellae, a molecular dynamics trajectory of approximately 1.3 ns (after equilibration) was obtained. It is found that the states involved in the shuttle mechanism are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfers involved in the shuttle. In contract, the fully biprotonated state, implicated in the gate opening mechanism, is not stable and the channel looses its structural integrity. If only two histidines are biprotonated the channel deforms but remains intact with the gate mostly closed. In summary, the results of this study lend support to the shuttle mechanism but not to the gate opening mechanism of proton gating in M2.

Schweighofer, Karl J.↗
NASA NTRS2018

Simulating the Effects of Carbon Nanotube Continuity and Interfacial Bonding on Composite Strength and Stiffness

Molecular dynamics simulations of carbon nanotube (CNT) composites, in which the CNTs are continuous across the periodic boundary, overestimate the experimentally measured mechanical properties of CNT composites along the fiber direction. Since the CNTs in these composites are much shorter than the composite dimensions, load must be transferred either directly between CNTs or through the matrix, a mechanism that is absent in simulations of effectively continuous CNTs. In this study, the elastic and fracture properties of high volume fraction discontinuous carbon nanotube/amorphous carbon composite systems were compared to those of otherwise equivalent continuous CNT composites using ReaxFF reactive molecular dynamics simulations. These simulations were used to show how the number of nanotube-matrix interfacial covalent bonds affect composite mechanical properties. Furthermore, the mechanical impact of interfacial bonding was decomposed to reveal its effect on the properties of the CNTs, the interfacial layer of matrix, and the bulk matrix. For the composites with continuous reinforcement, it was found that any degree of interfacial bonding has a negative impact on axial tensile strength and stiffness. This is due to disruption of the structure of the CNTs and interfacial matrix layer by the interfacial bonds. For the discontinuous composites, the modulus was maximized between 4%-7% interfacial bonding and the strength continues to increase up to the highest levels of interfacial bonding studied. Areas of low stress and voids were observed in the simulated discontinuous composites at the ends of the tubes, from which fracture was observed to initiate. Experimental carbon nanotube yarn composites were fabricated and tested. The results were used to illustrate knockdown factors relative to the mechanical performance of the tubes themselves.

Benjamin D Jensen↗
NASA NTRS2000

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗
NASA NTRS2013

Conductance of Ion Channels - Theory vs. Experiment

Transmembrane ion channels mediate a number of essential physiological processes in a cell ranging from regulating osmotic pressure to transmission of neural signals. Kinetics and selectivity of ion transport is of critical importance to a cell and, not surprisingly, it is a subject of numerous experimental and theoretical studies. In this presentation we will analyze in detail computer simulations of two simple channels from fungi - antiamoebin and trichotoxin. Each of these channels is made of an alpha-helical bundle of small, nongenomically synthesized peptides containing a number of rare amino acids and exhibits strong antimicrobial activity. We will focus on calculating ionic conductance defined as the ratio of ionic current through the channel to applied voltage. From molecular dynamics simulations, conductance can be calculated in at least two ways, each involving different approximations. Specifically, the current, given as the number of charges transferred through the channel per unit of time, can be obtained from the number of events in which ions cross the channel during the simulation. This method works well for large currents (high conductance values and/or applied voltages). If the number of crossing events is small, reliable estimates of current are difficult to achieve. Alternatively, conductance can be estimated assuming that ion transport can be well approximated as diffusion in the external potential given by the free energy profile. Then, the current can be calculated by solving the one-dimensional diffusion equation in this external potential and applied voltage (the generalized Nernst-Planck equation). To do so three ingredients are needed: the free energy profile, the position-dependent diffusion coefficient and the diffusive flux of ions into the channel. All these quantities can be obtained from molecular dynamics simulations. An important advantage of this method is that it can be used equally well to estimating large and small currents. In addition, once the free energy profile becomes available the full current-voltage dependence can be readily obtained. For both channels we carried out calculations using both approaches. We also tested the main assumptions underlying the diffusive model, such as uncorrelated nature of individual crossing events and Fickian diffusion. The accuracy and consistency of different methods will be discussed. Finally we will discuss how comparisons between calculated and measured ionic conductance and selectivity of transport can be used for determining structural models of the channels.

Pohorille, Andrew↗
NASA NTRS2015

Fracture of Carbon Nanotube - Amorphous Carbon Composites: Molecular Modeling

Carbon nanotubes (CNTs) are promising candidates for use as reinforcements in next generation structural composite materials because of their extremely high specific stiffness and strength. They cannot, however, be viewed as simple replacements for carbon fibers because there are key differences between these materials in areas such as handling, processing, and matrix design. It is impossible to know for certain that CNT composites will represent a significant advance over carbon fiber composites before these various factors have been optimized, which is an extremely costly and time intensive process. This work attempts to place an upper bound on CNT composite mechanical properties by performing molecular dynamics simulations on idealized model systems with a reactive forcefield that permits modeling of both elastic deformations and fracture. Amorphous carbon (AC) was chosen for the matrix material in this work because of its structural simplicity and physical compatibility with the CNT fillers. It is also much stiffer and stronger than typical engineering polymer matrices. Three different arrangements of CNTs in the simulation cell have been investigated: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. The SWNT and MWNT array systems are clearly idealizations, but the SWNT bundle system is a step closer to real systems in which individual tubes aggregate into large assemblies. The effect of chemical crosslinking on composite properties is modeled by adding bonds between the CNTs and AC. The balance between weakening the CNTs and improving fiber-matrix load transfer is explored by systematically varying the extent of crosslinking. It is, of course, impossible to capture the full range of deformation and fracture processes that occur in real materials with even the largest atomistic molecular dynamics simulations. With this limitation in mind, the simulation results reported here provide a plausible upper limit on achievable CNT composite properties and yield some insight on the influence of processing conditions on the mechanical properties of CNT composites.

Jensen, Benjamin D.↗
NASA NTRS1993

Molecular gas dynamics applied to low-thrust propulsion

The Direct Simulation Monte Carlo method is currently being applied to study flowfields of small thrusters, including both the internal nozzle and the external plume flow. The DSMC method is employed because of its inherent ability to capture nonequilibrium effects and proper boundary physics in low-density flow that are not readily obtained by continuum methods. Accurate prediction of both the internal and external nozzle flow is important in determining plume expansion which, in turn, bears directly on impingement and contamination effects.

Zelesnik, Donna↗
NASA NTRS1994

A comparison of dynamical and molecular gas masses in very luminous infrared galaxies

Observations of the 2.3 micron bands of CO have been used to measure the central stellar velocity dispersions in NGC 1614 and IC 694. K-band images were used to constrain models of the stellar mass distribution. With these data and a spectrum of Arp 220 from the literature, we have determined the total mass within the central regions of all three galaxies. The central molecular gas masses inferred from (12)CO J = 1-0 fluxes are larger than the total masses determined from stellar kinematics. Subtraction of the likely stellar mass from the total mass yields upper limits for the H2 mass that are 4-10 times less than the previously reported masses. This result supports other arguments that the I(sub CO)/M(sub H2) ratio is not the same in ultraluminous galaxies as it is in Milky Way giant molecular clouds.

Shier, L. M.↗
NASA NTRSDate not supplied

Multiscale Analysis of Thermo-Mechanical Behavior of Boron Nitride-Reinforced Epoxy Nanocomposites

The effect of size, shape, morphology, and arrangement of micro constituents of Boron Nitrite (BN) nanoplatelet/epoxy composites on their properties were investigated using a multi scale approach that includes Molecular Dynamics (MD) and micromechanics. The thermo-mechanical properties of the composites were evaluated using molecular theory and the analysis showed that the elastic constants of BN/epoxy composites were not severely affected by temperature. Also, the micromechanical analysis of Generalized Method of Cells (GMC) was utilized at higher length scale to evaluate elastic properties of the composites, for different geometries and arrangements of micro constituents. The predicted results of the analysis showed that the size, aspect ratios, morphology, and the arrangements of inclusions in BN/epoxy nanocomposites all have remarkable effect on the mechanical performances of the material systems.

nanocomposites↗
NASA NTRS1980

Clumpy molecular clouds - A dynamic model self-consistently regulated by T Tauri star formation

A model is proposed which can account for the longevity, energetics, and dynamical structure of dark molecular clouds. Recent observations of a high space density of T Tauri stars in some dark clouds provide the basis for the assertion that high-velocity winds from these low mass pre-main-sequence stars provide a continuous dynamic input into molecular clouds. It is postulated that as clumps are driven above the Jeans mass, both by coalescence and the enhancement of ram pressure through continuing acceleration by protostellar winds, collapse is followed by the formation of low-mass stars that generate additional protostellar winds. Finally, it is found that star formation occurs on a relatively slow time scale, comparable to the cloud lifetime

Norman, C.↗
NASA NTRS2004

Cross-Linked Nanotube Materials with Variable Stiffness Tethers

The constitutive properties of a cross-linked single-walled carbon nanotube material are predicted with a multi-scale model. The material is modeled as a transversely isotropic solid using concepts from equivalent-continuum modeling. The elastic constants are determined using molecular dynamics simulation. Some parameters of the molecular force field are determined specifically for the cross-linker from ab initio calculations. A demonstration of how the cross-linked nanotubes may affect the properties of a nanotube/polyimide composite is included using a micromechanical analysis.

Frankland, Sarah-Jane V.↗
NASA NTRS2001

Ranges of applicability for the continuum beam model in the mechanics of carbon nanotubes and nanorods

Limitations in the validity of the continuum beam model for carbon nanotubes (NTs) and nanorods are examined. Applicability of all assumptions used in the model is restricted by the two criteria for geometric parameters that characterize the structure of NTs. The key non-dimensional parameters that control the NT buckling behavior are derived via dimensional analysis of the nanomechanical problem. A mechanical law of geometric similitude for NT buckling is extended from continuum mechanics for different molecular structures. A model applicability map, where two classes of beam-like NTs are identified, is constructed for distinct ranges of non-dimensional parameters. Expressions for the critical buckling loads and strains are tailored for two classes of NTs and compared with the data provided by the molecular dynamics simulations. copyright 2001 Elsevier Science Ltd. All rights reserved.

Continuum beam model↗
NASA NTRSDate not supplied

Computationally Guided Design of Polymer-Coated Microparticles as Reusable Materials

Long-duration space exploration missions and sustained lunar or Martian surface operations present greater demands for multifunctional and reusable materials. By scaling down the amount of material to be launched from Earth, both mission cost and risk can be reduced. In this regard, leveraging in-space manufacturing capabilities with reusable feedstock materials is an attractive option, as it will allow for articles to be generated on demand, utilized, and then recycled for additional use. NASA’s Enabling Sustained Presence Using Recyclables (ESPUR) project aims to develop reusable materials using polymer-coated microparticles that are bonded via reversible Diels-Alder reactions, where only modest heat is needed to trigger the reverse reaction and enable reuse. For proof-of-concept demonstration, research is currently focused on the fabrication of epoxy microparticles that contain a copoly(carbonate urethane) coating with maleimide and furan functionalities. Here, we discuss the integration of computational materials modeling approaches to help navigate the large design space in this development effort. We perform molecular dynamics (MD) simulations with atomistic and coarse-grained models of the copolymer, which allow us to evaluate the effects of design parameters like the molecular weight and composition on the molecular interactions and chain dynamics. We show how the properties change with the reversible bonds. We also leverage discrete element method (DEM) simulations to assess how microparticle design parameters like the size ratio and volume fraction can be tuned to increase the packing density and number of microparticle contacts to improve the mechanical properties. Our results demonstrate how computational tools can be used in close collaboration with experimental efforts to accelerate material design.

reusable materials↗
NASA NTRSDate not supplied

Modeling Ring-Opening Polymerization with Machine-Learned Committors

Many important industrial and biological processes depend on ring-opening polymerization (ROP), including large-scale production polymers, such as nylon 6, and a variety of biodegradable polymers and high-performance resins. Large-scale simulations (>200K atoms) of ROP were performed using a deep neural network-based method that captured the effect of the local chemical environment of the reactive sites using training data from quantum chemical methods. The utility of this method was demonstrated by modeling the cationic ROP of bis-benzoxazine, a promising resin for use in advanced high-temperature application composites (>200°C). Density functional theory calculations were used to map out the various ring opening and polymerization reaction paths that occur during the thermal processing of bis-benzoxazine. The probabilities of these reactions, as a function of local atomic configurations, were used to construct committor functions using a neural network approach. The learned committor functions were then used to determine when a particular reaction occurred within the REACTER framework for reactive molecular dynamics simulations. The experimentally observed exotherm that occurs during bis-benzoxazine ROP was captured in an effective manner in the simulations by imposing a local heating of the system that corresponds to the enthalpy of reaction for each reaction. The dependence of the final polymer morphology on the ROP reaction temperature was analyzed at the molecular scale. The results presented in this work indicate that machine-learned committor functions are a promising approach for incorporating high-fidelity reactivity criteria that capture the effect of local chemical environments into large-scale polymerization simulations.

machine learning, polymer simulations, molecular d↗
NASA NTRS2008

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗
NASA NTRS2022

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗