Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic model data”

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 217 records · Page 12

Stiff neural ordinary differential equations

Neural Ordinary Differential Equations (ODEs) are a promising approach to learn dynamical models from time-series data in science and engineering applications. This work aims at learning neural ODEs for stiff systems, which are usually raised from chemical kinetic modeling in chemical and biological systems. We first show the challenges of learning neural ODEs in the classical stiff ODE systems of Robertson’s problem and propose techniques to mitigate the challenges associated with scale separations in stiff systems. We then present successful demonstrations in stiff systems of Robertson’s problem and an air pollution problem. The demonstrations show that the usage of deep networks with rectified activations, proper scaling of the network outputs as well as loss functions, and stabilized gradient calculations are the key techniques enabling the learning of stiff neural ODEs. The success of learning stiff neural ODEs opens up possibilities of using neural ODEs in applications with widely varying time-scales, such as chemical dynamics in energy conversion, environmental engineering, and life sciences.

97 MATHEMATICS AND COMPUTING↗

Surrogate Modeling of Landau Damping with Deep Operator Networks

Kinetic simulations excel at capturing microscale plasma physics phenomena with high accuracy, but their computational demands make them impractical for modeling large-scale space and astrophysical systems. In this context, we build a surrogate model, using Deep Operator Networks (DeepONets), based upon the Vlasov–Poisson simulation data to model the dynamical evolution of plasmas, focusing on the Landau damping process—a fundamental kinetic phenomenon in space and astrophysical plasmas. The trained DeepONets are able to capture the evolution of electric field energy in both linear and nonlinear regimes under various conditions. Extensive validation highlights DeepONets’ robust performance in reproducing complex plasma behaviors with high accuracy, paving the way for large-scale modeling of space and astrophysical plasmas.

plasma astrophysics↗

Fuel molecular structure effect on autoignition of highly branched iso-alkanes at low-to-intermediate temperatures: Iso-octane versus iso-dodecane

Highly branched iso-alkanes are an important class of hydrocarbons found in conventional petroleum-derived and alternative renewable fuels used for combustion applications. Recognizing that chemical kinetics for most of these iso-alkanes, especially at low-to-intermediate temperatures, has not been well studied, an experimental and modeling investigation of two selected iso-alkanes, iso-octane (2,2,4-trimethylpentane, iC8) and iso-dodecane (2,2,4,6,6-pentamethylheptane, iC12), is conducted to understand the fuel molecular structure effect on their autoignition characteristics. Using a rapid compression machine (RCM), the ignition responses of iC8 and iC12 at varying pressures, temperatures, and equivalence ratios are characterized and compared. The newly-acquired experimental ignition delay times have been compared with the literature RCM and shock tube data, demonstrating the complementary nature of the current dataset. Further comparison of the experimental pressure traces and ignition delay times illustrates the reactivity crossover between iC8 and iC12. Namely, there exists a temperature window in the negative temperature coefficient regime within which iC12 is less reactive than iC8, but iC12 becomes more reactive outside this temperature window. Furthermore, a chemical kinetic model of iso-alkanes including both iC8 and iC12 is developed. Simulated results using this model are then compared to the experimental data obtained in this study and available in the literature, showing its ability to predict the experimental trends. Chemical kinetic analyses have also been conducted to identify the important reaction pathways controlling autoignition at varying conditions, and to elucidate the underlying mechanism leading to different reactivity trends between iC8 and iC12.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring CO 2 reduction and crossover in membrane electrode assemblies

Electrochemical CO 2 reduction (CO 2 R) using renewable electricity is a key pathway toward synthesizing fuels and chemicals. In this study, multi-physics modeling is used to interpret experimental data obtained for CO 2 R to CO using Ag catalysts in a membrane electrode assembly. The one-dimensional model is validated using measured CO 2 crossover and product formation rates. The kinetics of CO formation are described by Marcus–Hush–Chidsey kinetics, which enables accurate prediction of the experimental data by accounting for the reorganization of the solvent during CO 2 R. Further, the results show how the performance is dictated by competing phenomena including ion formation and transport, CO 2 solubility, and water management. The model shows that increasing the ion-exchange capacity of the membrane and surface area of the catalyst increases CO formation rates by >100 mA cm –2 without negatively impacting CO 2 utilization. Here we provide insights into how to manage the trade-off between productivity and CO 2 utilization in CO 2 electrolyzers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical Kinetics of the Autoxidation of Poly(ethylenimine) in CO 2 Sorbents

