ECFM Models with LES

The ECFM-3Z and ECFM-CLEH models are both compatible with any of the LES turbulence models.

Using the ECFM-3Z or ECFM-CLEH model in conjunction with any LES turbulence model activates a new equation for modelling the subgrid scale (SGS) flame surface density.
Figure 1. EQUATION_DISPLAY
Σ t + T r e s T s g s + P = S r e s + S s g s + C r e s + C s g s
(3891)
T r e s is the resolved transport:
Figure 2. EQUATION_DISPLAY
T r e s = ∇⋅ ( v Σ )
(3892)
T s g s is the unresolved transport plus the molecular diffusion term:
Figure 3. EQUATION_DISPLAY
T s g s = ∇⋅ [ σ c ( D + v s g s S c t ) Σ ]
(3893)
P is the laminar propagation (where n is the unit vector normal to the flame front):
Figure 4. EQUATION_DISPLAY
P = ∇⋅ ( S d n Σ )
(3894)
S r e s is the resolved tangential strain:
Figure 5. EQUATION_DISPLAY
S r e s = ( ∇⋅ v n v n ) Σ
(3895)
C r e s is the resolved curvature:
Figure 6. EQUATION_DISPLAY
C r e s = S d ( ∇⋅ n ) Σ
(3896)
C s g s is the unresolved curvature:
Figure 7. EQUATION_DISPLAY
C s g s = β S l c * c ¯ c ¯ ( 1 c ¯ ) ( Σ Σ l a m ) Σ
(3897)
where β and c * are model constants. Σ l a m is the laminar part of the flame surface density:
Figure 8. EQUATION_DISPLAY
Σ l a m = | c ˜ | + ( c ¯ c ˜ ) ∇⋅ n
(3898)
S s g s is the unresolved tangential strain:
Figure 9. EQUATION_DISPLAY
S s g s = Γ { u ^ S l , Δ ^ δ l } u ^ Δ ^ Σ σ c
(3899)
which is a function of the non-dimensional parameters u ^ S l , Δ ^ δ l where the ^ symbol indicates the combustion filter size:
Figure 10. EQUATION_DISPLAY
Δ ^ = N r e s Δ
(3900)
the species, enthalpy and flame surface density equations are filtered at this filter size inside the flame. The Intermittent Turbulent Net Flame Stretch (ITNFS) function Γ in Eqn. (3899), models the effect of flame wrinkling by all vortices and is a function of the unresolved turbulence scales. Colin et al. [810] showed that vortices smaller than the turbulent flame brush l t are unable to wrinkle the flame. Therefore, it is necessary that the input values for Γ are filtered with a filter greater than or equal to l t . In other words, Δ l t .
Near the walls, the flames are quenched by q w a l l :
Figure 11. EQUATION_DISPLAY
q w a l l = q w f v
(3901)
where q w is the wall quenching factor and f v is the Van Driest damping function (see Eqn. (1389)).

For ECFM-3Z with LES the treatment at the wall is to multiply all the production terms (i.e. terms with a positive value on the RHS) in equation Eqn. (3891) by q w a l l .