Next Article in Journal
Stepping Out of Equilibrium: The Quest for Understanding the Role of Non-Equilibrium (Thermo-)Dynamics in Electronic and Electrochemical Processes
Previous Article in Journal
Designing Two Secure Keyed Hash Functions Based on Sponge Construction and the Chaotic Neural Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Implementation of the Onsager Theorem to Evaluate the Speed of the Deflagration Wave

Faculty of Aerospace Engineering, Technion-Israel Institute of Technology, Haifa 3200003, Israel
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(9), 1011; https://doi.org/10.3390/e22091011
Submission received: 29 July 2020 / Revised: 29 August 2020 / Accepted: 8 September 2020 / Published: 10 September 2020

Abstract

:
While considering the deflagration regime, the thermal theory of combustion proposes that the mechanism of heat transfer from the flame exothermic zone to the front neighborhood reactants layer dominates the flame behavior. The introduction of the Fourier law allows a closed solution of the continuity and energy conservation equations to yield the burning velocity. It is, however, clear that this classical solution does not conform to the momentum equation. In the present work, instead of introducing the Fourier law, we suggest the introduction of a simplified version of the Onsager relationship, which accounts for the entropy increase due to the heat transfer process from the front layer to its successive layer. Solving for the burning velocity yields a closed solution that also definitely conforms to the momentum equation. While it is realized that the pressure difference across the flame front in the deflagration regime is very small, we believe that violating the momentum equation is intolerable. Quite a good fitting, similarly to the classic theory predictions, has been obtained between our predictions and some experimentally observed values for the propagation flame deflagration velocity, while here, the momentum equation is strictly conserved.

1. Background

Flame propagation in a flammable quiescent homogeneous medium under steady-state, steady-flow (SSSF) conditions may occur under two different discrete regimes: deflagration (sub-sonic), or detonation (super-sonic). In 1886 Mallard and Le-Chateliet were the first to observe how a flammable gas/air mixture that is ignited at the open end of a long tube propagates and develops along the tube. Since then, a large body of experimental observations has been collected for a number of fuels under a range of initial conditions. The two different regimes were clearly noted and a number of successful attempts were made to explain them. Some useful expressions have also been developed to enable reasonable predictions.
When the three conservation equations (continuity, momentum and energy) are applied to a one-dimensional SSSF flame propagation process in a tube, it may soon be realized that in addition to a suitable equation of state, one more equation is needed to solve them [1]. The five unknowns are the inlet velocity, exit velocity, exit pressure, exit temperature and exit density. Of particular interest to us is the inlet velocity since it is defined as the basic flame propagation velocity. Rankine, Hugoniot and Rayleigh [1] showed that this set of equations (the three conservation equations combined with the equation of state) yields two simple algebraic equations (one quadratic and one linear). Customarily, the two equations are plotted on a pressure-density−1 coordinate system. The intersections of these two curves define the different regimes. It can be shown that the slope of the Rankine–Rayleigh equation (linear equation) depends on the flame propagation velocity, which is an unknown parameter. It was suggested by Chapman and Jouguet [1] that the two tangent points between these two curves represent the deflagration and the detonation conditions while the two different slopes of the Rankine–Rayleigh equation represent the two propagation velocities. When comparing the results with available experimental observations, it is clear that the predicted detonation propagation velocity fits well. The predicted propagation velocity in the deflagration regime is however far higher, about one or two orders of magnitude higher than the observed value. The thermal theory of combustion date back to Mallard and Le Chatelier [1], who proposed that unlike in the detonation region, where the hydrodynamic effects govern the flame propagation, when deflagration is considered, the mechanism of heat transfer from the flame exothermic zone to the front neighborhood reactants layer dominates the flame behavior. Zeldovich, Frank-Kamenetskii and Semenov [2] further improved the concept and presented it in detail. When the Fourier law is introduced to the above set of equations, the five equations allow a closed solution [2]. Some plausible assumptions may result in simple algebraic expression for the propagation velocity. It, however, can be shown that this classical solution does not conform to the momentum equation [[1], p. 150 line 42]. It is true that the pressure difference across the flame front in the deflagration regime is very small, but still, violating the momentum equation can hardly be accepted.
In the present work, we suggest a new approach to bridge this classic gap. Instead of introducing the Fourier law, we suggest introducing the Onsager relationship [3], which accounts for the entropy increase due to the heat transfer process from the front layer to its successive layer. Based on a few logical assumptions, a simplified Onsager relationship has been introduced to the set of the conservation equations, and a closed solution that definitely conforms to the momentum equation has been obtained. A very good fitting has been obtained between the predicted and the experimentally observed values of the propagation flame deflagration velocity.