The oxidative degradation rates of a CO 2 sorbent composed of a mesoporous alumina impregnated with poly- (ethylenimine) (PEI) are measured under systematically varied conditions and a reaction rate law is created. Good agreement is shown between the rate of oxidation obtained via in situ calorimetric heat measurement during oxidative degradation reactions and the loss of CO 2 capture performance presented as amine efficiency (mol CO 2 /mol amine). PEI mass loss and elemental composition are tracked over the course of the reaction and used in conjunction with the oxidation rate measurements to shed insight into the oxidation reaction(s). These data, in combination with measurements of the heat of reaction, suggest a common reaction set across the range of temperatures, oxygen concentrations, and sorbent compositions tested. The data are consistent with the basic autoxidation scheme (BAS), the accepted mechanism of autoxidation of aliphatic polymers. We propose a lumped kinetic model to describe the oxidation reaction set and estimate an activation energy of 105 kJ/mol and an oxygen reaction order of 0.5–0.7 from the data accordingly. Furthermore, these parameters can be incorporated into process cycle models to estimate the material lifetime, a critical uncertainty in the deployment of DAC technologies.

36 MATERIALS SCIENCE↗

Characterization of shale using Helium and Argon at high pressures

In order to estimate the shale gas in place and the eventual recovery during shale gas operations, one of the key requirements is to accurately characterize the shale’s petrophysical and transport properties such as porosity, permeability, diffusivity, and storage capacity. Despite the many efforts reported in the technical literature aiming to provide an improved understanding of the complex pore structures and the associated fluid flow in gas shales, complete characterization of organic-rich shale samples still poses a big challenge. Here, we have characterized mass transfer and sorption in shale at different length scales using Helium (He) and Argon (Ar) as probe gases. Thermogravimetric analysis (TGA) with a shale cube of ~1 cm 3 in size and gas expansion experiments with a full-diameter core (3.5” in diameter) were used to measure sorption kinetics/isotherms and mass transfer, respectively. Both samples are from the same depth/location in the Marcellus shale formation. The TGA steady-state technique was initially used to generate excess sorption isotherms for Ar, while dynamic TGA experiments were used to study its sorption kinetics. The TGA experiments demonstrate that Ar, which has a similar sorption potential as Methane, but is generally assumed to be inert, adsorbs onto the surfaces of the mesoporous and microporous regions of the shale samples according to a Langmuir-type behavior. Helium expansion experiments, on the full-diameter core, were used to measure the overall porosity, on the basis that He is a non-sorbing and inert gas as compared to Ar. The He expansion experiments, furthermore, allow us to delineate the mass transfer of gas across the inherent hierarchy of pore sizes, including macropores (macro- and microcracks), mesopores and micropores. Similar expansion experiments were also performed with Ar to study the combined impact of mass transfer and sorption. A triple-porosity model (TPM) was utilized to interpret the He expansion experiments with the shale core and to extract (estimate) relevant transport parameters. We report and compare here the diffusivities and permeabilities of the whole core for both He and Ar, as calculated from the modeling and fitting of the experimental data. On the premise that the shale cube is representative of the matrix region of the core, the Ar sorption kinetics from the cube experiments were subsequently combined with the transport parameters extracted from the He experiments to predict the behavior of the Ar expansion test with the full-diameter core. An excellent agreement is observed between the model predictions and the experimental data. The experimental observations and their interpretation indicate that one must be cautious when using Ar to estimate the true porosity and permeability of shales. In addition, we demonstrate that He and Ar probe gases, when used in tandem, can be employed effectively as a tool to characterize shales in terms of mass transfer and sorption dynamics across scales.

04 OIL SHALES AND TAR SANDS↗

Autoignition study of iso -cetane/tetralin blends at low temperature

