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 235 records · Page 13

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

The Earliest Ion Channels in Protocellular Membranes

Cellular membranes with their hydrophobic interior are virtually impermeable to ions. Bulk of ion transport through them is enabled through ion channels. Ion channels of contemporary cells are complex protein molecules which span the membrane creating a cylindrical pore filled with water. Protocells, which are widely regarded as precursors to modern cells, had similarly impermeable membranes, but the set of proteins in their disposal was much simpler and more limited. We have been, therefore, exploring an idea that the first ion channels in protocellular membranes were formed by much smaller peptide molecules that could spontaneously selfassemble into short-lived cylindrical bundles in a membrane. Earlier studies have shown that a group of peptides known as peptaibols is capable of forming ion channels in lipid bilayers when they are exposed to an electric field. Peptaibols are small, non-genetically encoded peptides produced by some fungi as a part of their system of defense against bacteria. They are usually only 14-20 residues long, which is just enough to span the membrane. Their sequence is characterized by the presence of non-standard amino acids which, interestingly, are also expected to have existed on the early earth. In particular, the presence of 2-aminoisobutyric acid (AIB) gives peptaibols strong helix forming propensities. Association of the helices inside membranes leads to the formation of cylindrical bundles, typically containing 4 to 10 monomers. Although peptaibols are excellent candidates for models of the earliest ion channels their structures, which are stabilized only by van der Waals forces and occasional hydrogen bonds between neighboring helices, are not very stable. Although it might properly reflect protobiological reality, it is also a major obstacle in studying channel behavior. For this reason we focused on two members of the peptaibol family, trichotoxin and antiamoebin, which are characterized by a single conductance level. This indicates that their structures are unique and stable. In addition, it is also believed that the trichotoxin channel displays some selectivity between potassium and chloride ions. This makes trichotoxin and antiamoebin ideal models of the earliest ion channels that could provide insight into the origins of ion conductance and selectivity. In the absence of crystal structure of the trichotoxin and antiamoebin channels, we propose their molecular models based on experimentally determined number of monomers forming the bundles. We use molecular dynamics simulations to validate the models in terms of their conductance and selectivity. On the basis of our simulations we show that the emergence of channels built of small, alpha-helical peptides was protobiologically plausible and did not require highly specific amino acid sequences, which is a convenient evolutionary trait. Despite their simple structure, such channels could possess properties that, at the first sight, appear to require markedly larger complexity. To this end, we will discuss how the amino acid sequence and structure of primitive channels give rise to the phenomena of ionic conductance and selectivity across the earliest cell walls, which were essential functions for the emergence and early evolution of protocells. Furthermore, we will argue that even though architectures of membrane proteins are not nearly as diverse as those of water-soluble proteins, they are sufficiently flexible to adapt readily to the functional demands arising during evolution.

Mijajlovic, Milan↗

Sputtering from a Porous Material by Penetrating Ions

Porous materials are ubiquitous in the universe and weathering of porous surfaces plays an important role in the evolution of planetary and interstellar materials. Sputtering of porous solids in particular can influence atmosphere formation, surface reflectivity, and the production of the ambient gas around materials in space, Several previous studies and models have shown a large reduction in the sputtering of a porous solid compared to the sputtering of the non-porous solid. Using molecular dynamics simulations we study the sputtering of a nanoporous solid with 55% of the solid density. We calculate the electronic sputtering induced by a fast, penetrating ion, using a thermal spike representation of the deposited energy. We find that sputtering for this porous solid is, surprisingly, the same as that for a full-density solid, even though the sticking coefficient is high.

Rodriguez-Nieva, J. F.↗

Formation, Structure and Properties of Amorphous Carbon Char from Polymer Materials in Extreme Atmospheric Reentry Environments

Amorphous carbonaceous char produced from the pyrolysis of polymer solids has many desirable properties for ablative heat shields for space vehicles. Molecular dynamics simulations are presented to study the transformation of the local atomic structure from virgin polymer to a dense, disordered char [1]. Release of polymer hydrogen is found to be critical to allow the system to collapse into a highly coordinated char structure. Mechanisms of the char formation process and the morphology of the resulting structures are elucidated. Thermal conductivity and mechanical response of the resulting char are evaluated [2]. During reenty, the optical response and oxidative reactivity of char are also important properties. Results of ab initio computations of char optical functions [3] and char reactivity [4] are also presented.

Lawson, John W.↗

The Effect of Surface Chemical Functionality Upon Ice Adhesion

In nature, anti-freeze proteins present in fish utilize specific organic functionalities to disrupt ice crystal formation and propagation. Based on these structures, surfaces with controlled chemical functionality and chain length were evaluated both experimentally and computationally to assess the effect of both parameters in mitigating ice formation. Linear aliphatic dimethylethoxysilanes terminated with methyl or hydroxyl groups were prepared, characterized, and used to coat aluminum. The effect upon icing using a microdroplet freezing apparatus and the Adverse Environment Rotor Test Stand found hydroxyl-terminated materials exhibited a greater propensity for ice formation and adhesion. Molecular dynamics simulations of a silica substrate bearing functionalized species of similar composition were brought into contact with a pre-equilibrated ice crystal. Several parameters including chain mobility were monitored to ascertain the size of a quasi-liquid layer. The studies suggested that chain mobility affected the interface between ice and the surface more than terminal group chemical composition.

