Modeling Nonequilibrium Reactive Systems with REACTER
Explore the source record for details and available documents.
machine learning, polymer simulations, molecular d↗
SEARCH · Engineering Papers
Results for “Reactive 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.
Lunar and Planetary Science XXXVI, Part 9
The following topics were discussed: Monitoring floods with NASA's ST6 autonomous spacecraft experiment; Dynamical cloud models constrained by high resolution spectroscopy of zodiacal light; The oxygen isotopic composition of the sun and implications for oxygen processing in molecular clouds; A nochian/hisperian hiatus and erosive reactivation of martian valley networks; Hard x-ray spectro-microscopy techniques; Thermoluminescence studies of carbonaceous chondrites, etc.
Source record↗
First Principles Based Reactive Atomistic Simulations to Understand the Effects of Molecular Hypervelocity Impact on Cassini's Ion and Neutral Mass Spectrometer
We report here on the predicted impact of species such as ice-water, CO2, CH4, and NH3, on oxidized titanium, as well as HC species on diamond surfaces. These simulations provide the dynamics of product distributions during and after a hypervelocity impact event, ionization fractions, and dissociation probabilities for the various species of interest as a function of impact velocity (energy). We are using these results to determine the relevance of the fragmentation process to Cassini INMS results, and to quantify its effects on the observed spectra.
hypervelocity collisions↗
Isomerization reaction dynamics and equilibrium at the liquid-vapor interface of water. A molecular-dynamics study
The gauche-trans isomerization reaction of 1,2-dichloroethane at the liquid-vapor interface of water is studied using molecular-dynamics computer simulations. The solvent bulk and surface effects on the torsional potential of mean force and on barrier recrossing dynamics are computed. The isomerization reaction involves a large change in the electric dipole moment, and as a result the trans/gauche ratio is considerably affected by the transition from the bulk solvent to the surface. Reactive flux correlation function calculations of the reaction rate reveal that deviation from the transition-state theory due to barrier recrossing is greater at the surface than in the bulk water. This suggests that the system exhibits non-Rice-Ramsperger-Kassel-Marcus behavior due to the weak solvent-solute coupling at the water liquid-vapor interface.
Benjamin, Ilan↗
Time-dependent Computational Studies of Premixed Flames in Microgravity
This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.
Kailasanath, K.↗
Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators
Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.
K Swaminathan Gopalan↗
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.↗
Spectral line shapes. Volume 4 - Proceedings of the Eighth International Conference, College of William and Mary, Williamsburg, VA, June 9-13, 1986
The present conference discusses topics in line broadening and shift phenomena in low-to-moderately dense plasmas, resonances in the continuum, hot and dense laser-produced plasmas, collisional broadening and intermolecular potentials, molecular broadening and collision-induced spectroscopy, collisional redistribution and laser-assisted collision processes, coherent and nonlinear spectroscopy, astrophysical line-shapes, and short-wavelength lasers. Attention is given to the widths and shifts of Be ion lines, the Stark broadening of Xe II and Fe I lines, perturbed and unperturbed series in the line-shapes of autoionizing Rydberg series, the line-broadening of hydrogenic ions in dense plasmas, Doppler-free line-shapes, the Franck-Condon principle, the dynamics of laser-assisted reactive collisions, the gas-kinetic effects of light, collision-induced processes in nonlinear optics, and line-shaping phenomena in cool stars.
Exton, Reginald J.↗
Time-dependent computational studies of flames in microgravity
The research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of the NASA Microgravity Science and Applications Program is described. The primary focus was on investigating fundamental questions concerning the propagation and extinction of premixed flames in Earth gravity and in microgravity environments. The approach was to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. A combination of one-dimensional and two-dimensional simulations was used to investigate the effects of curvature and dilution on ignition and propagation of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the mechanism leading to cellular instability, and to study the effects of gravity on the transition to cellular structure. A flame in a microgravity environment can be extinguished without external losses, and the mechanism leading to cellular structure is not preferential diffusion but a thermo-diffusive instability. The simulations have also lead to a better understanding of the interactions between buoyancy forces and the processes leading to thermo-diffusive instability.
Oran, Elaine S.↗
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↗
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↗
Analysis of leading edge and trailing edge cover glass samples before and after treatment with advanced satellite contamination removal techniques
Two samples from Long Duration Exposure Facility (LDEF) experiment M0003-4 were analyzed for molecular and particulate contamination prior to and following treatment with advanced satellite contamination removal techniques (CO2 gas/solid jet spray and oxygen ion beam). The pre- and post-cleaning measurements and analyses are presented. The jet spray removed particulates in seconds. The low energy reactive oxygen ion beam removed 5,000 A of photo polymerized organic hydrocarbon contamination in less than 1 hour. Spectroscopic analytical techniques were applied to the analysis of cleaning efficiency including: Fourier transform infrared, Auger, x ray photoemissions, energy dispersive x ray, and ultraviolet/visible. The results of this work suggest that the contamination studied here was due to spacecraft self-contamination enhanced by atomic oxygen plasma dynamics and solar UV radiation. These results also suggest the efficacy for the jet spray and ion beam contamination control technologies for spacecraft optical surfaces.
Hotaling, S. P.↗
Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels
This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.
Pellett, Gerald L.↗
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↗
Reactive control and reasoning assistance for scientific laboratory instruments
Scientific laboratory instruments that are involved in chemical or physical sample identification frequently require substantial human preparation, attention, and interactive control during their operation. Successful real-time analysis of incoming data that supports such interactive control requires: (1) a clear recognition of variance of the data from expected results; and (2) rapid diagnosis of possible alternative hypotheses which might explain the variance. Such analysis then aids in decisions about modifying the experiment protocol, as well as being a goal itself. This paper reports on a collaborative project at the NASA Ames Research Center between artificial intelligence researchers and planetary microbial ecologists. Our team is currently engaged in developing software that autonomously controls science laboratory instruments and that provides data analysis of the real-time data in support of dynamic refinement of the experiment control. the first two instruments to which this technology has been applied are a differential thermal analyzer (DTA) and a gas chromatograph (GC). coupled together, they form a new geochemicstry and microbial analysis tool that is capable of rapid identification of the organiz and mineralogical constituents in soils. The thermal decomposition of the minerals and organics, and the attendance release of evolved gases, provides data about the structural and molecular chemistry of the soil samples.
Thompson, David E.↗
Radiolytic Cryovolcanism Revisited
Location of Enceladus within the inner magnetosphere of Saturn makes it likely that energetic electron irradiation, ion-neutral chemistry of the emergent cryovolcanic plume gas, and resultant radiolytic modification of surface ice composition could impact long-term evolution of molecular composition in the global ice crust. Thermally-driven convection of ice in the south polar terrain would bring radiolytic oxidant products into contact with subsurface reservoirs of primordial organics on million-year time scales. The chemical reactions leading to CO2 gas production as a primary gas-piston driver of cryovolcanic activity would be exponentially elevated in the warm-ice margins of a heated fluid. The fluid temperature is typically assumed to be below 273 K but could be much higher in a gas-pressurized deep subsurface environment. The Perrier Ocean model has demonstrated how a CO2-loaded fluid could account for the observed jets, while the heat content of the fluid arising from the moon's deeper interior could support high levels of chemical reactivity in the thermal margins. Since mass loss and tidal dissipation arguments do not support continuous activity over billions of years, the activity is likely very episodic so that even low-level energy sources including irradiation-driven radiolysis of surface ices could substantially contribute to the chemical dynamics of the activity apparently now in high phase. A multi-phase thermochemical model, and supporting laboratory measurements of temperature-dependent reaction rates, are needed to investigate these potentially complex processes.
Cooper, John F.↗
Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations
A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.
Changwoon Jang↗
The dynamics of herpesvirus and polyomavirus reactivation and shedding in healthy adults: a 14-month longitudinal study
Humans are infected with viruses that establish long-term persistent infections. To address whether immunocompetent individuals control virus reactivation globally or independently and to identify patterns of sporadic reactivation, we monitored herpesviruses and polyomaviruses in 30 adults, over 14 months. Epstein-Barr virus (EBV) DNA was quantitated in saliva and peripheral blood mononuclear cells (PBMCs), cytomegalovirus (CMV) was assayed in urine, and JC virus (JCV) and BK virus (BKV) DNAs were assayed in urine and PBMCs. All individuals shed EBV in saliva, whereas 67% had >or=1 blood sample positive for EBV. Levels of EBV varied widely. CMV shedding occurred infrequently but occurred more commonly in younger individuals (P<.03). JCV and BKV virurias were 46.7% and 0%, respectively. JCV shedding was age dependent and occurred commonly in individuals >or=40 years old (P<.03). Seasonal variation was observed in shedding of EBV and JCV, but there was no correlation among shedding of EBV, CMV, and JCV (P>.50). Thus, adults independently control persistent viruses, which display discordant, sporadic reactivations.
Non-NASA Center↗