Next Article in Journal
Research on Optimal Scheduling Strategy of Differentiated Resource Microgrid with Carbon Trading Mechanism Considering Uncertainty of Wind Power and Photovoltaic
Next Article in Special Issue
Modeling of Hydrogen Combustion from a 0D/1D Analysis to Complete 3D-CFD Engine Simulations
Previous Article in Journal
Participation of Energy Communities in Electricity Markets and Ancillary Services: An Overview of Successful Strategies
Previous Article in Special Issue
A Comparative Study of the Hydrogen Auto-Ignition Process in Oxygen–Nitrogen and Oxygen–Water Vapor Oxidizer: Numerical Investigations in Mixture Fraction Space and 3D Forced Homogeneous Isotropic Turbulent Flow Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Water Mist on the Initial Evolution of Turbulent Premixed Hydrogen/Air Flame Kernels

by
Riccardo Concetti
1,*,
Josef Hasslberger
1,
Nilanjan Chakraborty
2 and
Markus Klein
1,*
1
Department of Aerospace Engineering, University of the Bundeswehr Munich, Werner-Heisenberg-Weg 39, 85577 Neubiberg, Germany
2
School of Engineering, Newcastle University, Newcastle-upon-Tyne NE1 7RU, UK
*
Authors to whom correspondence should be addressed.
Energies 2024, 17(18), 4632; https://doi.org/10.3390/en17184632
Submission received: 16 August 2024 / Revised: 4 September 2024 / Accepted: 6 September 2024 / Published: 16 September 2024
(This article belongs to the Special Issue Towards Climate Neutral Thermochemical Energy Conversion)

Abstract

:
In this study, a series of carrier-phase direct numerical simulations are conducted on spherical expanding premixed hydrogen/air flames with liquid water addition. An Eulerian–Lagrangian approach with two-way coupling is employed to describe the liquid–gas interaction. The impacts of preferential diffusion, the equivalence ratio, water loading, and the initial diameter of the water droplets are examined and analyzed in terms of flame evolution. It is observed that liquid water has the potential to influence flame propagation characteristics by reducing the total burning rate, flame area, and burning rate per unit area, attributed to flame cooling effects. However, these effects become discernible only under conditions where water evaporation is sufficiently intense. For the conditions investigated, the influence of preferential diffusion on flame evolution is found to be more significant than the interaction with liquid water. The results suggest that due to the slow evaporation rate of water, which is a result of its high latent heat of evaporation, the water droplets do not disturb the initial flame kernel growth significantly. This has implications for water injection concepts in internal combustion engines and for explosion mitigation.

1. Introduction

The utilization of liquid water injection in combustion applications has obtained considerable attention within the scientific community owing to its numerous practical advantages. The interaction between liquid water and the flame is based on three primary effects [1,2,3]:
  • Cooling effect: This results from the extraction of heat from the hot gases during the evaporation process.
  • Dilution effect: The presence of water acts as a diluent for the fuel and oxidizer, leading to a reduction in the heat released per unit of mass in the system.
  • Chemical effect: Water exhibits high efficiency in third-body reactions, influencing the flame structure (this effect is not analyzed in the current study, as it is considered a high-order effect).
The application of water injection in internal combustion engines (automotive and aeronautical engines) has been implemented since the last century (e.g., Pratt and Whitney J57 engine and the 1962 Oldsmobile Jetfire). Historical motivations for such applications include temporary increases in thrust and power output due to the augmented mass processed by the engine. Additional benefits encompass reduced maximum temperatures—in this way, the liquid water acts as thermal protection for structures—and decreased NOx emissions, as these are generally produced at elevated temperatures [4]. Previous research has demonstrated the impact of water injection on emissions in diesel engines and gas turbines [5,6,7] and also showed enhanced engine performance. Liquid water injection also holds significance in safety applications, as experimental and numerical studies have illustrated its ability to suppress or mitigate explosions and fires [8,9].
Nowadays, a primary focus of combustion research involves reducing CO2 emissions. In this context, investigations into the combustion of carbon-free fuels, particularly hydrogen, have intensified. Hydrogen combustion offers several advantages, including high availability and inherent renewability [10,11], which can contribute to the development of a future green and circular economy [12,13]. Additionally, hydrogen has the advantage of a high energy density per unit mass (141 MJ/kg). Despite these advantages, challenges such as low density, high volatility, wide flammability range, and safety concerns arise, which were previously encountered in rocket propulsion, for instance, and are now being considered for distribution and storage. The high diffusivity of hydrogen and tendency for preferential diffusion introduce complexities, leading to intrinsic instabilities, particularly thermodiffusive instability. This instability occurs when the effective Lewis number ( L e e f f ) of the mixture falls below a critical value, as characterized by the dispersion relation proposed by Matalon et al. [14]. Another criterion for thermodiffusive instability is the Markstein length, representing the effect of the stretch rate on the flame’s displacement speed; a negative value indicates susceptibility to thermodiffusive instability [15].
The propensity of mixtures with Lewis numbers below a critical threshold to exhibit intrinsic flame instability has been extensively studied [14,15,16,17,18,19,20,21,22,23] and remains a focal point of community interest. Spherical flames, which are typical configurations in internal combustion engines, enable the analysis of statistically non-zero stretch rates on flame velocities and the consideration of flame propagation characteristics (e.g., Markstein length). The evolution of expanding spherical flames becomes more intricate in the presence of preferential diffusion, as the close correlation between diffusion and topology introduces side effects [24], as the local focusing (or defocusing) of reactants can generate super-adiabatic temperature regions. Moreover, when spherical flames grow and reach a certain radius, unstable cellular structures begin to form without external forcing. However, the current study excludes consideration of this scenario due to the computational constraints of our setup. The primary emphasis of the current study is on examining the development of turbulent spherical flames in the presence of preferential diffusion and an evaporating liquid phase. The introduction of non-unity Lewis number conditions coupled with a non-zero mean curvature results in a distinct evolution compared to that of a statistically planar flame [25,26]. Moreover, the interaction between the flame kernel and turbulence amplifies the growth rate of the average flame radius, as evidenced in a prior investigation [27]. A growth rate proportional to the radius, rather than a constant rate over time, is observed. While there are instances in the literature where analyses have explored expanding flames interacting with turbulence [27,28] and have considered the impact of non-unity Lewis numbers [29,30], the authors are not aware of any studies comprehensively examining the simultaneous effects of turbulence, preferential diffusion, and the evaporation of liquid water on the evolution and propagation of spherical flames across various equivalence ratios, Lewis numbers, and water addition conditions.
The subsequent two sections will provide a presentation of the mathematical background and numerical implementation of our study. Section 4 will delve into the discussion of results and observations, followed by the presentation of conclusions in the final section.

2. Mathematical Background

