40.5.4. The Discrete Ordinates (DO) Radiation Model

The discrete ordinates (DO) radiation model solves the radiative transfer equation (RTE) for a finite number of discrete solid angles, each associated with a vector direction fixed in the global Cartesian system (x, y, z). The DO model solves for as many transport equations as there are directions . The solution method is identical to that used for the fluid flow and energy equations.

The DO Model Equations

The DO model considers the radiative transfer equation (RTE) in the direction as a field equation. Thus, Equation 40–107 is written as

(40–107)

Angular Discretization and Pixelation

Each octant of the angular space 4 at any spatial location is discretized into N x N solid angles of extent, called control angles. The angles are the polar and azimuthal angles respectively, and are measured with respect to the global Cartesian system (x, y, z) as shown in Figure 40.2: The Angular Coordinate System. The theta and phi extents of the control angle, are constant. In two-dimensional calculations, only four octants are solved due to symmetry, making a total of directions in all. In three-dimensional calculations, a total of directions are solved. In the case of the non-gray model equations are solved for each band.

Figure 40.2: The Angular Coordinate System

The Angular Coordinate System

When Cartesian meshes are used, it is possible to align the global angular discretization with the control volume face, as shown in Figure 40.3: The Face with No Control Angle Overhang. For generalized unstructured meshes however, control volume faces do not in general align with the global angular discretization, as shown in Figure 40.4: The Face with Control Angle Overhang leading to the problem of control angle overhang.

Figure 40.3: The Face with No Control Angle Overhang

The Face with No Control Angle Overhang

Figure 40.4: The Face with Control Angle Overhang

The Face with Control Angle Overhang

Essentially, control angles can straddle the control volume faces, so that they are partially outgoing to the face. Figure 40.4: The Face with Control Angle Overhang shows a 3D example of a face with control angle overhang.

Figure 40.5: The Face with Control Angle Overhang (3D)

The Face with Control Angle Overhang (3D)

The control volume face cuts the sphere representing the angular space at an arbitrary angle. The line of intersection is a great circle. Control angle overhang may also occur as a result of reflection and refraction. It is important in these cases to correctly account for the overhanging fraction. This is done through the use of pixelation.

Each overhanging control angle is divided into pixels, as shown in Figure 40.6: The Pixelation of Control Angle

Figure 40.6: The Pixelation of Control Angle

The Pixelation of Control Angle

The energy contained in each pixel is then treated as incoming or outgoing to the face. The influence of overhang can thus be accounted for within the pixel resolution. Ansys Icepak allows you to choose the pixel resolution. For problems involving gray-diffuse radiation, the default pixelation of 1 X 1 is usually sufficient. For problems involving symmetry, periodic, specular, or semi-transparent boundaries, a pixelation of 3 X 3 is recommended. You should be aware, however that increasing the pixelation adds to the cost of computation.