Coupled thermal-electrical analysis

Abaqus/Standard provides a fully coupled thermal-electrical procedure for analyzing Joule heating, which arises when the energy dissipated by an electrical current flowing through a conductor is converted into thermal energy.

The following topics are discussed:

Governing equations
Constitutive behavior
Thermal energy balance
Thermal energy due to electrical current
Surface conditions
Spatial discretization
Jacobian contributions

Governing equations

The electric field in a conducting material is governed by Maxwell's equation of conservation of charge. Assuming steady-state direct current, the equation reduces to

\begin{matrix} (1) & \int_{S} J \cdot n d S = \int_{V} r_{c} d V, \end{matrix}

where V is any control volume whose surface is S, $n$ is the outward normal to S, $J$ is the electrical current density (current per unit area), and $r_{c}$ is the internal volumetric current source per unit volume.

The divergence theorem is used to convert the surface integral into a volume integral:

\int_{V} [\frac{\partial}{\partial x} \cdot J - r_{c}] d V = 0,

and since the volume is arbitrary, this provides the pointwise differential equation

\frac{\partial}{\partial x} \cdot J - r_{c} = 0.

The equivalent weak form is obtained by introducing an arbitrary, variational, electrical potential field, $δ φ$ , and integrating over the volume:

\int_{V} δ φ [\frac{\partial}{\partial x} \cdot J - r_{c}] d V = 0.

Using first the chain rule and then the divergence theorem, this statement can be rewritten as

- \int_{V} \frac{\partial δ φ}{\partial x} \cdot J d V = \int_{S} δ φ J d S + \int_{V} δ φ r_{c} d V,

where $J \overset{def}{=} - J \cdot n$ is the current density entering the control volume across S.

Constitutive behavior

The flow of electrical current is described by Ohm's law:

J = σ^{E} \cdot E,

where $σ^{E} (θ, f^{α})$ is the electrical conductivity matrix; $θ$ is the temperature; and $f^{α}, α = 1, 2...$ are any predefined field variables. The conductivity can be isotropic, orthotropic, or fully anisotropic. $E (x)$ is the electrical field intensity defined as

\begin{matrix} (2) & E = - \frac{\partial φ}{\partial x} . \end{matrix}

Since a potential rise occurs when a charged particle moves against the electrical field, the direction of the gradient is opposite to that of the electrical field. Using this definition of the electrical field, Ohm's law is rewritten as

\begin{matrix} (3) & J = - σ^{E} \cdot \frac{\partial φ}{\partial x} . \end{matrix}

The constitutive relation is linear; that is, it assumes that the electrical conductivity is independent of the electrical field.

Introducing Ohm's law, the governing conservation of charge equation becomes

\begin{matrix} (4) & \int_{V} \frac{\partial δ φ}{\partial x} \cdot σ^{E} \cdot \frac{\partial φ}{\partial x} d V = \int_{V} δ φ r_{c} d V + \int_{S} δ φ J d S . \end{matrix}

Thermal energy balance

The heat conduction behavior is described by the basic energy balance relation

\begin{matrix} (5) & \int_{V} ρ \dot{U} δ θ d V + \int_{V} \frac{\partial δ θ}{\partial x} \cdot k \cdot \frac{\partial θ}{\partial x} d V = \int_{V} δ θ r d V + \int_{S} δ θ q d S, \end{matrix}

where V is a volume of solid material, with surface area S; $ρ$ is the density of the material; U is the internal energy; $k$ is the thermal conductivity matrix; q is the heat flux per unit area of the body, flowing into the body; and r is the heat generated within the body. The thermal problem is discussed in detail in Uncoupled heat transfer analysis.

Equation 4 and Equation 5 describe the electrical and thermal problems, respectively. Coupling arises from two sources: the conductivity in the electrical problem is temperature dependent, $σ^{E} = σ^{E} (θ)$ , and the internal heat generation in the thermal problem is a function of electrical current, $r = r_{e c} (J)$ , as described below.

Thermal energy due to electrical current

Joule's law describes the rate of electrical energy, $P_{e c}$ , dissipated by current flowing through a conductor as

P_{e c} = E \cdot J .

Using Equation 2 and Equation 3, Joule's law is rewritten as

P_{e c} = E \cdot σ^{E} \cdot E .

In a steady-state analysis $P_{e c}$ is evaluated at time $t + Δ t$ . In a transient analysis an averaged value of $P_{e c}$ is obtained over the increment