The method employed in this study to characterize the evolution of flame kernels in the presence of monodispersed water droplets utilizes a hybrid Eulerian–Lagrangian approach with two-way coupling. This approach entails representing the liquid phase as individual droplets treated as point sources and tracing their motion and evolution within the domain (Lagrangian phase). At the same time, the gas phase is characterized through standard conservation equations, with conditions for different fields computed at each grid point (Eulerian phase). Two-way coupling implies that the interaction between the two phases is considered bidirectionally, which means terms related to the other phase are present in the equations for both phases. It is noteworthy that, unlike four-way coupling, which considers particle–particle interaction, two-way coupling neglects this phenomenon. This omission is justifiable due to the sufficiently high average interdistance between the droplets ( s d / η = 1.45 1.93 , where s d is the average interdistance and η is the Kolmogorov length scale), which guarantees the low probability of droplet–droplet collision. Moreover, the Lagrangian description of the particles is possible because the initial diameter of the droplets ( a d ) is significantly smaller than the grid resolution ( a d / Δ x = 0.384 0.512 ) and the Kolmogorov length scale ( a d / η = 0.09 0.12 ). This allows us to avoid resolving the flow at the particle–gas interface [31,32,33]. It has been demonstrated that the present methodology accurately captures the evaporation characteristics in multiphase simulations [34] when the droplet dimension is below the Kolmogorov length scale. Furthermore, in a study by De Chaisemartin et al. [35] that analyzed the combustion of dispersed droplets with flames using both a Eulerian–Lagrangian approach and a fully droplet-resolved approach, similar behavior using the two methods was observed. However, the computational cost of the Eulerian–Lagrangian approach was shown to be 10 times lower than that of the fully Eulerian method.
In the Eulerian–Lagrangian formulation, for individual droplets, their positions, velocities, temperatures, and diameters ( x d , u d , T d , and a d , respectively) are calculated over time using the following expressions:
d x d d t = u d ; d u d d t = u ( x d , t ) u d τ d u ; d a d 2 d t = a d 2 τ d p ; d T ^ d d t = T ^ ( x d , t ) T ^ d B d L v / C p g τ d T .
In the various equations presented, terms associated with the gas field are evident, such as the gas velocity u ( x d , t ) and the gas temperature T ^ ( x d , t ) . The droplet temperature equation assumes infinite thermal conductivity for the liquid phase, resulting in a uniform temperature across the particle equal to T ^ d . Notably, this equation incorporates the latent heat of vaporization ( L v ), the gaseous specific heat ( C p g ), and the thermal Spalding number ( B d ), which is considered identical to the mass transfer Spalding number. The Spalding number is determined as B d = ( Y W S Y W ( x d , t ) ) / ( 1 Y W S ) . The water concentration at the interface Y S W is computed based on equilibrium conditions through the partial pressure derived from the Clausius–Clapeyron equation:
Y W S = 1 + W g W l P ( x d , t ) P W S 1 1 ; P W S = P r e f e x p L v R 1 T ^ r e f S 1 T ^ d S .
The atmospheric pressure ( P r e f ) and boiling temperature at that pressure ( T ^ r e f S ) serve as references in this context. The molar masses of the gas phase ( W g ) and liquid phase ( W l ) are denoted as such, with T ^ d S representing the particle temperature. Equation (1) incorporates three distinct relaxation times ( τ d u , τ d p , and τ d T ) associated with the evolution of particle velocity, diameter, and temperature, respectively. The definitions of these relaxation times are as follows:
τ d u = ρ d a d 2 18 C u μ ; τ d p = ρ d a d 2 4 μ S c S h c l n ( 1 + B d ) ; τ d T = ρ d a d 2 6 μ P r N u c B d l n ( 1 + B d ) C p l C p g .
The relaxation times depend on the properties of both the gaseous and liquid phases. In Equation (3), ρ d and C p l represent the mass density and specific heat of the droplets, respectively. The dynamic viscosity of the gaseous phase is denoted as μ , while S c and P r are the Schmidt and Prandtl numbers, respectively, for the gaseous phase. These values are considered identical in Equation (3). Similarly, the Sherwood and Nusselt numbers are identical and are calculated as S h c = N u c = 2 + 0.555 R e d S c / ( 1.232 + R e d S c 4 / 3 ) 1 / 2 . Finally, C u is an aerodynamic resistance coefficient computed as C u = 1 + R e d 2 / 3 / 6 , where R e d is the Reynolds number of the droplet ( R e d = ρ | u d u ( x d , t ) | a d / μ ). The general conservation equation for the gaseous phase is:
( ρ φ ) t + ( ρ u j φ ) x j = x j R φ φ 1 x j + ω ˙ φ + S ˙ g + S ˙ φ
The variable φ is a general variable and, depending on the equation, assumes a value of (1, u i , e, Y F , Y O , or Y W ), while φ 1 represents (1, u i , T ^ , Y F , Y O , or Y W ) as appropriate. Here, e represents the specific stagnation internal energy. The transport coefficients R φ in the equation are treated as constants with respect to temperature and composition, and their values are derived from the dynamic viscosity μ . The term S ˙ g is a source/sink term in the gas phase equation: for instance, representing the pressure gradient in the momentum equation. On the other hand, the term S ˙ φ serves as the source/sink term associated with the interaction with droplets and is calculated as follows:
S ˙ φ = 1 Δ V d d ( m d φ d ) d t
This term is computed individually for each droplet, and its impact is subsequently distributed among the nearest grid points. Mass diffusion of distinct species is addressed through the constant Lewis number approach, with unity Lewis numbers considered for oxygen and water. For hydrogen, the Lewis number ( L e H 2 ) is set at 0.2851 , which corresponds to unburned hydrogen’s Lewis number at a very low equivalence ratio as computed with Cantera [36]. To isolate the effects of preferential diffusion, additional simulations were conducted with an artificial hydrogen possessing identical characteristics as real hydrogen but assigned a Lewis number of unity.
While three-dimensional direct numerical simulations (DNSs) with detailed chemistry are currently feasible [37,38], such simulations remain computationally demanding, particularly in the context of an extensive parametric study. In this work, an irreversible single-reaction chemical mechanism is therefore considered [39]. The approach involves adjusting the Zeldovich number ( β ) and the heat of combustion (H) based on the local equivalence ratio ( ϕ ) [39]. The expressions for the Zeldovich number and heat release are as follows [39]:
β = β s t [ 1.0 + 8.250 ( ϕ 0.64 ) 2 ] if ϕ 0.64 β s t if 0.64 ϕ 1.07 β s t [ 1.0 + 1.443 ( ϕ 1.07 ) 2 ] if ϕ 1.07
H = H s t if ϕ 1.0 H s t [ 1.0 0.18 ( ϕ 1.0 ) ] if ϕ 1.0
The reaction rate is then calculated via an Arrhenius-type formulation:
ω ˙ F = ρ B * Y F Y O e x p β ( 1 T ) 1 α ( 1 T )
Here and in the following, T indicates the non-dimensional temperature defined as T = ( T ^ T ^ 0 ) / ( T ^ a d , ( ϕ = 1 ) T ^ 0 ) . In the equation of the chemical mechanism, the heat release parameter, denoted as α , is defined by the expression α = ( T ^ a d , ( ϕ = 1 ) T ^ 0 ) / T ^ a d , ( ϕ = 1 ) . With the selection of T ^ 0 = 300 K, the heat release parameter has a value of α = 0.874 . This single-step reaction mechanism has demonstrated its ability to capture the variation in laminar burning velocity with the gaseous equivalence ratio obtained from detailed chemistry simulations [1]. Several results presented in this article are based on a progress variable denoted as c. This variable increases monotonically from zero under unburned conditions (u) and unity in the burned gas under equilibrium conditions (b). In this context, the progress variable is specifically defined as the normalized oxygen mass fraction:
c = Y O u Y O Y O u Y O b
This definition has already been used in many previous works [1,2,3,40].

