Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemical dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Shock Hugoniot measurements of single-crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) compressed to 83 GPa

We present laser-driven shock Hugoniot measurements of single-crystal (SC) 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) between 15 and 83 GPa, spanning pressures below and well above the Chapman–Jouguet pressure of ~28 GPa for TATB formulations (TATB grains mixed with plastic binders at 5–10 wt. %). The new SC data are generally ~3% more compressible than previously published data on neat and formulated TATB measured in gas-gun and explosive-driven experiments. An exception is at compressions in the density of ~1.5 (~30–40 GPa), where our new SC data exhibit significantly lower pressures than previous results on overdriven TATB formulations, suggesting that our SC samples remain largely unreacted below 35 GPa over the short nanosecond-time scales inherent to our laser-driven experiments. We report these novel equation-of-state measurements are a critical step toward understanding TATB in its most fundamental form and improving predictive modeling of TATB-based explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Comparison of two pressure–temperature equilibration methods

We compare and contrast the traditionally used method of solving the pressure–temperature equilibration problem in hydrodynamics, where specific internal energy and density are considered independent variables, with a different method where pressure and temperature are independent variables. With the goal of examining the robustness of the two methods as the number of components increases, we examine 2-, 4-, 6-, and 8-component systems. After equilibrating more than 10 4 initial conditions for each system using both methods, we demonstrate that the latter method constrains the search space by lowering its dimensionality and forces a better initial guess, resulting in a higher probability of convergence to solution with fewer, cheaper iterations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Nucleation kinetics model for primary crystallization in Al–Y–Fe metallic glass

The high density of aluminum nanocrystals (>10 21 m -3 ) that develop during the primary crystallization in Al-based metallic glasses indicates a high nucleation rate (~10 18 m -3 s -1 ). Several studies have been advanced to account for the primary crystallization behavior, but none have been developed to completely describe the reaction kinetics. Recently, structural analysis by fluctuation electron microscopy has demonstrated the presence of the Al-like medium range order (MRO) regions as a spatial heterogeneity in as-spun Al 88 Y 7 Fe 5 metallic glass that is representative for the class of Al-based amorphous alloys that develop Al nanocrystals during primary crystallization. From the structural characterization, an MRO seeded nucleation configuration is established, whereby the Al nanocrystals are catalyzed by the MRO core to decrease the nucleation barrier. The MRO seeded nucleation model and the kinetic data from the delay time (τ) measurement provide a full accounting of the evolution of the Al nanocrystal density (N v ) during the primary crystallization under isothermal annealing treatments. Moreover, the calculated values of the steady state nucleation rates (J ss ) predicted by the nucleation model agree with the experimental results. Moreover, the model satisfies constraints on the structural, thermodynamic, and kinetic parameters, such as the critical nucleus size, the interface energy, and the volume-free energy driving force that are essential for a fully self-consistent nucleation kinetics analysis. In conclusion, the nucleation kinetics model can be applied more broadly to materials that are characterized by the presence of spatial heterogeneities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

In situ thermal conductivity measurement revealing kinetics of thermochemical reactions

Utilizing thermochemical reactions for thermal energy storage and solar fuel production has been an emerging research topic. Thermal transport properties of the materials are an important parameter that can determine the kinetics and efficiency of thermochemical reactions. With the increasing number of new thermochemical materials (TCMs); however, there is a lack of reliable techniques to monitor the thermal transport property of the materials and their changes as a function of reactions in real time. In this work, we report the in situ monitoring of thermochemical reactions using modulated photothermal radiometry (MPR). The thermal conductivities of two TCMs, namely, calcium hydroxide (Ca(OH) 2 ) and Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585), were measured as a function of temperature and time using the MPR technique. The measured thermal conductivities were correlated to the reaction. The work has two significant contributions to the research communities. First, it provides a non-invasive diagnostic tool for monitoring the thermal transport properties of TCMs that can potentially be a high-throughput measurement technique conducive to optimizing TCMs, reactors, and related thermal systems. Second, for TCMs that show observable changes in thermal transport properties, a correlation between the measured thermal conductivity and the conversion fraction of the reaction can be established for monitoring the reaction kinetics based on thermal characterization.

14 SOLAR ENERGY↗
DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗
DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Synthetic soils for ecological and synthetic biology applications

Soils are heterogeneous and dynamic systems characterized by complex physical, chemical, and biological interactions. Understanding these interactions is critical, as they influence plant productivity, global biogeochemical cycles, and ecosystem resilience. While ecologists have long studied soils in field, greenhouse, and laboratory settings, their complexity and heterogeneity make it challenging to pinpoint key properties driving biological processes and derive mechanistic insights. Advancements in synthetic biology, which seeks to engineer and control biological processes in soils, have increased the demand for standardized and controllable experimental platforms. These platforms, referred to here as ‘synthetic soils’, are systems designed to reproduce selected physicochemical characteristics of natural soils in a simplified and defined format, allowing scientists to systematically change soil physicochemical properties (i.e. texture, mineralogy, pH) to study how biological components (i.e. microbes, plants, soil fauna, etc.) respond to, modify, or interact within these controlled environments. This review explores existing synthetic soils, their advantages, limitations, and applications in ecology and synthetic biology, and discusses potential directions for their future development.

