Nonlinear Steady-State Heat Transfer Analysis

Calculates the temperature distribution in a system, in which material properties are a function of temperature.

(1)
( K C +H )T+R ( T+ T abs ) 4 =f MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaadaqadaqaaiaahUeadaWgaaWcbaGaam 4qaaqabaGccqGHRaWkcaWHibaacaGLOaGaayzkaaGaaCivaiabgUca RiaahkfadaqadaqaaiaahsfacaWHRaGaaCivamaaBaaaleaacaWHHb GaaCOyaiaahohaaeqaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacaaI 0aaaaOGaeyypa0JaaCOzaaaa@434B@
Where,
K C MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHlbWaaSbaaSqaaiaadoeaaeqaaa aa@3383@
Conductivity matrix
H
Boundary convection matrix, due to free convection.
T
Unknown nodal temperature
R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWGsbaaaa@3292@
Radiation exchange matrix
T abs MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHubWaaSbaaSqaaiaahggacaWHIb GaaC4Caaqabaaaaa@3595@
Absolute temperature scale defined via PARAM, TABS
f
Thermal loading vector

The system of equations is solved to find nodal temperature T .

Thermal load vector can be expressed as:(2)
f= f B + f H + f Q + f R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbGaeyypa0JaaCOzamaaBaaale aacaWGcbaabeaakiabgUcaRiaahAgadaWgaaWcbaGaamisaaqabaGc cqGHRaWkcaWHMbWaaSbaaSqaaiaadgfaaeqaaOGaey4kaSIaaCOzam aaBaaaleaacaWGsbaabeaaaaa@3E21@
Where,
f B MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
Power, due to heat flux at boundary specified by QBDY1.
f H MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
Boundary convection vector, due to convection specified by CONV (automatic free convection definition can be activated via CONVG Bulk/Subcase pair).
f Q MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
Power vector, due to internal heat generation specified by QVOL.
f R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
Boundary radiation vector, due to radiation specified by RADBC.

The matrix on the left-hand side of Equation 1 is singular, unless temperature boundary conditions are specified.

The system of equation (shown above) is solved using the Newton's method. Solution control is provided by defining parameters on the NLPARM Bulk Data Entry. TEMPERATURE(INITIAL) can be used to provide an initial temperature distribution. The temperature results from the nonlinear heat transfer analysis can be used in subsequent structural analysis.

Nonlinear Steady-State Heat Transfer Analysis Setup

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.

Input for Radiation

  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

Example: ANALYSIS and NLPARM Usage

SUBCASE  5
ANALYSIS=NLHEAT
SPC=10
LOAD=20
NLPARM=30
…
BEGIN BULK
…
NLPARM, 30
…
ENDDATA
Note:
  1. Optimization based on nonlinear heat transfer analysis is currently not supported.
  2. CNTNLSUB is currently not supported for nonlinear heat transfer analysis.
  3. Shell elements are membranes in Heat Transfer Analysis. Composite properties (PCOMP) are supported in NLHEAT by homogenization. The temperature distribution through the thickness of shell elements is not calculated. Only nodal temperature is determined (1 degree of freedom per grid).
  4. Flux for composite plate or shell elements are calculated with homogenized conductivity of the entire element.