Next Article in Journal
Numerical Investigation of Thermally Developing and Fully Developed Electro-Osmotic Flow in Channels with Rounded Corners
Next Article in Special Issue
Computational Study of Premixed Flame Propagation in Micro-Channels with Nonslip Walls: Effect of Wall Temperature
Previous Article in Journal
Low-Speed DSMC Simulations of Hotwire Anemometers at High-Altitude Conditions
Previous Article in Special Issue
Stability Analysis for an Interface with a Continuous Internal Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ultra-Lean Gaseous Flames in Terrestrial Gravity Conditions

Joint Institute for High Temperatures of Russian Academy of Sciences, 125412 Moscow, Russia
*
Author to whom correspondence should be addressed.
Submission received: 15 December 2020 / Revised: 28 December 2020 / Accepted: 30 December 2020 / Published: 3 January 2021

Abstract

:
Development of the combustion process in the gaseous mixtures of near-limit composition is of great interest for fundamental aspects of combustion theory and fire-safety applications. The dynamics of ultra-lean gaseous flames in near-limit mixtures is governed by many effects, such as buoyancy, preferential diffusion, radiation, and instability development. Though ultra-lean combustion was extensively studied in microgravity conditions, the influence of gravity on the ultra-lean flame structure and stability is still poorly understood. The paper is devoted to deepening the knowledge of ultra-lean flame dynamics in hydrogen-air mixtures under terrestrial gravity conditions. The spatial structures of the flame developing under the effect of buoyancy forces are investigated employing detailed numerical analysis. Different modes of near-limit flame evolution are observed depending on the mixture concentration. In particular, we registered and described three distinct spatial structures: individual kernels tending to extinguish in leanest compounds, complex multi-kernel structures in marginal compositions, and stable cap-shaped flames in more chemically active mixtures. We apply the flame-bubble analogy to interpret flame dynamics. On this basis, the diagram in the Re-Fr plane is developed. That allows classifying the emerging flame structures and determine flame stability. Additionally, different ignition modes are studied, and the mechanisms determining the impact of ignition mode on the flammability limits are distinguished. Obtained results provide useful insights into the processes of flame quenching and development in near-limit hydrogen-air mixtures under real gravity conditions and can be applied in the design of contemporary fire-safety systems.

1. Introduction

