Submit manuscript...
eISSN: 2576-4500

Aeronautics and Aerospace Open Access Journal

Research Article Volume 2 Issue 5

Assessment of deconvolution-based flamelet methods for progress variable rate modeling

ZM Nikolaou,1 L Vervisch2

1Computation-Based Science and Technology Research Centre (CaSToRC), The Cyprus Institute, Cyprus
2CORIA-CNRS, Normandie Universite, INSA de Rouen Normandie, Saint-Etiennedu-Rouvray, France

Correspondence: ZM Nikolaou, Computation-Based Science and Technology Research Centre (CaSToRC), The Cyprus Institute, Nicosia, 2121, Cyprus

Received: August 27, 2018 | Published: September 21, 2018

Citation: Nikolaou ZM, Vervisch L. Assessment of deconvolution-based flamelet methods for progress variable rate modeling. Aeron Aero Open Access J. 2018;2(5):274-281. DOI: 10.15406/aaoaj.2018.02.00061

Download PDF

Abstract

A novel approach for modeling the progress variable reaction rate in Large Eddy Simulations of turbulent and reacting flows is proposed. This is done in the context of two popular flamelet models which require the progress variable variance as input. The approach is based on using a recently proposed deconvolution method for modeling the variance. The deconvolution-modeled variance is used as an input in the flamelet models for modeling the filtered progress variable rate. The assessment of the proposed approach is conducted a priori using direct numerical simulation data of turbulent premixed flames. For the conditions tested in this study, deconvolution does not introduce a significant bias in the flamelet models' predictions, while a quantitatively good prediction of the progress variable rate is obtained for both flamelet models considered.

Introduction

Large Eddy Simulation (LES) is becoming the de facto approach both in industry and academia for modeling a wide range of flows. This includes a wide range of devices with complex geometry in the aerospace, energy, marine, and automotive industries to name but a few. The numerical modeling however of turbulent and reacting flows is computationally challenging and expensive, particularly in the context of LES. Reacting flows, and specifically combustion processes, include a wide range of length and time scales, all of which need to be adequately resolved during the simulation, on a mesh which for practical geometries can include millions of cells.1 An accurate description of the chemical processes involved in such flows, requires the use of fairly large detailed chemical kinetics mechanisms with hundreds of species and reactions which poses a significant computational workload. In an effort to reduce the dimensionality of the problem, various combustion modeling strategies have been developed throughout the years. The earliest modeling attempts employed simple, 1-step chemistry, for modeling the chemical kinetics.2 This, reduces the computational load substantially. However, with an ever increasing need for accuracy but at the same time detailed chemistry information (e.g. pollutants concentrations, CO2, NOx etc.), more elaborate combustion models were developed able to describe complex chemistry effects. These include, among others, Thickened Flames (TF), Flamelet Methods (FM), Conditional Moment Closure (CMC), transported-pdf methods etc. Each method has both its merits and drawbacks, and an in-depth discussion on the subject can be found in the reviews.3‒5

