## 13.3.7 Surface-to-Surface (S2S) Radiation Model Theory

The surface-to-surface radiation model can be used to account for the radiation exchange in an enclosure of gray-diffuse surfaces. The energy exchange between two surfaces depends in part on their size, separation distance, and orientation. These parameters are accounted for by a geometric function called a "view factor''.

The main assumption of the S2S model is that any absorption, emission, or scattering of radiation can be ignored; therefore, only "surface-to-surface'' radiation need be considered for analysis.

FLUENT's S2S radiation model assumes the surfaces to be gray and diffuse. Emissivity and absorptivity of a gray surface are independent of the wavelength. Also, by Kirchoff's law [ 248], the emissivity equals the absorptivity ( ). For a diffuse surface, the reflectivity is independent of the outgoing (or incoming) directions.

The gray-diffuse model is what is used in FLUENT. Also, as stated earlier, for applications of interest, the exchange of radiative energy between surfaces is virtually unaffected by the medium that separates them. Thus, according to the gray-body model, if a certain amount of radiant energy ( ) is incident on a surface, a fraction ( ) is reflected, a fraction ( ) is absorbed, and a fraction ( ) is transmitted. Since for most applications the surfaces in question are opaque to thermal radiation (in the infrared spectrum), the surfaces can be considered opaque. The transmissivity, therefore, can be neglected. It follows, from conservation of energy, that , since (emissivity), and .

The S2S Model Equations

The energy flux leaving a given surface is composed of directly emitted and reflected energy. The reflected energy flux is dependent on the incident energy flux from the surroundings, which then can be expressed in terms of the energy flux leaving all other surfaces. The energy reflected from surface is

 (13.3-70)

where is the energy flux leaving the surface, is the emissivity, is Boltzmann's constant, and is the energy flux incident on the surface from the surroundings.

The amount of incident energy upon a surface from another surface is a direct function of the surface-to-surface "view factor,'' . The view factor is the fraction of energy leaving surface that is incident on surface . The incident energy flux can be expressed in terms of the energy flux leaving all other surfaces as

 (13.3-71)

where is the area of surface and is the view factor between surface and surface . For surfaces, using the view factor reciprocity relationship gives

 (13.3-72)

so that

 (13.3-73)

Therefore,

 (13.3-74)

which can be written as

 (13.3-75)

where represents the energy that is given off (or radiosity) of surface , and represents the emissive power of surface . This represents equations, which can be recast into matrix form as

 (13.3-76)

where is an matrix, is the radiosity vector, and is the emissive power vector.

Equation  13.3-76 is referred to as the radiosity matrix equation. The view factor between two finite surfaces and is given by

 (13.3-77)

where is determined by the visibility of to . = 1 if is visible to and 0 otherwise.

Clustering

The S2S radiation model is computationally very expensive when there are a large number of radiating surfaces. To reduce the computational time as well as the storage requirement, the number of radiating surfaces is reduced by creating surface "clusters''. The surface clusters are made by starting from a face and adding its neighbors and their neighbors until a specified number of faces per surface cluster is collected.

A new algorithm has been implemented for the creation of surface clusters which is faster and supports non-conformal interfaces, hanging nodes, or grid adaption. This algorithm is now the default. If you wish to use the old algorithm, you may use the TUI command but adaption and non-conformal interfaces will not be supported.

The radiosity, , is calculated for the surface clusters. These values are then distributed to the faces in the clusters to calculate the wall temperatures. Since the radiation source terms are highly non-linear (proportional to the fourth power of temperature), care must be taken to calculate the average temperature of the surface clusters and distribute the flux and source terms appropriately among the faces forming the clusters.

The surface cluster temperature is obtained by area averaging as shown in the following equation:

 (13.3-78)

where is the temperature of the surface cluster, and and are the area and temperature of face . The summation is carried over all faces of a surface cluster.

Smoothing

Smoothing can be performed on the view factor matrix to enforce the reciprocity relationship and conservation.

The reciprocity relationship is represented by

 (13.3-79)

where is the area of surface , is the view factor between surfaces and , and is the view factor between surfaces and .

Once the reciprocity relationship has been enforced, a least-squares smoothing method [ 191] can be used to ensure that conservation is satisfied, i.e.,

 (13.3-80)

Previous: 13.3.6 Discrete Ordinates (DO)