Thermal Contact Resistance and Interfacial Conductance Models
When two solid surfaces touch, heat doesn’t flow across the gap as easily as it does through the solids themselves — that ‘resistance’ at the contact point is thermal contact resistance.
⚠️ Why It Matters
📘 Definition
Thermal contact resistance (R<sub>c</sub>) is the temperature drop per unit heat flux across the interface between two nominally solid, contacting bodies under steady-state conduction. It arises from surface roughness, interfacial voids, oxide layers, and imperfect mechanical contact, and is inversely related to interfacial conductance (h<sub>c</sub> = 1/R<sub>c</sub>). It is a non-intrinsic property dependent on contact pressure, surface finish, material pairing, and interstitial medium.
🎨 Concept Diagram
AI-generated illustration for visual understanding
💡 Engineering Insight
Thermal contact resistance is rarely dominated by bulk conduction — it’s governed by the *statistical distribution of asperities* and their elastic/plastic response. A polished surface (Ra < 0.1 μm) may underperform a moderately rough one (Ra ≈ 0.4 μm) under low pressure because the latter achieves more distributed, stable micro-contacts. Always measure — never assume — R<sub>c</sub>; even identical-looking assemblies can vary ±40% due to bolt-torque scatter and surface oxidation.
📖 Detailed Explanation
Deeper analysis reveals R<sub>c</sub> depends on three parallel mechanisms: (1) solid conduction through deformed asperities (plastic/elastic), (2) conduction through interstitial fluid/gas in the non-contact regions, and (3) radiation across gaps (significant above ~500 K). The Cooper–Mikic–Yovanovich (CMY) model combines these using statistical descriptions of surface roughness, material hardness, and contact mechanics — making it the industry standard for predictive design.
Advanced treatments incorporate time-dependent effects: thermal expansion mismatch causes contact degradation during cycling; interfacial chemical reactions (e.g., Al₂O₃ growth on aluminum) increase R<sub>c</sub> over time; and nanoscale phenomena (phonon mismatch, Kapitza resistance) dominate below 100 nm gap thicknesses. For cryogenic systems or ultra-high-flux laser diodes, molecular dynamics simulations and spectral phonon transmission models are now used alongside CMY to resolve wavelength-dependent scattering at atomically clean interfaces.
🔄 Engineering Workflow
📋 Decision Guide
| Rock/Field Condition | Recommended Design Action |
|---|---|
| High-power semiconductor (T<sub>j</sub> > 125°C) with bare copper-aluminum interface, Ra = 0.8 μm, P = 0.5 MPa | Apply conformal, low-viscosity phase-change TIM (k ≈ 6 W/m·K); target h<sub>c</sub> ≥ 15,000 W/m²·K |
| Vacuum environment, no TIM, stainless-steel–titanium interface, Ra = 2.5 μm, P = 2 MPa | Use indium foil interlayer or soft metal plating; avoid grease-based TIMs (outgassing risk) |
| High-cycle thermal cycling (>10⁴ cycles), aluminum–copper cold plate interface | Specify silver-filled electrically insulating elastomer TIM with creep recovery >90% and CTE matching |
📊 Key Properties & Parameters
Contact Pressure
0.1–100 MPaNormal compressive stress applied across the interface, influencing real contact area and deformation.
Doubling pressure typically reduces R<sub>c</sub> by 30–70% for metallic interfaces due to increased asperity deformation.
Surface Roughness (Ra)
0.02–10 μmArithmetic average deviation of surface profile from its mean line, quantifying micro-scale topography.
A tenfold increase in Ra can increase R<sub>c</sub> by 2–5× under identical pressure and material conditions.
Interfacial Conductance (h<sub>c</sub>)
100–100,000 W/m²·KReciprocal of thermal contact resistance per unit area; quantifies heat transfer efficiency at the interface.
Values < 1,000 W/m²·K indicate severe bottlenecking — often requiring thermal interface materials (TIMs) to meet junction temperature targets.
Interfacial Gap Conductivity (k<sub>gap</sub>)
0.026 (air) – 80 (solder) W/m·KEffective thermal conductivity of the material (e.g., air, grease, solder) filling microscopic gaps between surfaces.
Substituting air (k ≈ 0.026 W/m·K) with phase-change TIM (k ≈ 5–10 W/m·K) improves h<sub>c</sub> by up to 20×.
📐 Key Formulas
Cooper–Mikic–Yovanovich (CMY) Model
h_c = 0.46 \left( \frac{k_s k_g}{k_s + k_g} \right)^{0.5} \left( \frac{H}{\sigma} \right)^{0.95} \left( \frac{P_c}{H} \right)^{0.95}Empirical correlation for interfacial conductance based on surface hardness (H), roughness (σ), contact pressure (P<sub>c</sub>), and solid/gap conductivities (k<sub>s</sub>, k<sub>g</sub>).
| Symbol | Name | Unit | Description |
|---|---|---|---|
| h_c | interfacial conductance | W/(m^2·K) | Conductive heat transfer coefficient at the solid-solid interface |
| k_s | solid conductivity | W/(m·K) | Thermal conductivity of the solid material |
| k_g | gap conductivity | W/(m·K) | Thermal conductivity of the interfacial gap material (e.g., gas or vacuum) |
| H | surface hardness | MPa | Indentation hardness of the softer contacting surface |
| σ | surface roughness | m | Root-mean-square (RMS) surface roughness |
| P_c | contact pressure | Pa | Average normal pressure at the interface |
Thermal Contact Resistance (R_c)
R_c = \frac{1}{h_c A}Total resistance across an interface of area A, given interfacial conductance h_c.
| Symbol | Name | Unit | Description |
|---|---|---|---|
| R_c | Thermal Contact Resistance | K/W | Total resistance across an interface |
| h_c | Interfacial Conductance | W/(m^2·K) | Heat transfer conductance per unit area at the interface |
| A | Interface Area | m^2 | Contact area between two surfaces |
🏭 Engineering Example
NVIDIA H100 GPU Accelerator Module
Not applicable — electronic system🏗️ Applications
- High-performance computing cooling
- Power electronics thermal management
- Spacecraft radiator joints
- Nuclear fuel cladding interfaces
- Thermoelectric generator modules
📋 Real Project Case
Air-Cooled Condenser Retrofit for 600 MW Coal Power Plant
Retrofit of legacy water-cooled condenser at Midwest US plant