Warmteoverdracht door geleiding

De fysica van warmtegeleiding

Op atomair niveau ontstaat thermische geleiding door twee mechanismen. In niet-metallische vaste stoffen wordt warmte voornamelijk getransporteerd door roostertrillingen — gekwantiseerde pakketten trillingsenergie die fononen worden genoemd. Wanneer een gebied van een kristalrooster warmer is, trillen de atomen heviger, en deze trillingen planten zich voort naar naburige atomen en dragen daarbij energie mee. In metalen komt een tweede — en meestal dominant — mechanisme in het spel: vrije elektronen. De zee van gedelokaliseerde geleidingselektronen die metalen hun elektrische geleidbaarheid geeft, fungeert ook als een uiterst efficiënte drager van thermische energie. Dit verklaart waarom metalen zowel goede elektrische als goede thermische geleiders zijn, een verband dat geformaliseerd is in de wet van Wiedemann-Franz.

De thermische geleidbaarheid $k$ van een materiaal kwantificeert hoe gemakkelijk het warmte geleidt. Koper, met $k \approx 400 \; \mathrm{W/(m \cdot K)}$, behoort tot de beste geleiders. Diamant bereikt dankzij zijn extreem stijf rooster $k \approx 2000 \; \mathrm{W/(m \cdot K)}$. Aan het andere uiterste bereikt aerogel — een materiaal op basis van silica dat grotendeels uit ingesloten luchtholtes bestaat — $k \approx 0{,}015 \; \mathrm{W/(m \cdot K)}$, waarmee het een van de beste bekende thermische isolatoren is.

De wet van Fourier voor warmtegeleiding

De hoeksteen van conductieve warmteoverdracht is de wet van Fourier, voorgesteld door Jean-Baptiste Joseph Fourier in 1822. In één dimensie stelt deze dat de warmteflux — de snelheid van warmteoverdracht per eenheid van oppervlak — evenredig is met de negatieve temperatuurgradiënt:

$${\displaystyle q = -k \, \frac{dT}{dx}}$$

Waarin:

• $q$ de warmteflux is $\left[\mathrm{W/m^2}\right]$.
• $k$ de thermische geleidbaarheid van het materiaal is $\left[\mathrm{W/(m \cdot K)}\right]$.
• $dT/dx$ de temperatuurgradiënt in de richting van de warmtestroom is $\left[\mathrm{K/m}\right]$.

Het minteken weerspiegelt de tweede hoofdwet van de thermodynamica: warmte stroomt spontaan van warm naar koud. De totale warmteoverdrachtssnelheid door een oppervlak met oppervlakte $A$ is eenvoudigweg:

$${\displaystyle \dot{Q} = q \, A = -k \, A \, \frac{dT}{dx}}$$

In drie dimensies generaliseert de wet van Fourier tot een vectorvergelijking met de temperatuurgradiënt $\nabla T$:

$${\displaystyle \vec{q} = -k \, \nabla T}$$

Deze compacte uitdrukking vormt het startpunt voor vrijwel alle analytische en numerieke behandelingen van warmtegeleiding.

De warmtevergelijking

De combinatie van de wet van Fourier met energiebehoud levert de warmtevergelijking (ook wel de diffusievergelijking genoemd) op, de fundamentele partiële differentiaalvergelijking die transiënte warmtegeleiding in een vaste stof beschrijft:

$${\displaystyle \rho \, c_p \, \frac{\partial T}{\partial t} = \nabla \cdot \left( k \, \nabla T \right) + \dot{q}_{\text{gen}}}$$

Waarin:

• $\rho$ de materiaaldichtheid is $\left[\mathrm{kg/m^3}\right]$.
• $c_p$ de soortelijke warmtecapaciteit is $\left[\mathrm{J/(kg \cdot K)}\right]$.
• $\dot{q}_{\text{gen}}$ de volumetrische interne warmtegeneratiesnelheid is $\left[\mathrm{W/m^3}\right]$, bijvoorbeeld door elektrische weerstandsverwarming of chemische reacties.

Voor een homogeen materiaal met constante thermische geleidbaarheid vereenvoudigt dit tot:

$${\displaystyle \frac{\partial T}{\partial t} = \alpha \, \nabla^2 T + \frac{\dot{q}_{\text{gen}}}{\rho \, c_p}}$$

