Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced chemical kinetics”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Reduced chemical kinetics for propane combustion

It is pointed out that a detailed chemical kinetics mechanism for the combustion of propane consists of 40 chemical species and 118 elementary chemical reactions. An attempt is made to reduce the number of chemical species and elementary chemical reactions so that the computer run times and storage requirements may be greatly reduced in three-dimensional gas turbine combustion flow calculations, while maintaining accurate predictions of the propane combustion and exhaust emissions. By way of a sensitivity analysis, the species of interest and chemical reactions are classified in descending order of importance. Nineteen species are chosen, and their pressure, temperature, and concentration profiles are presented for the reduced mechanisms, which are then compared with those from the full 118 reactions. It is found that 45 reactions involving 27 species have to be kept for comparable agreement. A comparison of the results obtained from the 45 reactions to that of the full 118 shows that the pressure and temperature profiles and concentrations of C3H8, O2, N2, H2O, CO, and CO2 are within 10 percent of maximum change.

Ying, Shuh-Jing↗
DOE OSTI2026

Effect of NO on DME-Methanol HCCI Combustion Using a Reduced Chemical Kinetics Mechanism

Methanol is an attractive fuel for the maritime sector due to its wide availability. Its direct use as a fuel, however, is accompanied by challenges such as high latent heat of vaporization and low cetane number. A potential solution to overcome the ignition properties of methanol could be through on-board generation of dimethyl ether (DME) via catalytic dehydration of methanol. The resulting mixture from dehydration can be mixed in with the intake air to generate a homogenous charge compression ignition (HCCI) preburn for subsequent direct injection (DI) mixing controlled compression ignition (MCCI) of neat methanol. Within that context, complementary experimental work found that the influence of combustion residuals on the heat release rate (HRR) was significant, specifically for residual NO. This finding motivated the present computational and kinetic evaluation of the effects of NO on the low (LTHR) and high (HTHR) temperature heat release rates. The strong influence of small quantities of NO on the combustion process of a DME/methanol/H2O mixture (low catalyst or reactor efficiency) necessitated a kinetics-based investigation into this phenomenon. A mechanism sourced from the existing literature with NO had 172 species and 1375 reactions, making it computationally expensive for use. Hence, a mechanism reduction effort was implemented, and a rate constant (k) tuning effort based on sensitivity analysis was needed to validate experimental results using a zero-dimensional engine model in Cantera. The reduced mechanism was able to successfully capture the negligible influence of NO addition on DME HCCI combustion, whereas an advancement in LTHR and HTHR for a DME/methanol/H2O mixture was kinetically confirmed. Reaction pathway analysis showed that addition of NO chemically counteracted the OH sink created by alcohols like methanol, increasing the effectiveness of DME ignition.

Tyrewala, Daanish [ORNL] (ORCID:0000000208599324)↗
NASA NTRS2006

Development of Comprehensive Reduced Kinetic Models for Supersonic Reacting Shear Layer Simulations

Large-scale simulations of multi-dimensional unsteady turbulent reacting flows with detailed chemistry and transport can be computationally extremely intensive even on distributed computing architectures. With the development of suitable reduced chemical kinetic models, the number of scalar variables to be integrated can be decreased, leading to a significant reduction in the computational time required for the simulation with limited loss of accuracy in the results. A general MATLAB-based automated mechanism reduction procedure is presented to reduce any complex starting mechanism (detailed or skeletal) with minimal human intervention. Based on the application of the quasi steady-state (QSS) approximation for certain chemical species and on the elimination of the fast reaction rates in the mechanism, several comprehensive reduced models, capable of handling different fuels such as C2H4, CH4 and H2, have been developed and thoroughly tested for several combustion problems (ignition, propagation and extinction) and physical conditions (reactant compositions, temperatures, and pressures). A key feature of the present reduction procedure is the explicit solution of the concentrations of the QSS species, needed for the evaluation of the elementary reaction rates. In contrast, previous approaches relied on an implicit solution due to the strong coupling between QSS species, requiring computationally expensive inner iterations. A novel algorithm, based on the definition of a QSS species coupling matrix, is presented to (i) introduce appropriate truncations to the QSS algebraic relations and (ii) identify the optimal sequence for the explicit solution of the concentration of the QSS species. With the automatic generation of the relevant source code, the resulting reduced models can be readily implemented into numerical codes.

Zambon, A. C.↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
NASA NTRS2006

Opposed Jet Burner Approach for Characterizing Flameholding Potentials of Hydrocarbon Scramjet Fuels

