Product: ABAQUS/Standard
The ABAQUS/Standard capability for uncoupled heat transfer analysis is intended to model solid body heat conduction with general, temperature-dependent conductivity; internal energy (including latent heat effects); and quite general convection and radiation boundary conditions. This section describes the basic energy balance, constitutive models, boundary conditions, finite element discretization, and time integration procedures used.
Heat transfer in flowing materials (convection) is discussed in “Convection/diffusion,” Section 2.11.3. Radiation heat transfer in cavities is discussed in “Cavity radiation,” Section 2.11.4. All such heat transfer mechanisms can be present in a model.
The basic energy balance is (Green and Naghdi)
where V is a volume of solid material, with surface area S;
is the density of the material; 
is the material time rate of the internal energy; q is the heat flux per unit area of the body, flowing into the body; and r is the heat supplied externally into the body per unit volume.It is assumed that the thermal and mechanical problems are uncoupled in the sense that 
only, where 
is the temperature of the material, and q and r do not depend on the strains or displacements of the body. For simplicity a Lagrangian description is assumed, so “volume” and “surface” mean the volume and surface in the reference configuration.
This relationship is usually written in terms of a specific heat, neglecting coupling between mechanical and thermal problems:
Heat conduction is assumed to be governed by the Fourier law,
where
is the conductivity matrix, 
; 
is the heat flux; and 
is position. The conductivity can be fully anisotropic, orthotropic, or isotropic.Boundary conditions can be specified as prescribed temperature, 
; prescribed surface heat flux, 
per area; prescribed volumetric heat flux, 
per volume; surface convection: 
, where 
is the film coefficient and 
is the sink temperature; and radiation: 
, where A is the radiation constant (emissivity times the Stefan-Boltzmann constant) and 
is the absolute zero on the temperature scale used. Surfaces can also participate in cavity radiation effects. The cavity radiation formulation in ABAQUS is described in “Cavity radiation,” Section 2.11.4.
A variational statement of the energy balance, Equation 2.11.1–1, together with the Fourier law, Equation 2.11.1–2, is obtained directly by the standard Galerkin approach as
where
is an arbitrary variational field satisfying the essential boundary conditions. The body is approximated geometrically with finite elements, so the temperature is interpolated as 
are nodal temperatures. The Galerkin approach assumes that , the variational field, is interpolated by the same functions: 
. With these interpolations the variational statement, Equation 2.11.1–3, becomes 
are arbitrarily chosen, this gives the system of equations This set of equations is the “continuous time description” of the geometric approximation.ABAQUS/Standard uses the backward difference algorithm:
This operator is chosen for a number of reasons. First of all, we choose from one-step operators of the form
such operators are only conditionally stable for linear heat transfer problems. We prefer to work with unconditionally stable methods, because ABAQUS is most commonly applied to problems where the solution is sought over very long time periods (compared to the stability limit for the explicit form of the operator, 
), and so choose 
. Of these operators the central difference method, 
, has the highest accuracy. However, that form of the operator tends to produce oscillations in the early time solution that are not present in the backward difference form. Thus, we use 
: backward difference. Introducing the operator, Equation 2.11.1–5, into the energy balance Equation 2.11.1–4 gives This nonlinear system is solved by a modified Newton method. The method is modified Newton because the tangent matrix (the Jacobian matrix)—that is, the rate of change of the left-hand side of Equation 2.11.1–6 with respect to 
—is not formed exactly. The formation of the terms in this tangent matrix is now described.The internal energy term gives a Jacobian contribution:
is the specific heat, 
, outside the latent heat range, and is 
if 
at the integration point, where 
and 
are the liquidus and solidus temperatures and L is the latent heat associated with this phase change.In severe latent heat cases this term can result in numerical instabilities, as the stiffness term 
is small outside the solidus-liquidus temperature range and is very stiff inside that rather narrow range. To avoid such instabilities in those cases this term is modified to a secant term during the early iterations of the solution to a time step. Since the modification occurs only in cases involving latent heat, it affects only those problems.
The conductivity term gives a Jacobian contribution:
With film and radiation conditions, the surface flux term gives a Jacobian contribution:
and, hence, in 
, so a single equation solution provides the . Since the method usually is only a minor modification of Newton's method, convergence is rapid.ABAQUS/Standard uses an automatic (self-adaptive) time stepping algorithm to choose 
. This is based on a user-supplied tolerance on the maximum temperature change allowed in a time increment, and the increment is adjusted according to this parameter, as well as the convergence rate of Equation 2.11.1–7 in nonlinear cases.
The first-order heat transfer elements (such as 2-node link, 4-node quadrilateral, and 8-node brick) use a numerical integration rule with the integration stations located at the corners of the element for the heat capacitance terms. This means that the Jacobian term associated with the internal energy rate is diagonal. This approach is especially effective when strong latent heat effects are present. The second-order elements use conventional Gaussian integration. Thus, second-order elements are to be preferred for problems when the solution will be smooth (without latent heat effects), whereas the first-order elements should be used in nonsmooth cases (with latent heat).
The HEATCAP element is available for modeling lumped heat capacitance at a point. The associated concentrated film and concentrated radiation loading options are specified by the user. These loading options are also allowed in coupled temperature-displacement and coupled thermal-electrical analysis.