Even today, fire and explosion safety remains a relevant problem for energy, technology, and industrial applications. Among the most acute issues is the determination of conditions and mechanisms of non-stationary flame propagation, flame acceleration, and transition to fast subsonic and supersonic flame propagation regimes, including detonations. Another critical problem is related to the assessment of flammability limits and, in particular, to the investigation of combustion modes in lean and ultra-lean flammable gaseous mixtures. It is widely accepted that the most dangerous and devastating accidents are related to the occurrence of fast flames and detonations, so the problem of flame non-stationary development was a focus of intense studies worldwide [1]. In past decades, through real experimental studies and numerical analysis, many significant results have been obtained crucial for understanding non-stationary flame acceleration in confined vessels and combustion chambers filled with reactive gaseous mixtures. However, fast flame development, first of all, is observed in highly active mixtures with the composition near stoichiometric one. In the conditions intrinsic to the real technical systems and energy facilities operation, the formation and explosion of stoichiometric gaseous compounds are uncommon. When considering the scenarios related to the explosion safety, in most cases, combustible gases are gradually released as a result of decompression of fuel tanks or chemical generation of the reactive component. In large volumes characteristic for the real technical environments, very intense release and long periods are required for a flammable gas to mix with the oxidizer atmosphere and form a highly active gaseous mixture with nearly stoichiometric composition. In the real case scenario, the presence of venting also efficiently reduces the concentration of combustible components in the gaseous mixture inside of the volume. Besides, the mixture possesses a vertical concentration gradient due to natural convection. If the reactive component is heavier than the atmospheric air, then the composition near the floor is richer than near the ceiling and vice versa for the light fuels such as hydrogen. Thus, the problem of precise determination of the lean flammability limits and the evolution of the flame in near-limit gaseous mixtures is of primary importance for the safety analysis of real facilities. The problem of assessment of lean flammability limit is one of the fundamental ones and was thoroughly studied. For now, a large amount of experimental data is available on concentration limits for a wide range of combustible compounds and conditions. However, the dynamics of near-limit flames and physical mechanisms determining the development of combustion in real large-scale environments under terrestrial gravity conditions are still poorly understood.
Flame structures developing in ultra-lean premixed combustible mixtures with scarce fuel content possess a set of peculiar features. Theoretical analysis of the ultra-lean flame in zero gravity conditions have shown that the key mechanism of the combustion stability is the diffusion of the deficient reactant to the reaction zone [2]. In contrast, for the deflagrative combustion intrinsic to the more chemically active mixtures, the main mechanisms of flame propagation are heat transfer and diffusion of radicals [3]. Experimental evidence of stable ultra-lean kernels, so-called “flame-balls”, was obtained for the first time in microgravity experiments [4,5,6]. Stable ultra-lean flame kernels were also experimentally observed and numerically investigated in the co-flow setup at terrestrial gravity conditions [7]. The evolution of ultra-lean flames in large-scale volumes at terrestrial gravity conditions that is of paramount interest for real technical applications is little discussed in the literature. In classical papers by Bohm an Clusius [8] and Babkin et al. [9], the formation of rising cap-shaped flames and more complex flame structures in flammability tubes filled with the lean mixtures of the near limit compositions were experimentally observed. In recent papers [10,11], it was shown that the non-confined flame in the ultra-lean mixture evolves into a complex multi-kernel structure after ignition due to the influence of buoyancy forces and gas-dynamical flows associated with hot products upward motion. Herewith, gas-dynamic flows formed on the initial stages of the flame development may cause the breakup of the flame structure that defines the formation of small-scale flame kernels that are not able to sustain stable combustion and tend to extinguish [12]. Thus, it is of importance to take into account buoyancy-driven flows associated with the flame kernel rise to determine flammability limits in the terrestrial gravity conditions.
Ultra-lean flame kernel dynamics is very similar to the one of gaseous bubble rising in liquid or dynamics of liquid drops in immiscible fluids [13]. Along with the features of the upward motion, such as the character of rising velocity change with time and the tendency to approach the terminal velocity value, the flame kernel also possesses similar structural behavior to the propagating gaseous bubbles or liquid drops in immiscible fluids [14,15]. The flame-bubble analogy was fruitfully adopted by V. Bychkov to analyze curved flames propagation and associated Rayleigh–Taylor instability development in horizontal tubes [16,17] and astrophysical flames [18].
In the present paper, we aim to study the ultra-lean flames by means of numerical modeling and classify the flame structure by employing flame-bubble analogy. Here, the large-scale structure of the ultra-lean flame in terrestrial gravity conditions is analyzed numerically in a range of near-limit compositions. For the first time, the flame-bubble analogy is extended to describe spatial flame structure and stability by means of correlation diagrams commonly used for the classification of bubbles and droplets [19].

2. Problem Setup

The computational domain represents a two-dimensional semi-confined space with the free outlet boundary condition imposed at the top and sides of the domain and adiabatic wall condition at the bottom (see Figure 1a). The volume is initially filled with a hydrogen-air mixture of near-limit composition (4.5–6.0% H 2 in the air) at normal conditions ( T 0 = 300 K, p 0 = 1 atm). The mixture is ignited near the adiabatic wall by the instantaneous temperature increase up to 2000 K. The ignition region is a circle with a radius R 0 = 1 4 mm. After the ignition, the flame propagates upward due to buoyancy forces. The problem setup resembles the setups we studied previously in [10,11]. Here, however, the flame dynamics is studied in a range of mixture compositions, and the effect of ignition mode is also considered.
To investigate the relations between flame dynamics and the evolution of the gaseous bubble, we performed calculations in the approximation of “immiscible gases”. The idea was to model the upward motion of the light gas bubble in the heavier atmosphere neglecting the effects of diffusion to sustain the density gradient between the bubble and the atmosphere. The atmosphere, as well as the contents of the bubble, consisted of the hydrogen-air mixture with different compositions defining different density ratios. The density ratio between the bubble and the surrounding gas was varied in a range θ b = 1.5 4 . The bubble had an initial circular shape with a radius R 0 , b = 1 4 mm. A schematic view of this problem setup is presented in Figure 1b.

3. Numerical Methods

