SJMaths.com • UPSC ASO / Day 19: Heat Transfer II
DGCA Airworthiness & Engineering Stream UPSC Tier-2 Core Syllabus

Day 19: Fins & Transient Conduction

Comprehensive study module engineered for the UPSC Air Safety Officer (ASO) examination. Master the advanced thermodynamics and heat transfer mechanics governing aircraft thermal management, turbine blade cooling, avionics thermal dissipation, transient cryogenic fuel soaking, and critical insulation radii under severe flight operating envelopes.

Estimated Study Time: 4.5 Hours
6 Exhaustive Concept Cards + 22 Practice MCQs
DGCA CAR Section 2 & ICAO Annex 8 Compliant

Pedagogical Note for UPSC ASO Candidates

The following exhaustive study notes cover the exact subtopics mandated for Day 19. Every single subtopic below is presented as a standalone, rigorously structured concept card containing formal definitions, molecular/continuum mechanics, full mathematical derivations, dimensional analysis, aviation reference data, real aircraft applications, UPSC traps, and step-by-step worked numerics. Read every section thoroughly.

01

Critical Radius of Insulation

Subtopic 1 of 6

1. Formal Technical Definition

The critical radius of insulation ($r_{cr}$) is defined as the specific outer radius of a cylindrical or spherical insulation layer applied to a bare circular or spherical metal conduit (such as an aircraft bleed air duct or cryogenic fuel line) for which the total heat transfer rate from the system reaches an absolute maximum. Adding insulation up to $r_{cr}$ increases heat loss, whereas insulation added beyond $r_{cr}$ successfully suppresses heat loss.

2. Basic Concept & Physical Mechanism

At the continuum level, heat transfer from a cylindrical pipe to ambient air involves two competing thermal resistances in series: (i) Conduction resistance through the cylindrical insulation layer ($\ln(r_2/r_1) / 2\pi k L$), which increases as the outer radius $r_2$ increases; and (ii) Convection resistance at the outer surface boundary layer ($1 / (2\pi r_2 L h)$), which decreases as the outer radius $r_2$ increases due to the larger convective surface area $A = 2\pi r_2 L$. Initially, for very thin insulation ($r_2 < r_{cr}$), the addition of insulation causes the reduction in convection resistance to dominate over the increase in conduction resistance, leading to an increase in total heat dissipation. Only when $r_2 > r_{cr}$ does the logarithmic conduction term dominate, fulfilling the primary engineering intent of thermal insulation.

3. Mathematical Derivation & Equations

Consider a long cylindrical tube of inner radius $r_1$, outer metal radius $r_2$, insulated with an outer radius $r_3$ (letting $r_3 = r$), thermal conductivity $k$, exposed to convection coefficient $h$ and ambient temperature $T_\infty$. The steady-state heat transfer rate $\dot{Q}$ is given by:

$$\dot{Q} = \frac{T_1 - T_\infty}{R_{cond} + R_{conv}} = \frac{T_1 - T_\infty}{\frac{\ln(r/r_1)}{2\pi k L} + \frac{1}{2\pi r L h}}$$

To find the critical radius $r_{cr}$ where heat transfer is maximized, we differentiate $\dot{Q}$ with respect to the outer radius $r$ and set the derivative to zero ($\frac{d\dot{Q}}{dr} = 0$), which is equivalent to minimizing the total thermal resistance $R_{total}$:

$$\frac{d R_{total}}{dr} = \frac{d}{dr} \left[ \frac{\ln(r) - \ln(r_1)}{2\pi k L} + \frac{1}{2\pi h L r} \right] = 0$$

$$\frac{1}{2\pi k L r} - \frac{1}{2\pi h L r^2} = 0 \implies \frac{1}{k r} = \frac{1}{h r^2}$$

$$r_{cr,\text{cylinder}} = \frac{k}{h}, \quad \text{and for spheres:} \quad r_{cr,\text{sphere}} = \frac{2k}{h}$$

Where:

$k$ = Thermal conductivity of the insulation material ($\text{W}/(\text{m}\cdot\text{K})$)

$h$ = Convective heat transfer coefficient at the outer boundary ($\text{W}/(\text{m}^2\cdot\text{K})$)

$r_{cr}$ = Critical radius of insulation ($\text{m}$)

4. Dimensions & Units

Fundamental Dimensions:
[L] (Length)
SI Units:
Meters (m)
CGS Units:
Centimeters (cm)

5. Engineering Reference Data Table

