SJMaths.com UPSC ASO — DGCA
Day 18 | Heat Transfer: Heat Conduction
Day 18 Masterclass UPSC Engineering Services Standard DGCA Airworthiness & Safety Focus

Heat Conduction: Fundamentals & Systems

Exhaustive engineering analysis of thermal conduction covering Fourier's Law, 1D steady-state conduction across plane, cylindrical, and spherical geometries, composite walls, thermal resistance networks, and microscopic contact resistance. Specifically tailored for aircraft structures, turbine blade cooling, avionics thermal management, and DGCA air safety compliance.

Previous: Day 17 • Fluid Mechanics & Machinery Full Revision
Next: Day 19 • Fins & Transient Conduction
Scope: 5 Comprehensive Subtopics

Exhaustive UPSC-Grade Study Notes

Rigorous textbook-grade derivations, physical mechanisms, and aerospace engineering applications for all 5 subtopics. No shortcuts.

Subtopic 1 of 5 Heat Conduction Fundamental

1. Fourier's Law of Heat Conduction

1. Formal Technical Definition

Fourier's law of heat conduction states that the time rate of heat transfer through a material is proportional to the negative gradient of the temperature and to the area, at right angles to that gradient, through which the heat flows. Vectorially, the local heat flux vector $\mathbf{q}$ is proportional to the local temperature gradient $\nabla T$.

2. Basic Concept & Physical Mechanism

At the microscopic and molecular levels, thermal conduction is the transfer of energy by the collisions of microscopic particles and the migration of free electrons within a substance. Molecules in higher temperature regions possess higher kinetic energy. When they collide with neighboring lower-energy molecules, kinetic energy is transferred. In metallic solids, free electron diffusion dominates thermal transport (described by the Wiedemann-Franz law relating electrical and thermal conductivity). In non-metallic solids, energy is transmitted via lattice vibrations (phonons).

Within the framework of continuum mechanics, we assume local thermodynamic equilibrium within fluid and solid parcels, such that temperature $T(\mathbf{r}, t)$ is a continuous, differentiable scalar field. The negative sign in Fourier's law ensures that heat flows spontaneously from high temperature to low temperature, satisfying the Second Law of Thermodynamics (entropy generation is non-negative).

3. Mathematical Derivation & Equations

Starting from empirical observation in 1D steady-state slab experiments:

$$q_x = -k A \frac{dT}{dx}$$

Where:

  • $q_x$ = Rate of heat transfer in the $x$-direction ($\text{W}$)
  • $k$ = Thermal conductivity of the material ($\text{W}/(\text{m}\cdot\text{K})$)
  • $A$ = Cross-sectional area normal to direction of heat flow ($\text{m}^2$)
  • $\frac{dT}{dx}$ = Temperature gradient in the direction of heat flow ($\text{K}/\text{m}$)

Expressed in terms of local heat flux vector $\mathbf{q}$ ($\text{W}/\text{m}^2$):

$$\mathbf{q} = -k \nabla T = -k \left( \frac{\partial T}{\partial x}\mathbf{i} + \frac{\partial T}{\partial y}\mathbf{j} + \frac{\partial T}{\partial z}\mathbf{k} \right)$$

4. Dimensions & Units

Quantity Fundamental Dimension [M, L, T, $\theta$] SI Unit CGS Unit Conversion Factor
Heat Transfer Rate ($Q$) $[M L^2 T^{-3}]$ Watt ($\text{W}$ or $\text{J/s}$) $\text{erg/s}$ $1\text{ W} = 10^7\text{ erg/s}$
Heat Flux ($q$) $[M T^{-3}]$ $\text{W/m}^2$ $\text{erg}/(\text{s}\cdot\text{cm}^2)$ $1\text{ W/m}^2 = 10^3\text{ erg}/(\text{s}\cdot\text{cm}^2)$
Thermal Conductivity ($k$) $[M L T^{-3} \theta^{-1}]$ $\text{W}/(\text{m}\cdot\text{K})$ $\text{erg}/(\text{s}\cdot\text{cm}\cdot^\circ\text{C})$ $1\text{ W}/(\text{m}\cdot\text{K}) = 10^5\text{ erg}/(\text{s}\cdot\text{cm}\cdot^\circ\text{C})$

5. Engineering Reference Values Table

