SEARCH · Engineering Papers

Results for “forcing fields”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1960

On Ion Clusters in the Interstellar Gas

In a recent paper V.I. Krassovsky (1958) predicts the occurrence of clusters of large numbers of atoms and molecules around ions in the interstellar gas. He then proposes a number of physicochemical processes that would be considerably enhanced by the high particle density in such clusters. In particular, he suggests that absorption by negative ions formed in the clusters would account for the interstellar extinction without any necessity for the presence of grains. Because of the important consequences that ion clusters could have, it is necessary to examine their occurrence more fully. This note re-examines the formation of ion clusters in space and shows that even ion-molecule pairs are essentially non-existent. Ion clusters have been considered by Bloom and Margenau (1952) from the same point of view as that used by Krassovsky, whose basic reference (Joffe and Semenov 1933) unfortunately is not available. A different approach has been used by Eyring, Hirschfelder, and Taylor (1936) following the methods of chemical equilibrium. Both the references cited here enable one to conclude that clustering is negligible. Therefore, the treatment of Eyring et al. is more appropriate than the method of Bloom and Margenau, which depends on the statistical equilibrium of an atmosphere in a force field.

Donn, Bertram↗

NASA NTRS2001

Modeling Modern Methane Emissions from Natural Wetlands: Model Description and Results - 1

Methane is an important greenhouse gas which contributes about 22 percent to the present greenhouse effect. Natural wetlands currently constitute the biggest methane source and were the major source in preindustrial times. Wetland emissions depend highly on the climate, i.e., on soil temperature and water table. To investigate the response of methane emissions from natural wetlands to climate variations, a process-based model that derives methane emissions from natural wetlands as a function of soil temperature, water table, and net primary productivity is used. For its application on the global scale, global data sets for all model parameters are generated. In addition, a simple hydrologic model is developed in order to simulate the position of the water table in wetlands. The hydrologic model is tested against data from different wetland sites, and the sensitivity of the hydrologic model to changes in precipitation is examined. The global methane hydrology model constitutes a tool to study temporal and spatial variations in methane emissions from natural wetlands. The model is applied using high-frequency atmospheric forcing fields from European Center for Medium-range Weather Forecasts (ECMWF) re-analyses of the period from 1982 to 1993. We calculate global annual methane emissions from wetlands to be 260 teragrams per year. Twenty-five percent of these methane emissions originate from wetlands north of 30 degrees North Latitude. Only 60 percent of the produced methane is emitted, while the rest is re-oxidized. A comparison of zonal integrals of simulated global wetland emissions and results obtained by an inverse modeling approach shows good agreement. In a test with data from two wetlands the seasonality of simulated and observed methane emissions agrees well.

Biogeochemical processes↗

NASA NTRS2016

Li-Doped Ionic Liquid Electrolytes: From Bulk Phase to Interfacial Behavior

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.

molecular dynamics↗

NASA NTRS2016

Molecular Dynamics Modeling of Carbon Nanotube Composite Fracture Using ReaxFF

Carbon nanotube (CNT) fiber reinforced composites with specific tensile strengths and moduli approaching those of aerospace grade carbon fiber composites have recently been reported. This achievement was enabled by the emerging availability of high N/tex yarns in kilometer-scale quantities. While the production of this yarn is an impressive advance, its strength is still much lower than that of the individual CNTs comprising the yarn. Closing this gap requires understanding load transfer between CNTs at the nanometer dimensional scale. This work uses reactive molecular dynamics simulations to gain an understanding at the nanometer scale of the key factors that determine CNT nanocomposite mechanical performance, and to place more realistic upper bounds on the target properties. While molecular dynamics simulations using conventional force fields can predict elastic properties, the ReaxFF reactive forcefield can also model fracture behavior because of its ability to accurately describe bond breaking and formation during a simulation. The upper and lower bounds of CNT composite properties are investigated by comparing systems composed of CNTs continuously connected across the periodic boundary with systems composed of finite length CNTs. These lengths, effectively infinite for the continuous tubes and an aspect ratio of 13 for the finite length case, result from practical limitations on the number of atoms that can be included in a simulation. Experimentally measured aspect ratios are typically on the order of 100,000, so the calculated results should represent upper and lower limits on experimental mechanical properties. Finally, the effect of various degrees of covalent crosslinking between the CNTs and amorphous carbon matrix is considered to identify the amount of CNT-matrix covalent bonding that maximizes overall composite properties.

Jensen, Benjamin D.↗

NASA NTRS2017

The dependence of the oceans MOC on mesoscale eddy diffusivities: A model study