The governing equations represent the full Navier-Stokes model with account of thermal conductivity, multicomponent diffusion, energy release associated with the chemical transformations, and gravity forces. Due to the low intensity of the process of near-limit combustion, it is reasonable to use low-Mach approximation to avoid strict limitations on the time step of numerical integration. The following equations constitute the mathematical model of the considered process [20]:
ρ t + ρ u i x i = 0
ρ Y k t + ρ Y k u i x i = ρ Y k V k , i + ω ˙ k
u t u × ω + H p ^ 1 ρ = 1 ρ ρ ρ 0 g + · σ
ρ h s t + ρ h s u i x i = d p ¯ d t k = 1 N ω ˙ k Δ h f , k 0 x i ρ k = 1 N h s , k Y k V k , i x i κ T x i + σ i j u i x j
σ i j = μ Y k , T u i x j + u j x i 2 3 δ i j u l x l
p ¯ = ρ R T k Y k M k
d h s = C p Y k , T d T
here ρ is the mass density, u is the mass velocity, u i are the mass velocity vector components, Y k is the mass fraction of k-th component of gaseous mixture, M k is the molecular weight of k-th component of gaseous mixture, p ^ is the dynamic component of pressure fluctuations, which is much smaller by the order of magnitude compare with thermodynamic pressure p ¯ and p x , t = p ¯ t + p ^ x , t , σ is the viscous stresses tensor, σ i j are the components of the viscous stresses tensor, H = u 2 / 2 + p ^ / ρ is the stagnation energy per unit mass, ω is the vorticity vector, g = 0 , 9.8 is the acceleration of gravity, h s is the specific enthalpy of the mixture, h s , k is the specific enthalpy of k-th component of gaseous mixture, T is the temperature, κ Y k , T is the thermal conductivity coefficient, μ Y k , T is the viscosity coefficient, h f , k 0 is the specific enthalpy of formation of k-th component of the gaseous mixture, V k , i is the diffusion velocity vector component of k-th specie, C p ( Y k , T ) is the specific heat capacity at constant pressure of the mixture. Term ω ˙ k represents the change in mass fraction of k-th specie due to the chemical reactions.
Data form the JANAF tables [21] was used to calculate specific heat capacities and enthalpies of formation. Diffusion was modeled in the zeroth-order Hirshfelder-Curtiss approximation [22]. Mixture averaged transport coefficients were obtained from the first principles of the gas kinetics theory [23]. Correction velocity approach proposed in [24] was implemented to calculate diffusion velocities. Chemical kinetics of the hydrogen-air mixture oxidation was modeled via detailed oxidation scheme with 21 elemental reaction and 9 species [25].
The governing Equations (1)–(7) were integrated via the second-order predictor/corrector method described in [26]. Computational cell size was equal 0.4 mm that allows obtaining convergent solution (see comparison between solutions obtained in computational domains with cells 0.2 mm and 0.4 mm in Figure 2). Flame thickness for the 6% H 2 -Air mixture, which is the most reactive mixture among the considered ones, is about 10 mm. Thus, the spatial flame front structure is resolved with at least 25 cells for the cell size Δ x = 0.4 mm that is sufficient for the convergence.

4. Results and Discussion