Insulation Material Thermal Cond. $k$ ($\text{W}/\text{m}\cdot\text{K}$) Typical Convection $h$ ($\text{W}/\text{m}^2\cdot\text{K}$) Calculated $r_{cr}$ ($\text{mm}$)
Aircraft Glass Wool Blanket 0.038 10.0 (Free convection) 3.8 mm
Polyimide Foam (Avionics Ducting) 0.033 8.5 3.88 mm
Aerogel Blanket (Cryogenic Line) 0.021 15.0 (Forced airflow) 1.4 mm

6. Aircraft Application & Air Safety Officer Connection

In commercial aircraft, high-pressure pneumatic bleed air extracted from the compressor stages of gas turbine engines (reaching temperatures up to $450\,^\circ\text{C}$ to $600\,^\circ\text{C}$) is routed through titanium or stainless steel ducting embedded within the wing root and fuselage. According to DGCA CAR Section 2, Series 'I', Part III and FAA AC 25-something fire protection mandates, these ducts must be insulated with high-temperature blankets to protect adjacent composite structures, hydraulic lines, and electrical wiring bundles. Air Safety Implication: If maintenance engineers apply thin insulation layers whose outer radius falls below $r_{cr}$ on hot bleed air ducts, total heat rejection to the wing bay increases rather than decreases, leading to premature thermal degradation of nearby fuel tank baffling or electrical insulation, potentially causing catastrophic in-flight electrical arcs or uncontained wing bay fires.

7. UPSC Traps & Common Misconceptions

  • Trap 1: Assuming adding insulation always decreases heat transfer. (False; for $r < r_{cr}$, heat loss increases).
  • Trap 2: Confusing cylinder critical radius ($r_{cr} = k/h$) with sphere critical radius ($r_{cr} = 2k/h$). UPSC often swaps these formulas in objective matching questions.
  • Trap 3: Neglecting units when $k$ is given in $\text{W}/\text{m}\cdot\text{K}$ and $h$ in $\text{W}/\text{m}^2\cdot\text{K}$, resulting in conversion errors between meters and millimeters.

8. Step-by-Step Worked Numerical

Problem: An aircraft pneumatic bleed air steel tube having an outer radius of $r_1 = 12\,\text{mm}$ is insulated with asbestos-substitute material ($k = 0.18\,\text{W}/\text{m}\cdot\text{K}$). The tube is exposed to ambient aircraft bay air with a convective heat transfer coefficient $h = 12\,\text{W}/\text{m}^2\cdot\text{K}$. Calculate the critical radius of insulation and determine whether applying a $10\,\text{mm}$ thick layer of insulation will increase or decrease heat loss.

Step 1: Calculate Critical Radius ($r_{cr}$)

$$r_{cr} = \frac{k}{h} = \frac{0.18\,\text{W}/\text{m}\cdot\text{K}}{12\,\text{W}/\text{m}^2\cdot\text{K}} = 0.015\,\text{m} = 15\,\text{mm}$$

Step 2: Evaluate Proposed Insulation Thickness

Inner radius of insulation = tube outer radius = $12\,\text{mm}$. With a thickness of $10\,\text{mm}$, the new outer radius $r_2 = 12 + 10 = 22\,\text{mm}$.

Since the new outer radius $r_2 ($22\,\text{mm}$)$ is greater than the critical radius $r_{cr} ($15\,\text{mm}$)$, applying this $10\,\text{mm}$ insulation layer will successfully decrease the heat transfer rate from the duct.

Operational Takeaway: To ensure heat containment on bleed ducts, insulation thickness must always push the outer radius well beyond $r_{cr}$.

02

Extended Surfaces / Fins

Subtopic 2 of 6

1. Formal Technical Definition

An extended surface (or fin) is a solid metallic protrusion extending from a primary heat-transfer boundary surface into the surrounding fluid medium. Its primary engineering function is to enhance total convective heat dissipation from the base object by increasing the effective surface area available for heat transfer.

2. Basic Concept & Physical Mechanism

When a prime heat source (e.g., an aircraft piston engine cylinder head, avionics cooling chassis, or electronic control unit housing) has a relatively low convective heat transfer coefficient $h$ with ambient air, the temperature difference required to dissipate high thermal loads without fins would cause structural overheating. By attaching highly conductive metal fins (aluminum or copper), heat conducts along the fin length while simultaneously being stripped away from the fin flanks by convection. The temperature progressively drops from the base ($T_b$) to the fin tip ($T_{tip}$).

3. Mathematical Derivation & Equations

Consider a straight rectangular fin of length $L$, thickness $t$, width $W$, cross-sectional area $A_c = t W$, and perimeter $P = 2(t + W)$. Applying an energy balance on a differential fin element $dx$ under steady state:

