## 13.3.4 Rosseland Radiation Model Theory

The Rosseland or diffusion approximation for radiation is valid when the medium is optically thick ( ), and is recommended for use in problems where the optical thickness is greater than 3. It can be derived from the P-1 model equations, with some approximations. This section provides details about the equations used in the Rosseland model.

The Rosseland Model Equations

As with the P-1 model, the radiative heat flux vector in a gray medium can be approximated by Equation  13.3-4:

 (13.3-22)

where is given by Equation  13.3-3.

The Rosseland radiation model differs from the P-1 model in that the Rosseland model assumes that the intensity is the black-body intensity at the gas temperature. (The P-1 model actually calculates a transport equation for .) Thus , where is the refractive index. Substituting this value for into Equation  13.3-22 yields

 (13.3-23)

Since the radiative heat flux has the same form as the Fourier conduction law, it is possible to write

 (13.3-24) (13.3-25) (13.3-26)

where is the thermal conductivity and is the radiative conductivity. Equation  13.3-24 is used in the energy equation to compute the temperature field.

Anisotropic Scattering

The Rosseland model allows for anisotropic scattering, using the same phase function (Equation  13.3-7) described for the P-1 model in Section  13.3.3.

Boundary Condition Treatment for the Rosseland Model at Walls

Since the diffusion approximation is not valid near walls, it is necessary to use a temperature slip boundary condition. The radiative heat flux at the wall boundary, , is defined using the slip coefficient :

 (13.3-27)

where is the wall temperature, is the temperature of the gas at the wall, and the slip coefficient is approximated by a curve fit to the plot given in [ 332]:

 (13.3-28)

where is the conduction to radiation parameter at the wall:

 (13.3-29)

and .

Boundary Condition Treatment for the Rosseland Model at Flow Inlets and Exits

No special treatment is required at flow inlets and outlets for the Rosseland model. The radiative heat flux at these boundaries can be determined using Equation  13.3-24.