Abstract
The effects of turbulent Reynolds number on the statistical behaviour of the displacement speed have been studied using three-dimensional Direct Numerical Simulation of statistically planar turbulent premixed flames. The probability of finding negative values of the displacement speed is found to increase with increasing turbulent Reynolds number when the Damköhler number is held constant. It has been shown that the statistical behaviour of the Surface Density Function, and its strain rate and curvature dependence, plays a key role in determining the response of the different components of displacement speed. Increasing the turbulent Reynolds number is shown to reduce the strength of the correlations between tangential strain rate and dilatation rate with curvature, although the qualitative nature of the correlations remains unaffected. The dependence of displacement speed on strain rate and curvature is found to weaken with increasing turbulent Reynolds number when either Damköhler or Karlovitz number is held constant, but the qualitative nature of the correlation remains unaltered. The implications of turbulent Reynolds number effects in the context of Flame Surface Density (FSD) modelling have also been addressed, with emphasis on the influence of displacement speed on the curvature and propagation terms in the FSD balance equation.
1. Introduction
Displacement speed is an important quantity for the modelling of turbulent premixed flames using both level-set [1] and flame surface density (FSD) [2] approaches. The displacement speed is defined as the speed at which the flame surface moves normal to itself with respect to an initially coincident material surface. Statistics of displacement speed for turbulent premixed flames have been studied extensively in previous work based on two-dimensional direct numerical simulations (DNS) with detailed chemistry [3–13] and on three-dimensional DNS with simplified chemistry [14–22]. These studies have looked at the effects of local strain rate and curvature effect on [13–24] and have included the influence of root-mean-square (rms) turbulent velocity fluctuation magnitude and Lewis number Le [9, 15, 17, 22] on these effects. Peters [1] showed using an order of magnitude analysis that the curvature dependence of displacement speed becomes the leading order contributor in the thin reaction zones regime, whereas the curvature dependence remains relatively weak in the corrugated flamelets regime. This observation was supported using DNS data by Chakraborty [21] who also studied the individual response of the reaction, normal displacement, and curvature components of displacement speed. However, the effects of turbulent Reynolds number on displacement speed have yet to be addressed in detail. Hence, in this paper, the influence of on the statistics of has been analysed based on a three-dimensional compressible DNS database of statistically planar turbulent premixed flames. The main objectives of this study are as follows:(1)to demonstrate the influence of on the statistical distribution of (2)to demonstrate the influence of on the local curvature and strain rate dependencies of (3)to indicate the implications of the dependence of on in the context of flame surface density-based turbulent combustion modelling.
In order to address the above objectives, a series of DNS runs has been carried out for turbulent Reynolds numbers ranging from 20 to 100. In one group of simulations, the change in is obtained by modifying the Karlovitz number for a given value of Damköhler number , whereas in the other group of simulations the change in turbulent Reynolds number is induced by modifying the Damköhler number for a given value of Karlovitz number.
The rest of the paper is organised as follows. The mathematical background and the numerical details relevant to this work are discussed in the next two sections of this paper. Following this, results will be presented and subsequently discussed. In the final section, the main findings are summarised and conclusions are drawn.
2. Mathematical Background
Direct numerical simulations (DNS) of turbulent reacting flows ideally should be carried out in three-dimensions with detailed chemistry. However, such simulations remain extremely expensive and are not yet feasible for parametric studies. Previous work has indicated that displacement speed statistics obtained from three-dimensional DNS with simplified chemistry [14–22] are in agreement with those obtained from two-dimensional DNS with complex chemistry [3–13]. In the present study, three-dimensional DNS with simplified chemistry is adopted in order to allow for an extensive parametric study without excessive computational cost. The chemical reaction is represented in terms of single step Arrhenius type chemistry, and the species field is uniquely represented by a reaction progress variable, which may be defined based on a suitable product mass fraction according to where subscripts 0 and are used to refer to quantities in unburned reactants and fully burned products, respectively. The transport equation for reaction progress variable is given by where is the chemical reaction rate, is the fluid density, and is the diffusivity of reaction progress variable. The transport equation for reaction progress variable can be written in kinematic form for a given isosurface defined by [6, 7]
Comparing (2a) and (2b), the expression for the displacement speed of a given isosurface can be obtained as [1, 3–7] It is often useful to decompose the displacement speed into three separate components [3–7] where , , and are the reaction, normal diffusion, and tangential diffusion components which are given as where is the arithmetic mean of the two principal curvatures, which is defined as According to (6), flame surface that is curved convex to the unburned reactants has a positive curvature and vice versa. The definition of displacement speed and its components in (4)–(6) suggests that the statistics of the surface density function () [15, 25–27] and curvature are each likely to have a significant effect on the statistical behaviour of displacement speed. The transport equation for the SDF for an isosurface of is given by [15, 25–27] where is the tangential strain rate given by
It is evident from (7) and (8) that tangential strain rate is likely to affect the statistical behaviour of , which in turn can affect the local curvature and strain rate dependence of and . It is worth noting that in some previous studies [3–12], the statistics of the density-weighted displacement speed were analysed instead of the statistics of itself. For low Mach number unity Lewis number flames, and are related by for a given isosurface of c, and thus the statistics of are similar to those of . For this reason, and also since appears explicitly in the SDF and FSD transport equations [2, 15, 25–27], the statistics of are addressed in the present work. The local strain rate and curvature dependence of and its components for different values of turbulent Reynolds number will be discussed in Section 4 of this paper.
It is essential for the purposes of the present study to note the relationship between the turbulent Reynolds number and both the Damköhler and Karlovitz numbers. The turbulent Reynolds number is defined as , where is the density and is the dynamic viscosity of the unburned reactants, and is the integral length scale of the turbulence. For unity Lewis number flames, the turbulent Reynolds number can be scaled as [1] where is the Damköhler number and is the Karlovitz number. Here, is the thermal flame thickness with , and being the adiabatic flame temperature, unburned gas temperature, and instantaneous dimensional temperature, respectively, while the subscript “” refers to the unstrained planar laminar flame quantities. Equation (9) indicates that a change of either or will lead to a change in the turbulent Reynolds number . Thus, the influence of on the statistics of is essentially induced by changes in and/or .
3. Numerical Implementation
Three-dimensional compressible DNS runs for statistically planar turbulent premixed flames have been carried out using a DNS code called SENGA [28] in which the conservation equations of mass, momentum, energy, and species are solved in nondimensional form. The simulations have been carried out in a rectangular domain of size . The simulation domain has been discretised using a Cartesian mesh of size 345 × 230 × 230 with uniform grid spacing in each direction. The mean direction of flame propagation is assumed to align with the -direction and the transverse directions are assumed to be periodic. The boundaries in the direction of mean flame propagation are taken to be partially nonreflecting and the boundary conditions are specified using the Navier-Stokes characteristic boundary conditions (NSCBC) technique [29]. The spatial discretisation is carried out using a 10th order central difference scheme for the internal points and the order of discretisation drops gradually to a one-sided 2nd order scheme at nonperiodic boundaries. The time advancement is carried out using a low storage 3rd order Runge-Kutta scheme [30]. The turbulent velocity fluctuations are initialised using an initially homogeneous isotropic field using a pseudospectral method proposed by Rogallo [31]. The combustion is initailised by a steady unstrained planar laminar premixed flame solution. For the present study, the thermophysical properties such as dynamic viscosity , thermal conductivity , specific heats and and the density-weighted mass diffusivity are taken to be constant and independent of temperature following several previous studies [11–18, 32–34]. For the present study, standard values are taken for Prandtl number (Pr = 0.7) and the ratio of specific heats (). The flame Mach number is taken to be . The Lewis number Le is taken to be equal to unity for all the cases considered here. As and are often used to characterise the combustion situation on combustion regime diagram, the variation of in the present study is brought about by modifying and independently of each other. In order to achieve this, the values of both and are modified as both and are dependent on the values of and . The initial values of normalised rms turbulent velocity fluctuation , integral length scale to flame thickness ratio , heat release parameter , Damköhler number , Karlovitz number and turbulent Reynolds number are listed in Table 1. From Table 1, it is evident that remains greater than unity for all cases, which indicates that the combustion belongs to the thin reaction zones regime according to the regime diagram by Peters [1]. In cases B, C, and D the values of and were chosen in such a manner as to vary the turbulent Reynolds number by changing the Damköhler number Da while keeping the Karlovitz number constant (see (9)). On the other hand, in cases A, C, and E, and were chosen to vary by changing , while was kept constant. The range of turbulent Reynolds number values considered in this study remains modest although several previous studies [3–22] with comparable values of have made valuable contributions to the fundamental understanding of displacement speed in turbulent premixed flames. It is also worth noting that the range of considered here is comparable to that of previous laboratory-scale experiments (e.g., Kobayashi et al. [35]). Despite the limited range, a number of important effects on the displacement speed statistics have been identified in the present parametric study and are discussed in detail in Section 4 of this paper.
In all cases, the flame-turbulence interaction takes place under conditions of decaying turbulence, in which simulations should be carried out for a minimum time , where is the initial eddy turn over time and is the chemical time scale. All the cases considered here were run for one chemical time scale , which corresponds to a time equal to in case D, in cases A, C, and E, and for case B. The present simulation time is comparable to the simulation times used for previous DNS studies which focused on the analysis of displacement speed statistics [3–21, 23–27]. By the time statistics were extracted, the global turbulent kinetic energy and volume-averaged burning rate were no longer changing rapidly with time. The temporal variations of these quantities were presented in [36], and thus are not repeated here for the sake of conciseness. Analysis of the volume-averaged burning rate indicates that the effects of unsteadiness have become weak by the time when statistics were extracted. At this time, the global level of turbulent velocity fluctuation had decayed by 52.66%, 61.11%, 45%, 24%, and 34% in comparison to the initial values for cases A to E, respectively. By contrast, the integral length scale had increased by factors of between 1.5 to 2.25, ensuring that sufficient numbers of turbulent eddies were retained in each direction to obtain useful statistics. The values for , , and at the time when statistics were extracted have been presented in Table 2 of [36]. For all cases, the thermal flame thickness remains greater than the Kolmogorov length scale at that time, which indicates the combustion in all cases still belongs to the thin reaction zones regime. Moreover, the turbulent Reynolds number for all cases remains greater than unity (i.e., ) for the values of and at the time when statistics were extracted, which suggests that all the flames remain turbulent according to the regime diagram by Peters [1].
4. Results and Discussion
4.1. Flame-Turbulence Interaction
The contours of the reaction progress variable in the x1-x2 mid-plane of the domain at time are shown in Figures 1(a)–1(e) for cases A to E, respectively. It is evident that the extent of the wrinkling of the isosurfaces increases with increasing , that is, in going from case A through to case E (see Table 1). Moreover, the contours of representing the preheat zone (i.e., ) are much more distorted than those representing the reaction zone (i.e., ). This tendency is more prevalent for the case with the highest Karlovitz number (case E), since the scale separation between and increases with increasing , allowing more energetic eddies to enter into the preheat zone and causing greater distortion of the flame.
(a)
(b)
(c)
(d)
(e)
4.2. Probability Density Functions of Displacement Speed
The probability density functions (pdfs) of normalised displacement speed on the , and 0.9 isosurfaces are shown in Figures 2(a)–2(e), respectively. It is evident in all cases that there is a small but finite probability of finding negative displacement speed although the mean value of remains strongly positive throughout the flame brush. Negative values of displacement speed suggest that the progress variable isosurface in question is moving locally in a direction opposite to that of mean flame propagation. This situation arises when the negative contribution of the molecular diffusion rate dominates over the positive semi-definite chemical reaction rate, as discussed in detail by Gran et al. [4]. Moreover, Figure 2(e) indicates that the probability of finding negative values of is greatest for the highest Karlovitz number (case E). This behaviour can be explained in terms of the scaling analysis of Peters [1] for unity Lewis number flames, which suggested that the different components of displacement speed scale as where is the Kolmogorov velocity scale. Here, in accordance with Peters [1], () is taken to scale with and the curvature is taken to scale with the Kolmogorov length scale . These scalings of () and were confirmed by DNS data presented in [21]. Equation (10) clearly suggests that in the thin reaction zones regime (i.e., ) the effects of () are likely to weaken progressively in comparison to the contribution of with increasing Karlovitz number. This suggests that for large values of Karlovitz number the negative contribution of more readily overcomes the predominantly positive contribution of () to yield a negative value of . This eventually gives rise to an increased probability of finding negative values of for flames with increasing value of Karlovitz number, as suggested by Figure 2. In order to demonstrate that the observed negative values of originates principally due to , the pdfs of are shown in Figure 3, respectively, for all cases. Comparing Figures 2 and 3 reveals that the probability of finding negative values of is significantly smaller than that of finding negative values of .
(a)
(b)
(c)
(d)
(e)
(a)
(b)
(c)
(d)
(e)
Comparing Figures 2(a)–2(e) and Figures 3(a)–3(e), it can be seen that the most probable values of and remain of the order of and that the widths of the displacement speed pdf for a given isosurface increases with increasing (i.e., going from case A to case E). In order to explain this behavior, it is useful to examine the pdfs of the normalised tangential component of displacement speed which are depicted in Figure 4 and show that these pdfs remain almost symmetrical with a peak close to a value of zero. This behaviour is consistent with the statistically planar nature of these flames. However, the greater extent of flame wrinkling observed in Figure 1 for larger values of results in a greater spread of flame curvature, which in turn increases the width of the pdfs according to (5). The negligible mean value of essentially leads to the mean values of and becoming almost equal to each other, but the statistical variation of induced by leads to broader pdfs of than those of . These results clearly suggest that the distributions of and its components and generally broaden with , which essentially leads to broadening of the distributions with increasing when either or is held constant.
(a)
(b)
(c)
(d)
(e)
4.3. Strain Rate and Curvature Dependences of
For low Mach number unity Lewis number flames, both density and nondimensional temperature remain uniform on a given isosurface of and thus the statistical behaviours of and (see (5)) are mainly affected by the dependence of the SDF on strain rate and curvature. The values of conditional on normalised tangential strain rate are shown in Figure 5(a) for the isosurface at (close to the location of maximum reaction rate), which suggests clearly that and are positively correlated for all cases. The correlation coefficients between and are shown in Figure 5(b) for values of between 0.1 and 0.9, and these indicate that the positive correlation is maintained throughout the flame brush. This correlation can be explained by noting that the dilatation rate can be expressed as For unity Lewis number flames, the dilatation rate can be scaled as [21], whereas can be scaled as [37], which suggests that
(a)
(b)
Equation (11b) essentially shows that for low values of the Damköhler number, remains much greater than . This essentially suggests that an increase in tangential strain rate is associated with a decrease in the normal strain rate, since the normal strain rate scales as when dilatation rate effects are much weaker than those of turbulent straining. Under decreasing normal straining, the isoscalar lines are brought closer to each other and hence give rise to an increase in . The decreasing trend of normal strain rate with increasing tangential strain rate is reflected in the positive correlation between and throughout the flame brush. It can be seen from Figure 5 that the correlation between and does not show any monotonic trend with either Karlovitz number (Damköhler number) or turbulent Reynolds number variations when Damköhler number (Karlovitz number) is held constant.
The variation of conditional on normalised curvature is shown in Figure 6(a) for the isosurface at . Both positive and negative correlation trends can be observed, which is consistent with several previous studies [14–22, 26, 27]. As a result of these opposing trends, the net correlation between and turns out be weak throughout the flame brush, as evident from the correlation coefficients between and presented for different values of across the flame brush in Figure 6(b). Both the positively and negatively correlating branches can be explained using (11a). In turbulent premixed flames the dilatation rate is negatively correlated with curvature , because of focusing (defocusing) of heat in the negatively (positively) curved regions of the flame surface [18, 21]. The tangential strain rate is also found to be negatively correlated with curvature [18, 21]. The correlation coefficients for the and correlations for different values across the flame brush are shown in Figures 7(a) and 7(b), respectively. Comparing the correlation coefficients of the and correlations for different cases reveals that the strength of both of these negative correlations decreases with increasing and when either or is held constant.
(a)
(b)
(a)
(b)
The combination of the positive correlation between and and the negative correlation between and leads to small values of in regions of high positive curvature. This gives rise to the negative correlation trend between and , as observed in Figure 6. However, the dilatation rate may locally attain large values in the highly negatively curved regions of the flame due to strong focussing of heat. Hence, may locally exceed and induce an increase in according to (11a). As isoscalar lines move apart from each other under increasing normal strain rate, the increase in with increasing negative curvature leads to a positive correlation between and in negatively curved regions, as also observed in Figure 6.
4.4. Curvature Dependence of
The mean normalised displacement speed conditional on normalised curvature is shown in Figure 8(a) for , which indicates that is negatively correlated with curvature and that the correlation is nonlinear in nature. These observations are consistent with several previous DNS studies with both detailed and simplified chemistry [3–11, 13–18, 20–22]. The variation of the correlation coefficient between and for different isosurfaces across the flame brush is shown in Figure 8(b), which reveals that strength of the negative correlation decreases with increasing when either (see cases A, C, and E) or (see cases B, C, and D) is held constant (see Table 1). The tangential diffusion component of displacement speed is deterministically negatively correlated with curvature with correlation coefficient equal to unity for unity Lewis number flames (see (5)). Thus, the nonlinearity of the correlation between and originates due to the curvature dependence of and .
(a)
(b)
Mean values of the reaction component conditional on normalised curvature are shown in Figure 9(a) for the isosurface at . The plot shows both positive and negative correlation trends. For low Mach number unity Lewis number flames, the curvature dependence of is principally determined by the curvature dependence of . The negative (positive) correlation between and ultimately gives rise to the positive (negative) correlation between and as observed in Figure 9. As a result of the weak net correlation between and , the net correlation between and also turns out to be weak, as indicated by the correlation coefficients plotted in Figure 9(b) Note that Figure 9(b) includes isosurfaces of from 0.5 to 0.9 only, since the reaction rate is negligible in the preheat zone (see also (5)).
(a)
(b)
Corresponding mean values of the normal diffusion component conditional on normalised curvature are shown in Figure 10(a) for , while the correlation coefficients between and are shown in Figure 10(b) for values of c from 0.1 to 0.9. The magnitude of the normal molecular diffusion rate in the thin reaction zones regime can be scaled as where is the characteristic flame thickness which in turn can be scaled as . This suggests that scales as towards the reactant side of the flame brush and as towards the product side. As a result of this, the positively (negatively) correlating branch of the − correlation gives rise to a positive (negative) correlation between and towards the reactant side, where assumes predominantly positive values. By contrast, the positive (negative) correlation between and results in a negative (positive) correlation between and towards the product side where assumes predominantly negative values. Both positive and negative correlation trends can be discerned for all the isosurfaces shown in Figure 10(b). As a result of this, the net correlation between and turns out to be weak throughout the flame brush. The positive and negative correlating trends with observed for both and with give rise to positive and negative correlating branches in the variation of conditional on as shown in Figure 10(c) for . This nonlinear correlation between and ultimately induces a nonlinear curvature dependence of displacement speed (see Figure 8(a)), noting that the correlation coefficient arising from the deterministic relation between and remains close to −1.0 throughout the flame brush for all cases. The correlation coefficients between and are shown in Figure 10(d) for values of from 0.1 to 0.9, and demonstrate that the correlation is much weaker than the correlation in all cases. Hence, the net correlation between and turns out to be negative throughout the flame brush for all the cases considered here (see Figure 8(a)).
(a)
(b)
(c)
(d)
4.5. Tangential Strain Rate Dependence of
The mean values of normalised displacement speed conditional on normalised tangential strain rate are shown in Figure 11(a) for which suggests that and are positively correlated for all cases. The observed strain rate dependence of is found to be qualitatively consistent with previous DNS studies for flames with global Lewis number close to unity [3, 5, 12, 14, 15, 17–19, 21, 22]. The variations of the correlation coefficients between and for different isosurfaces are shown in Figure 11(b), which reveals that the strength of the correlation between and decreases with increasing when either or is held constant. This is consistent with the analytical treatment by Joulin [38] who indicated that the strain rate dependence of weakens with increasing frequency of straining. In order to explain the observed behaviour it is useful to look into the dependence of , , and on . For low Mach number unity, Lewis number flames the strain rate dependence of is essentially determined by the strain rate dependence of . It has been shown already that and are positively correlated (see Figure 5), which gives rise to a negative correlation between and according to (5). This is substantiated by the correlation coefficients between and as shown in Figure 11(c) for the isosurfaces from c = 0.5 to 0.9. The correlation between and does not show any monotonic trend with either Karlovitz number (Damköhler number) or turbulent Reynolds number variations when the Damköhler number (Karlovitz number) is held constant, due to the lack of any such trend in the correlation.
(a)
(b)
(c)
(d)
(e)
(f)
It has been discussed already that scales as on the reactant side of the flame brush, and as on the product side. Since there is a positive correlation between and , this gives rise to a positive (negative) correlation between and on the reactant (product) side. This behaviour can be seen in the correlation coefficients plotted in Figure 11(d), which also indicates that the correlation coefficient between and does not show any monotonic trend with changes in either or . However, it has been found that the correlation between and in general remains stronger for smaller values of when either or is held constant.
Since remains negligible towards the reactant side of the flame brush and is positively correlated with , the correlation coefficient between and remains positive towards the reactant side (see Figure 11(e)). By contrast, both and are negatively correlated with towards the product side, which leads to a locally negative correlation between and .
Since and are negatively correlated, the tangential diffusion component is positively correlated with because is proportional to the negative of curvature (i.e., ). This positive correlation is evident from Figure 11(f), in which the correlation coefficient between and is shown for isosurfaces of across the flame brush. The negative correlation between and (see Figure 7(a)) arises principally due to heat release, and thus, the strength of this correlation weakens with increasing Karlovitz number, since the effects of turbulent velocity fluctuations are likely to mask the heat release effects as the combustion tends more towards the broken reaction zones regime [1]. Moreover, with decreasing Damköhler number the effects of heat release weaken in comparison to those of turbulent straining, as suggested by (11b). As a result of this, the strength of the negative correlation between and weakens with decreasing (increasing) Damköhler number (Karlovitz number) when Ka (Da) is held constant (see Figure 7(a)), and this gives rise to a corresponding change in the correlation between and (see Figure 11(f)). This essentially suggests that the strength of the positive correlation between and decreases with increasing turbulent Reynolds number when either or is held constant. The effects of the positive and negative correlations remain strongest in the reaction zone and the strengths of these correlations decrease towards both the reactant and product sides of the flame brush (see Figures 11(f) and 7(b), resp.).
It is evident from Figures 11(e) and 11(f) that the positive correlations between and , and between and , aid each other towards the reactant side of the flame brush. By contrast, towards the product side the positive correlation between and overcomes the negative correlation between and to result in a net positive correlation between and in all the cases considered here (see Figure 11(a)). The weakening of the correlation between and with increasing ultimately leads to decrease in the strength of the correlation between and when either or is held constant (see Figures 11(a) and 11(f)).
4.6. Implications for FSD Modelling
In the context of Reynolds Averaged Navier Stokes (RANS) modelling, the generalised FSD is defined as [39] where indicates the Reynolds averaged value of a general quantity . It has been demonstrated in [15, 23, 26, 27] that the statistical behaviour of and the terms of its transport (7) in the thin reaction zones regime depend significantly on the choice of progress variable isosurface . Hence, it is more appropriate to use the generalised FSD (see(12a)) which does not depend on the choice of isosurface. Reynolds averaging of (7) yields the transport equation for the generalised FSD [2, 23, 24] where indicates a surface averaging operation for a general quantity [39]. The first term on the left hand side is the transient term while the second term is the advection term. The three terms on the right hand side represent the effects of strain rate, curvature, and propagation, respectively. It is evident from (12b) that the statistical behaviour of the curvature and propagation terms depends on the statistics of displacement speed. These terms can be rewritten in terms of the displacement speed components as
The variations of all of the terms appearing in these equations with Favre averaged reaction progress variable are shown in Figure 12. The Reynolds averages are evaluated by ensemble averaging the relevant quantities in the directions normal to the direction of mean flame propagation.
It can be seen from Figure 12 that the FSD curvature term (line with tick symbols) assumes mostly negative value throughout the flame brush for all cases although small positive contributions are found towards the reactant side of the flame brush for small Reynolds numbers (i.e., cases A and B). The magnitude of the term (line with filled triangles) remains comparable to the contribution of the term (line with open squares) for small values of (i.e., cases A and B). However, the contribution of remains smaller than that of for larger values of and the statistical behaviour of the FSD curvature term is then principally governed by the latter contribution. It is also evident from Figure 12 that the magnitude of () becomes increasingly more important than the magnitude of with increasing Karlovitz number (Damköhler number) when () is held constant. The statistical behaviour of is determined by the correlation between and and between and and the strength of both of these correlations tends to weaken with increasing (increasing) Karlovitz number (Damköhler number) when () is held constant (see Figures 6(b) and 10(d)), causing a corresponding reduction in the relative magnitude of the term .
The FSD propagation term (line with crosses) remains positive towards the reactant side, but it eventually becomes negative towards the product side. It can be seen that the term (line with filled squares) remains almost equal to the FSD propagation term itself, because the contribution of (line with open circles) remains negligible throughout the flame brush due to the statistically planar nature of the flames considered here. The contribution of remains positive towards the reactant side and assumes negative values towards the product side.
In some previous studies, the FSD curvature and propagation terms have been modelled by replacing with to yield the modelled expressions and , respectively; [40, 41]. Figure 12 further indicates that there are significant differences between the original terms and these modelled expressions, which are shown, respectively, by a line with filled circles and a line with open triangles. This essentially indicates that cannot be approximated simply by in the modelling of the FSD curvature and propagation terms. Moreover, Figure 12 suggests that modelling of is necessary for the closure of the FSD curvature term especially in the thin reaction zones regime since the term remains close to throughout the flame brush for high values of the turbulent Reynolds number. The effects of turbulent Reynolds number on the modelling of the FSD curvature and propagation terms are beyond of the scope of this study and will be addressed in future work.
5. Conclusions
The turbulent Reynolds number dependence of the statistics of displacement speed in turbulent premixed flames has been studied using a DNS database of statistically planar flames, in which the variation of turbulent Reynolds number from 20 to 100 has been brought about by modifying either Damköhler number or Karlovitz number independently of each other. As both and are functions of and , the variation of for the present database has been achieved by modifying the initial values of and simultaneously.
The mean value of displacement speed remains positive throughout the flame brush, but there is nonzero probability of finding negative values of displacement speed, in accordance with the findings of previous work. It has been shown here that the probability of finding negative displacement speed increases with increasing turbulent Reynolds number when the Damköhler number is held constant.
The dependences of tangential strain rate and dilatation rate on flame curvature are shown to have a significant influence on the strain rate and curvature dependences of the SDF , which in turn affects the statistical behaviour of displacement speed in response to strain rate and curvature. It has been found that the variation of turbulent Reynolds number does not alter the qualitative nature of the correlations between tangential strain rate and dilatation rate with curvature, but the strength of these correlations is found to weaken with increasing turbulent Reynolds number when either Damköhler or Karlovitz number is held constant. Similarly, the strain rate and curvature dependences of displacement speed are shown to weaken with increasing turbulent Reynolds number when either Damköhler or Karlovitz number is kept unaltered.
Detailed physical explanations are provided to explain the influences of turbulent Reynolds number on the strain rate and curvature dependences of displacement speed in terms of the individual response of the reaction, normal diffusion and tangential diffusion components of displacement speed.
The implications of the turbulent Reynolds number dependence of displacement speed for the modelling of the FSD transport equation has been explored. The statistical behaviour of the FSD curvature term contributions arising due to the combined reaction and normal diffusion components of displacement speed and to the tangential diffusion component of displacement speed have been analysed in detail. At low turbulence Reynolds numbers, the magnitudes of these contributions remain comparable. As the turbulence Reynolds number increases, the curvature dependence of the tangential component of displacement speed ensures that this contribution becomes increasingly important and is the predominant influence on the FSD curvature term for larger values of turbulent Reynolds number. Hence, the curvature dependence of the tangential diffusion component of displacement speed cannot be ignored for accurate modelling of the FSD curvature term for flames within the thin reaction zones regime.
It is worth noting that the effects of differential diffusion of heat and mass are not taken into account in the present study. However, some previous DNS studies with detailed chemistry [3–7, 9] suggested that intermediate species may play an important role even for the flames with global Lewis number close to unity. Moreover, only a moderate range of turbulent Reynolds number has been considered in this study, so three-dimensional DNS studies with detailed chemistry and at higher values of turbulent Reynolds number will be needed for deeper understanding and for the purpose of quantitative predictions.
Acknowledgment
N. Chakraborty gratefully acknowledges the financial assistance of the EPSRC, UK.