3. Numerical Implementation

Simulations were conducted utilizing a 3D compressible reactive flow solver: SENGA+ [1,2,3,40]. Spatial discretization employed a 10th-order central finite difference scheme, transitioning to a second-order one-sided finite difference scheme at the domain boundaries. Temporal advancement utilized an explicit third-order low-storage Runge–Kutta method [41]. In this study, a series of spherical flames with varying equivalence ratios ( ϕ = 0.8 and 1.0 ) were simulated, both with and without water addition ( Y W o v = 0 and 0.1 , where Y W o v is the ratio of added water to the total mass, Y W o v = Y W g + Y W l = m W / ( m 0 + m W ) , where m W is the added mass of water and m 0 = m A + m F is the sum of the fuel m F and air m A masses in the unburned part of the simulated volume). Two different types of water droplets, characterized by a d / δ s t = 0.03 and 0.04 , were used for water addition. The thickness of the unstretched stoichiometric flame, denoted as δ s t , is defined as:
δ s t = T ^ a d , ( ϕ = 1 ) T ^ 0 m a x ( | T ^ | L )
The subscript L designates laminar unstretched conditions. The initial droplet temperature is set equal to the unburned gas temperature T ^ 0 for all the cases. Although the initial temperature of the droplets affects the heat needed to reach the water saturation temperature, its effect is much weaker than the cooling effect caused by water evaporation due to large magnitude of the latent heat of evaporation for water. The different parameters for the simulations are presented in Table 1. The influence of a non-zero mean curvature, which is the case for spherical expanding flames, significantly affects the development and propagation of the flame, particularly under non-unity Lewis number conditions [25,26]. The method employed for initializing the flame kernel involves adopting a 1D laminar flame profile with a predetermined equivalence ratio and distributing it in a spherical shape. The kernel is positioned at the center of the domain, and conditions are set as fully burnt up to r = 0.2 δ s t . Beyond this point, the remaining 1D flame profile is mapped spherically. However, this procedure introduces discrepancies in thermochemical fields, especially in the presence of preferential diffusion, since the presence of a non-zero average curvature modifies the evolution of the local thermochemical fields due to the effects of diffusion. In the current dataset, the case with ϕ = 1 and L e H 2 < 1 is the one that evolves the fastest at the same stoichiometric chemical timescale ( t c h e m = α t h / S L , ( ϕ = 1 ) , with α t h as the thermal diffusivity of the mixture). The stoichiometric L e H 2 < 1 flame is allowed to evolve under laminar conditions until the thermochemical fields are relaxed, and this establishes a target flame radius of around 10 δ s t . In this way, all cases are restarted under laminar conditions. Turbulence is superimposed when the target flame radius is reached, and then the simulations are continued for 2 t c h e m . The turbulent field, which follows the Batchelor–Townsend energy spectrum, has an integral length scale L 11 / δ s t = 2.5 and intensity u / S L , ( ϕ = 1 ) = 4.0 , which yields a flame in the thin-reaction zone regime. Water addition is initiated at the beginning when the conditions are still laminar throughout the entire domain except for the central part containing the ignited flame. The decision to introduce droplets at the start of the simulations rather than when turbulence is initiated does not significantly impact the subsequent evolution of the flame or the evaporation process. This is a consequence of the fact that the primary mass and heat transfer processes between phases generally occur in the hotter regions of the flame (i.e., the flame and burned gas region), and the droplets are initially positioned outside of these areas. Boundary conditions are established as standard non-reflective inflow–outflow (NSCBC), as there are no pronounced flow directions in this case, in contrast to previous configurations with statistically planar flames [1,2,3]. The domain is cubic with dimensions of 40 δ s t × 40 δ s t × 40 δ s t , and the number of grid points is set at 512 × 512 × 512 . This ensures a minimum of 10 points to resolve the thermal flame thickness δ s t and at least 2 points for the Kolmogorov length scale η . Similar configurations have been explored in several previous studies [30,42].

4. Results

