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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗
NASA NTRS1997

Toroidal Single Wall Carbon Nanotubes in Fullerene Crop Circles

We investigate energetics and structure of circular and polygonal single wall carbon nanotubes (SWNTs) using large scale molecular simulations on NAS SP2, motivated by their unusual electronic and magnetic properties. The circular tori are formed by bending tube (no net whereas the polygonal tori are constructed by turning the joint of two tubes of (n, n), (n+1, n-1) and (n+2, n-2) with topological pentagon-heptagon defect, in which n =5, 8 and 10. The strain energy of circular tori relative to straight tube decreases by I/D(sup 2) where D is torus diameter. As D increases, these tori change from buckling to an energetically stable state. The stable tori are perfect circular in both toroidal and tubular geometry with strain less than 0. 03 eV/atom when D greater than 10, 20 and 40 nm for torus (5,5), (8,8) and (10, 10). Polygonal tori, whose strain is proportional to the number of defects and I/D are energetically stable even for D less than 10 nm. However, their strain is higher than that of perfect circular tori. In addition, the local maximum strain of polygonal tori is much higher than that of perfect circular tori. It is approx. 0.03 eV/atom or less for perfect circular torus (5,5), but 0.13 and 0.21 eV/atom for polygonal tori (6,4)/(5,5) and (7,3)/(5,5). Therefore, we conclude that the circular tori with no topological defects are more energetically stable and kinetically accessible than the polygonal tori containing the pentagon-heptagon defects for the laser-grown SWNTs and Fullerene crop circles.

Han, Jie↗
NASA NTRS2003

CO2 Acquisition Membrane (CAM)

The objective of CAM is to develop, test, and analyze thin film membrane materials for separation and purification of carbon dioxide (CO2) from mixtures of gases, such as those found in the Martian atmosphere. The membranes are targeted toward In Situ Resource Utilization (ISRU) applications that will operate in extraterrestrial environments and support future unmanned and human space missions. A primary application is the Sabatier Electrolysis process that uses Mars atmosphere CO2 as raw material for producing water, oxygen, and methane for rocket fuel and habitat support. Other applications include use as an inlet filter to collect and concentrate Mars atmospheric argon and nitrogen gases for habitat pressurization, and to remove CO2 from breathing gases in Closed Environment Life Support Systems (CELSS). CAM membrane materials include crystalline faujasite (FAU) zeolite and rubbery polymers such as silicone rubber (PDMS) that have been shown in the literature and via molecular simulation to favor adsorption and permeation of CO2 over nitrogen and argon. Pure gas permeation tests using commercial PDMS membranes have shown that both CO2 permeance and the separation factor relative to other gases increase as the temperature decreases, and low (Delta)P(Sub CO2) favors higher separation factors. The ideal CO2/N2 separation factor increases from 7.5 to 17.5 as temperature decreases from 22 C to -30 C. For gas mixtures containing CO2, N2, and Ar, plasticization decreased the separation factors from 4.5 to 6 over the same temperature range. We currently synthesize and test our own Na(+) FAU zeolite membranes using standard formulations and secondary growth methods on porous alumina. Preliminary tests with a Na(+) FAU membrane at 22 C show a He/SF6 ideal separation factor of 62, exceeding the Knudsen diffusion selectivity by an order of magnitude. This shows that the membrane is relatively free from large defects and associated non-selective (viscous flow) transport mechanisms. The Membrane Test Facility (MTF) has been developed to measure membrane permeance over a wide range of temperature and pressure. The facility uses two volume compartments separated by the membrane that are instrumented to measure temperature, delta pressure across the membrane, and gas composition. A thermal shroud supports and encloses the membrane, and provides temperature control. Methods were developed to determine membrane permeance using the first order decay of the pressure difference between the sealed compartments, using the total pressure for pure gases, and partial pressure of each species in gas mixtures. The technique provides an end-to-end measurement of gas permeance that includes concentration polarization effects. Experiments have shown that in addition to membrane permeance properties, the geometry and design of associated structures play an important role in how membrane systems will function on Mars.

Mason, Larry W.↗
NASA NTRS2004

Novel Microstructures for Polymer-Liquid Crystal Composite Materials

There are a number of interface-dominated composite materials that contain a liquid crystalline (LC) phase in intimate contact with an isotropic phase. For example, polymer- dispersed liquid crystals, used in the fabrication of windows with switchable transparency, consist of micron size LC droplets dispersed in an isotropic polymer matrix. Many other types of liquid crystal composite materials can be envisioned that might have outstanding optical properties that could be exploited in novel chemical sensors, optical switches, and computer displays. This research project was based on the premise that many of these potentially useful LC composite materials can only be fabricated under microgravity conditions where gravity driven flows are absent. In the ground-based research described below, we have focused on a new class of LC composites that we call thermotropic- lyotropic liquid crystal systems (TLLCs). TLLCs consist of nanosize droplets of water dispersed in an LC matrix, with surfactants at the interface that stabilize the structure. By varying the type of surfactant one can access almost an infinite variety of unusual LC composite microstructures. Due to the importance of the interface in these types of systems, we have also developed molecular simulation models for liquid crystals at interfaces, and made some of the first measurements of the interfacial tension between liquid crystals and water.

Magda, Jules J.↗
NASA NTRS2008

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspect of out understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to to simulate these process in numerical models of chemistry and transport. These models depend on laboratory-measured kinetic reaction rates and photlysis cross section to simulate molecular interactions. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluation are applied in random combinations. We determine the key reaction and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Kawa, S. Randolph↗
NASA NTRS2008

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspects of our understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to simulate these processes in numerical models of chemistry and transport. The fidelity of the models is assessed in comparison with a wide range of observations. These models depend on laboratory-measured kinetic reaction rates and photolysis cross sections to simulate molecular interactions. A typical stratospheric chemistry mechanism has on the order of 50- 100 species undergoing over a hundred intermolecular reactions and several tens of photolysis reactions. The rates of all of these reactions are subject to uncertainty, some substantial. Given the complexity of the models, however, it is difficult to quantify uncertainties in many aspects of system. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluations are applied in random combinations. We determine the key reactions and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Kawa, S. Randolph↗
NASA NTRS2018

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is a suite of tools for the analysis of porous materials and generation of material microstructures. From microstructural data, often obtained through X-ray microtomography, PuMA can determine a number of effective material properties and perform material response simulations. Version 2.2 includes capabilities for computing volume fractions, porosity, specific surface area, effective thermal and electrical conductivities, and continuum and rarefied diffusive tortuosity. PuMA can also simulate competitive diffusion/reaction processes at the micro-scale, such as surface oxidation. In this poster, recent advancements to the PuMA software are detailed, including the full refactoring of PuMA into v3.0, a new module to compute heat conduction in anisotropic materials, a particle method for simulating molecular beam experiments, a new finite-volume Laplace solver, complex fibrous material generation, woven material generation, and a coupling of PuMA with the DAKOTA software for advanced statistics.

PuMA↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗
NASA NTRSDate not supplied

REACTER 2.0: Quantum-Informed Reaction Constraints and Automated Interaction Typing

REACTER is a heuristic method for modeling chemical reactions in classical molecular dynamics simulations, implemented in LAMMPS as fix bond/react. The authors recently extended LAMMPS to support alphanumeric labels for atom types, bond types etc., which enables the pre- and post-reaction templates required by the REACTER protocol to be portable between different simulations and greatly simplifies the task of creating simulation-ready reaction templates. To further increase the generality of reaction templates, support for wildcard characters within atom types has been added, along with the automatic assignment of interaction types for new bonds, angles, etc. based on the involved atom types. In some cases, this feature can express a class of reactions with one pair of reaction templates, where previously dozens may have been required. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, as well as a custom constraint that leverages LAMMPS’ powerful built-in variable framework. Other new features include variable support for various inputs (e.g., to allow reaction rates or cutoffs to be dependent on overall conversion), on-the-fly update of molecule IDs, and the ability to create new atoms positioned with respect to the reaction site. The new features are applied to modeling polymeric, thermosetting and composite materials, and advanced applications of the new reaction constraints are demonstrated. For example, REACTER is shown to accurately reproduce mechanically-induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a constraint on the total potential energy of the involved atoms.

polymer simulations↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗
NASA NTRS1993

Simulations of a molecular plasma in collisional-radiative nonequilibrium

A code for the simulation of nonequilibrium plasmas is being developed, with the capability to couple the plasma fluid-dynamics for a single fluid with a collisional-radiative model, where electronic states are treated as separate species. The model allows for non-Boltzmann distribution of the electronic states. Deviations from the Boltzmann distributions are expected to occur in the rapidly ionizing regime behind a strong shock or in the recombining regime during a fast expansion. This additional step in modeling complexity is expected to yield more accurate predictions of the nonequilibrium state and the radiation spectrum and intensity. An attempt at extending the code to molecular plasma flows is presented. The numerical techniques used, the thermochemical model, and the results of some numerical tests are described.

Cambier, Jean-Luc↗
NASA NTRS2001

Ranges of Applicability for the Continuum-beam Model in the Constitutive Analysis of Carbon Nanotubes: Nanotubes or Nano-beams?

Ranges of validity for the continuum-beam model, the length-scale effects and continuum assumptions are analyzed in the framework of scaling analysis of NT structure. Two coupled criteria for the applicability of the continuum model are presented. Scaling analysis of NT buckling and geometric parameters (e.g., diameter and length) is carried out to determine the key non-dimensional parameters that control the buckling strains and modes of NT buckling. A model applicability map, which represents two classes of NTs, is constructed in the space of non-dimensional parameters. In an analogy with continuum mechanics, a mechanical law of geometric similitude is presented for two classes of beam-like NTs having different geometries. Expressions for the critical buckling loads and strains are tailored for the distinct groups of NTs and compared with the data provided by the molecular dynamics simulations. Implications for molecular dynamics simulations and the NT-based scanning probes are discussed.

Harik, Vasyl Michael↗
NASA NTRS1984

A low Earth orbit molecular beam space simulation facility

A brief synopsis of the low Earth orbit (LEO) satellite environment is presented including neutral and ionic species. Two ground based atomic and molecular beam instruments are described which are capable of simulating the interaction of spacecraft surfaces with the LEO environment and detecting the results of these interactions. The first detects mass spectrometrically low level fluxes of reactively and nonreactively surface scattered species as a function of scattering angle and velocity while the second ultrahigh velocity (UHV) molecular beam, laser induced fluorescence apparatus is capable of measuring chemiluminescence produced by either gas phase or gas-surface interactions. A number of proposed experiments are described.

Cross, J. B.↗
NASA NTRS1978

Ultraviolet-photoproduced organic solids synthesized under simulated Jovian conditions - Molecular analysis

In an earlier paper, Khare and Sagan reported the production of a brownish polymeric material from the near-ultraviolet irradiation of simulated jovian atmospheres with a low hydrogen abundance. Examination of this product indicates that hydrogen sulfide is the initial photon acceptor; the powder resulting after extraction with benzene is 84 percent sulfur, largely S8. In results reported here, the remaining 16 percent was pyrolyzed and then examined by gas chromatography-mass spectrometry. Pyrolysis at 450 C yielded a series of alkanes, alkenes, C3-alkylbenzenes, aromatics, thiophenes, alkylthiophenes, alkylmercaptans, alkyldisulfides, together with the nitrogenous compounds hydrogen cyanide, methyl cyanide, alkylisothiocyanates, acrylonitrile, and allylisothiocyanates. Some of these compounds might be sought on Jupiter and Saturn and their satellites by remote infrared and ultraviolet spectroscopy and directly by entry probes.

Khare, B. N.↗
NASA NTRS1999

Design Through Simulation of a Molecular Sieve Column for Treatment of MON-3

The presence of water in propellant-grade MON-3 is a concern in the Aerospace Industry. NASA Johnson Space Center (JSC), White Sands Test Facility (WSTF) Propulsion Department has evaluated many types of molecular sieves for control of iron, the corrosion product of water in Mixed Oxides of Nitrogen (MON-3). In 1995, WSTF initiated laboratory and pilot-scale testing of molecular sieve type 3A for removal of water and iron. These tests showed sufficient promise that a series of continuous recycle tests were conducted at WSTF. Periodic samples of the circulating MON-3 solution were analyzed for water (wt %) and iron (ppm, wt). This test column was modeled as a series of transfer units; i. e., each unit represented the height equivalent of a theoretical plate. Such a model assumes there is equilibrium between the adsorbent material and the effluent stream from the unit. Operational and design parameters were derived based on the simulation results. These parameters were used to predict the design characteristics of a proposed molecular sieve column for removal of water and iron from MON-3 at the NASA Kennedy Space Center (KSC). In addition, these parameters were used to simulate a small, single-pass operation column at KSC currently used for treating MON-3. The results of this work indicated that molecular sieve type 3A in 1/16 in. diameter pellets, in a column 2.5 ft. in diameter, 18 ft. in height, and operated at 25 gpm is adequate for the required removal of water and iron from MON-3.

Swartz, A. Ben↗
NASA NTRS1991

Molecular dynamics and quasidynamics simulations of the annealing of bulk and near-surface interstitials formed in molecular-beam epitaxial Si due to low-energy particle bombardment during deposition

The relaxation, diffusion, and annihilation of split and hexagonal interstitials resulting from 10 eV Si irradiation of (2x1)-terminated Si(100) are investigated. Molecular dynamics and quasidynamics simulations, utilizing the Tersoff many-body potential are used in the investigation. The interstitials are created in layers two through six, and stable atomic configurations and total potential energies are derived as a function of site symmetry and layer depth. The interstitial Si atoms are allowed to diffuse, and the total potential energy changes are calculated. Lattice configurations along each path, as well as the starting configurations, are relaxed, and minimum energy diffusion paths are derived. The results show that the minimum energy paths are toward the surface and generally involved tetrahedral sites. The calculated interstitial migration activation energies are always less than 1.4 eV and are much lower in the near-surface region than in the bulk.

Kitabatake, M.↗