Skip to content
Ahmet Çelik
← Back to Study Guides

Chapter 5: Transient Conduction

MECH302

Source files used: Heat Transfer (13), Heat Transfer (14), related review comments in later combined-mode lectures; textbook Chapter 5 as support. No worked examples are included.

1. Big Picture

Chapter 5 adds time to conduction.

In Chapters 2 and 3, many problems were steady: temperature did not change with time. In Chapter 5, temperature changes as an object heats up or cools down.

The big decision in Biddle’s lecture is the Biot number:

Bi=hLck\boxed{Bi=\frac{hL_c}{k}}

The Biot number compares internal conduction resistance to external convection resistance.

  • Small BiBi: temperature inside the body is nearly uniform.
  • Large BiBi: temperature varies significantly inside the body.

So the first question in a transient conduction problem is not “which equation do I use?” It is:

Can I treat the object temperature as spatially uniform?\boxed{\text{Can I treat the object temperature as spatially uniform?}}

If yes, use the lumped heat capacity model. If no, use spatial transient conduction methods such as one-term approximations or Heisler/Grober charts.


2. Core Ideas

2.1 Lumped Heat Capacity Model

The lumped model assumes the whole object has one temperature at a given time:

T=T(t)T=T(t)

Temperature changes with time but not with position inside the solid.

This is valid when internal conduction is fast compared with surface convection. In that case, the inside of the object does not have enough thermal resistance to develop a major temperature gradient.

Biddle’s decision rule:

Bi<0.1lumped model is acceptableBi<0.1 \quad \Rightarrow \quad \text{lumped model is acceptable}

2.2 Spatial Effects

If Bi>0.1Bi>0.1, the object cannot be treated as isothermal. Then temperature depends on position and time:

T=T(x,t)orT=T(r,t)T=T(x,t) \quad \text{or} \quad T=T(r,t)

For a plane wall, centerline and surface temperatures are different. For a cylinder, center and surface temperatures are different. The center usually responds more slowly than the surface.


2.3 Characteristic Length

For lumped capacitance:

Lc=VAsL_c=\frac{V}{A_s}

where:

  • VV: object volume,
  • AsA_s: surface area exposed to convection.

For one-term transient charts/equations, the characteristic length depends on geometry:

GeometryOne-Term BiBi and FoFo lengthV/AsV/A_s (lumped Bi check only)
Plane wall of total thickness 2L2LLL, half-thicknessLL
Long cylinderror_oro/2r_o/2
Sphereror_oro/3r_o/3

⚠️ Common error: ro/2r_o/2 (cylinder) and ro/3r_o/3 (sphere) are the V/AsV/A_s values used only when checking lumped capacitance (Bi=hLc/k<0.1Bi = hL_c/k < 0.1). Once you move to the one-term approximation, use ror_o for both BiBi and FoFo. Confirmed in textbook Ch. 5 and Solutions Manual (Problem 5.72 comments).

Do not mix the lumped V/AsV/A_s length with the Heisler/one-term geometry length.


2.4 Fourier Number

The Fourier number measures dimensionless time:

Fo=αtL2Fo=\frac{\alpha t}{L^2}

or for cylinders/spheres:

Fo=αtro2Fo=\frac{\alpha t}{r_o^2}

Large FoFo means more time has passed for heat to diffuse through the object.


2.5 Energy Transferred

In transient conduction, you may be asked not only for temperature but also for how much energy has been gained or lost.

For lumped systems:

Q=ρVcp(TiT(t))Q=\rho V c_p\left(T_i-T(t)\right)

for cooling, with sign chosen according to the problem.

Biddle explicitly corrected a lecture note mistake: for energy stored in a cylinder, use volume πr2L\pi r^2L, not surface area.


3. Main Governing Equations and Formulas

3.1 Biot Number

Bi=hLckBi=\frac{hL_c}{k}

where:

  • hh: convection coefficient,
  • LcL_c: characteristic length,
  • kk: solid thermal conductivity.

Use to decide whether lumped capacitance is valid.

Decision:

Bi<0.1lumped heat capacity modelBi<0.1 \Rightarrow \text{lumped heat capacity model} Bi>0.1spatial transient methodBi>0.1 \Rightarrow \text{spatial transient method}

3.2 Characteristic Length for Lumped Model

Lc=VAsL_c=\frac{V}{A_s}

