Nonlinear Steady-State Heat Transfer Analysis

(1)

$$\left(K_C+H\right)T+R{\left(T+T_{abs}\right)}^4=f$$

Where,

$K_C$

Conductivity matrix

$H$

Boundary convection matrix, due to free convection.

$T$

Unknown nodal temperature

$R$

Radiation exchange matrix

$T_{abs}$

Absolute temperature scale defined via PARAM, TABS

$f$

Thermal loading vector

Thermal load vector can be expressed as:(2)

$$f=f_B+f_H+f_Q+f_R$$

Where,

$f_B$

Power, due to heat flux at boundary specified by QBDY1.

$f_H$

Boundary convection vector, due to convection specified by CONV (automatic free convection definition can be activated via CONVG Bulk/Subcase pair).

$f_Q$

Power vector, due to internal heat generation specified by QVOL.

$f_R$

Boundary radiation vector, due to radiation specified by RADBC.

The following steps are a guide to setup a Nonlinear Steady-State Heat Transfer Analysis.

1. Use the solution sequence identifier (ANALYSIS) in the Subcase Information Entry section to select the nonlinear steady-state heat transfer analysis using: ANALYSIS=NLHEAT.
2. The likely initial temperature distribution can be defined using the TEMPERATURE Subcase Information Entry (type=INITIAL). A good initial temperature estimate improves the convergence of the solver.
3. The MATT4 Bulk Data Entry can be used to define temperature dependent thermal material properties.
4. To indicate that a nonlinear solution is required for any subcase, a NLPARM Subcase Information Entry is required. This subcase entry points to a NLPARM Bulk Data Entry that specifies convergence tolerances and other nonlinear parameters.
5. Loads and boundary conditions are defined in the Bulk Data Entry section of the input deck. These should be referenced in the Subcase Information Entry section using SPC and LOAD entries in a subcase. Each Subcase defines a load vector.

1. Radiation is considered as a boundary condition (similar to Convection).
2. Radiant exchange is between the surface element and a black body space node.
3. The ambient space node is defined by an SPC Bulk Data Entry and is identified via the NODAMB field of the RADBC Bulk Data Entry.
4. RADBC is only used when applied to a surface element (CHBDYE Bulk Data).
5. The ambient node is treated as a black body with its own temperature for radiation exchange between the surface element and space.
6. The radiation view factor (FAMB on RADBC entry) defines the geometric association between the element surface and the ambient node.
7. The emissivity (EMIS1 field) and absorptivity (ABSORP field) can be defined on the RADM Bulk Data Entry.
8. The RADM entry is directly referenced by a surface element entry (CHBDYE entry).
9. The two parameters are required for radiation calculations:
   • PARAM, TABS: Absolute temperature scale factor
   • PARAM, SIGMA: Stefan-Boltzmann constant

SUBCASE  5
ANALYSIS=NLHEAT
SPC=10
LOAD=20
NLPARM=30
…
BEGIN BULK
…
NLPARM, 30
…
ENDDATA