Stationäre Wärmeleitung berechnen – Modul EA

Das Modul berechnet die stationäre Wärmeleitung nach dem VDI-Wärmeatlas (12. Auflage 2019) für die klassischen Grundgeometrien: ebene Wände und Platten, Hohlzylinder (Rohre) und Hohlkugeln, jeweils auch mehrschichtig, sowie Rippen, Stäbe und Drähte.

Modul EARegelwerk VDI-Wärmeatlas, 12. Auflage 2019Lesedauer 7 minDE / EN

Aufgabenstellung und Berechnungsziel

Das Modul berechnet die stationäre Wärmeleitung nach dem VDI-Wärmeatlas (12. Auflage 2019) für die klassischen Grundgeometrien: ebene Wände und Platten, Hohlzylinder (Rohre) und Hohlkugeln, jeweils auch mehrschichtig, sowie Rippen, Stäbe und Drähte. Zusätzlich lassen sich zweidimensionale Leitungsprobleme über Formkoeffizienten behandeln, etwa erdverlegte Rohre oder Körper mit ovalen und rechteckigen Querschnitten.

Stationäre Wärmeleitung berechnen ist eine der Grundaufgaben des Apparate- und Anlagenbaus: Temperaturverläufe in Behälter- und Ofenwänden, Wärmeströme durch Rohrwände und Dämmschichten, die Wirksamkeit von Rippen auf Heiz- und Kühlflächen oder die Wandtemperaturen mehrschichtiger Aufbauten. Aus Wärmeleitfähigkeit, Schichtdicke und Fläche ergibt sich der Wärmeleitwiderstand jeder Schicht; die Reihenschaltung der Widerstände liefert Wärmestrom und Zwischentemperaturen.

Das Modul folgt dabei den Beziehungen aus Abschnitt E1 des VDI-Wärmeatlas und stellt die Ergebnisse so bereit, dass sie direkt in Wärmedurchgangs- und Verlustrechnungen (etwa mit dem Modul zum Wärmeverlust von Wänden und Rohrleitungen) weiterverwendet werden können.

Regelwerk und Berechnungsbasis: VDI-Wärmeatlas, 12. Auflage 2019

So läuft die Berechnung ab

  1. Geometrie festlegen: Zuerst wird die Grundgeometrie gewählt: ebene Wand, Hohlzylinder oder Hohlkugel, ein- oder mehrschichtig, alternativ Rippe, Stab oder ein zweidimensionaler Fall mit Formkoeffizient. Die Geometrie bestimmt, welche Leitungsgleichung angesetzt wird.
  2. Abmessungen und Stoffwerte vorgeben: Für jede Schicht werden Dicke bzw. Innen- und Außendurchmesser sowie die Wärmeleitfähigkeit bei der maßgebenden Mitteltemperatur eingegeben. Bei temperaturabhängiger Leitfähigkeit wird mit dem Wert bei der mittleren Schichttemperatur gerechnet, was ggf. eine kurze Iteration erfordert.
  3. Randbedingungen setzen: Als thermische Randbedingungen dienen die Oberflächen- bzw. Wandtemperaturen auf beiden Seiten oder ein vorgegebener Wärmestrom. Damit ist das stationäre Leitungsproblem vollständig bestimmt.
  4. Widerstände und Wärmestrom berechnen: Für jede Schicht wird der Wärmeleitwiderstand gebildet – bei der ebenen Wand proportional zur Dicke, beim Hohlzylinder mit dem Logarithmus des Durchmesserverhältnisses, bei der Hohlkugel mit der Differenz der Kehrwerte der Radien. Die Reihenschaltung liefert den Gesamtwiderstand und daraus den Wärmestrom.
  5. Temperaturverlauf und Sonderfälle auswerten: Aus Wärmestrom und Einzelwiderständen folgen die Temperaturen an den Schichtgrenzen. Bei Rippen und Stäben werden zusätzlich Rippenwirkungsgrad und Temperaturverlauf längs der Rippe berechnet; bei zweidimensionalen Fällen wird der Wärmestrom über den Formkoeffizienten der Anordnung bestimmt.
Eingabegrößen24 / 34 Größen
GrößeSymbolEinheit
FlächeA
Wandtemperatur 1ϑ0°C
Wandtemperatur 2ϑn°C
Anzahl Schichtenn
Schicht 1s1 λ1m
Schicht 2s2 λ2m
Schicht 3s3 λ3m
Schicht 4s4 λ4m
Schicht 1r1 λ1W/(m·K)
Schicht 2r2 λ2W/(m·K)
Schicht 3r3 λ3W/(m·K)
Schicht 4s4 λ4W/(m·K)
Schicht 1ϑ1°C
Schicht 2ϑ2°C
Schicht 3ϑ3°C
Schicht 4ϑ4°C
X-Koordinatexm
Beginn der Schicht Alle Koordinaten beziehen sich auf x = 0xkm
WärmeleitfähigkeitλkW/(m·K)
Temperatur am Beginnϑk°C
Längelm
Innenr0m
Schicht 1r1 λ1m
Schicht 2r2 λ2m
Berechnungsergebnisse24 / 28 Größen
GrößeSymbolEinheit
Wärmestrom je LängeneinheitQ/l = QlW/m
WärmeleitfähigkeitλW/(m·K)
FormfaktorSl-
Temperaturϑ2°C
Temperaturϑ1°C
Berechnung nach VDI-Wärmeatlas Tabelle Nr.Tafel-
Anordnung Nr.Anordnung-
Parameter b bzw b1bm
Parameter Deltaδm
Parameter a bzw. a1a cm
Parameter ca cm
Parameter nn-
Parameter hhm
Parameter eem
Parameter ddm
Parameter b2b1 b2m
Parameter a2a1 a2m
FormfaktorSr-
Parameter rrm
Parameter r1r1 r2m
Parameter r2r1 r2m
Kreisscheibe / KugelKreisscheibe-
_Bemerkung_Bemerkung-
_Bemerkung_Bemerkung-

