Wärmeübertragung durch Wärmeleitung

Die Physik der Wärmeleitung

Auf atomarer Ebene entsteht Wärmeleitung durch zwei Mechanismen. In nichtmetallischen Festkörpern wird Wärme hauptsächlich durch Gitterschwingungen transportiert — quantisierte Pakete von Schwingungsenergie, genannt Phononen. Wenn ein Bereich eines Kristallgitters heißer ist, schwingen seine Atome stärker, und diese Schwingungen breiten sich zu Nachbaratomen aus und tragen Energie mit sich. In Metallen kommt ein zweiter — und meist dominierender — Mechanismus ins Spiel: freie Elektronen. Das Meer delokalisierter Leitungselektronen, das Metallen ihre elektrische Leitfähigkeit verleiht, dient auch als äußerst effizienter Träger thermischer Energie. Deshalb sind Metalle sowohl gute elektrische als auch gute thermische Leiter, ein Zusammenhang, der durch das Wiedemann-Franz-Gesetz formalisiert wird.

Die Wärmeleitfähigkeit $k$ eines Werkstoffs quantifiziert, wie leicht er Wärme leitet. Kupfer mit $k \approx 400 \; \mathrm{W/(m \cdot K)}$ gehört zu den besten Leitern. Diamant erreicht dank seines extrem starren Gitters $k \approx 2000 \; \mathrm{W/(m \cdot K)}$. Am anderen Ende des Spektrums erreicht Aerogel — ein siliziumbasiertes Material, das hauptsächlich aus eingeschlossenen Lufttaschen besteht — $k \approx 0{,}015 \; \mathrm{W/(m \cdot K)}$ und ist damit einer der besten bekannten Wärmeisolatoren.

Fouriersches Gesetz der Wärmeleitung

Der Eckpfeiler der konduktiven Wärmeübertragung ist das Fouriersche Gesetz, das 1822 von Jean-Baptiste Joseph Fourier aufgestellt wurde. In einer Dimension besagt es, dass die Wärmestromdichte — die Wärmeübertragungsrate pro Flächeneinheit — proportional zum negativen Temperaturgradienten ist:

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

Dabei gilt:

• $q$ ist die Wärmestromdichte $\left[\mathrm{W/m^2}\right]$.
• $k$ ist die Wärmeleitfähigkeit des Materials $\left[\mathrm{W/(m \cdot K)}\right]$.
• $dT/dx$ ist der Temperaturgradient in Richtung des Wärmeflusses $\left[\mathrm{K/m}\right]$.

Das negative Vorzeichen spiegelt den zweiten Hauptsatz der Thermodynamik wider: Wärme fließt spontan von heiß nach kalt. Die gesamte Wärmeübertragungsrate durch eine Fläche $A$ ist einfach:

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

In drei Dimensionen verallgemeinert sich das Fouriersche Gesetz zu einer Vektorgleichung mit dem Temperaturgradienten $\nabla T$:

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

Dieser kompakte Ausdruck ist der Ausgangspunkt für praktisch alle analytischen und numerischen Behandlungen der Wärmeleitung.

Die Wärmeleitungsgleichung

Die Kombination des Fourierschen Gesetzes mit der Energieerhaltung ergibt die Wärmeleitungsgleichung (auch Diffusionsgleichung), die fundamentale partielle Differentialgleichung der instationären Wärmeleitung in einem Festkörper:

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

Dabei gilt:

• $\rho$ ist die Materialdichte $\left[\mathrm{kg/m^3}\right]$.
• $c_p$ ist die spezifische Wärmekapazität $\left[\mathrm{J/(kg \cdot K)}\right]$.
• $\dot{q}_{\text{gen}}$ ist die volumetrische interne Wärmeerzeugungsrate $\left[\mathrm{W/m^3}\right]$, z. B. durch elektrische Widerstandserwärmung oder chemische Reaktionen.

Für ein homogenes Material mit konstanter Wärmeleitfähigkeit vereinfacht sich dies zu:

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

