Heat transfer

In a heat transfer analysis, triggered by the *HEAT TRANSFER procedure card, the temperature is the independent degree of freedom. In essence, the energy equation is solved subject to temperature and flux boundary conditions ([19]). For steady-state calculations it leads to a Laplace-type equation.

The governing equation for heat transfer reads

$\displaystyle \nabla \cdot (- \boldsymbol{\kappa} \cdot \nabla{T}) + \rho c \dot{T} = \rho h$ (357)

where $ \boldsymbol{\kappa}$ contains the conduction coefficients, $ \rho$ is the density, $ h$ the heat generation per unit of mass and $ c$ is the specific heat.

The temperature can be defined using the *BOUNDARY card using degree of freedom 11. Flux type boundary conditions can consist of any combination of the following:

  1. Concentrated flux, applied to nodes, using the *CFLUX card (degree of freedom 11)
  2. Distributed flux, applied to surfaces or volumes, using the *DFLUX card
  3. Convective flux defined by a *FILM card. It satisfies the equation

    $\displaystyle q = h(T-T_0)$ (358)

    where $ q$ is the a flux normal to the surface, $ h$ is the film coefficient, $ T$ is the body temperature and $ T_0$ is the environment temperature (also called sink temperature). CalculiX can also be used for forced convection calculations, in which the sink temperature is an unknown too. This applied to all kinds of surfaces cooled by fluids or gases.

  4. Radiative flux defined by a *RADIATE card. The equation reads

    $\displaystyle q = \epsilon (\theta^4-\theta_0^4)$ (359)

    where $ q$ is a flux normal to the surface, $ \epsilon$ is the emissivity, $ \theta$ is the absolute body temperature (Kelvin) and $ \theta_0$ is the absolute environment temperature (also called sink temperature). The emissivity takes values between 0 and 1. A zero value applied to a body with no absorption nor emission and 100 $ \%$ reflection. A value of 1 applies to a black body. The radiation is assumed to be diffuse (independent of the direction of emission) and gray (independent of the emitted wave length). CalculiX can also be used for cavity radiation, simulating the radiation interaction of several surfaces. In that case, the viewfactors are calculated, see also [35] for the fundamentals of heat transfer and [6] for the calculation of viewfactors.

    The calculation of viewfactors involves the solution of a four-fold integral. By using analytical formulas derived by Lambert this integral can be reduced to a two-fold integral. This is applied in CalculiX right now: the interacting surfaces are triangulated and the viewfactor between two triangles is calculated by taking a one-point integration for the base triangle (in the center of gravity) and the analytical formula for the integration over the other triangles covering a hemisphere about the base triangle. One can switch to a more accurate integration over the base triangle by increasing the variable “factor” in subroutine radmatrix, look at the comments in that subroutine. This, however, will increase the computational time.

For a heat transfer analysis the conductivity coefficients of the material are needed (using the *CONDUCTIVITY card) and for transient calculations the heat capacity (using the *SPECIFIC HEAT card). Furthermore, for radiation boundary conditions the *PHYSICAL CONSTANTS card is needed, specifying absolute zero in the user's temperature scale and the Boltzmann constant.

Notice that a phase transition can be modeled by a local sharp maximum of the specific heat. The energy $ U$ per unit of mass needed to complete the phase transition satisfies

$\displaystyle U = \int_{T_0}^{T_1} C dT,$ (360)

where $ C$ is the specific heat and $ [T_0,T_1]$ is the temperature interval within which the phase transition takes place.