Flamelet methods in particular are popular since they are relatively straightforward to implement, and have been applied in a number of LES studies both for premixed and for non-premixed flames, with overall good results.6‒16 These studies include ows in complex geometries and for different flow configurations (freely-propagating, recirculating, etc.). The filtered progress variable c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ which is a main LES solution variable, and it's variance σ 2 = c ˜ 2 c ˜ c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHdp WCjuaGdaahaaWcbeqaaKqzadGaaGOmaaaajugibiabg2da9iqadoga gaacaKqbaoaaCaaaleqabaqcLbmacaaIYaaaaKqzGeGaeyOeI0Iabm 4yayaaiaGabm4yayaaiaaaaa@4388@ as obtained from a suitable model, are used as parameters in prebuilt tables for obtaining important information and for providing closures for un-closed terms in the governing equations. This includes species mean (filtered) mass fractions, mean (filtered) reaction rates etc. Such methods offer a simple, yet effective way of providing detailed chemistry information (a posteriori), while keeping the computational requirements to a minimum. A key unclosed term in the transport equation for the progress variable, c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ , is the filtered progress variable reaction rate, w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ , where the overbar denotes a spatial filtering operation as defined in the context of  LES. This is a dominant term, particularly in the reaction zone of the flame,17 and an accurate model is required for this term. Two popular amelet methods for modeling w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ , and which involve the progress variable variance, are the presumed-pdf approach and the Filtered Laminar Flame (FLF) approach.7,13 In the presumed-pdf approach, the progress variable c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ and it's variance σ2, are used to parameterise the progress variable pdf p ˜ ( ζ; c ˜ , σ 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGWb GbaGaajuaGdaqadaqaaKqzGeGaeqOTdONaai4oaiqadogagaacaiaa cYcacqaHdpWCjuaGdaahaaqabeaajugWaiaaikdaaaaajuaGcaGLOa Gaayzkaaaaaa@433D@ , where ζ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaH2o GEaaa@3842@ is the sample-space variable for c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ . This pdf is usually taken to be a β-function. A usually steady, 1D, laminar flame solution is used for integrating with the pdf thus obtaining a closure for the rate.

In the FLF approach, a 1D laminar flame solution is pre-filtered, and a table is constructed which contains values of the filtered reaction rate (or other variables of interest). This is done for discrete filter-width values, and the table contains the filtered values of the variables at all spatial positions across the 1D flame. The progress variable as obtained from the LES, and its variance, are used as parameters for accessing the table, and obtaining a value of the filtered rate. It is clear that in both methods the progress variable variance is an important parameter for determining the filtered rate. In turn, the value of the rate affects the flame propagation characteristics (flame displacement speed), and the interaction of the flame with the flow field.17 In the traditional modeling approach, an algebraic model is used or a transport equation is solved for in order to obtain an estimate of the variance, and there exist many studies in the literature using these two classic modeling approaches.18‒24 Algebraic based models are generally computationally efficient, whereas an additional transport equation requires solving in the case of transport equation-based models. Both of these approaches however, usually include additional unknown model parameters whose value significantly affects the predictive ability of the model.

 In classic, gradient-based methods for example, a filter-size dependent constant is involved.20 In transport-equation based methods, additional un-closed terms appear in the transport equation for the variance which also require modeling-an important un-closed term is the scalar dissipation rate. In addition, the value of the models’ parameters depends on many factors such as the flow configuration, Reynolds number Re, Damkohler number Da etc. This fact, may limit the generality of classic modeling approaches, and as a result the applicability range of flamelet methods which require the variance as input. An alternative modeling approach in the context of LES is based on deconvolution. Deconvolution-based methods aim to recover a good approximation ϕ * MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHvp GzjuaGdaahaaqabeaajugWaiaacQcaaaaaaa@3AD9@ of the un-filtered field ϕ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHvp Gzaaa@384D@ from knowledge of the filtered field ϕ ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacuaHvp Gzgaqeaaaa@3865@ on the LES mesh. Deconvolution was recently found to be a powerful tool to complement physical modeling by directly exploiting the information contained in the resolved signals.25, 26 In reacting ows, deconvolution methods were used in,27,28 and more recently in29‒31 in order to recover variables such as the temperature, with overall promising results in both a priori and a posteriori studies.

In recent work,32,33 an iterative deconvolution method was proposed which was shown to provide quantitatively good predictions of the un-filtered progress variable. Furthermore, it was shown that Iterative Deconvolution combined with Explicit Filtering (IDEF) provides quantitatively good predictions of the progress variable variance and the scalar flux.32,33 Given that diffusive terms in the transport equation for the progress variable can be adequately modeled with standard practices, the only remaining term is the filtered progress variable reaction rate. A hybrid deconvolution-flamelet formulation is thus proposed for the above two flamelet approaches, for modeling the filtered rate. The main advantages of such an approach are twofold:

  1. Detailed chemistry information is provided a posteriori while still solving a single transport equation for the progress variable.
  2. A parameter-free method, namely deconvolution, which is not tied to a specific flow regime or combustion mode, is used for obtaining the variance.

These two advantages combined, have the potential to render a more robust flamelet strategy for modeling the reaction rate in turbulent premixed flames in the context of LES. Since the deconvolution method described in32, 33 is local i.e. the method provides a value of the variance at every mesh point, it is still unclear what the overall effect of modeling the progress variable reaction rate will be when using a flamelet method such as a presumed-pdf approach, or FLF. The aim of this study is to evaluate the sensitivity in modeling the filtered rate using these two approaches, where IDEF is employed for modeling the variance. In contrast to most a priori assessments which are conducted on the fine DNS mesh, the assessment in this study is conducted on a much coarser LES mesh in order to simulate an actual LES. Even though a priori assessments in general using DNS data do not guarantee functionality of the model in actual LES, the results in32,33 indicate that the simulated LES a priori assessment procedure is more suitable for evaluating LES sub-grid models. This is particularly true for sub-grid scale models which include gradients of the resolved fields.32 In addition, such a priori assessments are useful in clearly evaluating the performance of a combustion model, free from the influences of turbulence models and of the often dissipative numerical schemes which are used in actual LES. In the text which follows, section 2 lists details of the DNS database used for the assessment, section 3 describes the hybrid deconvolution-flamelet approach used, and results are presented/discussed in section 4.

Description of the DNS database

Direct simulations were conducted using the SENGA2 code.34 The code solves the compressible Navier-Stokes equations using a 10th order finite difference scheme for interior points, and a 4th order Runge-Kutta scheme for the time-stepping,

ρ t + ρ u k x k =0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaSaaae aacqGHciITcqaHbpGCaeaacqGHciITcaWG0baaaiabgUcaRmaalaaa baGaeyOaIyRaeqyWdiNaamyDamaaBaaabaqcLbmacaWGRbaajuaGbe aaaeaacqGHciITcaWG4bWaaSbaaeaajugWaiaadUgaaKqbagqaaaaa cqGH9aqpcaaIWaaaaa@4AE8@ (1)

ρ u i t + ρ u k u i x k = ρ x i + T ki x k MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaSaaae aacqGHciITcqaHbpGCcaWG1bWaaSbaaeaajugWaiaadMgaaKqbagqa aaqaaiabgkGi2kaadshaaaGaey4kaSYaaSaaaeaacqGHciITcqaHbp GCcaWG1bWaaSbaaeaajugWaiaadUgaaKqbagqaaiaadwhadaWgaaqa aKqzadGaamyAaaqcfayabaaabaGaeyOaIyRaamiEamaaBaaabaqcLb macaWGRbaajuaGbeaaaaGaeyypa0ZaaSaaaeaacqGHciITcqaHbpGC aeaacqGHciITcaWG4bWaaSbaaeaajugWaiaadMgaaKqbagqaaaaacq GHRaWkdaWcaaqaaiabgkGi2kaadsfadaWgaaqaaKqzadGaam4Aaiaa dMgaaKqbagqaaaqaaiabgkGi2kaadIhadaWgaaqaaKqzadGaam4Aaa qcfayabaaaaaaa@6638@ (2)

ρE t + ρ u k E x k = ρ u k x k q k x k + T km u m x k MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaSaaae aacqGHciITcqaHbpGCcaWGfbaabaGaeyOaIyRaamiDaaaacqGHRaWk daWcaaqaaiabgkGi2kabeg8aYjaadwhadaWgaaqaaKqzadGaam4Aaa qcfayabaGaamyraaqaaiabgkGi2kaadIhadaWgaaqaaKqzadGaam4A aaqcfayabaaaaiabg2da9iabgkHiTmaalaaabaGaeyOaIyRaeqyWdi NaamyDamaaBaaabaqcLbmacaWGRbaajuaGbeaaaeaacqGHciITcaWG 4bWaaSbaaeaajugWaiaadUgaaKqbagqaaaaacqGHsisldaWcaaqaai abgkGi2kaadghadaWgaaqaaKqzadGaam4AaaqcfayabaaabaGaeyOa IyRaamiEamaaBaaabaqcLbmacaWGRbaajuaGbeaaaaGaey4kaSYaaS aaaeaacqGHciITcaWGubWaaSbaaeaajugWaiaadUgacaWGTbaajuaG beaacaWG1bWaaSbaaeaajugWaiaad2gaaKqbagqaaaqaaiabgkGi2k aadIhadaWgaaqaaKqzadGaam4Aaaqcfayabaaaaaaa@741B@ (3)

ρ Y α t + ρ u k Y α x k = w ˙ α ρ V α,k Y α x k MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaSaaae aacqGHciITcqaHbpGCcaWGzbWaaSbaaeaajugWaiabeg7aHbqcfaya baaabaGaeyOaIyRaamiDaaaacqGHRaWkdaWcaaqaaiabgkGi2kabeg 8aYjaadwhadaWgaaqaaKqzadGaam4AaaqcfayabaGaamywamaaBaaa baqcLbmacqaHXoqyaKqbagqaaaqaaiabgkGi2kaadIhadaWgaaqaaK qzadGaam4Aaaqcfayabaaaaiabg2da9iqadEhagaGaamaaBaaabaqc LbmacqaHXoqyaKqbagqaaiabgkHiTmaalaaabaGaeyOaIyRaeqyWdi NaamOvamaaBaaabaqcLbmacqaHXoqycaGGSaGaam4AaaqcfayabaGa amywamaaBaaabaqcLbmacqaHXoqyaKqbagqaaaqaaiabgkGi2kaadI hadaWgaaqaaKqzadGaam4Aaaqcfayabaaaaaaa@6B07@ (4)

Where the symbols have their usual meanings. Table 1 lists details of the DNS database used to conduct the analysis. This is a freely-propagating premixed multi-component fuel-air flame. A detailed chemical mechanism which was developed specifically for such fuels was used with 49 reactions and 15 species.35 urms is the rms value of the fluctuating component of the incoming velocity field, with integral length scale lint on the reactant side. The turbulence Reynolds number is Re= u rms l int v r MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaciOuai aacwgacqGH9aqpdaWccaqaaiaadwhadaWgaaqaaKqzadGaamOCaiaa d2gacaWGZbaajuaGbeaacaWGSbWaaSbaaeaajugWaiGacMgacaGGUb GaaiiDaaqcfayabaaabaGaamODamaaBaaabaqcLbmacaWGYbaajuaG beaaaaaaaa@488B@ , the Damkohler number is Da= ( l int u rms ) ( δ S L ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamirai aadggacqGH9aqpkmaaliaabaWaaeWaaeaajuaGdaWccaqaaiaadYga daWgaaqaaKqzadGaciyAaiaac6gacaGG0baajuaGbeaaaeaacaWG1b WaaSbaaeaajugWaiaadkhacaWGTbGaam4CaaqcfayabaaaaaGccaGL OaGaayzkaaaabaWaaeWaaeaajuaGdaWccaGcbaqcLbsacqaH0oazaO qaaKqzGeGaam4uaKqbaoaaBaaaleaajugWaiaadYeaaSqabaaaaaGc caGLOaGaayzkaaaaaaaa@4F8D@ and the Karlovitz number is Ka= ( δ η k ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacaWGlb GaamyyaOGaeyypa0ZaaeWaaeaajuaGdaWccaGcbaqcLbsacqaH0oaz aOqaaKqzGeGaeq4TdGwcfa4aaSbaaSqaaKqzadGaam4AaaWcbeaaaa aakiaawIcacaGLPaaadaahaaWcbeqaaKqzadGaaGOmaaaaaaa@44FB@ , where sL is the laminar flame speed, and the thickness δ=( v r S L ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaH0o azcqGH9aqpjuaGdaqadaGcbaqcfa4aaSGaaOqaaKqzGeGaamODaKqb aoaaBaaaleaajugWaiaadkhaaSqabaaakeaajugibiaadofajuaGda WgaaWcbaqcLbmacaWGmbaaleqaaaaaaOGaayjkaiaawMcaaaaa@44AE@ . Note that the laminar flame thickness is defined as δ L = ( T p T r ) max( dT dx ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaH0o azjuaGdaWgaaqaaKqzadGaamitaaqcfayabaqcLbsacqGH9aqpjuaG daWccaqaamaabmaabaqcLbsacaWGubqcfa4aaSbaaeaajugWaiaadc haaKqbagqaaKqzGeGaeyOeI0IaamivaKqbaoaaBaaabaqcLbmacaWG YbaajuaGbeaaaiaawIcacaGLPaaaaeaajugibiGac2gacaGGHbGaai iEaKqbaoaabmaabaWaaSGaaeaajugibiaadsgacaWGubaajuaGbaqc LbsacaWGKbGaamiEaaaaaKqbakaawIcacaGLPaaaaaaaaa@5511@ where Tr, Tp are the reactant and product temperatures respectively. These conditions place the flame in the distributed/broken reaction zones regime according to the classic combustion diagram by Peters.36 Further details of the simulations are given in.37,38 In the sections which follow, the progress variable c is based on temperature.

Case

urms/sl

lint

Ret

Da

Ka

A

3.18

16.54

52.66

5.19

1.39

B

9

16.66

150.05

1.85

6.62

C

14.04

16.43

230.69

1.17

12.97

Table 1 Turbulent flame parameters for the DNS studies

Mathematical background

Filtering and deconvolution operations

In the context of LES, filtering is a convolution operation between the signal ϕ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHvp Gzaaa@384D@ and the filter G,

ϕ ¯ ( x _ ,t )= x _ = G( x _ x _ ) ϕ( x _ ,t )d x _ =G*ϕ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacuaHvp GzgaqeaKqbaoaabmaabaqcLbsaceWG4bGba0bacaGGSaGaamiDaaqc faOaayjkaiaawMcaaKqzGeGaeyypa0tcfa4aa8qCaeaajugibiaadE eajuaGdaqadaqaaKqzGeGabmiEayaaDaGaeyOeI0IabmiEayaafyaa DaaajuaGcaGLOaGaayzkaaaabaqcLbmaceWG4bGbauGba0bacqGH9a qpcqGHsislcqGHEisPaKqbagaajugWaiabg6HiLcqcLbsacqGHRiI8 aiabew9aMLqbaoaabmaabaqcLbsaceWG4bGbauGba0bacaGGSaGaam iDaaqcfaOaayjkaiaawMcaaKqzGeGaamizaiqadIhagaqbgaqhaiab g2da9iaadEeacaGGQaGaeqy1dygaaa@62D7@ (5)

A Gaussian filter is used in the present study,

G( x _ )= ( 6 π Δ 2 ) 3 2 e 6 x _ x _ Δ 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacaWGhb qcfa4aaeWaaeaajugibiqadIhagaqhaaqcfaOaayjkaiaawMcaaKqz GeGaeyypa0tcfa4aaeWaaeaadaWcaaqaaKqzGeGaaGOnaaqcfayaaK qzGeGaeqiWdaNaeuiLdqucfa4aaWbaaeqabaqcLbmacaaIYaaaaaaa aKqbakaawIcacaGLPaaadaahaaqabeaalmaalaaajuaGbaqcLbmaca aIZaaajuaGbaqcLbmacaaIYaaaaaaajugibiaadwgajuaGdaahaaWc beqaamaalaaabaqcLbmacqGHsislcaaI2aGabmiEayaaDaGaeyyXIC TabmiEayaaDaaaleaajugWaiabfs5aeTWaaWbaaWqabeaajugWaiaa ikdaaaaaaaaaaaa@5BB4@ (6)

Where Δ is the filter width. Favre-filtered variables are defined as,

ϕ ¯ ( x _ ,t )= ρϕ ¯ ρ ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacuaHvp GzgaqeaKqbaoaabmaabaqcLbsaceWG4bGba0bacaGGSaGaamiDaaqc faOaayjkaiaawMcaaKqzGeGaeyypa0tcfa4aaSaaaeaadaqdaaqaaK qzGeGaeqyWdiNaeqy1dygaaaqcfayaaKqzGeGafqyWdiNbaebaaaaa aa@47B3@ (7)

 

The DNS data ρ and ρc are first filtered using Eq. 5 on the fine DNS mesh. This is done for filter widths Δ+=Δ/δL of 1.0, 2.0 and 3.0 respectively. In order to simulate an actual LES, the filtered data as obtained on the DNS mesh are then sampled onto a much coarser LES mesh. This is done using high-order Lagrange polynomials for the interpolation. The filtered variables as obtained by sampling on the LES mesh are then ρ ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacuaHbp GCgaqeaaaa@385D@

and ρc ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aa0aaae aacqaHbpGCcaWGJbaaaaaa@393D@ . The size of the LES mesh is determined based on the criterion described in32: the deconvolution method requires the filtered fields to be adequately resolved on the LES mesh, and the criterion derived in32 ensures this. This requirement is translated to an LES mesh spacing ratio h/Δ=0.25 where h is the LES mesh spacing. Table 2 and Table 3 list the DNS and LES meshes for each of the three DNS databases used in this study. In comparison to the Kolmogorov length scale ηk, the LES mesh is much coarser: h/ηk for cases A, B and C for the finest LES mesh (Δ+=1.0) is 6.8, 14.9 and 20.7 respectively, and much larger for the coarsest LES mesh considered in this study. Therefore, small-scale information of the order of ηk is lost on the simulated LES mesh.

+

Nx

Ny

Nz

h/lT

DNS

768

384

384

0.03

1

74

37

37

0.35

2

37

18

18

0.69

3

24

12

12

1.04

Table 2 DNS and LES meshes for cases A and B with h/Δ=0.25

+

Nx

Ny

Nz

h/ lT

DNS

1632

544

544

0.02

1

112

37

37

0.21

2

56

18

18

0.42

3

37

12

12

0.64

Table 3 DNS and LES meshes for case C with h/Δ=0.25

The filtered variables ρ ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacuaHbp GCgaqeaaaa@385D@ and ρc ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aa0aaae aacqaHbpGCcaWGJbaaaaaa@393D@ as obtained by sampling on the LES mesh are deconvoluted using an iterative algorithm in order to obtain estimates of the un-filtered fields namely ρ* and (ρc)* (on the LES mesh). The algorithm (Van-Cittert) reads,39‒41

ϕ *n+1 = ϕ *n +b( ϕ ¯ G* ϕ *n ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHvp GzjuaGdaahaaWcbeqaaKqzadGaaiOkaiaad6gacqGHRaWkcaaIXaaa aKqzGeGaeyypa0Jaeqy1dywcfa4aaWbaaSqabeaajugWaiaacQcaca WGUbaaaKqzGeGaey4kaSIaamOyaKqbaoaabmaakeaajugibiqbew9a MzaaraGaeyOeI0Iaam4raiaacQcacqaHvpGzjuaGdaahaaWcbeqaaK qzadGaaiOkaiaad6gaaaaakiaawIcacaGLPaaaaaa@5306@ (8)

With ϕ *0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHvp GzjuaGdaahaaqabeaajugWaiaacQcacaaIWaaaaaaa@3B93@ = ϕ ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacuaHvp Gzgaqeaaaa@3865@ , and the progress variable variance is calculated by explicitly filtering the deconvoluted fields as obtained on the LES mesh,

σ 2 IDEF = ( ρc )*( ρc )*/ ρ* ¯ ρ ¯ c ˜ c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHdp WCjuaGdaahaaqabeaajugWaiaaikdaaaqcfa4aaSbaaeaajugWaiaa dMeacaWGebGaamyraiaadAeaaKqbagqaaKqzGeGaeyypa0tcfa4aaS aaaeaadaqdaaqaamaalyaabaWaaeWaaeaajugibiabeg8aYjaadoga aKqbakaawIcacaGLPaaajugibiaacQcajuaGdaqadaqaaKqzGeGaeq yWdiNaam4yaaqcfaOaayjkaiaawMcaaKqzGeGaaiOkaaqcfayaaKqz GeGaeqyWdiNaaiOkaaaaaaaajuaGbaqcLbsacuaHbpGCgaqeaaaacq GHsislceWGJbGbaGaaceWGJbGbaGaaaaa@59D2@ (9)

In practice, the deconvolution algorithm is implemented with an error controller in order to ensure that the deconvolved fields lie within correct physical limits. Further details of the exact implementation of the method are given in.32,33 The estimate of the variance as given by Eq. 9 is then used for parameterising the reaction rate models which are presented in the next two sections.

Unstrained flamelet presumed-pdf (UF) model

In the classical UF modeling approach the filtered rate is calculated using,

w ˙ ¯ c ( x _ ,t )= ρ ¯ ( x _ ,t ) 0 1 w ˙ cL ( ζ ) ρL( ζ ) p ˜ ( ζ; c ˜ ( x _ ,t ), σ 2 ( x _ ,t ) )dζ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaWaaeWaaeaa jugibiqadIhagaqhaiaacYcacaWG0baajuaGcaGLOaGaayzkaaqcLb sacqGH9aqpcuaHbpGCgaqeaKqbaoaabmaabaqcLbsaceWG4bGba0ba caGGSaGaamiDaaqcfaOaayjkaiaawMcaamaapehabaWaaSaaaeaaju gibiqadEhagaGaaKqbaoaaBaaabaqcLbmacaWGJbGaamitaaqcfaya baWaaeWaaeaajugibiabeA7a6bqcfaOaayjkaiaawMcaaaqaaKqzGe GaeqyWdiNaamitaKqbaoaabmaabaqcLbsacqaH2oGEaKqbakaawIca caGLPaaaaaaabaqcLbmacaaIWaaajuaGbaqcLbmacaaIXaaajugibi abgUIiYdGabmiCayaaiaqcfa4aaeWaaeaajugibiabeA7a6jaacUda ceWGJbGbaGaajuaGdaqadaqaaKqzGeGabmiEayaaDaGaaiilaiaads haaKqbakaawIcacaGLPaaajugibiaacYcacqaHdpWCjuaGdaahaaqa beaajugWaiaaikdaaaqcfa4aaeWaaeaajugibiqadIhagaqhaiaacY cacaWG0baajuaGcaGLOaGaayzkaaaacaGLOaGaayzkaaqcLbsacaWG KbGaeqOTdOhaaa@7F34@ (10)

where ζ is the sample space variable for c, w ˙ cL MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiaajuaGdaWgaaqaaKqzadGaam4yaiaadYeaaKqbagqaaaaa@3BAE@ is the laminar progress variable rate, and ρL is the laminar density. The progress variable pdf is taken to be a β-function in accordance to usual practice in presumed pdf methods i.e.,

p ˜ ( ζ )= 1 C ζ a1 ( 1ζ ) b1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGWb GbaGaajuaGdaqadaqaaKqzGeGaeqOTdOhajuaGcaGLOaGaayzkaaqc LbsacqGH9aqpjuaGdaWcaaqaaKqzGeGaaGymaaqcfayaaKqzGeGaam 4qaaaacqaH2oGEjuaGdaahaaWcbeqaaKqzadGaamyyaiabgkHiTiaa igdaaaqcfa4aaeWaaOqaaKqzGeGaaGymaiabgkHiTiabeA7a6bGcca GLOaGaayzkaaqcfa4aaWbaaSqabeaajugWaiaadkgacqGHsislcaaI Xaaaaaaa@52A7@ (11)

where C is a normalisation constant and where the parameters a, b are chosen so that the filtered progress variable c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ and the variance σ2 are recovered. These are given by: a= c ˜ ( 1/g 1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacaWGHb Gaeyypa0Jabm4yayaaiaqcfa4aaeWaaeaadaWcgaqaaKqzGeGaaGym aaqcfayaaKqzGeGaam4zaaaacqGHsislcaaIXaaajuaGcaGLOaGaay zkaaaaaa@411E@ , and b=( 1 c ˜ )( 1/g 1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacaWGIb Gaeyypa0tcfa4aaeWaaeaajugibiaaigdacqGHsislceWGJbGbaGaa aKqbakaawIcacaGLPaaadaqadaqaamaalyaabaqcLbsacaaIXaaaju aGbaqcLbsacaWGNbaaaiabgkHiTiaaigdaaKqbakaawIcacaGLPaaa aaa@456D@ , where g= σ 2 c ˜ ( 1 c ˜ ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacaWGNb Gaeyypa0tcfa4aaSGaaeaajugibiabeo8aZLqbaoaaCaaabeqaaKqz adGaaGOmaaaaaKqbagaajugibiqadogagaacaaaajuaGdaqadaqaaK qzGeGaaGymaiabgkHiTiqadogagaacaaqcfaOaayjkaiaawMcaaaaa @45EA@ .

Note that in practice, the integration as specified by Eq. 10 may be problematic since the progress variable pdf takes non-finite values in the case the pdf is bimodal. A way around this issue is to expand the integral in Eq. 10 and use the cumulative distribution function P ˜ ( ζ; c ˜ , σ 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGqb GbaGaajuaGdaqadaqaaKqzGeGaeqOTdONaai4oaiqadogagaacaiaa cYcacqaHdpWCjuaGdaahaaqabeaajugWaiaaikdaaaaajuaGcaGLOa Gaayzkaaaaaa@431D@ instead for calculating the mean of a variable y as, y ¯ ( x _ ,t )= ρ ¯ ( x _ ,t ) ρ L ( 1 ) y L ( 1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOabmyEay aaraWaaeWaaeaajugibiqadIhagaqhaiaacYcacaWG0baajuaGcaGL OaGaayzkaaGaeyypa0ZaaSaaaeaacuaHbpGCgaqeamaabmaabaGabm iEayaaDaGaaiilaiaadshaaiaawIcacaGLPaaaaeaajugibiabeg8a YLqbaoaaBaaabaqcLbmacaWGmbaajuaGbeaadaqadaqaaKqzGeGaaG ymaaqcfaOaayjkaiaawMcaaaaacaWG5bWaaSbaaeaajugWaiaadYea aKqbagqaamaabmaabaGaaGymaaGaayjkaiaawMcaaaaa@5327@

ρ ¯ ( x _ ,t ) 0 1 1 ρ L ( ζ ) ( ρ L ( ζ ) d yL ( ζ ) dζ y L ( ζ ) d ρL ( ζ ) dζ ) P ˜ ( ζ; c ˜ , σ 2 )dζ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqGHsi slcuaHbpGCgaqeaKqbaoaabmaakeaajugibiqadIhagaqhaiaacYca caWG0baakiaawIcacaGLPaaajuaGdaWdXbGcbaqcfa4aaSaaaOqaaK qzGeGaaGymaaGcbaqcLbsacqaHbpGCjuaGdaWgaaWcbaqcLbmacaWG mbaaleqaaKqbaoaabmaakeaajugibiabeA7a6bGccaGLOaGaayzkaa aaaaWcbaqcLbmacaaIWaaaleaajugWaiaaigdaaKqzGeGaey4kIipa juaGdaqadaGcbaqcLbsacqaHbpGCjuaGdaWgaaWcbaqcLbmacaWGmb aaleqaaKqbaoaabmaakeaajugibiabeA7a6bGccaGLOaGaayzkaaqc fa4aaSaaaOqaaKqzGeGaamizaKqbaoaaBaaaleaajugWaiaadMhaca WGmbaaleqaaKqbaoaabmaakeaajugibiabeA7a6bGccaGLOaGaayzk aaaabaqcLbsacaWGKbGaeqOTdOhaaiabgkHiTiaadMhajuaGdaWgaa qaaKqzadGaamitaaqcfayabaWaaeWaaeaajugibiabeA7a6bqcfaOa ayjkaiaawMcaamaalaaakeaajugibiaadsgajuaGdaWgaaWcbaqcLb macqaHbpGCcaWGmbaaleqaaKqbaoaabmaakeaajugibiabeA7a6bGc caGLOaGaayzkaaaabaqcLbsacaWGKbGaeqOTdOhaaaGccaGLOaGaay zkaaqcLbsaceWGqbGbaGaajuaGdaqadaqaaKqzGeGaeqOTdONaai4o aiqadogagaacaiaacYcacqaHdpWCjuaGdaahaaqabeaajugWaiaaik daaaaajuaGcaGLOaGaayzkaaqcLbsacaWGKbGaeqOTdOhaaa@929E@ (12)

where the derivatives of y with respect to ζ are typically well-defined for quantities of interest, and so is the cumulative distribution function P ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGqb GbaGaaaaa@3769@ .

Filtered laminar flame model

In the FLF approach,7,13 the filtered rate is calculated from look-up tables constructed using filtered laminar profiles of a canonical flame solution e.g. a 1D unstrained flame. This process essentially corresponds to a filtered laminar flame pdf approach,7,13

w ˙ ¯ c ( x _ ,t )= ρ ¯ ( x _ ,t ) 0 1 w ˙ cL ( ζ ) ρL( ζ ) p ˜ L ( ζ; c ˜ ( x _ ,t ); σ 2 ( x _ ,t ) )dζ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzGeGaam4yaaqcfayabaWaaeWaaeaa jugibiqadIhagaqhaiaacYcacaWG0baajuaGcaGLOaGaayzkaaqcLb sacqGH9aqpcuaHbpGCgaqeaKqbaoaabmaabaqcLbsaceWG4bGba0ba caGGSaGaamiDaaqcfaOaayjkaiaawMcaamaapehabaWaaSaaaeaaju gibiqadEhagaGaaKqbaoaaBaaabaqcLbsacaWGJbGaamitaaqcfaya baWaaeWaaeaajugibiabeA7a6bqcfaOaayjkaiaawMcaaaqaaKqzGe GaeqyWdiNaamitaKqbaoaabmaabaqcLbsacqaH2oGEaKqbakaawIca caGLPaaaaaaabaqcLbsacaaIWaaajuaGbaqcLbsacaaIXaaacqGHRi I8aiqadchagaacaKqbaoaaBaaabaqcLbmacaWGmbaajuaGbeaadaqa daqaaKqzGeGaeqOTdONaai4oaiqadogagaacaKqbaoaabmaabaqcLb saceWG4bGba0bacaGGSaGaamiDaaqcfaOaayjkaiaawMcaaiaacUda jugibiabeo8aZLqbaoaaCaaabeqaaKqzGeGaaGOmaaaajuaGdaqada qaaKqzGeGabmiEayaaDaGaaiilaiaadshaaKqbakaawIcacaGLPaaa aiaawIcacaGLPaaajugibiaadsgacqaH2oGEaaa@7E47@ (13)

Where the filtered laminar flame pdf is given by,

p ˜ L ( ζ; c ˜ ; σ 2 )= ρ L ( ζ ) ρ L ΔL | dζ dx | 1 G( x L ΔL x( ζ ) ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGWb GbaGaajuaGdaWgaaqaaKqzadGaamitaaqcfayabaWaaeWaaeaajugi biabeA7a6jaacUdaceWGJbGbaGaacaGG7aGaeq4Wdmxcfa4aaWbaae qabaqcLbmacaaIYaaaaaqcfaOaayjkaiaawMcaaiabg2da9maalaaa baGaeqyWdi3aaSbaaeaacaWGmbaabeaadaqadaqaaiabeA7a6bGaay jkaiaawMcaaaqaaiabeg8aYnaaBaaabaqcLbmacaWGmbaajuaGbeaa daahaaqabeaajugWaiabgkHiTiabfs5aejaadYeaaaaaaKqbaoaaem aabaWaaSaaaeaacaWGKbGaeqOTdOhabaGaamizaiaadIhaaaaacaGL hWUaayjcSdWaaWbaaeqabaqcLbmacqGHsislcaaIXaaaaKqbakaadE eadaqadaqaaiaadIhadaWgaaqaaKqzadGaamitaaqcfayabaWaaWba aeqabaqcLbmacqqHuoarcaWGmbaaaKqbakabgkHiTiaadIhadaqada qaaiabeA7a6bGaayjkaiaawMcaaaGaayjkaiaawMcaaaaa@7101@ (14)

And the filter function G in this case is in 1D. In the FLF approach, the laminar flame filter size ΔL which is used to parameterize the laminar flame pdf is generally smaller than the actual LES filter size Δ. The laminar flame filter ΔL, and the laminar flame position xL, are chosen so that the progress variable and its variance as obtained from the LES, match the corresponding 1D-filtered laminar flame values i.e. c ˜ ( x _ ,t;Δ )= c ˜ L ( x L Δ L ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaajuaGdaqadaGcbaqcLbsaceWG4bGba0bacaGGSaGaamiDaiaa cUdacqqHuoaraOGaayjkaiaawMcaaiabg2da9KqzGeGabm4yayaaia GcdaWgaaWcbaqcLbmacaWGmbaaleqaaOWaaeWaaeaajugibiaadIha kmaaBaaaleaajugWaiaadYeaaSqabaGccqGHsisljugibiabfs5aeP WaaSbaaSqaaKqzadGaamitaaWcbeaaaOGaayjkaiaawMcaaaaa@4E86@ , and σ 2 ( x _ ,t;Δ )= σ 2 L ( x L Δ L ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHdp WCjuaGdaahaaqabeaajugWaiaaikdaaaqcfa4aaeWaaOqaaKqzGeGa bmiEayaaDaGaaiilaiaadshacaGG7aGaeuiLdqeakiaawIcacaGLPa aacqGH9aqpjugibiabeo8aZLqbaoaaCaaabeqaaKqzadGaaGOmaaaa kmaaBaaaleaajugWaiaadYeaaSqabaGcdaqadaqaaKqzGeGaamiEaO WaaSbaaSqaaKqzadGaamitaaWcbeaakiabgkHiTKqzGeGaeuiLdqKc daWgaaWcbaqcLbmacaWGmbaaleqaaaGccaGLOaGaayzkaaaaaa@5552@ . This ensures that the filtered-flame pdf recovers the LES mean and variance values.

In practice, an unstrained laminar flame is filtered for a wide range of filter widths in order to obtain a table of values of c ˜ L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaakmaaBaaaleaajugWaiaadYeaaSqabaaaaa@39BC@ , and σ 2 L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHdp WCjuaGdaahaaqabeaajugWaiaaikdaaaGcdaWgaaWcbaqcLbmacaWG mbaaleqaaaaa@3D22@ . The LES values c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ and σ 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqaHdp WCjuaGdaahaaqabeaajugWaiaaikdaaaaaaa@3AE2@ are then used in order to determine the corresponding filter width Δ L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqqHuo arkmaaBaaaleaajugWaiaadYeaaSqabaaaaa@3A2B@ and spatial position xL in the table which match the filtered laminar flame values. The mean rate is then calculated by filtering the laminar flame profile at Δ L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqqHuo arkmaaBaaaleaajugWaiaadYeaaSqabaaaaa@3A2B@ and obtaining its value at xL, as discussed in.7,13

Results and discussion

Modeled pdfs comparison

Irrespective of the flamelet model used, small discrepancies in the modeled variance using deconvolution (or any other variance model) may cause large discrepancies in the modeled pdf and as a result large discrepancies in the filtered progress variable reaction rate. The aim of this section is to investigate this point for the UF and FLF models. In order to do that, the progress variable pdfs, as defined for the UF and FLF models in sections 3.2 and 3.3 respectively, are calculated at every point in the domain on the LES mesh. The pdfs are parameterised in two different ways:

  1. Using the actual variance as obtained by filtering the DNS data (and sampling onto the LES mesh) to obtain p(ζ) and pL(ζ).
  2. Using the modeled variance as obtained using deconvolution to obtain p*( ζ) and p*L(ζ).

Thus, for every grid-point on the LES mesh, any changes in the modeled pdfs are a result of differences in the variance only, since c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ is always the same. The pdfs as obtained at every grid-point and for every time-step, are then averaged in the homogeneous directions (y; z) and for all time-steps, in order to obtain the averaged x-direction pdfs < p(ζ; x) >, < p*(ζ; x) > and so on. The Favre-filtered progress variable c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ is also averaged in the same way to obtain < c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ > as a function of x. This averaging procedure allows us to make a direct comparison between the two modeled pdfs, in order to investigate whether the deconvolution-based modeling for the variance induces any significant bias.

Figure 2(a)-2(c) shows the averaged pdfs < p(ζ) > and < p*(ζ) > for the UF model for cases A-C respectively. The pdfs at a spatial location corresponding to < c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ >=0.5 are shown, where there is a significant contribution in the burning mode part of the pdf. The agreement between <p> and <p*> is best for the lowest filter width,

Δ+=1, as expected for all three cases. This result is consistent with the results in32 where the variance estimate is most accurate for the lowest filter widths. For increasing filter widths, the agreement somewhat deteriorates. This is also to be expected since deconvolution accuracy is inversely proportional to filter width size [32]. This results in the variance being somewhat over-estimated, hence the generally broader pdfs observed in Figure 2(a)-2(c) for Δ+=2, and, 3. Overall though, the agreement between the two pdfs remains good. The same process is repeated for the FLF model. The results for the filtered laminar flame pdfs are shown in Figure 3(a)-3(c) for cases A-C respectively. These pdfs also correspond to a spatial position having < c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ >=0.5. A similarly good agreement between the modeled pdfs is observed in the case of the FLF model. These results indicate that the deconvolution-based approach for modeling the variance does not introduce a significant bias in the modeled pdfs for either amelet model used.

Figure 1 Comparison of instantaneous progress variable rates for Δ+=1.0: (a) w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ , (b) w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : UF, (c) w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : FLF, (d) w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : UF MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0caaa@36E3@ IDEF, (e) w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : FLF MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0caaa@36E3@ IDEF.

  • Figure 2 Progress variables pdfs for the unstrained amelet model at < c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ >=0.5 for (a) case A, (b) case B, and (c) case C. p is obtained using the actual variance, and p* is obtained using deconvolution for modeling the variance.

  • Figure 3 Progress variables pdfs for the FLF model at < c ˜ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWGJb GbaGaaaaa@377C@ >=0.5 for (a) case A, (b) case B, and (c) case C. pL is obtained using the actual variance, and p*L is obtained using deconvolution for modeling the variance.

Modeled rates comparison

The results of the previous section are encouraging, and the aim of this section is to examine whether the hybrid amelet-deconvolution approach introduces a significant bias in the modeled rates. The pdfs as obtained on the LES mesh are used in order to calculate the progress variable rates for the UF and FLF models (on the LES mesh) as described in sections 3.2 and 3.3. The pdfs are calculated using both the actual variance as obtained

by filtering the DNS data, and the modeled variance using deconvolution, as explained in the previous section. In the case of the UF model, this leads to filtered rate predictions which are denoted here by w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : UF, and w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : UF MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0caaa@36E3@ IDEF. In the case of the FLF model, the corresponding rate predictions are denoted as w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : FLF, and w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : FLF MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0caaa@36E3@ IDEF. Figure 1 shows a visual comparison between the actual rate (a) and the four modeled rates (b)-(e) as defined above. This is done for the highest turbulence level, case C, for Δ+=1.0. Figure 1 shows instantaneous contours of the progress variable rate in the direction of mean flame propagation. These contours are equally normalized using the maximum instantaneous value of the DNS-filtered rate in order to elucidate any differences in the modeled rates. Clearly, the flame is highly convoluted due to the intense turbulence level as one may observe from 1(a). Nevertheless both the UF (1(b)) and FLF models (1(c)) provide reasonably good estimates of the progress variable rate. In addition, comparing the results in Figure 1(b) and 1(d) for the UF model it is clear that the deconvolution method provides equally well predictions of the progress variable rate. The same applies for the FLF models by comparing Figures (c) and (e).

Another important point to note is the effect of the filter on the rate: for increasing filter widths the filtered rates generally decrease in comparison with the un-filtered laminar flame result. Simply evaluating the rate using the filtered values would not result in a correct evaluation of the filtered rate since for a non-zero filter width  (ρ; T; Yk) . Hence the need for modeling the filtered rate. In order to warrant a more quantitative evaluation, conditional averages based on the Favre-filtered progress variable, , are also calculated for each case, and for each flamelet model. The conditional rates for the UF and FLF models for each filter width are shown in Figure 4(a)-4(c) for case A, Figure 5(a)-5(c) for case B, and Figure 6(a)-6(c) for case C. Also shown in each figure is the laminar flame result, and the actual conditional rate as obtained by filtering the data on the DNS mesh (and sampling onto the LES mesh). The assessment is also complemented by using a suitable local error measure between the modeled rates. This error is defined as for two variables x and y. In the case of the UF model, the errors e(: UF; : UFIDEF) are calculated in order to elucidate any deviations from the UF model predictions. In the case of the FLF model, a similar error is calculated e (: FLF; : FLFIDEF). For each error measure, the threshold is taken to be at 10% of the maximum instantaneous value of the progress variable rate, so that regions of significant heat release are considered while regions with low heat release do not affect the statistics. These errors are shown in Table 4 and Table 5 for the UF and FLF models respectively.

  • Figure 4 Conditionally averaged filtered progress variable rate for case A.

Case

+=1         ∆+=2        ∆+=3

A

1.2

5.7

8.8

B

1.1

3.5

5.3

C

1.7

6.6

10.4

Table 4 Percentage error between modeled rates for the UF model, using IDEF for the pdf variance, and using the actual variance for the pdf, e( w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : UF MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0caaa@36E3@ IDEF; w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : UF).

The results in Figure 4, Figure 5 and Figure 6 show a reasonably good agreement between the modeled rates both for the UF and FLF approaches. This is also supported by the results in Table 4 and Table 5. For both models the bias in using deconvolution for modeling the variance is lowest for the smallest filter width which is to be expected. The bias generally increases for increasing filter width for each case, but remains within reasonable bounds for the largest filter width. It is interesting to note that the largest error for the UF model occurs for case C at Δ+=3.0, while in the case of the FLF model this occurs for case A instead, at Δ+=3.0. There does not appear to be a correlation in the error with increasing turbulence level for either the UF or FLF models-further parametric DNS studies are required for determining such behaviour in order to cover a wider space of possible progress variable variance values. What these results clearly indicate however is that IDEF produces quantitatively good estimates of the progress variable rate, with no significant bias in the models' predictions.

  • Figure 5 Conditionally averaged filtered progress variable rate for case B.

  • Figure 6 Conditionally averaged filtered progress variable rate for case C.

Case

+=1         ∆+=2         ∆+=3

A

6.1

9.3

11.5

B

4.9

7.6

9.2

C

4.5

7.8

7.4

Table 5 Percentage error between modeled rates for the FLF model, using IDEF for the pdf variance, and using the actual variance for the pdf, e( w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : FLF MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0caaa@36E3@ IDEF; w ˙ ¯ c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaceWG3b GbaiGbaebajuaGdaWgaaqaaKqzadGaam4yaaqcfayabaaaaa@3AF4@ : FLF).

Conclusion

A direct numerical simulation database of a freely-propagating turbulent premixed flame in a canonical inflow-outflow configuration is used in order to investigate whether a deconvolution-based amelet approach for modeling the filtered progress variable reaction rate in the context of LES, introduces any significant bias in the predictions of the flamelet models used. The assessment of the method is conducted a priori, on a simulated LES mesh which is constructed using the DNS-filtered data. On the LES mesh, the filtered fields ρ ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacuaHbp GCgaqeaaaa@385D@ and ρc ¯ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aa0aaae aacqaHbpGCcaWGJbaaaaaa@393D@ are deconvoluted, and the progress variable variance is obtained by explicitly filtering the deconvoluted fields on the same mesh. The variance as obtained is used to parameterize the progress variable pdfs in two popular flamelet methods namely the unstrained flamelet model using a β-pdf, and the filtered laminar flame model.
The results of this study show that the deconvolution-based approach for modeling the variance does not introduce any significant bias in the shape of the modeled progress variable pdfs or on the values of the modeled reaction rates for both flamelet models tested. As a result, deconvolution can be used in conjuction with such flamelet methods for modeling the progress variable reaction rate, thus rendering a more parameter-free flamelet-based modeling strategy.

Acknowledgements

This work made use of the facilities of HECToR, the UK's national high-performance computing service, which is provided by UoE HPCx Ltd at the University of Edinburgh, Cray Inc and NAG Ltd, and funded by the Office of Science and Technology through EPSRC's High End Computing Programme.

Conflicts of Interest

The author declares that they have no conflict of interest.

References

  1. Sagaut P. Large Eddy Simulation for Incompressible Flows: An Introduction. 2nd ed. Newyork: Springer-Verlag; 2001.
  2. Bray KNC, Moss JB. A unified statistical model of the premixed turbulent flame. Acta Astron. 1977;4(3-4):291‒319.
  3. Gicquel LYM, Staffelbach G, Poinsot T. Large eddy simulations of gaseous flames in gas turbine combustion chambers. Prog En Combust Sci. 2012;38(6):782‒817.
  4. Veynante D, Vervisch L. Turbulent combustion modeling. Prog Energy Combust Sci. 2002;28(3):193‒266.
  5. Pitsch H. Large eddy simulation of turbulent combustion. Ann Rev Fluid Mech. 2006;38:453‒482.
  6. Domingo P, Vervisch L, Veynante D. Large-eddy simulation of a lifted methane jet flame in a vitiated co-flow. Combust Flame. 2008;152(3):415‒432.
  7. Moureau V, Domingo P, Vervisch L. From large-eddy simulation to direct numerical simulation of a lean premixed swirl flame: filtered laminar flame-PDF modeling. Combust Flame. 2011;158(7):1340‒1357.
  8. Galpin J, Naudin A, Vervisch L, et al. Large-eddy simulation of a fuel-lean premixed turbulent swirl-burner. Combust Flame. 2008;155(1-2):247‒266.
  9. Kempf A, Malalasekera W, Ranga-Dinesh KJ, et al. Large eddy simulations of swirling non-premixed flames with amelet models: a comparison of numerical methods. Flow Turbulence Combust. 2008;81(4):523‒561.
  10. Malalasekera W, Ranga-Dinesh KJ, Ibrahim SS, et al. LES of recirculation and vortex breakdown in swirling flames. Combust Sci Techn. 2008;180(5):809‒832.
  11. Domingo P, Vervisch L, Veynante D. Large-eddy simulation of a lifted methane jet flame in a vitiated co-flow. Combust Flame. 2008;152(3):415‒432.
  12. Knudsen E, Kolla H, Hawkes ER, et al. LES of a premixed jet flame DNS using a strained flamelet model. Combust Flame. 2013;160(12):2911‒2927.
  13. Nambully S, Domingo P, Moureau V, et al. A Filtered-Laminar-Flame PDF sub-grid scale closure for LES of premixed turbulent flames. Part I: Formalism and application to a bluff-body burner with differential diffusion. Combust Flame. 2014;161(7):1756‒1774.
  14. Domingo P, Vervisch L. Large eddy simulation of premixed turbulent combustion using approximate deconvolution and explicit flame filtering. Proc Combust Inst. 2015;35(2):1349‒1357.
  15. Langella I, Swaminathan N. Unstrained and strained flamelets for LES of premixed combustion. Combust Th Model. 2016;20(3):410‒440.
  16. Langella I, Swaminathan N, Pitz RW. Application of unstrained flamelet SGS closure for multi-regime premixed combustion. Combust Flame. 2016;173:161‒178.
  17. Boger M, Veynante D, Boughanem H, et al. Direct Numerical Simulation analysis of flame surface density concept for Large Eddy Simulation of turbulent premixed combustion. Proc Combust Inst. 1998;27(1):917‒925.
  18. Cook AW, Riley JJ. A sub-grid model for equilibrium chemistry in turbulent flows. Phys Fluids. 1994;6(8):2868‒2870.
  19. Cook AW. Determination of the constant coefficient in scale similarity models of turbulence. Phys Fluids. 1997;9(5):1485‒1487.
  20. Girimaji S, Zhou Y. Analysis and modeling of subgrid scalar mixing using numerical data. Phys Fluids. 1996;8(5):1224‒1236.
  21. Pierce CD, Moin P. A dynamic model for subgrid-scale variance and dissipation rate of a conserved scalar. Phys Fluids. 1998;10(12):3041‒3044.
  22. Veynante D, Knikker R. Comparison between LES results and experimental data in reacting flows. J Turbulence. 2006;7:1‒20.
  23. Nikolaou ZM, Swaminathan N. Assessment of FSD and SDR closures for turbulent flames of alternative fuels. Flow Turb Combust. 2018; 1‒16 p.
  24. Balarac G, Pitsch H, Raman V. Development of a dynamic model for the sub filter scalar variance using the concept of optimal estimators. Phys Fluids. 2008;20(3):1‒9.
  25. Bose S, Moin P. A dynamic slip boundary condition for wall-modeled large-eddy simulation. Phys Fluids. 2014;26(1):1‒18.
  26. Locci C, Vervisch L. Eulerian scalar projection in Lagrangian point source context: An approximate inverse filtering approach. Flow Turbulence Combust. 2016;97(1):363‒368.
  27. Mathew J. Large Eddy Simulation of a premixed flame with approximate deconvolution modeling. Proc Combust Inst. 2002;29(2):1995‒2000.
  28. Vreman AW, Bastiaans RJ, Geurts BJ. A similarity sub-grid model for premixed turbulent combustion. Flow Turbulence Combust. 2009;82(2):233‒248.
  29. Domingo P, Vervisch L. DNS and approximate deconvolution as a tool to analyse one-dimensional filtered flame sub-grid scale modeling. Combust Flame. 2017;177:109‒122.
  30. Mehl C, Idier J, Fiorina B. Evaluation of deconvolution modeling applied to numerical combustion. Combust Th Model. 2018;22(1):38‒70.
  31. Wang Q, Ihme M. Regularized deconvolution method for turbulent combustion modeling. Combust Flam. 2017;176:125‒142.
  32. Nikolaou ZM, Vervisch L. A priori assessment of an iterative deconvolution method for LES sub-grid scale variance modeling. Flow Turbulence Combust. 2018;101(1):33‒53.
  33. Nikolaou ZM, Cant RS, Vervisch L. Scalar flux modeling in turbulent flames using iterative deconvolution. Phys Rev Fluids. 2018;3(4):1‒19.
  34. Cant RS. SENGA2 User Guide. CUED/A-THERMO/TR67; 2012. 46 p.
  35. Nikolaou Z, Swaminathan N. A 5-step reduced mechanism for combustion of CO/H2/H2O/CH4/CO2 mixtures with low hydrogen/methane and high H2O content. Combust Flame. 2013;160(1):56‒75.
  36. Peters N. Laminar flamelet concepts in turbulent combustion. Symposium on Combustion. 1988;21(1):1231‒1250.
  37. Nikolaou ZM, Swaminathan N, Chen JY. Evaluation of a reduced mechanism for turbulent premixed combustion. Combust Flame. 2014;161(12):3085‒3099.
  38. Nikolaou ZM, Swaminathan N. Direct numerical simulation of complex fuel combustion with detailed chemistry: physical insight and mean reaction rate modeling. Comb Sci Tech. 2015;187(11):1759‒1789.
  39. Jansson PA. Deconvolution with applications in spectroscopy. New York: Academic Press; 1984. 342 p.
  40. Van Cittert PH. On the influence of the gap width on the intensity distribution in spectral lines. Magazine for Physics. 1931;69(5-6):298‒308.
  41. Benjamin P. A quantitative evaluation of various iterative deconvolution algorithms. IEEE Transactions on Instrumentation and Measurement. 1991;40(3):558‒562.
Creative Commons Attribution License

©2018 Nikolaou, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.