iso-Cetane and tetralin are the two fuel components commonly considered in literature diesel surrogates, and are hydrocarbon classes representative of iso-alkanes and naphthoaromatics, respectively. Since both surrogate components are involved in the key reactions controlling the autoignition reactivity of diesel surrogates, autoignition studies of iso-cetane, tetralin, and their blending behavior at low temperatures can be helpful in developing/refining chemical kinetic models of diesel surrogates, in order to better predict diesel ignition response under low temperature combustion conditions. In this investigation, autoignition experiments of iso-cetane, tetralin, and their binary mixtures at varying blending ratios have been conducted in a rapid compression machine (RCM) in the temperature/pressure range of 630–930 K and 10–20 bar. For the neat fuel components, the newly-acquired RCM data have been compared with the literature experimental results at some overlapping conditions. In addition, a chemical kinetic model describing the oxidation of iso-cetane and tetralin compiled/updated from the diesel surrogate model developed by Lawrence Livermore National Laboratory is used to simulate the experimental datasets. It is shown that this model exhibits good agreement with the experimental results of tetralin and also predicts the total ignition delay times of iso-cetane at lower temperatures reasonably well. However, the current model fails to fully capture the increase of total ignition delay time with increasing amount of tetralin in the binary blends observed in the present RCM experiments. Therefore, chemical kinetic analyses of the current model are conducted to help identify possible reasons leading to the discrepancies between simulations and experiments. Furthermore, the chemical kinetic interactions between iso-cetane and tetralin are explored by including additional cross-reactions among iso-cetane, tetralin, and their fuel fragments in the current model to assess their effects on model predictions. As a result, based on the chemical kinetic analysis results, future direction for model refinements is discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Vertical effective stress and temperature as controls of quartz cementation in sandstones: Evidence from North Sea Fulmar and Gulf of Mexico Wilcox sandstones

Here we present quantitative petrographic data, high spatial resolution oxygen isotope analyses of quartz cement, basin modelling and a kinetic model for quartz precipitation for two Paleocene-Eocene Wilcox Group sandstones from Texas and two Jurassic Fulmar Formation sandstones from the Central North Sea. At each location, one sandstone has been buried to ca. 145 °C and one to ca. 185 °C. A key difference between the Wilcox and Fulmar burial histories is that the Wilcox sandstones are currently at higher vertical effective stresses and, from basin modelling studies, have been subjected to generally higher vertical effective stresses through their burial history. The amounts of quartz cement in the Wilcox sandstones are between 12 and 18%, and between 2 and 6% in the Fulmar sandstones. High-spatial-resolution oxygen isotope data obtained from the quartz cements suggest temperature ranges for quartz precipitation from 60 to 80 °C to values approaching maximum burial temperature. Factors such as grain coatings or the timing of petroleum emplacement cannot explain the differences in the amounts of quartz cement. Petrographic data show that most of the silica for quartz cement can be derived from intergranular pressure dissolution. Although the sample set is small, we interpret the results to suggest that the differences in quartz cementation in Fulmar and Wilcox sandstones can be explained better by differences in their vertical effective stress history than their temperature history; in this case, the supply of silica rather than the precipitation of quartz becomes an important control on the rate and extent of cementation.

58 GEOSCIENCES↗

The role of collisional energy transfer in the thermal and prompt dissociation of 1-methyl allyl

Collisional energy transfer plays a pivotal role in a-priori calculations of the pressure-dependent rate constants obtained from a master equation. However, accurate determinations of collisional energy transfer parameters are rare and so most kinetic studies rely on best-fits of these parameters to experimental measurements (if available) or estimates guided from literature studies. In this work, we have quantified the effect of the uncertainty in energy transfer parameters on the thermal and prompt dissociation kinetics of a resonance stabilized radical, 1-methyl allyl (1MA), of relevance to the combustion of 1- and 2-butene isomers. Simulations using literature kinetics models were performed to assess the impact of these uncertainties on flame propagation and speciation data in laminar flames of 1- and 2-butene. Analyses of the uncertainty propagated by the energy transfer parameters for 1MA dissociation to the flame simulations, in particular the laminar flame speed, indicate an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Ab-initio based theoretical calculations were also performed to obtain pressure-dependent kinetics for the reaction of 1MA with H-atoms. Lastly, theoretically calculated energy transfer parameters were used to best characterize the kinetics and branching between the chain propagating 1MA dissociation and the chain-terminating reaction 1MA + H.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A study of silver acetate under extreme conditions