wobei $\alpha = k / (\rho \, c_p)$ die Temperaturleitfähigkeit $\left[\mathrm{m^2/s}\right]$ ist — eine einzelne Eigenschaft, die angibt, wie schnell sich Temperaturänderungen durch ein Material ausbreiten. Kupfer mit $\alpha \approx 1{,}17 \times 10^{-4} \; \mathrm{m^2/s}$ reagiert thermisch weit schneller als Ziegel mit $\alpha \approx 5{,}2 \times 10^{-7} \; \mathrm{m^2/s}$.

Stationäre Wärmeleitung: Wärmewiderstand

Wenn sich die Bedingungen zeitlich nicht ändern ($\partial T / \partial t = 0$) und keine interne Wärmeerzeugung vorliegt, reduziert sich die Wärmeleitungsgleichung auf die Laplace-Gleichung: $\nabla^2 T = 0$. Für den einfachen, aber äußerst nützlichen Fall eindimensionaler, stationärer Wärmeleitung durch eine ebene Wand der Dicke $L$ ist das Temperaturprofil linear und die Wärmeübertragungsrate beträgt:

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

wobei der Wärmewiderstand für die Wärmeleitung durch die Wand definiert ist als:

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

Diese Analogie zum Ohmschen Gesetz ($V = IR$) ist außerordentlich leistungsfähig. Genau wie elektrische Widerstände in Reihen- und Parallelschaltung kombiniert werden können, lassen sich Wärmewiderstände aneinanderketten, um zusammengesetzte Wände, Dämmschichten, Kontaktwiderstände und konvektive Randbedingungen in einem einheitlichen Rahmen zu modellieren.

Mehrschichtige Wände

Betrachten wir eine Wand aus drei Schichten (z. B. Putz, Ziegel, Dämmung) mit konvektiven Randbedingungen auf beiden Seiten. Der gesamte Wärmewiderstand beträgt:

$${\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}}$$

und die Gesamt-Wärmeübertragungsrate ergibt sich direkt als $\dot{Q} = \Delta T_{\text{gesamt}} / R_{\text{total}}$. Dieser Ansatz ist das Rückgrat der Gebäudeenergieberechnungen und wird täglich von Architekten und Maschinenbauingenieuren zur Bewertung der Dämmleistung und U-Werte eingesetzt.

Radiale Wärmeleitung: Zylinder und Kugeln

Viele Ingenieurbauteile — Rohre, Kabel, Reaktorbehälter — haben zylindrische oder sphärische Geometrie. Für stationäre Wärmeleitung durch einen Hohlzylinder (Innenradius $r_1$, Außenradius $r_2$, Länge $L_{\text{cyl}}$) ist das Temperaturprofil logarithmisch und die Wärmeübertragungsrate beträgt:

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

Der entsprechende Wärmewiderstand für radiale Wärmeleitung durch eine zylindrische Schale ist:

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

Eine interessante Konsequenz ergibt sich beim Aufbringen von Isolierung auf ein dünnes Rohr oder einen Draht: Es existiert ein kritischer Isolationsradius $r_{\text{cr}} = k_{\text{ins}} / h$, unterhalb dessen das Hinzufügen von Isolierung den Wärmeverlust tatsächlich erhöht, da die größere Außenfläche die Konvektion stärker fördert, als die Isolierung die Wärmeleitung reduziert. Dieses kontraintuitive Ergebnis hat praktische Auswirkungen auf die Isolierung von Elektrokabeln und Rohrleitungen mit kleinem Durchmesser.

Instationäre Wärmeleitung und die Biot-Zahl

Wenn sich Temperaturen mit der Zeit ändern — ein Bauteil wird abgeschreckt, ein Ofen heizt vor, eine Bremsscheibe nimmt Reibungswärme auf — betreten wir den Bereich der instationären Wärmeleitung. Die erste Frage ist, ob die Innentemperatur des Körpers als räumlich gleichmäßig behandelt werden kann oder ob sich signifikante interne Gradienten entwickeln. Die Biot-Zahl liefert die Antwort:

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