Berechnungsoptionen

Kreisscheibe / Kugel

Kreisscheibe · Kugel

Gelöstes Beispiel

Eine Rohrleitung ist mit einer 50 mm dicken Mineralwolleschale gedämmt. Die Dämmung beginnt am Außendurchmesser des Rohres von 88,9 mm; ihre Oberflächentemperaturen sind innen 150 °C und außen 40 °C. Die Wärmeleitfähigkeit der Mineralwolle beträgt bei der Mitteltemperatur 0,05 W/(m·K). Gesucht ist der längenbezogene Wärmestrom durch die Dämmschicht.

Gegebene Werte

Innendurchmesser der Dämmschicht d188,9 mm
Außendurchmesser der Dämmschicht d2188,9 mm
Innere Oberflächentemperatur ϑ1150 °C
Äußere Oberflächentemperatur ϑ240 °C
Wärmeleitfähigkeit λ0,05 W/(m·K)

Lösungsweg

1

Formel für den Hohlzylinder

Für die stationäre Wärmeleitung durch eine zylindrische Schicht gilt der längenbezogene Wärmestrom:

ql = 2 · π · λ · (ϑ1 − ϑ2) / ln(d2/d1)

2

Durchmesserverhältnis einsetzen

ln(d2/d1) = ln(188,9 / 88,9) = ln(2,125) = 0,7537

3

Wärmestrom berechnen

ql = 2 · π · 0,05 W/(m·K) · (150 − 40) K / 0,7537

ql = 34,56 / 0,7537 W/m ≈ 45,9 W/m

Je Meter Rohrleitung strömen also rund 46 W durch die Dämmschicht. Der Temperaturverlauf in der Schicht ist dabei logarithmisch über dem Radius.

Ergebnis

Längenbezogener Wärmestrom q_l≈ 45,9 W/m

Alle Werte sind beispielhaft. Maßgeblich sind das jeweils gültige Regelwerk und die projektspezifischen Randbedingungen.

Häufige Fragen

Wann darf ich stationär rechnen und wann brauche ich eine instationäre Betrachtung?

Stationäre Rechnung setzt voraus, dass sich die Randbedingungen zeitlich nicht oder nur langsam gegenüber der thermischen Zeitkonstante des Bauteils ändern. Für Dauerbetriebszustände von Öfen, Rohrleitungen und Behältern ist das erfüllt. An- und Abfahrvorgänge, Chargenprozesse oder Störfallszenarien erfordern dagegen eine instationäre Rechnung, weil dort die Wärmespeicherung in der Wand das Verhalten dominiert.

Warum wird beim Rohr mit dem Logarithmus des Durchmesserverhältnisses gerechnet?

Beim Hohlzylinder wächst die durchströmte Fläche mit dem Radius, die Wärmestromdichte nimmt also nach außen ab. Die Integration der Fourierschen Gleichung über den Radius liefert deshalb den Widerstand R = ln(d_a/d_i)/(2·π·λ·L) je Rohrlänge. Die Näherung über eine mittlere Fläche ist nur bei dünnwandigen Rohren zulässig; bei dickwandigen Schichten wie Dämmungen führt sie zu deutlichen Fehlern.

Wie gehe ich mit Kontaktwiderständen zwischen den Schichten um?

Die ideale Reihenschaltung setzt vollflächigen Kontakt der Schichten voraus. Reale Grenzflächen – etwa zwischen Rohr und aufgeschobener Dämmschale oder zwischen Wand und Bekleidung – haben Luftspalte, die als zusätzlicher Widerstand wirken. Bei Dämmberechnungen wirkt das meist günstig (mehr Widerstand), bei Kühlaufgaben aber gefährlich, weil die Wärmeabfuhr überschätzt wird. Im Zweifel sind Kontaktwiderstände konservativ anzusetzen.

Was sagt der Rippenwirkungsgrad aus?

Der Rippenwirkungsgrad ist das Verhältnis des tatsächlich über die Rippe abgeführten Wärmestroms zu dem Wärmestrom, den die Rippe bei durchgehend Fußtemperatur abführen würde. Er sinkt mit zunehmender Rippenhöhe, steigendem Wärmeübergangskoeffizienten und abnehmender Wärmeleitfähigkeit des Rippenmaterials. Lange, dünne Rippen aus schlecht leitendem Material bringen deshalb kaum Zusatzleistung – die Rippenauslegung ist immer ein Optimum aus Fläche und Wirkungsgrad.

Verwandte Berechnungen