# 3.2.5. Turbulent Flow Equations, Favre-Averaged

The Favre-averaged turbulent transport equations are derived from the laminar equations of Laminar Flow Equations by passing the equations through either the RANS temporal filter of (3.68) or the LES spatial filter of (3.72). The mathematical form of the equations are essentially identical between the two filtering methods, so only a single set of equations will be presented. Care should be taken to interpret the filters as either temporal or spatial, depending on the closure models selected. While it is the Favre-averaged form of the equations that are solved, a comparison of the simple Reynolds-averaged and the Favre-averaged form is given at the end of this page for reference.

The approach most commonly used in turbulence modeling is called the Boussinesq eddy viscosity approximation, which relates the turbulent stress tensor to the filtered strain rate tensor through a modeled turbulent eddy viscosity. This general modeling approach has shown remarkable success for a broad range of problems (Wilcox [13]), and is the approach used in SIERRA/Fuego. A similar approach is used for scalar transport, where the scalar flux vector is related to scalar gradients through a modeled diffusion coefficient.

The following subsections describe the turbulent transport equations expressed in terms of a turbulent eddy viscosity or turbulent diffusion coefficient through the Boussinesq approximation. The treatment of these coefficients is dependent upon which of the many closure models are selected, and will be described in Turbulence Closure Models.

## 3.2.5.1. Conservation of Mass

The integral form of the Favre-filtered continuity equation used for turbulent transport is

(3.77)

This equation is in closed form, and no additional modeling is required.

## 3.2.5.2. Conservation of Momentum

The integral form of the Favre-filtered momentum equations used for turbulent transport are

(3.78)

where the turbulent stress is defined as

(3.79)

### 3.2.5.2.1. RANS Modeling

For RANS simulations, represents the Reynolds stress tensor and can be reduced to the form by substitution of the Favre decomposition of each variable and simplifying. The deviatoric (trace-free) part of the stress tensor is defined as

(3.80)

(3.81)

where the turbulent kinetic energy is defined as . The deviatoric part of the Reynolds stress tensor is modeled by the Boussinesq approximation which relates the Reynolds stresses to the filtered strain rate tensor through a modeled turbulent viscosity , resulting in

(3.82)

(3.83)

where the filtered strain rate tensor is defined by

(3.84)

Substituting this into (3.81) yields the modeled form of the full Reynolds stress tensor (Kuo [12], p. 445)

(3.85)

The Favre-filtered momentum equations then become

(3.86)

(3.87)

where RANS closure models for the turbulent viscosity are presented in Turbulence Closure Models.

### 3.2.5.2.2. LES Modeling

For LES, in (3.78) represents the subgrid stress tensor. The deviatoric part of the subgrid stress tensor is defined as

(3.88)

(3.89)

where the subgrid turbulent kinetic energy is defined as . The deviatoric part of the subgrid stress tensor is then modeled similar to RANS closures as (Moin, et al. [14])

(3.90)

Substituting this into (3.89) yields the modeled form of the full subgrid stress tensor

(3.91)

For low Mach-number flows, a vast majority of the turbulent kinetic energy is contained at resolved scales (Erlebacher, et al. [15]). For this reason, the subgrid turbulent kinetic energy will not be directly treated and will instead be included in the pressure as an additional normal stress. The Favre-filtered momentum equations then become

(3.92)

(3.93)

where LES closure models for the subgrid turbulent eddy viscosity are presented in Turbulence Closure Models.

## 3.2.5.3. Conservation of Energy

The integral form of the Favre-filtered energy equation used for turbulent transport is

(3.94)

(3.95)

The simple Fickian diffusion velocity approximation, (3.46), is assumed, so that the mean diffusive heat flux vector is

(3.96)

If Sc = Pr, i.e., unity Lewis number (Le = 1), then the diffusive heat flux vector simplifies to . The viscous dissipation term is closed by

(3.97)