$$\dot{Q}_{cond,x} = \dot{Q}_{cond,x+dx} + \dot{Q}_{conv}$$

$$-k A_c \frac{dT}{dx} = \left( -k A_c \frac{dT}{dx} + \frac{d}{dx}\left(-k A_c \frac{dT}{dx}\right)dx \right) + h(P\,dx)(T - T_\infty)$$

Defining excess temperature $\theta(x) = T(x) - T_\infty$ and fin parameter $m^2 = \frac{hP}{k A_c}$, we obtain the governing 1D fin differential equation:

$$\frac{d^2\theta}{dx^2} - m^2\theta = 0$$

General solution: $\theta(x) = C_1 \cosh(mx) + C_2 \sinh(mx)$. For an infinitely long fin ($L \to \infty$, $\theta(L) = 0$):

$$\theta(x) = \theta_b e^{-mx}, \quad \text{and total heat loss:} \quad \dot{Q}_{fin} = \sqrt{h P k A_c} \theta_b$$

4. Dimensions & Units

Fin Parameter ($m$):
$[L^{-1}]$ | $\text{m}^{-1}$
Heat Transfer ($\dot{Q}$):
$[M L^2 T^{-3}]$ | Watts (W)
Thermal Cond. ($k$):
$[M L T^{-3} \theta^{-1}]$ | $\text{W}/\text{m}\cdot\text{K}$

5. Engineering Reference Data Table

Fin Material Thermal Conductivity $k$ ($\text{W}/\text{m}\cdot\text{K}$) Typical Application Max Operating Temp ($^\circ\text{C}$)
Copper C11000 386 Avionic Heat Sinks 250 $^\circ\text{C}$
Aluminum Alloy 6061-T6 167 Gearbox & Oil Coolers 200 $^\circ\text{C}$
Inconel 718 (Turbine shroud fins) 11.4 Turbine Gas Paths 700 $^\circ\text{C}$

6. Aircraft Application & Air Safety Officer Connection

Finned surfaces are ubiquitous in aviation: cooling fins on piston-engine aircraft (e.g., Cessna 172 Lycoming engines), oil coolers (matrix heat exchangers), avionics cooling racks in transport category aircraft, and internal cooling passages in high-pressure turbine blades. Air Safety Implication: In turbine blade cooling, internal serpentine cooling passages use pin-fins and turbulators. If foreign object debris (FOD) or thermal oxidation blocks internal fin channels, local metal temperatures exceed the material creep limit, resulting in blade untwist, high-energy turbine blade containment failure, and catastrophic engine uncontained rotor failure (investigated under DGCA and ICAO Annex 13 guidelines).

7. UPSC Traps & Common Misconceptions

  • Trap 1: Assuming temperature gradient is uniform along the fin length. (False; temperature decays exponentially or via hyperbolic functions from base to tip).
  • Trap 2: Confusing perimeter $P$ and cross-sectional area $A_c$ in the fin parameter $m = \sqrt{hP/kA_c}$. For rectangular fins, $P = 2(w+t)$ and $A_c = wt$.

8. Step-by-Step Worked Numerical

Problem: An aluminum electronic control unit cooling fin ($k = 200\,\text{W}/\text{m}\cdot\text{K}$) is $L = 50\,\text{mm}$ long, $t = 2\,\text{mm}$ thick, and $W = 100\,\text{mm}$ wide. It is exposed to ambient air at $T_\infty = 30\,^\circ\text{C}$ with $h = 25\,\text{W}/\text{m}^2\cdot\text{K}$. The base temperature is $T_b = 120\,^\circ\text{C}$. Assuming an infinitely long fin approximation, calculate the heat dissipation rate.

Step 1: Compute Cross-sectional Area ($A_c$) and Perimeter ($P$)

$A_c = t \times W = 0.002\,\text{m} \times 0.1\,\text{m} = 2 \times 10^{-4}\,\text{m}^2$

$P = 2(t + W) = 2(0.002 + 0.1) = 0.204\,\text{m}$

Step 2: Calculate Fin Parameter ($m$)

$m = \sqrt{\frac{h P}{k A_c}} = \sqrt{\frac{25 \times 0.204}{200 \times 2 \times 10^{-4}}} = \sqrt{\frac{5.1}{0.04}} = \sqrt{127.5} = 11.29\,\text{m}^{-1}$

Step 3: Calculate Heat Loss for Infinite Fin