2. Theory

Consider a combustion wave propagating in SSSF conditions through a combustible mixture in a two open-sided, open-ended adiabatic stationary tube. The wave is opposed by the reactants (unburned gases) flowing at a velocity exactly equal to the wave propagation velocity, so that the wave is fixed with respect to the tube. Following the classical theory, while neglecting any transport phenomena in the fluid (no viscosity, heat transfer and matter diffusivity), we may apply the following conservation equations between the tube entrance (reactants) and the tube exit (products):
Continuity:
ρ R u R = ρ P u P
Momentum:
ρ R u R 2 + P R = ρ P u P 2 + P P
Energy:
u R 2 2 + h R = u P 2 2 + h P
where ρ , P , and u are the fluid density, pressure and velocity, respectively, and h its enthalpy. The indices R and P indicate reactant and product, respectively. Further we assume that both the reactants and products follow the ideal gas law, thus:
P = ρ R T
We have now four equations with five unknowns: the temperature, pressure, density and velocity of the products, and the reactants’ velocity. The classical theory suggests that if we further assume invariant gas constant across the flame and identical heat capacity for the reactants and products, Equations (1)–(5) may be easily combined to yield the following two relationships [1]:
The Rayleigh equation:
P P = P R ( 1 ρ P 1 ρ R ) ρ R 2 u R 2
The Hugoniot equation:
k k 1 ( P P ρ P P R ρ R ) 1 2 ( P P P R ) ( 1 ρ R + 1 ρ P ) = h c
where k is the specific heat ratio and h c the heat of combustion. Now we have two equations with three unknowns. The Rayleigh and Hugoniot equations may be plotted on a P 1 / ρ coordinate system, while the intersections of the Rayleigh lines with the Hugoniot curve (four intersection points) define four possible solutions. Note that the two different families of the Rayleigh lines slopes (two families of beams) that intersect with the Hugoniot curve signify the reactants’ velocities in the detonation and in the deflagration regimes, which are still unknown parameters. When considering the slopes of these two families, it is clear that there are only two discrete slopes that generate two tangent beams to the Hugoniot curve, the steeper beam designates the detonation wave and the milder beam designates the deflagration wave. Stability and entropy considerations reveal that the only stable detonation wave is represented by the tangent point between the Rayleigh and the Hugoniot curves (the upper Chapman–Jouguet point) [1]. The suitable algebraic equation for the tangent conditions provides the additional equation that is needed to solve Equations (1)–(4) (or (5) and (6)). A wide range of experimental observations confirmed the velocity of the detonation wave as predicted by this set of equations. The prediction of the deflagration velocity has however been found far higher (1–2 order of magnitudes) than the observed value.
It was proposed by Mallard and Le Chatelier that the structure of the deflagration wave (which is essentially different from the structure of the detonation wave) suggests that the flame propagation under the deflagration regime is controlled mainly by the heat transfer from the flame front exothermic layer to the fresh reactants’ layer ahead of it [1]. Semenov, Zeldovich, and Frank-Kamenetskii adopted this concept and considered the Fourier law of heat conduction to take account of it [2]:
J t = K t d T d x
where J t is the heat flux and K t is the thermal conductivity of the reactants.
While assuming a linear temperature profile and ignoring the mass diffusion, the classical thermal theory [1] suggests equating the heat conduction to the energy increase from reactants to products, thus obtaining an algebraic equation that provides the additional equation that is needed to solve Equations (1)–(4) (or (5) and (6)) for a deflagration wave. Further simplifications and the introduction of the flame thickness concept [1] yield the classical relationship between the laminar burning velocity, S L (=the reactants’ velocity), the flame thickness, δ , and the thermal diffusivity of the reactants, α :
S L = α δ
Experimental observations confirmed the velocity of the deflagration wave as predicted by Equation (8). It can be shown however that this classic simple solution does not conform to the momentum equation [[1], p. 150 line 42]. Although the pressure difference across the flame front in the deflagration regime is very small, the violation of the momentum equation is improper.