waar $\alpha = k / (\rho \, c_p)$ de thermische diffusiviteit is $\left[\mathrm{m^2/s}\right]$ — één enkele eigenschap die aangeeft hoe snel temperatuurveranderingen zich door een materiaal voortplanten. Koper, met $\alpha \approx 1{,}17 \times 10^{-4} \; \mathrm{m^2/s}$, reageert thermisch vele malen sneller dan baksteen met $\alpha \approx 5{,}2 \times 10^{-7} \; \mathrm{m^2/s}$.

Stationaire geleiding: thermische weerstand

Wanneer de omstandigheden niet in de tijd veranderen ($\partial T / \partial t = 0$) en er geen interne warmtegeneratie is, reduceert de warmtevergelijking tot de vergelijking van Laplace: $\nabla^2 T = 0$. Voor het eenvoudige maar bijzonder nuttige geval van eendimensionale, stationaire geleiding door een vlakke wand met dikte $L$ is het temperatuurprofiel lineair en de warmteoverdrachtssnelheid:

$${\displaystyle \dot{Q} = \frac{k \, A}{L} \left(T_1 - T_2\right) = \frac{T_1 - T_2}{R_{\text{cond}}}}$$

waar de thermische weerstand voor geleiding door de wand gedefinieerd is als:

$${\displaystyle R_{\text{cond}} = \frac{L}{k \, A}}$$

Deze analogie met de wet van Ohm ($V = IR$) is buitengewoon krachtig. Net zoals elektrische weerstanden in serie en parallel geschakeld kunnen worden, kunnen thermische weerstanden aan elkaar gekoppeld worden om samengestelde wanden, isolatielagen, contactweerstanden en convectieve randvoorwaarden in één enkel geünificeerd kader te modelleren.

Samengestelde wanden

Beschouw een wand bestaande uit drie lagen (bv. pleisterwerk, baksteen, isolatie) met convectieve randvoorwaarden aan beide zijden. De totale thermische weerstand is:

$${\displaystyle R_{\text{total}} = \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}}$$

en de totale warmteoverdrachtssnelheid volgt rechtstreeks als $\dot{Q} = \Delta T_{\text{overall}} / R_{\text{total}}$. Deze aanpak vormt de ruggengraat van gebouwenergetische berekeningen en wordt dagelijks gebruikt door architecten en werktuigbouwkundig ingenieurs om isolatieprestaties en U-waarden te evalueren.

Radiale geleiding: cilinders en bollen

Veel technische componenten — leidingen, kabels, reactorvaten — hebben een cilindrische of bolvormige geometrie. Voor stationaire geleiding door een holle cilinder (binnenradius $r_1$, buitenradius $r_2$, lengte $L_{\text{cyl}}$) is het temperatuurprofiel logaritmisch en de warmteoverdrachtssnelheid:

$${\displaystyle \dot{Q} = \frac{2\pi \, k \, L_{\text{cyl}} \left(T_1 - T_2\right)}{\ln\left(r_2 / r_1\right)}}$$

De bijbehorende thermische weerstand voor radiale geleiding door een cilindrische schaal is:

$${\displaystyle R_{\text{cyl}} = \frac{\ln\left(r_2 / r_1\right)}{2\pi \, k \, L_{\text{cyl}}}}$$

Een interessant gevolg treedt op wanneer isolatie wordt aangebracht op een dunne leiding of draad: er bestaat een kritische isolatieradius $r_{\text{cr}} = k_{\text{ins}} / h$ waaronder het aanbrengen van isolatie het warmteverlies feitelijk vergroot, omdat het grotere buitenoppervlak de convectie meer bevordert dan de isolatie de geleiding vermindert. Dit contra-intuïtieve resultaat heeft praktische gevolgen voor de isolatie van elektrische kabels en leidingen met een kleine diameter.

Transiënte geleiding en het Biot-getal

Wanneer temperaturen in de tijd veranderen — een onderdeel dat wordt afgeschrikt, een oven die opwarmt, een remschijf die wrijvingswarmte absorbeert — betreden we het domein van de transiënte geleiding. De eerste vraag is of de interne temperatuur van het lichaam als ruimtelijk uniform kan worden beschouwd of dat er significante interne gradiënten ontstaan. Het Biot-getal geeft het antwoord:

$${\displaystyle \mathrm{Bi} = \frac{h \, L_c}{k_s}}$$