At first, let us consider the flame structures and dynamics of flame kernels in near-limit hydrogen-air mixtures. In Figure 3 spatial flame structures formed in the mixtures with compositions 5.0 6.0 % H 2 and ignited by the heated region with R 0 = 1 mm radius are presented. One can see that a small variation of the composition, less than 2% of hydrogen content in the mixture, leads to the drastic change in the ultra-lean flame structure. For the mixture with 4.5% hydrogen content in the air, successful flame propagation was not observed in the case of ignition via the heated region with R 0 = 1 mm. In the compound with 5% hydrogen content, a cap-shaped flame structure, formed on the early stage of the flame evolution, is almost immediately being broken up by the vortical flows generated by the rising combustion products [12]. That leads to the formation of two separate flame kernels, which continue to propagate upward under the buoyancy effect, and the combustion area only slightly expands with time. Increasing hydrogen content up to 5.5% results in a more complex flame structure. Here, similarly to the case of 5% H 2 in air, the initially formed cap-shaped flame is subjected to the stretching forces determined by the vortical flows. Here flame is broken up as well. However, the combustion intensity is high enough to produce a larger amount of combustion products and to maintain three separate combustion zones—on the leading tip of the cap and the periphery. The middle kernel demonstrates similar behavior as after the ignition. Its upward motion generates vortical flow in the wake. As a result, the secondary break-up proceeds in the leading flame kernel. This time the volume of the accumulated combustion products is not enough, and the flame breaks up in a similar scenario as in the 5% H 2 mixture. The overall evolution of the flame is very sensitive to gas-dynamics. Emerging gas flows stretch and easily tear apart the flame surface into smaller flame kernels. Nevertheless, due to higher burning rate, such break-up events do not cause the global quenching of the flame, and it continues to rise upward in the form of a separated multi-kernel structure. Exactly this multi-kernel combustion was presumably registered in [8]. Finally, the stable cap-shaped flame structure is formed in the hydrogen-air mixture with 6% hydrogen content. Here one can also observe break-up events, although they happen only on the periphery of the cap-shaped flame. Detailed analysis of this combustion mode was carried out in [10].
Analysis of the ultra-lean combustion dynamics under terrestrial gravity conditions allows concluding that gas-dynamical flows associated with natural convection and buoyancy of hot combustion products in cold fresh mixture atmosphere can inhibit flame propagation. The inhibition mechanism here is similar to that of the turbulent flows, where the excessive flow-induced stretch leads to the flame quenching [27]. However, turbulence can also cause flame intensification, especially in the mixtures with a pronounced preferential diffusion effects [28]. Thus, in ultra-lean hydrogen mixtures with Lewis number sufficiently lower than 1, one can expect an increase in burning velocity due to flame stretch and flame surface break-up.
Let us continue with the analysis of the ignition effect on the ultra-lean flame dynamics. In Figure 4, the comparison between flame dynamics in a 6% hydrogen-air mixture initiated by three different ignition modes ( R 0 = 1 , 2 , 4 mm) is demonstrated. One can see that general trends are quite similar. In all the cases the maximal temperature of the flame approaches super adiabatic value T m a x 1150 K. Additional 1D calculations in spherical symmetry provide value of maximal flame temperature T m a x , 1 D = 1104 K. Maximal flame temperature in two-dimensional flame can exceed values obtained in 1D setup due to gas-dynamical stretch and curvature of the flame front. The burning velocity can be assessed as the time derivative of the radius of the circular area equivalent to the area of the combustion products d R e q / d t . This evaluation provides an almost identical time history of the burning velocity throughout the process in each case. It can be noted that after ignition burning velocity tends to the value of ≈30 mm/s that is close to the normal burning velocity obtained in 1D simulations. However, after peripheral break-up and formation of additional small kernels, burning velocity increases. That is associated with the positive effect of curvature due to preferential diffusion. The main difference between three ignition modes is related to the initial cooling rate of the flame kernel and the initial value of the equivalent radius. The higher amount of the initially heated gas cools down slower and supports the increased rate of combustion for a longer time. That results in higher starting values of R e q . In the case of a 6% H 2 -air mixture that mechanism does not play a significant role. The situation changes dramatically closer to limit compositions. The mixture with 4.5% H 2 content, which is not flammable by the ignition source with R 0 = 1 mm, occurs to be able to burn when the area of ignition is increased. In Figure 5a, the flame structure in the 4.5% H 2 -air mixture is demonstrated. The flame consists of two separate kernels that do not grow in time but travel upwards in a zigzagging trajectory. Such rising mode can be also observed in case of small gaseous bubbles [29]. In Figure 5b, time histories of the maximum flame temperature and R e q are depicted. One can see that flame temperature changes in an oscillating way, approaching the maximal temperature obtained in 1D spherical calculations ( T m a x , 1 D ) in the lowest position. R e q shows similar behavior and oscillates near the average value of ≈1.5 mm. The excessive amount of combustion products generated in the early stage of the flame development due to the increased ignition region allows sustaining combustion in the 4.5% H 2 -air mixture.
To analyze flame dynamics from the flame-bubble analogy point of view, we performed calculations of two dimensionless factors, namely Reynolds number Re = ρ a m b U d e q μ a m b and Froude number Fr = U g d e q . The Reynolds number reflects the effect of inertial and viscous forces, while the Froude number indicates the ratio of inertial and gravity forces. Such a set of dimensionless parameters is commonly used for the analysis of bubble dynamics in liquid [14]. Characteristic velocity U was calculated as a difference between the rising velocity U f , L and burning velocity estimated as d R e q / d t . For bubbles, characteristic velocity was equal to the rising velocity of the bubble. In Figure 6, the diagrams obtained in the Re-Fr plane are given. At first, we performed calculations of dimensionless quantities for bubbles with different initial radii R 0 , b and density ratios θ b . It occurs that the trajectory of the bubble in the Re-Fr plane represents a line segment during the initial accelerated upward motion (see Figure 6a). Convective mixing smooths out the density gradient, and the bubble area decreases until it ceases to be distinguishable. A convective decrease in bubble area results in a deviation from the linear trajectory in the Re-Fr plane as the Reynolds number increases with the bubble area, while the Froude number, on the contrary, decreases. That process is inevitable due to the absence of surface tension. Convective mixing is more pronounced with θ b increased, as the density ratio defines the initial buoyant acceleration of the bubble. Thus, the length of the linear segment decreases with increasing θ b . At θ b = 5 , linear segment ceases to exist as the convective mixing starts to dominate the structure of the bubble. Therefore, we plotted linear segments of the initial bubble dynamics for density ratios θ b = 1.5 4 . Linear segments of bubbles with the same initial radius R 0 , b occur to have the same slope, while density ratio defines the starting point of the linear segment and its length. All the bubbles with the same initial radius R 0 , b develop similar spatial structures that allowed us to divide the Re-Fr plane into areas defined by the slopes of linear segments (see Figure 6) where one expects to observe one or another structure of the bubble. The obtained diagram can be used to predict and classify the spatial structure of bubbles and gaseous flame kernels. Comparing flame structures in Figure 3 and bubble structures in Figure 6b, one can note a clear resemblance. In the region near the R 0 , b = 1 mm line, bubbles loose connectivity and acquire a two-kernel structure, similar to the flame modes observed in 5.0% and 4.5% H 2 -Air mixtures. On the other hand, line R 0 , b = 4 mm defines the region with stable cap-shaped bubbles. Lines R 0 , b = 2 mm and R 0 , b = 3 mm correspond to unstable bubbles susceptible to the development of Rayleigh-Taylor instability on the interface [30].
For the flame kernels, trajectories on the Re-Fr plane are not linear (see arrows in Figure 6b). Flame trajectories possess two distinct attraction limits, namely low Fr < 1 and high Re and high Fr > 1 and low Re . The low value of the Froude number Fr < 1 is characteristic to stable modes of flame kernel evolution insusceptible to Rayleigh-Taylor instability development in terrestrial gravity conditions such as multi-kernel (5.5% H 2 -Air) or cap-shaped (6.0% H 2 -Air) flames. For that kind of flame, effective combustion area rises together with the rising velocity, and rising velocity increases slower or with the same increment as the square of combustion products area. Apparently, that is the case for the mixtures with high burning velocity. On the other hand, the second limit of high Fr and low Re attracts flames of leaner compositions characterized by modes with small kernels (5.0% H 2 -Air and 4.5% H 2 -Air with R 0 = 2 mm). In such flames, convection causes a strong quenching effect. Therefore, each break-up event results in a major loss in the effective flame area. So the equivalent radius of the flame rapidly reduces, and the Froude number increases. It is interesting to note that the combustion process in 5.5% H 2 -Air has a loop in trajectory in the Re-Fr plane. The trajectory enters the region with the possibility of flame break-up due to Rayleigh-Taylor instability. That behavior reflects the formation of a multi-kernel structure on the initial stages of flame evolution, during which the convective quenching effect is manifested most strongly, and the effective combustion area substantially decreases. After the initial impact, the flame recovers and proceeds the propagation without significant opposition.
It should be noted that the rising velocity of the flame kernel is determined by the expansion ratio of the flame, which is almost constant in the considered range of hydrogen content. Thus, it is the burning velocity and the ability to withstand convective cooling during the break-up events that determine the difference in trajectories in the Re-Fr plane for different mixtures. That explains the shift in the flammability limit with the increase of the initial radius of the ignition area R 0 , b . The higher flame area is less affected by the convective cooling during the initial break-up of the flame in 4.5% hydrogen-air mixture. Therefore, combustion is still possible after the break-up. The same effect can be obtained with the increase of ignition temperature (intensity of the ignition source) instead of enlargement of the area of the ignition zone.