$\dot{Q} = \sqrt{h P k A_c}\theta_b = \sqrt{25 \times 0.204 \times 200 \times 2 \times 10^{-4}} \times (120 - 30)$

$\dot{Q} = \sqrt{1.02} \times 90 = 1.0099 \times 90 = \mathbf{90.89\,\text{W}}$

Operational Takeaway: Extended surfaces dramatically increase heat dissipation, maintaining avionics junction temperatures well within safe military/aviation operational limits.

03

Fin Efficiency & Effectiveness

Subtopic 3 of 6

1. Formal Technical Definition

Fin efficiency ($\eta_f$) is the ratio of actual heat transferred by the fin to the ideal heat that would be dissipated if the entire fin were maintained at the base temperature ($T_b$). Fin effectiveness ($\epsilon_f$) is the ratio of the total heat transfer from the finned surface to the heat transfer from the same surface if no fins were present.

2. Basic Concept & Physical Mechanism

Because finite thermal conductivity creates a temperature drop along a fin, the tip is cooler than the base, meaning the average temperature difference driving convection is less than $\theta_b$. Fin efficiency quantifies this thermal penalty ($\eta_f < 1$). Fin effectiveness ($\epsilon_f$) determines whether adding fins is even thermodynamically justified. For an effective thermal design, engineering standards dictate that $\epsilon_f$ must be significantly greater than 2.0; otherwise, the weight and manufacturing complexity of adding fins outweigh their thermal benefit.

3. Mathematical Derivation & Equations

Fin Efficiency ($\eta_f$):

$$\eta_f = \frac{\dot{Q}_{fin}}{\dot{Q}_{max}} = \frac{\dot{Q}_{fin}}{h A_{fin} (T_b - T_\infty)}$$

For a rectangular profile fin with an adiabatic tip, $\eta_f = \frac{\tanh(mL)}{mL}$.

Fin Effectiveness ($\epsilon_f$):

$$\epsilon_f = \frac{\dot{Q}_{finned}}{\dot{Q}_{unfinned}} = \frac{\dot{Q}_{fin}}{h A_b (T_b - T_\infty)} = \frac{\eta_f A_{fin}}{A_b}$$

$A_{fin}$ = Total surface area of the fin

$A_b$ = Base cross-sectional area uncovered by the fin

4. Dimensions & Units

Fin Efficiency ($\eta_f$):
Dimensionless (Ratio, typically 0 to 1 or %)
Fin Effectiveness ($\epsilon_f$):
Dimensionless (Ratio > 1 for valid design)

5. Engineering Reference Data Table

Fin Configuration Typical Efficiency ($\eta_f$) Typical Effectiveness ($\epsilon_f$) Design Assessment
Optimized Al Aircraft Oil Cooler Fin 0.85 - 0.92 5.5 - 8.0 Highly Efficient
Thick Steel Engine Cylinder Fin 0.60 - 0.70 2.1 - 3.0 Marginal (Weight penalty)
Poorly Conductive Plastic Fin < 0.20 < 1.0 Defeats purpose (Acts as insulation)

6. Aircraft Application & Air Safety Officer Connection

Engineers designing oil coolers and fuel-oil heat exchangers for turbofan engines must verify that fin effectiveness ($\epsilon_f$) remains high across all flight regimes (takeoff, climb, cruise, and descent). Air Safety Implication: If inadequate fin design or fouling by ingested runway dust reduces fin effectiveness below 2.0, lubricating oil temperature spikes during high-thrust climb. This leads to oil coking in bearing compartments, bearing seizure, and catastrophic in-flight engine failure—a major audit point during DGCA Airworthiness inspections of reciprocating and turbine powerplants.

7. UPSC Traps & Common Misconceptions

  • Trap 1: Confusing $\eta_f$ (always $\le 1.0$) with $\epsilon_f$ (can be $>> 1.0$).
  • Trap 2: Believing that making a fin infinitely long maximizes heat transfer. While total heat $\dot{Q}$ increases with length, the incremental gain approaches zero, and efficiency plummets.

8. Step-by-Step Worked Numerical

Problem: A straight rectangular aluminum fin has length $L = 20\,\text{mm}$, thickness $t = 1.5\,\text{mm}$, and thermal conductivity $k = 180\,\text{W}/\text{m}\cdot\text{K}$. It operates in an airflow where $h = 40\,\text{W}/\text{m}^2\cdot\text{K}$. Calculate the fin efficiency $\eta_f$ assuming an adiabatic tip boundary condition.

Step 1: Compute Fin Parameter $mL$

