Closures for the chemical source term

The CFD model developed above is an example of a moment closure. Unfortunately, when applied to reacting scalars such as those considered in Section III, moment closures for the chemical source term are not usually accurate (Fox, 2003). An alternative approach that yields the same moments can be formulated in terms of a presumed PDF method (Fox, 1998). Here we will consider only the simplest version of a multi-environment micromixing model. Readers interested in further details on other versions of the model can consult Wang and Fox (2004). 

A theoretical framework based on the one-point, one-time joint probability density function (PDF) is developed. It is shown that all commonly employed models for turbulent reacting flows can be formulated in terms of the joint PDF of the chemical species and enthalpy. Models based on direct closures for the chemical source term as well as transported PDF methods, are covered in detail. An introduction to the theory of turbulence and turbulent scalar transport is provided for completeness. 

In other closures for the chemical source term, a model for the conditional scalar dissipation rate (e

scalar Laplacian, the conditional scalar... 

Another conditional expectation that frequently occurs in closures for the chemical source term is the conditional mean of the composition variables given the mixturefraction. The latter, defined in Chapter 5, is an inert scalar formed by taking a linear combination of the components of 0 ... 

In developing closures for the chemical source term and the PDF transport equation, we will also come across conditional moments of the derivatives of a field conditioned on the value of the field. For example, in conditional-moment closures, we must provide a functional form for the scalar dissipation rate conditioned on the mixture fraction, i.e.,... 

The turbulence models discussed in this chapter attempt to model the flow using low-order moments of the velocity and scalar fields. An alternative approach is to model the one-point joint velocity, composition PDF directly. For reacting flows, this offers the significant advantage of avoiding a closure for the chemical source term. However, the numerical methods needed to solve for the PDF are very different than those used in standard CFD codes. We will thus hold off the discussion of transported PDF methods until Chapters 6 and 7 after discussing closures for the chemical source term in Chapter 5 that can be used with RANS and LES models. 

Figure 5.1. Closures for the chemical source term can be understood in terms of their relationship to the joint composition PDF. The simplest methods attempt to represent the joint PDF by its (lower-order) moments. At the next level, the joint PDF is expressed in terms of the product of the conditional joint PDF and the mixture-fraction PDF. The conditional joint PDF can then be approximated by invoking the fast-chemistry or flamelet limits, by modeling the conditional means of the compositions, or by assuming a functional form for the PDF. Similarly, it is also possible to assume a functional form for the joint composition PDF. The best method to employ depends strongly on the functional form of the chemical source term and its characteristic time scales.

The simplest closure for the chemical source term is to assume that the joint composition PDF can be represented by its moments. In general, this assumption is of limited validity. Nevertheless, in this section we review methods based on moment closures in order to illustrate their limitations. 

In (5.297), the interpolation parameter is defined separately for each component. Note, however, that unlike the earlier examples, there is no guarantee that the interpolation parameters will be bounded between zero and one. For example, the equilibrium concentration of intermediate species may be negligible despite the fact that these species can be abundant in flows dominated by finite-rate chemistry. Thus, although (5.297) provides a convenient closure for the chemical source term, it is by no means guaranteed to produce accurate predictions A more reliable method for determining the conditional moments is the formulation of a transport equation that depends explicitly on turbulent transport and chemical reactions. We will look at this method for both homogeneous and inhomogeneous flows below. 

Hence, it is not necessary to solve a separate LES transport equation for the filtered compositions. Indeed, in the limit of one environment, (5.422) reduces to the LES model for the filtered compositions with the simplest possible closure for the chemical source term (i.e., one that neglects all SGS fluctuations). 

---