\begin{matrix} P_{e c} & = \frac{1}{Δ t} \int_{Δ t} P_{e c} d t \\ = E \cdot σ^{E} \cdot E - E \cdot σ^{E} \cdot Δ E + \frac{1}{3} Δ E \cdot σ^{E} \cdot Δ E, \end{matrix}

where $E$ and $σ^{E}$ are values at time $t + Δ t$ . The amount of this energy released as internal heat is

r = η_{v} P_{e c},

where $η_{v}$ is an energy conversion factor.

Surface conditions

The surface—S—of the body consists of parts on which boundary conditions can be prescribed— $S_{p}$ —and parts that can interact with nearby surfaces of other bodies— $S_{i}$ . Prescribed boundary conditions include the electrical potential, $φ = φ (x, t)$ ; temperature, $θ = θ (x, t)$ ; electrical current density, $J = J (x, t)$ ; heat flux, $q = q (x, t)$ ; and surface convection and radiation conditions. The surface interaction model includes heat conduction and radiation effects between the interface surfaces and electrical current flowing across the interface. Heat conduction and radiation are modeled by

q_{c} = k_{g} (θ_{B} - θ),

and

q_{r} = F_{B} {(θ_{B} - θ^{z})}^{4} - F {(θ - θ^{z})}^{4},

respectively, where $θ$ is the temperature on the surface of the body under consideration, $θ_{B}$ is the temperature on the surface of the other body, $θ^{z}$ is the value of absolute zero on the temperature scale being used, $k_{g} (\bar{θ}, {\bar{f}}^{α})$ is the gap thermal conductance, $\bar{θ} = \frac{1}{2} (θ + θ_{B})$ is the average interface temperature, ${\bar{f}}^{α} = \frac{1}{2} (f_{A}^{α} + f_{B}^{α})$ is the average of any predefined field variables at A and B, and F and $F_{B}$ are constants.

The electrical current flowing between the interface surfaces is modeled as

J = σ_{g} (φ_{B} - φ),

where $φ$ is the electrical potential on the surface of the body under consideration, $φ_{B}$ is the electrical potential on the surface of the other body, and $σ_{g} (\bar{θ}, {\bar{f}}^{α})$ is the gap electrical conductance. The electrical energy dissipated by the current flowing across the interface,

P_{e c} = J (φ_{B} - φ) = σ_{g} {(φ_{B} - φ)}^{2},

is released as heat on the surfaces of the bodies:

q_{e c} = f η_{g} P_{e c}, and q_{e c}^{B} = (1 - f) η_{g} P_{e c},

where $η_{g}$ is an energy conversion factor and f specifies how the total heat is distributed between the interface surfaces. $P_{e c}$ is evaluated at the end of the time increment in a steady-state analysis, and an averaged value over the time increment is used in a transient analysis. This is described in detail in Heat generation caused by electrical current.

Introducing the surface interaction effects and electrical energy released as thermal energy, the governing electric and thermal equations become

\begin{matrix} (6) & \int_{V} \frac{\partial δ φ}{\partial x} \cdot σ^{E} \cdot \frac{\partial φ}{\partial x} d V = \int_{V} δ φ r_{c} d V + \int_{S_{p}} δ φ J d S + \int_{S_{i}} δ φ σ_{g} (φ_{B} - φ) d S, \end{matrix}

and

\begin{matrix} (7) & \begin{matrix} \int_{V} ρ \dot{U} δ θ d V & + \int_{V} \frac{\partial δ θ}{\partial x} \cdot k \cdot \frac{\partial θ}{\partial x} d V = \int_{V} δ θ r d V + \int_{V} δ θ η_{v} P_{e c} d V \\ + \int_{S_{p}} δ θ q d S + \int_{S_{i}} δ θ (q_{c} + q_{r} + q_{e c}) d S . \end{matrix} \end{matrix}

Spatial discretization

In a finite element model equilibrium is approximated as a finite set of equations by introducing interpolation functions. Discretized quantities are indicated by uppercase superscripts (for example, $φ^{N}$ ). The summation convention is adopted for the superscripts. The discretized quantities represent nodal variables, with nodes shared between adjacent elements and appropriate interpolation chosen to provide adequate continuity of the assumed variation.

The virtual electrical potential field is interpolated by

δ φ = N^{N} δ φ^{N},

where $N^{N}$ are the interpolation functions. The discretized electrical equation is then written as

δ φ^{N} {\int_{V} \frac{\partial N^{N}}{\partial x} \cdot σ^{E} \cdot \frac{\partial φ}{\partial x} d V = \int_{V} N^{N} r_{c} d V + \int_{S_{p}} N^{N} J d S + \int_{S_{i}} N^{N} σ_{g} (φ_{B} - φ) d S} .