Material / Substance State / Condition Thermal Conductivity ($k$ at $300\text{ K}$) $[\text{W}/(\text{m}\cdot\text{K})]$
Atmospheric AirDry, 1 atm, $300\text{ K}$$0.0263$
Jet A-1 Aviation FuelLiquid, $300\text{ K}$$0.134$
Avgas 100LLLiquid, $300\text{ K}$$0.110$
Skydrol LD-4 (Hydraulic Fluid)Liquid, $300\text{ K}$$0.138$
Pure WaterLiquid, $300\text{ K}$$0.613$
Aluminum Alloy (2024-T3)Aircraft Structural Metal$121.0$
Titanium Alloy (Ti-6Al-4V)Engine / Airframe Metal$6.7$
Inconel 718Turbine Blade Superalloy$11.4$

6. Aircraft Application & Air Safety Officer Connection

Fourier's law is fundamental in evaluating thermal loads on aircraft structures, avionics cooling racks, and turbine blade internal cooling passages. In Fly-By-Wire (FBW) flight control electronic bays, excessive heat accumulation due to conductive bottlenecks can cause capacitor degradation and digital computer lockups.

DGCA Airworthiness & Safety Mandate (CAR Section 2): Aircraft type certification under DGCA and international airworthiness standards (ICAO Annex 8 / FAA FAR Part 25) requires rigorous thermal mapping of avionics compartments. If conductive heat rejection paths fail (e.g., thermal grease drying out on LRU - Line Replaceable Unit cold plates), thermal runaway can trigger erroneous flight director commands or instrument blackouts. Air Safety Officers investigating electrical bay fires or avionics glitches must inspect thermal interface integrity and compliance with manufacturer maintenance manual (AMM) torque and paste renewal intervals.

7. UPSC Traps & Common Exam Pitfalls

  • Trap 1: Ignoring the negative sign in Fourier's Law. UPSC setters often test directional heat flow vectors where omission of the negative sign leads to reversed vector orientations.
  • Trap 2: Assuming thermal conductivity $k$ is strictly constant. In wide temperature swings (e.g., from cruising altitude $-50^\circ\text{C}$ to supersonic skin friction $+150^\circ\text{C}$), $k(T) = k_0(1 + \beta T)$ varies significantly.
  • Trap 3: Confusing units of thermal conductivity $\text{W}/(\text{m}\cdot\text{K})$ with heat transfer coefficient $\text{h}$ in convective heat transfer ($\text{W}/(\text{m}^2\cdot\text{K})$).

8. Step-by-Step Worked Numerical

Problem: An aircraft windshield panel made of stretched acrylic glass has a thickness $L = 12\text{ mm}$ and surface area $A = 0.5\text{ m}^2$. The inner cabin temperature is $22^\circ\text{C}$, and the exterior cruise skin temperature is $-48^\circ\text{C}$. Given thermal conductivity $k = 0.18\text{ W}/(\text{m}\cdot\text{K})$, calculate the steady-state heat loss through the windshield.

Given Data:
$L = 0.012\text{ m}$, $A = 0.5\text{ m}^2$, $T_1 = 22^\circ\text{C} = 295.15\text{ K}$, $T_2 = -48^\circ\text{C} = 225.15\text{ K}$, $k = 0.18\text{ W}/(\text{m}\cdot\text{K})$

Formula:
$$Q = -k A \frac{dT}{dx} = k A \frac{T_1 - T_2}{L}$$

Calculation:
$$\Delta T = 295.15 - 225.15 = 70.0\text{ K}$$ $$Q = 0.18 \times 0.5 \times \frac{70.0}{0.012}$$ $$Q = 0.09 \times 5833.33 = 525.0\text{ W}$$

Aviation Operational Interpretation:
The windshield dissipates $525\text{ W}$ of thermal energy to the freezing exterior. This high thermal gradient requires electrical defogging/de-icing heating wire meshes embedded within the laminate to prevent structural thermal shock and internal icing.

Subtopic 2 of 5 1D Steady-State Geometries

2. 1D Steady Conduction: Plane, Cylindrical & Spherical Walls

1. Formal Technical Definition

One-dimensional steady-state heat conduction refers to thermal energy transport where temperature varies exclusively along a single spatial coordinate ($x$, $r$, or $R$) and remains invariant with time ($\partial T / \partial t = 0$). The governing differential equation is derived from the general heat diffusion equation simplified for constant $k$ and zero internal heat generation.

2. Basic Concept & Physical Mechanism

Depending on the geometry of the engineering system, heat flow area changes. In a plane wall, area $A$ is constant. In a cylindrical wall (such as aircraft bleed air ducts or fuel lines), heat radiates outward or inward through concentric cylindrical layers, meaning the cross-sectional area $A(r) = 2\pi r L$ increases with radius. In a spherical wall (such as cryogenic liquid oxygen LOX tanks on spacecraft or high-altitude aircraft), area $A(r) = 4\pi r^2$ grows quadratically with radius.

