2.11.5 Viewfactor calculation

Product: ABAQUS/Standard

Cavity radiation occurs when surfaces of the model can see each other and, thus, exchange heat with each other by radiation (Figure 2.11.5–1). Such exchange depends on viewfactors that measure the relative interaction between the surfaces composing the cavity. Viewfactor calculation is rather complicated for anything but the most trivial geometries. ABAQUS offers an automatic viewfactor calculation capability for two- and three-dimensional cases as well as for axisymmetric situations. This capability can take into account general surface blocking (or shadowing) as well as the most common forms of radiation symmetry. The viewfactor calculation can also be automatically repeated a number of times throughout the analysis history (this is user-controlled) if cavity surfaces are moved in space causing the viewfactors to change.

Figure 2.11.5–1 Heat exchange between surfaces by radiation.

ABAQUS, Inc., is pleased to acknowledge that the viewfactor calculation technology implemented in ABAQUS was derived from technology originally developed by the Atomic Energy Authority of the United Kingdom; see, for example, Johnson (1987). The remainder of this section contains a general description of this technology.

Viewfactor between elementary areas

The viewfactor between two elementary areas, and , can be generally written as

where R is the distance between the two areas and

are the angles between R and the normals to the surfaces of the areas (Figure 2.11.5–2).

Figure 2.11.5–2 Schematic of viewfactor calculation.

This formula applies when the areas

and

are small compared with the distance R. We compute viewfactors in an area-lumped fashion, so that the integral of Equation 2.11.5–1 would become

However, this expression is singular in

, so it is not well behaved for general surfaces. Thus, we use

The limit of Equation 2.11.5–3 for large

is Equation 2.11.5–2.

Discretization

Cavities are composed of surfaces; and surfaces, in turn, are made up of finite element faces. For the purpose of viewfactor calculations, one can then think of a cavity as a collection of element faces (or facets) corresponding to the finite element discretization around the cavity.In the two- and three-dimensional cases the element faces composing cavities can be treated as elementary areas and, thus, Equation 2.11.5–1 applies. In axisymmetric cases the element faces represent rings so that the viewfactors involve two ring surfaces looking at each other. This requires an integration over where it is important to account for “horizon” effects (Johnson, 1987).

In so far as the viewfactor calculations are concerned, first- and second-order element faces are treated similarly in the sense that the midside nodes of the faces in the second-order elements are ignored. This means that a pair of four-noded faces looking at each other will produce the same viewfactor as a pair of eight-noded faces with corner nodes coinciding with the nodes of the four-noded faces.

Radiation blocking

Radiation within a cavity implies that every surface exchanges heat with every other surface. The problem is made more complex when solid bodies are interposed between radiating surfaces blocking (or shadowing) off some but not all the possible paths along which heat can be radiated from the facets of one surface to the facets of another surface (Figure 2.11.5–3).

Figure 2.11.5–3 Blocking or shadowing example.

It is inconceivable that the user could handle this complexity in all but the simplest situations. Therefore, by default, ABAQUS automatically checks if blocking takes place for every possible radiation path in a cavity. This requires that the program check if the ray joining the centers of each pair of facets intersects any other facet. For cavities with a large number of facets this can be very time consuming. For this reason ABAQUS allows the user to guide its blocking algorithm by accepting input of which surfaces cause blocking, thus significantly reducing the computational effort required. If a ray between two facets intersects any other facet, then in the two- and three-dimensional cases that ray is eliminated and no radiative heat transfer takes place between the facets. In the axisymmetric case blocking is much more complicated since each element face in the finite element model represents a ring. This is handled automatically and requires that the program calculate which part of the extent of the ring is blocked.

Radiation symmetries

Use of symmetry can greatly reduce the size of a problem, but—in the case of cavity radiation—it requires that special facilities for definition and handling of symmetries be available. ABAQUS provides capabilities for three different kinds of symmetries: simple reflection symmetry, periodic symmetry, and cyclic symmetry. Reflection symmetry allows one additional image of the model to be created by reflection through a line in two dimensions or reflection through a plane in three dimensions. Periodic symmetry can be used to create multiple images of the model by periodic repetition in two- or three-dimensional space according to a periodic distance vector. Cyclic symmetry creates multiple images of the model by cyclic repetition about a point in two dimensions or by cyclic repetition about an axis in three dimensions. Combinations of the different types of symmetry are supported.

To illustrate the handling of symmetries during viewfactor calculation, consider the case of a simple reflection symmetry in two-dimensional space (Figure 2.11.5–4).

Figure 2.11.5–4 Reflection symmetry example.

Radiation between facet i (with its centroid at point A) and facet j (with its centroid at point B) has two contributions: one arising from the ray between points A and B and the other coming from the ray between points A and

, where

is the mirror image of B. The length of ray

is defined directly in the model. The definition of the length of ray

requires that point C on the reflection symmetry line be located such that

and

make equal angles to it; ray

then has length

. Similar logic can be extended to the three-dimensional case.

In axisymmetric cases symmetry about the axis of symmetry of the model is always implied, and the only other symmetries allowed are simple reflection through a plane normal to the axis of symmetry or periodic repetition in the direction of the axis of symmetry.

Viewfactor checking

The viewfactor is a purely geometrical quantity, and it has some special properties. One property that allows us to check the accuracy of the calculation is that for a completely enclosed cavity:

indicating that any ray from surface i in whatever direction it leaves the surface will reach another surface in the enclosed cavity.

The quantity in Equation 2.11.5–4 is calculated for every facet of each cavity, and its value is used to provide a check to control the accuracy of viewfactor calculation.

Radiation to ambient

The quantity calculated in Equation 2.11.5–4 can deviate from unity so long as the cavity is not fully enclosed. The user can define such an open cavity by giving the value of the ambient temperature in the cavity definition. In this case the difference between one and the quantity calculated in Equation 2.11.5–4 for each facet of the open cavity is considered to be the fraction radiating from that facet to the surrounding medium.