For a rectangular profile, $A_c = t \times 1$ (per unit width) and $P = 2(t + 1) \approx 2\,\text{m}$ (per unit width). Thus $m = \sqrt{h P / k A_c} = \sqrt{\frac{40 \times 2}{180 \times 0.0015}} = \sqrt{\frac{80}{0.27}} = \sqrt{296.3} = 17.21\,\text{m}^{-1}$.

$mL = 17.21\,\text{m}^{-1} \times 0.02\,\text{m} = 0.3442$

Step 2: Calculate Fin Efficiency ($\eta_f$)

$$\eta_f = \frac{\tanh(mL)}{mL} = \frac{\tanh(0.3442)}{0.3442} = \frac{0.3314}{0.3442} = \mathbf{0.963} \text{ (or } 96.3\% \text{)}$$

Operational Takeaway: A high efficiency of 96.3% indicates that the fin operates nearly isothermally with minimal temperature drop along its short length.

04

Lumped-Capacitance Method

Subtopic 4 of 6

1. Formal Technical Definition

The lumped-capacitance method is a transient heat transfer approximation in which the internal thermal resistance of a solid body is assumed to be negligible compared to the external surface convective resistance, allowing the temperature of the entire solid to be treated as uniform spatially at any given instant in time $T(t)$.

2. Basic Concept & Physical Mechanism

When a solid object with high thermal conductivity (like a solid metal turbine blade stub or steel landing gear pin) is suddenly exposed to a fluid at a different temperature, heat flows between the surface and the fluid. If internal conduction is extremely rapid, heat distributes instantly throughout the volume, eliminating internal temperature gradients. The solid behaves like a "lump" of uniform temperature, and transient analysis simplifies from a complex partial differential equation (PDE) to a straightforward ordinary differential equation (ODE).

3. Mathematical Derivation & Equations

Applying conservation of energy to a solid body of volume $V$, density $\rho$, specific heat $c$, and surface area $A_s$ undergoing convective cooling/heating:

$$-h A_s (T - T_\infty) = \rho V c \frac{dT}{dt}$$

Defining excess temperature $\theta(t) = T(t) - T_\infty$ and integrating from $t = 0$ (where $\theta(0) = \theta_i$) to $t$:

$$\frac{d\theta}{\theta} = -\frac{h A_s}{\rho V c} dt \implies \ln\left(\frac{\theta(t)}{\theta_i}\right) = -\left(\frac{h A_s}{\rho V c}\right)t$$

$$\frac{T(t) - T_\infty}{T_i - T_\infty} = e^{-(\frac{h A_s}{\rho V c})t} = e^{-t / \tau_t}$$

Where time constant $\tau_t = R_{therm} C_{therm} = \left(\frac{1}{h A_s}\right)(\rho V c)$ represents the thermal time constant of the system.

4. Dimensions & Units

Time Constant ($\tau_t$):
[T] | Seconds (s)
Specific Heat ($c$):
$[L^2 T^{-2} \theta^{-1}]$ | $\text{J}/(\text{kg}\cdot\text{K})$
Density ($\rho$):
$[M L^{-3}]$ | $\text{kg}/\text{m}^3$

5. Engineering Reference Data Table

Component / Material Density $\rho$ ($\text{kg}/\text{m}^3$) Specific Heat $c$ ($\text{J}/\text{kg}\cdot\text{K}$) Thermal Cond. $k$ ($\text{W}/\text{m}\cdot\text{K}$)
Titanium Alloy Ti-6Al-4V 4430 526.3 6.7
Aluminum 2024-T3 2770 875.0 177.0
Inconel 718 8190 435.0 11.4

6. Aircraft Application & Air Safety Officer Connection

Lumped-capacitance analysis is used by aeronautical engineers to model the thermal soak of aircraft structural fuel tanks during rapid altitude changes (e.g., supersonic cruise heating or sudden high-altitude descent into sub-freezing air) and to design thermocouple temperature probes (exhaust gas temperature - EGT probes). Air Safety Implication: EGT thermocouple response time is governed by its lumped thermal time constant. If an EGT probe has excessive thermal mass or scale buildup, its time lag during rapid throttle transients (such as go-around maneuver execution) masks engine over-temperature spikes, leading to turbine blade melting and uncontained engine failure.

7. UPSC Traps & Common Misconceptions

  • Trap 1: Applying lumped capacitance without checking the Biot number. If $Bi > 0.1$, the method is invalid and yields massive errors.
  • Trap 2: Confusing characteristic length $L_c = V/A_s$ for general bodies with $L_c = r_0/2$ for spheres and $r_0/2$ for cylinders.

