Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic modeling”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-throughput approach for investigating interdiffusion in medium- and high-entropy alloys

Interdiffusion experiments are usually time-consuming and tedious since diffusion couples must be annealed at several temperatures for a long time. The efforts required to study interdiffusion in multicomponent alloys increase dramatically as multiple diffusion couples are required to cover broad composition ranges and determine the diffusivities of individual elements in different chemical environments. To circumvent this challenge, we present a high-throughput approach applicable to single-phase and compositionally complex alloys, which are assumed to approximate ideal solid solutions. Here, a simple diffusion-multiple experiment combined with a physically based kinetic model is proposed to efficiently determine the diffusion coefficients of the constituent elements in quaternary CrFeCoNi alloys. Compared with tracer diffusivities reported in the literature, the results, thus, obtained do not differ by more than a factor of 2 and were obtained from a single interdiffusion experiment. In contrast, the diffusivities simulated with commercial mobility and thermodynamic databases are strongly overestimated by a factor ranging from 1 to 16. Therefore, our approach enables high-throughput determination of diffusivities and can help in the design of alloys for high-temperature applications where diffusion plays a key role.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Kinetics of antigen binding to arrays of antibodies in different sized spots

A fluorescence-based array biosensor has been developed which can measure the binding kinetics of an antigen to an immobilized antibody in real time. A patterned array of antibodies immobilized on the surface of a planar waveguide was used to capture a Cy5-labeled antigen present in a solution that was continuously flowed over the surface. The CCD image of the waveguide was monitored continuously for 25 min. The resulting exponential rise in fluorescence signal was determined by image analysis software and fitted to a reaction-limited kinetics model, giving a kf of 3.6 x 10(5) M(-1) s(-1). Different spot sizes were then patterned on the surface of the waveguide using either a PDMS flow cell or laser exposure, producing width sizes ranging from 80 to 1145 microm. It was demonstrated that under flow conditions, the reduction of spot size did not alter the association rate of the antigen with immobilized antibody; however, as the spot width decreased to < 200 nm, the signal intensity also decreased.

Non-NASA Center↗
DOE OSTI2020

Chemical Modeling of the Capture of Radioiodine Species

Iodine is a common fission product resulting from the transmutation of uranium fuel in nuclear power reactors. The short-lived radioactive isotope, 131 I, with a half-life of around eight days, would be a particular health concern if released into the environment. Capture materials, such as activated carbon, are used routinely in nuclear power plants and medical isotope production facilities to remove this hazard during normal operations to reduce emissions to acceptable levels. Iodine, unlike other common fission products such as noble gases, is challenging because of the wide range of molecular forms and oxidation states (from +VII in IF 7 to -I in iodide) that it displays, which can interconvert both in the gas phase and in the solution phase for example. The work in this project is focused on three computational tasks. In the first, we consider the binding of iodine-containing molecules of varying oxidation states to the major component of activated carbon – pristine graphene-like sheets. The second task focused on molecular binding and possible chemical transformations resulting from the reaction of molecular iodine with chemical functional groups on the edge of the graphene sheets or surface defects in the sheet itself. The final task was associated with radiolytic breakdown of 131 I to 131 Xe in iodine containing molecules and investigated the resulting molecular fragmentation process both in the gas phase and on the graphene-like surface. For each of the tasks, the intent was to learn about the thermodynamics and kinetics of the chemical transformations and use the calculations to provide parameters for input into kinetic models describing these processes for interpretation of experimental data in future work.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Anomalous conductivity due to relativistic Landau quantization

Here we use a recently developed kinetic model derived from the Dirac equation, in order to study electromagnetic wave propagation in superstrong magnetic fields, such as in magnetars, where relativistic Landau quantization is prominent. The leading contribution to the conductivity tensor in such a plasma is calculated. It is found that the electron Hall current has an anomalous contribution, in the quantum relativistic regime, where the effective particle energy has a significant contribution from the diamagnetic and Zeeman energy. As a result, a new quantum resonance frequency appears, and the dispersion relation for the left- and right-hand polarized modes are strongly modified for long and moderate wavelengths. The implications for magnetar physics are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Influence of functional groups on low-temperature combustion chemistry of biofuels

