Mathematical Description of Two-Dimensional Steady-State Heat Conduction Problem

While one-dimensional steady-state heat conduction models offer an intuitive grasp of thermal physics, they are often in

While one-dimensional steady-state heat conduction models offer an intuitive grasp of thermal physics, they are often insufficient for addressing the complexities of modern engineering. In practical applications—such as the thermal management of integrated circuits, the design of high-performance heat sinks, or the analysis of building envelopes—temperature variations occur across multiple spatial dimensions. To accurately predict temperature distributions in these scenarios, we must transition to a two-dimensional (2D) mathematical framework.

This article provides a rigorous mathematical description of the 2D steady-state heat conduction problem, deriving the governing equations from first principles and categorizing the essential boundary conditions required for a complete solution.
The mathematical modeling of heat transfer begins with Fourier’s Law of Heat Conduction. This empirical law states that the local heat flux density is proportional to the negative gradient of the temperature. In a two-dimensional Cartesian coordinate system $(x, y)$, the heat flux vector $\mathbf{q}$ is expressed as:

$$\mathbf{q} = -k \nabla T$$

Where:

  • $k$ represents the thermal conductivity of the material ($\text{W/m}\cdot\text{K}$).
  • $\nabla T$ is the temperature gradient.

In a 2D plane, the components of the heat flux in the $x$ and $y$ directions are defined as:

  1. $q_x = -k \frac{\partial T}{\partial x}$
  2. $q_y = -k \frac{\partial T}{\partial y}$

Physically, this implies that heat energy naturally flows from regions of higher temperature to regions of lower temperature, with the magnitude of the flow determined by both the material's ability to conduct heat and the steepness of the temperature change.

Derivation of the Governing Equation

To describe the temperature field $T(x, y)$ across a continuous domain, we apply the Principle of Conservation of Energy to an infinitesimal control volume.

1. Energy Balance on a Control Volume

Consider a microscopic rectangular element located at $(x, y)$ with dimensions $dx$ and $dy$. Under steady-state conditions, the temperature at any given point does not change over time. Therefore, the energy balance within this element simplifies to:

$$\text{Net Heat Flux Inflow} + \text{Internal Heat Generation} = 0$$

2. Mathematical Formulation

If we assume the presence of an internal heat source $\dot{q}$ (heat generated per unit volume), we analyze the energy exchange as follows:

  • X-direction balance: The net heat entering via the $x$-direction is $-\frac{\partial q_x}{\partial x} dx dy$.
  • Y-direction balance: The net heat entering via the $y$-direction is $-\frac{\partial q_y}{\partial y} dy dx$.
  • Generation term: The total heat generated within the element is $\dot{q} dx dy$.

Summing these components and setting them to zero yields:

$$-\frac{\partial q_x}{\partial x} dx dy - \frac{\partial q_y}{\partial y} dx dy + \dot{q} dx dy = 0$$

By dividing through by the area $dx dy$, we obtain the continuity equation for heat flux:

$$\frac{\partial q_x}{\partial x} + \frac{\partial q_y}{\partial y} = \dot{q}$$

3. The Partial Differential Equation (PDE)

By substituting Fourier’s Law into the continuity equation, we derive the governing PDE for temperature:

$$\frac{\partial}{\partial x} \left( k \frac{\partial T}{\partial x} \right) + \frac{\partial}{\partial y} \left( k \frac{\partial T}{\partial y} \right) + \dot{q} = 0$$

For an isotropic material where the thermal conductivity $k$ is constant, the equation simplifies to:

$$\nabla^2 T = -\frac{\dot{q}}{k}$$

Depending on the presence of internal heat generation, this equation takes two fundamental forms:

  • Laplace Equation ($\nabla^2 T = 0$): Used when there is no internal heat generation ($\dot{q} = 0$).
  • Poisson Equation ($\nabla^2 T = -\frac{\dot{q}}{k}$): Used when internal heat generation is present ($\dot{q} \neq 0$).

Boundary Conditions

A partial differential equation defines the physics within a domain, but it does not yield a unique solution without specified constraints at the boundaries. In 2D heat conduction, three types of boundary conditions are typically employed:

1. Dirichlet Condition (First Kind)

This condition specifies a fixed temperature along the boundary. For instance, if one edge of a plate is held in contact with a constant-temperature reservoir:
$$T(x, y) = T_{surface}$$

2. Neumann Condition (Second Kind)

This condition specifies a prescribed heat flux across the boundary. Mathematically, this involves the temperature gradient. A common special case is the adiabatic (insulated) boundary, where no heat flows through the surface:
$$-k \frac{\partial T}{\partial n} = 0 \implies \frac{\partial T}{\partial n} = 0$$
(where $n$ is the direction normal to the boundary).

3. Robin Condition (Third Kind)

The Robin condition is a combination of the temperature and its gradient, used to model convective heat transfer between the surface and a surrounding fluid. It is expressed as:
$$-k \frac{\partial T}{\partial n} = h(T - T_\infty)$$
Where $h$ is the convective heat transfer coefficient and $T_\infty$ is the ambient fluid temperature. This is the most frequent condition used in engineering to simulate cooling via air or liquid.

Practical Application Example

To illustrate these concepts, consider a rectangular metal plate of length $L$ and width $W$ with a uniform internal heat source $\dot{q}$. The thermal environment is defined as follows:

  • The bottom edge ($y=0$) is maintained at a constant temperature $T_{bottom}$.
  • The top edge ($y=W$) is exposed to ambient air via convection ($h, T_\infty$).
  • The left and right edges ($x=0, x=L$) are perfectly insulated.

The resulting mathematical model is:

  1. Governing Equation: $\frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\dot{q}}{k} = 0$
  2. Boundary Conditions:
    • $T(x, 0) = T_{bottom}$ (Dirichlet)
    • $-k \frac{\partial T}{\partial y} \big|{y=W} = h(T(x, W) - T\infty)$ (Robin)
    • $\frac{\partial T}{\partial x} \big|_{x=0} = 0$ (Neumann/Adiabatic)
    • $\frac{\partial T}{\partial x} \big|_{x=L} = 0$ (Neumann/Adiabatic)

Solving this system allows engineers to map the entire temperature field, identifying "hot spots" that could lead to structural failure or component degradation.

Conclusion

The mathematical description of 2D steady-state heat conduction provides the essential bridge between physical phenomena and computational analysis. By synthesizing Fourier's Law with the principle of energy conservation, we arrive at the Laplace and Poisson equations. When coupled with appropriate Dirichlet, Neumann, or Robin boundary conditions, these equations form the bedrock of thermal simulation, enabling the development of sophisticated numerical tools like Finite Element Analysis (FEA) used in modern industry.