waar $h$ de convectieve oppervlaktecoëfficiënt is, $L_c$ een karakteristieke lengte (doorgaans volume/oppervlak), en $k_s$ de thermische geleidbaarheid van de vaste stof. Wanneer $\mathrm{Bi} < 0{,}1$, zijn interne temperatuurgradiënten verwaarloosbaar en is de lumped capacitance methode van toepassing — het gehele lichaam wordt behandeld als ware het op één enkele, tijdsafhankelijke temperatuur, beschreven door de afkoelingswet van Newton:

$${\displaystyle T(t) = T_{\infty} + \left(T_0 - T_{\infty}\right) e^{-t/\tau}}$$

met tijdsconstante $\tau = \rho \, V \, c_p / (h \, A_s)$. Wanneer $\mathrm{Bi} > 0{,}1$, worden ruimtelijke temperatuurvariaties binnen het lichaam belangrijk en moet de volledige warmtevergelijking worden opgelost — analytisch voor eenvoudige geometrieën (via scheidingsvariabelen en Heisler-diagrammen) of numeriek met de eindige-elementenmethode voor complexe vormen.

Engineering toepassingen

1. Gebouwisolatie en energie-efficiëntie

Geleiding door wanden, daken en vloeren is de primaire weg voor warmteverlies in gebouwen. Ingenieurs gebruiken de thermische weerstandsnetwerkmethode om de totale U-waarde van constructie-elementen te evalueren en te bepalen waar aanvullende isolatie het beste rendement oplevert. Materialen zoals geëxpandeerd polystyreen ($k \approx 0{,}035 \; \mathrm{W/(m \cdot K)}$), minerale wol en vacuümisolatiepanelen worden specifiek geselecteerd vanwege hun lage thermische geleidbaarheid.

2. Thermisch beheer van elektronica

