Calculator D2

Fourier’s Law of Conduction and Thermal Resistance Networks

Heat flows through solid materials like a river flowing downhill — faster when there’s a bigger temperature difference and slower when the material is thicker or less conductive.

⚠️ Why It Matters

1
Inaccurate thermal resistance modeling
2
Over- or under-sized heat sinks
3
Localized component overheating
4
Premature electronic failure
5
System-level reliability degradation
6
Costly redesigns and warranty claims

📘 Definition

Fourier’s Law of Conduction states that the conductive heat flux vector q (W/m²) is proportional to the negative temperature gradient, expressed as q = −k∇T, where k is the thermal conductivity (W/m·K). For one-dimensional steady-state conduction across a homogeneous slab, it simplifies to Q = kΔT/δ, where Q is heat transfer rate per unit area, ΔT is the temperature difference across thickness δ. Thermal resistance networks model composite conduction paths analogously to electrical resistances in series and parallel.

🎨 Concept Diagram

T_H = 85°CT_C = 25°CQ = k·(T_H−T_C)/δδ = 0.02 m, k = 200 W/m·K

AI-generated illustration for visual understanding

💡 Engineering Insight

Never treat contact resistance as a fixed value — it varies inversely with interfacial pressure and exponentially with surface roughness. In power modules, a 20% reduction in clamping force can increase R_c by 200%, making mechanical preload as critical as thermal interface material selection. Always specify TIM application method (dispense volume, cure profile, bond line thickness) in procurement docs — not just k-value.

📖 Detailed Explanation

Fourier’s Law is the foundational principle for all conduction analysis: it defines heat flow as driven solely by temperature gradients, independent of bulk motion or radiation. The negative sign indicates heat flows from hot to cold — a directional constraint essential for correct network polarity. In simple slabs, this yields linear temperature profiles and enables direct analogy to electrical circuits, where temperature difference (ΔT) acts like voltage, heat flow (Q) like current, and resistance (R) like ohmic resistance.

Beyond planar solids, Fourier’s Law generalizes to multidimensional and transient cases via the heat diffusion equation ∂T/∂t = α∇²T, where α = k/ρc_p is thermal diffusivity. Resistance networks remain valid only when spatial gradients are predominantly one-dimensional and time constants exceed system response requirements. At interfaces, the law breaks down locally — requiring empirical correction via contact resistance models (e.g., Cooper-Mikic-Yovanovich) that incorporate surface topography, hardness, and load.

Advanced applications include anisotropic materials (e.g., graphite heat spreaders with k_∥ ≈ 400 W/m·K, k_⊥ ≈ 5 W/m·K), where resistance networks must resolve directional conductivities separately; microscale effects (Kapitza resistance at phonon boundaries); and coupled conduction-convection-radiation systems where radiation dominates above ~500°C and must be linearized (h_rad ≈ 4εσT³) to preserve network linearity. Modern practice embeds resistance networks within digital twins, updating R-values in real time using embedded thermistor data and Bayesian inference.

🔄 Engineering Workflow

Step 1
Step 1: Identify geometry, materials, and boundary conditions (temperatures, convection coefficients, heat sources)
Step 2
Step 2: Segment domain into homogeneous layers/regions and define conduction paths
Step 3
Step 3: Compute individual conductive resistances (R = δ/kA) and convective resistances (R = 1/hA)
Step 4
Step 4: Construct equivalent thermal resistance network (series/parallel) and validate topology against physics
Step 5
Step 5: Solve for unknown temperatures or heat flows using Ohm’s-law analog (Q = ΔT/R_total)
Step 6
Step 6: Verify assumptions (steady-state, 1D dominance, negligible radiation) via order-of-magnitude check or CFD screening
Step 7
Step 7: Iterate with measured interface pressures, TIM performance data, or infrared thermography validation

📋 Decision Guide

Rock/Field Condition Recommended Design Action
High-power semiconductor module (P > 100 W), air-cooled heatsink Model contact resistance explicitly using TIM thickness and pressure-dependent k; use R_total = R_silicon + R_solder + R_baseplate + R_contact + R_heatsink + 1/(hA)
Building envelope with layered insulation, structural framing, and air films Use parallel-series network: R_total = R_out + [1/(1/R_framing + 1/R_insulation) + R_sheathing + R_drywall] + R_in
Transient thermal event (e.g., motor startup, pulsed laser diode) Augment conduction-only network with lumped capacitance (C_th) elements; apply RC time constant analysis before assuming steady-state validity

📊 Key Properties & Parameters

Thermal Conductivity (k)

0.024–400 W/m·K (air: 0.024, copper: 401, concrete: 1.4, fiberglass: 0.04)

Material property quantifying its ability to conduct heat; ratio of heat flux to temperature gradient.

⚡ Engineering Impact:

Directly governs required thickness for insulation or heat spreader design — low k demands greater thickness to achieve target R-value.

Thermal Resistance (R)

0.05–10 m²·K/W (single-glazed window: ~0.17, 100 mm mineral wool: ~2.5, high-performance wall assembly: ~6.0)

Measure of a material’s or assembly’s opposition to heat flow; R = δ/k for conduction, units m²·K/W.