Despite the simulations being conducted for a predetermined duration after the start of turbulence, regardless of the simulation parameters, the duration varies due to the differing time required for each kernel to reach the target radius. Consequently, the concentration level of steam differs slightly between the cases upon the activation of turbulence. Nonetheless, steam concentrations remain low across all cases at the time of turbulence activation because the droplets reside outside the burned gas region, and the differences in steam concentration are minimal for this reason.
The behavior of the progress variable and non-dimensional temperature field are equal in adiabatic and equidiffusive flames (i.e., L e e f f = 1 ). However, our simulations encompass non-adiabatic conditions owing to droplet evaporation. Consequently, we expect a notable difference between the progress variable and the non-dimensional temperature, particularly in proximity to evaporating droplet locations due to their cooling and preferential diffusion effects.
Figure 1 shows iso-surfaces of the progress variable and non-dimensional temperature at c = 0.9 and T = 0.9 . The two iso-surfaces exhibit similar topological characteristics, albeit with differences that are mainly located in regions where droplet evaporation has occurred. The formation of small dimples is attributed to the cooling effect induced by the phase transition of water going from the liquid to the gaseous state. Notably, this cooling effect manifests more prominently in the temperature field, while the progress variable is only indirectly influenced by evaporation and the effects are less discernible. Moreover, the impact of a non-unity Lewis number is less relevant in comparison to non-adiabatic effects in this figure, especially because at stoichiometric conditions, preferential diffusion is less significant compared to leaner conditions, such as ϕ = 0.8 . An essential consequence of preferential diffusion is the occurrence of regions where the temperature surpasses the adiabatic flame temperature corresponding to the specific equivalence ratio. The examination of the non-dimensional temperature field underscores this characteristic trait of lean flames with L e H 2 < 1 .
Figure 2 illustrates the X Y mid-plane depicting the non-dimensional temperature θ = ( T ^ T ^ 0 ) / ( T ^ a d , ( ϕ ) T ^ 0 ) fields from the various simulated cases in this study at t = 1.5   t c h e m after the initiation of turbulence. It is immediately clear that the maximum temperature is much higher in the cases with ϕ = 0.8 and L e H 2 < 1 than in their counterparts without preferential diffusion. In the equidiffusive cases, θ is mostly around 1, while there is much more occurrence of θ > 1 with preferential diffusion, as is also confirmed later on in Figure 3. Furthermore, it is clear that the case with ϕ = 0.8 and without preferential diffusion exhibits the slowest flame propagation. Moreover, the temperature field in the case with ϕ = 0.8 and L e H 2 < 1 demonstrates significant distortion induced by turbulence. This pronounced wrinkling of the flame surface is attributed to the relatively stronger impact of preferential diffusion under lean conditions, which leads to an increased flame surface area, which is consistent with previous findings [20,43,44]. The importance of preferential diffusion as a function of the equivalence ratio can be evaluated through the effective Lewis number L e e f f . This parameter considers the overall behavior of the mixture concerning mass and temperature diffusion. The formulation suggested by Bechtold and Matalon for lean conditions is as follows [43]:
L e e f f = 1 + ( L e O 2 1 ) + ( L e H 2 1 ) A 1 + A ; A = 1 + β ( 1 ϕ 1 ) .
This definition comprises a weighted mean of the Lewis number of the primary reactants, contingent upon the equivalence ratio via the function A, which is influenced by the Zeldovich number β .
In Table 2, the effective Lewis numbers are reported for cases A-F for the two equivalence ratios examined in this study, revealing a discernible decrease in L e e f f as ϕ diminishes. Additionally, in Figure 2, a cooling effect stemming from droplet evaporation is observed, which is dependent on the maximum gas temperature and residence time of the droplets. Consequently, the most pronounced cooling effects are observed in the case with ϕ = 1 and L e H 2 = 1 , where the maximum temperature is only slightly less than the case with L e H 2 < 1 and the duration of the interaction of the droplets with the flame and the burned gasses is higher than that of the stoichiometric case with L e H 2 < 1 . Conversely, despite a significantly longer droplet–flame interaction time in the case with ϕ = 0.8 and no preferential diffusion, evaporation remains minimal due to the relatively lower burned gas temperature.
The probability density functions depicted in Figure 3 substantiate in a quantitative manner the earlier observations based on slices of the temperature field in Figure 2. Specifically, cases characterized by preferential diffusion under lean conditions manifest regions featuring super-adiabatic temperatures. Analysis of the probability density functions indicates that the cases with ϕ = 0.8 generally yield flames exhibiting higher occurrences of super-adiabatic temperatures than the stoichiometric cases when L e H 2 < 1 . This discrepancy arises due to the much lower effective Lewis number of the mixture, which causes more favorable conditions for this particular preferential diffusion effect. Moreover, the influence of turbulence-induced stretch enhances the consequences of preferential diffusion, which is particularly pronounced in lean scenarios compared to stoichiometric ones due to the favorable effects on topology, as observed by Concetti et al. [45]. Notably, temperature uniformity is enhanced in instances where preferential diffusion is absent. Probability density functions in the right column attain higher values and lack the extended tails observed in the results in the left column due to the lower thermal non-uniformity in cases with L e H 2 = 1 than in the cases with L e H 2 < 1 . Across nearly all cases, temperature distributions demonstrate only moderate sensitivity to the addition of water, even for the case with the smallest initial diameter. However, the stoichiometric case without preferential diffusion exhibits the most substantial change of the probability density function, as evidenced by the peak shifting to lower temperature values and a reduction in the peak height attributable to thermal non-uniformity arising from water droplet evaporation. Nonetheless, all droplet diameters fail to induce markedly distinct behavior in the temperature distributions, even in instances where evaporation is the most intense among all cases considered in this study. Examination of the Y W g fields offers a qualitative assessment of evaporation levels and their anticipated dilution effects on the flame.
Figure 4 illustrates the fields of gaseous water concentration resulting from evaporation across the various cases. The concentrations remain comparable across different initial droplet dimensions and are contingent solely upon variations in the equivalence ratio and flame propagation speed. The latter factor is particularly significant, as it influences the residence time of particles within the flame’s hot region (i.e., flame and burned gas region), as explained at the beginning of the present section. It can be observed that for the majority of cases, evaporation is insufficient to substantially impact temperature fields in a significant manner. The sole exception is noted in the case with ϕ = 1 and the unity Lewis number. Nevertheless, minor cooling effects, which are also discernible in Figure 3, may still influence flame speed and propagation characteristics.
It is well known that the flame speed is influenced by stretch effects (=strain+curvature effects) via the Markstein relation [46]. Hence, all propagation-related properties are now plotted as functions of the average radius r V . This approach aims to isolate the dependence of the propagation speed on the curvature similarly across all analyzed cases. The average radius is computed from the volume using the following formula:
r V = 3 4 π V c d V 1 / 3
The normalized total burning rate S relative to the time of turbulence initiation is defined as follows:
S = V ω ˙ c d V ρ 0 S L , ( ϕ ) A 0 = V ω ˙ c d V ρ 0 S L , ( ϕ ) 4 π ( r V 0 ) 2
In this way, S can be simply interpreted as the relative evolution of the total burning rate. In the definition of S, the stretched laminar flame speed S L , ( ϕ ) at the target radius for turbulence initiation and the initial projected flame area A 0 , computed through the initial radius based on volume r V 0 , are used for the estimation of the initial total burning rate.
Figure 5 depicts the relative evolution of the total burning rate for cases with and without preferential diffusion, which are presented on the left and right sides, respectively. The presence of preferential diffusion notably enhances the growth of this quantity, particularly under lean conditions. The higher values of S in lean conditions when L e H 2 < 1 can be elucidated from the effective Lewis number values presented in Table 2. At the same time, the lean case without preferential diffusion exhibits a notably slower evolution compared to all other cases, which also justifies the limited span of data along the abscissa in the plots. Furthermore, slight differences can be observed between the trends associated with evaporation when comparing cases with and without water addition. However, these deviations remain marginal due to the modest degree of evaporation witnessed in the simulated cases. This is attributed to the limited residence time of droplets within the flame and the limited time spent in the burned gas region. The former directly influences the flame structure, while the latter primarily has indirect effects related to the decrease in thermal expansion. Nonetheless, despite the observed effects being limited, discernible physical trends and phenomena are already apparent. The evolution of the total burning rate is influenced by variations in reactivity, which are reflected in the burning rate per unit area Ω = ( S T / S L , ( ϕ ) 0 ) / ( A c / A ) = V ω c ˙ d V / ( ρ 0 S L , ( ϕ ) 0 A c ) , where S L , ( ϕ ) 0 is the unstretched laminar burning velocity, A is the projected flame area at the instant of the evaluation, the turbulent flame area A c = V | c | d V , and the turbulent flame speed S T = V ω c ˙ d V / ( ρ 0 A ) . In line with the findings of several previous studies [23,44], instances characterized by pronounced preferential diffusion induced by L e e f f < 1 exhibit significantly greater values of both Ω and A c compared to cases where this phenomenon is absent. Consequently, these characteristics are now observed in order to discern their responses to a non-zero curvature and the introduction of liquid water.
According to Damköhler’s first hypothesis, the value of Ω , as it is defined in the present work, is expected to equal 1.0 under adiabatic and equidiffusive premixed statistically planar flame conditions [47,48]. Statistically planar turbulent flames with a Lewis number smaller (larger) than one are characterized by a value of Ω > 1 ( Ω < 1 ) [49]. Unity Lewis number flames with a globally negative mean curvature (e.g., Bunsen flames) are characterized by a value of Ω > 1 [48], whereas flames with a globally positive mean curvature (e.g., flame kernels) are characterized by a value of Ω < 1 [50]. For curved flames with a non-unity Lewis number, both effects compete with each other [50]. In addition, as previously noted, the conditions under investigation in this study are non-adiabatic due to water evaporation, and notably, the effects of preferential diffusion are also considered. The values of burning rate per unit area, depicted in Figure 6, illustrate that when L e H 2 < 1 , the quantity Ω > 1 , reaching values up to Ω 1.5 . The cases with L e H 2 = 1 settle to values marginally larger than unity for ϕ = 0.8 and smaller than unity for ϕ = 1.0 . In order to explain these effects, we introduce the Markstein relation:
S d = S L , ( ϕ ) 0 ( 1 L K / S L , ( ϕ ) 0 ) ; K = 1 A c d A c d t = 2 S d κ m + a t
The displacement speed S d is the velocity with which the flame propagates with respect to an initial coincident material surface [51]. In this context, L denotes the Markstein length, which is typically positive in the absence of thermodiffusive instability, and K signifies the stretch rate, which describes the variation of flame area over time. The stretch rate is a function of flame curvature, denoted as κ m = 0.5 · n , and tangential strain rate, represented by a t = n n : u + · u . In the definitions of curvature and tangential strain rate, n signifies the normal direction of the flame front and is defined as n = c / | c | . The displacement speed, defined in Equation (13), is a local quantity dependent on the reaction rate, diffusive fluxes, concentration gradients, and evaporation effects, where the latter two are generally negligible compared to the reaction rate and diffusive fluxes [52].
ρ S d | c | ω ˙ c + · ( ρ D c ) .
Upon integration over the computational volume, the contribution of the diffusive flux vanishes [48], leading to the following equation:
V ρ S d | c | d V = V ω ˙ c d V .
The volume integral of the reaction rate can also be expressed as Ω ρ 0 S L , ( ϕ ) 0 A c = V ω ˙ c d V . Using the surface averaging introduced by Boger et al. [53], Q s = Q | c | / | c | , where Q = Q d V / V denotes the volume average of a general quantity Q, the left-hand side of Equation (15) can be rewritten as:
Ω ρ 0 S L , ( ϕ ) 0 A c = ρ S d s A c
By substituting the density-weighted Markstein relation (Equation (13)) under the assumption of a constant Markstein length, the following expression is obtained, indicating that the sign of ( Ω 1 ) depends on the product of the sign of the Markstein length and the sign of the surface-averaged stretch rate:
Ω = 1 L K s / S L , ( ϕ ) 0 .
According to the definition of the Markstein length associated with the displacement speed given in [18], L attains the values L = 0.31 , 0.14 for ϕ = 0.8 , 1.0 , respectively, for L e H 2 < 1 , whereas it remains larger than unity for L e H 2 = 1 , i.e., L = 0.34 , 0.29 for ϕ = 0.8 , 1.0 , respectively. This explains the qualitative trends observed in Figure 6. A quantitative assessment is complicated by the fact that Equation (17) is only valid for small stretch rates; the Markstein length is not a constant (see, e.g., [51]) and is very sensitive to the definition of the flame speed [18].
The influence of water on the normalized burning rate per unit area is nearly negligible, as evaporation remains minimal. However, as demonstrated in prior works by the authors [1,3], significant cooling effects of water under substantial evaporation rates can considerably impact Ω , leading to values much lower than unity for both hydrogen and hydrocarbon combustion. Another factor for investigating the influence on flame propagation characteristics is the flame area. Figure 7 illustrates the temporal evolution of the flame area for cases with preferential diffusion on the left and without preferential diffusion on the right. The trends of the normalized flame area are strongly correlated to those of the normalized total burning rate profiles, highlighting its significant influence on propagation speed.
Additionally, water evaporation marginally impacts the trend of flame area evolution, resulting in diminished area generation due to the cooling effect, which mitigates flame wrinkling induced by turbulence. Another noteworthy characteristic is observed in cases where preferential diffusion is active, where the slope of the curves grows faster with the radius than in the cases without preferential diffusion. In a spherical expanding flame configuration, area increase is influenced by two factors: flame outward expansion, evident in the increase in radius based on volume, and flame surface wrinkling. The intensity of the latter can be estimated by the ratio between the radius estimated from the volume and the radius estimated from the flame area ( r A = [ A c / ( 4 π ) ] 1 / 2 ); when this ratio substantially exceeds unity, it suggests intense surface wrinkling.
Figure 8 illustrates the deviation of this ratio from unity. As anticipated, cases without preferential diffusion generally exhibit slower growth compared to those with preferential diffusion. However, it is noteworthy that the ratio r A / r V is higher for lean cases than for stoichiometric cases when L e H 2 = 1 . This is due to the data being plotted over the radial evolution rather than time, with the duration of flame–turbulence interaction being the determining factor for flame wrinkling. Water influences the ratio r A / r V in a stronger way in cases with preferential diffusion than in cases with L e H 2 = 1 .
Initial flame kernels are prone to quenching caused by curvature stretch effects and turbulence. Water injection has the potential to further reduce the reaction rate, as found by earlier studies [1,45]. Nevertheless, the present results suggest that due to the low evaporation rate of water, water droplets are unlikely to affect the initial flame kernel growth significantly and have a comparably small effect on the conditions for sustained combustion. Nevertheless, for longer interaction times, water droplets are likely to have the potential to reduce flame acceleration or turbulent flame speed for both the preferential diffusion and equidiffusive cases, as has been observed in previous works [1,45].