Use only for the lumped-capacitance Biot number.


3.3 Lumped Capacitance Temperature Response

T(t)TTiT=exp(hAsρVcpt)\frac{T(t)-T_\infty}{T_i-T_\infty}=\exp\left(-\frac{hA_s}{\rho Vc_p}t\right)

or

θθi=et/τ\frac{\theta}{\theta_i}=e^{-t/\tau}

where:

θ=TT\theta=T-T_\infty θi=TiT\theta_i=T_i-T_\infty τ=ρVcphAs\tau=\frac{\rho Vc_p}{hA_s}

τ\tau is the thermal time constant.

Use when Bi<0.1Bi<0.1.


3.4 Lumped Heat Rate

Instantaneous convection heat rate:

q˙(t)=hAs[T(t)T]\dot{q}(t)=hA_s\left[T(t)-T_\infty\right]

Heat rate decreases as the object temperature approaches the fluid temperature.


3.5 Lumped Energy Transferred

Energy transferred to/from the surroundings by convection over time:

Q=0thAsθdtQ=\int_{0}^{t} hA_s\theta\,dt

Using the lumped result θ=θiet/τ\theta=\theta_i e^{-t/\tau}:

Q=(ρcV)θi[1exp(tτ)]Q=(\rho cV)\theta_i\left[1-\exp\left(-\frac{t}{\tau}\right)\right]

For cooling:

Q=ρVcp(TiT(t))Q=\rho Vc_p\left(T_i-T(t)\right)

Maximum possible energy transfer:

Q0=ρVcp(TiT)Q_0=\rho Vc_p\left(T_i-T_\infty\right)

Fraction transferred:

QQ0=1T(t)TTiT\frac{Q}{Q_0}=1-\frac{T(t)-T_\infty}{T_i-T_\infty}

3.6 Fourier Number

Plane wall:

Fo=αtL2Fo=\frac{\alpha t}{L^2}

Long cylinder or sphere:

Fo=αtro2Fo=\frac{\alpha t}{r_o^2}

Use with one-term approximations or transient charts.


3.7 One-Term Approximation: Plane Wall

For a plane wall of total thickness 2L2L, with convection on both sides:

θθi=C1eζ12Focos(ζ1xL)\frac{\theta}{\theta_i}=C_1e^{-\zeta_1^2Fo}\cos\left(\zeta_1\frac{x}{L}\right)

where:

θ=TT\theta=T-T_\infty Fo=αtL2Fo=\frac{\alpha t}{L^2} x=xLx^*=\frac{x}{L}

C1C_1 and ζ1\zeta_1 come from the transient conduction table as functions of Bi=hL/kBi=hL/k.

At the centerline, x=0x=0, so cos(0)=1\cos(0)=1:

θ0θi=C1eζ12Fo\frac{\theta_0}{\theta_i}=C_1e^{-\zeta_1^2Fo}

Use when Bi>0.1Bi>0.1 and the one-term approximation is valid, commonly for sufficiently large FoFo such as Fo>0.2Fo>0.2.


3.8 One-Term Approximation: Long Cylinder

For a long cylinder:

θθi=C1eζ12FoJ0(ζ1rro)\frac{\theta}{\theta_i}=C_1e^{-\zeta_1^2Fo}J_0\left(\zeta_1\frac{r}{r_o}\right)

where:

Fo=αtro2Fo=\frac{\alpha t}{r_o^2}

and J0J_0 is the Bessel function of the first kind.

At the centerline, r=0r=0, J0(0)=1J_0(0)=1:

θ0θi=C1eζ12Fo\frac{\theta_0}{\theta_i}=C_1e^{-\zeta_1^2Fo}

Use table values for C1C_1 and ζ1\zeta_1.


3.9 One-Term Approximation: Sphere

For a sphere:

θθi=C1eζ12Fosin(ζ1r/ro)ζ1r/ro\frac{\theta}{\theta_i}=C_1e^{-\zeta_1^2Fo}\frac{\sin(\zeta_1r/r_o)}{\zeta_1r/r_o}

At the center:

θ0θi=C1eζ12Fo\frac{\theta_0}{\theta_i}=C_1e^{-\zeta_1^2Fo}

Use sphere table values for C1C_1 and ζ1\zeta_1.


3.10 Multidimensional Product Solution

For shapes that can be built from intersections of one-dimensional shapes, approximate:

(θθi)multiD=(θθi)1D\left(\frac{\theta}{\theta_i}\right)_{multi-D} =\prod \left(\frac{\theta}{\theta_i}\right)_{1D}

Examples:

  • finite cylinder = infinite cylinder solution × plane wall solution,
  • rectangular solid = plane wall solution in xx × plane wall solution in yy × plane wall solution in zz.

Use only when the textbook charts/tables support the geometry.


4. Problem-Solving Workflow

4.1 Lumped-Capacitance Workflow

  1. Identify the body, volume VV, and exposed surface area AsA_s.
  2. Compute Lc=V/AsL_c=V/A_s.
  3. Compute Bi=hLc/kBi=hL_c/k.
  4. If Bi<0.1Bi<0.1, use lumped capacitance.
  5. Compute τ=ρVcp/(hAs)\tau=\rho Vc_p/(hA_s).
  6. Use θ/θi=et/τ\theta/\theta_i=e^{-t/\tau}.
  7. Solve for unknown time, temperature, or heat transfer.

4.2 Spatial Transient Workflow

  1. Compute BiBi using the geometry-specific length.
  2. Since Bi>0.1Bi>0.1, do not use lumped capacitance.
  3. Identify geometry: plane wall, cylinder, sphere.
  4. Compute Fo=αt/L2Fo=\alpha t/L^2 or αt/ro2\alpha t/r_o^2.
  5. Get C1C_1 and ζ1\zeta_1 from the proper table for the geometry and Bi.
  6. Use the one-term expression to find center temperature.
  7. Use position factor to find off-center or surface temperature if needed.
  8. Use energy charts/equations if asked for total energy transferred.

5. Decision Rules / Decision Trees

5.1 Main Transient Decision

Compute Bi=hLc/kBi=hL_c/k.

Bi<0.1Bi<0.1? → lumped heat capacity: T=T(t)T=T(t) only.

Bi>0.1Bi>0.1? → spatial effects matter: T=T(x,t)T=T(x,t) or T=T(r,t)T=T(r,t).

5.2 Which Length?

Checking lumped capacitance? → Lc=V/AsL_c=V/A_s

Using plane-wall transient chart/equation? → L = half-thickness

Using long-cylinder chart/equation? → L=roL=r_o

Using sphere chart/equation? → L=roL=r_o

5.3 What Is Being Asked?

Temperature after time t? → use temperature ratio equation/chart

Time to reach temperature? → solve temperature ratio for tt or FoFo

Energy transferred? → use Q/Q0Q/Q_0 relation or lumped energy balance

Surface temperature? → find center temperature first, then use surface/center relation

5.4 Lumped Heating vs Cooling

Object hotter than fluid? → cooling, TT decreases toward TT_\infty

Object colder than fluid? → heating, TT increases toward TT_\infty

In both cases: θ/θi=exp(t/τ)\theta/\theta_i=\exp(-t/\tau)


6. Important Tables / Correlations Needed

6.1 Biot Number Rule

Biot NumberModel
Bi<0.1Bi<0.1Lumped capacitance acceptable
Bi>0.1Bi>0.1Spatial gradients important

6.2 Geometry Lengths

GeometryLumped LcL_cOne-Term/Chart Length
Plane wall, total thickness 2L2LV/AsV/A_shalf-thickness LL
Long cylinderV/AsV/A_sradius ror_o
SphereV/AsV/A_sradius ror_o

6.3 Table 5.1 Type Data

For one-term approximations, use tables that give:

  • ζ1\zeta_1,
  • C1C_1,
  • as functions of BiBi,
  • separately for plane wall, cylinder, and sphere.

Do not use plane wall constants for a cylinder or sphere.


7. Key Takeaways

  • Chapter 5 is about temperature changing with time.
  • Always start with the Biot number.
  • Bi<0.1Bi<0.1 means the lumped model is valid.
  • Lumped capacitance assumes uniform temperature inside the object.
  • The lumped response is exponential.
  • τ=ρVcp/(hAs)\tau=\rho Vc_p/(hA_s) controls how fast the object responds.
  • Bi>0.1Bi>0.1 means spatial gradients matter.
  • For spatial effects, use plane wall/cylinder/sphere one-term solutions or charts.
  • Use the correct characteristic length for the method.
  • For energy stored in a cylinder, use cylinder volume πr2L\pi r^2L, not surface area.