artificial soils↗
DOE OSTI2020

Modified atomic layer deposition of MoS 2 thin films

As one of the most attractive transition metal dichalcogenides (TMDs), the growth of molybdenum disulfide (MoS 2 ) with industrial compatibility is of great importance. Atomic layer deposition (ALD) has been shown to be a promising method to achieve the growth of high-quality TMD materials. However, MoS 2 films deposited by ALD often are amorphous with nonideal stoichiometry and require high-temperature post-deposition annealing. Here, we introduce a modified ALD recipe using Mo(CO) 6 and H 2 S, resulting in controllable linear growth behavior, a S-to-Mo ratio of 2:1, and crystalline films at a temperature as low as 190 °C. The growth mechanisms and key factors leading to this improvement are proposed and complemented by kinetics calculations. This newly developed methodology relies on aligning the process time with the reaction kinetics of carbonyl disassociation. The MoS 2 films prepared herein were shown to be active hydrogen evolution reaction catalysts.

2D materials↗
DOE OSTI2025

Nonpyrophoric alternative to trimethylaluminum for the atomic layer deposition of Al 2 O 3

We have examined the atomic layer deposition (ALD) of Al 2 O 3 using a nonpyrophoric precursor, which possesses only Al–N bonds and no Al–C bonds: Al(N(CH 3 ) 2 ) 2 (–N(C 2 H 5 )–C 2 H 4 –N(C 2 H 5 ) 2 ), which we refer to as bis-dimethylamino-diamino-aluminum (BDMADA-Al). We employed a quartz crystal microbalance (QCM) to monitor ALD in situ and in real time, and the deposited thin films have been characterized using x-ray photoelectron spectroscopy, spectroscopic ellipsometry, x-ray reflectivity, and atomic force microscopy. Films deposited with BDMADA-Al result in near-stoichiometric Al2O3 at temperatures ranging from T = 120 to 285 °C using H 2 O as the coreactant and at T = 285 °C using t-BuOH as the coreactant. The properties of the films (density, C incorporation, stoichiometry, growth rates) are comparable using either BDMADA-Al or trimethylaluminum (TMA) as the precursor under similar reaction conditions. The N content in the films deposited with BDMADA-Al as the precursor is <1% at both T = 120 and 285 °C. Additionally, both BDMADA-Al|t-BuOH and TMA|t-BuOH processes deposit near-stoichiometric Al 2 O 3 at T = 285 °C, while neither process chemistry produces a constant growth rate per cycle at T = 120 °C. Close examination of single cycle QCM data from BDMADA-Al ALD indicates unique ligand exchange reaction mechanisms during each of the half-reactions. Experimentally, we have evidence supporting the mechanism in which upon the chemisorption of BDMADA-Al, the species on the surface retains the –N(CH 3 ) 2 ligands, while the diamine ligand is mostly lost on dissociative adsorption, except possibly at the higher temperature (285 °C). As a result, our work demonstrates that BDMADA-Al can be used as an effective alternative precursor to TMA for the ALD of Al 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic STEM-EELS for single-atom and defect measurement during electron beam transformations

This study introduces the integration of dynamic computer vision–enabled imaging with electron energy loss spectroscopy (EELS) in scanning transmission electron microscopy (STEM). This approach involves real-time discovery and analysis of atomic structures as they form, allowing us to observe the evolution of material properties at the atomic level, capturing transient states traditional techniques often miss. Rapid object detection and action system enhances the efficiency and accuracy of STEM-EELS by autonomously identifying and targeting only areas of interest. This machine learning (ML)–based approach differs from classical ML in that it must be executed on the fly, not using static data. We apply this technology to V-doped MoS 2 , uncovering insights into defect formation and evolution under electron beam exposure. This approach opens uncharted avenues for exploring and characterizing materials in dynamic states, offering a pathway to increase our understanding of dynamic phenomena in materials under thermal, chemical, and beam stimuli.

47 OTHER INSTRUMENTATION↗
NASA NTRS1980

Iron and titanium distribution on the moon from orbital gamma ray spectrometry with implications for crustal evolutionary models

A set of Fe and Ti maps and regional values are obtained from the Apollo 15 and 16 orbital gamma ray data by energy band analysis. High-Ti basalts predominate the early and late stages of mare volcanism with high-Fe basaltic volcanism in the interim. The first evidence of a high-Ti-KREEP basalt association is found in the Aristarchus region. A N-S asymmetry for Fe and Ti in the east limb and farside highlands complicates the E-W asymmetry for Th but substantiates crustal inhomogeneity. The observed crustal inhomogeneity adds an additional objection to the primitive source model for crustal evolution. The high-Ti-KREEP basalt association and the general trend of decreasing mare basalt Ti with time lend support to the cumulate source model; however, this model cannot account for young, high-Ti maria. The dynamic assimilation model better accounts for chemical variations observed on the moon.

