StartseiteThermodynamikWärmegleichung — 2D-Temperaturverteilung

🌡 Wärmegleichung — 2D-Temperaturverteilung

Simulieren Sie die 2D-Wärmeleitung mit der Finite-Differenzen-Methode. Zeichnen Sie Wärmequellen und -senken, beobachten Sie, wie die Temperatur diffundiert, und vergleichen Sie stationäre mit instationären Lösungen.

Thermodynamik3DMittel60 FPS
heat-equation ↗ Eigenständig öffnen

Über diese Simulation

Diese Simulation löst die zweidimensionale Wärmegleichung ∂T/∂t = α∇²T auf einem 80×80-Gitter mit einer expliziten Finite-Differenzen-Methode. Malen Sie heiße Quellen, kalte Senken und isolierende Wände direkt auf das Canvas und beobachten Sie, wie sich die Temperatur nach außen ausbreitet, während der Löser wiederholt den fünf-Punkte-Laplace-Stern anwendet. Eine separate Gauß-Seidel-Relaxationsroutine lässt Sie direkt zur stationären Verteilung springen — dem harmonischen Feld, das durch die Laplace-Gleichung ∇²T = 0 beschrieben wird, sobald alle Einschwingvorgänge abgeklungen sind.

🔬 Was es zeigt

Jedes Animationsbild berechnet die Temperatur jeder freien Zelle aus ihren vier Nachbarn neu, mit T ← T + (α·Δt/Δx²)·(Tlinks+Trechts+Toben+Tunten−4T), farbcodiert von tiefem Blau (0 °C) über Grün und Gelb bis Weiß (100 °C). Heiße Zellen bleiben bei 100 °C fixiert, kalte Zellen bei 0 °C, und Wände liegen bei fest 20 °C, während sie die Diffusion über sich hinweg blockieren. Die Anzeigen für Max, Min und Mittelwert sowie der Schrittzähler werden jedes Bild aktualisiert.

🎮 Anwendung

Wählen Sie Heiß, Kalt oder Wand und klicken-ziehen Sie dann auf dem Canvas, um Quellen in der gewählten Pinselgröße zu malen; Löschen setzt Zellen zurück zur freien Diffusion. Der Diffusivität-α-Schieberegler (0,01–0,48) legt fest, wie schnell sich Wärme ausbreitet, während der Code den effektiven Zeitschritt so begrenzt, dass α·Δt/Δx² nie etwa 0,24 überschreitet. Stationärer Zustand führt sofort 3000 Gauß-Seidel-Durchläufe aus; Löschen setzt die Platte auf einheitliche 20 °C zurück; die Schaltflächen Stab und Raum laden fertige Voreinstellungen.

💡 Wussten Sie schon?

Joseph Fourier führte diese Gleichung 1822 in seiner Théorie analytique de la chaleur ein. Dieselbe Mathematik, mit Temperatur umbenannt in Wahrscheinlichkeitsdichte, beschreibt die Brownsche Bewegung und liegt der in der modernen Finanzwelt verbreiteten Black-Scholes-Optionspreisformel zugrunde.

Häufig gestellte Fragen

Was steuert der Diffusivität-α-Schieberegler?

Alpha legt die thermische Diffusivität in ∂T/∂t = α∇²T fest und steuert, wie schnell sich Wärme in jedem Schritt durch die Platte ausbreitet. Die Simulation skaliert ihren internen Zeitschritt automatisch so, dass α·Δt/Δx² nie etwa 0,24 überschreitet, wodurch das explizite Finite-Differenzen-Update selbst beim Maximalwert des Reglers von 0,48 numerisch stabil bleibt.

Wie unterscheiden sich Wände von heißen und kalten Quellen?

Heiße und kalte Zellen sind auf feste Temperaturen (100 °C und 0 °C) fixiert und wirken als kontinuierliche Quellen oder Senken, die in jedem Schritt weiter Wärme zuführen oder abführen. Wände werden bei konstant 20 °C gehalten, aber wichtiger ist, dass das Rechenschema eine benachbarte Wandzelle als gleich der eigenen Temperatur der aktuellen Zelle behandelt, sodass sie den Wärmefluss über die Grenze hinweg effektiv wie ein Isolator blockiert.

Was passiert, wenn ich auf Stationärer Zustand klicke?

Die Schaltfläche Stationärer Zustand stoppt die animierte Zeitschrittberechnung und führt stattdessen 3000 Iterationen der Gauß-Seidel-Relaxation aus, wobei jede freie Zelle wiederholt mit ihren vier Nachbarn gemittelt wird. Dies konvergiert gegen die Lösung der Laplace-Gleichung ∇²T = 0 — die Verteilung, in die sich die Platte schließlich einpendeln würde, wenn man sie unendlich lange laufen ließe — fast augenblicklich berechnet statt Bild für Bild.