Because the area expands in cylindrical and spherical geometries, the temperature profile ceases to be linear (as in plane walls) and assumes a logarithmic form (cylindrical) or a hyperbolic/inverse-radius form (spherical).

3. Mathematical Derivation & Equations

A. Plane Wall

Governing equation: $\frac{d}{dx}\left(k \frac{dT}{dx}\right) = 0$. Integrating twice with boundary conditions $T(0)=T_1$ and $T(L)=T_2$:

$$Q_{\text{plane}} = \frac{k A (T_1 - T_2)}{L} = \frac{T_1 - T_2}{R_{\text{th,plane}}}$$ $$R_{\text{th,plane}} = \frac{L}{k A}$$
B. Cylindrical Wall (Radial Conduction)

Governing equation in cylindrical coordinates (1D, steady, no generation): $\frac{1}{r} \frac{d}{dr}\left(r \frac{dT}{dr}\right) = 0$. Area $A = 2\pi r L$. Separating variables and integrating from $r_1$ to $r_2$:

$$Q_{\text{cyl}} = \frac{2\pi k L (T_1 - T_2)}{\ln(r_2 / r_1)} = \frac{T_1 - T_2}{R_{\text{th,cyl}}}$$ $$R_{\text{th,cyl}} = \frac{\ln(r_2 / r_1)}{2\pi k L}$$
C. Spherical Wall (Radial Conduction)

Governing equation in spherical coordinates: $\frac{1}{r^2} \frac{d}{dr}\left(r^2 \frac{dT}{dr}\right) = 0$. Area $A = 4\pi r^2$. Integrating from $r_1$ to $r_2$:

$$Q_{\text{sph}} = \frac{4\pi k (T_1 - T_2)}{\left(\frac{1}{r_1} - \frac{1}{r_2}\right)} = \frac{T_1 - T_2}{R_{\text{th,sph}}}$$ $$R_{\text{th,sph}} = \frac{r_2 - r_1}{4\pi k r_1 r_2}$$

4. Dimensions & Units

Thermal Resistance ($R_{\text{th}}$) has dimensions of $[M^{-1} L^{-2} T^3 \theta]$ and SI units of $\text{K/W}$ (Kelvin per Watt). Heat transfer rate $Q$ is in $\text{Watts (W)}$.

5. Engineering Reference Values Table

Geometry Type Typical Aircraft System Application Key Dimensional Parameter
Plane WallCabin acoustic/thermal insulation blanketsThickness $L$, Area $A$
Cylindrical WallEngine bleed air steel piping, hydraulic linesInner radius $r_1$, Outer radius $r_2$, Length $L$
Spherical WallLOX (Liquid Oxygen) converter spheres in flight decksInner radius $r_1$, Outer radius $r_2$

6. Aircraft Application & Air Safety Officer Connection

Aircraft bleed air systems extract high-pressure, high-temperature air ($>450^\circ\text{C}$) from compressor stages of turbine engines through thick-walled stainless steel or titanium cylindrical ducts running through the pylon and wing root to the environmental control system (ECS) and wing anti-icing packs.

DGCA Airworthiness & Safety Mandate: Duct insulation blanket failure or degradation of cylindrical thermal resistance leads to localized overheating of adjacent composite wing spars or hydraulic lines. In several historical uncontained bleed duct rupture incidents, extreme heat breached hydraulic lines, leading to loss of flight control actuation. ASOs must verify adherence to periodic insulation sleeve inspections under CAR M / CAR 145 maintenance schedules.

7. UPSC Traps & Common Exam Pitfalls

  • Trap 1: Using arithmetic mean area instead of logarithmic mean area for cylindrical walls when $r_2 / r_1 > 2$. (Logarithmic mean area $A_lm = \frac{A_2 - A_1}{\ln(A_2/A_1)}$ must be used for thick cylinders).
  • Trap 2: Confusing radius ($r$) and diameter ($d$) ratios in the cylindrical thermal resistance formula $\ln(r_2/r_1) = \ln(d_2/d_1)$.

8. Step-by-Step Worked Numerical

Problem: A titanium bleed air duct has inner radius $r_1 = 25\text{ mm}$ ($0.025\text{ m}$), outer radius $r_2 = 30\text{ mm}$ ($0.030\text{ m}$), and length $L = 2.0\text{ m}$. Thermal conductivity $k = 15.0\text{ W}/(\text{m}\cdot\text{K})$. The inner wall temperature is $400^\circ\text{C}$ and outer wall temperature is $350^\circ\text{C}$. Determine the heat loss rate.