3. Present Suggested Approach

As an alternative to the Semenov group’s methodology, we propose in the present work a different approach to bridge this gap. Instead of introducing the Fourier law, we suggest introducing the Onsager relationship, which accounts for the entropy increase due to the heat transfer process from the front layer to its successive layer (see [3,4] for the essentials of the Onsager relationship). Based on a few logical assumptions, a simplified Onsager relationship can be derived and introduced to the set of conservation equations. We start from a basic postulation that suggests, in general, when a system is close to the state of equilibrium [4], any flux J i (heat, molecules, momentum, etc.) is expected to be proportional to the sum of linear functions of any driving force X j (gradient of temperature−1, concentration, velocity, etc.). Thus:
J i = j = 1 n L i j X j
where the coefficients L i j are the Onsager phenomenological coefficients. Based on statistical thermodynamics, Onsager [4] showed that whenever the total entropy production per unit time per unit volume resulting from all the irreversible processes can be represented as a linear sum of products of forces and fluxes, then L i j = L j i . Thus:
d S d t = i J i X i > 0
We assume that in the case of combustion, the contribution of the molecular diffusivity to the entropy production is marginal, and thus the entropy is generated due to two pronounced processes: heat generation and heat conduction, thus:
d S d t = ( d S d t ) h e a t   g e n e r a t i o n + ( d S d t ) h e a t   c o n d u c t i o n
Assuming that the main heat production occurs at the flame front at T P , and thus the entropy production due to the heat generation (here the driving force is the gradient of 1 / T , while the flux is the heat), is:
( d S d t ) h e a t   g e n e r a t i o n = Q ˙ ε δ T P = m ˙ h c ε δ T P = ρ R S L h c ε δ T P
where m ˙ and ε δ are the mass flow rate per unit area and the reactions’ zone thickness ( ε is an arbitrary number smaller but close to 1), respectively. Q ˙ is the heat production per unit area. Assuming d T d x = c o n s t , and K = c o n s t (here K is the thermal conductivity of the gases in [ k J m · s · K ] units), yields:
( d S d t ) h e a t   c o n d u c t i o n = K ( T P T R ) 2 T P T R δ 2
Assuming an ideal gas with c p = c o n s t , and that the entropy production equals to the entropy increase from reactants to products results:
K ( T P T R ) 2 T P T R δ 2 + ρ R S L h c ε δ T P = ρ R S L δ [ c p l n ( T P T R ) R l n ( P P P R ) ]
Note: the right term in Equation (14) represents the entropy increase from reactants to products. This term is obtained while integrating the thermodynamic relationship T d s = d h v d p from the inlet to the exit, while assuming an ideal gas law with c p = c o n s t .
Solving for S L :
S L = K ρ R c p δ ( T P T R ) 2 T P T R [ l n ( T P T R ) k 1 k l n ( P P P R ) h c ε c p T P ]
Equation (15) correlates the flame laminar velocity with three unknowns: T P , P P and the flame thickness δ. Following the classic thermal theory, in which the heat generation is compared to the heat convection, the flame thickness can be easily expressed in terms of the rate of heat generation, which include an Arrhenius-type expression [1]. It is noted that the result presented in Equation (15) resembles the classical result obtained from the energy balance between heat conduction and heat convection, for which the right term = 1. Equations (1)–(4) and (15) thus form a complete set of five equations with five unknowns: the temperature, pressure, density and velocity of the products, and the reactants’ velocity, and thus conforming to the three conservation laws, which include the momentum equation.