wobei $h$ der konvektive Oberflächenkoeffizient, $L_c$ eine charakteristische Länge (typischerweise Volumen/Oberfläche) und $k_s$ die Wärmeleitfähigkeit des Festkörpers ist. Wenn $\mathrm{Bi} < 0{,}1$, sind interne Temperaturgradienten vernachlässigbar und das Lumped-Capacitance-Verfahren ist anwendbar — der gesamte Körper wird als bei einer einzigen, zeitabhängigen Temperatur befindlich behandelt, beschrieben durch das Newtonsche Abkühlungsgesetz:

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

mit der Zeitkonstante $\tau = \rho \, V \, c_p / (h \, A_s)$. Wenn $\mathrm{Bi} > 0{,}1$, werden räumliche Temperaturvariationen innerhalb des Körpers bedeutsam und die vollständige Wärmeleitungsgleichung muss gelöst werden — analytisch für einfache Geometrien (mittels Variablentrennung und Heisler-Diagrammen) oder numerisch mit der Finite-Elemente-Methode für komplexe Formen.

Anwendungen im Ingenieurwesen

1. Gebäudedämmung und Energieeffizienz

Wärmeleitung durch Wände, Dächer und Böden ist der Hauptpfad für Wärmeverluste in Gebäuden. Ingenieure verwenden den Wärmewiderstandsnetzwerk-Ansatz, um den Gesamt-U-Wert von Baukonstruktionen zu bewerten und zu bestimmen, wo zusätzliche Dämmung die beste Rendite bringt. Materialien wie expandiertes Polystyrol ($k \approx 0{,}035 \; \mathrm{W/(m \cdot K)}$), Mineralwolle und Vakuumdämmplatten werden gezielt wegen ihrer niedrigen Wärmeleitfähigkeit ausgewählt.

2. Thermisches Management in der Elektronik

In der Mikroelektronik muss die von Transistoren erzeugte Wärme durch den Silizium-Chip, durch thermische Schnittstellenmaterialien (TIMs) und in den Wärmespreizer und Kühlkörper geleitet werden, bevor sie durch Konvektion abgeführt werden kann. Jede Schicht stellt einen Wärmewiderstand dar, und der Gesamtwiderstand vom Sperrschicht bis zur Umgebung bestimmt die Betriebstemperatur des Chips. Die Reduzierung jedes Glieds in dieser Kette — zum Beispiel durch Verwendung hochleitfähiger Kupfer- oder Diamant-Wärmespreizer — senkt direkt die Sperrschichttemperatur und verbessert die Zuverlässigkeit.

3. Metallurgische Wärmebehandlung

Prozesse wie Glühen, Abschrecken und Anlassen hängen entscheidend von der Steuerung der Geschwindigkeit und Gleichmäßigkeit der Temperaturänderung innerhalb eines Metallbauteils ab. Die interne Wärmeleitung bestimmt, wie schnell der Kern eines dicken Stahlschmiedestücks die Zieltemperatur erreicht im Vergleich zu seiner Oberfläche. Fehler hierbei führen zu Eigenspannungen, Verzug oder unerwünschten Gefügestrukturen. Prädiktive Thermalanalyse mittels FEA wird breit eingesetzt, um Wärmebehandlungszyklen zu entwerfen, die die gewünschten Eigenschaften bei minimalen Defekten erreichen.

4. Thermische Schutzsysteme

Raumfahrzeuge, die in die Atmosphäre eintreten, erfahren extreme Aufheizung — Oberflächentemperaturen können 1500 °C überschreiten. Ablative Hitzeschilde und keramische Kachelsysteme nutzen Materialien mit niedriger Wärmeleitfähigkeit, um die Struktur des Fahrzeugs kühl zu halten und Zeit zu gewinnen, bis der Wärmeimpuls vorübergeht, bevor er zu kritischen Komponenten vordringt. Die Auslegung dieser Systeme umfasst das Lösen instationärer Wärmeleitungsprobleme mit temperaturabhängigen Materialeigenschaften und Oberflächenablation.

Fallstudie: Wärmeverlust durch ein isoliertes Rohr