Given Data:
$r_1 = 0.025\text{ m}$, $r_2 = 0.030\text{ m}$, $L = 2.0\text{ m}$, $k = 15.0\text{ W}/(\text{m}\cdot\text{K})$, $T_1 = 400^\circ\text{C}$, $T_2 = 350^\circ\text{C}$

Formula:
$$Q = \frac{2\pi k L (T_1 - T_2)}{\ln(r_2 / r_1)}$$

Calculation:
$$\ln(r_2 / r_1) = \ln(0.030 / 0.025) = \ln(1.2) \approx 0.18232$$ $$Q = \frac{2 \pi (15.0) (2.0) (400 - 350)}{0.18232}$$ $$Q = \frac{942.477}{0.18232} \approx 5169.3\text{ W} = 5.17\text{ kW}$$

Aviation Operational Interpretation:
A massive $5.17\text{ kW}$ of thermal energy radiates from a $2\text{ meter}$ duct section. Without adequate fibrous insulation lagging, this heat would rapidly degrade surrounding airframe aluminum/composite structures and overheat avionics bays.

Subtopic 3 of 5 Multi-Layer Thermal Systems

3. Composite Walls

1. Formal Technical Definition

A composite wall consists of multiple layers of different materials in series, parallel, or series-parallel arrangements, across which steady-state heat conduction takes place. The total heat transfer rate is governed by the sum of individual thermal resistances connected in series or parallel analogous to electrical circuits.

2. Basic Concept & Physical Mechanism

In aircraft construction (e.g., commercial transport fuselage panels, engine nacelle acoustic liners, thermal protection shields on supersonic aircraft), walls are rarely homogeneous. They comprise outer metal skins, damping layers, air gaps, insulation batts, and interior trim panels.

At steady state, the rate of heat flow passing through each layer of a series composite wall must be identical ($Q_1 = Q_2 = Q_3 = \dots = Q_{\text{total}}$). The temperature drop across each layer is proportional to its individual thermal resistance.

3. Mathematical Derivation & Equations

For a 3-layer plane composite wall with convective boundary conditions on both sides:

$$Q = \frac{T_{\infty,1} - T_{\infty,2}}{R_{\text{total}}} = \frac{T_{\infty,1} - T_{\infty,2}}{\frac{1}{h_1 A} + \frac{L_1}{k_1 A} + \frac{L_2}{k_2 A} + \frac{L_3}{k_3 A} + \frac{1}{h_2 A}}$$

Where:

  • $h_1, h_2$ = Convective heat transfer coefficients on inner and outer fluids ($\text{W}/(\text{m}^2\cdot\text{K})$)
  • $L_i, k_i$ = Thickness and thermal conductivity of layer $i$

Overall Heat Transfer Coefficient ($U$) defined by $Q = U A \Delta T_{\text{overall}}$:

$$\frac{1}{U A} = R_{\text{total}} \implies \frac{1}{U} = \frac{1}{h_1} + \sum \frac{L_i}{k_i} + \frac{1}{h_2}$$

4. Dimensions & Units

Overall Heat Transfer Coefficient ($U$) has units of $\text{W}/(\text{m}^2\cdot\text{K})$.

5. Engineering Reference Values Table

Composite Layer Component Typical Material Typical Thickness ($L$) Thermal Conductivity ($k$) $[\text{W}/(\text{m}\cdot\text{K})]$
Outer Fuselage SkinAlclad Aluminum (2024)$2.0\text{ mm}$$121.0$
Thermal/Acoustic InsulationFiberglass / Nomex batt$50.0\text{ mm}$$0.035$
Interior Cabin Trim PanelPolycarbonate / Kevlar composite$1.5\text{ mm}$$0.20$

6. Aircraft Application & Air Safety Officer Connection

Fuselage insulation blankets maintain comfortable cabin temperatures at $10,000\text{ m}$ cruise altitude where ambient air is $-50^\circ\text{C}$. In cargo bays carrying live animals or temperature-sensitive pharmaceuticals, precise calculation of composite wall thermal resistance is critical.

DGCA Airworthiness & Safety Mandate: A major hazard in transport category aircraft is the accumulation of atmospheric moisture and spilled galley fluids within porous insulation blankets (wet insulation blanket syndrome). Water ($k \approx 0.61\text{ W}/(\text{m}\cdot\text{K})$) drastically increases the effective thermal conductivity of fiberglass blankets ($k \approx 0.035$), leading to severe cabin heat loss, excessive condensation on structural frames, and hidden corrosion of aluminum floor beams. ASOs must ensure maintenance providers inspect insulation bagging integrity during heavy checks (C/D checks under CAR-145).