Was zeigen die Voreinstellungen Stab und Raum?

Voreinstellung: Stab platziert eine senkrechte Spalte aus 100 °C heißen Zellen nahe dem linken Rand und eine 0 °C kalte Wand entlang des rechten Randes, mit isolierenden oberen und unteren Wänden, was einfache Leitung durch eine Platte zwischen einer Wärmequelle und einer Wärmesenke veranschaulicht. Voreinstellung: Raum platziert einen 80 °C heißen „Heizkörper"-Streifen entlang des Bodens und einen 0 °C kalten „Fenster"-Streifen oben rechts, umschlossen von isolierenden Wänden, und zeigt, wie sich Wärme ungleichmäßig in einem Raum ausbreitet.

Warum begrenzt die Simulation die Wirkung des Diffusivitätsreglers auf die Geschwindigkeit?

Explizite Finite-Differenzen-Schemata für die Wärmegleichung sind nur numerisch stabil, wenn α·Δt/Δx² unter 0,25 bleibt; oberhalb dieser Schwelle verstärken sich kleine Fehler mit jedem Schritt und das Temperaturfeld oszilliert und divergiert, anstatt sich zu glätten. Der Code erzwingt dies, indem er das effektive Verhältnis unabhängig davon, wie hoch der Diffusivitätsregler eingestellt ist, auf 0,24 begrenzt.

Häufig gestellte Fragen

Was stellt der Laplace-Operator in der Wärmegleichung dar?

Der Laplace-Operator ∇²u misst den Unterschied zwischen dem Wert von u an einem Punkt und seinem Mittelwert in einer umgebenden Nachbarschaft. Wo ∇²u > 0 ist (ein lokales Minimum), fließt Wärme hinein und die Temperatur steigt; wo ∇²u < 0 ist (ein lokales Maximum), fließt Wärme heraus und die Temperatur sinkt.

Was ist die stationäre Lösung der Wärmegleichung?

Für t → ∞ verschwindet die Zeitableitung und die Wärmegleichung reduziert sich auf die Laplace-Gleichung ∇²u = 0. Die stationäre Temperaturverteilung ist harmonisch — glatt, ohne lokale Maxima oder Minima innerhalb des Gebiets.

Wie hängt die Wärmegleichung mit Zufallsbewegungen zusammen?

Die Wärmegleichung ist der Kontinuumsgrenzwert einer Zufallsbewegung auf einem Gitter. Wenn Teilchen zufällig in gleichen Schritten diffundieren, erfüllt die Wahrscheinlichkeitsdichte, ein Teilchen an Position x zur Zeit t zu finden, dieselbe Gleichung wie das Temperaturfeld — eine Verbindung zwischen Diffusion, Wahrscheinlichkeitstheorie und Wärmeübertragung.

Welche numerischen Methoden lösen die Wärmegleichung?

Gängige Methoden sind explizite finite Differenzen (Vorwärts-Euler, einfach, aber bedingt stabil), implizite Methoden (Rückwärts-Euler, unbedingt stabil) und das Crank-Nicolson-Schema (zweiter Ordnung genau und unbedingt stabil). Finite-Elemente-Methoden bewältigen unregelmäßige Gebiete.

Kann die Wärmegleichung die Abkühlung der Erde modellieren?

Ja. Geophysiker verwenden die Wärmegleichung, um zu modellieren, wie ursprüngliche Wärme durch die Erdkruste und den Erdmantel entweicht. Lord Kelvin nutzte sie berühmterweise, um das Alter der Erde auf 20-100 Millionen Jahre zu schätzen, obwohl ihm die radioaktive Erwärmung, die heute die Hauptquelle der inneren Wärme ist, nicht bekannt war.

⚙ Unter der Haube

Simulieren Sie die 2D-Wärmeleitung mit der Finite-Differenzen-Methode. Zeichnen Sie Wärmequellen und -senken, beobachten Sie, wie die Temperatur diffundiert, und vergleichen Sie stationäre mit instationären Lösungen.

heat equation simulation2d temperature distributionfinite difference methodthermal conductionsteady statetransient heat transfer

3D · Three.js / WebGL-Renderer · Ziel 60 FPS · läuft vollständig clientseitig, keine Installation nötig

Was haben Sie gefunden?

Schritte zur Reproduktion hinzufügen (optional)