Opposed Jet Burner (OJB) tools have been used extensively by the authors to measure Flame Strength (FS) extinction limits of laminar H2/N2 air and (recently) hydrocarbon (HC) air Counterflow Diffusion Flames (CFDFs) at one atm. This paper details normalization of FSs of N2- diluted H2 and HC systems to account for effects of fuel composition, temperature, pressure, jet diameter, inflow Reynolds number, and inflow velocity profile (plug, contoured nozzle; and parabolic, straight tube). Normalized results exemplify a sensitive accurate means of validating, globally, reduced chemical kinetic models at approx. 1 atm and the relatively low temperatures approximating the loss of non-premixed idealized flameholding, e.g., in scramjet combustors. Laminar FS is defined locally as maximum air input velocity, U(sub air), that sustains combustion of a counter-jet of g-fuel at extinction. It uniquely characterizes a fuel. And global axial strain rate at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or (sub t)) can be compared directly with computed extinction limits, determined using either a 1-D Navier Stokes stream-function solution, using detailed transport and finite rate chemistry, or (better yet) a detailed 2-D Navier Stokes numerical simulation. The experimental results define an idealized flameholding reactivity scale that shows wide ranging (50 x) normalized FS s for various vaporized-liquid and gaseous HCs, including, in ascending order: JP-10, methane, JP-7, n-heptane, n-butane, propane, ethane, and ethylene. Results from H2 air produce a unique and exceptionally strong flame that agree within approx. 1% of a recent 2-D numerically simulated FS for a 3 mm tube-OJB. Thus we suggest that experimental FS s and/or FS ratios, for various neat and blended HCs w/ and w/o additives, offer accurate global tests of chemical kinetic models at the Ts and Ps of extinction. In conclusion, we argue the FS approach is more direct and fundamental, for assessing, e.g., idealized scramjet flameholding potentials, than measurements of laminar burning velocity or blowout in a Perfectly Stirred Reactor, because the latter characterize premixed combustion in the absence of aerodynamic strain. And FS directly measures a chemical kinetic characteristic of non-premixed combustion at typical flameholding temperatures. It mimics conditions where gfuels are typically injected into a subsonic flameholding recirculation zone that captures air, where the effects of aerodynamic strain and associated multi-component diffusion become important.

Pellett, Gerald L.↗
NASA NTRS1990

Reduced and simplified chemical kinetics for air dissociation using Computational Singular Perturbation

The Computational Singular Perturbation CSP methods is employed (1) in the modeling of a homogeneous isothermal reacting system and (2) in the numerical simulation of the chemical reactions in a hypersonic flowfield. Reduced and simplified mechanisms are constructed. The solutions obtained on the basis of these approximate mechanisms are shown to be in very good agreement with the exact solution based on the full mechanism. Physically meaningful approximations are derived. It is demonstrated that the deduction of these approximations from CSP is independent of the complexity of the problem and requires no intuition or experience in chemical kinetics.

Goussis, D. A.↗
NASA NTRS2014

Assessment of Reduced-Kinetics Mechanisms for Combustion of Jet Fuel in CFD Applications

A computational effort was undertaken to analyze the details of fluid flow in Lean-Direct Injection (LDI) combustors for next-generation LDI design. The National Combustor Code (NCC) was used to perform reacting flow computations on single-element LDI injector configurations. The feasibility of using a reduced chemical-kinetics approach, which optimizes the reaction rates and species to model the emissions characteristics typical of lean-burning gas-turbine combustors, was assessed. The assessments were performed with Reynolds- Averaged Navier-Stokes (RANS) and Time-Filtered Navier Stokes (TFNS) time-integration, with a Lagrangian spray model with the NCC code. The NCC predictions for EINOx and combustor exit temperature were compared with experimental data for two different single-element LDI injector configurations, with 60deg and 45deg axially swept swirler vanes. The effects of turbulence-chemistry interaction on the predicted flow in a typical LDI combustor were studied with detailed comparisons of NCC TFNS with experimental data.

Ajmani, Kumud↗
DOE OSTI2025

Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks

Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.

autoencoder neural networks↗
NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗
NASA NTRS1995

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗
NASA NTRS1992

Asymptotic analysis with reduced chemistry for the burning of n-heptane droplets

The method of rate-ratio asymptotics is used with reduced chemistry to analyze the flame structure and extinction of an isolated n-heptane droplet burning under quasisteady, spherically symmetrical conditions. The outer transport zones are described by the classical Burke-Schumann solution. The inner reaction zone consists of a thin layer, on the rich side of the flame, where the fuel is consumed, and on the lean side, a broader but still asymptotically thin oxidation layer, where H2 and CO are consumed. Special attention is given to differences in predictions of extinction conditions, caused by different chemical-kinetic approximations in the reduced chemistry, including fuel-chemistry effects through molecules containing more than one carbon atom. From the analysis, extinction diameters for n-heptane droplets are estimated for different pressures and ambient oxygen concentrations. The results show that extinction diameters are extremely sensitive to the number of radicals consumed in breaking down each fuel molecule.

