Vol.:(0123456789) Flow, Turbulence and Combustion (2022) 109:101–123 https://doi.org/10.1007/s10494-021-00315-5 1 3 Turbulence Intensity and Length Scale Effects on Premixed Turbulent Flame Propagation Shrey Trivedi1 · R. S. Cant1 Received: 2 August 2021 / Accepted: 15 December 2021 / Published online: 29 December 2021 © The Author(s) 2021 Abstract The effects of varying turbulence intensity and turbulence length scale on premixed tur- bulent flame propagation are investigated using Direct Numerical Simulation (DNS). The DNS dataset contains the results of a set of turbulent flame simulations based on separate and systematic changes in either turbulence intensity or turbulence integral length scale while keeping all other parameters constant. All flames considered are in the thin reaction zones regime. Several aspects of flame behaviour are analysed and compared, either by varying the turbulence intensity at constant integral length scale, or by varying the integral length scale at constant turbulence intensity. The turbulent flame speed is found to increase with increasing turbulence intensity and also with increasing integral length scale. Changes in the turbulent flame speed are generally accounted for by changes in the flame surface area, but some deviation is observed at high values of turbulence intensity. The probability density functions (pdfs) of tangential strain rate and mean flame curvature are found to broaden with increasing turbulence intensity and also with decreasing integral length scale. The response of the correlation between tangential strain rate and mean flame curvature is also investigated. The statistics of displacement speed and its components are analysed, and the findings indicate that changes in response to decreasing integral length scale are broadly similar to those observed for increasing turbulence intensity, although there are some interesting differences. These findings serve to improve current understanding of the role of turbulence length scales in flame propagation. Keywords  Turbulent premixed flames · Turbulence intensity · Turbulence integral length scale 1  Introduction The response of a premixed flame when subjected to increasingly intense turbulence has been the subject of many different experimental and theoretical studies (Bradley 1992; Driscoll 2008; Peters 1999; Poludnenko and Oran 2011, 2010). Turbulence causes flames to stretch, resulting in an increase of the flame surface area and hence * Shrey Trivedi st634@cam.ac.uk 1 Department of Engineering, University of Cambridge, Cambridge, UK http://crossmark.crossref.org/dialog/?doi=10.1007/s10494-021-00315-5&domain=pdf 102 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 the burning rate (Damköhler 1940; Pope 1988; Candel and Poinsot 1990; Law 1988). In addition, it has been theorised that small scales of turbulence may enhance the trans- port due to turbulence and cause further enhancement in the flame speed (Damköhler 1940; Nivarti et al. 2019; Driscoll 2008). A majority of related studies has focused on the variations in turbulence intensity (Driscoll 2008; Lipatnikov and Chomiak 2002), including Direct Numerical Simulation (DNS) studies which have been applied to investigate the effect of increasing the turbulence intensity u′ while holding all other parameters constant as far as possible (Nivarti and Cant 2017). There have been fewer studies which have produced information about the effect on the flame of varying the turbulence integral length scale l0 (Peters 1999; Lipatnikov and Chomiak 1999; Song et al. 2021). Several experimental investigations dealing with variations in �0 have been discussed in Sect.  3.3.2 in Ref. (Lipatnikov and Chomiak 2002). This is harder to achieve in DNS, due to the need to resolve the full range of tur- bulence length scales from the integral scale down to the Kolmogorov scale � . There- fore, given the high computational cost, the scope for variation of the integral length scale l0 is necessarily limited. Nevertheless, fundamental study of the effect of length scales is important in order to explore the combustion physics, to define the constraints on flamelet-based modelling (Driscoll et al. 2020; Skiba et al. 2017, 2021; Hawkes and Chen 2006) and to suggest model improvements. Recent papers have analysed the vari- ations in flame speed with �0 both experimentally (Kim et al. 2020; Mohammadnejad et al. 2021) and using DNS (Yu and Lipatnikov 2017a, b; Kulkarni and Bisetti 2021a; Kulkarni et al. 2021; Kulkarni and Bisetti 2021b; Varma et al. 2021). The dependence of flame speed on the turbulence Reynolds number was investi- gated experimentally (Smith and Gouldin 1979; Shy et  al. 2015) and showed any increase in flame speed with an increasing turbulent Reynolds number, consistent with the previous predictions (Damköhler 1940; Peters 1999). Recent analysis based on DNS of turbulent slot jet flames (Attili et al. 2021; Luca et al. 2019) investigated flames at nearly constant Karlovitz number (Ka) but with variations in Reynolds num- ber (Re). The key finding of this work was that thickening of the inner reaction layer occurs as Re increases, causing the flame speed to increase faster than the flame area as Re increases (Attili et al. 2021). Furthermore, the statistics of flame stretch for this dataset (Luca et al. 2019) indicated that tangential strain rate and mean flame curvature normalised by the Kolmogorov length � were found to have a weak dependence on Re. The present study aims to investigate the response of a turbulent premixed flame to variations in u′ and l0 , treating each of these quantities separately as an independent variable while holding all the other variables constant. Unlike the previous analyses (Attili et  al. 2021; Luca et  al. 2019), the Komogorov scale and hence the Karlovitz number changes across the dataset. The rest of the paper is organised as follows. Section 2 presents the relevant theory to outline and define the flame properties of interest. A matrix of test cases is described in Sect. 3 covering a range of both of the independent variables of interest. Section 4 presents the results for the flame properties including flame speed and flame area, sta- tistics of tangential strain rate and mean flame curvature, and statistics of displacement speed and its components. The effects of varying u′ at constant l0 are compared with previous findings (Nivarti and Cant 2017), while results illustrating the effects of vary- ing l0 at constant u′ are also presented and discussed. 103Flow, Turbulence and Combustion (2022) 109:101–123 1 3 2 � Theory The increase in flame surface area AT due to the action of turbulence (Law 1988; Pope 1988) and the subsequent increase in the turbulent flame speed sT has been studied extensively both experimentally (Driscoll 2008) and computationally (Chen et al. 1999; Poludnenko and Oran 2011; Nivarti and Cant 2017), and is governed by the expression where I0 is the stretch factor (Driscoll 2008; Bray and Cant 1991). Note that the above rela- tionship is valid under the Klimov-Williams limit (Klimov 1963; Williams 1985), i.e. when all turbulence length scales are larger than �L . In the above expression, the turbulent flame speed sT can be evaluated using the integral where �R and YFR are the unburnt mixture density and fuel mass fraction respectively and A0 is the cross sectional area of the domain. It is convenient to introduce a reaction pro- gress variable c (based for example on normalised fuel mass fraction) which rises montoni- cally from zero in fresh reactants to unity in fully-burned products. The relation between flame area and flame surface density is stated as The flame surface density Σ is important in premixed combustion modelling (Pope 1988; Peters 1999; Candel and Poinsot 1990; Cant et al. 1990), and can be defined as the expected value of the flame surface area per unit volume for an isosurface at c = c∗ , i.e. where the “fine–grained” flame surface density Σ� can be defined by using the gradient of progress variable (Boger et al. 1998) according to The balance equation for Σ takes the form Pope (1988), Candel and Poinsot (1990) where the flame normal vector ni is given by and the displacement speed Sd of the flame is obtained using the expression (1) sT sL ∼ AT AL => sT sL = AT AL I0 (2)sT = − 1 𝜌RYFRA0 ∫V 𝜔̇FdV (3)AT = ∫ ∞ −∞ �dV (4)Σ = Σ� (5)Σ� = �A �V = |∇c|�(c − c ∗) (6) 𝜕Σ 𝜕t + 𝜕 𝜕x i (u i + S d n i )Σ = ṠΣ (7)ni = − ∇c |∇c| , (8)𝜌Sd𝜎 = 𝜔̇ + 𝜕 𝜕xk ( 𝜌D 𝜕c 𝜕xk ) 104 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 where � is the generalised FSD (Boger et al. 1998) defined as � = |∇c| . The flame stretch rate Ṡ is defined as the fractional rate of increase of flame area where at is the rate of hydrodynamic strain tangential to the flame and � is the mean curva- ture of the flame. An expression for at can be stated as where �ij is the Kronecker delta, and the mean curvature � is defined as Flame stretch can have significant effects on the internal structure of the flame (Law 1988). These effects can be accounted for using the stretch factor I0 (Bray and Cant 1991; Driscoll 2008), defined as the ratio of local instantaneous propagation speed s∗ L to the unstrained laminar flame speed sL: The quantity I0 has been found to be ∼ 1 for thermo-diffusively neutral mixtures, i.e. Le = 1 , but the value of I0 has been found to deviate from unity for non-unity Lewis num- bers (Hawkes and Chen 2006; Bradley 1992; Lipatnikov and Chomiak 2002) or for very high values of Karlovitz number (Driscoll 2008). Moreover, thickening of the flame can enhance the flame speed over and above that estimated by change in flame area alone (Attili et al. 2021; Nivarti et al. 2019) leading to values of I0 greater than unity. 3 � DNS Dataset for the Current Study The simulations were carried out using the well-established in-house DNS code SENGA2 (Jenkins and Cant 1999; Cant 2012). This code solves the Navier–Stokes momentum equa- tions in their compressible form, augmented by conservation equations for mass and energy and by balance equations for the relevant chemical species mass fractions. The spatial dis- cretisation scheme uses 10th order explicit centred finite differences, reducing smoothly to 4th order one-sided differences at non-periodic boundaries. Time-stepping is effected using a 4th-order low–storage explicit Runge–Kutta algorithm with adaptive step-size control (Kennedy and Carpenter 2003). The code is fully parallelised using a domain decomposi- tion technique. Boundary conditions are imposed using the Navier–Stokes Characteristic Boundary Condition (NSCBC) formalism (Poinsot and Lele 1992). For the purposes of the present study, the chemistry is treated using a single-step Arrhe- nius formulation tuned to replicate the flame propagation speed sL of a stoichiometric meth- ane—air flame, i.e. 39 cm/s, and also to capture the corresponding thermal flame thickness �L given as �L = (Tb − Tu)∕(dT∕dx)max , where Tb and Tu are the temperature of the burnt and unburnt gas respectively. The density ratio of the burned to unburned gases is ∼ 0.13. (9)Ṡ = 1 A dA dt = at + 2Sd𝜅 (10)at = (�ij − ninj) �ui �xj (11)� = 1 2 �ni �xi (12)I0 = s∗ L sL 105Flow, Turbulence and Combustion (2022) 109:101–123 1 3 The Lewis number of the reacting species is set to unity. It is acknowledged that single— step chemistry may not capture all the details of the flame—turbulence interaction, and this will be investigated in future work using a more detailed chemical reaction mechanism. The computational domain consists of a cuboid with periodic boundary conditions in the spanwise y and z directions. The streamwise x direction has NSCBC outflow boundary conditions on opposite faces. The thermochemical pressure is set to 1 bar and the reactant temperature is set to 300 K. The dataset created for analysis is based on a pair of statisti- cally planar flames facing each other with fresh reactants between them. The flames are simulated in a cuboidal domain of cross section 5 mm × 5 mm. Turbulence of the required intensity and integral length scale is generated using a spectral algorithm (Orszag 1972) initialised with the Bachelor–Townsend spectrum (Batchelor and Townsend 1948) and is imposed on the domain. The turbulence in the unburned reactants is forced linearly (Lun- dgren 2003) in order to counteract the process of decay and hence to maintain the required properties. A known shortcoming of Lundgren’s method (Lundgren 2003) is that the tur- bulent length scale changes over time until it reaches a steady value (Klein et  al. 2017; Rosales and Meneveau 2005). Klein et  al. (2017) attempted to fix this issue by making the forcing term proportional to a high pass filtered value of velocity. However, no such correction has been applied in this work and it is likely that the associated length scales of turbulence change in time. In future, this issue will be addressed.. Within the dataset, both turbulence intensity u′ and integral length scale �0 are var- ied systematically. For values of u�∕sL = 10, 20 and 30, simulations are carried out with integral length scales corresponding to �0∕�L = 5.0 , �0∕�L = 2.0 and �0∕�L = 1.25 . For u�∕sL = 40 only the case of �0∕�L = 2.0 is simulated. The mesh spacing �x is varied in each case to match the Kolmogorov scale � and the simulation time was kept at 6 eddy turnover times, i.e. t∕� = 6 for all cases. A summary of the parameters for each simulation is provided in Table 1. The separation between the flames is varied for each individual case such that the flames do not collide even at the end of the run. The forcing region was kept between the two flames and of the forced region varied between cases since each case had different separa- tion between flames. The width of the forced region was chosen such that it does not inter- act with the flames during the run time. Figure 1 shows results from the cases with three different �0 values for u�∕sL = 10 . For the highest �0 value (Fig. 1a), the large scale wrin- kling requires larger separation between the two flames whereas for the smaller �0 values (Fig. 1b and c) a smaller separation is sufficient to keep the two flames from interacting. Table 1   Main parameters of the different cases in the dataset �0∕�L Ka Da Re � ( �m) �� ( �s) �x ( �m) u �∕s L = 10 5.0 14.14 0.5 50 39.27 11.03 41.66 2.0 22.36 0.2 20 24.85 6.97 26.04 1.25 28.28 0.125 12.5 20.33 5.51 20.833 u �∕s L = 20 5.0 40.0 0.25 100 28.76 3.9 26.04 2.0 63.25 0.1 40 18.19 2.46 14.38 1.25 80.0 0.0625 25 14.38 1.95 13.02 u �∕s L = 30 5.0 73.5 0.166 150 23.48 2.12 22.72 2.0 116.2 0.0667 60 14.85 1.34 13.02 1.25 147.0 0.0417 37.5 11.90 1.06 11.90 u �∕s L = 40 2.0 179 0.05 80 12.86 0.87 11.90 106 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 The dataset can be used to analyse the effects of integral length scale �0 while keep- ing the turbulence intensity u′ fixed. Alternatively, a fixed value of �0 can be chosen and the variation in u′ can be studied. Note that the Kolmogorov scale � also decreases as �0 decreases and/or u′ increases. The data points corresponding to the dataset are presented on a Borghi diagram as shown in Fig. 2. All the cases are either in, or very close to, the thin reaction zones regime. It has been noted that broken reaction zones are not observed even well beyond Ka=100 (Aspden et al. 2019; Song et al. 2021), which was also found to be true for all of the present cases. In all cases, the turbulent flame approached a statisti- cally stationary state, albeit with considerable fluctuations of sT about its mean value, as observed in other work (Poludnenko and Oran 2011; Lipatnikov et al. 2015). 4 � Results In this section, the flame properties are presented as functions of the systematic variations in both u′ and �0 . The results for increasing u′ at a fixed value of �0 are found to be consistent with those from the literature (Nivarti and Cant 2017; Lipatnikov and Chomiak 2002). The (a) �0/δL = 5.0 (b) �0/δL = 2.0 (c) �0/δL = 1.25 Fig. 1   Two dimensional slices of the 3-D snapshots of product mass fraction for cases �0∕�L = 5.0 , �0∕�L = 2.0 and �0∕�L = 1.25 (top to bottom) for fixed u�∕s L = 10 . The separation between the flames is adjusted to keep them from interacting 107Flow, Turbulence and Combustion (2022) 109:101–123 1 3 results for varying �0 at constant u′ are presented for analysis, and comparisons are drawn between the known trends for increasing u′ versus those observed for varying �0. Turbulent flame speed and flame area are analysed first, followed by the probability den- sity functions (pdfs) of mean curvature, tangential strain rate and displacement speed. Subse- quently, the correlations between strain, curvature, displacement speed and its components are evaluated and presented. 4.1 � Turbulent Flame Speed The turbulent flame speed sT and flame area AT (for the isosurface at c = 0.8) are evaluated using Eqs. 2 and 3 respectively and are normalised with the unstrained laminar flame speed sL and the domain cross—sectional area AL for the full dataset. The turbulent flame speed sT is plotted against u′ (also normalised with sL ) in Fig. 3 with different lines representing different values of �0 . It can be seen that sT increases with u′ for all values of �0 . Furthermore, the val- ues of sT are higher for larger values of �0 , i.e. This behaviour is consistent with previous simulation results for increasing u′ (Nivarti and Cant 2017; Lipatnikov and Chomiak 2002) as well as with previous analyses (Peters 1999; Lipatnikov and Chomiak 1999). Note that the values of sT∕sL obtained here are similar to sT ,�0∕𝛿L=5.0 > sT ,�0∕𝛿L=2.0 > sT ,�0∕𝛿L=1.25 Fig. 2   Representation of the current DNS dataset on a Borghi diagram Fig. 3   Turbulent flame speed s T ∕s L plotted against turbulence intensity u�∕s L for three different values of integral length scale �0∕�L equal to 5.0, 2.0 and 1.25 108 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 those obtained in the constant density flows (for example in Sabelnikov et al. (2019)) as well as in variable density flow (Nivarti and Cant 2017), consistent with the discussion in Lipatnikov et al. (2018). Values of sT∕sL are also presented in Table 2 along with the corresponding values of AT∕AL and the ratio I0 = (sT∕sL)∕(AT∕AL) (see Eq. 1). There are several important obser- vations to be made from this table. First, the area ratio AT∕AL follows the same trend as sT∕sL in that values corresponding to �0∕�L = 5.0 are higher than those corresponding to �0∕�L = 2.0 and �0∕�L = 1.25 for each u′ value. It can be seen that sT∕sL (and also AT∕AL ) differs by around a factor of two between �0∕�L = 5.0 and �0∕�L = 1.25 . This indicates that the principal factor affecting the turbulent flame speed is the increase of flame sur- face area, due to either higher turbulence intensity or larger turbulence integral length scale causing increased large-scale flame wrinkling. Following Refs. (Kulkarni and Bisetti 2021a; Chaudhuri et al. 2012), a power—law fit for the dependence of ST on Reynolds number in the form ST ∼ Ren was found to yield an exponent n ∼ 0.4 which is similar but lower than the findings of Chaudhuri et al. (2012) ( n ∼ 0.5 ) and Kulkarni et al. (2021) ( n ∼ 0.56 ). This could be a result of the differences in flow configuration, fuel type or the fact that the Reynolds numbers used here are lower than in these studies. This merits further investigation using a DNS dataset covering a greater range of Reynolds number than is accessible in the present work. At the same time, the value of I0 becomes greater than unity as u′ increases for the two cases with larger �0 , and especially for the case with the largest value of �0∕�L = 5.0 . Like- wise, the values of I0 also increase slightly as u′ increases for any fixed value of �0 . This indicates that factors other than the flame area, such as changes to the internal structure of the local flame sheet due to strain and curvature may affect the observed increase in the turbulent flame speed. However, the values of I0 remain very close to unity for all the cases and hence require further investigation. An increase in I0 with Reynolds number was reported by Attili et al. (2021) caused by the thickening of the inner flame layer. In the pre- sent study, the values of I0 for any given u′ do not change significantly as the values of �0 change. For each value of u′ , it can be seen from Table 2 (4th column) that the decrease in AT∕AL as �0 decreases is monotonic but is not linear. The ratio (AT∕AL)∕(�0∕�L) is presented in the rightmost column of the table. This ratio is lowest for the cases with �0∕�L = 5.0 and increases for the cases with smaller length scales. This provides an indication that an Table 2   Turbulent flame speed s T and flame surface area A T (evaluated at c = 0.8 ) normalised by unstrained laminar flame speed s L and domain cross section area A L respectively for all cases Values of I0 and the ratio (A T ∕A L )∕(�0∕�L) are also presented u �∕u L �0∕�L s T ∕s L A T ∕A L I0 (A T ∕A L )∕(�0∕�L) 10 5.0 3.237 3.296 0.982 0.659 2.0 2.287 2.299 0.995 1.149 1.25 1.758 1.776 0.989 1.421 20 5.0 3.805 3.830 0.993 0.766 2.0 2.753 2.704 1.018 1.352 1.25 2.256 2.194 1.028 1.755 30 5.0 5.023 4.907 1.024 0.981 2.0 3.571 3.398 1.051 1.699 1.25 2.685 2.619 1.024 2.095 40 2.0 4.628 4.160 1.112 2.080 109Flow, Turbulence and Combustion (2022) 109:101–123 1 3 increasing amount of small-scale turbulence plays a role in increasing the flame surface area, and is more efficient in doing so at lower values of �0. The flame surface density Σ conditioned on progress variable c is presented in Fig. 4. Figure 4a shows the plots for fixed �0∕�L = 2.0 and increasing u′ while Fig. 4b represents those for a fixed u�∕SL ( = 20) but decreasing �0 . The conditional value of Σ decreases as u′ increases. This is because even though the flame area increases, the volume correspond- ing to the wrinkled flame surface also increases. It is interesting to note that the same trend of decreasing Σ is observed also for decreasing �0 at fixed u′ , and indeed the similarity between Figs. 4a and b is striking. Note that there are no significant changes in the profiles of Σ for c > 0.6 and hence the changes in Σ are mostly restricted to the preheat zone, con- sistent with the previous findings (Nivarti and Cant 2017). Fig. 4   Flame surface density Σ conditioned on progress variable c for a varying u′ at fixed �0 = 2.0� L and b varying �0 at fixed u � = 20 s L (a) Flame surface density for different u′ at a fixed �0 = 2.0δL. (b) Flame surface density for different �0 at a fixed u′ = 20 sL. 110 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 4.2 � Curvature and Strain Rate In this subsection, the pdfs of mean curvature � and tangential strain rate at are plotted for different values of u′ and �0 . These pdfs are analysed for the isosurface correspond- ing to c = 0.8 . First, the pdfs of curvature and tangential strain rate are presented in Fig.  5 for varying u′ at fixed values of �0 . The curvature pdfs are shown in the left column while the tangential strain rate pdfs are shown in the right column. The pdfs (a) Pdfs of mean curvature (left) and tangential strain rate (right) at �0/δL = 5.0. (b) Pdfs of mean curvature (left) and tangential strain rate (right) at �0/δL = 2.0. (c) Pdfs of mean curvature (left) and tangential strain rate (right) at �0/δL = 1.25. Fig. 5   Pdfs of mean curvature (left column) and strain rate (right column) for different values of u�∕s L and fixed �0∕�L = 5.0 (top row), 2.0 (middle row) and 1.25 (bottom row) 111Flow, Turbulence and Combustion (2022) 109:101–123 1 3 corresponding to �0∕�L = 5.0, 2.0 and 1.25 are shown in the top, middle and bottom row respectively. Looking at the mean curvature pdfs (Fig.  5, left column), as u′ is increased from 10 (red line) to 20 (blue line) to 30 (yellow line), the pdfs become wider and their peak value decreases, and this is true for all values of �0∕�L . On physical grounds, the increase in the width of the mean curvature pdf is expected since at higher u′ the eddies have more energy available to wrinkle the flame surface. Furthermore, at higher u′ there is an addition of Kolmogorov—scale eddies since 𝜂u�∕sl=10 > ... > 𝜂u�∕sL=30  . These additional scales lie in the range 𝜂u�∕sL=10 > l > 𝜂u�∕sL=20 for the u�∕sL = 20 case and likewise for u�∕sL = 30 case. These additional smaller scales result in higher values of � ( ∼ 1∕R ) and hence broader pdfs. The behaviour of these pdfs is consistent with the results in a previous analysis by Nivarti and Cant (2017) and for the experimental results for V-flames (Kheirkhah and Gül- der 2013). It is also interesting to note that the height and width of the mean curvature pdfs for u�∕sL = 20 and 30 are similar to each other, suggesting that there is some tendency towards saturation of flame curvature at higher turbulence intensity. For the tangential strain rate pdfs (Fig.  5, right column) the height of the peaks decreases and their width increases as u�∕sL increases, for all values of �0∕�L . The effect is most pronounced for the highest u′ cases for which the pdf becomes significantly broader which was also observed from the experimental results of spherically expanding flames (Liu et al. 2021). The peak of the pdf for all cases is observed at a positive value of at , and this value becomes more positive as u�∕sL increases which is consistent with previous DNS findings (Nivarti and Cant 2017). Again, these effects result from greater energy within the eddies at higher u�∕sL . The additional smaller scales strain the flame surfaces more effectively due to their higher velocity gradients. Note that the flame thickness �L and flame time �c(= �L∕L) are used to normalise � and at respectively. This is done so that the nor- malisation still shows the differences between the results across the dataset, since �L and �c are constant for all cases presented here. It is expected that other scales are better suited for such normalisation, and to justify the above phenomenological explanations. These will be explored in future work. Results showing the variation in the pdfs of mean curvature � and tangential strain rate at with changes in �0∕�L for fixed values of u�∕sL are presented in Fig. 6. As in Fig. 5 the mean curvature pdfs are shown in the left column while the tangential strain rate pdfs are shown in the right column. The pdfs corresponding to �0∕�L = 5.0, 2.0 and 1.25 are repre- sented by blue, red and yellow coloured lines respectively (see legend). The values of u�∕sL are 10 (top row), 20 (middle row) and 30 (bottom row). In each plot, the available mean kinetic energy of the turbulence is the same and hence changes in the pdf are due only to the different length scales of turbulence. It can be seen by comparison with Fig.  5 that the trends observed for decreasing �0∕�L are similar to those for increasing u�∕sL . The pdf of � (left column) for all values of u′ has the greatest peak value but smallest overall width for the highest value of �0∕�L = 5.0 (blue line). The pdfs become lower and wider as the value of �0∕�L decreases to 2.0 (red line) and to 1.25 (yellow line). As the value of �0∕�L decreases, the corresponding Kolmogorov length scale � also decreases, i.e. 𝜂 �0∕𝛿L=5.0 > ... > 𝜂 �0∕𝛿L=1.25  . This introduces an additional range of smaller length scales between 𝜂 �0∕𝛿L = 5.0 > l > 𝜂 �0∕𝛿L = 2.0 for the �0∕�L = 2.0 case. These additional small scales for each case correspond to higher � ( ∼ 1∕R ) and result in higher probability of high curvature events as �0 decreases. Hence, the width of the mean cur- vature pdfs increases. Likewise, since the cases at lower �0 contain fewer large scale eddies, the probability of low curvature events decreases, and thus the peak height of 112 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 the pdf decreases. The height and width of the pdfs for �0∕�L = 2.0 and �0∕�L = 1.25 remains roughly the same, again suggesting a tendency towards saturation in curvature for values of �0∕�L less than about 2.0. Note that the peak of the � pdfs occurs close to zero curvature for all values of �0. (a) Pdf of mean curvature (left) and tangential strain rate (right) at u′/sL = 10. (b) Pdfs of mean curvature (left) and tangential strain rate (right) at u′/sL = 20. (c) Pdfs of mean curvature (left) and tangential strain rate (right) at u′/sL = 30. Fig. 6   Pdfs of mean curvature (left column) and strain rate (right column) for different values of �0∕�L and fixed u�∕s L = 10 (top row), 20 (middle row) and 30 (bottom row) 113Flow, Turbulence and Combustion (2022) 109:101–123 1 3 Turning to the changes in the pdfs of tangential strain rate for different values of �0∕�L (Fig.  6, right column), for all values of u′ the peak value of each pdf is again the high- est for the largest integral length scale, i.e. �0∕�L = 5.0 . The pdfs become lower and their width increases as �0∕�L decreases, but the mean tangential strain rate remains positive for all cases. Also, the value of at at the peak becomes increasingly positive with decreas- ing �0∕�L . Again, the smaller values of Kolmogorov scale � for smaller values of �0∕�L result in greater positive straining of the flame due to their large velocity gradients. Once again, on comparing Figs. 5 and 6 it can be seen that for both � and at there is consider- able similarity in the changes observed in the pdfs in response to either increasing u�∕sL or decreasing �0∕�L . As before, proper scaling will be needed in future for greater clarity on the above explanations. The correlation between mean curvature � and tangential strain rate at can be observed from the joint pdfs of these quantities in the form of iso-probability curves shown in Fig. 7. The subfigures on the top row show the correlation between � and at for u�∕sL = 10, the middle row shows u�∕sL = 20 and bottom row shows u�∕sL = 30. The value of �0∕�L (a) Joint pdfs of mean curvature κ and tangential strain rate at at u′/sL = 10 and �0/δL = 5.0, 2.5 and 1.25 from left to right respectively. (b) Joint pdfs of mean curvature κ and tangential strain rate at at u′/sL = 20 and �0/δL = 5.0, 2.5 and 1.25 from left to right respectively. (c) Joint pdfs of mean curvature κ and tangential strain rate at at u′/sL = 30 and �0/δL = 5.0, 2.5 and 1.25 from left to right respectively. Fig. 7   joint pdfs of mean curvature and tangential strain rate for different values of �0∕�L (left to right) and for different values of u�∕s L (top to bottom) 114 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 decreases from 5.0 to 2.0 and then to 1.25 in the left, middle and right subfigures respec- tively in each row. In all cases, a negative correlation between � and at is observed. Similar plots have been presented in the past for increasing u′ (and fixed �0 ) (Nivarti and Cant 2019), and these showed a negative correlation between � and at which weakened as u′ was increased. This is a consequence of greater wrinkling of the flame due to higher energy in all scales of tur- bulent eddies. The current results also show the same trend with u′ at fixed �0 (Fig. 7, top to bottom in any column). By contrast, the strength of the �–at correlation tends to increase slightly at each fixed value of u′ as �0∕�L is decreased. This is consistent with the broaden- ing of the strain rate pdf (see Fig 6 right column) at fixed u′ , which indicates that the avail- able kinetic energy (for a fixed u′ ) is spread over a greater range of length scales giving a better match with the length scales of flame curvature. 4.3 � Displacement Speed Statistics It is interesting to consider the pdf of the displacement speed Sd and its dependence on flame stretch quantities. The pdfs of Sd are presented in Fig. 8a for different values of u�∕sL at fixed �0∕�L = 2.0 , and in Fig. 8b for different values of �0∕�L at fixed u�∕sL = 20. It can be seen that the displacement speed pdfs exhibit similar trends to those observed for the mean curvature and tangential strain rate pdfs. For increasing values of u�∕sL at fixed �0∕�L , the pdfs become shorter and wider, and the same pattern is observed for decreasing �0∕�L at fixed u�∕sL . In all cases, the mean value of Sd was found to be ∼ �sL∕�R where �R is the reactant density. The relatively strong correlation between Sd and flame stretch, i.e. both strain and curva- ture, ensures that the pdfs of Sd (Fig. 8) also become shorter and broader in a similar man- ner to the of � and at pdfs with either increasing u′ (Fig. 5) or decreasing �0 (Fig. 6). Joint pdfs of displacement speed Sd and mean curvature � in the form of iso-probability curves are shown in Fig. 9a (top row) while those for displacement speed and tangential strain rate at are shown in Fig. 9a (bottom row). For both rows, the value of �0∕�L = 2.0 is fixed while u�∕sL increases from 10 (left) to 20 (middle) to 30 (right). A distinct negative correlation can be seen between Sd and � (top row). The negative correlation is stronger for positive curvature than for negative curvature, as observed in previous studies (Chakraborty and Cant 2005b; Nivarti and Cant 2019) and explained in terms of the combined effects of both curvature and strain rate on the different components of displacements speed. As u′ increases, the strength of the correlation between Sd and � remains broadly unchanged. There is a positive correlation between Sd and at (bottom row), and this is consistent with the negative correlation observed between � and at (see Fig. 7). As u′ increases the underly- ing positive correlation between Sd and at weakens somewhat, and there is a tendency to favour Sd∕sL close to unity at low mean curvature for a broad range of strain rates, again as previously observed (Chakraborty and Cant 2005a; Nivarti and Cant 2019). Further joint pdfs of Sd and � are shown in Fig. 9b (top row) while joint pdfs of Sd and at are shown in Fig. 9b (bottom row), but this time the value of u�∕sL = 20 is held fixed while �0∕�L is varied from 5.0 (left) to 2.0 (middle) to 1.25 (right). Comparison of Fig. 9a and b shows that the trends for the Sd–� joint pdfs with decreasing �0 are similar to those observed for increasing u′ . This is less apparent for the Sd-at joint pdfs (bottom row). Here the tendency for Sd∕sL to cluster close to unity weakens as �0 decreases, due to the lower probability of finding low mean curvatures (see Fig. 5 left column). 115Flow, Turbulence and Combustion (2022) 109:101–123 1 3 The correlation coefficients corresponding to the joint pdfs presented above are summa- rised in Tables 3 and 4. First, correlation coefficients for increasing u′ at a fixed �0∕�L = 2.0 are shown in Table 3. Consistent with Fig. 7, the correlation coefficient for the �-at correla- tion becomes less negative, i.e. the correlation becomes weaker as the turbulence intensity increases and the same is true for the Sd-at correlation, while the Sd-� correlation coeffi- cient does not change much. Similar results were found for all values of �0∕�L. The corresponding correlation coefficients for a fixed value of u�∕sL at different values of �0∕�L are shown in Table 4. The Sd–at correlation coefficient is positive and increases (a) Pdf of displacement speed Sd for a range of values of u′/sL at fixed �0/δL = 2.0. (b) Pdf of displacement speed Sd for a range of values of �0/δL at fixed u′/sL = 20. Fig. 8   Probability density functions of displacement speed S d for a different values of u′ at fixed �0∕�L = 2.0 and b different values of �0∕�L at fixed u�∕s L = 20 116 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 slightly with decreasing �0∕�L . The Sd–� correlation coefficient, however, is negative and broadly constant for these cases. (a) Joint pdfs of displacement speedSd with κ (top row) and at (bottom row) for values of u′/sL = 10 (left), 20 (middle) and 30 (right) at fixed �0/δL = 2.0. (b) Joint pdfs of displacement speedSd with κ and at for a values of �0/δL = 5.0 (left), 2.0 (middle) and 1.25 (right) at fixed u′/sL = 20. Fig. 9   Joint pdfs of displacement speed S d with � and a t for a different values of u′ at fixed �0∕�L = 2.0 and b different values of �0∕�L at fixed u�∕s L = 20 Table 3   Correlation coefficients over a range of u�∕s L at fixed �0∕�L = 2.0 u �∕s L � − a t S d − a t S d − � �0∕�L = 2.0 10 − 0.732 0.514 − 0.454 20 − 0.553 0.390 − 0.525 30 − 0.450 0.241 − 0.530 117Flow, Turbulence and Combustion (2022) 109:101–123 1 3 The different components of displacement speed, i.e. the reaction component Sr , normal diffusion component Sn and tangential diffusion component St are defined by The joint pdfs for the reaction and normal diffusion components of Sd with � and at are pre- sented in Fig. 10 for fixed �0 and varying u′ . The first row shows the joint pdfs of the nor- mal diffusion component Sn with � while the second row shows the reaction component Sr with � . Similarly, the third row shows Sn with at while the fourth row shows Sr with at . The left, middle and right columns show the joint pdfs for u �∕sL = 10, 20 and 30 respectively. Joint pdfs are not shown for the tangential component St since this component is determin- istically related to � (see eqn. 13). (13)Sr = 𝜔̇ 𝜌𝜎 ; Sn = − nk 𝜌𝜎 𝜕 𝜕xk (𝜌D𝜎); St = −2D𝜅 Table 4   Correlation coefficients over a range of �0∕�L at fixed u �∕s L = 20 �0∕�L � − a t S d − a t S d − � u �∕s L = 20 5.0 − 0.476 0.271 − 0.489 2.0 − 0.553 0.390 − 0.525 1.25 − 0.565 0.379 − 0.552 κ− Sn κ− Sr at − Sn at − Sr Fig. 10   Joint pdfs of different components of the displacement speed with � and a t for different values of u �∕s L at fixed �0∕�L = 2.0 118 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 The relevant correlation coefficients are presented in Table  5 for a fixed value of �0∕�L = 2.0 and increasing u�∕sL . The correlation between St and at (fourth column) is of the same value but opposite sign to that between � and at (see Table 4) since St is directly proportional to negative � . This also results in a correlation coefficient between St and � equal to negative unity (right column). The normal component of displacement speed Sn shows a very weak correlation with at (third column) but a negative correlation with � (second column from the right). On the other hand, the reaction component Sr shows a positive correlation with at (second column) and a weak correlation with � (third column from the right). These results are consistent with the findings of several previous studies (Chakraborty and Cant 2005a, b; Nivarti and Cant 2019). Corresponding joint pdfs for the reaction and normal diffusion components of Sd with � and at are presented in Fig. 11 for fixed u�∕sL = 20 and varying �0∕�L . Again the first row Table 5   Correlation coefficients for the correlations between the reaction ( S r  ), normal diffusion ( S n  ) and tangential diffusion ( S t  ) components of displacement speed with mean curvature � and tangential strain rate a t over a range of u�∕s L at fixed �0∕�L = 2.0 u �∕s l S r − a t S n − a t S t − a t S r − � S n − � S t − � �0∕�L = 2.0 10 0.190 0.031 0.732 0.067 − 0.297 − 1.0 20 0.194 − 0.012 0.553 0.110 − 0.428 − 1.0 30 0.166 − 0.039 0.450 0.092 − 0.484 − 1.0 κ− Sn κ− Sr at − Sn at − Sr Fig. 11   Joint pdfs of different components of the displacement speed with � and a t for different values of �0∕�L at fixed u�∕s L = 20 119Flow, Turbulence and Combustion (2022) 109:101–123 1 3 shows normal diffusion component Sn with � while the second row Sr with � . The third row shows Sn with at while the fourth row shows Sr with at . The left, middle and right columns show the joint pdfs for �0∕�L = 5.0, 2.0 and 1.25 respectively. Comparison of Figs. 11 with 10 indicates a strong similarity between the joint pdfs for Sn and Sr with � (top two rows in each Figure), as well as similar trends for the plots of Sn and Sr with at (bottom two rows in each Figure). The correlation coefficients for the components of Sr , Sn and St with mean curvature and tangential strain rate are provided in Table 6 for fixed u�∕sL = 20 and a range of values of �0∕�L. Again, the correlation between St and at (fourth column) is the same in value but oppo- site in sign to that of � and at (see Fig. 5). There is a perfect negative correlation between St and � (right column), i.e. the correlation coefficient is equal to negative unity. As in the previous cases, Sn has a weak correlation with at (third column) while Sr has a positive correlation with at (second column). The correlation between Sn and � is negative (second column from the right) while the correlation between Sr and � is weakly positive (third column from the right). For the displacement speed components, the trends with decreasing �0 at fixed u′ are broadly similar to those observed for increasing u′ at fixed �0 . The correlation between Sn and at remains weak for all cases, and this is true also for the correlation between Sr and � . The negative correlation between Sn and � becomes stronger for variations in either u′ or �0 . The correlation between Sr and at is positive and weakens with increasing u′ , whereas for decreasing �0 the positive correlation strengthens slightly, consistent with the trends observed for the correlation between � and at . The correlation between St and at in all cases simply reflects the correlation between � and at and the deterministic dependence of St on �. 5 � Conclusions A set of Direct Numerical Simulations has been carried out in which the integral length scale �0∕�L as well as the turbulence intensity u�∕sL is varied in a systematic manner while keeping all other parameters constant. A twin flame configuration has been used in which the separation between the flames is kept sufficiently large so that the flames do not collide. Linear forcing of the turbulence has been applied between the two flames so that the turbu- lence does not decay rapidly as the flames propagate towards each other. The turbulent flame speed sT and flame area AT have been evaluated for all the cases. The turbulent flame speed is found to be greatest for the highest value of �0 at a fixed value of u′ , consistent with previous predictions (Peters 1999; Damköhler 1940) and DNS data (Nivarti et al. 2019). The ratio of sT∕sL to AT∕AL is found to be close to unity but increases Table 6   Correlation coefficients for the correlations between the reaction ( S r  ), normal diffusion ( S n  ) and tangential diffusion ( S t  ) components of displacement speed with mean curvature � and tangential strain rate a t over a range of �0∕�L at fixed u�∕s L = 20 �0∕�L S r − a t S n − a t S t − a t S r − � S n − � S t − � u �∕s L = 20 5.0 0.186 − 0.124 0.476 0.135 − 0.382 − 1.0 2.0 0.194 − 0.012 0.553 0.110 − 0.428 − 1.0 1.25 0.200 0.028 0.565 0.063 − 0.451 − 1.0 120 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 slightly as u�∕sL increases for fixed values of the length scale ratio �0∕�L . The flame sur- face density � has been evaluated for variations in both u′ and �0 . Changes observed for decreasing values of �0 at fixed u′ are found to be remarkably similar to those observed for increasing values of u′ at fixed �0. The probability density functions (pdfs) of the mean flame curvature and tangential strain rate have been analysed, again with the focus on varying �0∕�L for fixed u�∕sL . It has been found that as �0 decreases, the pdfs of both mean curvature and tangential strain rate become broader and flatter. The peak of the curvature pdfs shifts to a slightly more nega- tive value whereas that of the strain rate shifts to a more positive value. Analysis the displacement speed Sd indicates similar trends. The pdfs of Sd become broader and flatter with decreasing �0 at fixed u′ , just as they do with increasing u′ at fixed �0 . Correlations between displacement speed and mean curvature, as well as those between displacement speed and tangential strain rate, have been investigated using joint pdfs and correlation coefficients. The corresponding correlations for the reaction, normal diffusion and tangential diffusion components of displacement speed have been investigated also. The results for increasing u′ at fixed �0 are consistent with previous observations, and the variations observed for decreasing �0 at fixed u′ are found to be broadly similar. Future work will focus on extending the range of values of �0∕�L that can be explored. This will require larger computational domains and correspondingly larger computational resources. Proper scaling for the quantities presented in this paper and their dependence on the Reynolds number will be explored. It will also be necessary also to extend the study to incorporate any effects of chemistry. More detailed investigation of the physics of small- scale turbulence-flame interactions will be carried out, and the implications for modelling will be explored. Acknowledgements  The authors are grateful to EPSRC (Grant Number: EP/R029369/1) and ARCHER for financial and computational support as a part of their funding to the UK Consortium on Turbulent Reacting Flows (www.ukctrf.com). Special thanks are due to Dr. Girish Nivarti for valuable discussions on this topic. Declarations  Conflict of interest  The authors declare that they have no conflict of interest. Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Com- mons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://​creat​iveco​mmons.​org/​licen​ses/​by/4.​0/. References Aspden, A.J., Day, M.S., Bell, J.B.: Towards the distributed burning regime in turbulent premixed flames. J. Fluid Mech. 871, 1–21 (2019). https://​doi.​org/​10.​1017/​jfm.​2019.​316 Attili, A., Luca, S., Denker, D., Bisetti, F., Pitsch, H.: Turbulent flame speed and reaction layer thickening in premixed jet flames at constant Karlovitz and increasing Reynolds numbers. Proc. Combust. Inst. 38(2), 2939–2947 (2021). https://​doi.​org/​10.​1016/j.​proci.​2020.​06.​210 http://creativecommons.org/licenses/by/4.0/ https://doi.org/10.1017/jfm.2019.316 https://doi.org/10.1016/j.proci.2020.06.210 121Flow, Turbulence and Combustion (2022) 109:101–123 1 3 Batchelor, G.K., Townsend, A.A.: Decay of turbulence in the final period. Proc. R. Soc. Lond. A194, 527– 543 (1948) Boger, M., Veynante, D., Boughanem, H., Trouve, A.: Direct numerical simulation analysis of flame surface density concept for large eddy simulation of turbulent premixed combustion. In: twenty seventh sym- posium (international) on combustion 27, 917–925 (1998) Bradley, D.: How fast can we burn?. In: twenty-fourth symposium (international) on combustion, the com- bustion institute pp. 247–262 (1992) Bray, K.N.C., Cant, R.S.: Some applications of Kolmogorov’s turbulence research in the field of combus- tion. Proc. R. Soc. A 434, 217–240 (1991) Candel, S.M., Poinsot, T.J.: Flame stretch and the balance equation for the flame area. Combust. Sci. Tech- nol. 70, 1–15 (1990) Cant, R.S.: Senga2 user guide. Technical Report CUED-THERMO (2012) Cant, R.S., Pope, S.B., Bray, K.N.C.: Modelling of flamelet surface-to-volume ratio in turbulent premixed combustion. In: Twenty-third symposium (international) on combustion, the combustion institute pp. 809–815 (1990) Chakraborty, N., Cant, R.S.: Effects of strain rate and curvature on surface density function transport in tur- bulent premixed flames in the thin reaction zones regime. Phys. Fluids 17(6), 065108 (2005). https://​ doi.​org/​10.​1063/1.​19230​47 Chakraborty, N., Cant, R.S.: Influence of Lewis number on curvature effects in turbulent premixed flame propagation in the thin reaction zones regime. Phys. Fluids 17(10), 105105 (2005). https://​doi.​org/​10.​ 1063/1.​20842​31 Chaudhuri, S., Wu, F., Zhu, D., Law, C.K.: Flame speed and self-similar propagation of expanding turbu- lent premixed flames. Phys. Rev. Lett. 108, 044503 (2012). https://​doi.​org/​10.​1103/​PhysR​evLett.​108.​ 044503 Chen, J.H., Echekki, T., Kollman, W.: The mechanism of two-dimensional pocket formation in lean pre- mixed methane-air flames with implications to turbulent combustion. Combust. Flame 116, 15–48 (1999) Damköhler, G.: Der einfluss der turbulenz auf die flammengeschwindigkeit in gasgemischen. Z. Elektro- chem. Angew. Phys. Chem. 46(11), 601–626 (1940). https://​doi.​org/​10.​1002/​bbpc.​19400​461102 Driscoll, J.F.: Turbulent premixed combustion: flamelet structure and its effect on turbulent burning veloci- ties. Progress Energy Combust. Sci. 34, 91–134 (2008) Driscoll, J.F., Chen, J.H., Skiba, A.W., Carter, C.D., Hawkes, E.R., Wang, H.: Premixed flames subjected to extreme turbulence: some questions and recent answers. Progress Energy Combust. Sci. 76, 100802 (2020). https://​doi.​org/​10.​1016/j.​pecs.​2019.​100802 Hawkes, E.R., Chen, J.H.: Comparison of direct numerical simulation of lean premixed methane-air flames with strained laminar flame calculations. Combust. Flame 144, 112–125 (2006) Jenkins, K.W., Cant, R.S.: Direct Numerical Simulation of Turbulent Flame Kernels. In: Knight, D., Sakell, L. (eds.) Recent Advances in DNS and LES, pp. 191–202. Springer, Dordrecht (1999) Kennedy, C.A., Carpenter, M.H.: Additive Runge-Kutta schemes for convection-diffusion-reaction equa- tions. Appl. Numerical Math. 11(1), 139–181 (2003). https://​doi.​org/​10.​1016/​S0168-​9274(02)​00138-1 Kheirkhah, S., Gülder, L.: Turbulent premixed combustion in v-shaped flames: characteristics of flame front. Phys. Fluids 25(5), 055107 (2013). https://​doi.​org/​10.​1063/1.​48070​73 Kim, J., Satija, A., Lucht, R.P., Gore, J.P.: Effects of turbulent flow regime on perforated plate stabilized piloted lean premixed flames. Combust. Flame 211, 158–172 (2020). https://​doi.​org/​10.​1016/j.​combu​ stfla​me.​2019.​09.​027 Klein, M., Chakraborty, N., Ketterl, S.: A comparison of strategies for direct numerical simulation of tur- bulence chemistry interaction in generic planar turbulent premixed flames. Flow Turbul. Combust. 99, 955–971 (2017) Klimov, A.M.: Laminar flame in turbulent flow. Prikladnoy Mekhaniki i Tekhnicheskoy Fiziki Zhurnal (English translation AD-A200 241 Foreign Technology Division. Air Force Syst. Command 3, 49–58 (1963) Kulkarni, T., Bisetti, F.: Analysis of the development of the flame brush in turbulent premixed spherical flames. Combusti. Flame 234, 111640 (2021). https://​doi.​org/​10.​1016/j.​combu​stfla​me.​2021.​111640 Kulkarni, T., Bisetti, F.: Surface morphology and inner fractal cutoff scale of spherical turbulent premixed flames in decaying isotropic turbulence. Proc. Combust. Inst. 38(2), 2861–2868 (2021). https://​doi.​org/​ 10.​1016/j.​proci.​2020.​06.​117 Kulkarni, T., Buttay, R., Kasbaoui, M.H., Attili, A., Bisetti, F.: Reynolds number scaling of burning rates in spherical turbulent premixed flames. J. Fluid Mech. 906, 25–89 (2021). https://​doi.​org/​10.​1017/​jfm.​ 2020.​784 Law, C.K.: Dynamics of stretched flames. Proc. Combust. Inst. 22, 1381–1402 (1988) https://doi.org/10.1063/1.1923047 https://doi.org/10.1063/1.1923047 https://doi.org/10.1063/1.2084231 https://doi.org/10.1063/1.2084231 https://doi.org/10.1103/PhysRevLett.108.044503 https://doi.org/10.1103/PhysRevLett.108.044503 https://doi.org/10.1002/bbpc.19400461102 https://doi.org/10.1016/j.pecs.2019.100802 https://doi.org/10.1016/S0168-9274(02)00138-1 https://doi.org/10.1063/1.4807073 https://doi.org/10.1016/j.combustflame.2019.09.027 https://doi.org/10.1016/j.combustflame.2019.09.027 https://doi.org/10.1016/j.combustflame.2021.111640 https://doi.org/10.1016/j.proci.2020.06.117 https://doi.org/10.1016/j.proci.2020.06.117 https://doi.org/10.1017/jfm.2020.784 https://doi.org/10.1017/jfm.2020.784 122 Flow, Turbulence and Combustion (2022) 109:101–123 1 3 Lipatnikov, A., Chomiak, J.: Effects of turbulence length scale on flame speed: a modelling study. In: Rodi, W., Laurence, D. (eds.) Engineering Turbulence Modelling and Experiments, vol. 4, pp. 841–850. Elsevier Science Ltd, Oxford (1999). https://​doi.​org/​10.​1016/​B978-​00804​3328-8/​50081-3 Lipatnikov, A., Chomiak, J.: Turbulent flame speed and thickness: phenomenology, evaluation, and applica- tion in multi-dimensional simulations. Progress Energy Combust. Sci. 28, 1–74 (2002) Lipatnikov, A.N., Chomiak, J., Sabelnikov, V.A., Nishiki, S., Hasegawa, T.: Unburned mixture fingers in premixed turbulent flames. Proc. Combust. Inst. 35(2), 1401–1408 (2015). https://​doi.​org/​10.​1016/j.​ proci.​2014.​06.​081 Lipatnikov, A.N., Chomiak, J., Sabelnikov, V.A., Nishiki, S., Hasegawa, T.: A DNS study of the physical mechanisms associated with density ratio influence on turbulent burning velocity in premixed flames. Combust. Theory Model. 22(1), 131–155 (2018). https://​doi.​org/​10.​1080/​13647​830.​2017.​13902​65 Liu, Z., Unni, V.R., Chaudhuri, S., Law, C.K., Saha, A.: Local statistics of laminar expanding flames sub- jected to Darrieus–Landau instability. Proc. Combust. Inst. 38(2), 1993–2000 (2021). https://​doi.​org/​ 10.​1016/j.​proci.​2020.​06.​118 Luca, S., Attili, A., Lo Schiavo, E., Creta, F.: On the statistics of flame stretch in turbulent premixed jet flames in the thin reaction zone regime at varying Reynolds number. Proc. Combust. Inst. 37(2), 2451– 2459 (2019). https://​doi.​org/​10.​1016/j.​proci.​2018.​06.​194 Lundgren, T.S.: Linearly Forced Isotropic Turbulence. Annual Research Briefs, pp. 461–473. CTR, Stanford (2003) Mohammadnejad, S., An, Q., Vena, P., Yun, S., Kheirkhah, S.: Contributions of flame thickening and extinctions to a heat release rate marker of intensely turbulent premixed hydrogen-enriched methane- air flames. Combust. Flame 231, 111481 (2021). https://​doi.​org/​10.​1016/j.​combu​stfla​me.​2021.​111481 Nivarti, G.V., Cant, R.S.: Direct numerical simulation of the bending effect in turbulent premixed flames. Proc. Combust. Inst. 36, 1903–1910 (2017) Nivarti, G.V., Cant, R.S.: Stretch rate and displacement speed correlations for Increasingly–Turbu- lent premixed flames. Flow Turbul. Combust. 102, 957–971 (2019). https://​doi.​org/​10.​1007/​ s10494-​018-​9990-7 Nivarti, G.V., Cant, R.S., Hochgreb, S.: Reconciling turbulent burning velocity with flame surface area in small-scale turbulence. J. Fluid Mech. 858, 25–87 (2019) Orszag, S.A.: Numerical methods for the simulation of turbulence. Phys. Fluids Suppl. 2, 250–257 (1972) Peters, N.: The turbulent burning velocity for large-scale and small-scale turbulence. J. Fluid Mech. 384, 107–132 (1999) Poinsot, T., Lele, S.: Boundary conditions for direct simulations of compressible viscous flows. J. Comput. Phys. 101, 104–129 (1992) Poludnenko, A.Y., Oran, E.S.: The interaction of high-speed turbulence with flames: global properties and internal flame structure. Combust. Flame 157, 995–1011 (2010) Poludnenko, A.Y., Oran, E.S.: The interaction of high-speed turbulence with flames: turbulent flame speed. Combust. Flame 158, 301–326 (2011) Pope, S.B.: The evolution of surfaces in turbulence. Int. J. Eng. Sci. 26, 445–469 (1988) Rosales, C., Meneveau, C.: Linear forcing in numerical simulations of isotropic turbulence: physical space implementations and convergence properties. Phys. Fluids 17, 095106 (2005) Sabelnikov, V.A., Yu, R., Lipatnikov, A.N.: Thin reaction zones in constant-density turbulent flows at low Damköhler numbers: theory and simulations. Phys. Fluids 31(5), 055104 (2019). https://​doi.​org/​10.​ 1063/1.​50901​92 Shy, S., Liu, C., Lin, J., Chen, L., Lipatnikov, A., Yang, S.: Correlations of high-pressure lean methane and syngas turbulent burning velocities: effects of turbulent Reynolds, Damköhler, and Karlovitz numbers. Proc. Combust. Inst. 35(2), 1509–1516 (2015). https://​doi.​org/​10.​1016/j.​proci.​2014.​07.​026 Skiba, A.W., Carter, C.D., Hammack, S.D., Driscoll, J.F.: High-fidelity flame-front wrinkling measurements derived from fractal analysis of turbulent premixed flames with large Reynolds numbers. Proc. Com- bust. Insti. 38(2), 2809–2816 (2021). https://​doi.​org/​10.​1016/j.​proci.​2020.​06.​041 Skiba, A.W., Wabel, T.M., Carter, C.D., Hammack, S.D., Temme, J.E., Lee, T., Driscoll, J.F.: Reaction layer visualization: a comparison of two plif techniques and advantages of khz-imaging. Proc. Combust. Inst. 36, 4593–4601 (2017) Smith, K.O., Gouldin, F.C.: Turbulence effects on flame speed and flame structure. AIAA J. 17(11), 1243– 1250 (1979). https://​doi.​org/​10.​2514/3.​61305 Song, W., Hernández Pérez, F.E., Tingas, E., Im, H.G.: Statistics of local and global flame speed and struc- ture for highly turbulent H 2/air premixed flames. Combust. Flame 232, 111523 (2021). https://​doi.​org/​ 10.​1016/j.​combu​stfla​me.​2021.​111523 https://doi.org/10.1016/B978-008043328-8/50081-3 https://doi.org/10.1016/j.proci.2014.06.081 https://doi.org/10.1016/j.proci.2014.06.081 https://doi.org/10.1080/13647830.2017.1390265 https://doi.org/10.1016/j.proci.2020.06.118 https://doi.org/10.1016/j.proci.2020.06.118 https://doi.org/10.1016/j.proci.2018.06.194 https://doi.org/10.1016/j.combustflame.2021.111481 https://doi.org/10.1007/s10494-018-9990-7 https://doi.org/10.1007/s10494-018-9990-7 https://doi.org/10.1063/1.5090192 https://doi.org/10.1063/1.5090192 https://doi.org/10.1016/j.proci.2014.07.026 https://doi.org/10.1016/j.proci.2020.06.041 https://doi.org/10.2514/3.61305 https://doi.org/10.1016/j.combustflame.2021.111523 https://doi.org/10.1016/j.combustflame.2021.111523 123Flow, Turbulence and Combustion (2022) 109:101–123 1 3 Varma, A.R., Ahmed, U., Klein, M., Chakraborty, N.: Effects of turbulent length scale on the bending effect of turbulent burning velocity in premixed turbulent combustion. Combust. Flame 233, 111569 (2021). https://​doi.​org/​10.​1016/j.​combu​stfla​me.​2021.​111569 Williams, F.A.: Combustion Theory. Westview Press, Boulder (1985) Yu, R., Lipatnikov, A.N.: Direct numerical simulation study of statistically stationary propagation of a reac- tion wave in homogeneous turbulence. Phys. Rev. E 95, 063101 (2017). https://​doi.​org/​10.​1103/​PhysR​ evE.​95.​063101 Yu, R., Lipatnikov, A.N.: DNS study of dependence of bulk consumption velocity in a constant-density reacting flow on turbulence and mixture characteristics. Phys. Fluids 29(6), 065116 (2017). https://​doi.​ org/​10.​1063/1.​49908​36 https://doi.org/10.1016/j.combustflame.2021.111569 https://doi.org/10.1103/PhysRevE.95.063101 https://doi.org/10.1103/PhysRevE.95.063101 https://doi.org/10.1063/1.4990836 https://doi.org/10.1063/1.4990836 Turbulence Intensity and Length Scale Effects on Premixed Turbulent Flame Propagation Abstract 1 Introduction 2 Theory 3 DNS Dataset for the Current Study 4 Results 4.1 Turbulent Flame Speed 4.2 Curvature and Strain Rate 4.3 Displacement Speed Statistics 5 Conclusions Acknowledgements References