4. Discussion

Some important results of the present calculations are depicted in Figure 1. This schematic chart shows the range of the two regimes (classic map), together with the variation of the products’ entropy for reactants that initially are at standard conditions. As expected, the stable deflagration wave is located on the Hugoniot curve (black curve) on the deflagration regime side ( 1 / ρ > 1 / ρ R ) , where the products’ density is lower than that of the Chapman–Jouguet (C-J) lower point, while the pressure is lower than that of the reactants. The increase in the entropy at this point reflects the entropy gain due to the heat generation and the heat conduction as suggested in Equation (11).
For typical values of a stoichiometric hydrocarbon/air mixture, the right term of Equation (15), is in the order of one. Moreover, since ε is an arbitrary number that is smaller but close to one, its effect on the burning velocity (Equation (15)) is practically insignificant. For the following brief comparison, we further assume that ε 1 . The effect of the reactants’ temperature, pressure and equivalence ratio on the burning velocity are therefore mainly attributed to their effects on the flame thickness and to a smaller extent on the reactants’ density. This is in a similar manner as suggested by the classical result, which is obtained from the energy balance between heat conduction and heat convection while ignoring the pressure difference across the flame. The ratio between the predicted burning wave velocity, S L , of the present approach and that of the classical thermal theory is:
( S L ) P r e s e n t   a p p r o a c h ( S L ) T h e r m a l   t h e o r y = ( T P T R ) 2 T P T R [ l n ( T P T R ) k 1 k l n ( P P P R ) h c ε c p T P ]
While the present approach ensures that the momentum basic law is conserved, Equation (16) suggests that the end results for the burning velocity do not noticeably differ from the value as predicted by the classical theory. For a typical hydrocarbon flame, the ratio is close to one.
Figure 2 shows some experimental results and the present predictions (solid line) for methane/air mixtures. As expected, the predicted results show quite a good fitting to the observed values of several experimental data, as collected by Andrews and Bradley [5] in a wide range of initial conditions. It is however noted that since the Arrhenius-type equation for a single-step reaction does not maximize at or close to the stoichiometric conditions, the predictions for fuel rich mixtures ( φ > 1 in Figure 2c) are inadequate.

5. In Summary

The thermal theory of combustion proposes that for deflagrating flames the mechanism of heat transfer from the flame exothermic zone to the front neighborhood reactants layer dominates the flame behavior. The introduction of the Fourier law provides the additional required equation that after some plausible assumptions allows a closed solution to yield the burning velocity. It is however clear that this classical solution does not conform to the momentum equation. In the present work, instead of introducing the Fourier law, we suggest introducing a simplified version of the Onsager relationship, which accounts for the entropy increase due to the heat transfer process from the front layer to its successive layer. In this way, a closed solution that definitely conforms to the momentum equation has been obtained. It is anticipated that the pressure difference across the flame front in the deflagration regime is very small, but still, violating the momentum equation can hardly be accepted. Unsurprisingly, quite a good fitting has been obtained between the predicted and the experimentally observed values of the propagation flame deflagration velocity, while the momentum equation is strictly conserved.

Author Contributions

Data curation, Y.M. and I.M.-M.; methodology, E.S. and B.C. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the financial support of Minerva Research Center (Max Planck Society Contract No. AZ5746940764).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

c p [ k J k g · K ] Heat capacity at constant pressure
h [ k J k g ] Enthalpy
J [ k J m 2 s ] Flux
K t [ k J m · s · K ] Thermal conductivity
k [ ] Specific heat ratio
L i j [ k J · K m · s ] The Onsager phenomenological coefficients
m ˙ [ k g m 2 s ] Mass flow rate per unit area
Q ˙ [ k J m 2 s ] Heat production per unit area
P [ k P a ] Pressure
R [ k J k g · K ] Gas constant
S L [ m s ] Laminar burning velocity
S [ k J m 3 · K ] Entropy per unit volume
T [ K ] Temperature
t [ s ] Time
u [ m s ] Velocity
X j [ 1 m · K ] Driving Force
x [ m ] Coordinate
Greek symbols
α [ m 2 s ] Thermal diffusivity
δ [ m ] Flame thickness
ε [ ] An arbitrary number smaller but close to 1
ρ [ k g m 3 ] Density
Indices
c Combustion
P Products
R Reactants
t Thermal

