## 13.3.8 Radiation in Combusting Flows

The Weighted-Sum-of-Gray-Gases Model

The weighted-sum-of-gray-gases model (WSGGM) is a reasonable compromise between the oversimplified gray gas model and a complete model which takes into account particular absorption bands. The basic assumption of the WSGGM is that the total emissivity over the distance can be presented as

 (13.3-81)

where are the emissivity weighting factors for the th fictitious gray gas, the bracketed quantity is the th fictitious gray gas emissivity, is the absorption coefficient of the th gray gas, is the sum of the partial pressures of all absorbing gases, and is the path length. For and FLUENT uses values obtained from [ 67] and [ 344]. These values depend on gas composition, and also depend on temperature. When the total pressure is not equal to 1 atm, scaling rules for are used (see Equation  13.3-87).

The absorption coefficient for is assigned a value of zero to account for windows in the spectrum between spectral regions of high absorption ( ) and the weighting factor for is evaluated from [ 344]:

 (13.3-82)

The temperature dependence of can be approximated by any function, but the most common approximation is

 (13.3-83)

where are the emissivity gas temperature polynomial coefficients. The coefficients and are estimated by fitting Equation  13.3-81 to the table of total emissivities, obtained experimentally [ 67, 79, 344].

The absorptivity of the radiation from the wall can be approximated in a similar way [ 344], but, to simplify the problem, it is assumed that  [ 247]. This assumption is justified unless the medium is optically thin and the wall temperature differs considerably from the gas temperature.

Since the coefficients and are slowly varying functions of and , they can be assumed constant for a wide range of these parameters. In [ 344] these constant coefficients are presented for different relative pressures of the CO and H O vapor, assuming that the total pressure is 1 atm. The values of the coefficients shown in [ 344] are valid for  atm-m and  K. For  K, coefficient values suggested by [ 67] are used. If for all , Equation  13.3-81 simplifies to

 (13.3-84)

Comparing Equation  13.3-84 with the gray gas model with absorption coefficient , it can be seen that the change of the radiation intensity over the distance in the WSGGM is exactly the same as in the gray gas model with the absorption coefficient

 (13.3-85)

which does not depend on . In the general case, is estimated as

 (13.3-86)

where the emissivity for the WSGGM is computed using Equation  13.3-81. as defined by Equation  13.3-86 depends on , reflecting the non-gray nature of the absorption of thermal radiation in molecular gases. In FLUENT, Equation  13.3-85 is used when m and Equation  13.3-86 is used for m. Note that for m, the values of predicted by Equations  13.3-85 and  13.3-86 are practically identical (since Equation  13.3-86 reduces to Equation  13.3-85 in the limit of small ).

FLUENT allows you to specify as the mean beam length or the characteristic cell size. The model based on the mean beam length is the recommended approach, especially when you have a nearly homogeneous medium and you are interested in the radiation exchange between the walls of the enclosure. You can specify the mean beam length or have FLUENT compute it. If you do decide to use the WSGGM based on the characteristic cell size, note that the predicted values of will be grid dependent (this is a known limitation of the model). See Section  8.8.1 for details about setting properties for the WSGGM.

 The WSGGM cannot be used to specify the absorption coefficient in each band when using the non-gray DO model. If the WSGGM is used with the non-gray DO model, the absorption coefficient will be the same in all bands.

When  atm

The WSGGM, as described above, assumes that --the total (static) gas pressure--is equal to 1 atm. In cases where is not unity (e.g., combustion at high temperatures), scaling rules suggested in [ 94] are used to introduce corrections. When  atm or  atm, the values for in Equations  13.3-81 and  13.3-85 are rescaled:

 (13.3-87)

where is a non-dimensional value obtained from [ 94], which depends on the partial pressures and temperature of the absorbing gases, as well as on .

The Effect of Soot on the Absorption Coefficient

When soot formation is computed, FLUENT can include the effect of the soot concentration on the radiation absorption coefficient. The generalized soot model estimates the effect of the soot on radiative heat transfer by determining an effective absorption coefficient for soot. The absorption coefficient of a mixture of soot and an absorbing (radiating) gas is then calculated as the sum of the absorption coefficients of pure gas and pure soot:

 (13.3-88)

where is the absorption coefficient of gas without soot (obtained from the WSGGM) and

 (13.3-89)

with

is the soot density in kg/m .

The coefficients and were obtained [ 318] by fitting Equation  13.3-89 to data based on the Taylor-Foster approximation  [ 370] and data based on the Smith et al. approximation [ 344].

See Sections  8.8 and 20.3.3 for information about including the soot-radiation interaction effects.

The Effect of Particles on the Absorption Coefficient

FLUENT can also include the effect of discrete phase particles on the radiation absorption coefficient, provided that you are using either the P-1 or the DO model. When the P-1 or DO model is active, radiation absorption by particles can be enabled. The particle emissivity, reflectivity, and scattering effects are then included in the calculation of the radiative heat transfer. See Section  22.14 for more details on the input of radiation properties for the discrete phase.