5. Conclusions

In this study, direct numerical simulations of turbulent spherical expanding premixed flames are conducted under varying equivalence ratios, water loading, and initial droplet diameters. The simulations are conducted with and without preferential fuel diffusion to elucidate the effects arising from this characteristic. The primary objective of this study is to identify qualitative trends in the interaction between water, flame, turbulence, and preferential diffusion, where the level of detail obtained through direct numerical simulation goes beyond that obtained by existing experimental results. Initial observations reveal that under non-adiabatic conditions and in the presence of preferential diffusion, the temperature and progress variable fields exhibit subtle differences. Notably, cases with preferential diffusion manifest regions with temperatures surpassing the adiabatic temperature under the prescribed conditions. In this context, the burned gas temperature of the system, the residence time within the flame, and the lifetime in the burned gas region of droplets are crucial parameters for determining the evaporation rate when the initial droplet diameters are identical. However, in the configuration considered in this study, the evaporation rate of water remains moderate, and the differences in propagation characteristics are not significantly altered between cases with and without water addition. Nevertheless, a reduction in total burning rate, flame area, and burning rate per unit area is observed.
In the presence of preferential diffusion, the total burning rate increases with radius much more rapidly compared to the thermodiffusively neutral (i.e., unity Lewis number) cases. The burning rate per unit area exhibits a different response to flame stretching depending on whether the condition of a non-unity Lewis number is considered or not. In the former case, the burning rate increases with time and flame stretching, while in the latter case, a slight decrease is observed. This phenomenon is linked to the differing sign of the Markstein length. Additionally, the trend of the flame area varies between cases with and without preferential diffusion, primarily due to flame wrinkling, which can be estimated by the ratio of the radius computed with the flame area to the flame volume. The maximum value of this quantity is attained for cases with preferential diffusion and an equivalence ratio of 0.8 .
While the present results suggest a weak influence of water droplets on the initial growth of a spherical flame, further analyses will be necessary so that the current findings can be validated and more comprehensive physical insights can be extracted in the presence of detailed chemical and transport mechanisms. Alternatively, it might be of interest to study the long-term interaction of the flame and the water, which will require considerably larger computational domains and might limit the analysis to two-dimensional cases.

Author Contributions