Ongoing progress in synthetic biology, metabolic engineering, and catalysis continues to produce a diverse array of advanced biofuels with complex molecular structure and functional groups. In order to integrate biofuels into existing combustion systems, and to optimize the design of next-generation combustion systems, understanding connections between molecular structure and ignition at low-temperature conditions (< 1000 K) remains a priority that is addressed in part using chemical kinetics modeling. The development of predictive models relies on detailed information, derived from experimental and theoretical studies, on molecular structure and chemical reactivity, both of which influence the balance of chain reactions that occur during combustion – propagation, termination, and branching. In broad context, three main categories of reactions affect ignition behavior: (i) initiation reactions that generate a distribution of organic radicals, ; (ii) competing unimolecular decomposition of and bimolecular reaction of with O2; (iii) decomposition mechanisms of peroxy radical adducts (RO), including isomerization via RO ⇌ OOH. Furthermore, all three categories are influenced by functional groups in different ways, which causes a shift in the balance of chain reactions that unfold over complex temperature- and pressure-dependent mechanisms.

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

On the Rate Constant for NH 2 +HO 2 and Third-Body Collision Efficiencies for NH 2 +H(+M) and NH 2 +NH 2 (+M)

In low-temperature flash photolysis of NH 3 /O 2 /N 2 mixtures, the NH 2 consumption rate and the product distribution is controlled by the reactions NH 2 + HO 2 → products (R1), NH 2 + H (+M) → NH 3 (+M) (R2), and NH 2 + NH 2 (+M) → N 2 H 4 (+M) (R3). In the present work, published flash photolysis experiments by, among others, Cheskis and co-workers, are re-interpreted using recent direct measurements of NH 2 + H (+N 2 ) and NH 2 + NH 2 (+N 2 ) from Altinay and Macdonald. To facilitate analysis of the FP data, relative third-body collision efficiencies compared to N 2 for R2 and R3 were calculated for O 2 and NH 3 as well as for other selected molecules. We report results were in good agreement with the limited experimental data. Based on reported NH 2 decay rates in flash photolysis of NH 3 /O 2 /N 2 , a rate constant for NH 2 + HO 2 → NH 3 + O 2 (R1a) of $k_{1\text{a}}$ = 1.5(±0.5) × 10 14 cm 3 mol –1 s –1 at 295 K was derived. This value is higher than earlier determinations based on the FP results but in good agreement with recent theoretical work. Kinetic modeling of reported N 2 O yields indicates that NH 2 + HO 2 → H 2 NO + O (R1c) is competing with R1a, but perturbation experiments with addition of CH4 indicate that it is not a dominating channel. Measured HNO profiles indicate that this component is formed directly by NH 2 + HO 2 → HNO + H 2 O (R1b), but theoretical work indicates that R1b is only a minor channel. Based on this analysis, we estimate $k_{1\text{c}}$ = 2.5 × 10 13 cm 3 mol –1 s –1 and $k_{1\text{b}}$ = 2.5 × 10 12 cm 3 mol –1 s –1 at 295 K, with significant uncertainty margins.

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

The effect on soot and its gas precursors of doping ethylene with 2,2,4,6,6-pentamethyl-heptane in the nitrogen-fuel stream of a laminar non-premixed Planar Mixing Layer Flame (PMLF)

Synthetic Aviation Turbine Fuels (SATFs) are promising for reducing soot emissions from the aviation sector and diversifying Jet Fuel (JF) sources. Accurately predicting the combustion and emissions behavior of SATFs (and other JFs) necessitates robust experimental databases to elucidate the chemistry of long-chain iso-paraffins, which can compose up to two-thirds of SATF blends and whose behavior is considered to be well-represented by that of iso-dodecane isomers. Here, this study characterizes two laminar non-premixed Planar Mixing Layer Flames (PMLFs) with mild soot loads fueled by nitrogen-diluted ethylene, pure and doped with 2,2,4,6,6-pentamethyl-heptane, respectively. The two PMLFs have the same stoichiometric mixture fraction and total hydrocarbon mole fraction in the fuel stream (X F,F =X C2H4,F +X C12H26,F = 0.260), resulting in nearly the same maximum temperature (T max ≈1800 K) and simple identification of the effects of doping. Importantly, any horizontal PMLF cross-section has a self-similar structure that can be modeled as an equivalent One-Dimensional Counterflow Flame (1D-CF) with vanishingly small strain rate (a). The cross-section at a Height Above the Burner (HAB) of 50 mm is characterized in terms of C 0 -C 18 gas species using capillary sampling followed by GC-MS analyses. Laser-Induced Emission Spectroscopy (LIES) quantifies the soot volume fraction (ƒ v ) profiles at HAB=25 and 50 mm where Elastic Laser Light Scattering (E-LLS) is performed to determine the a of the equivalent 1D-CFs and the profile of the E-LLS equivalent diameter ( d 6,3 ) of soot. The substitution of 1500 ppm of ethylene with 2,2,4,6,6-pentamethyl-heptane causes an increase of ≈1.5 in the concentrations of several polycyclic aromatic hydrocarbons and fv. Concurrently, the measured d 6,3 doubles in the oxidizer stream, yet remains the same in the fuel stream, at HAB=50 mm. Instead, at HAB= 25 mm, the iso-dodecane doping does not affect the d 6,3 profile in either stream. The experimental results partially validate the chemical reactions and soot formation kinetic model developed at Lawrence Livermore National Laboratory and provide directions to further improve its predictions.