The dependence of the depth and strength of the ocean's global meridional overturning cells (MOC) on the specification of mesoscale eddy diffusivity (K) is explored in two ocean models. The GISS and MIT ocean models are driven by the same prescribed forcing fields, configured in similar ways, spun up to equilibrium for a range of K 's and the resulting MOCs mapped and documented. Scaling laws implicit in modern theories of the MOC are used to rationalize the results. In all calculations the K used in the computation of eddy-induced circulation and that used in the representation of eddy stirring along neutral surfaces, is set to the same value but is changed across experiments. We are able to connect changes in the strength and depth of the Atlantic MOC, the southern ocean upwelling MOC, and the deep cell emanating from Antarctica, to changes in K.

Global ocean↗

NASA NTRS2017

Small Ion Channel Linking Molecular Simulations and Electrophysiology

Ion channels are pore-forming protein assemblies that mediate the transport of small ions across cell membranes. Otherwise, membrane bilayers would be almost impermeable to ions incapable to traverse the low dielectric constant, hydrophobic membrane core. Ion channels are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses, cardiac functions, muscle contraction and apoptosis. On the other extreme of biological complexity, viral ion channels (viroporins) influence many stages of the virus infection cycle either through regulating virus replication, such as entry, assembly and release or modulating the electrochemical balance in the subcellular compartments of host cells. Ion channels were crucial components of protocells. Their emergence facilitated adaptation of nascent life to different environmental conditions. The earliest ion channels must have been much simpler than most of their modern ancestors. Viral channels are among only a few naturally occurring models to study the structure, function and evolution of primordial channels. Experimental studies of these properties are difficult and often unreliable. In principle, computational methods, and molecular dynamics (MD) simulations in particular, can aid in providing information about both the structure and the function of ion channels. However, MD suffers from its own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. It is, therefore, essential to determine the reliability of MD simulations. We propose to do so on the basis of two criteria. One is channel stability on time scales that extend for several microseconds or longer. The other is the ability to reproduce the measured ionic conductance as a function of applied voltage. If both the stability and the calculated ionic conductance are satisfactory it will greatly increase our confidence that the structure and the function of a channel are described sufficiently accurately. To our knowledge, long time scale stability (approx.10 micro-sec) and the correct electrophysiology have been shown so far for only one channel - the synthetic LS3 hexamer). In this presentation, this approach will be discussed in application to two viral channels - Vpu, encoded by the HIV-1 genome and p7 of hepatitis C.

Pohorille, Andrzej↗

NASA NTRS2016

Interfacial Structure and Capacitance of Li-doped Ionic Liquid Electrolytes from Molecular Simulations

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries, supercapacitors, and hybrid energy storage devices. Though Li-salt is often present in these systems, its influence on interfacial properties is largely uncharacterized. We, thereby, present an extensive computational analysis, supported by experimental comparisons, of the properties of a representative set of these electrolytesat an ideal carbon interface as a function of Li-salt doping and voltage. We have performed polarizable molecular (MD) dynamics simulations, using the APPLEP force field, to evaluate electric double layer (EDL) capacitance and distribution of Li+ in the EDL. Differential capacitance exhibits the characteristic camel profile and is insensitive to Li-doping. Li+ localizes in the second molecular layer of the EDL, which is a result of confinement from free energy barriers associated with ion layering. Joint MDelectronic structure computations show the electrochemical window of the electrolytes to be a weak function of Li-doping. Estimates of supercapacitor specific energy are made using the computed window and capacitance. The magnitude and trends in specific energy are in good agreement with experiment.

batteries↗

NASA NTRS2015

The Structure, Anharmonic Vibrational Frequencies, and Intensities of NNHNN+

A semi-global potential energy surface (PES) and quartic force field (QFF) based on fitting high-level electronic structure energies are presented to describe the structures and spectroscopic properties of NNHNN+. The equilibrium structure of NNHNN+ is linear with the proton equidistant between the two nitrogen groups and thus of D(sub ∞h) symmetry. Vibrational second-order perturbation theory (VPT2) calculations based on the QFF fails to describe the proton “rattle” motion, i.e., the antisymmetric proton stretch, due to the very flat nature of PES around the global minimum, but performs properly for other modes with sharper potential wells. Vibrational self-consistent field/virtual state configuration interaction (VSCF/VCI) calculations using a version of MULTIMODE without angular momentum terms successfully describe this motion and predict the fundamental to be at 759 cm(exp -1). This is in good agreement with the value of 746 cm(exp -1) from a fixed-node diffusion Monte Carlo calculation and the experimental Ar-tagged result of 743 cm(exp -1). Other VSCF/VCI energies are in good agreement with other experimentally reported ones. Both double-harmonic intensity and rigorous MULTIMODE intensity calculations show the proton transfer fundamental has a very strong intensity.