In de micro-elektronica moet de warmte die door transistors wordt gegenereerd, worden geleid door de silicium-chip, door thermische interfacematerialen (TIM's), en naar de warmtespreider en het koellichaam voordat ze door convectie kan worden afgevoerd. Elke laag vertegenwoordigt een thermische weerstand, en de totale junction-naar-omgeving weerstand bepaalt de bedrijfstemperatuur van de chip. Het verlagen van elke schakel in deze keten — bijvoorbeeld door warmtespreiders van hooggeleideind koper of diamant te gebruiken — verlaagt rechtstreeks de junctietemperatuur en verbetert de betrouwbaarheid.

3. Metallurgische warmtebehandeling

Processen als gloeien, afschrikken en ontlaten zijn sterk afhankelijk van de beheersing van de snelheid en uniformiteit van temperatuurverandering in een metalen onderdeel. Interne geleiding bepaalt hoe snel de kern van een dik stalen smeedstuk de doeltemperatuur bereikt in vergelijking met het oppervlak. Foutieve beheersing resulteert in residuele spanningen, vervorming of ongewenste microstructuren. Voorspellende thermische analyse met FEA wordt op grote schaal ingezet om warmtebehandelingscycli te ontwerpen die de gewenste eigenschappen bereiken en tegelijkertijd defecten minimaliseren.

4. Thermische beschermingssystemen

Ruimtevaartuigen die de atmosfeer binnenkomen, ondervinden extreme verhitting — oppervlaktetemperaturen kunnen 1500 °C overschrijden. Ablatieve hitteschilden en keramische tegelsystemen benutten materialen met een lage thermische geleidbaarheid om de structuur van het voertuig koel te houden, waardoor de warmtepuls de tijd krijgt om te passeren voordat deze kritieke componenten bereikt. Het ontwerp van deze systemen omvat het oplossen van transiënte geleidingsproblemen met temperatuurafhankelijke materiaaleigenschappen en oppervlakte-ablatie.

Casestudy: warmteverlies door een geïsoleerde leiding

Een stalen leiding ($k_{\text{steel}} = 50 \; \mathrm{W/(m \cdot K)}$) transporteert stoom op 200 °C. De leiding heeft een binnenradius $r_1 = 25 \; \mathrm{mm}$, buitenradius $r_2 = 30 \; \mathrm{mm}$, en is bedekt met 40 mm minerale wolisolatie ($k_{\text{ins}} = 0{,}04 \; \mathrm{W/(m \cdot K)}$, buitenradius $r_3 = 70 \; \mathrm{mm}$). De omgevingstemperatuur is 20 °C met een convectieve buitencoëfficiënt $h_o = 10 \; \mathrm{W/(m^2 \cdot K)}$. Per meter leidinglengte ($L = 1 \; \mathrm{m}$):

$${\displaystyle R_{\text{steel}} = \frac{\ln(r_2/r_1)}{2\pi \, k_{\text{steel}} \, L} = \frac{\ln(30/25)}{2\pi \times 50 \times 1} = \frac{0.1823}{314.16} \approx 5.8 \times 10^{-4} \; \mathrm{K/W}}$$

$${\displaystyle R_{\text{ins}} = \frac{\ln(r_3/r_2)}{2\pi \, k_{\text{ins}} \, L} = \frac{\ln(70/30)}{2\pi \times 0.04 \times 1} = \frac{0.8473}{0.2513} \approx 3.371 \; \mathrm{K/W}}$$

$${\displaystyle R_{\text{conv}} = \frac{1}{h_o \, 2\pi \, r_3 \, L} = \frac{1}{10 \times 2\pi \times 0.07 \times 1} = \frac{1}{4.398} \approx 0.227 \; \mathrm{K/W}}$$

De totale weerstand is $R_{\text{total}} = 5{,}8 \times 10^{-4} + 3{,}371 + 0{,}227 \approx 3{,}60 \; \mathrm{K/W}$, en het warmteverlies per meter bedraagt:

$${\displaystyle \dot{Q} = \frac{200 - 20}{3.60} \approx 50 \; \mathrm{W/m}}$$

Merk op dat de stalen wand vrijwel niets bijdraagt aan de totale weerstand ($0{,}016\%$), terwijl de isolatielaag $93{,}6\%$ van de thermische weerstand levert. Zonder isolatie zou het warmteverlies stijgen tot ruwweg 800 W/m — een zestienvoudige toename. Dit type thermische weerstandsanalyse is routine in het ontwerp van industriële leidingsystemen, procestechniek en energie-audits.

Numerieke methoden voor complexe geleidingsproblemen

Analytische oplossingen bestaan enkel voor relatief eenvoudige geometrieën en randvoorwaarden. In de praktijk hebben technische componenten onregelmatige vormen, temperatuurafhankelijke materiaaleigenschappen, interne warmtebronnen, contactweerstanden aan interfaces en complexe samenstellingen van meerdere materialen. Voor deze problemen is de eindige-elementenmethode (FEM) het standaard numerieke gereedschap. FEM discretiseert de geometrie in een mesh van elementen, benadert het temperatuurveld met stuksgewijze polynomiale functies, en lost het resulterende stelsel algebraïsche vergelijkingen op om een gedetailleerde temperatuurkaart van het volledige onderdeel te produceren.

Moderne FEA-software kan stationaire en transiënte geleiding behandelen, gecombineerd met convectieve en straling randvoorwaarden, niet-lineaire materiaaleigenschappen, faseovergangen, en zelfs thermo-mechanische koppeling waarbij thermische spanningen en vervormingen terugkoppelen naar het thermische probleem. Wanneer de vloeistofstroming zelf deel uitmaakt van het probleem — gedwongen of natuurlijke convectie, geconjugeerde warmteoverdracht door een warmtewisselaar — neemt Computational Fluid Dynamics (CFD) het over en lost de stromings- en energievergelijkingen gelijktijdig op. Samen vormen FEA en CFD de ruggengraat van moderne warmteoverdrachtsanalyse, waarmee ingenieurs temperaturen, warmtestromen en thermische spanningen in componenten van elke complexiteit kunnen voorspellen.

Conclusie

Warmtegeleiding is de meest intuïtieve en wiskundig hanteerbare van de drie warmteoverdrachtsmodi, maar het technische belang ervan is enorm. De wet van Fourier en het concept van thermische weerstand bieden elegante analytische hulpmiddelen voor snelle schattingen en ontwerprichtlijnen, terwijl de warmtevergelijking — numeriek opgelost door de eindige-elementenmethode — de volledige complexiteit van problemen uit de praktijk aankan. Of het doel nu is om warmteverlies door een gebouwschil te minimaliseren, temperaturen in een microchip te beheersen, de afkoelsnelheid tijdens een afschrikproces te regelen, of een thermisch beschermingssysteem voor terugkeer in de atmosfeer te ontwerpen: geleiding maakt steeds deel uit van het verhaal. Voor technische uitdagingen waar nauwkeurigheid ertoe doet en de belangen hoog zijn, bieden professionele thermische analysediensten de grondigheid en het inzicht die nodig zijn om thermisch begrip om te zetten in betrouwbare, geoptimaliseerde ontwerpen.