Iso-dodecane (2,2,4,6,6-Pentamethyl-Heptane)↗
DOE OSTI2025

Water‐mediated Decomposition Pathways of ε‐Hexanitrohexaazaisowurtzitane

The dominant initial decomposition mechanism of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (HNIW) is understood from the literature to be the N─NO 2 bond scission unimolecular mechanism. Here we have found a water-mediated HONO release initial decomposition mechanism energetically competitive with the N─NO 2 bond scission mechanism modeled using ab initio nudged elastic band simulations. The activation energy of the water-mediated HONO mechanism was calculated to be 138.9 kJ/mol with the Perdew, Burke, and Ernzerhof (PBE) functional (177.8 with PBE0), while the double NO 2 release mechanism is 171.9 kJ/mol (203.3 with PBE0). Branching secondary, tertiary, and quaternary decomposition steps were also discovered including oxidation of HONO or H 2 O causing ring-opening that leads to C─N bond breaking and N 2 O or NO release. The reaction rate of HONO oxidation is much faster than HONO release, making the HONO concentration low. The Helmholtz free energy barrier of the double NO 2 release mechanism is lower than the water-mediated HONO release above 438 K (165°C) due to the vibrational contribution to the free energy. This helps explain why the literature reports more NO 2 release at higher temperatures. The present study reveals a physical mechanism for how water can catalyze the decomposition of HNIW at low enough temperatures, and reveals secondary, tertiary, and quaternary reaction mechanisms leading to NO x release. In conclusion, as water is realistically always present, the reported barriers are important for building kinetic models needed to ensure the safety and reliability of HNIW.