With the aim of exploring chemical systems that may undergo metallization when irradiated with hard X-rays, we selected silver acetate (AgC 2 H 3 O 2 ) for the subject of this study. X-ray-induced decomposition of silver acetate under ambient and high-pressure conditions was observed in a diamond anvil cell (DAC), leading to the formation of metallic nanograins of silver at ambient pressure and 1.65 GPa. At 4 GPa, no decomposition was observed. The Avrami kinetics equation also provides information about novel structural formation at ambient pressure and 1.65 GPa. By modeling of the XRD data, it was found that the size of the silver nanocrystallites formed at 1.65 GPa pressure steadily increased to ∼5 nm after 600 min of X-ray irradiation as determined by applying the Scherrer equation to the diffraction peak widths. Time-resolved X-ray diffraction (XRD) revealed pressure-dependent kinetics, demonstrating that coupling pressure with irradiation enables controlled photochemical pathways in this model system. Concurrent with previous studies, the application of high pressure (HP) can be considered as a means of controlling X-ray synthetic photochemistry.

36 MATERIALS SCIENCE↗

Probing High-Temperature Amine Chemistry: Is the Reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 Important?

The reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 (R1) has been identified as a key step to explain experimental results for pyrolysis and oxidation of ammonia. However, no direct experimental or theoretical evidence for the reaction has been reported. In the present work, the reaction was studied by ab initio theory and by reinterpretation of experimental data. We could not locate a transition state for R1 occurring as a direct process, but alternative mechanisms yielded an upper bound to k 1 of 1.5 × 10 13 exp(–58.9 kcal mol –1 /RT) cm 3 mol –1 s –1 over 1000–2500 K, several orders of magnitude lower than values applied in modeling. Consistent with the theoretical work, re-evaluation of NH 3 pyrolysis data supported a very low value of k 1 . However, this finding opens up a novel unresolved issue. Current kinetic models cannot capture the NH 3 oxidation behavior in a number of laminar flow reactor and jet-stirred reactor experiments without adopting an improbably high value for k 1 . Furthermore, important oxidation steps might be underestimated or missing from mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lessons Learned in Employing Data Analytics to Predict Oxidation Kinetics and Spallation Behavior of High-Temperature NiCr-Based Alloys

Machine learning (ML) can offer many advantages in predicting material properties over traditional materials development methods based solely on limited experimental investigations or physical-based simulations with the capability to reduce development cost, risk, and time. However, so far, limited efforts have been made to predict alloy oxidation kinetics and spallation behavior via ML due to the lack of consistently measured and sufficient experimental data and the inherent complexity in oxidation behavior of multicomponent high-temperature alloys. A previous study reported the ability of ML to predict oxidation kinetics of NiCr-based alloys as a function of alloy composition and operating conditions. Here, the performance of a ML model in predicting rate constants and spallation probability was evaluated in light of the roles of the data distribution of the experimental dataset (data analytics), the alloy composition, the exposure environment and the chosen oxidation approach to extracting kinetic values from the measured mass changes (but using either a simple parabolic law or a statistical cyclic oxidation model). Potential strategies to improve the predictions and enhance the extrapolative capability of the previously trained model will be discussed.

36 MATERIALS SCIENCE↗

A Coupled Small and Wide-Angle X-Ray Scattering Study of Phase Transformation Mechanisms in U-6Wt Pct Nb

The elevated temperature phase transformation kinetics of U-6wt%Nb in both wrought and cast conditions were investigated using in-situ high energy synchrotron X-ray scattering. The evolution of the phase content, γ-phase composition, and microstructure morphology was determined through the analysis of simultaneously collected wide-angle diffraction and small-angle scattering (SAXS) data. The results of the wide-angle scattering analysis, especially the phase fractions and average composition of the γ-phase, provide model parameters to enable a quantitative analysis using the SAXS data. Greater clarity regarding the kinetics of different precipitation mechanisms, namely the continuous precipitation (CP) of α-U within the γ-phase matrix and the discontinuous (DP) grain boundary precipitation of a lamellar microstructure, consisting of α-U and Nb-rich γ-phase are obtained from the coupled measurement modes. These results indicate that the CP mechanism is responsible for the initial precipitation of α-U at all temperatures considered. Near the nose in the time-temperature-transformation (TTT) diagram at ~798 K, both mechanisms are active during the early stages of phase transformation. However, the DP mechanism has an increasingly long incubation period with decreasing temperature below the TTT nose. In conclusion, the average interlamellar spacing obtained from the SAXS analysis increases with time, providing further evidence of the so-called divergent DP behavior previously observed in U-6wt%Nb alloys.

