Lumped Capacitance

In the study of transient heat conduction, determining how temperature evolves within a solid body over time and space i

In the study of transient heat conduction, determining how temperature evolves within a solid body over time and space is a formidable task. Mathematically, this typically requires solving complex partial differential equations (PDEs) where temperature $T$ is a function of spatial coordinates and time, expressed as $T = f(x, y, z, t)$.

However, in many engineering contexts, we encounter scenarios where the internal thermal conductivity of a material is so high that the temperature remains nearly uniform throughout the body during a heating or cooling process. When the internal resistance to heat flow is negligible compared to the external convective resistance, we can simplify the problem significantly. This simplification is known as the Lumped Capacitance Method.
The fundamental premise of the lumped capacitance method is the assumption of spatial uniformity. Instead of tracking temperature at every specific point within a volume, we treat the entire object as a single "lump" with a uniform temperature $T(t)$ that changes only with time.

This approach relies on several key assumptions:

  • Negligible Temperature Gradients: The temperature difference between the center and the surface of the object is assumed to be zero ($\nabla T \approx 0$).
  • Dominant Convection/Conduction Ratio: The rate of heat transfer via internal conduction is much faster than the rate of heat transfer via surface convection.
  • Dimensional Reduction: The problem is reduced from a complex partial differential equation (PDE) to a much more manageable first-order ordinary differential equation (ODE).

Mathematical Derivation

The derivation of the lumped capacitance model is rooted in the Principle of Conservation of Energy. For a solid body undergoing a temperature change, the rate of change of its internal energy must equal the rate of heat transfer at its surface.

Let us define the following parameters:

  • $\rho$: Density of the body ($\text{kg/m}^3$)
  • $V$: Volume of the body ($\text{m}^3$)
  • $C_p$: Specific heat capacity ($\text{J/kg}\cdot\text{K}$)
  • $A$: Surface area of the body ($\text{m}^2$)
  • $h$: Convective heat transfer coefficient ($\text{W/m}^2\cdot\text{K}$)
  • $T(t)$: Instantaneous temperature of the body at time $t$
  • $T_\infty$: Ambient temperature of the surrounding fluid

The energy balance equation is expressed as:
$$\rho V C_p \frac{dT(t)}{dt} = -h A [T(t) - T_\infty]$$

This is a first-order linear differential equation. By employing the method of separation of variables and integrating from the initial state ($t=0, T=T_i$) to a general state, we arrive at the analytical solution:

$$\frac{T(t) - T_\infty}{T_i - T_\infty} = \exp\left( -\frac{h A}{\rho V C_p} t \right)$$

To simplify this expression, we often introduce a time constant, $\tau$:
$$\tau = \frac{\rho V C_p}{h A}$$

The equation then takes a classic exponential decay form:
$$\theta(t) = \theta_i \exp\left( -\frac{t}{\tau} \right)$$
where $\theta$ represents the dimensionless temperature difference.

Determining Validity: The Biot Number

The most critical question in applying this method is: Is the assumption of uniform temperature actually valid? To answer this, we use a dimensionless parameter called the Biot Number ($Bi$).

1. Definition and Physical Interpretation

The Biot number is defined as the ratio of the internal conductive resistance to the external convective resistance:
$$Bi = \frac{h L_c}{k}$$

Where:

  • $h$ is the convective heat transfer coefficient.
  • $k$ is the thermal conductivity of the solid.
  • $L_c$ is the characteristic length.

Physically, $Bi$ tells us which process dominates the thermal response:

  • Small $Bi$ ($Bi \ll 1$): Conduction within the solid is extremely efficient. The heat spreads through the body much faster than the surface can exchange it with the environment, resulting in a nearly uniform temperature.
  • Large $Bi$ ($Bi \gg 1$): Convection at the surface is very rapid, but the internal conduction is slow. This creates significant temperature gradients, where the surface temperature changes much faster than the core temperature.

2. The Characteristic Length ($L_c$)

The characteristic length $L_c$ is a geometric parameter defined as the ratio of the volume to the surface area:
$$L_c = \frac{V}{A}$$

The value of $L_c$ varies depending on the geometry of the object:

  • Infinite Flat Plate (thickness $2L$): $L_c = L$
  • Long Cylinder (radius $r$): $L_c = r/2$
  • Sphere (radius $r$): $L_c = r/3$

3. Engineering Criterion

In practical engineering applications, the lumped capacitance method is considered a reliable approximation if:
$$Bi < 0.1$$

If $Bi$ exceeds this threshold, the temperature gradients within the body become too significant to ignore, and one must resort to more sophisticated methods, such as Heisler charts, infinite series solutions, or numerical techniques like Finite Element Analysis (FEA).

Worked Example: Cooling of a Copper Sphere

Problem Statement:
A copper sphere with a radius of $r = 0.02,\text{m}$ is initially at $100^\circ\text{C}$. It is submerged in a fluid at $20^\circ\text{C}$ with a convective heat transfer coefficient of $h = 500,\text{W/m}^2\cdot\text{K}$.
Given the properties of copper:

  • Thermal conductivity $k = 400,\text{W/m}\cdot\text{K}$
  • Density $\rho = 8900,\text{kg/m}^3$
  • Specific heat $C_p = 385,\text{J/kg}\cdot\text{K}$

Calculate:

  1. Whether the lumped capacitance method is applicable.
  2. The temperature of the sphere after 60 seconds.

Solution:

Step 1: Calculate the characteristic length ($L_c$)
For a sphere:
$$L_c = \frac{r}{3} = \frac{0.02}{3} \approx 0.00667,\text{m}$$

Step 2: Calculate the Biot Number ($Bi$)
$$Bi = \frac{h L_c}{k} = \frac{500 \times 0.00667}{400} \approx 0.0083$$
Since $0.0083 < 0.1$, the lumped capacitance method is valid.

Step 3: Determine the time constant ($\tau$)
$$\tau = \frac{\rho L_c C_p}{h} = \frac{8900 \times 0.00667 \times 385}{500} \approx 45.67,\text{s}$$

Step 4: Calculate the temperature at $t = 60,\text{s}$
Using the exponential decay formula:
$$\frac{T(60) - 20}{100 - 20} = \exp\left( -\frac{60}{45.67} \right)$$
$$\frac{T(60) - 20}{80} = \exp(-1.3137) \approx 0.2688$$
$$T(60) = 20 + (80 \times 0.2688) \approx 41.5^\circ\text{C}$$

Result: After 60 seconds, the temperature of the copper sphere is approximately $41.5^\circ\text{C}$.

Summary

The lumped capacitance method is an indispensable tool for rapid thermal analysis. By utilizing the Biot Number to assess the competition between internal conduction and external convection, engineers can bypass complex spatial calculations when the thermal gradients are minimal. While it is a simplification, it provides highly accurate results for high-conductivity materials and small geometries, serving as a vital first step in any transient thermal investigation.