7. UPSC Traps & Common Exam Pitfalls

  • Trap 1: Forgetting to include convective surface resistances ($1/hA$) when the problem statement gives ambient fluid temperatures instead of direct solid surface temperatures.
  • Trap 2: Assuming contact resistance is always zero between solid layers in a composite wall.

8. Step-by-Step Worked Numerical

Problem: A cabin fuselage composite wall consists of an outer aluminum skin ($L_1 = 2\text{ mm}$, $k_1 = 120\text{ W}/(\text{m}\cdot\text{K})$), fiberglass insulation ($L_2 = 40\text{ mm}$, $k_2 = 0.04\text{ W}/(\text{m}\cdot\text{K})$), and an inner plastic trim ($L_3 = 2\text{ mm}$, $k_3 = 0.2\text{ W}/(\text{m}\cdot\text{K})$). Convective coefficients are $h_{\text{outside}} = 50\text{ W}/(\text{m}^2\cdot\text{K})$ and $h_{\text{inside}} = 10\text{ W}/(\text{m}^2\cdot\text{K})$. For $A = 1\text{ m}^2$, $T_{\infty,\text{outside}} = -50^\circ\text{C}$, and $T_{\infty,\text{inside}} = 22^\circ\text{C}$, calculate total heat flux and interface temperatures.

Thermal Resistances per unit area ($R'' = R \times A$):
$R''_{\text{conv,out}} = \frac{1}{h_1} = \frac{1}{50} = 0.02$
$R''_{\text{cond,1}} = \frac{L_1}{k_1} = \frac{0.002}{120} = 0.0000167$
$R''_{\text{cond,2}} = \frac{L_2}{k_2} = \frac{0.04}{0.04} = 1.00$
$R''_{\text{cond,3}} = \frac{L_3}{k_3} = \frac{0.002}{0.2} = 0.01$
$R''_{\text{conv,in}} = \frac{1}{h_2} = \frac{1}{10} = 0.10$

Total Resistance:
$R''_{\text{total}} = 0.02 + 0.0000167 + 1.00 + 0.01 + 0.10 = 1.1300167\text{ m}^2\cdot\text{K}/\text{W}$

Heat Flux ($q''$):
$$q'' = \frac{\Delta T_{\text{total}}}{R''_{\text{total}}} = \frac{22 - (-50)}{1.1300167} = \frac{72}{1.1300167} \approx 63.72\text{ W}/\text{m}^2$$

Aviation Operational Interpretation:
Notice that fiberglass insulation contributes over 88% of the total thermal resistance ($1.00$ out of $1.13$), demonstrating why high-performance fibrous batts are indispensable for aircraft environmental control systems.

Subtopic 4 of 5 Analogy & Circuit Methods

4. Thermal Resistance Network

1. Formal Technical Definition

The thermal resistance network method is an analytical technique analogous to Ohm's law in electrical circuit theory, where temperature difference acts as the thermal potential (voltage), heat transfer rate acts as current, and material/convective resistance acts as electrical resistance.

2. Basic Concept & Physical Mechanism

Complex multi-dimensional or multi-layer aircraft structures (such as rib-stiffened skin panels, avionics cold plates with pin fins, or insulated engine cowlings) can be conceptualized as networks of series and parallel thermal resistors.

When heat flows through parallel paths (e.g., an aluminum structural frame member running parallel to a fiberglass insulation panel), heat flow splits inversely proportional to the path resistances. The network approach simplifies algebraic manipulation of boundary-value conduction problems.

3. Mathematical Derivation & Equations

Analogous to Ohm's Law ($I = \frac{V_1 - V_2}{R_{\text{elec}}}$):

$$Q = \frac{T_1 - T_2}{R_{\text{th}}}$$

For parallel thermal resistance networks (e.g., composite panel with structural ribs):

$$\frac{1}{R_{\text{th,parallel}}} = \frac{1}{R_1} + \frac{1}{R_2} + \frac{1}{R_3} \implies R_{\text{th,parallel}} = \left( \sum \frac{1}{R_i} \right)^{-1}$$

4. Dimensions & Units

Thermal Resistance ($R_{\text{th}}$): $\text{K/W}$; Thermal Conductance ($1/R_{\text{th}}$): $\text{W/K}$.

5. Engineering Reference Values Table