The turbulent diffusive flux vector in (3.95) is defined as

(3.98)

For RANS simulations, represents the turbulent energy diffusive flux vector and is simplified to the form by substitution of the Favre decomposition of each variable. It is then modeled by

(3.99)

where is the turbulent Prandtl number and is the modeled turbulent eddy viscosity from momentum closure. For LES, represents the subgrid turbulent energy diffusive flux vector, and is modeled in the same way as

(3.100)

where is the subgrid turbulent Prandtl number and is the modeled subgrid turbulent eddy viscosity from momentum closure.

The resulting filtered and modeled turbulent energy equation for both RANS and LES is given in Libby and Williams [11], p. 25, as

(3.101)

This equation is also given in Gran et al. [16] (without the transient and radiation source terms and the additional term for laminar transport). The turbulent Prandtl number must have the same value as the turbulent Schmidt number for species transport to maintain unity Lewis number.

## 3.2.5.4. Conservation of Species

The integral form of the Favre-filtered species equation used for turbulent transport is

(3.102)

where the form of diffusion velocities (see (3.46)) assumes the Fickian approximation with a constant value of diffusion velocity for consistency with the turbulent form of the energy equation, (3.95).

The turbulent diffusive flux vector is defined as

(3.103)

For RANS simulations, represents the turbulent species diffusive flux vector and is simplified to the form by substitution of the Favre decomposition of each variable. It is then modeled as

(3.104)

where is the turbulent Schmidt number for all species and is the modeled turbulent eddy viscosity from momentum closure. For LES, represents the subgrid turbulent species diffusive flux vector, and is modeled identically as

(3.105)

where is the subgrid turbulent Schmidt number for all species and is the subgrid modeled turbulent eddy viscosity from momentum closure.

The Favre-filtered and modeled turbulent species transport equation for both RANS and LES then becomes (Gran et al. [16])

(3.106)

If transporting both energy and species equations, the laminar Prandtl number must be equal to the laminar Schmidt number and the turbulent Prandtl number must be equal to the turbulent Schmidt number to maintain unity Lewis number. Although there is a species conservation equation for each species in a mixture of species, only species equations need to be solved since the mass fractions sum to unity and

(3.107)

## 3.2.5.5. Radiation Transport

The Favre-averaged energy equation, (3.101), requires the time-averaged radiative flux divergence. From (3.63), the time-averaged radiative flux divergence is given by

(3.108)

For optically thin turbulent eddies, which is the case for many combustion applications, fluctuations in the absorption coefficient and the scalar flux are weakly correlated [7] so (3.108) may be simplified to

(3.109)

The time averaged scalar flux is obtained from the time averaged Boltzmann radiative transport equation

(3.110)

where the correlation between the turbulent fluctuations in the absorption coefficient and the intensity is assumed small to simplify the absorption term.

Both (3.109) and (3.110) include the time averaged emission term, , which may significantly increase the radiative emission from a turbulent flame above what would be estimated from the mean temperature and absorption coefficient values. The details of the closure used for this term are discussed in the turbulent combustion model section.

## 3.2.5.6. Comparison Between Time-Averaging and Favre-Averaging

The time-averaged (Reynolds-averaged) and Favre-averaged transport equations are given in the following section.

### 3.2.5.6.1. Conservation of Mass

The continuity equation:

> - Time averaged:

>   (3.111)

> - Favre averaged:

>   (3.112)

### 3.2.5.6.2. Conservation of Momentum

The momentum transport equations:

> - Time averaged:

>   (3.113)

> - Favre averaged:

>   (3.114)

### 3.2.5.6.3. Conservation of Energy

The energy transport equation, assume Lewis number is one:

> - Favre averaged:

>   (3.115)

### 3.2.5.6.4. Conservation of Species

The species transport equation:

> - Time averaged:

>   > (3.116)

> > - Favre averaged:

> >   (3.117)