Smith, Joseph G., Jr.↗

Multiscale Modeling of Carbon/Phenolic Composite Thermal Protection Materials: Atomistic to Effective Properties

Next generation ablative thermal protection systems are expected to consist of 3D woven composite architectures. It is well known that composites can be tailored to achieve desired mechanical and thermal properties in various directions and thus can be made fit-for-purpose if the proper combination of constituent materials and microstructures can be realized. In the present work, the first, multiscale, atomistically-informed, computational analysis of mechanical and thermal properties of a present day - Carbon/Phenolic composite Thermal Protection System (TPS) material is conducted. Model results are compared to measured in-plane and out-of-plane mechanical and thermal properties to validate the computational approach. Results indicate that given sufficient microstructural fidelity, along with lowerscale, constituent properties derived from molecular dynamics simulations, accurate composite level (effective) thermo-elastic properties can be obtained. This suggests that next generation TPS properties can be accurately estimated via atomistically informed multiscale analysis.

Materials Science↗

Monte Carlo Simulation of a Knudsen Effusion Mass Spectrometer Sampling System

Knudsen flow is easily simulated with a Monte Carlo method. In this study we develop a Visual Basic for Excel (VBA) code to simulate the molecular beam from a vaporizing solid. The system at NASA Glenn uses the restricted collimation method of Chatillon and colleagues, which consists of two apertures between the effusion cell and the ionizer. The diameter of the first aperture is smaller than the diameter of the effusion cell orifice, so the ionizer effectively only sees inside the effusion cell. The code is able to calculate the transmission coefficient through the cell orifice, through the cell orifice and the first aperture, and through the cell orifice and first and second aperatures. Calculated transmission coefficients through the cell orifice are compared to tabulated values to validate the code. Then transmission coefficients are calculated through the cell orifice and both apertures to the ionizer. This allows the geometry (aperture spacing and diameters) of the sampling system the sampling system to be optimized. Calculated transmission factors are also compared to literature values calculated via an analytic method.

Radke, Michael J.↗

Toward Ultralight High Strength Structural Materials via Collapsed Carbon Nanotube Bonding

The growing commercial availability of carbon nanotube (CNT) macro-assemblies such as sheet and yarn is making their use in structural composite components increasingly feasible. However, the mechanical properties of these materials continue to trail those of state-of-the-art carbon fiber composites due to relatively weak inter-tube load transfer. Forming covalent links between adjacent CNTs promises to mitigate this problem, but it has proven difficult in practice to introduce them chemically within densified and aligned CNT materials due to their low permeability. To avoid this limitation, this work explores the combination of pulsed electrical current, temperature, and pressure to introduce inter-CNT bonds. Reactive molecular dynamics simulations identify the most probable locations, configurations, and conditions for inter-nanotube bonds to form. This process is shown to introduce covalent linkages within the CNT material that manifest as improved macroscale mechanical properties. The magnitude of this effect increases with increasing levels of prealignment of the CNT material, promising a new synthesis pathway to ultralight structural materials with specific strengths and stiffnesses exceeding 1 and 100 GPa/(g/cu.cm), respectively.

Jensen, Benjamin D.↗

A Multiscale Cohesive Law for Carbon Fiber Networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. Here, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

W. SCHILL↗

Free Energy Landscapes for Elucidating the Structural Consequences of Exon-20 mutations on the ErbB Family of Protein Kinases

The ErbB family of protein kinases plays an important role in major cellular functions and consequently mutations in the functional regions of these proteins are implicated in several types of cancer growths. To envision rational design of small molecule drugs that target the diseased proteins it is important to quantify the structural effects of the mutations, as some of these mutants render the protein resistant to tyrosine kinase inhibitors (TKIs). Herein we use advanced sampling techniques and long-timescale molecular dynamics simulations to predict the effect of major exon 20 mutations on the ErbB family, specifically EGFR and HER2 proteins. Exon 20 mutations have been clinically known to induce TKI resistance, though the mechanisms of such an effect is poorly understood. By mapping out the free energy landscape of the mutants and comparing them against the wild-type, we elucidate the structural differences in the binding pocket region that alter the nature of drug-protein interactions. We believe that these insights will play a pivotal role in developing small molecule drugs that overcome the TKI resistance.

Ashwin Ravichandran↗

Twin nucleation and growth mechanism in Ni-based superalloys

While micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 C, many aspects of twin nucleation and growth remain unexplored. Kolbe mechanism for micro-twinning, based on thermally activated reordering, is probably the only concept currently widely accepted by the scientific community. We propose a qualitatively different mechanism for nucleation and growth of twins. The mechanism can be briefly described as follows. Penetration of a gamma prime precipitate by two 1/2(110) edge dislocations travelling on adjacent {111} glide planes triggers nucleation (at the interface of the precipitate and the matrix) and emission of Shockley partial of screw character on the glide plane of the edge dislocation, which entered the precipitate first (generating trailing high-energy anti-phase boundary (APB)). Propagation of this Shockley partial into the precipitate converts the APB into super intrinsic stacking fault (SISF). The recurring arrival of additional edge dislocations on the glide planes adjacent to the configuration described above leads to formation of super extrinsic stacking fault SESF, subsequent micro-twin formation, and growth of the twinned region. We demonstrate the proposed mechanism via molecular dynamics simulations.

Ni-superalloys↗

Predicting Char Yield of High-Temperature Resins

A simulation protocol is developed to predict the char yield of organic resins during high-temperature processing. Such in silico methods can help screen promising new formulations for advanced materials, but previously no chemistry-sensitive technique existed to predict the important experimental value of char yield. The method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to processing temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small byproduct molecules to mimic outgassing, and the application of high pressure to eliminate porosity and encourage graphitization. More than ten different resin chemistries are investigated, including arylacetylenes, cyanate esters, phthalonitriles, and polyimides, representing a diverse group of precursors with respect to initial cyclic content, heteroatoms and reactive groups. The protocol correctly predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. The properties of the resins during processing are compared, including outgassing products, morphology of the final chemical configurations, cyclic content and mechanical properties.

polymer simulations, molecular dynamics↗

Twin nucleation and growth mechanism in Ni-based superalloys

While micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 C, many aspects of twin nucleation and growth remain unexplored. Kolbe mechanism for micro-twinning, based on thermally activated reordering, is probably the only concept currently widely accepted by the scientific community. We propose a qualitatively different mechanism for nucleation and growth of twins. The mechanism can be briefly described as follows. Penetration of a   precipitate by two 1/2(110) edge dislocations travelling on adjacent {111} glide planes triggers nucleation (at the interface of the precipitate and the matrix) and emission of Shockley partial of screw character on the glide plane of the edge dislocation, which entered the precipitate first (generating trailing high-energy anti-phase boundary (APB)). Propagation of this Shockley partial into the precipitate converts the APB into super intrinsic stacking fault (SISF). The recurring arrival of additional edge dislocations on the glide planes adjacent to the configuration described above leads to formation of super extrinsic stacking fault SESF, subsequent micro-twin formation, and growth of the twinned region. We demonstrate the proposed mechanism via molecular dynamics simulations.

Ni-based superalloys↗

Predicting Char Yield of High-Temperature Resins

A simulation technique has been developed for predicting the char yield of organic resins during high-temperature processing. In silico methods can aid in the screening of new advanced materials for a number of important properties, but no chemistry-sensitive protocol currently exists for predicting the important experimental value of char yield. The proposed method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to pyrolysis temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small by-product molecules to mimic outgassing, and the application of high pressure to encourage the formation of a dense, glassy network. Six different resin chemistries were investigated: an ethynyl, a phenylethynyl, a cyanate ester, a phthalonitrile, acrylonitrile, and adamantane. These candidates represent a diverse group of precursors with respect to initial cyclic content, presence of heteroatoms, and types of reactive groups. The protocol developed accurately predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. Several simulated properties of the carbonized structures are compared with experimental results, including outgassing products, morphology of the final chemical configurations, cyclic content, and mechanical properties.

polymer simulations, molecular dynamics↗

Predicting Char Yield of High-Temperature Resins

A simulation protocol is developed to predict the char yield of organic resins during high-temperature processing. Simulation protocol in silico methods can help screen promising new formulations for advanced materials, but previously no chemistry-sensitive technique existed to predict the important experimental value of char yield. The method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to processing temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small byproduct molecules to mimic outgassing, and the application of high pressure to eliminate porosity and encourage graphitization. More than ten different resin chemistries are investigated, including arylacetylenes, cyanate esters, phthalonitriles, and polyimides, representing a diverse group of precursors with respect to initial cyclic content, heteroatoms and reactive groups. The protocol correctly predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. The properties of the resins during processing are compared, including outgassing products, morphology of the final chemical configurations, cyclic content and mechanical properties.

polymer simulations, molecular dynamics↗

Effect of Cr on Twinning in Ni-Based Superalloys

Micro-twinning is the major creep deformation mechanism in Ni-based superalloys at temperatures above 700 °C. Recent experiments indicate that formation and extension of micro-twins is reliant on the segregation of gamma phase formers (e.g., Co and Cr) along the twin interface in gamma prime phase. In this study we use molecular dynamics simulation to investigate the effect of site preference of Cr in Ni 3 Al gamma prime phase on deformation twinning. We demonstrate that, depending on the site preference (Cr on Ni sites, or Cr on Al sites), two drastically different deformation behaviors can be expected.

Ni-based superalloys↗