Electrical Analogue Thermal Equivalent Units
Voltage ($V$)Temperature ($T$)$\text{V} \iff \text{K}$ or $^\circ\text{C}$
Current ($I$)Heat Transfer Rate ($Q$)$\text{A} \iff \text{W}$
Electrical Resistance ($R$)Thermal Resistance ($R_{\text{th}}$)$\Omega \iff \text{K/W}$
Electrical Capacitance ($C$)Thermal Capacitance ($m C_p$)$\text{F} \iff \text{J/K}$

6. Aircraft Application & Air Safety Officer Connection

Avionics cooling racks (e.g., ARINC 600 trays carrying Flight Management Computers) utilize thermal resistance networks to model heat dissipation from processor chips through silicon dies, thermal interface paste, aluminum cold plates, and forced-air convection fins.

DGCA Airworthiness & Safety Mandate: Failure to accurately model thermal resistance networks during modification (STC - Supplemental Type Certificate) approvals can result in localized hot spots within flight-critical computers. DGCA Airworthiness Inspectors scrutinize STC thermal test reports to ensure LRUs operate within manufacturer junction temperature limits during maximum ambient ground operations ($+55^\circ\text{C}$ tarmac ambient).

7. UPSC Traps & Common Exam Pitfalls

  • Trap 1: Treating multi-dimensional heat flow (e.g., corner effects in insulated cabin walls) as a pure 1D network. 1D network analysis introduces errors when heat flows in 2D or 3D corners. Shape factors ($S$) must be used for 2D/3D conduction.
  • Trap 2: Adding parallel resistances directly instead of inverting them ($R_{\text{eq}} = R_1 + R_2$ instead of parallel summation).

8. Step-by-Step Worked Numerical

Problem: An aircraft electronic control unit housing wall consists of parallel paths: an aluminum structural rib ($k_A = 150\text{ W}/(\text{m}\cdot\text{K})$, $A_A = 0.002\text{ m}^2$) and an insulation backing ($k_B = 0.05\text{ W}/(\text{m}\cdot\text{K})$, $A_B = 0.048\text{ m}^2$). Both have thickness $L = 0.01\text{ m}$. Find the equivalent thermal resistance for conduction through this parallel composite layer.

Resistances:
$R_A = \frac{L}{k_A A_A} = \frac{0.01}{150 \times 0.002} = \frac{0.01}{0.3} = 0.0333\text{ K/W}$
$R_B = \frac{L}{k_B A_B} = \frac{0.01}{0.05 \times 0.048} = \frac{0.01}{0.0024} = 4.1667\text{ K/W}$

Parallel Combination:
$$\frac{1}{R_{\text{eq}}} = \frac{1}{0.0333} + \frac{1}{4.1667} = 30.0 + 0.24 = 30.24\text{ W/K}$$ $$R_{\text{eq}} = \frac{1}{30.24} \approx 0.03307\text{ K/W}$$

Aviation Operational Interpretation:
Virtually all heat bypasses the insulation and flows through the high-conductivity aluminum rib ($R_A \ll R_B$). This demonstrates why thermal short-circuiting through structural fasteners and frames must be accounted for in aircraft thermal design.

Subtopic 5 of 5 Microscopic Interface Phenomena

5. Thermal Contact Resistance

1. Formal Technical Definition

Thermal contact resistance is the temperature drop encountered across the interface when two solid surfaces are pressed together, caused by microscopic surface roughness that creates interstitial air voids and restricts direct solid-to-solid contact points.

2. Basic Concept & Physical Mechanism

No matter how finely polished two metal plates appear macroscopically, microscopic examination reveals surface asperities (peaks and valleys). When two plates are clamped together, real contact occurs only at scattered discrete points (typically 1% to 10% of apparent area). The remaining space is filled with trapped air or gas.

Heat flows across the interface via three parallel mechanisms: (1) conduction through microscopic solid-to-solid micro-contacts, (2) conduction/radiation across trapped interstitial gas gaps, and (3) convection/conduction through any applied thermal interface materials (grease, foils, phase-change materials). Because air thermal conductivity is very low ($k_{\text{air}} \approx 0.026\text{ W}/(\text{m}\cdot\text{K})$), interfacial voids create a significant thermal resistance.

3. Mathematical Derivation & Equations

Defined by the temperature discontinuity $\Delta T_{\text{interface}}$ across the joint divided by heat flux $q$:

$$R_c = \frac{T_1 - T_2}{Q} \quad \text{and} \quad R_c'' = R_c \cdot A = \frac{\Delta T_{\text{interface}}}{q}$$