Since $δ φ$ is arbitrary,

\begin{matrix} (8) & I_{φ}^{N} = \int_{V} \frac{\partial N^{N}}{\partial x} \cdot σ^{E} \cdot \frac{\partial φ}{\partial x} d V - \int_{V} N^{N} r_{c} d V - \int_{S_{p}} N^{N} J d S - \int_{S_{i}} N^{N} σ_{g} (φ_{B} - φ) d S = 0. \end{matrix}

The temperature field in the thermal problem is approximated by the same set of interpolation functions:

δ θ = N^{P} δ θ^{P} .

Using these interpolation functions and a backward difference operator to integrate the internal energy rate, $\dot{U}$ , the thermal energy balance relation is obtained:

\begin{matrix} (9) & \begin{matrix} I_{θ}^{P} & = \frac{1}{Δ t} \int_{V} N^{P} ρ (U_{t + Δ t} - U_{t}) d V + \int_{V} \frac{\partial N^{P}}{\partial x} \cdot k \cdot \frac{\partial θ}{\partial x} d V - \int_{V} N^{P} r d V \\ - \int_{V} N^{P} η_{v} P_{e c} d V - \int_{S_{p}} N^{P} q d S - \int_{S_{i}} N^{P} (q_{c} + q_{r} + q_{e c}) d S = 0. \end{matrix} \end{matrix}

Jacobian contributions

The Jacobian contributions are obtained by taking variations of Equation 8 and Equation 9 with respect to the electrical potential, $φ$ , and the temperature, $θ$ , at time $t + Δ t$ . This yields

\begin{matrix} K_{φ φ}^{N M} = \frac{\partial I_{φ}^{N}}{\partial φ^{M}} & = \int_{V} \frac{\partial N^{N}}{\partial x} \cdot σ^{E} \cdot \frac{\partial N^{M}}{\partial x} d V + \int_{S_{i}} N^{N} σ_{g} N^{M} d S, \\ K_{φ θ}^{N Q} = \frac{\partial I_{φ}^{N}}{\partial θ^{Q}} & = - \int_{V} \frac{\partial N^{N}}{\partial x} \cdot \frac{\partial σ^{E}}{\partial θ} \cdot E N^{Q} d V - \frac{1}{2} \int_{S_{i}} N^{N} \frac{\partial σ_{g}}{\partial \bar{θ}} (φ_{B} - φ) N^{Q} d S, \end{matrix}

\begin{matrix} K_{θ φ}^{P M} = \frac{\partial I_{θ}^{P}}{\partial φ^{M}} & = \int_{V} N^{P} η_{v} (J - \frac{1}{3} Δ J) \cdot \frac{\partial N^{M}}{\partial x} d V + 2 \int_{S_{i}} N^{P} f η_{g} σ_{g} (φ_{B} - φ) N^{M} d S, \\ K_{θ θ}^{P Q} = \frac{\partial I_{θ}^{P}}{\partial θ^{Q}} & = \frac{1}{Δ t} \int_{V} N^{P} ρ \frac{d U}{d θ} N^{Q} d V + \int_{V} \frac{\partial N^{P}}{\partial x} \cdot k \cdot \frac{\partial N^{Q}}{\partial x} d V + \int_{V} \frac{\partial N^{P}}{\partial x} \cdot \frac{\partial k}{\partial θ} \cdot \frac{\partial θ}{\partial x} N^{Q} d V \\ - \int_{V} N^{P} η_{v} (E \cdot \frac{\partial σ^{E}}{\partial θ} \cdot E - E \cdot \frac{\partial σ^{E}}{\partial θ} \cdot Δ E + \frac{1}{3} Δ E \cdot \frac{\partial σ^{E}}{\partial θ} \cdot Δ E) N^{Q} d V \\ - \int_{S_{p}} N^{P} \frac{\partial q}{+ \partial θ} N^{Q} d S - \int_{S_{i}} N^{P} (\frac{\partial q_{c}}{\partial θ} + \frac{\partial q_{r}}{\partial θ} + \frac{\partial q_{e c}}{\partial θ}) N^{Q} d S . \end{matrix}

The term $\partial q / \partial θ$ in the $K_{θ θ}^{P Q}$ component includes prescribed surface convection and radiation conditions. The surface interaction terms $\partial q_{c} / \partial θ$ , $\partial q_{r} / \partial θ$ , and $\partial q_{e c} / \partial θ$ are evaluated in Heat generation caused by electrical current.

The Jacobian contributions give rise to an unsymmetric system of equations, requiring the use of the nonsymmetric matrix storage and solution scheme.