Astrochemistry↗

NASA NTRS2016

Computational and Experimental Study of Li-doped Ionic Liquids at Electrified Interfaces

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Psi). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Psi. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Psi-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Psi and suggest a source of impedance for Li+ to reach the electrode.

Haskins, Justin B.↗

NASA NTRS2017

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗

NASA NTRS2020

Computationally Efficient Modeling and Data Assimilation of Near-Surface Variability

Near-surface (< 20m) ocean exhibits high variability due to coupled interactions, for e.g., with the atmosphere, sea ice, land, etc. Here we focus on atmospheric heat and momentum (wind) forcing, which are known to cause diurnal variability within the mixed layer. Only recently with a combination of sufficiently high vertical/horizontal resolution (75L, 1/4deg) and sub-daily atmospheric forcing fields, ocean models are starting to resolve this diurnal variability. However, the computation expense of such a high vertical resolution is burdensome in the context of coupled modeling and data assimilation. An alternative approach is to parameterize this diurnal variability with a prognostic model, that is embedded into the ocean model.In the first part of this presentation, we will demonstrate results with the above two approaches, by comparing them to profiles of near-surface temperature and salinity. In the context of data assimilation and reanalysis, this modeling capability opens the door to re-examine and perhaps improve specification of background (or, ensemble) error characteristics. The second half of this talk will focus on illustrating diurnally varying errors within an ensemble DA, and possible approaches to improve localization (horizontal/vertical) to extract maximum possible observational information content from in-situ and satellite observations of sea surface temperature.

Akella, Santha↗

NASA NTRSDate not supplied

REACTER: A Heuristic Method for Reactive Molecular Dynamics

REACTER (www.reacter.org) is a heuristic protocol that enables the simulation of complex reactions using atomistic molecular dynamics (MD) with a fixed-valence force field. Incorporating reactions into classical MD with this approach allows modeling of reactive systems over greatly-increased time scales, enabling systems to be modeled with MD that would not otherwise be feasible. One or more competing multi-step reactions or series of reactions can be invoked simultaneously. Special treatment can be applied to neighboring atoms to relax high energy configurations while the simulation progresses. The original version of REACTER, which was implemented in the open-source LAMMPS simulation package as fix bond/react, was only available for serial simulations. This work describes the expansion of fix bond/react for use in parallel simulations, as well as the addition of various new options, including deletion of reaction by-products, reversible reactions, and custom reaction constraints. These new capabilities are demonstrated through large-scale simulations (200,000+ atoms) of the polymerization of polystyrene and nylon 6,6. The morphologies of both polymers are analyzed after reaching >99% extent of polymerization. Finally, the newly-added reversible reactions feature is demonstrated by rupturing these highly-entangled systems under uniaxial strain by defining a chain scission reaction.

polymer simulations, molecular dynamics↗

NASA NTRS2020

Lambert’s Problem – A Geometric Approach

A fundamental problem in spacecraft mission design is to find a free flight path from one place to another in a chosen travel time. Lambert studied this problem for free flight in an inverse square central force field, and Lagrange produced a solution in 1778. Although this is an old problem, a new approach may be of some value. There are two steps to the new solution. First, find every ellipse with focus at the origin that intersects two circles both centered at the origin. Then select only those ellipses that have a specified angle between the intersection points of the ellipse with the inner and outer circles. Second, find the travel times between intersections for each ellipse. These are all the possible travel times from which one may choose. Standard solutions to Lambert’s problem can be found in references [1-4].

Easton, Robert↗

NASA NTRSDate not supplied

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↗

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

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗

NASA NTRSDate not supplied

Hierarchical screening for Li-based solid electrolytes using fast, interpretable machine-learned potentials

Li-based solid-state electrolyte materials enable safer, all-solid-state batteries but the computational search for candidates with favorable stability and Li-ion conductivity is challenging due to the size of the search space and the cost of evaluating transport properties with ab initio methods. The prohibitive cost of high-throughput screening with DFT has lead to the development of surrogate models using geometric analysis, empirical potentials, and descriptors for ionic transport. Here, I will discuss a hierarchical screening approach for identifying promising materials using a combination of density functional theory, bond-valence methods, and machine learning potentials generated with the Ultra-Fast Force Fields (UF3) framework. We show how the inexpensive bond-valence method can be used to guide the generation of training samples for machine learning, in addition to filtering candidates. Finally, we apply the hierarchical workflow to screen for ionic conductivity across a database of Li-containing compounds.

Materials discovery↗

NASA NTRSDate not supplied

Hierarchical Screening for Li-Based Solid Electrolytes Using Fast, Interpretable Machine-Learned Potentials

Materials discovery↗