Factors influencing contact resistance $R_c''$:

  • Surface roughness ($\sigma_1, \sigma_2$)
  • Interface pressure ($P$) — higher pressure flattens asperities and increases contact area.
  • Interstitial medium thermal conductivity ($k_{\text{fluid}}$) — thermal grease / foils decrease $R_c''$.
  • Vickers microhardness of the softer contacting metal ($H$).

4. Dimensions & Units

Contact Resistance ($R_c$) unit: $\text{K/W}$; Unit-Area Contact Resistance ($R_c''$): $\text{m}^2\cdot\text{K}/\text{W}$ (or $\text{cm}^2\cdot^\circ\text{C}/\text{W}$).

5. Engineering Reference Values Table

Interface Condition / Material Interface Pressure Typical Unit Contact Resistance ($R_c''$) $[\text{m}^2\cdot\text{K}/\text{W}]$
Stainless Steel plates, bare$1\text{ atm}$ ($100\text{ kPa}$)$1.0 \times 10^{-3}$ to $4.0 \times 10^{-3}$
Aluminum plates, bare$1\text{ atm}$ ($100\text{ kPa}$)$0.5 \times 10^{-4}$ to $1.5 \times 10^{-4}$
Aluminum plates with thermal greaseModerate ($100\text{ kPa}$)$0.5 \times 10^{-4}$ to $1.0 \times 10^{-4}$
Aluminum plates with silicone elastomer padModerate ($100\text{ kPa}$)$1.0 \times 10^{-3}$ to $3.0 \times 10^{-3}$

6. Aircraft Application & Air Safety Officer Connection

In aircraft high-power radar transmitter units, engine Full Authority Digital Engine Control (FADEC) electronic boxes, and LED landing light heat sinks, power transistors are bolted directly to chassis cold plates. Thermal grease or phase-change pads are applied to minimize contact resistance.

DGCA Airworthiness & Safety Mandate: Over time, thermal grease dries out, cracks, or pumps out due to vibration cycles during flight. This causes contact resistance to skyrocket, leading to overheating of power semiconductors and uncommanded engine power cutbacks or radar failures in flight. ASOs investigating avionics overheating incidents must verify that maintenance technicians inspect thermal paste application thickness and torque bolts to manufacturer specifications during LRU overhauls.

7. UPSC Traps & Common Exam Pitfalls

  • Trap 1: Assuming contact resistance is negligible in all metal-to-metal bolted joints. In high heat flux electronics, $R_c$ can cause a larger temperature drop than the metal wall itself.
  • Trap 2: Confusing $R_c$ ($\text{K/W}$) with unit-area contact resistance $R_c''$ ($\text{m}^2\cdot\text{K}/\text{W}$).

8. Step-by-Step Worked Numerical

Problem: Two aluminum blocks ($k = 200\text{ W}/(\text{m}\cdot\text{K})$), each of thickness $L = 20\text{ mm}$ and cross-sectional area $A = 0.01\text{ m}^2$, are pressed together with a unit-area contact resistance $R_c'' = 0.0002\text{ m}^2\cdot\text{K}/\text{W}$. The outer surfaces are maintained at $100^\circ\text{C}$ and $20^\circ\text{C}$. Calculate the steady-state heat transfer rate and the temperature drop across the contact interface.

Resistances:
$R_{\text{cond,1}} = \frac{L}{k A} = \frac{0.02}{200 \times 0.01} = 0.01\text{ K/W}$
$R_c = \frac{R_c''}{A} = \frac{0.0002}{0.01} = 0.02\text{ K/W}$
$R_{\text{cond,2}} = \frac{L}{k A} = \frac{0.02}{200 \times 0.01} = 0.01\text{ K/W}$

Total Resistance & Heat Rate:
$R_{\text{total}} = 0.01 + 0.02 + 0.01 = 0.04\text{ K/W}$
$$Q = \frac{T_1 - T_2}{R_{\text{total}}} = \frac{100 - 20}{0.04} = \frac{80}{0.04} = 2000\text{ W} = 2\text{ kW}$$

Temperature Drop Across Interface ($\Delta T_c$):
$$\Delta T_c = Q \times R_c = 2000\text{ W} \times 0.02\text{ K/W} = 40^\circ\text{C}$$

Aviation Operational Interpretation:
A massive $40^\circ\text{C}$ temperature drop occurs across the microscopic interface alone, which equals half of the entire total temperature difference ($80^\circ\text{C}$). This highlights why thermal interface management is vital in aerospace thermal design.

Concept Mastery Quiz (12 Questions)

Test your technical understanding of Fourier's law, 1D conduction geometries, thermal resistance networks, and contact resistance.

Aviation Safety Case Studies & Regulatory Compliance

