The Eddy Dissipation Concept—Analysis of Different Fine Structure Treatments for Classical Combustion

The Eddy Dissipation Concept (EDC) is common in modeling turbulent combustion. Several model improvements have been proposed in literature; recent modifications aim to extend its validity to Moderate or Intense Low oxygen Dilution (MILD) conditions. In general, the EDC divides a fluid into a reacting and a non-reacting part. The reacting part is modeled as perfectly stirred reactor (PSR) or plug flow reactor (PFR). EDC theory suggests PSR treatment, while PFR treatment provides numerical advantages. Literature lacks a thorough evaluation of the consequences of employing the PFR fine structure treatment. Therefore, these consequences were evaluated by employing tests to isolate the effects of the EDC variations and fine structure treatment and by conducting a Sandia Flame D modeling study. Species concentration as well as EDC species consumption/production rates were evaluated. The isolated tests revealed an influence of the EDC improvements on the EDC rates, which is prominent at low shares of the reacting fluid. In contrast, PSR and PFR differences increase at large fine fraction shares. The modeling study revealed significant differences in the EDC rates of intermediate species. Summarizing, the PFR fine structure treatment might be chosen for schematic investigations, but for detailed investigations a careful evaluation is necessary.


Introduction
Nowadays, the Eddy Dissipation Concept (EDC) [1][2][3] is very common in turbulent combustion modeling. Since the initial publication of the EDC in the 1970s, several modifications of the EDC were proposed in literature [4][5][6][7][8], all of them claiming to improve the modeling results. In recent years, emphasis is on the extension of the EDC for Moderate or Intense Low oxygen Dilution (MILD) combustion conditions [7][8][9]. In contrast to classical turbulent combustion, MILD combustion features widely distributed reaction zones and slower reactions due to lower oxygen concentrations as well as slightly lower turbulence [10][11][12][13]. In classical turbulent combustion, chemical reactions occur only in small areas, the so-called fine structures, where educts are mixed at molecular scale. These fine structures are of Kolmogorov length scale in at least two spatial dimensions and of larger scale in the third [3,4,14,15]. A fine structure share (γ * ) and a fine structure residence time (τ * ) are defined for characterization based on turbulence parameters. The non-reacting fluid is called the surroundings in the EDC context. At sufficiently high Reynolds numbers, combustion reactions are limited by the

Theory
The EDC splits the reacting fluid into a reacting zone (fine structures) and a non-reacting zone (surroundings). The fine structures are of vortex-sheet or tubelike shape and their size is of Kolmogorov scale in two dimensions but at bigger scale in the third dimension [3,4,14,15]. Reactions occur in those spatially isolated regions, where educt gas phase species mix at the molecular scale. This is done by turbulent mixing in the Kolmogorov structures where fluid motion is dissipated to thermal energy due to viscous forces. Ertesvåg and Magnussen [25] compared the EDC theoretical turbulence energy cascade presumed for the EDC with turbulence energy spectra of reacting flows and showed that cascade reproduces the energy spectrum well. The different EDC modifications are discussed in more detail in the next paragraphs. There, the fine structure, surrounding and fluid mean properties are denoted with * , • and − , respectively.

'EDC1981'
In the original EDC [3], the fine structure reaction rate R * i is equal to a mass transfer rate between the fine structures and the surroundings and is calculated as follows for species i: τ * is the fine structure residence time or EDC mixing time, ρ * and ρ • are the fine structure and surroundings densities, and c * i and c • i are the fine structure and surrounding concentrations of species i. The inverse fine structure residence time is proportional to the dimensionless fine structure mass transfer rate (ṁ * ). The fine structure and surrounding quantities are linked by the fluid average properties and the fine structure mass share γ * : A dimensionless fine structure length fraction is defined as the cubic root of the fine structure share: The fine structure length fraction γ L and the fine structure residence time τ * are calculated as ratios of the turbulent kinetic energy (k), the turbulent dissipation rate ( ), and the model constants C τ or C γ . An alternative definition based on the turbulent Reynolds number (Re t , see Equation (4)) is also possible. The model constants are proportional to the C D1 and C D2 constants of the turbulent energy cascade [25]: By substituting Equations (3) and (2) into Equation (1) and replacing the fine structure mass transfer rate with the mean mass transfer rate (ṁ =ṁ * γ 3 L = γ 3 L /τ * ), the mean reaction rate is obtained as: ρ, Y i , and Y * denote the mean density, and the mean and the fine structure mass fractions of species i. χ represents the reacting fraction of the fine structures. and is assumed to be unity, although this might be wrong, especially for low turbulent flows as indicated by Lewandowski et al. [26]. They also speculate a relation of χ and Re t at least for MILD combustion conditions [26,27].
If the reacting parts of the fine structures are assumed to be homogeneously mixed, they can be modeled as perfectly stirred reactors (PSR) [4]. The PSR has to be solved to steady-state in order to obtain the fine structure composition and, thus, the EDC reaction rates. Since computing the steady-state is numerically expensive, the fine structures are alternatively treated as plug flow reactors (PFR) and solved for the fine structure residence time [17].