Card, J. M.↗
NASA NTRS2016

The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure-temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

planets and satellites↗
NASA NTRS2012

AdapChem

AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.

Oluwole, Oluwayemisi O.↗
NASA NTRS2018

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

This work presents a general framework for model reduction of non-equilibrium energy transfer and dissociation processes. The multi-group maximum-entropy method is coupled with the quasi-classical trajectory method to directly construct a reduced order model for chemical non-equilibrium. Kinetic data is calculated by applying detailed balance at a microscopic level, overcoming the limitation of calculating recombination kinetic data. This approach enables the construction of a reduced order model for kinetics which bypasses the need to compute state-to-state kinetic data. This physics based reduced order model ensures that an equilibrium distribution is reached given infinite time, while allowing for non-equilibrium distributions during the relaxation and dissociation processes. A proof-of-concept test case demonstrates the applicability of this model by comparison with state-to-state kinetic data for the N2(Chi (sup 1) Sigma ((sub g), (sup plus sign)) minus N ((exp 4) S (sub u)) system with excellent agreement.

Maximum Entropy↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
NASA NTRS2011

A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane

The reduction of elementary or skeletal oxidation kinetics to a subgroup of tractable reactions for inclusion in turbulent combustion codes has been the subject of numerous studies. The skeletal mechanism is obtained from the elementary mechanism by removing from it reactions that are considered negligible for the intent of the specific study considered. As of now, there are many chemical reduction methodologies. A methodology for deriving a reduced kinetic mechanism for alkane oxidation is described and applied to n-heptane. The model is based on partitioning the species of the skeletal kinetic mechanism into lights, defined as those having a carbon number smaller than 3, and heavies, which are the complement of the species ensemble. For modeling purposes, the heavy species are mathematically decomposed into constituents, which are similar but not identical to groups in the group additivity theory. From analysis of the LLNL (Lawrence Livermore National Laboratory) skeletal mechanism in conjunction with CHEMKIN II, it is shown that a similarity variable can be formed such that the appropriately non-dimensionalized global constituent molar density exhibits a self-similar behavior over a very wide range of equivalence ratios, initial pressures and initial temperatures that is of interest for predicting n-heptane oxidation. Furthermore, the oxygen and water molar densities are shown to display a quasi-linear behavior with respect to the similarity variable. The light species ensemble is partitioned into quasi-steady and unsteady species. The reduced model is based on concepts consistent with those of Large Eddy Simulation (LES) in which functional forms are used to replace the small scales eliminated through filtering of the governing equations; in LES, these small scales are unimportant as far as the overwhelming part of dynamic energy is concerned. Here, the scales thought unimportant for recovering the thermodynamic energy are removed. The concept is tested by using tabular information from the LLNL skeletal mechanism in conjunction with CHEMKIN II utilized as surrogate ideal functions replacing the necessary functional forms. The test reveals that the similarity concept is indeed justified and that the combustion temperature is well predicted, but that the ignition time is over-predicted, a fact traced to neglecting a detailed description of the processes leading to the heavies chemical decomposition. To palliate this deficiency, functional modeling is incorporated into this conceptual reduction in addition to the modeling the evolution of the global constituent molar density, the enthalpy evolution of the heavies, the contribution to the reaction rate of the unsteady lights from other light species and from the heavies, the molar density evolution of oxygen and water, and the mole fractions of the quasisteady light species. The model is compact in that there are only nine species-related progress variables. Results are presented showing the performance of the model for predicting the temperature and species evolution. The model reproduces the ignition time over a wide range of equivalence ratios, initial pressure, and initial temperature.

Harstad, Kenneth G.↗
DOE OSTI2024

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗
NASA NTRS2007

Combustion

An overview of the emissions related research being conducted as part of the Fundamental Aeronautics Subsonics Fixed Wing Project is presented. The overview includes project metrics, milestones, and descriptions of major research areas. The overview also includes information on some of the emissions research being conducted under NASA Research Announcements. Objective: Development of comprehensive detailed and reduced kinetic mechanisms of jet fuels for chemically-reacting flow modeling. Scientific Challenges: 1) Developing experimental facilities capable of handling higher hydrocarbons and providing benchmark combustion data. 2) Determining and understanding ignition and combustion characteristics, such as laminar flame speeds, extinction stretch rates, and autoignition delays, of jet fuels and hydrocarbons relevant to jet surrogates. 3) Developing comprehensive kinetic models for jet fuels.

Bulzan, Dan↗