Conceptualization, J.H. and N.C.; methodology, J.H. and N.C.; software, J.H. and N.C.; formal analysis, R.C. and J.H.; investigation, R.C., J.H. and M.K.; resources, M.K.; data curation, R.C.; writing—original draft preparation, R.C., J.H. and M.K.; writing—review and editing, R.C., J.H., M.K. and N.C.; visualization, R.C.; supervision, J.H., M.K. and N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by dtec.bw—Digitalization and Technology Research Center of the Bundeswehr; dtec.bw is funded by the European Union—NextGenerationEU.

Data Availability Statement

The data used for the present analysis are available from the corresponding author under reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hasslberger, J.; Ozel-Erol, G.; Chakraborty, N.; Klein, M.; Cant, S. Physical effects of water droplets interacting with turbulent premixed flames: A Direct Numerical Simulation analysis. Combust. Flame 2021, 229, 111404. [Google Scholar] [CrossRef]
  2. Hasslberger, J.; Concetti, R.; Chakraborty, N.; Klein, M. Inertial effects on the interaction of water droplets with turbulent premixed flames: A direct numerical simulation analysis. Proc. Combust. Inst. 2023, 39, 2575–2586. [Google Scholar] [CrossRef]
  3. Concetti, R.; Hasslberger, J.; Chakraborty, N.; Klein, M. Analysis of Water Droplet Interaction with Turbulent Premixed and Spray Flames Using Carrier Phase Direct Numerical Simulations. Combust. Sci. Technol. 2023, 195, 1411–1433. [Google Scholar] [CrossRef]
  4. Zeldovich, Y.B. The oxidation of nitrogen in combustion and explosions. J. Acta Physicochim. 1946, 21, 577. [Google Scholar]
  5. Tesfa, B.; Mishra, R.; Gu, F.; Ball, A. Water injection effects on the performance and emission characteristics of a CI engine operating with biodiesel. Renew. Energy 2012, 37, 333–344. [Google Scholar] [CrossRef]
  6. Mingrui, W.; Sa, N.T.; Turkson, R.F.; Jinping, L.; Guanlun, G. Water injection for higher engine performance and lower emissions. J. Energy Inst. 2017, 90, 285–299. [Google Scholar] [CrossRef]
  7. Lellek, S. Pollutant Formation in Premixed Natural Gas Swirl Flames with Water Injection. Ph.D. Thesis, Technische Universität München, Munich, Germany, 2017. [Google Scholar]
  8. Thomas, G. On the conditions required for explosion mitigation by water sprays. Process. Saf. Environ. Prot. 2000, 78, 339–354. [Google Scholar] [CrossRef]
  9. Kim, S.C.; Ryou, H.S. An experimental and numerical study on fire suppression using a water mist in an enclosure. Build. Environ. 2003, 38, 1309–1316. [Google Scholar] [CrossRef]
  10. Ohi, J. Hydrogen energy cycle: An overview. J. Mater. Res. 2005, 20, 3180–3187. [Google Scholar] [CrossRef]
  11. Ahmed, M.R.; Barua, T.; Das, B.K. A comprehensive review on techno-environmental analysis of state-of-the-art production and storage of hydrogen energy: Challenges and way forward. Energy Sources Part Recover. Util. Environ. Eff. 2023, 45, 5905–5937. [Google Scholar] [CrossRef]
  12. Oliveira, A.M.; Beswick, R.R.; Yan, Y. A green hydrogen economy for a renewable energy society. Curr. Opin. Chem. Eng. 2021, 33, 100701. [Google Scholar] [CrossRef]
  13. Falcone, P.M.; Hiete, M.; Sapio, A. Hydrogen economy and sustainable development goals: Review and policy insights. Curr. Opin. Green Sustain. Chem. 2021, 31, 100506. [Google Scholar] [CrossRef]
  14. Matalon, M.; Cui, C.; Bechtold, J. Hydrodynamic theory of premixed flames: Effects of stoichiometry, variable transport coefficients and arbitrary reaction orders. J. Fluid Mech. 2003, 487, 179–210. [Google Scholar] [CrossRef]
  15. Keppeler, R.; Pfitzner, M. Modelling of Landau–Darrieus and thermo-diffusive instability effects for CFD simulations of laminar and turbulent premixed combustion. Combust. Theory Model. 2015, 19, 1–28. [Google Scholar] [CrossRef]
  16. Frankel, M.; Sivashinsky, G. The effect of viscosity on hydrodynamic stability of a plane flame front. Combust. Sci. Technol. 1982, 29, 207–224. [Google Scholar] [CrossRef]
  17. Class, A.G.; Matkowsky, B.; Klimenko, A. Stability of planar flames as gasdynamic discontinuities. J. Fluid Mech. 2003, 491, 51–63. [Google Scholar] [CrossRef]
  18. Lipatnikov, A.; Chomiak, J. Molecular transport effects on turbulent flame propagation and structure. Prog. Energy Combust. Sci. 2005, 31, 1–73. [Google Scholar] [CrossRef]
  19. Altantzis, C.; Frouzakis, C.; Tomboulides, A.; Matalon, M.; Boulouchos, K. Hydrodynamic and thermodiffusive instability effects on the evolution of laminar planar lean premixed hydrogen flames. J. Fluid Mech. 2012, 700, 329–361. [Google Scholar] [CrossRef]
  20. Berger, L.; Kleinheinz, K.; Attili, A.; Pitsch, H. Characteristic patterns of thermodiffusively unstable premixed lean hydrogen flames. Proc. Combust. Inst. 2019, 37, 1879–1886. [Google Scholar] [CrossRef]
  21. Creta, F.; Lapenna, P.E.; Lamioni, R.; Fogla, N.; Matalon, M. Propagation of premixed flames in the presence of Darrieus–Landau and thermal diffusive instabilities. Combust. Flame 2020, 216, 256–270. [Google Scholar] [CrossRef]
  22. Rasool, R.; Chakraborty, N.; Klein, M. Effect of non-ambient pressure conditions and Lewis number variation on direct numerical simulation of turbulent Bunsen flames at low turbulence intensity. Combust. Flame 2021, 231, 111500. [Google Scholar] [CrossRef]
  23. Berger, L.; Grinberg, M.; Jürgens, B.; Lapenna, P.E.; Creta, F.; Attili, A.; Pitsch, H. Flame fingers and interactions of hydrodynamic and thermodiffusive instabilities in laminar lean hydrogen flames. Proc. Combust. Inst. 2023, 39, 1525–1534. [Google Scholar] [CrossRef]
  24. Rutland, C.; Trouvé, A. Direct simulations of premixed turbulent flames with nonunity Lewis numbers. Combust. Flame 1993, 94, 41–57. [Google Scholar] [CrossRef]
  25. Klein, M.; Chakraborty, N.; Jenkins, K.W.; Cant, R.S. Effects of initial radius on the propagation of premixed flame kernels in a turbulent environment. Phys. Fluids 2006, 18, 055102. [Google Scholar] [CrossRef]
  26. Klein, M.; Chakraborty, N.; Cant, R. Effects of turbulence on self-sustained combustion in premixed flame kernels: A direct numerical simulation (DNS) study. Flow Turbul. Combust. 2008, 81, 583–607. [Google Scholar] [CrossRef]
  27. Ashurst, W.T. A simple illustration of turbulent flame ball growth. Combust. Sci. Technol. 1995, 104, 19–32. [Google Scholar] [CrossRef]
  28. Jenkins, K.; Cant, R. Curvature effects on flame kernels in a turbulent environment. Proc. Combust. Inst. 2002, 29, 2023–2029. [Google Scholar] [CrossRef]
  29. Jomaas, G.; Law, C.K.; Bechtold, J. On transition to cellularity in expanding spherical flames. J. Fluid Mech. 2007, 583, 1–26. [Google Scholar] [CrossRef]
  30. Ozel Erol, G.; Hasslberger, J.; Klein, M.; Chakraborty, N. A direct numerical simulation investigation of spherically expanding flames propagating in fuel droplet-mists for different droplet diameters and overall equivalence ratios. Combust. Sci. Technol. 2019, 191, 833–867. [Google Scholar] [CrossRef]
  31. Reveillon, J.; Vervisch, L. Analysis of weakly turbulent dilute-spray flames and spray combustion regimes. J. Fluid Mech. 2005, 537, 317–347. [Google Scholar] [CrossRef]
  32. Balachandar, S.; Eaton, J.K. Turbulent dispersed multiphase flow. Annu. Rev. Fluid Mech. 2010, 42, 111–133. [Google Scholar] [CrossRef]
  33. Elghobashi, S. Direct numerical simulation of turbulent flows laden with droplets or bubbles. Annu. Rev. Fluid Mech. 2019, 51, 217–244. [Google Scholar] [CrossRef]
  34. Haruki, Y.; Pillai, A.L.; Kitano, T.; Kurose, R. Numerical investigation of flame propagation in fuel droplet arrays. At. Sprays 2018, 28, 357–388. [Google Scholar] [CrossRef]
  35. De Chaisemartin, S.; Fréret, L.; Kah, D.; Laurent, F.; Fox, R.O.; Reveillon, J.; Massot, M. Eulerian models for turbulent spray combustion with polydispersity and droplet crossing. C. R. Mécanique 2009, 337, 438–448. [Google Scholar] [CrossRef]
  36. Goodwin, D.G.; Moffat, H.K.; Schoegl, I.; Speth, R.L.; Weber, B.W. Cantera: An Object-Oriented Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport Processes. Version 3.0.0. 2023. Available online: https://zenodo.org/records/8137090 (accessed on 15 August 2024).
  37. Göktolga, M.U.; van Oijen, J.A.; de Goey, L.P.H. 3D DNS of MILD combustion: A detailed analysis of heat loss effects, preferential diffusion, and flame formation mechanisms. Fuel 2015, 159, 784–795. [Google Scholar] [CrossRef]
  38. Xin, Y.X.; Yoo, C.S.; Chen, J.H.; Law, C.K. A DNS study of self-accelerating cylindrical hydrogen–air flames with detailed chemistry. Proc. Combust. Inst. 2015, 35, 753–760. [Google Scholar] [CrossRef]
  39. Fernández-Tarrazo, E.; Sánchez, A.L.; Liñán, A.; Williams, F.A. A simple one-step chemistry model for partially premixed hydrocarbon combustion. Combust. Flame 2006, 147, 32–38. [Google Scholar] [CrossRef]
  40. Abo-Amsha, K.; Chakraborty, N. Surface Density Function and Its Evolution in Homogeneous and Inhomogeneous Mixture n-Heptane MILD Combustion. Combust. Sci. Technol. 2023, 195, 1483–1508. [Google Scholar] [CrossRef]
  41. Wray, A.A. Minimal Storage Time Advancement Schemes for Spectral Methods; Report No. MS; NASA Ames Research Center: Moffett Field, CA, USA, 1990; Volume 202. [Google Scholar]
  42. Ozel Erol, G.; Hasslberger, J.; Klein, M.; Chakraborty, N. A direct numerical simulation analysis of spherically expanding turbulent flames in fuel droplet-mists for an overall equivalence ratio of unity. Phys. Fluids 2018, 30, 086104. [Google Scholar] [CrossRef]
  43. Bechtold, J.; Matalon, M. The dependence of the Markstein length on stoichiometry. Combust. Flame 2001, 127, 1906–1913. [Google Scholar] [CrossRef]
  44. Chakraborty, N.; Klein, M. Influence of Lewis Number on curvature effects in turbulent premixed flame propagation in the thin reaction zones regime. Phys. Fluids 2005, 17, 105105. [Google Scholar] [CrossRef]
  45. Concetti, R.; Hasslberger, J.; Chakraborty, N.; Klein, M. Effects of liquid water addition on turbulent premixed hydrogen/air combustion. Fuel 2024, 373, 132314. [Google Scholar] [CrossRef]
  46. Clavin, P.; Joulin, G. Premixed flames in large scale and high intensity turbulent flow. J. Phys. Lett. 1983, 44, 1–12. [Google Scholar] [CrossRef]
  47. Damköhler, G. Der einfluss der turbulenz auf die flammengeschwindigkeit in gasgemischen. Zeitschrift Elektrochem. Angew. Phys. Chem. 1940, 46, 601–626. [Google Scholar] [CrossRef]
  48. Chakraborty, N.; Alwazzan, D.; Klein, M.; Cant, R.S. On the validity of Damköhler’s first hypothesis in turbulent Bunsen burner flames: A computational analysis. Proc. Combust. Inst. 2019, 37, 2231–2239. [Google Scholar] [CrossRef]
  49. Chakraborty, N.; Wang, L.; Klein, M. Streamline segment statistics of premixed flames with nonunity Lewis numbers. Phys. Rev. E 2014, 89, 033015. [Google Scholar] [CrossRef]
  50. Ozel-Erol, G.; Klein, M.; Chakraborty, N. Lewis Number Effects on Flame Speed Statistics in Spherical Turbulent Premixed Flames. Flow Turbul. Combust. 2021, 106, 1043–1063. [Google Scholar] [CrossRef]
  51. Herbert, A.; Ahmed, U.; Chakraborty, N.; Klein, M. Applicability of extrapolation relations for curvature and stretch rate dependences of displacement speed for statistically planar turbulent premixed flames. Combust. Theory Model. 2020, 24, 1021–1038. [Google Scholar] [CrossRef]
  52. Ozel-Erol, G.; Hasslberger, J.; Chakraborty, N.; Klein, M. Effects of water droplet injection on turbulent premixed flame propagation: A direct numerical simulation investigation. Flow Turbul. Combust. 2023, 110, 105–123. [Google Scholar] [CrossRef]
  53. Boger, M.; Veynante, D.; Boughanem, H.; Trouvé, A. Direct numerical simulation analysis of flame surface density concept for large eddy simulation of turbulent premixed combustion. In Symposium (International) on Combustion; Elsevier: Amsterdam, The Netherlands, 1998; Volume 27, pp. 917–925. [Google Scholar]
