SEARCH · Search NASA

Results for “reaction rate equations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Maximum Sample Temperature for Mars Sample Return: A Historical Perspective

Since the first Mars Sample Return (MSR) report published by the Jet Propulsion Laboratory (JPL) in 1974 [1], a series of panels, reports, and white papers have recognized the importance of sample temperature and offered an informed sample maximum temperature (henceforth SMT) limit for returning martian samples to Earth. The Mars Sample Handling and Requirements Panel (MSHARP, 1999) stated that "[t]he main issue in sample preservation is temperature" [2]. More recently, the Mars Exploration Program Analysis Group (MEPAG)'s "Science Priorities for Mars Sample Return" report (2008), declared that "[s]ignificant loss, particularly to biological studies, occurs if samples reach +50C for three hours", whereby "scientific objectives related to life goals could be seriously compromised" [3]. By contrast, the Mars 2020 mission has adopted a SMT of +60C as spelled out in Beaty et al., 2016 [4]. Samples will be collected and then deposited on the surface in sealed tubes for possible retrieval and return to Earth. Beaty et al. [4] calculates that the samples will experience maximum temperatures of ~+30 to +60C, depending on latitude. At present, there is no mission requirement for the measurement/data logging of sample temperature during this period. We will explore the history of martian SMTs, as they have been recorded since 1974 [1], effectively representing input across multiple generations of Mars scientists. Ten separate publications present SMTs for MSR samples [1-10]. One report [10] is for a mission concept specifically designed to exclude life detection investigations, and recommended an SMT of 50C. Another did not specify a temperature, recommending "Mars ambient temperature" [5]. Of the remaining eight, SMTs are given as: -30C [1], -20C [3], 60C [4], -73 to 41C depending on sample type [6], -40C [7], -43 to 13C depending on type [2,8], and -33C [9]. If we restrict the temperatures to samples highlighted in the Mars 2020 mission goals, i.e. organics-bearing and sedimentary rocks, then the average SMT is -28+/-39C (n=8). Applying a Dixon's Q Test at P=0.05 (two-tailed), the 60C SMT [4] fails with Q=0.602 versus Qcrit=0.526. Excluding the outlier produces an average SMT of -40+/-17C (n=7). Therefore, the average SMT expressed by the Mars science community over the past 44 years (two generations) is a sample temperature no greater than -40C. The difference in chemical reaction rates between this average SMT and Beaty et al [4] can be estimated using the Arrhenius equation. Assuming a generic chemical reaction with an activation energy of 50 kJ/mol and a pre-exponential factor invariant with temperature, this reaction will proceed 2300x faster at 60C than at -40C. To illustrate the effects of the increased reaction rate, consider 10 ppb of alanine in a Mars 2020 cache, and assume that it becomes unmeasurable if it degrades to 1 ppb, as per the Mars 2020 Organic Contamination Panel contamination limits [11]. If we illustrate the effect with an arbitrary degradation rate such that the alanine will become undetectable in ten years at -40C, then the same 10 ppb alanine degrades beyond detectability in only 38 days at 60C. Further research is required to quantify expected analyte losses in the cached samples due to thermal processing.

Fries, Marc↗

DOE OSTI2020

Formation of phenanthrene via H–assisted isomerization of 2–ethynylbiphenyl produced in the reaction of phenyl with phenylacetylene

Model chemistry G3(MP2,CC)//B3LYP/6–311G(d,p) calculations of the potential energy surface for the reaction of phenyl radical (C 6 H 5 ) with phenylacetylene (C 8 H 6 ) have been carried out and combined with Rice–Ramsperger–Kassel–Marcus/Master Equation calculations of temperature– and pressure–dependent rate constants. The results showed that the reaction can serve as a viable source for the formation of phenanthrene via an indirect route involving a primary reaction of phenyl addition to the ortho carbon in the ring of phenylacetylene and H elimination producing 2–ethynylbiphenyl followed by secondary H–assisted isomerization of 2–ethynylbiphenyl to phenanthrene. In the secondary reaction, the H atom adds to the α carbon of the ethynyl side chain, then a six–member ring closure takes place followed by aromatization via an H loss. The channel of H addition to the side chain of 2–ethynylbiphenyl appears to be much faster than H addition to the ortho carbon in the ethynyl–substituted ring leading back to the initial C 6 H 5 + C 8 H 6 reactants. Rate constants for the primary C 6 H 5 + C 8 H 6 ⇌2–ethynylbiphenyl (p1) + H and secondary p1 + H⇌phenanthrene (p2) + H reactions have been computed in the temperature range of 500–2500 K at pressures of 30 Torr, 1, 10, and 100 atm and fitted to modified Arrhenius expressions. Here, the suggested kinetic scheme and rate constants are proposed as a prototype for the modeling of the growth of polycyclic aromatic hydrocarbons via the phenyl addition–dehydrocyclization (PAC) mechanism involving an addition of a PAH radical to an ethynyl–substituted PAH molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Direct Discontinuous Galerkin methods for the reacting multi-component flow equations

The Direct Discontinuous Galerkin (DDG (Liu and Yan, 2008)) method and a counterpart with Interface Correction (DDGIC (Danis and Yan, 2022)) are extended to compute diffusion terms that arise when solving the compressible multi-component flow equations in thermochemical nonequilibrium. Thermodynamic properties, transport properties, chemical reaction rates, and energy exchange terms are computed using Mutation++ (Scoggins et al., 2020). The DG method is applied on unstructured grids, where the accuracy and convergence rates can be sensitive to the numerical method chosen for parabolic terms. A method for determining the homogeneity tensor of the flow equations required for DDGIC is shown. The convergence properties of the DDG methods are studied and compared to the Interior Penalty (IP) method. A number of numerical experiments are conducted to assess the accuracy and performance of the method. The numerical results and convergence studies indicate that DDG and DDGIC provide accurate solutions and perform well for general flows in thermochemical nonequilibrium.

Diffusion↗

NASA NTRS2014

Calculation of State Specific Rate Coefficients for Non-Equilibrium Hypersonics Applications: from H(Psi) = E(Psi) to k(T) = A *exp(-E(sub a)/RT)

Development of High-Fidelity Physics-Based Models to describe hypersonic flight through the atmospheres of Earth and Mars is underway at NASA Ames Research Center. The goal is to construct chemistry models of the collisional and radiative processes that occur in the bow shock and boundary layers of spacecraft during atmospheric entry that are free of empiricism. In this talk I will discuss our philosophy and describe some of our progress. Topics to be covered include thermochemistry, internal energy relaxation, collisional dissociation and radiative emission and absorption. For this work we start by solving the Schrodinger equation to obtain accurate interaction potentials and radiative properties. Then we invoke classical mechanics to compute state-specific heavy particle collision cross sections and reaction rate coefficients. Finally, phenomenological rate coefficients and relaxation times are determined from master equation solutions.

aerothermodynamics↗

DOE OSTI2023

Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software

The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

A transport equation for the scalar dissipation in reacting flows with variable density: First results

Although the different regimes of premixed combustion are not well defined, most of the recent developments in turbulent combustion modeling are led in the so-called flamelet regime. The goal of these models is to give a realistic expression to the mean reaction rate (w). Several methods can be used to estimate (w). Bray and coworkers (Libby & Bray 1980, Bray 1985, Bray & Libby 1986) express the instantaneous reaction rate by means of a flamelet library and a frequency which describes the local interaction between the laminar flamelets and the turbulent flowfield. In another way, the mean reaction rate can be directly connected to the flame surface density (Sigma). This quantity can be given by the transport equation of the coherent flame model initially proposed by Marble & Broadwell 1977 and developed elsewhere. The mean reaction rate, (w), can also be estimated thanks to the evolution of an arbitrary scalar field G(x, t) = G(sub O) which represents the flame sheet. G(x, t) is obtained from the G-equation proposed by Williams 1985, Kerstein et al. 1988 and Peters 1993. Another possibility proposed in a recent study by Mantel & Borghi 1991, where a transport equation for the mean dissipation rate (epsilon(sub c)) of the progress variable c is used to determine (w). In their model, Mantel & Borghi 1991 considered a medium with constant density and constant diffusivity in the determination of the transport equation for (epsilon(sub c)). A comparison of different flamelet models made by Duclos et al. 1993 shows the realistic behavior of this model even in the case of constant density. Our objective in this present report is to present preliminary results on the study of this equation in the case of variable density and variable diffusivity. Assumptions of constant pressure and a Lewis number equal to unity allow us to significantly simplify the equation. A systematic order of magnitude analysis based on adequate scale relations is performed on each term of the equation. As in the case of constant density and constant diffusivity, the effects of stretching of the scalar field by the turbulent strain field, of local curvature, and of chemical reactions are predominant. In this preliminary work, we suggest closure models for certain terms, which will be validated after comparisons with DNS data.