Ein Stahlrohr ($k_{\text{steel}} = 50 \; \mathrm{W/(m \cdot K)}$) führt Dampf bei 200 °C. Das Rohr hat einen Innenradius $r_1 = 25 \; \mathrm{mm}$, Außenradius $r_2 = 30 \; \mathrm{mm}$ und ist mit 40 mm Mineralwolledämmung umhüllt ($k_{\text{ins}} = 0{,}04 \; \mathrm{W/(m \cdot K)}$, Außenradius $r_3 = 70 \; \mathrm{mm}$). Die Umgebungstemperatur beträgt 20 °C mit einem äußeren konvektiven Koeffizienten $h_o = 10 \; \mathrm{W/(m^2 \cdot K)}$. Pro Meter Rohrlänge ($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}}$$

Der Gesamtwiderstand beträgt $R_{\text{total}} = 5{,}8 \times 10^{-4} + 3{,}371 + 0{,}227 \approx 3{,}60 \; \mathrm{K/W}$, und der Wärmeverlust pro Meter ist:

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

Beachten Sie, dass die Stahlwand nahezu nichts zum Gesamtwiderstand beiträgt ($0{,}016\%$), während die Dämmschicht $93{,}6\%$ des Wärmewiderstands liefert. Ohne Dämmung würde der Wärmeverlust auf etwa 800 W/m springen — eine sechzehnfache Steigerung. Diese Art der Wärmewiderstandsanalyse ist Routine in der industriellen Rohrleitungsplanung, Verfahrenstechnik und Energieberatung.

Numerische Methoden für komplexe Wärmeleitungsprobleme

Analytische Lösungen existieren nur für relativ einfache Geometrien und Randbedingungen. In der Praxis haben Ingenieurbauteile unregelmäßige Formen, temperaturabhängige Materialeigenschaften, interne Wärmequellen, Kontaktwiderstände an Grenzflächen und komplexe Mehrwerkstoff-Baugruppen. Für diese Probleme ist die Finite-Elemente-Methode (FEM) das numerische Standardwerkzeug. FEM diskretisiert die Geometrie in ein Netz von Elementen, approximiert das Temperaturfeld mit stückweisen Polynomfunktionen und löst das resultierende System algebraischer Gleichungen, um eine detaillierte Temperaturkarte des gesamten Bauteils zu erzeugen.

Moderne FEA-Software kann stationäre und instationäre Wärmeleitung handhaben, gekoppelt mit konvektiven und Strahlungs-Randbedingungen, nichtlinearen Materialeigenschaften, Phasenänderungen und sogar thermo-mechanischer Kopplung, bei der thermische Spannungen und Verformungen auf das thermische Problem zurückwirken. Wenn die Fluidströmung selbst Teil des Problems ist — erzwungene oder natürliche Konvektion, konjugierter Wärmeübergang durch einen Wärmetauscher — übernimmt Computational Fluid Dynamics (CFD), das Strömungs- und Energiegleichungen gleichzeitig löst. Zusammen bilden FEA und CFD das Rückgrat der modernen Wärmeübertragungsanalyse und ermöglichen es Ingenieuren, Temperaturen, Wärmeströme und Thermospannungen in Bauteilen beliebiger Komplexität vorherzusagen.

Fazit

Wärmeleitung ist der intuitivste und mathematisch handhabbarste der drei Wärmeübertragungsmechanismen, doch seine ingenieurstechnische Bedeutung ist immens. Das Fouriersche Gesetz und das Wärmewiderstandskonzept bieten elegante analytische Werkzeuge für schnelle Abschätzungen und Auslegungsrichtlinien, während die Wärmeleitungsgleichung — numerisch mit der Finite-Elemente-Methode gelöst — die volle Komplexität realer Probleme bewältigt. Ob das Ziel die Minimierung von Wärmeverlusten durch eine Gebäudehülle, das Temperaturmanagement in einem Mikrochip, die Steuerung der Abkühlrate beim Abschrecken oder die Auslegung eines Wärmeschutzsystems für den Wiedereintritt ist — Wärmeleitung ist stets Teil der Geschichte. Für Ingenieurherausforderungen, bei denen Genauigkeit zählt und viel auf dem Spiel steht, liefern professionelle Thermalanalyse-Dienstleistungen die Rigorosität und Einsicht, die nötig sind, um thermisches Verständnis in zuverlässige, optimierte Designs umzusetzen.