⚡ Engineering Impact:

Determines total heat loss/gain in building envelopes and junction temperatures in electronics packaging — errors >10% cause non-compliance with energy codes.

Contact Resistance (R_c)

1×10⁻⁵–5×10⁻³ m²·K/W (bare aluminum-aluminum: ~10⁻⁵, unlubricated steel-steel: ~10⁻³)

Additional thermal resistance at interfaces due to surface roughness, air gaps, and interfacial pressure.

⚡ Engineering Impact:

Dominates total resistance in high-power semiconductor modules — unaccounted R_c causes junction temperature overestimation by >20°C.

Convection Coefficient (h)

2–25,000 W/m²·K (natural convection air: 2–25, forced air: 10–500, water cooling: 500–10,000, boiling water: 2,500–25,000)

Empirical parameter relating surface heat flux to temperature difference between surface and fluid (q = hΔT).

⚡ Engineering Impact:

Controls boundary condition fidelity in resistance networks — misestimating h by factor of 2 shifts predicted case temperature by 30–50°C in power electronics.

📐 Key Formulas

Conductive Resistance (Slab)

R_{cond} = \frac{\delta}{k A}

Thermal resistance of a uniform solid layer of thickness δ, conductivity k, and cross-sectional area A.

Variables:
Symbol Name Unit Description
R_{cond} Conductive Resistance K/W Thermal resistance of a uniform solid layer
δ Thickness m Thickness of the slab
k Thermal Conductivity W/(m·K) Material's ability to conduct heat
A Cross-sectional Area Area perpendicular to heat flow
Typical Ranges:
Power electronics baseplate (Cu, 2 mm)
0.0005–0.002 m²·K/W
Roof insulation (fiberglass, 200 mm)
2.0–3.0 m²·K/W
⚠️ For electronics: R_cond < 0.001 m²·K/W per 100 W; for buildings: R_total ≥ 4.0 m²·K/W (IECC 2021 Climate Zone 5)

Convective Resistance

R_{conv} = \frac{1}{h A}

Thermal resistance at a solid-fluid interface governed by convection coefficient h and surface area A.

Variables:
Symbol Name Unit Description
R_{conv} Convective Resistance K/W Thermal resistance at a solid-fluid interface
h Convection Heat Transfer Coefficient W/(m2·K) Coefficient governing heat transfer rate between solid surface and fluid
A Surface Area m2 Area of the solid surface in contact with the fluid
Typical Ranges:
Natural convection on PCB (h ≈ 8 W/m²·K, A = 0.002 m²)
62.5 m²·K/W
Forced air heatsink (h ≈ 120 W/m²·K, A = 0.03 m²)
0.28 m²·K/W
⚠️ Ensure R_conv ≤ 0.3 m²·K/W for silicon devices operating >100 W to avoid exceeding T_j,max

Total Resistance (Series)

R_{total} = R_1 + R_2 + \dots + R_n

Sum of individual resistances for heat flow through sequential layers or interfaces.

Variables:
Symbol Name Unit Description
R_{total} Total Resistance K/W Total thermal resistance for heat flow through sequential layers or interfaces
R_1 Resistance of Layer 1 K/W Thermal resistance of the first layer or interface
R_2 Resistance of Layer 2 K/W Thermal resistance of the second layer or interface
R_n Resistance of Layer n K/W Thermal resistance of the nth layer or interface
Typical Ranges:
SiC module (die–solder–baseplate–TIM–heatsink)
0.15–0.45 cm²·K/W
Residential wall assembly (exterior film–sheathing–insulation–framing–drywall–interior film)
3.2–6.8 m²·K/W
⚠️ Electronics: R_total < 0.3 cm²·K/W per 100 W; Buildings: R_total ≥ prescribed by ASHRAE 90.1 or IECC for climate zone

🏭 Engineering Example

Tesla Gigafactory Berlin Power Inverter Assembly Line

N/A — applied to SiC MOSFET module thermal management
k_TIM
6.5 W/m·K
TIM_thickness
0.05 mm
Case_temp_target
95 °C
Clamping_pressure
1.2 MPa
Max_junction_temp
175 °C
R_contact_measured
0.18 cm²·K/W

🏗️ Applications

  • Power electronics thermal management
  • Building energy modeling and code compliance
  • Aerospace thermal control systems
  • Battery pack thermal design

📋 Real Project Case

Air-Cooled Condenser Retrofit for 600 MW Coal Power Plant

Retrofit of legacy water-cooled condenser at Midwest US plant

Challenge: Water scarcity forcing shift to dry cooling; risk of summer turbine backpressure rise
Read full case study →

🎨 Technical Diagrams

HotColdk=1.4δ=0.15 mQ
R₁R₂R₃ΔT₁ΔT₂ΔT₃R_total = R₁ + R₂ + R₃
R_AR_BR_C1/R_total = 1/R_A + 1/R_B + 1/R_C

📚 References

[1]
ASHRAE Handbook—Fundamentals — American Society of Heating, Refrigerating and Air-Conditioning Engineers
[2]
Thermal Management of Microelectronic Equipment — IEEE Press / John Wiley & Sons