Steele, Brad A. [Lawrence Livermore National Labor↗
DOE OSTI2021

Effects of non-thermal termolecular chemistry on detonation development in hydrogen (H2) / methane (CH4) - air mixtures

In combustion simulations, it is usually assumed that highly-energized collision complexes (either radical or stable species) formed in exothermic reactions are always under thermal equilibrium as they evolve through competing networks of reactions. However, in practical flames, non-trivial amounts of reactive radicals such as H, O and OH are present apart from O2. As a result, collisions of these reactive species with the energized collision complexes have been shown in recent studies to induce non-thermal reactivity. These studies have also demonstrated that such non-thermal reactions can be suitably represented in macroscopic kinetics models as chemically termolecular reactions. The present work was focused on identifying and quantifying the effects of including such reactions on the evolution of an initial deflagration front to a developing detonation in H2/CH4-air mixtures under boosted internal combustion (IC) engine conditions. Specifically, fully resolved simulations, with and without non-thermal reactivity, were performed for a constant volume reactor containing stoichiometric H2/CH4-air mixture. It was found that inclusion of termolecular chemistry resulted in a delayed onset of end-gas auto-ignition. Concurrently, the developing detonation intensity was observed to be significantly higher. Chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions responsible for the observed variation in the results.

Detonation↗
NASA NTRS1995

The mosaic structure of plasma bulk flows in the Earth's magnetotail

Moments of plasma distributions observed in the magnetotail vary with different time scales. In this paper we attempt to explain the observed variability on intermediate timescales of approximately 10-20 min that result from the simultaneous energization and spatial structuring of solar wind plasma in the distant magnetotail. These processes stimulate the formation of a system of spatially disjointed. highly accelerated filaments (beamlets) in the tail. We use the results from large-scale kinetic modeling of magnetotail formation from a plasma mantle source to calculate moments of ion distribution functions throughout the tail. Statistical restrictions related to the limited number of particles in our system naturally reduce the spatial resolution of our results, but we show that our model is valid on intermediate spatial scales Delta(x) x Delta(z) equal to approximately 1 R(sub E) x 1000 km. For these spatial scales the resulting pattern, which resembles a mosaic, appears to be quite variable. The complexity of the pattern is related to the spatial interference between beamlets accelerated at various locations within the distant tail which mirror in the strong near-Earth magnetic field. Global motion of the magnetotail results in the displacement of spacecraft with respect to this mosaic pattern and can produce variations in all of the moments (especially the x-component of the bulk velocity) on intermediate timescales. The results obtained enable us to view the magnetotail plasma as consisting of two different populations: a tailward-Earthward system of highly accelerated beamlets interfering with each other, and an energized quasithermal population which gradually builds as the Earth is approached. In the near-Earth tail, these populations merge into a hot quasi-isotropic ion population typical of the near-Earth plasma sheet. The transformation of plasma sheet boundary layer (PSBL) beam energy into central plasma sheet (CPS) quasi-thermal energy occurs in the absence of collisions or noise. This paper also clarifies the relationship between the global scale where an MHD description might be appropriate and the lower intermediate scales where MHD fails and large-scale kinetic theory should be used.

Ashour-Abdalla, M.↗
DOE OSTI2021

Efficient and reversible electron bifurcation with either normal or inverted potentials at the bifurcating cofactor

A longstanding mystery surrounding electron bifurcation is the significance of inverted (or “crossed”) reduction potentials of the two-electron bifurcating cofactor. Using a many-electron open-system kinetic model, we show that reversible and efficient electron bifurcation is possible without inverted reduction potentials at the bifurcating site if the absolute value of the difference between first and second reduction potentials of the bifurcating species is sufficiently large (on the scale of the redox-potential span of the high- and low-potential branches). Surprisingly, the case with strong, normally ordered potentials at the bifurcating cofactor can produce electron bifurcation that is just as effective as the case with strongly inverted potentials. Lastly, this finding amplifies the puzzle surrounding the recruitment of inverted potentials in the few well-characterized bifurcating systems of nature and suggests that electron bifurcating cofactors without strongly inverted potentials may yet be discovered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2022

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗
DOE OSTI2025

Intrinsic activation energies for ring contraction of allylic cations in zeolites

Cyclic carbocations are important intermediates in zeolite-catalyzed chemistries such as methanol-to-hydrocarbons conversion, naphthenes ring opening, or coke formation. While information about their thermodynamic stability exists, little is known about the kinetics of formation and transformation of cyclic carbocations in zeolites. To fill this knowledge gap, ring contraction of the 1,3,5,5-tetramethylcyclohexenyl cation (C 10 H 17 + ), a representative of 6-membered ring allylic cations, was investigated by in situ UV–vis and IR spectroscopy. Protonic forms of zeolites served as catalysts, at temperatures from 80 °C to 135 °C. Significant oligomerization and hydride transfer in BEA and FAU hampered kinetics analysis, whereas ring contraction dominated in the channels of MOR. The reactant cation, characterized by an electronic absorption at 314 nm and an allylic stretch at 1549 cm −1 , contracted to both a 1,3-alkyl-substituted cyclopentenyl cation (287 nm and 1506 cm −1 ) and a 1,2,3-alkyl-substituted cyclopentenyl cation (297 nm and 1489 cm −1 ). Collection of time-resolved IR spectra and fitting of the intensities with various kinetic models revealed a third transformation, which is expected from thermodynamics: the 1,3-alkyl-substituted cyclopentenyl cation isomerizes to the 1,2,3-alkyl-substituted cyclopentenyl cation. Series of IR spectra recorded at different temperatures delivered intrinsic activation enthalpies (entropies) in MOR of 67 ± 2 kJ mol −1 (−130 ± 6 J mol −1 K −1 ) and 90 ± 3 kJ mol −1 (−70 ± 8 J mol −1 K −1 ) for the contraction to 1,3- and 1,2,3-substituted species, and of 88 ± 5 kJ mol −1 (−84 ± 12 J mol −1 K −1 ) for the isomerization. The findings characterize one path – via contraction of larger rings – to different cyclopentenyl species in zeolites; and the associated, moderate activation energies suggest such transformations contribute to many complex hydrocarbon reaction networks.

Acid catalysis↗
DOE OSTI2026

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗
DOE OSTI2023

A projection-based analytical Jacobian framework for chemical kinetics applications

A major challenge in simulating complex combustion systems with detailed chemical kinetic models is the cost of integrating the chemical source terms, often done using stiff ODE solvers that require frequent Jacobian evaluations. Using analytically derived Jacobian matrices instead of divided-difference-based numerical Jacobian approximations can significantly reduce the associated computational cost. However, ambiguities arise in the formulation of analytical Jacobians because the chemical state of the system, or state vector, can be expressed in multiple ways, involving variables that are typically not independent from one another. Here, in this work, the consequences of those ambiguities on practical calculations are characterized in detail, and a generalized, projection-based framework is proposed as a mitigation strategy. Performances are assessed in a series of test cases involving a variety of configurations and numerical solution approaches. Results show that with proper treatment, commonly used analytical Jacobian formulations can be considered as equivalent for practical purposes, thereby alleviating concerns that the state vector chosen to express the governing equations and corresponding analytical Jacobian may significantly impact the accuracy of the simulations.

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

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

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

Influence of gasoline fuel formulation on lean autoignition in a mixed-mode-combustion (deflagration/autoignition) engine

We report stoichiometric spark-ignition engines suffer efficiency penalties due to throttling losses at low loads, a low specific-heat ratio of the stoichiometric working fluid, and limits on compression ratio due to end-gas autoignition leading to undesirable knocking. Mixed-Mode Combustion (MMC) mitigates these shortcomings by using a lean working fluid where a spark-initiated pilot-stabilized deflagrative flame front is followed by controlled end-gas autoignition. This MMC study investigates the effects of initial conditions (intake air temperature, intake pressure, equivalence ratio, and intake oxygen fraction) on autoignition tendency of four gasoline-range fuels with varying properties and composition. The use of fuels with varying octane sensitivity (S) allowed exploring the importance of low-temperature heat release in triggering autoignition. Fuels with high S were less reactive for conditions that promote low-temperature chemistry (operation at high intake air pressure or without N2 dilution). Conversely, an Alkylate fuel with low S showed a greater autoignition resistance at operating conditions that were unfavorable for low-temperature chemistry. Next, the effect of residual gas composition on autoignition tendency of fuels was examined with a chemical-kinetics model. Among the various molecules in the residual gas, nitric oxide (NO) enhanced the low-temperature chemistry and increased the autoignition tendency most significantly. The fuels’ autoignition response to increasing NO amount corroborates the experimental observations. Next, the sequential autoignition of the end-gas was assessed to be less impacted by thermal stratification because of lean mixtures showing relatively less low-temperature chemistry, when compared to stoichiometric mixtures. Next, the effect of changing equivalence ratio on the autoignition was found to be similar for all fuels, regardless of their S. With changing intake air temperature, the response of fuels’ autoignition tendency depended on the dilution level used. At high dilution (i.e. low intake [O2]), fuels’ reactivity increased with increasing intake air temperature. In contrast, for operation without dilution, the autoignition tendency of the low-S Alkylate fuel decreased with increasing intake air temperature, while that of high-S High Cycloalkane fuel still increased with increasing intake air temperature. In conclusion, conventional octane metrics (RON and MON) have utility in assessing the autoignition tendency under lean MMC operation. Moreover, the fuel requirements for MMC align with that of stoichiometric operation: i.e., high RON and high S fuels are desirable for stable non-knocking operation.

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

A workflow for automatic generation and efficient refinement of individual pressure-dependent networks

Manual chemical kinetic model construction often requires modelists to at least implicitly guess all possible decay pathways and their relative fluxes for all important species. We present a workflow and associated tool for automatic generation and re finement of pressure-dependent networks that should enable modelists to efficiently and comprehensively identify decay pathways and estimate respective parameters for chemical species. This tool combines the capabilities of the Reaction Mechanism Generator (RMG) software to generate possible reaction paths, determine thermochemistry, approximate rate coefficients and estimate frequencies with the capabilities of the Arkane software to use quantum chemical parameters to compute thermochemistry and pressure dependent rates. A flux-based algorithm is used to decide which isomers to add to the network. Isomers added to the network are reacted to form other channels. When enough of the flux is accounted for by the isomers and bimolecular product channels, the generation process terminates to yield a comprehensive network. Network sensitivity analysis is applied to the network to identify the most important wells and barriers that can be refined using quantum chemistry calculations. Iterative refinement can then be used to achieve accuracy. The comprehensive network can be easily reduced using energy and flux-based algorithms to the most important channels and isomers in the network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