Figure 1. Iso-surfaces of progress variable ( c = 0.9 ) on the left and non-dimensional temperature ( T = 0.9 ) on the right. The data are taken for the case with ϕ = 1 and L e H 2 < 1 with water droplet initial diameter a d / δ s t = 0.04 at t = 1.5 t c h e m after superimposing the turbulence.
Figure 1. Iso-surfaces of progress variable ( c = 0.9 ) on the left and non-dimensional temperature ( T = 0.9 ) on the right. The data are taken for the case with ϕ = 1 and L e H 2 < 1 with water droplet initial diameter a d / δ s t = 0.04 at t = 1.5 t c h e m after superimposing the turbulence.
Energies 17 04632 g001
Figure 2. X Y mid-plane of the non-dimensional temperature θ = ( T ^ T ^ 0 ) / ( T ^ a d , ( ϕ ) T ^ 0 ) fields at t = 1.5 t c h e m after superimposing the turbulence for the cases with ϕ = 1 (1st and 3rd rows) and ϕ = 0.8 (2nd and 4th rows) without droplets (left column), with droplets of initial diameter a d / δ s t = 0.04 (central column), and with droplets of initial diameter a d / δ s t = 0.03 (right column). The cases with L e H 2 < 1 are in the first two rows, while the cases in which L e H 2 = 1 are in the last two rows. The white iso-lines refer to different progress variable values ( c = 0.1 , 0.5 , and 0.9 ).
Figure 2. X Y mid-plane of the non-dimensional temperature θ = ( T ^ T ^ 0 ) / ( T ^ a d , ( ϕ ) T ^ 0 ) fields at t = 1.5 t c h e m after superimposing the turbulence for the cases with ϕ = 1 (1st and 3rd rows) and ϕ = 0.8 (2nd and 4th rows) without droplets (left column), with droplets of initial diameter a d / δ s t = 0.04 (central column), and with droplets of initial diameter a d / δ s t = 0.03 (right column). The cases with L e H 2 < 1 are in the first two rows, while the cases in which L e H 2 = 1 are in the last two rows. The white iso-lines refer to different progress variable values ( c = 0.1 , 0.5 , and 0.9 ).
Energies 17 04632 g002
Figure 3. Probability density function of the non-dimensional temperature θ within the reaction-dominated region of the flame c = 0.7 0.99 . The data are collected at t = 1.5 t c h e m after superimposing the turbulence. The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Figure 3. Probability density function of the non-dimensional temperature θ within the reaction-dominated region of the flame c = 0.7 0.99 . The data are collected at t = 1.5 t c h e m after superimposing the turbulence. The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Energies 17 04632 g003
Figure 4. X Y mid-plane of the gaseous water concentration coming only from the evaporation Y W g field at t = 1.5 t c h e m after superimposing the turbulence for the cases with ϕ = 1 (1st and 3rd rows) and ϕ = 0.8 (2nd and 4th rows) with droplets with initial diameter a d / δ s t = 0.04 (left column) and with droplets with initial diameter a d / δ s t = 0.03 (right column). The cases with L e H 2 < 1 are in the first two rows, while the cases in which L e H 2 = 1 are in the last two rows. The white iso-lines refer to different progress variable values ( c = 0.1 , 0.5 , and 0.9 ).
Figure 4. X Y mid-plane of the gaseous water concentration coming only from the evaporation Y W g field at t = 1.5 t c h e m after superimposing the turbulence for the cases with ϕ = 1 (1st and 3rd rows) and ϕ = 0.8 (2nd and 4th rows) with droplets with initial diameter a d / δ s t = 0.04 (left column) and with droplets with initial diameter a d / δ s t = 0.03 (right column). The cases with L e H 2 < 1 are in the first two rows, while the cases in which L e H 2 = 1 are in the last two rows. The white iso-lines refer to different progress variable values ( c = 0.1 , 0.5 , and 0.9 ).
Energies 17 04632 g004
Figure 5. Evolution of the normalized total burning rate S versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Figure 5. Evolution of the normalized total burning rate S versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Energies 17 04632 g005
Figure 6. Burning rate per unit area normalized by the unstretched laminar equivalent Ω versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Figure 6. Burning rate per unit area normalized by the unstretched laminar equivalent Ω versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Energies 17 04632 g006
Figure 7. Flame area normalized by the flame area when the turbulence is superimposed A c / A 0 versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right it is L e H 2 = 1 .
Figure 7. Flame area normalized by the flame area when the turbulence is superimposed A c / A 0 versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right it is L e H 2 = 1 .
Energies 17 04632 g007
Figure 8. Flame-area-based radius normalized by volume-based radius r A / r V versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Figure 8. Flame-area-based radius normalized by volume-based radius r A / r V versus the radius computed from the flame kernel volume r V . The blue continuous lines and markers refer to the stoichiometric cases, while the red continuous lines and markers refer to the cases with ϕ = 0.8 . The continuous lines represent the cases without water addition, and the markers refer to the cases with water addition. The different markers indicate the initial diameter of the droplets. On the left, the fuel Lewis number is L e H 2 < 1 , while on the right, L e H 2 = 1 .
Energies 17 04632 g008
Table 1. Overview of the parametric study in terms of equivalence ratio ϕ , overall water loading Y W o v = Y W g + Y W l , normalized initial droplet diameter a d / δ s t , and fuel Lewis number L e H 2 .
Table 1. Overview of the parametric study in terms of equivalence ratio ϕ , overall water loading Y W o v = Y W g + Y W l , normalized initial droplet diameter a d / δ s t , and fuel Lewis number L e H 2 .
Case ϕ Y W ov a d / δ st Le H 2
A 1.0 0.0 0.2851
B 1.0 0.1 0.04 0.2851
C 1.0 0.1 0.03 0.2851
D 0.8 0.0 0.2851
E 0.8 0.1 0.04 0.2851
F 0.8 0.1 0.03 0.2851
G 1.0 0.0 1.0
H 1.0 0.1 0.04 1.0
I 1.0 0.1 0.03 1.0
J 0.8 0.0 1.0
K 0.8 0.1 0.04 1.0
L 0.8 0.1 0.03 1.0
Table 2. Effective Lewis numbers L e e f f at different equivalence ratios ϕ for the present thermochemistry (cases A–F).
Table 2. Effective Lewis numbers L e e f f at different equivalence ratios ϕ for the present thermochemistry (cases A–F).
ϕ Le eff
0.8 0.51
1.0 0.64
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Concetti, R.; Hasslberger, J.; Chakraborty, N.; Klein, M. Effects of Water Mist on the Initial Evolution of Turbulent Premixed Hydrogen/Air Flame Kernels. Energies 2024, 17, 4632. https://doi.org/10.3390/en17184632

AMA Style

Concetti R, Hasslberger J, Chakraborty N, Klein M. Effects of Water Mist on the Initial Evolution of Turbulent Premixed Hydrogen/Air Flame Kernels. Energies. 2024; 17(18):4632. https://doi.org/10.3390/en17184632

Chicago/Turabian Style

Concetti, Riccardo, Josef Hasslberger, Nilanjan Chakraborty, and Markus Klein. 2024. "Effects of Water Mist on the Initial Evolution of Turbulent Premixed Hydrogen/Air Flame Kernels" Energies 17, no. 18: 4632. https://doi.org/10.3390/en17184632

APA Style

Concetti, R., Hasslberger, J., Chakraborty, N., & Klein, M. (2024). Effects of Water Mist on the Initial Evolution of Turbulent Premixed Hydrogen/Air Flame Kernels. Energies, 17(18), 4632. https://doi.org/10.3390/en17184632

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