5. Conclusions

In the present paper, we carried out numerical analysis of the ultra-lean hydrogen-air flame dynamics in terrestrial gravity conditions. It is revealed that there is a wide range of possible spatial flame structures forming during the combined process of convective motion and burning of the gaseous mixture. The flame structure, as well as flammability limit, are strongly affected by gas-dynamical flows developed in the process of the flame upward motion under the action of buoyancy in terrestrial gravity conditions. In near-limit compositions, gas-dynamical stretch associated with vortical flows formed in the wake of rising combustion products can result in flame break-up and quenching. On the other hand, due to the effects of preferential diffusion, distortion of the flame front in lean hydrogen-air mixture by gas-dynamical flows can lead to the increase in combustion temperature and amplification of burning velocity. The developed approach of bubble-flame analogy allowed performing classification of different flame modes in near-limit mixtures. Based on computations of non-reactive gaseous bubbles, regions of stable and unstable flame modes were distinguished on the Re-Fr plane of dimensionless Reynolds and Froude number. The convection and ignition mode effects on the flammability limits are interpreted through the analysis of the flame trajectories in the Re-Fr plane.

Author Contributions

Conceptualization, A.K. and I.Y.; methodology, A.K., K.M., I.Y.; software, I.Y.; validation, K.M., I.Y.; formal analysis, A.K., K.M., I.Y.; investigation, K.M., I.Y.; resources, I.Y.; data curation, K.M., I.Y.; writing—original draft preparation, I.Y.; writing—review and editing, A.K., K.M., and I.Y.; visualization, K.M. and I.Y.; supervision, I.Y.; project administration, I.Y.; funding acquisition, I.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The work was financially supported by the state support of young Russian scientists grant MK-3473.2019.2.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. Obtained amount of data is too large to make it publicly available.