8. Step-by-Step Worked Numerical

Problem: A stainless steel spherical component ($k = 15\,\text{W}/\text{m}\cdot\text{K}$, $\rho = 7800\,\text{kg}/\text{m}^3$, $c = 500\,\text{J}/\text{kg}\cdot\text{K}$) of diameter $D = 10\,\text{mm}$ is initially at $300\,^\circ\text{C}$ and is suddenly immersed in an oil bath maintained at $T_\infty = 40\,^\circ\text{C}$ with $h = 250\,\text{W}/\text{m}^2\cdot\text{K}$. Calculate the time required for the sphere to cool to $100\,^\circ\text{C}$.

Step 1: Check Biot Number Validity

$L_c = \frac{V}{A_s} = \frac{\frac{4}{3}\pi r^3}{4\pi r^2} = \frac{r}{3} = \frac{0.0025}{3} = 8.33 \times 10^{-4}\,\text{m}$

$Bi = \frac{h L_c}{k} = \frac{250 \times 8.33 \times 10^{-4}}{15} = 0.0139 < 0.1$. Lumped capacitance is valid!

Step 2: Calculate Time Constant ($\tau_t$)

$\tau_t = \frac{\rho V c}{h A_s} = \frac{7800 \times \left(\frac{0.005}{3}\right) \times 500}{250} = \frac{6500}{250} = 26.0\,\text{s}$

Step 3: Calculate Time $t$ for $T = 100\,^\circ\text{C}$

$\frac{T(t) - T_\infty}{T_i - T_\infty} = e^{-t / \tau_t} \implies \frac{100 - 40}{300 - 40} = e^{-t / 26} \implies \frac{60}{260} = 0.2308 = e^{-t / 26}$

$-t / 26 = \ln(0.2308) = -1.466 \implies t = 1.466 \times 26 = \mathbf{38.12\,\text{s}}$

Operational Takeaway: The component reaches $100\,^\circ\text{C}$ in approximately 38.1 seconds under forced oil quenching.

05

Biot & Fourier Numbers

Subtopic 5 of 6

1. Formal Technical Definition

The Biot number ($Bi$) is a dimensionless parameter representing the ratio of internal conductive thermal resistance of a solid body to the external convective thermal resistance at its surface. The Fourier number ($Fo$) is a dimensionless time parameter representing the ratio of heat conducted through a solid to the thermal energy stored within it (dimensionless transient time).

2. Basic Concept & Physical Mechanism

In transient heat conduction, these two numbers govern similarity and regime selection. If $Bi \le 0.1$, internal resistance is negligible compared to surface convection, validating the lumped-capacitance model. If $Bi > 0.1$, internal temperature gradients cannot be ignored, requiring full transient analytical solutions (Heisler charts or Bessel function series). The Fourier number acts as a dimensionless clock for transient diffusion, where systems with identical $Fo$ and $Bi$ experience identical fractional temperature distributions regardless of their physical size or material.

3. Mathematical Derivation & Equations

Biot Number ($Bi$):

$$Bi = \frac{R_{cond}}{R_{conv}} = \frac{L_c / (k A_s)}{1 / (h A_s)} = \frac{h L_c}{k}$$

Fourier Number ($Fo$):

$$Fo = \frac{\text{Heat Conduction Rate}}{\text{Thermal Energy Storage Rate}} = \frac{\alpha t}{L_c^2} = \frac{k t}{\rho c L_c^2}$$

$\alpha = \frac{k}{\rho c}$ = Thermal diffusivity ($\text{m}^2/\text{s}$)

$L_c = \frac{V}{A_s}$ = Characteristic length ($\text{m}$)

$t$ = Time ($\text{s}$)

4. Dimensions & Units

Biot Number ($Bi$):
Dimensionless ($[M^0 L^0 T^0 \theta^0]$)
Fourier Number ($Fo$):
Dimensionless ($[M^0 L^0 T^0 \theta^0]$)

5. Engineering Reference Data Table

Transient Regime Biot Number Condition Analytical Solution Method Aviation Example
Lumped System $Bi \le 0.1$ Exponential decay (ODE) Thin thermocouple probe in duct
Spatial Gradients Required $Bi > 0.1$ Heisler Charts / Series Solution Thick turbine rotor disc during start
Quasi-Steady Transient $Fo > 0.2$ First-term series approximation Wing spar heat soak during climb

6. Aircraft Application & Air Safety Officer Connection