References

  1. Glassman, I.; Yetter, R.A. Combustion, 4th ed.; Elsevier: Amsterdam, The Netherlands, 2008. [Google Scholar]
  2. Semenov, N.N. Some Problems in Chemical Kinetics and Reactivity; Pergamon Press: Oxford, UK, 1959. [Google Scholar]
  3. Callen, H.B. Thermodynamics and an Introduction to Thermostatistics, 2nd ed.; John Wiley and Sons: Hoboken, NJ, USA, 1985. [Google Scholar]
  4. Onsager, L. Reciprocal Relations in Irreversible Processes. I. Phys. Rev. 1931, 37, 405. [Google Scholar] [CrossRef]
  5. Andrews, G.E.; Bradley, D. The burning velocity of methane-air mixture. Combust. Flame 1972, 19, 275–288. [Google Scholar] [CrossRef]
  6. Johnston, W.; Johnston, W.C. Measures flame velocity of fuels at low pressures. Soc. Automot. Eng. J. 1947, 55, 62–65. [Google Scholar]
  7. Dugger, G.L. Effect of Initial Mixture Temperature on Flame Speed of Methane-Air, Propane-Air, and Ethylene-Air Mixtures (NACA Report of Investigations No. 1061); Lewis Flight Propulsion Laboratory: Cleveland, OH, USA, 1952. [Google Scholar]
  8. Halpern, C. Measurement of flame speeds by a nozzle burner method. J. Res. Natl. Bur. Stand. 1958, 60, 535. [Google Scholar] [CrossRef]
  9. Babkin, V.; Kozachenko, L. Study of normal burning velocity in methane-air mixtures at high pressures. Combust. Explos. Shock Waves 1966, 2, 46–52. [Google Scholar] [CrossRef]
  10. Barassin, A.; Lisbet, R.; Combourieu, J.; Laffitte, P. Bull. Soc. Chim. France 1967, 7, 2521. [Google Scholar]
  11. Gu, X.J.; Haq, M.Z.; Lawes, M.; Woolley, R. Laminar burning velocity and Markstein lengths of methane–air mixtures. Combust. Flame 2000, 121, 41–58. [Google Scholar] [CrossRef]
  12. Egerton, A.C.; Lefebvre, A. Flame propagation: The effect of pressure variation on burning velocities. Proc. R. Soc. Lond. Ser. Math. Phys. Sci. 1954, 222, 206–223. [Google Scholar]
  13. Clingman, W.H.; Pease, R.N. Critical Considerations in the Measurement of Burning Velocities of Bunsen Burner Flames and Interpretation of the Pressure Effect. Measurements and Calculations for Methane 1. J. Am. Chem. Soc. 1956, 78, 1775–1780. [Google Scholar] [CrossRef]
  14. Diederichsen, J.; Wolfhard, H. The burning velocity of methane flames at high pressure. Trans. Faraday Soc. 1956, 52, 1102–1109. [Google Scholar] [CrossRef]
  15. Manton, J.; Milliken, B.; Proc, B. Gas Dynamics Symposium on Aerothermochemistry; Northwestern University: Evanton, Ill, USA, 1955. [Google Scholar]
  16. Singer, J.; Grumer, J.; Cook, E. Burning Velocities by the Bunsen Burner Method. I. Hydrocarbon-Oxyen Mixtures at One Atmosphere. II. Hydrocarbon-Air Mixtures at Subatmospherc Presusres; James Forrestal Research Center, Princeton University: Princeton, NJ, USA, 1955. [Google Scholar]
  17. Strauss, W.A.; Edse, R. Gas Dynamics Symposium on Aerothermochemistry; Northwestern University: Evanston, IL, USA, 1956. [Google Scholar]
  18. Gilbert, M. The influnence of pressure on flame speed. Symp. (Int.) Combust. 1957, 6, 74–83. [Google Scholar] [CrossRef]
  19. Agnew, J.; Graiff, L. The pressure dependence of laminar burning velocity by the spherical bomb method. Combust. Flame 1961, 5, 209–219. [Google Scholar] [CrossRef]
  20. Babkin, V.; Kozachenko, L.; Kuznetsov, I. The effect of pressure on the normal burning velocity of a methane–air mixture. Zhurnal Prikl. Mekhaniki Tekhnicheskoi Fiz. 1964, 3, 145–149. [Google Scholar]
  21. Bradley, D.; Hundy, G. Burning velocities of methane-air mixtures using hot-wire anemometers in closed-vessel explosions. Symp. Int. Combust. 1971, 13, 575–583. [Google Scholar] [CrossRef]
  22. Van Maaren, A.; Thung, D.; de Goey, L.P.H. Measurement of flame temperature and adiabatic burning velocity of methane/air mixtures. Combust. Sci. Technol. 1994, 96, 327–344. [Google Scholar] [CrossRef] [Green Version]
  23. Vagelopoulos, C.M.; Egolfopoulos, F.N. Direct experimental determination of laminar flame speeds. Symp. Int. Combust. 1998, 27, 513–519. [Google Scholar] [CrossRef]
  24. Bosschaart, K.J.; de Goey, L.R.H. The laminar burning velocity of flames propagating in mixtures of hydrocarbons and air measured with the heat flux method. Combust. Flame 2004, 136, 261–269. [Google Scholar] [CrossRef]
