Charged Species Concentration in Combusting Mixtures Using Equilibrium Chemistry.
Ionization in flames has been a topic of research since the late 1940s [1-10]. Early investigations in this area were motivated primarily by defense applications such as the need to define the size and extent of radar targets, predict electromagnetic wave attenuation in propulsion and re-entry systems, and understand various communication-related problems. Another important application of ionized gas flows was in the area of power generation. Magnetohydrodynamic (MHD) power generation received considerable attention in the 1960s. MHD generators produce electricity directly from a fast moving stream of ionized gases without the need for any moving mechanical parts such as turbines and rotary generators. One of the main technical challenges in operating MHD generators is the creation and sustenance of the ions in the working fluid. Candidate working fluids included flow streams derived from combustion, noble gases such as argon and gas streams seeded with alkali metal vapors such as Potassium and Cesium. The design and development of these MHD generators required the accurate computation of the charged species concentrations to determine the conductivity of the working fluid.) In recent years, research on ionization in flames has been spurred by the need for advanced ignition methods for engineering applications. Ultra-lean combustion in internal combustion engines and high-speed flows in aerospace applications (ramjets and scramjets) are two notable examples [1, 2]. Internal combustion engines represent applications characterized by low temperatures and high pressures (T ~ 600-700 K, P ~ 15-40 atm at ignition) and low speeds (gas velocity 10-50 m/s), whereas aerospace applications represent a case of moderate temperature, low pressure (P ~ 0.1-1 atm) and high velocity (V ~ 500 m/s). Given the emphasis on reduction of emissions in the power generation and transportation sectors, the role of ions in soot formation is another important area of study. Soot formation is believed to be influenced by charged species in fuel-rich hydrocarbon-air combustion. The main precursor for soot formation is assumed to be the chemi-ion, [C.sub.3][H.sub.3.sup.+], which is the dominant ion in fuel-rich hydrocarbon flames. Understanding the process of initiation and development of the combustion kernel plays a key role in design, development, and optimization of engineering applications using advanced ignition methods. Active radicals and charged species (electrons and ions) are believed to be important in the early flame development phase. Determining the identity and concentration of important species using both, experimental and numerical techniques, has been a topic of research since the late nineteen forties [5-11]. More recent studies on the topic can be found in [12-19]. Charged species concentrations in reacting mixtures can be evaluated using equilibrium chemistry or finite-rate chemistry (FRC) assumptions. Most of the numerical studies aimed at evaluating charged species concentrations in flames use finite-rate chemistry (see [12-16] and references therein) in low-pressure, laboratory flames using simple hydrocarbons such as methane. In power generating combustion devices such as internal combustion engines and gas turbines, combustion takes place under high pressure conditions with complex fuels and hence the mechanisms used to study combustion of simple hydrocarbons at low pressures would likely have to be extended considerably to produce reliable and accurate results.
In some of the earliest studies by Calcote and his coworkers, the concentration of positive ions in the combustion zone of a laboratory hydrocarbon flame (P = 1 atm) was estimated to be about [10.sup.12] ions/cc using a simplistic mechanical model [5, 6] (although later experiments have shown the ion concentrations to be in the range of 109-1010 ions/cc [7, 11].) Calcote and King  sought to explain the positive ion concentration deduced from experimental measurements of current in stoichiometric propane-air flames by using equilibrium assumptions. They assumed that the combustion mixture consisted of a set of neutral species and an adiabatic flame temperature of 2200 K for a near-stoichiometric propane-air mixture at a pressure of 1 atmosphere. They computed the ion concentration under thermal equilibrium assumptions using the Saha equation (this technique will henceforth be referred to as the Calcote method). The following neutral species were considered in their computations: [H.sub.2], [N.sub.2], CO, C[O.sub.2], H, OH, O, [H.sub.2]O, [O.sub.2], NO, Na, and C. Only reactions of the kind M = [M.sup.+] + [e.sup.-] were considered. Furthermore, a given ion concentration was computed assuming that only a single ionic species was present in the combustion mixture. For instance, to compute the concentration of [H.sup.+] ions, they assumed that the combustion products contained only [H.sub.2], [N.sub.2], CO, C[O.sup.2], H, OH, O, [H.sub.2]O, [O.sub.2], NO, Na, C, [H.sup.+], and [e.sup.-]. Calcote and King reported that at 2200 K, N[O.sup.+] has the highest ion concentration at 1.7 x [10.sup.7] ions/cc (excluding [Na.sup.+] and [C.sup.+]), as shown in Table 1 of . The concentration of the dominant ionic specie, N[O.sup.+], is well below the value of [10.sup.12] ions/cc estimated from experiments on stoichiometric propane-air flames. Given the large discrepancy between ion concentrations estimated from experimental data and equilibrium calculations, Calcote and King concluded that "ions are produced not by thermal processes but by chemi-ionization or events similar to those using abnormal electronic excitation in flames. The results have also been taken as evidence that thermal equilibrium does not exist between the species in the flame front" .
The large discrepancy between the maximum ion concentration reported by Calcote's computations and the estimates inferred from experimental data could stem from the restrictive assumptions and computational methodology used in Calcote's calculations. Specifically, Calcote and King considered only positive ions and electrons in their equilibrium calculations (produced by thermal ionization). The presence of negative ions such as O[H.sup.-], [O.sup.-], or [O.sub.2.sup.-] and positive ions such as [H.sub.3][O.sup.+] were ignored. While Calcote and King did conclude that thermal equilibrium does not exist between the ions and their parent species, and that chemi-ions could be responsible for the higher ion concentrations deduced from experiments, they did not extend their computations to include these effects. Goodings et al.  showed that one of the dominant charged species in rich C[H.sub.4]-[O.sup.2] flames is [H.sub.3][O.sup.+] and that its measured value is identical to the expected equilibrium value. Goodings et al. also showed that the measured concentration of [H.sub.3][O.sup.+] is ~1.4 x [10.sup.10] ions/cc, which is two orders of magnitude less than the ~[10.sup.12] ions/cc estimated by Calcote. Additionally, Goodings et al. showed that the concentration of the second most dominant species, CH[O.sup.+], is about three orders of magnitude less than that of [H.sub.3][O.sup.+] (see Table 1 in ). The importance of [H.sub.3][O.sup.+] in hydrocarbon-air flames is well known and has been reported in more recent studies using finite-rate chemical kinetic calculations. Similarly, recent studies have shown that the ion current in a constant volume chamber and HCCI engines can be explained well by using equilibrium assumptions and that [H.sub.3][O.sup.+] is indeed the dominant positive ion [17, 18].
While there are a considerably large number of studies using finite-rate chemistry to investigate charged species concentrations in laboratory and low-pressure flames, there is relatively little work done on the use of equilibrium chemistry methods to evaluate charged species composition in reacting mixtures. In this work, a more in-depth study of the equilibrium concentration in stoichiometric hydrocarbon flames is conducted by including the effects of chemi-ions and elevated electron temperatures, thus complementing Calcote's initial work. Revisiting Calcote's calculations to determine the appropriateness of using equilibrium assumptions for computing charged species concentrations in reacting mixtures can also benefit the development of fast and robust design/optimization tools for modern applications of plasma-based combustion systems which typically operate in the lean/stoichiometric regime (equivalence ratio < 1). Use of equilibrium chemistry computations for determining the concentration of combustion products is advantageous to using finite-rate chemistry calculations. Finite-rate chemistry calculations require the use of an appropriate mechanism that can be large (a few hundred reactions and species) even for simple hydrocarbons and can present serious computational challenges. For instance, Prager et al.  used a mechanism consisting of 275 reactions and 49 species to model lean/stoichiometric methane-air flames while Starik and Titova  proposed an ionic mechanism containing 392 reactions and 59 species to predict the ion concentrations in a perfectly stirred reactor. Use of such large mechanisms requires considerable computational resources, in terms of both hardware and wall-time, and hence precludes their use in multi-dimensional reacting flow simulations for design/developmental studies. Equilibrium chemistry computations assume that the reactions occur at time scales much shorter than the flow (mixing and/or diffusion) time scale; and the concentration of the final products depends only on the temperature, pressure, and initial elemental mixture composition of the system. Approximately 30 species (including ions and electrons) are usually adequate for describing the combustion products required for engineering design and analyses of most fuel-air combustion systems. The concentrations of species in the mixture are computed using thermophysical properties (enthalpy, entropy, Gibbs free energy) of the system. These thermophysical properties can be computed using constants that are easily available (for instance, in the CHEMKIN and NASA-CEA databases). Finite-rate chemistry computations, on the other hand, require reaction rate constants for each elementary reaction in the mechanism. These reaction rates constants are not readily available, especially for elementary reactions involving electrons and charges species. Equilibrium chemistry computations, especially for hydrocarbonair combustions, are also "extendible" in that equilibrium concentrations for any hydrocarbon-air combustion can be obtained by merely replacing the thermophysical constants of one hydrocarbon with another with all other species remaining the same. Finite-rate chemistry computations, on the other hand, would require a different mechanism with the inclusion of more intermediate species to compute the species concentrations. Given these considerations, equilibrium calculations can provide valuable estimates of charged species concentrations in plasma-based engineering systems characterized by large Damkohler numbers at greatly reduced computational cost and effort.
To this end, we conducted a systematic study of the problem of computing equilibrium concentrations of various neutrals and ions in combusting stoichiometric hydrocarbon-air mixtures under various conditions of temperature, pressure, and thermal disequilibrium. While the identity of dominant charged species in hydrocarbon flames is well-established in laboratory and low-pressure flames using finite-rate chemistry calculations, our goal is to answer the following questions with regard to the equilibrium calculations conducted by Calcote and King and assess the suitability of the use of equilibrium chemistry computations for engineering applications in power generating equipment.
(1) Can a more accurate method to compute equilibrium concentrations of all ions and neutrals considered in Calcote's original analyses yield a value of the magnitude of the dominant charged specie which is closer to ion concentrations estimated from experimental data?
(2) What is the impact of mixture composition (considering a different set of neutral and charged species as combustion products) on the magnitude of the dominant charged specie?
(3) What is the impact of elevated electron temperature (thermal non-equilibrium) on the equilibrium concentration of charged species in hydrocarbon combustion?
Our computations seek to determine whether equilibrium chemistry calculations conducted with the inclusion of [H.sub.3][O.sup.+] ions and the effect of elevated electron temperatures can indeed yield results consistent with experimental observations, as postulated by Calcote and King. Furthermore, the validity of equilibrium computations at higher pressures and thermal disequilibrium is also investigated. To the best of the author's knowledge, such a systematic study has not been conducted.
This paper is organized as follows. Section 2 describes the method of solution, and Section 3 discusses the results. Section 4 briefly summarizes the main findings of this work.
2. Method of Solution
The equilibrium concentrations of charged and neutral species were obtained by using the Chemical Equilibrium Analysis (CEA) code. The CEA code computes equilibrium concentrations by minimizing the Gibbs free energy of the system containing all the ions and neutral species assumed to be present in the system. Details of the computational methodology can be obtained from . The CEA code is well-validated and widely used for equilibrium computations and hence used in this study. The neutral and charged species originally considered by Calcote (namely, [H.sub.2], [N.sub.2], CO, C[O.sub.2], H, OH, O, [H.sub.2]O, [O.sub.2], NO, Na, C, [H.sub.2.sup.+], [N.sub.2.sup.+], C[O.sup.+] C[O.sub.2.sup.+], [H.sup.+], O[H.sup.+], [O.sup.+], [H.sub.2][O.sup.+], [O.sub.2.sup.+], N[O.sup.+],[ Na.sup.+], [C.sup.+], and [e.sup.-]) were used to obtain the equilibrium concentrations of a stoichiometric propane-air mixture at T = 2200 K and P = 1 atm using the CEA code. This set of neutral and charged species used by Calcote will be referred to as mixture I henceforth. The equilibrium concentrations obtained by using the CEA code were compared with those obtained by using the Calcote method . Equilibrium computations of two other propane-air mixtures consisting of different sets of ions and neutrals than those used by Calcote were conducted to assess the impact of mixture composition on the dominant charged species concentration. The first propane-air mixture considered in this study, referred to as mixture II, had 18 neutral and 7 charged species. The neutral species in mixture II were [C.sub.3][H.sub.8] (propane), [O.sub.2], C[O.sub.2], [H.sub.2]O, [N.sub.2], N, O, NO, OH, H, [N.sub.2]O, CO, [H.sub.2], N[O.sub.2], H[O.sub.2], [C.sub.2][H.sub.2], C, and CH. The seven charged species in mixture II were CH[O.sup.+], [e.sup.-], [H.sub.3][O.sup.+], N[O.sup.+], [O.sub.2.sup.-], [O.sup.-], and O[H.sup.-]. The neutral and charged species assumed to be present in mixture II have shown good agreement with measured ion currents in constant volume cylinders and internal combustion engines, as discussed in [17, 18] and hence used in this study. Validation of the ion current predictions in a constant volume methane-air mixture using mixture II (as discussed in ) is shown in Section 3.4 for the benefit of the reader.
For the sake of completeness, equilibrium composition of a stoichiometric propane-air mixture including all the species in mixture I and mixture II, termed as mixture, III was also studied. The species present in mixture III were [C.sub.3][H.sub.8] [H.sub.2] [N.sub.2] CO C[O.sub.2] H OH O [H.sub.2]O [O.sub.2] CH NO [H.sub.2.sup.+] [N.sub.2.sup.+] C[O.sup.+] C[O.sub.2.sup.+] [H.sup.+] O[H.sup.+] [O.sup.+] [H.sub.2][O.sup.+] [O.sub.2.sup.+] C[H.sup.+] N[O.sup.+] [e.sup.-] N [N.sub.2]O N[O.sub.2] H[O.sub.2] [C.sub.2][H.sub.2], C HC[O.sup.+] [H.sub.3][O.sup.+] [O.sub.2.sup.-] [O.sup.-] O[H.sup.-]. Positively charged species, namely, [H.sub.2.sup.+] [N.sub.2.sup.+] C[O.sup.+] C[O.sub.2.sup.+] [H.sup.+] O[H.sup.+] [O.sup.+] [H.sub.2][O.sup.+] [O.sub.2.sup.+] C[H.sup.+] from Calcote's original calculations, were included in this mixture. The effect of elevated electron temperature was studied by using the method of equilibrium constants as described in detail in Section 3.4.
3. Results and Discussions
This section discusses the equilibrium concentration of ions and neutrals in a stoichiometric premixed propane-air mixture using different computations methods (CEA and Calcote), mixture compositions (mixture I, mixture II, and mixture III), and electron temperatures. Additionally, the effects of gas temperature and pressure for engine-like conditions are also discussed.
3.1. Comparison of CEA Calculations with the Calcote Method (Code Validation). There are no detailed experiments reporting the temperature and concentrations of ions and neutral species in stoichiometric propane-air flames at ambient pressure (Calcote's conditions) and hence it is not possible to compare the numerical results presented by Calcote (see Table 1 in ) with experimental data. For the sake of validation, equilibrium concentrations were computed using the Calcote method with the CEA code. Following the Calcote method, the equilibrium mixture composition of the most dominant ion, N[O.sup.+], was computed considering only N[O.sup.+] and [e.sup.-] as the charged species, along with other neutral species considered by Calcote ([H.sub.2], [N.sub.2], CO, C[O.sub.2], H, OH, O, [H.sub.2]O, [O.sub.2], and NO). These computations were conducted with a stoichiometric propane-air mixture at 2200 K and 1 atm, as reported in . Table 1 shows a comparison of the important charged and neutral species. The equilibrium concentrations computed with the CEA code compares very well with those reported in  and hence serves as a validation of the computational methodology used in this work. It was verified that the concentrations of other charged species (such as [H.sub.2][O.sup.+], [H.sub.2.sup.+] etc.) were the same as those reported in  using the Calcote method.
3.2. Equilibrium Computations considering All of Calcote's Ions and Neutrals Simultaneously. Equilibrium concentrations of a stoichiometric propane-air mixture at T = 2200 K and P = 1 atm were obtained by minimizing the Gibbs free energy of a system consisting of all the ions and neutrals considered in Calcote's study (mixture I). This method of computing the equilibrium concentrations is more accurate than the Calcote method. Table 2 shows a comparison between the equilibrium concentrations of the dominant ions using the Calcote method and the CEA code considering all species together. It is seen that the peak concentration of the dominant ion, namely, N[O.sup.+] predicted by Calcote is nearly the same as that predicted by the CEA code. Moreover, both methods predict the same four most dominant ions, namely, C[O.sub.2.sup.+], [O.sub.2.sup.+], [H.sub.2][O.sup.+] and N[O.sup.+]; and both predict that N[O.sup.+] >> [H.sub.2][O.sup.+] (which is the second most dominant positive ion). Hence the total ion concentration is determined solely by N[O.sup.+] and the contribution of other ions to the total ion concentration is insignificant in results obtained by both methods. Table 2 also shows that the ion concentrations of C[O.sub.2.sup.+], [O.sub.2.sup.+], and [H.sub.2][O.sup.+] predicted by the CEA code are two to three orders of magnitude less than those predicted by Calcote's method. This difference is due to the method of computation of the equilibrium concentrations. As explained earlier, the Calcote method computed equilibrium concentrations by assuming that only a single ion was present in the mixture while the CEA method obtains the equilibrium concentrations using a more rigorous and accurate method, by minimizing the Gibbs free energy of the entire system under consideration. There are no detailed experiments reporting the measured concentrations of ions such as [H.sub.2][O.sup.+], [O.sub.2.sup.+], and C[O.sub.2.sup.+] in stoichiometric propane-air flames at ambient pressure and hence it is not possible to compare the numerical results with experimental data. However, for most engineering applications, the peak ion concentration is of primary importance. The more accurate calculations using the CEA code show that the main conclusion in , namely, that, in a stoichiometric propane-air mixture at T = 2200 K and P = 1 atm, N[O.sup.+] is the dominant ion with an ion concentration of ~[10.sup.7] ions/cc, is still valid.
3.3. Equilibrium Computations including the Effects of Negative Ions and [H.sub.3][O.sup.+]. Calcote and King considered only positive ions formed by thermal ionization in their equilibrium calculations. In contrast, Green and Sugden  proposed that the elementary reactions responsible for ionization in hydrocarbon flames were as follows:
CH + O = CH[O.sup.+] + [e.sup.-] (1)
CH[O.sup.+] + [H.sub.2]O = [H.sub.3][O.sup.+] + CO (2)
[H.sub.3] [O.sup.+] + [e.sup.-] = [H.sub.2]O + H (3)
Chemi-ions such as CH[O.sup.+] and [H.sub.3][O.sup.+] were not considered in Calcote's computation. The impact of the equilibrium mixture composition on the identity and magnitude of the most dominant charged species was assessed by using different sets of ions and neutrals, namely, mixtures II and III. The equilibrium composition of mixtures II and III were also computed at T = 2200 K and 1 atm for a stoichiometric propane-air mixture using the CEA code. Table 3 shows a comparison of the equilibrium concentrations of the dominant ions in mixtures I, II, and III. It is seen that N[O.sup.+] is the most dominant ion in mixtures II and III as reported by Calcote (for mixture I). Furthermore, the concentration of N[O.sup.+] for all three mixtures was between 1.7 x [10.sup.7] and 1.8 x [10.sup.7]. The ion concentration of N[O.sup.+] in mixtures II and II compares well (~5.5% higher) with the ion concentration of 1.7 x [10.sup.7] ions/cc reported by Calcote and King.
[H.sub.3][O.sup.+] was the second most dominant ion in mixtures II and III (mixture I did not include chemi-ions) and its concentration was ~ 3.06 x [10.sup.6] ions/cc, which is considerably lower than the ion concentration of N[O.sup.+]. The results obtained by using the CEA code demonstrate that neither the method of solution nor the mixture composition has a significant impact on the identity of the most dominant ion (N[O.sup.+]) and its concentration (~1.7 x [10.sup.7]) in laboratory stoichiometric propane-air flames at ambient pressure (P = 1 atm). Results presented in Sections 3.1-3.3 show that, even with the inclusion of chemi-ions (such as [H.sub.3][O.sup.+] and HC[O.sup.+]), the peak ion concentrations observed experimentally cannot be explained under conditions of thermal equilibrium between the electrons and heavy particles (ions and neutrals).
3.4. Effect of Thermal Non-Equilibrium (Elevated Electron Temperature) on the Charged Species Concentrations. The effect of elevated electron temperature on increasing the ion concentration in hydrocarbon flames is well known . The CEA code does not compute equilibrium concentrations under conditions of thermal disequilibrium hence an in-house solver based on the equilibrium constant method was used to study the effect of elevated electron temperature on the equilibrium concentration of charged species in reacting mixture II. Furthermore, the equilibrium constant method makes it easier to explain the impact of electron temperature on charged species concentrations since the equilibrium constants are functions of temperature.
Chemical equilibrium can be described by two equivalent formulations, namely, the equilibrium constant method and minimization of the free energy . The methodology for computing equilibrium concentrations using the minimization of the free energy is described in detail in . Details of the equilibrium constant method are described briefly next for the benefit of the reader (additional details are described in [17, 18]).
3.4.1. Equilibrium Constant Method. Consider a mixture of N distinct species in equilibrium in a reacting mixture under a given set of thermodynamic conditions. Let the N species be composed on M chemical elements. To obtain the equilibrium concentrations of these N species, a system of equations consisting of M linear element-conservation relations and N-M nonlinear equilibrium equations is required. For systems with charged species, an additional charge conservation equation is also needed; hence N-M-1 non-linear equilibrium equations are needed.
An elementary reaction (reversible or irreversible) involving the N chemical species can be written in the general form as
[N.summation over (k=1)] [a.sub.ks][X.sub.k] [??] [N.summation over (k=1)] [b.sub.ks][X.sub.k] (4)
where [a.sub.ks] and [b.sub.ks] are integer numbers representing stoichiometric coefficients of species k in reaction s (1, ..., N-M) and [X.sub.k] is the chemical symbol for the kth species.
The equilibrium relation for reaction s can be written as
[k.sub.s][R.sub.s] = [P.sub.s], (5)
[mathematical expression not reproducible], (6)
[k.sub.s] is the equilibrium constant of reaction s and [[c.sub.k]] is the molar concentration of species k.
The linear element-conservation equations are of the form (i = 1, M):
[N.summation over (k=1)] [d.sub.ki][X.sub.k] = [C.sub.i] (7)
where [C.sub.i] is the initial number of atoms of element i and [d.sub.ki] is the number of atoms of element i in species [X.sub.k]. For instance, in species C[O.sub.2], d = 1 for C and d = 2 for [O.sub.2].
The equilibrium constant for a given reaction s is computed using the treatment outlined in standard thermodynamic textbooks as
ln [k.sub.s] = [DELTA][S.sub.s.sup.[omicron]](T)/R - [DELTA][H.sub.s.sup.[omicron]](T)/RT (8)
where [DELTA] refers to the change in entropy and enthalpy in passing from reactants to products in a given reaction s, R is the universal gas constant, and T is the temperature. Thus, the terms in the RHS of (8) can be expressed as
[DELTA][S.sub.s.sup.[omicron]](T)/R = [N.summation over (k=1)] [v.sub.ks] [S.sub.s.sup.[omicron]](T)/RT (9)
[DELTA][H.sub.s.sup.[omicron]](T)/RT = [N.summation over (k=1)] [v.sub.ks] [H.sub.s.sup.[omicron]](T)/RT (10)
where [v.sub.ks] is the difference between the stoichiometric coefficients of products and reactants of a given species in reaction s.
[v.sub.ks] = [b.sub.ks] - [a.sub.ks] (11)
The specific enthalpy and entropy of species k required for evaluating (9) and (10) were computed by using polynomials as shown in (12) and (13) (as discussed in ). The polynomial coefficients in (12) and (13) were obtained from the CHEMKIN database.
[mathematical expression not reproducible] (12)
[mathematical expression not reproducible] (13)
In this work, the number of species, N = 25 (including ions) and number of elements, M = 4. Hence, to compute the equilibrium concentrations of the 25 species, one requires 4 element balance equations, 1 charge balance equation, and 20 nonlinear equilibrium equations. The set of species, and nonlinear equilibrium equations used in this work are as shown in Tables 6 and 7, respectively, in Appendix. The linear equations describing element and charge conservation are also shown in Appendix.
3.4.2. Validation of the In-House Solver Using the Equilibrium Constant Method. The solver used to compute the equilibrium concentrations of species in mixture II has been well-validated with results from the CEA code and experimental data. Table 4 compares the equilibrium mole fractions of important charged and neutral species in a stoichiometric propane-air mixture at T = 2200 K and P = 1 atm. It is seen that the species mole fractions predicted by the in-house solver based on the equilibrium constant method compares well with those predicted by the CEA code based on the method of minimization of the Gibbs free energy.
As further validation, the solver results were also compared with experimental data. Figure 1 shows a comparison of the computed ion current with experimental data in a constant volume cylinder consisting of a methane-air mixture for two different air-fuel ratios (AFR) .
It is seen that the location and the magnitude of the peak current is captured well. The ion current is overpredicted for t > 40 milliseconds since the thermal model describing the temperature variation in the cylinder was assumed to be adiabatic (for the sake of simplicity). Thus, the higher cylinder temperature due to the adiabatic assumption leads to a higher predicted ion current for t > 40 milliseconds, whereas in the actual experiment the heat loss from the cylinder would likely lead to a faster drop in the cylinder temperature and hence the ion current. The results presented in Table 4 and Figure 1 validate the in-house solver based on the equilibrium constant method used in this work. Reference  compares the ion current in an HCCI engine for various fuel-air mixtures with experimental data presented in . It is seen that the equilibrium chemistry computations agree well with finite-rate chemistry computations presented in reference .
3.4.3. Effect of Elevated Electron Temperature. The effect of elevated electron temperature on the equilibrium concentration of charged species in a reacting mixture was studied by computing the equilibrium constants of reactions involving electrons at the elevated electron temperature whereas the equilibrium constants of all the other reactions were computed at the gas temperature. Hence, equilibrium constants of reactions 15-18 and reaction 20 shown in Table 6 were computed at the assumed electron temperature using Eq. (8). This procedure is equivalent to the formulation described in Ref.  to compute the equilibrium concentration for a system with multiple gas temperatures.
Figure 2 shows the effect of electron temperature ([T.sub.e]) on the concentration (ions/cc) of [H.sub.3][O.sup.+] and N[O.sup.+]. One can see that as the degree of disequilibrium increases from [T.sub.e] = [T.sub.g] = 2200 K to [T.sub.e] = 5,000 K and [T.sub.g] = 2200 K, the concentration of N[O.sup.+] increases by 30%, whereas that of [H.sub.3][O.sup.+] increases by about 4 orders of magnitude. At [T.sub.e] ~ 3700K, or [T.sub.e] ~ 1.7 [T.sub.g], the concentration of [H.sub.3][O.sup.+] is approximately 1.4 x [10.sup.10] ions/cc, which is close to the experimental value reported by Goodings et al. .
Figure 2 shows that inclusion of [H.sub.3][O.sup.+] in the equilibrium mixture and elevated electron temperature can yield positive ion concentration comparable to those observed in laboratory flame experiments. At T = 2200 K, the [H.sub.3][O.sup.+] ion concentration is below that of N[O.sup.+] as stated earlier and its concentration is ~[10.sup.7] as noted by Calcote. Figure 2 also shows that, in the absence of [H.sub.3][O.sup.+], the high values of ion concentration in flames cannot be explained by N[O.sup.+] concentration alone, even at elevated electron temperatures. Figure 2 thus shows the equilibrium calculations which include chemi-ions ([H.sub.3][O.sup.+]) at slightly elevated electron temperature ([T.sub.e] ~ 1.7 [T.sup.g]) can yield ion concentrations comparable to experiments conducted on laboratory flames, thus confirming the conclusions by Calcote and King.
The sharp increase in [H.sub.3][O.sup.+] with elevated electron temperature as compared to N[O.sup.+] can be explained by examining the effect of temperature on the equilibrium constants of the reactions N + O [??] N[O.sup.+] + [e.sup.-] and [H.sub.2]O + H [??] [H.sub.3][O.sup.+] + [e.sup.-] (reactions 16 and 17 in Table 7). Figure 3 shows the variation of the equilibrium constant for reactions 16 and 17 as a function of temperature. It is seen that the equilibrium constant of reaction 17 drops by about 10 orders of magnitude whereas that of reaction 16 increases by about 4 orders of magnitude for 2000 K < [T.sub.e] < 5000 K. The stronger temperature dependence of reaction 17 yields a markedly higher concentration of [H.sub.3][O.sup.+] compared to N[O.sup.+] with increase in electron temperatures.
3.5. Effect of Fame Temperature and Pressure on the Charged Species Computations. The discussion in Section 3.4 showed that ion concentrations in laboratory flames can be explained by assuming the presence of chemi-ions ([H.sub.3][O.sup.+]) and thermal disequilibrium. The concentration of charged species, however, is also strongly dependent on the flame temperature (or gas temperature) and pressure. Calcote's equilibrium calculation were for a typical laboratory flame such as a Bunsen burner operating at 1 atmosphere with [T.sub.g] = 2200K. In typical energy conversion devices such as automotive engines and gas turbines, the combustion reactants (fuel/air mixture) are preheated and combustion takes place at higher temperatures and pressures. For instance, in an internal combustion engine, the fuel-air mixture in the engine cylinder is heated to ~ 600-650 K during the isentropic compression before the ignition begins. The unburned fuel-air mixture in the cylinder is typically between 20 and 40 atm (depending on the compression ratio and boost pressure) just prior to ignition. The burned gas temperature and pressure near top-dead-center (TDC) can be close to 2800 K and 80 atm (or higher), respectively. Under these conditions, the electron and gas temperatures can be in local thermodynamic equilibrium (LTE) or [T.sub.e] ~ [T.sub.g]. The higher collision frequency at these temperature and pressure conditions reduces the thermal disequilibrium (locally) between the electron temperature and gas temperature (see Figure 2 in ). Table 5 shows the effect of [T.sub.e], [T.sub.g], and pressure on equilibrium concentration of the dominant ions, N[O.sup.+] and [H.sub.3][O.sup.+], for mixture II. One can see a marked increase in the charged species concentration under engine-like conditions as compared to the laboratory flame under thermal equilibrium conditions ([T.sub.g] = [T.sub.e] = 2200 K and P = 1 atm). The increase in the chemi-ion concentration of [H.sub.3][O.sup.+] is pronounced even with a slight degree of thermal disequilibrium between electron and gas temperatures. The equilibrium composition of charged species under thermal equilibrium conditions are consistent with experimentally measured ion currents in a constant volume chamber and internal combustion engines as shown in [17, 18].
These results suggest that under engine-like conditions at higher temperatures and pressures, ion concentrations at the flame front can be explained by the presence of [H.sub.3][O.sup.+] under thermal equilibrium or slight thermal disequilibrium.
Figures 4 and 5 show the effect of gas temperature and pressure on N[O.sup.+] and [H.sub.3][O.sup.+], respectively, for a range of engine-like temperatures under conditions of thermal equilibrium ([T.sub.e] = [T.sub.g]). The temperature range, 1800 K < [T.sub.g] < 2800 K was chosen since it is representative of the burned gas temperature in engines during the combustion process. The pressures, P = 40 atm and P = 80 atm, are representative values at start of ignition and peak pressure in an engine. Results for P = 1 atm represent the engine pressure at the bottom dead center (BDC) and is shown as a baseline case.
Figures 4 and 5 show that gas temperature has a strong impact on the ion concentration for both N[O.sup.+] and [H.sub.3][O.sup.+] over the temperature and pressure range considered in this work. It is seen that ion concentrations increase by 4-5 orders of magnitude at T = 2800 K as compared to their values at T = 1800 K. It is also seen that a given temperature, there is about 1-2 orders of magnitude increase in ion concentration between P = 1 atm and P = 40 atm, whereas the increase in ion concentrations between P = 40 atm and P = 80 atm is not as pronounced. Figure 3 showed the marked effect of temperature on the equilibrium constants for the reactions N + O [??] N[O.sup.+] + [e.sup.-] and [H.sub.2]O + H [??] [H.sub.3][O.sup.+] + [e.sup.-] which also explains the trends shown in Figures 4 and 5.
The equilibrium computations conducted by Calcote and King for charged species concentrations in laboratory hydrocarbon flames were extended to include the effects of chemiions and thermal disequilibrium. This study also evaluated the effect of computational methodology and mixture composition on the charged species computations. The study verified that more detailed and accurate equilibrium computations with the inclusion of chemi-ions yield ion concentrations similar to those reported by Calcote when the electron temperature was assumed equal to the gas temperature for a laboratory flame with [T.sub.g] = 2200 and P = 1 atm. Equilibrium compositions of a stoichiometric propane-air mixture showed a marked increase in the concentration of [H.sub.3][O.sup.+] ions at elevated electron temperatures. However, no significant increase in the N[O.sup.+] concentration was noted with increased electron temperatures. These equilibrium computation results validate the conclusions of Calcote and King, namely, that for a laboratory flame, thermal equilibrium does not exist between species in the flame front and that ions are produced not by thermal process but by chemi-ionization. The equilibrium computations under engine-like conditions showed a marked increase in the ion concentration of [H.sub.3][O.sup.+] with and without the assumption of LTE. Equilibrium computations using a reduced set of charged and neutral species (as in mixture II) which include the dominant chemi-ion, namely, [H.sub.3][O.sup.+], can provide useful estimates of charged species concentrations in modern combustion devices.
See Tables 6 and 7.
The element balance equations of C, H, O, and N can be written as follows:
[mathematical expression not reproducible] (A.1)
Here, [N.sub.s] represents the number of moles of species [N.sub.s] and [N.sub.c], [N.sub.h], [N.sub.o], and [N.sub.n] are the total number of moles of C, H, O, and N, respectively, in the system under consideration.
The charge conservation equation can be written as
[mathematical expression not reproducible] (A.2)
Reaction 11 in Table 7 identifies the fuel as [C.sub.x][H.sub.y]; hence one can study propane by specifying x = 3, y = 8.
The data used to support the findings of this study are available from the corresponding author upon request.
The research did not receive any specific funding, but was performed as part of the employment of the author at Argonne National Laboratory.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
This material was based upon work supported by the US Department of Energy, Office of Science, under Contract No. DE-AC02-06CH11357.
 C. Cardin, B. Renou, G. Cabot, and A. M. Boukhalfa, "Experimental analysis of laser-induced spark ignition of lean turbulent premixed flames: New insight into ignition transition," Combustion and Flame, vol. 160, no. 8, pp. 1414-1427, 2013.
 T. Mathur, M. Gruber, K. Jackson et al., "Supersonic combustion experiments with a cavity-based fuel injector," Journal of Propulsion and Power, vol. 17, no. 6, pp. 1305-1312, 2001.
 Y. Xiong, D. G. Park, B. J. Lee, S. H. Chung, and M. S. Cha, "DC field response of one-dimensional flames using an ionized layer model," Combustion and Flame, vol. 163, pp. 317-325, 2016.
 H. R. N. Jones and A. N. Hayhurst, "Measurements of the concentrations of positive and negative ions along premixed fuel-rich flames of methane and oxygen," Combustion and Flame, vol. 166, pp. 86-97, 2016.
 H. F. Calcote, "Electrical properties of flames: Burner flames in transverse electric fields," in Proceedings of the 3rd Symposium on Combustion and Flame and Explosion Phenomena, pp. 245-253, USA, September 1948.
 H. F. Calcote and N. R. Pease, "Electrical properties of flames. Burner flames in longitudinal electric fields," Industrial & Engineering Chemistry, vol. 43, no. 12, pp. 2726-2731, 1951.
 J. Lawton and F. J. Weinberg, Electrical Aspects of Combustion, Clarendon Press, 1969.
 I. G. McWilliam and R. A. Dewar, "Flame ionization detector for gas chromatography ," Nature, vol. 181, no. 4611, p. 760, 1958.
 J. A. Green and T. M. Sugden, "Some observations on the mechanism of ionization in flames containing hydrocarbons," Symposium (International) on Combustion, vol. 9, no. 1, pp. 607-621, 1963.
 H. F. Calcote and I. R. King, "Studies of ionization in flames by means of langmuir probes," Symposium (International) on Combustion, vol. 5, no. 1, pp. 423-434, 1955.
 J. M. Goodings, D. K. Bohme, and T. M. Sugden, "Positive ion probe of methane-oxygen combustion," Symposium (International) on Combustion, vol. 16, no. 1, pp. 891-902, 1977
 A. B. S. Alquaity, B. Chen, J. Han et al., "New insights into methane-oxygen ion chemistry," Proceedings of the Combustion Institute, vol. 36, no. 1, pp. 1213-1221, 2017
 J. Prager, U. Riedel, and J. Warnatz, "Modeling ion chemistry and charged species diffusion in lean methane-oxygen flames," Proceedings of the Combustion Institute, vol. 31, pp. 1129-1137, 2007.
 A. M. Starik and N. S. Titova, "Kinetics of ion formation in the volumetric reaction of methane with air," Combustion, Explosion and Shock Waves, vol. 38, no. 3, pp. 253-268, 2002.
 V. Rao and D. Honnery, "A comparison of two NOx prediction schemes for use in diesel engine thermodynamic modelling," Fuel, vol. 107, pp. 662-670, 2013.
 M. Belhi, P. Domingo, and P. Vervisch, "Direct numerical simulation of the effect of an electric field on flame stability," Combustion and Flame, vol. 157, no. 12, pp. 2286-2297, 2010.
 S. M. Aithal, "Analysis of the current signature in a constant-volume combustion chamber," Combustion Science and Technology, vol. 185, no. 2, pp. 336-349, 2013.
 S. M. Aithal, "Prediction of voltage signature in a homogeneous charge compression ignition (HCCI) engine fueled with propane and acetylene," Combustion Science and Technology, vol. 185, no. 8, pp. 1184-1201, 2013.
 P. Mehresh, D. Flowers, and R. W. Dibble, "Experimental and numerical investigation of effect of fuel on ion sensor signal to determine combustion timing in homogeneous charge compression ignition engines," International Journal of Engine Research, vol. 6, no. 5, pp. 465-474, 2005.
 S. Gordon and B. J. McBride, "Computer program for calculation of complex equilibrium composition, rocket performance, incident and reflected shocks and Chapman-Jouguet detonations," NASA, vol. SP-273, p. 168, 1971.
 W. J. Miller, "Ions in flames. Evaluation and prognosis," Symposium (International) on Combustion, vol. 14, no. 1, pp. 307-320, 1973.
 R. Kee, F. Rupley, and J. Miller, "CHEMKIN-III: A FORTRAN chemical kinetics package for the analysis of gas-phase chemical and plasma kinetics," No. SAND-96-8216 SAND--96-8216, Sandia National Labs., Livermore, CA (United States), 1996.
 J. C. Keck, "Rate-controlled constrained-equilibrium theory of chemical reactions in complex systems," Progress in Energy and Combustion Science, vol. 16, no. 2, pp. 125-154, 1990.
 D. W. Cho, S. H. Lee, and S. J. Na, "Characterization of welding arc and weld pool formation in vacuum gas hollow tungsten arc welding," Journal of Materials Processing Technology, vol. 213, no. 2, pp. 143-152, 2013.
S. M. Aithal (iD)
Argonne National Laboratory, 9700 S. Cass Ave, Lemont, IL 60439, USA
Correspondence should be addressed to S. M. Aithal; email@example.com
Received 30 May 2018; Accepted 16 September 2018; Published 4 October 2018
Academic Editor: Kalyan Annamalai
Caption: Figure 1: Comparison of the computed ion current with experimental data in a constant volume cylinder (C[H.sub.4]-air mixture, air-fuel-ratio (AFR) =15 & 17) .
Caption: Figure 2: Effect of electron temperature ([T.sub.e]) on the concentration (ions/cc) of [H.sub.3][O.sup.+] and N[O.sup.+] under laboratory flame conditions (P = 1 atm) for a stoichiometric propane-air mixture.
Caption: Figure 3: Effect of electron temperature ([T.sub.e]) on the equilibrium constants of [H.sub.3][O.sup.+] and N[O.sup.+].
Caption: Figure 4: Effect of gas temperature on N[O.sup.+] ion concentration for P = 1, 40 and 80 atm under conditions of thermal equilibrium ([T.sub.e] = [T.sub.g]) in a stoichiometric propane-air mixture.
Caption: Figure 5: Effect of gas temperature on [H.sub.3][O.sup.+] ion concentration for P = 1, 40 and 80 atm under conditions of thermal equilibrium ([T.sub.e] = [T.sub.g]) in a stoichiometric propane-air mixture.
Table 1: Comparison of equilibrium concentrations (Calcote and CEA). Species Calcote (#/cc) CEA (#/cc) [H.sub.2] 8.2 x [10.sup.15] 8.4 x [10.sup.15] [N.sub.2] 2.4 x [10.sup.18] 2.4 x [10.sup.18] CO 3.3 x [10.sup.16] 3.1 x [10.sup.16] C[O.sub.2] 3.5 x [10.sup.17] 3.5 x [10.sup.17] H 9.2 x [10.sup.15] 9.4 x [10.sup.15] OH 7.9 x [10.sup.15] 8.6 x [10.sup.15] O 6.6 x [10.sup.14] 6.0 x [10.sup.14] [H.sub.2]O 5.0 x [10.sup.17] 5.0 x [10.sup.17] [O.sub.2] 1.5 x [10.sup.15] 1.5 x [10.sup.15] NO 6.3 x [10.sup.15] 5.9 x [10.sup.15] N[O.sup.+] 1.7 x [10.sup.7] 1.8 x [10.sup.7] Table 2: Comparison of the dominant equilibrium ion concentrations using the Calcote method and minimization of the Gibbs free energy method. Species Calcote CEA (all species) (ions/cc) (ions/cc) N[O.sup.+] 1.7 x [10.sup.7] 1.8 x [10.sup.7] [H.sub.2][O.sup.+] 4.4 x [10.sup.4] 4.02 x [10.sup.2] [O.sub.2.sup.+] 2.2 x [10.sup.4] 8.4 x [10.sup.1] C[O.sub.2.sup.+] 1.5 x [10.sup.3] 2 Table 3: Comparison of equilibrium concentrations of the dominant ions in mixtures I-III. Species Mixture I Mixture II Mixture III N[O.sup.+] 1.7 x [10.sup.7] 1.8 x [10.sup.7] 1.78 x [10.sup.7] [H.sub.3] -- 3.06 x [10.sup.6] 3.05 x [10.sup.6] [O.sup.+] HC[O.sup.+] -- 1.03 x [10.sup.3] 1.03 x [10.sup.3] Table 4: Validation of equilibrium constant method with CEA code. No. Species Minimization of Gibbs Equilibrium Constant Free Energy (CEA code) method Mole fractions Mole fractions 1 [O.sub.2] 4.46 x [10.sup.-3] 4.45 x [10.sup.-3] 2 C[O.sub.2] 1.05 x [10.sup.-1] 1.06 x [10.sup.-1] 3 NO 1.76 x [10.sup.-3] 1.85 x [10.sup.-3] 4 [H.sub.2] 2.54 x [10.sup.-3] 2.55 x [10.sup.-3] 5 N[O.sup.+] 5.35 x [10.sup.-12] 5.06 x [10.sup.-12] 6 O[H.sup.-] 3.91 x [10.sup.-13] 3.83 x [10.sup.-13] 7 [O.sup.-] 1.05 x [10.sup.-14] 1.14 x [10.sup.-14] Table 5: Effect of [T.sub.e], [T.sub.g] and pressure on equilibrium concentration of ions (mixture II) in a stoichiometric propane-air mixture. [T.sub.g] [T.sub.e] (K) Pressure [H.sub.3][O.sup.+] (K) (atm) (ions/cc) 1 2200 2200 1 5.66 x [10.sup.6] 2 2800 2800 40 1.53 x [10.sup.10] 3 2800 4200 ([T.sub.e] = 40 2.81 x [10.sup.12] 1.5[T.sub.g]) N[O.sup.+] Total Application (ions/cc) (ions/cc) 1 1.69 x [10.sup.7] 2.2 x [10.sup.7] Laboratory flame (Bunsen burner) with thermal equilibrium 2 5.16 x [10.sup.10] 6.69 x [10.sup.10] Burned gas conditions (engines) with thermal equilibrium 3 7.05 x [10.sup.10] 2.81x1012 Burned gas conditions (engines) with thermal dis-equilibrium Table 6: List of species. Species 1 [C.sub.3][H.sub.8] 2 [O.sub.2] 3 C[O.sub.2] 4 [H.sub.2]O 5 [N.sub.2] 6 N 7 O 8 NO 9 OH 10 H 11 [N.sub.2]O 12 CO 13 [H.sub.2]O 14 N[O.sub.2] 15 H[O.sub.2] 16 [C.sub.2[H.sub.2] 17 C 18 CH 19 HC[O.sup.+] 20 e- 21 [H.sub.3][O.sup.+] 22 N[O.sup.+]] 23 [O.sub.2.sup.-] 24 [O.sup.-] 25 O[H.sup.-] Table 7: List of reactions. 1 1/2 [H.sub.2] = H [k.sub.1] = [H]/[[H.sub.2]].sup.0.5] [square root of P] 2 1/2 [O.sub.2] = O [k.sub.2] = [O]/[[O.sub.2]].sup.0.5] [square root of P] 3 1/2 [N.sub.2] = N [k.sub.3] = [N]/[[N.sub.2]].sup.0.5] [square root of P] 4 1/2 [H.sub.2] + 1/2 [k.sub.4] = [O.sub.2] = OH [OH]/[[H.sub.2]].sup.0.5] [[O].sup.0.5] 5 1/2 [N.sub.2] + 1/2 [k.sub.5] = [O.sub.2] = NO [NO]/[[N.sub.2]].sup.0.5] [[O].sup.0.5] 6 [H.sub.2] + 1/2 [O.sub.2] [k.sub.6] = = [H.sub.2]O [[H.sub.2]O]/[[H.sub.2]] [[[O.sub.2]].sup.0.5] [(P).sup.-0.5] 7 CO + 1/2 [O.sub.2] = [k.sub.7] = [C[O.sub.2]]/[CO] C[O.sub.2] [[[O.sub.2]].sup.0.5] [(P).sup.-0.5] 8 NO + 1/2 [O.sub.2] = [k.sub.8] = [N[O.sub.2]]/[NO] N[O.sub.2] [[[O.sub.2]].sup.0.5] [(P).sup.-0.5] 9 [O.sub.2] + 1/2 [H.sub.2] [k.sub.9] = = H[O.sub.2] [H[O.sub.2]]/[[O.sub.2]] [[[H.sub.2]].sup.0.5] [(P).sup.-0.5] 10 [N.sub.2] + 1/2 [O.sub.2] [k.sub.10] = = [N.sub.2]O [[N.sub.2]O]/[[N.sub.2]] [[[O.sub.2]].sup.0.5] [(P).sup.(y/2-1)] 11 [C.sub.x][H.sub.y] + (x + [mathematical expression not 2) [O.sub.2] = reproducible] xC[O.sub.2] + y/2 [H.sub.2]O 12 [C.sub.x] [H.sub.y] = [k.sub.12] = [C.sub.x] [H.sub.y-1] + H [[C.sub.x][H.sub.y-1]] [H]/[[C.sub.x][H.sub.y]] P 13 [C.sub.2] [H.sub.2] = 2C [k.sub.13] = [[C].sup.2] + [H.sub.2] [[H.sub.2]]/[[C.sub.2][H.sub.2]] P 14 CO + OH = CH + [O.sub.2] [k.sub.14] = [CH] [[O.sub.2]]/[CO][HO] 15 CH + O = HC[O.sup.+] + [k.sub.15] = [HC[O.sup.+]] [e.sup.-] [[e.sup.-]]/[CH][O] 16 N + O = N[O.sup.+] + [k.sub.16] = [N[O.sup.+]] [e.sup.-] [[e.sup.-]]/[N][O] 17 [H.sub.3] [O.sup.+] + [k.sub.17] = [[H.sub.2]O] [e.sup.-] = [H.sub.2]O + [H]/ H [[H.sub.3][O.sup.+]][[e.sup.-]] 18 [O.sub.2] + [e.sup.-] = [k.sub.18] = [[O.sub.2.sup.-]]/ [O.sub.2.sup.-] [[O.sub.2]][[e.sup.-]] [p.sup.-1] 19 O + [O.sub.2.sup.-] = [k.sub.19] = [O.sub.2] + [O.sub.-] [[O.sub.2]][[O.sup.-]]/ [O][[O.sub.2.sup.-]] 20 OH + [e.sup.-] = [k.sub.20] = O[H.sup.-] [O[H.sup.-]]/[OH][[e.sup.-]] [p.sup.-1]
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Publication:||Journal of Combustion|
|Date:||Jan 1, 2018|
|Previous Article:||Prediction of Pollutant Emissions from Bluff-Body Stabilised Nonpremixed Flames.|
|Next Article:||Coal and Biomass Combustion.|