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 asABAQUS/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 formThe internal energy term gives a Jacobian contribution:
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:
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.