Acknowledgments

We acknowledge high-performance computing support from the Joint Supercomputer Center of the Russian Academy of Sciences and Supercomputing Center of Lomonosov Moscow State University.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
JANAFJoint Army, Navy and Air Force

References

  1. Ciccarelli, G.; Dorofeev, S. Flame acceleration and transition to detonation in ducts. Prog. Energy Combust. Sci. 2008, 34, 499–550. [Google Scholar] [CrossRef]
  2. Zel’dovich, Y.B.; Branblatt, G.; Librovich, V.B.; Makhviladze, G. The Mathematical Theory of Combustion and Explosions, 1st ed.; Consultants Bureau: New York, NY, USA, 1985. [Google Scholar]
  3. Warnatz, J.; Maas, U.; Dibble, R. Combustion, 4th ed.; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2005. [Google Scholar]
  4. Ronney, P.D. Near-limit flame structures at low Lewis number. Combust. Flame 1990, 82, 1–14. [Google Scholar] [CrossRef]
  5. Ronney, P.D.; Wu, M.S.; Pearlman, G.H.; Karen, J.W. Experimental Study of Flame Balls in Space: Preliminary Results from STS-83. AIAA J. 1998, 36, 1361–1368. [Google Scholar] [CrossRef]
  6. Kagan, L.; Sivashinsky, G. Self-fragmentation of nonadiabatic cellular flames. Combust. Flame 1997, 108, 220–226. [Google Scholar] [CrossRef]
  7. Shoshin, Y.; van Oijen, J.; Sepman, A.; de Goey, L. Experimental and computational study of the transition to the flame ball regime at normal gravity. Proc. Combust. Inst. 2011, 33, 1211–1218. [Google Scholar] [CrossRef]
  8. Böhm, G.; Clusius, K. Die Struktur aufsteigender H2-O2-Flammens. Z. Naturforsch. 1948, A3, 386. [Google Scholar] [CrossRef]
  9. Babkin, V.S.; Zamashchikov, V.V.; Badalyan, A.M.; Krivulin, V.N.; Kudryavtsev, E.A.; Baratov, A.N. Effect of tube diameter on homogeneous gas flame propagation limits. Combust. Explos. Shock Waves 1982, 18, 164–171. [Google Scholar] [CrossRef]
  10. Yakovenko, I.S.; Ivanov, M.F.; Kiverin, A.D.; Melnikova, K.S. Large-scale flame structures in ultra-lean hydrogen-air mixtures. Int. J. Hydrogen Energy 2018, 43, 1894–1901. [Google Scholar] [CrossRef] [Green Version]
  11. Volodin, V.V.; Golub, V.V.; Kiverin, A.D.; Melnikova, K.S.; Mikushkin, A.Y.; Yakovenko, I.S. Large-scale Dynamics of Ultra-lean Hydrogen-air Flame Kernels in Terrestrial Gravity Conditions. Combust. Sci. Technol. 2020. [Google Scholar] [CrossRef]
  12. Kiverin, A.D.; Yakovenko, I.S.; Melnikova, K.S. On the structure and stability of ultra-lean flames. J. Phys. Conf. Ser. 2019, 1147, 012048. [Google Scholar] [CrossRef]
  13. Bychkov, V.; Liberman, M. Dynamics and stability of premixed flames. Phys. Rep. 2000, 325, 115–237. [Google Scholar] [CrossRef]
  14. Ohta, M.; Sussman, M. The buoyancy-driven motion of a single skirted bubble or drop rising through a viscous liquid. Phys. Fluids 2012, 24. [Google Scholar] [CrossRef] [Green Version]
  15. Ray, B.; Prosperetti, A. On skirted drops in an immiscible liquid. Chem. Eng. Sci. 2014, 108, 213–222. [Google Scholar] [CrossRef]
  16. Bychkov, V.V. Bubble motion in a horizontal tube and the velocity estimate for curved flames. Phys. Rev. E 1997, 55, 6898–6901. [Google Scholar] [CrossRef]
  17. Modestov, M.; Bychkov, V.; Betti, R.; Eriksson, L.E. Bubble velocity in the nonlinear Rayleigh–Taylor instability at a deflagration front. Phys. Plasmas 2008, 15, 042703. [Google Scholar] [CrossRef]
  18. Bychkov, V.; Modestov, M.; Akkerman, V.; Eriksson, L.E. The Rayleigh—Taylor instability in inertial fusion, astrophysical plasma and flames. Plasma Phys. Control. Fusion 2007, 49, B513–B520. [Google Scholar] [CrossRef]
  19. Bhaga, D.; Weber, M.E. Bubbles in viscous liquids: Shapes, wakes and velocities. J. Fluid Mech. 1981, 105, 61–85. [Google Scholar] [CrossRef] [Green Version]
  20. Kuo, K.K.; Acharya, R. Fundamentals of Turbulent and Multiphase Combustion, 1st ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012. [Google Scholar]
  21. Chase, M.W. NIST-JANAF Themochemical Tables. J. Phys. Chem. Ref. Data Monogr. 1998, 9, 1–1963. [Google Scholar]
  22. Hirschfelder, J.O.; Curtiss, C.F.; Bird, R.B. The Molecular Theory of Gases and Liquids; Revised ed.; Wiley-Interscience: Hoboken, NJ, USA, 1964. [Google Scholar]
  23. Kee, R.J.; Coltrin, M.E.; Glarborg, P. Chemically Reacting Flow: Theory and Practice, 1st ed.; Wiley-Interscience: Hoboken, NJ, USA, 2003. [Google Scholar]
  24. Coffee, T.; Heimerl, J. Transport algorithms for premixed, laminar steady-state flames. Combust. Flame 1981, 43, 273–289. [Google Scholar] [CrossRef]
  25. Kéromnès, A.; Metcalfe, W.K.; Heufer, K.A.; Donohoe, N.; Das, A.K.; Sung, C.J.; Herzler, J.; Naumann, C.; Griebel, P.; Mathieu, O.; et al. An experimental and detailed chemical kinetic modeling study of hydrogen and syngas mixture oxidation at elevated pressures. Combust. Flame 2013, 160, 995–1011. [Google Scholar] [CrossRef] [Green Version]
  26. McGrattan, K.; McDermott, R.; Hostikka, S.; Floyd, J.; Vanella, M. Fire Dynamics Simulator Technical Reference Guide; Volume 1: Mathematical Model; Technical Report NIST Special Publication 1018-1; U.S. Department of Commerce, National Institute of Standards and Technology: Gaithersburg, MD, USA, 2019. [CrossRef]
  27. Chomiak, J.; Jarosiński, J. Flame quenching by turbulence. Combust. Flame 1982, 48, 241–249. [Google Scholar] [CrossRef]
  28. Yang, S.; Saha, A.; Liang, W.; Wu, F.; Law, C.K. Extreme role of preferential diffusion in turbulent flame propagation. Combust. Flame 2018, 188, 498–504. [Google Scholar] [CrossRef]
  29. Tripathi, M.K.; Sahu, K.C.; Govindarajan, R. Dynamics of an initially spherical bubble rising in quiescent liquid. Nat. Commun. 2015, 6, 6268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Kull, H. Theory of the Rayleigh-Taylor instability. Phys. Rep. 1991, 206, 197–325. [Google Scholar] [CrossRef]