'EDC1996'
In 1996, Gran and Magnussen [5] proposed modifying the mean reaction rate to include mixing effects of fresh gas and fine structures. The conversion ratio of the mass exchange rate is changed from γ 3 L /τ * to γ 2 L /τ * . Accordingly, the mean reaction rate changes to:

'EDC2005'
In 2005, Magnussen [6] proposed a modification of the fine structure definition (γ * ) to incorporate the fact that reactions can extend outside the fine structures (Equation 10). Therefore, the mean reaction rate is modified to: Moreover, the relation between the mean fluid properties and the surrounding and fine structures changes:

'EDC2016'
Recently, emphasis is put on the applicability of the EDC to non-classical combustion conditions, e.g., MILD combustion.
First approaches tried to modify the EDC model constants to improve results [9,17,20,21]. However, these modifications were only applicable for specific cases. Parente et al. [7] proposed a more sophisticated approach to adapt the EDC. Since Kolmogorov eddies enter and thicken the reaction zone, the characteristic fine structure velocity is proportional to the flame speed in MILD combustion [12]. Therefore, the turbulent energy cascade is related to the Damköhler number at Kolmogorov scale (Da η ) and the turbulent Reynolds number (Re t ). The modified energy cascade gives alternative definitions of the model constants C γ and C τ , which can be seen as local qualitative estimates: The mean reaction rate is calculated according to Equation (8) by changing the model constants when calculating γ L Equation (5) and τ * Equation (6). Both EDC model constants are limited to their default values, e.g., C γ below 2.1377 and C τ above 0.4082.

'EDC2017'
Bao [8] modified Parente's idea and eliminated the flame speed from the derivations of the model constants by substituting it with the characteristic length and the chemical time scale Since the unknown proportional factor in the flame speed definition is eliminated, a quantitative relation for the local model constants C τ and C γ was derived. The definitions of Bao and Parente et al. show similar structures: The procedure for calculating the mean EDC rate is taken from Parente et al.
The review of the different EDC modifications shows three approaches for the model adaption: conversion from fine structure to mean reaction rates [5,6], modified model constants, e.g., static [17,20,21] or dynamic based on local properties [7,8], and fine structure treatment [5,16,17]. Figure 1 summarizes the differences of the EDC modifications available in literature based on a schematic illustration of a representative fine structure. The arrows surrounding the fine structures indicate the relative intensity of the mass transfer between non-reacting surroundings and the fine structures, which differs between the EDC versions. In case of the EDC-MILD, the colored arrows shall indicate that the fine fraction share and mass transfer rates can vary. indicates the fine structures, while represents the mass transfer between them and the surroundings.
illustrates the variable fine structure fractions and mass transfer intensities occurring in the EDC modifications for Moderate or Intense Low oxygen Dilution (MILD) combustion. Symbol sizes indicate their relative relations.

The Mean Reaction Rate Expression
Converting fine structure reaction rates to fluid mean rates by a characteristic volumetric ratio would be a logical choice Equation (16). However, different conversion ratios based on adapted model assumptions have been proposed in literature. Gran and Magnussen [5] modified the conversion of the mass transfer rate toṁ =ṁ * γ 2 L Equation (17). Actually, this means reaction rates are converted by a characteristic surface to volume ratio instead of a characteristic volume ratio. Keeping in mind that γ L should be between 0 and 1 to satisfy γ * ≤ 1, it is obvious that this modification increases the average reaction rate compared to the characteristic volume ratio conversion. In 2005, Magnussen [6] proposed a new modification Equation (18) based on a characteristic surface ratio. Therefore, both the mass transfer rate and the surrounding species mass fraction are converted by γ 2 L and 1 − γ 2 L , respectively. This surface ratio conversion gives the highest average reaction rates. In general, the conversion expressions differ only by the exponents of the characteristic fine structure length fraction γ L in the numerator and denominator. The following equations summarize the conversion factors of the EDC1981, EDC1996 and EDC2005 modification: EDC2005 : Obviously, all three expressions have a singularity near γ L = 1. Moreover, they can be approximated by γ 2 L or γ 3 L if γ L << 1, respectively. However, since the fine structure residence time (τ * ) should be smaller than or equal to a characteristic mixing time (τ mix ), the fine structure fraction should be limited in order to ensure a ratio (R τ ) of the fine structure to the mixing time below unity, e.g., for the original EDC [16][17][18]28]: As a result, the maximum γ L is different for the three conversion ratios and has values between 0.70 and 0.79. In contrast, Lewandowski et al. [26] proposed the trivial condition based on the mean and fine structure reaction rates, which give τ mix = τ * /γ 3 L : Requiring γ L < 1 allows significant higher species reaction rates. According to Lewandowski et al., a limitation of the fine structure length fraction is only necessary if the employed mixing time scale deviates from the EDC mixing time scale. Figure 2 shows the influence of the different conversion ratios on the average reaction rate. The top diagram depicts the conversion factor for all three conversion ratios, the dashed line represents the maximum conversion ratio to obey Equation (19). The bottom diagram depicts the conversion factors normalized by the volume ratio conversion factor Equation (7) over the whole γ L range. Both the surface to volume and the surface ratio give similar values for γ L < 0.4. Starting at γ L = 0.4, the surface to volume ratio approaches the volume conversion ratio.

Model Constant Modification
The C γ and C τ EDC constants were adapted to reproduce MILD combustion conditions by Parente et al. [7] and Bao [8]. The Damköhler number at Kolmogorov scale (Da η ) is of order unity in this combustion regime. Thus, the influence of the different model formulations is investigated for Da η < 2 and turbulent Reynolds numbers (Re t ) up to 1000. Figure 3 shows γ L (top) and τ * (bottom) for the modified EDC constants normalized by the results for the standard EDC constants on a double-log plot. The fine structure length fraction (γ L ) increases with increasing Damköhler number The normalized conversion ratios reveal that the surface conversion ratio results in approximately 1.5 times higher reaction source terms than the volume conversion ratio at high γ L . In contrast, at low γ L values (especially if smaller than 0.1), the surface to volume and the surface conversion ratios give several times higher reaction rates than the volume conversion ratio. The surface to volume conversion ratio is a hybrid version of the other two ratios: it is similar to the surface conversion ratio for low and similar to the volume conversion ratio for high γ L values. In general, the relative differences arising due to the different conversion ratios are higher at low γ L values (exceeding a factor of 100) and reduce to a factor below 2 for γ L > 0.5. However, the absolute differences are small for low γ L , since the conversion ratios are small compared to unity (γ L << 1).

Model Constant Modification
The C γ and C τ EDC constants were adapted to reproduce MILD combustion conditions by Parente et al. [7] and Bao [8]. The Damköhler number at Kolmogorov scale (Da η ) is of order unity in this combustion regime. Thus, the influence of the different model formulations is investigated for Da η < 2 and turbulent Reynolds numbers (Re t ) up to 1000. Figure 3 shows γ L (top) and τ * (bottom) for the modified EDC constants normalized by the results for the standard EDC constants on a double-log plot. The fine structure length fraction (γ L ) increases with increasing Damköhler number at a given Kolmogorov scale (Da η ) and increasing turbulent Reynolds number (Re t ). EDC2016 and EDC2017 modifications reach the original fine structure length fraction below a Damköhler number of around unity. If C γ is clipped by the EDC1981 value, the modifications affect the results for slow reactions in less turbulent regions. The slope of EDC2016 is higher, but the intercept is significantly lower than for the EDC2017.
The fine structure residence time (τ * ) behaves opposite to the fine structure length fraction. The EDC2016 and EDC2017 residence times approach the EDC1981 values between Damköhler numbers of 0.1 and 1 for the investigated turbulent Reynolds numbers (Re t ). The slopes of the EDC2016 and EDC2017 modifications are identical in the double-log representation, but the intercepts of the EDC2016 are significantly lower, indicating smaller fine structure residence times.
Farokh and Birouk [20,21] showed that the influences of C τ and C γ on temperature and species concentrations depend on the turbulence level. They performed a modeling study on Sandia Flames A and B; Flame A is a low turbulent flame (Re jet ≈ 1100), while Flame B is more turbulent (Re jet ≈ 17500). In case of low turbulence, the time scale ratio according to Equation (19) would exceed unity and, thus, fine structure shares are limited to ensure a proper ratio between mixing and fine structure residence time. If C τ is decreased in these regions, reaction rates increase and vice versa. In contrast, if the critical Reynolds number is exceeded, an increase of C τ gives higher reaction rates. For highly turbulent flames, e.g., Flame B, higher C γ values result in higher temperatures and reaction rates. The influence of C τ is less distinct than the one of C γ . However, increased C τ constants gave lower temperatures before and after the temperature peak, while the peak temperature was virtually the same. At peak temperatures, reactions are mixing controlled, since chemistry time scales are fast there. Up-and downstream, temperatures decrease and reactions become either mixing or chemical controlled.