36 MATERIALS SCIENCE↗

AWSD Reactive Burn Model for the HMX‐Based High Explosive LX‐04

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive burn model is applied to describe shock initiation and detonation properties of the HMX-based high explosive LX-04. The parameters in the model are calibrated to data from multiple sources. The thermodynamic equations of state used in the model are calibrated to a combination of thermochemical calculations for HMX and LX-04 as well as experimentally-measured cylinder expansion results for LX-04. The kinetic parameters are calibrated to velocity data from gas gun experiments performed using EDC-32—a high explosive with the same chemical composition as LX-04 but different structural properties, and scaled rate stick data for PBX 9012. The AWSD model is shown to accurately describe the shock initiation and propagation of LX-04. Very good agreement is observed between the available experimental data and the AWSD model output. The presented results constitute an accurate LX-04 reactive burn model for use in engineering-scale models and simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A Semi-Detailed Pyrolytic Gas-Phase Kinetic Model for the Volatiles of Polyethylene Thermal Degradation

This work presents a semi-detailed kinetic model to address the pyrolytic gas-phase reactivity of volatiles formed during thermal degradation of polyethylene (PE). The model builds on a validated multi-step condensed-phase model and employs validated lumping approaches. Short-chain compounds are modelled with high detail, while long-chain ones are described by surrogate species representative of diesel-cuts (NC16H32) and waxes (NC30H60). The reactivity of short chains is described through the comprehensive CRECK kinetic model, updated to align C5-C7 olefins based on recent literature experimental data. Due to the lack of experimental data for longer olefins, their reactivity is modeled by analogy to the shorter ones, ensuring an asymptotic behavior with increasing carbon numbers. The semi-detailed model is validated through experimental data on PE pyrolysis, assuming an instantaneous mixing of the inert inlet flow with released volatiles, followed by a segregated plug-flow behavior. Validation across different reactor setups confirms the model’s capability to predict detailed product distributions. Despite minor discrepancies, the proposed model effectively captures experimental trends. Further work will address modelling the reactivity in oxygen-containing environments.

kinetics↗

A Semi-Detailed Pyrolytic Gas-Phase Kinetic Model for the Volatiles of Polyethylene Thermal Degradation

This work presents a semi-detailed kinetic model to address the pyrolytic gas-phase reactivity of volatiles formed during thermal degradation of polyethylene (PE). The model builds on a validated multi-step condensed-phase model and employs validated lumping approaches. Short-chain compounds are modelled with high detail, while long-chain ones are described by surrogate species representative of diesel-cuts (NC16H32) and waxes (NC30H60). The reactivity of short chains is described through the comprehensive CRECK kinetic model, updated to align C5-C7 olefins based on recent literature experimental data. Due to the lack of experimental data for longer olefins, their reactivity is modeled by analogy to the shorter ones, ensuring an asymptotic behavior with increasing carbon numbers. The semi-detailed model is validated through experimental data on PE pyrolysis, assuming an instantaneous mixing of the inert inlet flow with released volatiles, followed by a segregated plug-flow behavior. Validation across different reactor setups confirms the model’s capability to predict detailed product distributions. Despite minor discrepancies, the proposed model effectively captures experimental trends. Further work will address modelling the reactivity in oxygen-containing environments.

kinetics↗

Material analysis of metal hydrides for bulk hydrogen storage

The United States Department of Energy (DOE) Hydrogen Storage Engineering Center of Excellence (HSECoE), 2009–2014, was established to evaluate the viability of material-based hydrogen storage for light duty vehicles. The HSECoE concluded that there were no metal hydrides that simultaneously met the DOE targets for gravimetric and volumetric capacity, and rate of onboard refueling. However, it was noted by the DOE that potential applications existed for stationary hydrogen storage systems. An application of particular interest is using a metal-hydride to store a hydrogen supply for fuel cell backup power to a 5 MWe data center. A generalized numerical model that coupled mass and energy conservation, chemical kinetics and thermodynamics was developed to simulate the discharge process for a given metal-hydride. The model was applied to an interstitial and a complex metal-hydride. In conclusion, it was found that the interstitial metal-hydride, which is discharged by a pressure swing is preferrable in a practical sense to the complex metal-hydride, which requires a temperature swing discharge process.

08 HYDROGEN↗