Figure 1. Schematic chart of the range of the two regimes (classic map), together with the variation of the products’ entropy.
Figure 1. Schematic chart of the range of the two regimes (classic map), together with the variation of the products’ entropy.
Entropy 22 01011 g001
Figure 2. Some experimental results [5] and the present prediction (solid line) for methane/air mixtures. Here the equivalence ratio is defined as the ratio between the air-fuel mass ratio and the stoichiometric fuel-air mass ratio. (a) Flame velocity vs. reactants’ temperature [6,7,8,9,10,11]. (b) Flame velocity vs. reactants’ pressure [9,12,13,14,15,16,17,18,19,20,21]. (c) Flame velocity vs. equivalence ratio [11,22,23,24].
Figure 2. Some experimental results [5] and the present prediction (solid line) for methane/air mixtures. Here the equivalence ratio is defined as the ratio between the air-fuel mass ratio and the stoichiometric fuel-air mass ratio. (a) Flame velocity vs. reactants’ temperature [6,7,8,9,10,11]. (b) Flame velocity vs. reactants’ pressure [9,12,13,14,15,16,17,18,19,20,21]. (c) Flame velocity vs. equivalence ratio [11,22,23,24].
Entropy 22 01011 g002

Share and Cite

MDPI and ACS Style

Sher, E.; Moshkovich-Makarenko, I.; Moshkovich, Y.; Cukurel, B. Implementation of the Onsager Theorem to Evaluate the Speed of the Deflagration Wave. Entropy 2020, 22, 1011. https://doi.org/10.3390/e22091011

AMA Style

Sher E, Moshkovich-Makarenko I, Moshkovich Y, Cukurel B. Implementation of the Onsager Theorem to Evaluate the Speed of the Deflagration Wave. Entropy. 2020; 22(9):1011. https://doi.org/10.3390/e22091011

Chicago/Turabian Style

Sher, Eran, Irena Moshkovich-Makarenko, Yahav Moshkovich, and Beni Cukurel. 2020. "Implementation of the Onsager Theorem to Evaluate the Speed of the Deflagration Wave" Entropy 22, no. 9: 1011. https://doi.org/10.3390/e22091011

APA Style

Sher, E., Moshkovich-Makarenko, I., Moshkovich, Y., & Cukurel, B. (2020). Implementation of the Onsager Theorem to Evaluate the Speed of the Deflagration Wave. Entropy, 22(9), 1011. https://doi.org/10.3390/e22091011

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop