The key to the premixed combustion model is the prediction of , the turbulent flame speed normal to the mean surface of the flame. The turbulent flame speed is influenced by the following:
In FLUENT, the turbulent flame speed is computed using a model for wrinkled and thickened flame fronts [ 415]:
|=||RMS (root-mean-square) velocity (m/s)|
|=||laminar flame speed (m/s)|
|=||molecular heat transfer coefficient of unburnt|
|mixture (thermal diffusivity) (m /s)|
|=||turbulence length scale (m)|
|=||turbulence time scale (s)|
|=||chemical time scale (s)|
The turbulence length scale, , is computed from
where is the turbulence dissipation rate.
The model is based on the assumption of equilibrium small-scale turbulence inside the laminar flame, resulting in a turbulent flame speed expression that is purely in terms of the large-scale turbulent parameters. The default value of 0.52 for is recommended by [ 415], and is suitable for most premixed flames. The default value of 0.37 for should also be suitable for most premixed flames.
The model is strictly applicable when the smallest turbulent eddies in the flow (the Kolmogorov scales) are smaller than the flame thickness, and penetrate into the flame zone. This is called the thin reaction zone combustion region, and can be quantified by Karlovitz numbers, Ka, greater than unity. Ka is defined as
|=||characteristic flame time scale|
|=||smallest (Kolmogorov) turbulence time scale|
Lastly, the model is valid for premixed systems where the flame brush width increases in time, as occurs in most industrial combustors. Flames that propagate for a long period of time equilibrate to a constant flame width, which cannot be captured in this model.
Turbulent Flame Speed for LES
For simulations that use the LES turbulence model, the Reynolds-averaged quantities in the turbulent flame speed expression (Equation 16.2-4) are replaced by their equivalent sub-grid quantities. In particular, the large eddy length scale is modeled as
where is the Smagorinsky constant and is the cell characteristic length.
The RMS velocity in Equation 16.2-4 is replaced by the sub-grid velocity fluctuation, calculated as
where is the sub-grid scale mixing rate (inverse of the sub-grid scale time scale), given in Equation 14.1-28.
Laminar Flame Speed
The laminar flame speed ( in Equation 16.2.1) can be specified as constant, or as a user-defined function. A third option appears for non-adiabatic premixed and partially-premixed flames and is based on the correlation proposed by Meghalchi and Keck [ 241],
In Equation 16.2-10, and are the unburnt reactant temperature and pressure ahead of the flame, and .
The reference laminar flame speed, , is calculated from
where is the equivalence ratio ahead of the flame front, and , and are fuel-specific constants. The exponents and are calculated from,
The Meghalchi-Keck laminar flame speeds are available for fuel-air mixtures of methane, methanol, propane, iso-octane and indolene fuels.
Unburnt Density and Thermal Diffusivity
The unburnt density ( in Equation 16.2.1) and unburnt thermal diffusivity ( in Equation 16.2-5) are specified constants that are set in the Materials panel. However, for compressible cases, such as in-cylinder combustion, these can change significantly in time and/or space. When the ideal gas model is selected for density, the unburnt density and thermal diffusivity are calculated as volume averages ahead of the flame front.
Flame Stretch Effect
Since industrial low-emission combustors often operate near lean blow-off, flame stretching will have a significant effect on the mean turbulent heat release intensity. To take this flame stretching into account, the source term for the progress variable ( in Equation 16.2-1) is multiplied by a stretch factor, [ 417]. This stretch factor represents the probability that the stretching will not quench the flame; if there is no stretching ( ), the probability that the flame will be unquenched is 100%.
The stretch factor, , is obtained by integrating the log-normal distribution of the turbulence dissipation rate, :
where erfc is the complementary error function, and and are defined below.
is the standard deviation of the distribution of :
where is the stretch factor coefficient for dissipation pulsation, is the turbulent integral length scale, and is the Kolmogorov micro-scale. The default value of 0.26 for (measured in turbulent non-reacting flows) is recommended by [ 415], and is suitable for most premixed flames.
is the turbulence dissipation rate at the critical rate of strain [ 415]:
By default, is set to a very high value ( ) so no flame stretching occurs. To include flame stretching effects, the critical rate of strain should be adjusted based on experimental data for the burner. Numerical models can suggest a range of physically plausible values [ 415], or an appropriate value can be determined from experimental data. A reasonable model for the critical rate of strain is
where is a constant (typically 0.5) and is the thermal diffusivity. Equation 16.2-16 can be implemented in FLUENT using a property user-defined function. See the separate UDF Manual. for details about user-defined functions.
Preferential diffusion accounts for the effect of variations in fuel molecular diffusion coefficients on heat release intensity in premixed turbulent combustion. Inclusion of this effect is important for simulation of combustion with light fuels (e.g., hydrogen) or heavy fuels (e.g., evaporated oil). The model for preferential diffusion is based on the concept of leading points, formulated in [ 189]. The authors of [ 189] derived formulas for the changes in mixture composition within the combustion zone due to the difference in the molecular diffusivities of fuel, , and oxidizer, . These formulas are rewritten in [ 417] as
|=||mass stoichiometric coefficient|
|=||stoichiometric ratio of unburnt mixture composition|
|=||stoichiometric ratio of leading-point composition|
The concept of leading points is applied to the FLUENT model by using instead of for the formulation of the laminar flame speed, , or the molecular heat transfer coefficient, . This simple approach results in reasonable agreement with the measurements of mass combustion rates in stirred bombs [ 417], without the need for additional empirical parameters.
Volume expansion at the flame front can cause counter-gradient diffusion. This effect becomes more pronounced when the ratio of the reactant density to the product density is large, and the turbulence intensity is small. It can be quantified by the ratio , where , , , and are the unburnt and burnt densities, laminar flame speed, and turbulence intensity, respectively. Values of this ratio greater than one indicate a tendency for counter-gradient diffusion, and the premixed combustion model may be inappropriate. Recent arguments for the validity of the turbulent-flame-speed model in such regimes can be found in Zimont et al. [ 416].