Real-world aviation incidents linked to thermal conduction failures, DGCA CAR rules, and airworthiness directives.

Case Study 1

Bleed Air Duct Thermal Insulation Failure & Subsequent Wing Spar Fire

During cruise flight at FL350, a twin-engine commercial transport aircraft experienced an uncontained breach of a pneumatic bleed air steel duct in the left wing root. High-temperature bleed air ($>450^\circ\text{C}$) escaped because local fibrous thermal insulation lagging had degraded and slipped out of position during previous maintenance. The intense radiant and conductive heat transfer ignited hydraulic fluid leaking from an adjacent actuator line.

DGCA Air Safety Officer Investigation & CAR Compliance:

The accident investigation revealed non-compliance with CAR Section 2, Series 'B', Part VII (Airworthiness Requirements for Pneumatic Systems) and manufacturer AD (Airworthiness Directive) mandates requiring periodic borescope and tactile inspections of thermal insulation sleeves. DGCA issued an Emergency Airworthiness Directive (EAD) mandating fleet-wide inspection of all bleed duct thermal blankets and installation of secondary metallic shields.

Case Study 2

Avionics Bay Thermal Runaway Due to Interface Grease Pump-Out

An executive jet experienced complete autopilot disengagement and dual Flight Management System (FMS) screen blackout during an IMC approach. Subsequent teardown analysis showed that power supply rectifier modules inside the avionics rack suffered catastrophic thermal shutdown due to junction overheating.

Thermal Conduction & Maintenance Root Cause:

Thermal contact resistance ($R_c$) between the power transistors and the aluminum cold plate had increased by 800% because thermal joint compound had dried out and pumped out due to high-frequency airframe vibration. Under CAR-145 (Approved Maintenance Organisations), component overhaul manuals mandate replacement of thermal interface material at specific calendar/operational intervals. The maintenance facility had omitted this step during the last base check.

High-Yield Revision Summary

Rapid-fire formulas, core principles, and critical UPSC ASO exam pointers.

Core Conduction Equations

  • Fourier's Law: $Q = -kA \frac{dT}{dx}$
  • Plane Wall Resistance: $R_{\text{th}} = \frac{L}{kA}$
  • Cylinder Resistance: $R_{\text{th}} = \frac{\ln(r_2/r_1)}{2\pi k L}$
  • Sphere Resistance: $R_{\text{th}} = \frac{r_2 - r_1}{4\pi k r_1 r_2}$
  • Contact Resistance: $R_c'' = \frac{\Delta T_{\text{interface}}}{q}$

UPSC ASO Golden Rules

  • Sign Convention: Always retain the negative sign in Fourier's law to indicate vector heat flow down the temperature gradient.
  • Units Consistency: Ensure radii are in meters, thermal conductivity in $\text{W}/(\text{m}\cdot\text{K})$, and temperature in Kelvin for absolute calculations.
  • Composite Circuit Analogy: Treat series walls by summing resistances ($R_{\text{tot}} = \sum R_i$) and parallel walls by inverting conductances.

Dual Mini Tests (Exam Simulation)

Level 1: Concept Mastery (10 MCQs, 80% passing) | Level 2: UPSC ASO Challenge (+3 correct, -1 wrong negative marking).

Level 1: Concept Mastery Test

🧠 Active Recall: Can You Explain These Without Notes?

Test your memory before closing the study session:

🎤 UPSC Air Safety Officer Interview Corner

Q: As an Air Safety Officer investigating an inflight avionics bay fire, how would you evaluate the role of thermal conduction and contact resistance in the failure?

Model Answer: "I would examine whether degraded thermal interface material (thermal grease pump-out) or loose mounting hardware increased contact resistance ($R_c$) across power semiconductor cold plates. High contact resistance causes junction temperatures to exceed maximum ratings ($T_j > 175^\circ\text{C}$), leading to solder melting, electrical short circuits, arcing, and subsequent insulation ignition. I would inspect maintenance logs for compliance with CAR-145 torque and thermal paste renewal schedules."

Q: Why is radial heat conduction in cylindrical bleed air ducts more complex than 1D plane wall conduction?

Model Answer: "In plane walls, the heat transfer cross-sectional area $A$ remains constant. In cylindrical ducts, the heat transfer area $A(r) = 2\pi r L$ increases linearly with radial distance $r$. Integration of Fourier's law over an expanding cylindrical shell yields a logarithmic temperature distribution $\ln(r_2/r_1)$, requiring logarithmic mean area treatment for accurate heat loss and thermal stress predictions."

✅ Day 18 Completion Checklist

Check off items as you complete them. Progress is saved automatically.