During high-speed aircraft acceleration (e.g., military fighter jets transitioning to supersonic Mach numbers), the skin friction heating creates severe transient thermal stresses in thick titanium wing spars and leading-edge structures. Air Safety Implication: Air Safety Officers investigating structural fatigue failures must analyze the cumulative thermal fatigue cycles ($Fo$ and $Bi$ variations) experienced by load-bearing metal structures. High thermal gradients ($Bi > 0.1$) coupled with cyclic aerodynamic heating induce low-cycle thermal fatigue, leading to stress corrosion cracking and catastrophic structural delamination in flight.

7. UPSC Traps & Common Misconceptions

  • Trap 1: Confusing the numerator and denominator of the Biot number. Remember: $Bi = h L_c / k$ (Internal Conduction resistance denominator vs External Convection numerator). If inverted, it is known as the Nusselt number ($Nu = h L_c / k_f$).
  • Trap 2: Assuming Biot number uses fluid thermal conductivity $k_f$. It uses the solid body's thermal conductivity $k_s$.

8. Step-by-Step Worked Numerical

Problem: A large stainless steel plate ($k = 15\,\text{W}/\text{m}\cdot\text{K}$, $\alpha = 3.95 \times 10^{-6}\,\text{m}^2/\text{s}$) of thickness $2L = 40\,\text{mm}$ ($L = 0.02\,\text{m}$) is initially at $20\,^\circ\text{C}$ and subjected to hot exhaust gases at $500\,^\circ\text{C}$ on both sides with $h = 100\,\text{W}/\text{m}^2\cdot\text{K}$. Calculate the Biot number and the Fourier number after $t = 300\,\text{seconds}$.

Step 1: Calculate Biot Number ($Bi$)

For a plane wall of half-thickness $L$, characteristic length $L_c = L = 0.02\,\text{m}$.

$$Bi = \frac{h L_c}{k} = \frac{100 \times 0.02}{15} = \frac{2}{15} = \mathbf{0.1333}$$

Since $Bi > 0.1$, lumped capacitance is not valid; internal temperature gradients are significant.

Step 2: Calculate Fourier Number ($Fo$)

$$Fo = \frac{\alpha t}{L^2} = \frac{3.95 \times 10^{-6} \times 300}{(0.02)^2} = \frac{1.185 \times 10^{-3}}{0.0004} = \mathbf{2.9625}$$

Operational Takeaway: $Fo = 2.96$ indicates that transient heat diffusion has progressed deeply into the plate interior after 5 minutes.

06

Transient Fundamentals

Subtopic 6 of 6

1. Formal Technical Definition

Transient heat conduction (unsteady-state conduction) refers to thermal energy propagation within a solid medium where temperatures vary with both spatial coordinates and time ($\frac{\partial T}{\partial t} \neq 0$), governed by the generalized Fourier-Kirchhoff heat diffusion equation without steady-state assumptions.

2. Basic Concept & Physical Mechanism

Unlike steady-state heat transfer where energy entering equals energy leaving every volume element instantly, transient conduction involves energy storage. As thermal waves penetrate a solid, a portion of the heat flux increases the internal energy of the material lattice (manifesting as temperature rise), while the remainder is conducted further. The transient response depends fundamentally on thermal diffusivity $\alpha = k / (\rho c)$, which measures the material's ability to conduct thermal energy relative to its ability to store thermal energy.

3. Mathematical Derivation & Equations

Starting from the general 3D unsteady heat diffusion equation with internal heat generation $\dot{q}$:

$$\nabla \cdot (k \nabla T) + \dot{q} = \rho c \frac{\partial T}{\partial t}$$

For constant thermal conductivity $k$ and no internal heat generation ($\dot{q} = 0$), this simplifies to the classical diffusion equation:

$$\nabla^2 T = \frac{1}{\alpha} \frac{\partial T}{\partial t}$$

In 1D Cartesian coordinates: $\frac{\partial^2 T}{\partial x^2} = \frac{1}{\alpha} \frac{\partial T}{\partial t}$. Solutions require one initial condition ($T(x,0) = T_i$) and two boundary conditions (e.g., convective surface condition $-k \frac{\partial T}{\partial x}\big|_{x=L} = h(T(L,t) - T_\infty)$).

4. Dimensions & Units

Thermal Diffusivity ($\alpha$):
$[L^2 T^{-1}]$ | $\text{m}^2/\text{s}$
Temperature Gradient ($\frac{\partial T}{\partial x}$):
$[\theta L^{-1}]$ | $\text{K}/\text{m}$
Unsteady Rate ($\frac{\partial T}{\partial t}$):
$[\theta T^{-1}]$ | $\text{K}/\text{s}$