Fine Structure Treatment-Perfectly Stirred Reactor (PSR) vs. Plug Flow Reactor (PFR)
According to Magnussen [4], fine structures can be seen as perfectly stirred reactors (PSR), since they provide a perfect analogy: the fine structure residence time τ * corresponds to the PSR residence time. The PSR inlet concentrations are equal to the species concentrations of the non-reacting surroundings and the PSR outlet concentrations are the fine structure concentrations. Neglecting radiative or diffusive thermal losses of the fine structures, the PSR is described by the following set of equations for species mass fraction (Y i ), enthalpy (h), and pressure (p): ω * i is the chemical species reaction rate, while the second term on the right hand side (RHS) represents the mixing between fine structures and surroundings. The PSR is solved to steady-state to obtain the fine structure composition. Solving this set of ordinary differential equations (ODE) is numerically expensive, since the chemical source terms are highly nonlinear. Furthermore, under certain conditions, periodic oscillating states might occur in PSRs, although steady input parameters are employed, due to chemical kinetic effects. An example of such a periodic oscillating state is given in Figure 4. A stoichiometric hydrogen mixture at 800 K and 0.1 bar is combusted in a PSR with a residence time of 0.75 s employing the kinetic mechanism of Connaire et al. [29]. The open source tool Cantera [30] was employed for the PSR calculation. Moreover, Cantera tutorials inspired the authors for the periodic oscillation example. A convenient way to get rid of the instabilities and to reduce the numerical effort, is to neglect the mixing term between the fine structures and the surroundings [9,[16][17][18]. This approach is referred to as the plug flow reactor (PFR) approach, since back-mixing is ignored in the species conservation equations and only chemistry is integrated over time. The integration time reduces to the fine structure residence time (τ * ). Using the same assumptions as for the PSR, the PFR is described by the following set of equations: Obviously, the results of these two reactor types should differ when using the same starting values. Since the fine structure species concentrations directly influence the mean reaction rates, the results using the PSR and PFR approach in flame modeling are expected to differ. Figure 5 gives a conceptional illustration of the PSR and the PFR fine structure treatment. In case of the PSR treatment, small eddies forming at the surface or diffusion effects enable mass exchange with the surroundings. This is indicated by tiny eddies surrounding the fine structures. In contrast, the PFR approach assumes isolated, possibly smooth structures lacking mass transfer with the surroundings. Treating the fine structures as PFR should increase the EDC mean rates compared to the PSR approach for fast reactions, since no educts are back-mixed to the PFR.  For a better illustration, a case study for different γ * and τ * values was conducted. A simple test case representing a single fluid element was used. The employed chemical mechanism for all calculations was the GRI3.0, which is an optimized methane combustion mechanism [32]. In order to investigate the influence of the reactor type, the ratio of the PFR and PSR mean reaction rate is taken. If the initial values are identical, the ratio reduces to the ratio of the difference between the mean and fine structure mass fractions: This ratio is investigated for the main combustion educts and products (CH 4 , CO, CO 2 , H 2 , H 2 O, O 2 ) and the main combustion radical species (H, O, OH) to quantify the effect of the fine structure treatment on the mean reaction rate.
Three sets of initial mole fractions were employed for the simple test case: lean and stoichiometric mixtures of main species, and a high radical mixture representing the reaction zone of a methane flame. Initial temperatures were set to 1000 K, and 2000 K for the main species and to 2000 K, and 2200 K for the high radical mixture, respectively. The fine structure residence times were chosen rather arbitrarily, but should reflect typical residence times. A PSR model based on Cantera [30] was employed for the test cases. Tables 1 and 2 summarize the initial conditions of the main species and radical species test cases, respectively.  Figures 6 and 7 show representative results of the total 192 data sets. In general, the differences in the mean reaction rates increase at bigger fine structure scales (γ * ), while they decrease for longer fine structure residence times (τ * ). For γ * = 0.1, mean reaction rates of most species are in the same range. In other words, PSR and PFR approaches give similar results for the mean rates in the case of small fine structure length fractions. Since classical turbulent combustion occurs in small isolated structures [3,12,13], both fine fraction treatments should give similar results. In contrast, in case of MILD combustion, reactions occur in a significant portion of the fluid and, therefore, the PFR approach will give different EDC reaction rates than the PSR approach.
Comparing the main and radical test case results reveals that the differences in the mean reaction rates are in general higher for the radical species cases. This is caused by the pronounced radical pool and the resulting small characteristic chemical time scales. The results confirm the findings of De et al. [17] that both approaches give similar results only for slow and slightly exothermic reactions. The differences in the reaction rates between the PSR and PFR fine structure treatment depend on the fine structure properties (γ L , τ * ), the fine structure state (temperature, pressure, chemical composition), and the employed reaction mechanism. However, under certain conditions, e.g., small fine structure shares, both reactor types can give comparable results.

Influences of the Model Variations in Flame Modeling
The influences of the proposed model variations are tested using Sandia Flame D, piloted under a stoichiometric premixed jet flame [22,23] surrounded by an air co-flow. The jet Reynolds number is approx. 22,400 and the visible flame length approx. 0.48 m. Table 3 shows the properties of the three gas streams (methane jet, pilot and co-flow stream), which are used as input parameters for the Flame D simulations. The flame is modeled using OpenFOAM R [24] and in-house libraries for the EDC model. Differential diffusion was neglected, since turbulent diffusion is dominant at high turbulence conditions. Although the modeling results are compared to the measured flame data, the main focus is on identifying the model differences. Therefore, a relatively coarse wedge grid with approx. 6090 cells is used. It has a 0.1 m upstream part for the methane jet and pilot inlet; the downstream part is 0.5 m in length and 0.15 m in radial direction. All five EDC modifications were tested on this grid using the PSR and PFR fine structure treatment, summing up ten test cases for the evaluation of the model differences and the fine structure treatment. The EDC2016 and EDC2017 were included in these tests to verify their validity for classical combustion conditions. Obviously, there should be no differences to the EDC1996, since the EDC2016 and EDC2017 are exactly the same as the EDC1996, if the model constants are limited to their standard values in the high turbulent regions. The chemical time scale (τ chem ) of Ren and Goldin [33] is employed to determine the EDC model constants in the EDC2016 and EDC2017 modification: Figure 8 compares the main species (CH 4 , O 2 , CO 2 , and H 2 O) mass fractions along the center line. In general, all variations give similar results except for the original EDC employing the PFR fine structure treatment, which did not ignite the mixture. Although no detailed investigation was conducted, the EDC reaction rate conversion factor, which is the lowest of all the EDC versions, might be the reason. Moreover, the other results show reasonable agreement to the experimental data. Due to almost identical results, lines in the subsequent diagrams overlap and only some are visible.  [22]. Not all EDC versions are visible due to overlapping profiles. Figure 9 shows the species mass fraction profiles of intermediate species (CO, H 2 , and OH) as well as the temperature profile along the center line. Compared to the main species, results for intermediates deviate more among the EDC modifications. However, only the OH radical profiles show distinct differences between the PSR and PFR fine structure treatment. The PFR OH peak mass fraction is approx. 1.5 times the PSR peak mass fraction, while the peak position and the profile shapes are similar. In this region, the PSR treatment is closer to the measured values. Temperature profiles confirm that the EDC1981 employing the PFR fine structure treatment failed to ignite the methane jet; therefore, this variant is excluded from further investigations. The other profiles resemble each other and show reasonable agreement with the measured temperatures. For a better quantification of the fine structure treatment influence, the PFR results are normalized using the PSR results of the corresponding EDC modification. In contrast to the species mass fraction profiles, these normalized profiles reveal more pronounced deviations between the PSR and PFR treatment. Especially for low species mass fractions, the PFR approach predicts up to several times higher or lower concentrations for consumed and produced species, respectively. Methane (CH 4 ) and oxygen (O 2 ) deviate the most 20 x/D upstream from the peak temperature, which is located at x/D ∼ 45. The PFR first predicts higher mass fractions than the PSR followed by lower ones (see Figure 10).
Carbon dioxide (CO 2 ) and water vapor (H 2 O) mass fractions are higher for the PFR downstream from the burner nozzle, but similar further downstream from the PSR values. In contrast, CO and H 2 mass fractions are higher for the PSR approach upstream from the peak temperature, while PFR and PSR results are similar further downstream (see Figure 11). OH radical mass fractions are significantly overpredicted by the PFR treatment compared to the PSR treatment directly above the burner tip (located at x/D = 0) and around the peak temperature location. Both give similar mass fractions further downstream, while OH mass fractions are higher for the PSR upstream. The deviations of the different profiles in the inlet and outlet regions of the intermediate species and CH 4 profiles are caused by numerical fluctuations for these amounts which are close to zero. The normalized temperature profile reveals that the PFR approach predicts slightly higher temperatures downstream from the burner tip. Starting at x/D ∼ 25, the PSR treatment predicts higher temperatures, but they converge further downstream. The maximum temperature deviation between the two fine structure treatments is around 10%. Comparing the fine structure share (γ * ) and the fine structure residence time (τ * ) along the centerline reveals that the values of the PFR and PSR treatment differ in the range between x/D ∼ 5 and x/D ∼ 45 for all EDC modifications (see Figure 12). Up to x/D ∼ 35, the fine structure share of the PFR computations is lower than using the PSR treatment. Further downstream, both fine fraction treatments predict similar fine structure shares. Assessing the species and relative species mass fraction profiles in this context shows that their deviations correlate with the variations in the fine structure shares (γ * ). Moreover, the EDC2005 fine structure shares are considerably higher than the ones of the other model variations, e.g., EDC1981 and EDC1996, which are virtually similar. That is obvious, since the fine structure share definition differs, e.g., Equation (10) for EDC2005 and Equation (3) for the other EDC versions. The EDC2016 and EDC2017 modifications give identical results as the EDC1996, thus the EDC1996 and EDC2016 values are overlapping in Figure 12. The fine structure deviations appear at the axial positions between x/D ∼ 5 and x/D ∼ 45, while the residence times deviate all over the investigated profile. Along the centerline, the normalized fine structure shares and residence times are identical for all EDC modifications, except the EDC2005 fine structure share, which deviates per definition. The fine structure residence times are similar for all investigated cases, although the cases employing the PFR treatment predict higher fine structure residence times than the PSR around x/D ∼ 40. The original EDC employing the PFR treatment (which failed to ignite the mixture) gives different residence times: lower ones up to x/D ∼ 55 and higher ones further downstream. For a better understanding of the variations in the PSR and PFR treatment, the species reaction rates along the center axis were investigated. Figure 13 shows the main species (e.g., CH 4 , O 2 , CO 2 , and H 2 O) rates for all EDC variations. The fluctuating rate profiles of the PSR treatment are probably caused by periodic oscillating states or numerical effects occurring when solving the PSR to steady-state. In contrast, the PFR rate profiles are smooth. The CH 4 , O 2 , and H 2 O consumption peaks of the PSR treatment are shifted approx. 5 x/D upstream compared to the PFR peaks. Moreover, the PSR peak magnitudes are almost twice as high as the PFR peak magnitudes. In contrast to the other main species rates, CO 2 is both produced and consumed along the flame D centerline. An initial production at x/D ∼ 30 is followed by a consumption zone and a second production zone between x/D ∼ 35 and 40 and x/D > 40, respectively. The predicted consumption peak and the second production peak are similar in position and magnitude for both the PSR and PFR fine structure treatment. The rates of the EDC1981 employing the PFR fine structure treatment shows species rates around zero, since the model failed to ignite the mixture.
The species consumption rate profiles of the intermediate species are given in Figure 14. The PSR fine structure treatment predicts CO and H 2 production peaks approx. 5 x/D upstream from the PFR rates. Peak magnitudes are about the same for CO and higher for H 2 in case of the PSR treatment.
Starting at x/D ∼ 40, CO and H 2 are consumed at almost similar rates for both fine structure treatments. In case of OH, the PFR fine structure treatment predicts a consumption peak shortly upstream that of the PSR treatment and has twice the magnitude. The PFR treatment gives smooth profiles, which shows an area of OH production between x/D ∼ 40 and 45 and a second consumption zone further downstream. After the initial production peak, the rates predicted by PSR treatment are alternating between consumption and production up to x/D ∼ 50, eventually reaching consumption far downstream. In general, the species consumption rates predicted by the PSR fine structure treatment indicate an earlier combustion start and higher rates than the PFR treatment. Furthermore, the CO, CO 2 and OH rate profiles project the inner and outer reaction zone occurring in Sandia Flame D. The PFR profiles are smooth, while the PSR ones are volatile in this investigation. The PSR behavior is likely caused by oscillating internal states during the solution procedure of the fine structures due to effects caused by the reaction system.

Conclusions
The focus of the present work is the comparison of the classical formulation of the Eddy Dissipation Concept (EDC) and its proposed modifications and variations in literature. A special focus is to extend the knowledge on the consequences of the plug flow reactor (PFR) simplification compared to the perfectly stirred reactor (PSR) fine structure treatment required by the EDC. Only schematic investigations of the PSR vs. PFR approach have been published in literature so far [16][17][18]. Thus, for the first time, a thorough evaluation of the consequences of the PFR approach is presented. The investigations are based on simple tests and a Sandia Flame D modeling [22,23] study. On top,

Conclusions
The focus of the present work is the comparison of the classical formulation of the Eddy Dissipation Concept (EDC) and its proposed modifications and variations in literature. A special focus is to extend the knowledge on the consequences of the plug flow reactor (PFR) simplification compared to the perfectly stirred reactor (PSR) fine structure treatment required by the EDC. Only schematic investigations of the PSR vs. PFR approach have been published in literature so far [16][17][18]. Thus, for the first time, a thorough evaluation of the consequences of the PFR approach is presented. The investigations are based on simple tests and a Sandia Flame D modeling [22,23] study. On top, the proposed EDC improvements [3][4][5][6][7][8] are discussed in detail and their influence on the EDC results is also evaluated. The findings of the present study are summarized in the following points: • The fine structures are treated as perfectly stirred (PSR) or plug flow (PFR) reactors in literature [4,6,16,17], although only the PSR treatment complies with the original EDC theory. A simple test series employing a single fluid element is used to study the fine structure treatment influence on the reaction rates. The results indicate increasing deviations between the PSR and PFR approaches for increasing fine structure shares (γ * ) and decreasing fine structure residence times (τ * ).
• Sandia Flame D [22,23] was modeled with an EDC extension for OpenFOAM R [24] to identify the influences of EDC variations with special focus on PSR and PFR fine structure treatment. The thorough investigation of the simulation results along the flame D center axis reveals that there is a difference between the PSR and PFR fine structure treatment for classical combustion conditions, e.g., high turbulence and fast chemistry (mixing controlled). The difference strongly depends on the local EDC parameters γ * and τ * as well as the characteristic chemical time scale. Nevertheless, taking into account the coarse grid, all investigated EDC variations (e.g., fine structure treatment, EDC formulation) give reasonable results. In general, the PFR treatment gives better results for CO and H 2 mass fraction profiles, while the PSR treatment gives better results for the OH mass fractions. The species mass fraction ratios reveal deviations between the PSR and PFR approach of more than 50% for the main and intermediate species, while temperatures deviate by up to 10%. Comparing the species consumption rates indicates an earlier combustion start and higher peak reaction rates for the PSR treatment. Both the fine structure share (γ * ) and the fine structure residence time (τ * ) vary up to 15%. Additionally, the consumption rate peaks are generally higher for the PSR approach. The original EDC employing the PFR approach failed to ignite Sandia Flame D in the simulation; the reasons were not investigated in detail.
The results indicate that the PFR fine structure treatment can give reasonable results for classical combustion conditions. Due to its significantly lower numerical effort compared to the PSR fine structure treatment, the PFR approach might be preferred for a priori investigations. However, a thorough selection should be done before choosing the fine structure treatment, if emphasis is put on detailed information of the reaction zone, e.g., intermediate species profiles.

Funding:
The authors gratefully acknowledge the funding support of K1-MET GmbH (Linz, Austria), metallurgical competence center. The research program of the competence center K1-MET is supported by COMET (Competence Center for Excellent Technologies), the Austrian program for competence centers. COMET is funded by the Federal Ministry for Transport, Innovation and Technology, the Federal Ministry for Science, Research and Economy, the province of Upper Austria, Tyrol, and Styria, the Styrian Business Promotion Agency.