Figure 1. Schematic of the problem setups. (a) Ultra-lean combustion of hydrogen-air mixture initiated near the bottom wall of the domain via instantaneous local temperature rise. (b) Rise of the light gaseous bubble under terrestrial gravity conditions in the “immiscible gases” approximation.
Figure 1. Schematic of the problem setups. (a) Ultra-lean combustion of hydrogen-air mixture initiated near the bottom wall of the domain via instantaneous local temperature rise. (b) Rise of the light gaseous bubble under terrestrial gravity conditions in the “immiscible gases” approximation.
Fluids 06 00021 g001
Figure 2. Grid sensitivity analysis for the combustion process in 6% H 2 -Air mixture. Temperature isolines T = 1000 K are presented for computational domains with cell size Δ x = 0.2 mm (solid black lines) and Δ x = 0.4 mm (orange dashed lines).
Figure 2. Grid sensitivity analysis for the combustion process in 6% H 2 -Air mixture. Temperature isolines T = 1000 K are presented for computational domains with cell size Δ x = 0.2 mm (solid black lines) and Δ x = 0.4 mm (orange dashed lines).
Fluids 06 00021 g002
Figure 3. Flame structures in near-limit ultra-lean hydrogen-air mixtures. Top—5.0% H 2 -Air; middle—5.5% H 2 -Air; bottom—6.0% H 2 -Air. Color palette indicates gas temperature.
Figure 3. Flame structures in near-limit ultra-lean hydrogen-air mixtures. Top—5.0% H 2 -Air; middle—5.5% H 2 -Air; bottom—6.0% H 2 -Air. Color palette indicates gas temperature.
Fluids 06 00021 g003
Figure 4. Time histories of the maximal temperature, equivalent radius and d R e q / d t in 6.0% hydrogen-air mixture. Solid lines—ignition region with R 0 = 1 mm, dash-dotted lines—ignition region with R 0 = 2 mm, dashed lines—ignition region with R 0 = 4 mm.
Figure 4. Time histories of the maximal temperature, equivalent radius and d R e q / d t in 6.0% hydrogen-air mixture. Solid lines—ignition region with R 0 = 1 mm, dash-dotted lines—ignition region with R 0 = 2 mm, dashed lines—ignition region with R 0 = 4 mm.
Fluids 06 00021 g004
Figure 5. (a) Flame structure in ultra-lean 4.5% hydrogen-air mixture ignited by the heated region with R 0 = 2 mm. (b) Time histories of the maximal temperature and equivalent radius in 4.5% hydrogen-air mixture.
Figure 5. (a) Flame structure in ultra-lean 4.5% hydrogen-air mixture ignited by the heated region with R 0 = 2 mm. (b) Time histories of the maximal temperature and equivalent radius in 4.5% hydrogen-air mixture.
Fluids 06 00021 g005
Figure 6. The diagram in Re-Fr plane. (a) Linear segments obtained via non-reactive bubbles computations. Linear approximations of those segments divide the plane into five sectors. Initial radius of the bubble defines the slope of the linear segment, density ratio defines starting point and length of the linear segment. Arrows indicate the direction of density ratio θ b increase. Linear segments colors: white— θ b = 1.5 , yellow— θ b = 2 , blue— θ b = 3 , red— θ b = 4 . (b) Flame trajectories in Re-Fr plane for different mixture compositions: green—6% H 2 -Air, R 0 = 1 mm, black—5.5% H 2 -Air, R 0 = 1 mm, orange—4.5% H 2 -Air, R 0 = 2 mm, blue—5.0% H 2 -Air, R 0 = 1 mm. Characteristic bubble structures for R 0 , b = 1 4 mm for each sector are presented outside the Re-Fr plane.
Figure 6. The diagram in Re-Fr plane. (a) Linear segments obtained via non-reactive bubbles computations. Linear approximations of those segments divide the plane into five sectors. Initial radius of the bubble defines the slope of the linear segment, density ratio defines starting point and length of the linear segment. Arrows indicate the direction of density ratio θ b increase. Linear segments colors: white— θ b = 1.5 , yellow— θ b = 2 , blue— θ b = 3 , red— θ b = 4 . (b) Flame trajectories in Re-Fr plane for different mixture compositions: green—6% H 2 -Air, R 0 = 1 mm, black—5.5% H 2 -Air, R 0 = 1 mm, orange—4.5% H 2 -Air, R 0 = 2 mm, blue—5.0% H 2 -Air, R 0 = 1 mm. Characteristic bubble structures for R 0 , b = 1 4 mm for each sector are presented outside the Re-Fr plane.
Fluids 06 00021 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yakovenko, I.; Kiverin, A.; Melnikova, K. Ultra-Lean Gaseous Flames in Terrestrial Gravity Conditions. Fluids 2021, 6, 21. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6010021

AMA Style

Yakovenko I, Kiverin A, Melnikova K. Ultra-Lean Gaseous Flames in Terrestrial Gravity Conditions. Fluids. 2021; 6(1):21. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6010021

Chicago/Turabian Style

Yakovenko, Ivan, Alexey Kiverin, and Ksenia Melnikova. 2021. "Ultra-Lean Gaseous Flames in Terrestrial Gravity Conditions" Fluids 6, no. 1: 21. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6010021

Article Metrics

Back to TopTop