Mantel, T.↗

NASA NTRS1977

Model for premixed turbulent flames

Flame speed and mean degree of reaction profiles are determined for an ensemble mean one-dimensional premixed turbulent flame. Elementary closure models are employed to treat turbulence closure problems arising from the nonlinear convective and chemical reaction rate terms in the time-dependent fuel conservation equation; a constant eddy diffusivity is useed to model correlation terms derived from the convective term. An expression developed for fuel consumption is solved to obtain the turbulent flame speed and the variation of mean fuel mass fraction across the flame. Comparison of the numerical analysis with experimental data indicates the importance of the quenching of the combustion reaction in high-wavenumber eddies.

Gouldin, F. C.↗

NASA NTRS1988

Numerical modeling of chemical vapor deposition (CVD) in a horizontal reactor

In the present numerical prediction of the deposition rate of silicon from silane in a CVD process, the conservation equations for mass, momentum, energy, and chemical species are solved on a staggered grid using the SIMPLE algorithm, while the rate of chemical reactions in the gas phase and on the susceptor surface is obtained from an Arrhenius rate equation. Predicted deposition rates as a function of position along the susceptor with and without the gas phase chemical reaction are compared with the available experimental and numerical data; agreement is excellent except at the leading edge of the susceptor, where the deposition rate is overpredicted.

Sheikholeslami, M. Z.↗

DOE OSTI2023

Learning effective stochastic differential equations from microscopic simulations: Linking stochastic numerics to deep learning

We identify effective stochastic differential equations (SDEs) for coarse observables of fine-grained particle- or agent-based simulations; these SDEs then provide useful coarse surrogate models of the fine scale dynamics. We approximate the drift and diffusivity functions in these effective SDEs through neural networks, which can be thought of as effective stochastic ResNets. The loss function is inspired by, and embodies, the structure of established stochastic numerical integrators (here, Euler–Maruyama and Milstein); our approximations can thus benefit from backward error analysis of these underlying numerical schemes. They also lend themselves naturally to “physics-informed” gray-box identification when approximate coarse models, such as mean field equations, are available. Existing numerical integration schemes for Langevin-type equations and for stochastic partial differential equations can also be used for training; we demonstrate this on a stochastically forced oscillator and the stochastic wave equation. Our approach does not require long trajectories, works on scattered snapshot data, and is designed to naturally handle different time steps per snapshot. We consider both the case where the coarse collective observables are known in advance, as well as the case where they must be found in a data-driven manner.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Models to Incorporate Reaction Mechanisms into DG-OSPREY: Fixed-Bed Simulations for Organic Iodides Capture Using Ag 0 Z

As one of the most potent radioisotopes released during spent nuclear fuel reprocessing, 129 I is strictly regulated and must be removed before discharge. Organic iodides (primarily alkyl iodides with different chain lengths, i.e., CH 3 I, C 4 H 9 I, and C 12 H 25 I) comprise ~2% of the total iodine in the reprocessing off-gases and are primarily present in vessel off-gas (VOG). Reduced silver mordenite (Ag 0 Z) is predominantly considered for the removal of radioiodine; however, its capture performance and underlying interaction processes with long-chain organic iodides are not fully understood. Two major tasks were accomplished in this study. First, to improve upon the previous experimental studies where Ag 0 Z was used to capture CH 3 I, C 4 H 9 I, and C 12 H 25 I at different concentrations, we comprehensively investigated the corresponding capture mechanisms by characterizing fully loaded Ag 0 Z samples. Second, computational codes were implemented to perform fixed-bed simulations that account for transport and reaction mechanisms. Scanning electron microscopy with energy dispersive X-ray analysis (SEM-EDX), powder X-ray diffraction (PXRD), UV-visible diffuse reflectance spectroscopy (UV-vis DRS), and thermogravimetric analysis (TGA) were conducted on Ag 0 Z samples that are saturated with I 2 , CH 3 I, C 4 H 9 I, and C 12 H 25 I), respectively. Results indicate that AgI is the predominant adsorption product regardless of the adsorbed iodine species, yet alkyl iodides with different carbon chain lengths may have different compositions of α- and γ-AgI. Synchrotron pair distribution function (PDF) measurements and TGA coupled with a Fourier transformed infrared detector (TGA FTIR) have been performed, and experimental data are currently analyzed. Results are expected to provide further insights into the adsorption mechanisms. The fixed-bed capture performance for CH 3 I was successfully simulated by solving mathematical equations that describe the underlying transport processes and adsorption reactions. The computational framework, Catalytic After Treatment System (CATS), that was originated in our research group, was used to solve the governing equations. Kinetic parameters, including the pore diffusivity and reaction rate constant were obtained by optimization techniques using data from thin-bed experiments performed at Oak Ridge National Laboratory. The performance of a deep bed predicted using the optimized parameters showed promising agreement with the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1990

Enhanced thermochemical nonequilibrium computations of flow around the aeroassist flight experiment vehicle

A three-dimensional explicit, finite-rate, shock-capturing numerical algorithm is used to calculate thermochemical nonequilibrium flowfields about the Aeroassist Flight Experiment vehicle at one of its flight trajectory points. The full Navier-Stokes equations and an eleven species chemical model with the latest reaction rates are incorporated into the code. Results are compared against experimental data and other numerical solutions. The effects of changes to the physical model on the stagnation line and base region flow and surface quantities on the forebody are investigated.

Palmer, Grant↗

NASA NTRS1994

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗

DOE OSTI2024

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2014

State-to-State Internal Energy Relaxation Following the Quantum-Kinetic Model in DSMC

A new model for chemical reactions, the Quantum-Kinetic (Q-K) model of Bird, has recently been introduced that does not depend on macroscopic rate equations or values of local flow field data. Subsequently, the Q-K model has been extended to include reactions involving charged species and electronic energy level transitions. Although this is a phenomenological model, it has been shown to accurately reproduce both equilibrium and non-equilibrium reaction rates. The usefulness of this model becomes clear as local flow conditions either exceed the conditions used to build previous models or when they depart from an equilibrium distribution. Presently, the applicability of the relaxation technique is investigated for the vibrational internal energy mode. The Forced Harmonic Oscillator (FHO) theory for vibrational energy level transitions is combined with the Q-K energy level transition model to accurately reproduce energy level transitions at a reduced computational cost compared to the older FHO models.

Liechty, Derek S.↗

DOE OSTI2025

Not-quite-transcendental Functions for Logarithmic Interpolation of Tabulated Data

From tabulated nuclear and degenerate equations of state to photon and neutrino opacities and nuclear reaction rates, tabulated data is ubiquitous in computational astrophysics. The dynamic range that must be covered by these tables typically spans many orders of magnitude. Here we present a novel strategy for accurately and performantly interpolating tabulated data that spans these large dynamic ranges. We demonstrate the efficacy of this strategy in tabulated lookups for nuclear and terrestrial equations of state. We show that this strategy is a faster drop-in replacement for linear interpolation of logarithmic grids.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

Kinetics of hexacelsian to celsian phase transformation in SrAl2Si2O8

The kinetics of hexacelsian to celsian phase transformation in SrAl2Si2O8 have been investigated. Phase pure hexacelsian was prepared by heat treatment of glass flakes at 990 C for 10 h. Bulk hexacelsian was isothermally heat treated at 1026, 1050, 1100, 1152, and 1200 C for various times. The amounts of monoclinic celsian formed were determined using quantitative X-ray diffraction. Values of reaction rate constant, k, at various temperatures were evaluated from the Avrami equation. The Avrami parameter was determined to be 1.1, suggesting a diffusionless, one-dimensional transformation mechanism. From the temperature dependence of k, the activation energy for this reaction was evaluated to be 527 plus or minus 50 kJ/mole (126 plus or minus 12 kcal/mole). This value is consistent with a mechanism involving the transformation of the layered hexacelsian structure to a three-dimensional network celsian structure which necessitates breaking of the strongest bonds, the Si-O bonds.

Bansal, Narottam P.↗

NASA NTRS1993

Kinetics of hexacelsian-to-celsian phase transformation in SrAl2Si2O8

Bansal, Narottam P.↗

NASA NTRS2005

Control of Combustion-Instabilities Through Various Passive Devices

It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.

Frendi, Kader↗