5. Engineering Reference Data Table

Material Thermal Diffusivity $\alpha$ ($\text{m}^2/\text{s} \times 10^{-6}$) Transient Thermal Response
Pure Silver 149.0 Extremely rapid thermal propagation
Aluminum 2024-T3 73.4 Rapid thermal equalization
Stainless Steel 304 3.95 Slow thermal diffusion (high thermal inertia)
Aircraft Epoxy Composite 0.52 Severe thermal barrier / slow soak

6. Aircraft Application & Air Safety Officer Connection

Transient heat conduction is critical in analyzing aircraft fire penetration resistance (e.g., FAA/DGCA firewall certification tests where engine nacelle firewalls must withstand $1100\,^\circ\text{C}$ flame penetration for 15 minutes without flame breakthrough) and de-icing/anti-icing boot thermal cycles on aircraft wings. Air Safety Implication: If composite wing structures undergo rapid transient heating during lightning strikes or ground de-icing fluid chemical reactions, internal thermal expansion mismatches cause subsurface delamination. Air Safety Officers must verify non-destructive testing (NDT) inspection intervals for composite aircraft following severe lightning strike events.

7. UPSC Traps & Common Misconceptions

  • Trap 1: Confusing thermal conductivity $k$ with thermal diffusivity $\alpha$. Conductivity measures steady-state heat flow ability; diffusivity measures how fast transient temperature changes propagate through a medium.
  • Trap 2: Assuming initial conditions can be ignored in transient PDEs. Initial uniform temperature $T_i$ is mandatory for non-homogeneous solution convergence.

8. Step-by-Step Worked Numerical

Problem: A semi-infinite aircraft titanium firewall slab ($\alpha = 6.5 \times 10^{-6}\,\text{m}^2/\text{s}$) is initially at $25\,^\circ\text{C}$. Suddenly, the surface is exposed to a jet engine fire torch at $800\,^\circ\text{C}$. Calculate the temperature at a depth of $x = 10\,\text{mm}$ into the firewall after $t = 60\,\text{seconds}$. (Use the error function solution: $\frac{T(x,t) - T_i}{T_s - T_i} = 1 - \text{erf}\left(\frac{x}{2\sqrt{\alpha t}}\right)$, given $\text{erf}(0.327) \approx 0.35$).

Step 1: Calculate Similarity Variable ($\eta = \frac{x}{2\sqrt{\alpha t}}$)

$$\sqrt{\alpha t} = \sqrt{6.5 \times 10^{-6} \times 60} = \sqrt{3.9 \times 10^{-4}} = 0.01975\,\text{m}$$

$$\eta = \frac{0.010}{2 \times 0.01975} = \frac{0.010}{0.0395} = 0.2532 \approx 0.253$$

Step 2: Evaluate Temperature via Error Function

For $\eta \approx 0.253$, $\text{erf}(0.253) \approx 0.280$.

$$\frac{T(x,t) - 25}{800 - 25} = 1 - \text{erf}(0.253) = 1 - 0.280 = 0.720$$

$$T(x,t) - 25 = 0.720 \times 775 = 558\,\text{deg}$$

$$T(10\,\text{mm}, 60\,\text{s}) = 558 + 25 = \mathbf{583\,^\circ\text{C}}$$

Operational Takeaway: After just 60 seconds of fire exposure, the titanium firewall reaches $583\,^\circ\text{C}$ at $10\,\text{mm}$ depth, highlighting the severity of nacelle fire containment requirements.

Active Recall Prompts: Can You Explain These Without Notes?

Prompt 1

Why does adding insulation to a pipe initially increase heat loss if the outer radius is less than the critical radius $r_{cr}$?

Prompt 2

What physical condition must the Biot number satisfy for the lumped-capacitance method to be valid, and why do internal temperature gradients matter?

UPSC Interview Corner: Board Viva Model Answers

Q: As an Air Safety Officer investigating an engine uncontained failure, how would you evaluate turbine blade cooling fin degradation? ▼

Model Answer: "I would examine the metallurgical microstructure of the failed turbine blade sections using scanning electron microscopy (SEM) to check for high-temperature creep voiding and thermal fatigue cracking. Simultaneously, I would inspect the internal serpentine cooling passages for particulate clogging or oxidation scale buildup that would impair fin effectiveness ($\epsilon_f$). A reduction in internal convective cooling efficiency causes local metal temperatures to exceed superalloy limits, resulting in loss of yield strength, blade untwist, and ultimate structural failure."

Day 19 Study Completion Checklist