Davis, P. A., Jr.↗
NASA NTRS1992

A new numerical model of the middle atmosphere. I - Dynamics and transport of tropospheric source gases

Attention is given to a new model of the middle atmosphere which includes, in addition to the equations governing the zonal mean state, a potential vorticity equation for a single planetary-scale Rossby wave, and an IR radiative transfer code for the stratosphere and lower mesosphere, which replaces the Newtonian cooling parameterization used previously. It is shown that explicit computation of the planetary-scale wave field yields a more realistic representation of the zonal mean dynamics and the distribution of trace chemical species. Wave breaking produces a well-mixed 'surf zone' equatorward of the polar night vortex and drives a meridional circulation with downwelling on the poleward side of the vortex. This combination of mixing and downwelling produces shallow meridional gradients of trace gases in the subtropics and middle latitudes, and very steep gradients at the edge of the polar vortex. Computed distributions of methane and nitrous oxide are shown to agree well with observations.

Garcia, Rolando R.↗
NASA NTRS1993

Upwind MacCormack Euler solver with non-equilibrium chemistry

A computer code, designated UMPIRE, is currently under development to solve the Euler equations in two dimensions with non-equilibrium chemistry. UMPIRE employs an explicit MacCormack algorithm with dissipation introduced via Roe's flux-difference split upwind method. The code also has the capability to employ a point-implicit methodology for flows where stiffness is introduced through the chemical source term. A technique consisting of diagonal sweeps across the computational domain from each corner is presented, which is used to reduce storage and execution requirements. Results depicting one dimensional shock tube flow for both calorically perfect gas and thermally perfect, dissociating nitrogen are presented to verify current capabilities of the program. Also, computational results from a chemical reactor vessel with no fluid dynamic effects are presented to check the chemistry capability and to verify the point implicit strategy.

Sherer, Scott E.↗
NASA NTRS1993

Interpretation of NO(x)/NO(y) observations from AASE-2 using a model of chemistry along trajectories

In situ measurements of NO and NO(y) are used to derive the ratio NO(x)/NO(y) along the flight track of the NASA ER-2 aircraft. Data are presented for two flights at midlatitudes in October 1991 during the Airborne Arctic Stratospheric Expedition-2 (AASE-2). Aerosol particle surface area was concurrently measured. The observations are compared with a photochemical model integrated along back trajectories from the aircraft flight track. Comparison of observations with the model run along trajectories and at a fixed position clearly and quantitatively demonstrates the importance of an air parcel's dynamic history in interpretation of local chemical observations. Comparison of the data with model runs under different assumptions regarding heterogeneous chemistry further reinforces the case for occurrence of the reaction of N2O5 + H2O on sulfate aerosol surfaces in the atmosphere. Finally, comparisons for which relative changes in the model and the data are not consistent caution that our ability to resolve all the observations is not yet complete.

Kawa, S. R.↗
NASA NTRS1995

Galacti chemical evolution: Hygrogen through zinc

Using the output from a grid of 60 Type II supernova models (Woosley & Weaver 1995) of varying mass (11 approx. less than (M/solar mass) approx. less than 40) and metallicity (0, 10(exp -4), 0.01, and 1 solar metallicity), the chemical evolution of 76 stable isotopes, from hydrogen to zinc, is calculated. The chemical evolution calculation employs a simple dynamical model for the Galaxy (infall with a 4 Gyr e-folding timescale onto a exponential dsk and 1/r(exp 2) bulge), and standard evolution parameters, such as a Salpeter initial mass function and a quadratic Schmidt star formation rate. The theoretical results are compared in detail with observed stellar abundances in stars with metallicities in the range -3.0 approx. less than (Fe/H) approx. less than 0.0 dex. While our discussion focuses on the solar neighborhood where there are the most observations, the supernova rates, an intrinsically Galactic quality, are also discussed.

Timmes, F. X.↗
NASA NTRS1997

Automated Simplification of Full Chemical Mechanisms

A code has been developed to automatically simplify full chemical mechanisms. The method employed is based on the Intrinsic Low Dimensional Manifold (ILDM) method of Maas and Pope. The ILDM method is a dynamical systems approach to the simplification of large chemical kinetic mechanisms. By identifying low-dimensional attracting manifolds, the method allows complex full mechanisms to be parameterized by just a few variables; in effect, generating reduced chemical mechanisms by an automatic procedure. These resulting mechanisms however, still retain all the species used in the full mechanism. Full and skeletal mechanisms for various fuels are simplified to a two dimensional manifold, and the resulting mechanisms are found to compare well with the full mechanisms, and show significant improvement over global one step mechanisms, such as those by Westbrook and Dryer. In addition, by using an ILDM reaction mechanism in a CID code, a considerable improvement in turn-around time can be achieved.

Norris, A. T.↗