Laplace-Differentialgleichung lösen mit FEM

Mit JavaScript wird die Methode der finiten Elemente (FEM) für online-Berechnungen der stationären Wärmeleitung angewendet.
Die stationäre Wärmeleitung wird durch die folgende partielle Differentialgleichung beschrieben:
λ(x,y) (∂²T/∂x² + ∂²T/∂y²) = 0.
Die Lösung wird mit Hilfe von 3-Knoten-Dreieckelementen und isoparametrischen 4-Knoten-Viereckselementen bestimmt.
Die Ansatzfunktionen beider Elemente haben entlang ihrer Kanten lineare Verläufe. Sie lassen sich daher problemlos kombinieren.
Als Lösung werden die Temperaturen an den Knoten bestimmt und zusätzlich die Wärmeströme an den Knoten, die einen Wärmeaustausch nach draußen haben.
Ferner werden alle Elementmatrizen und das Gesamtgleichungssystem ausgegeben.

 

Zu den Modelldaten

Bei der Elementeingabe wird pro Element eine Liste von 3 oder 4 Punkten erwartet. Dadurch wird indirekt der Typ des Elements festgelegt.
Beim Vierknotenelement ist darauf zu achten, dass diese Liste eine umlaufende Reihenfolge hat. Wo man beginnt ist dabei egal.
Ferner wird die Wärmeleitfähigkeit λ erwartet.
Elemente mit unterschiedlicher Wärmeleitfähigkeit werden hier durch unterschiedlich starke rot-Färbung kenntlich gemacht.
Hinweis: Bei der Elementgenerierung ist darauf zu achten, dass Knoten entlang der gemeinsamen Kante zweier Elemente auch Knoten beider Elemente sind.
Bei der Temperatureingabe und der Wärmequelleneingabe werden alle Knoten farblich hervorgehoben, die eine der beiden Vorgaben besitzen.
Rot sind Knoten mit Wärmequelle, blau sind Knoten mit Temperaturvorgabe.

Zu den Ergebnissen

Dargestellt wird der Verlauf der Temperatur T(x,y) im Gebiet.
Als Tabelle ausgegeben wird ferner die Temperatur an allen Knoten und die Wärmeströme zwischen dem untersuchten Gebiet und der Umgebung.
Falls dem Gebiet Energie zufließt (Wärmequellen) ist der Wärmestrom positiv.
Falls aus dem Gebiet Energie abfließt (Wärmesenken) ist der Wärmestrom negativ.

Zu den Einheiten

Die Punkt-Koordinaten sind in der Einheit m zu verstehen.
Die Wärmeleitfähigkeit λ hat die Einheit W/mK.
Die gleiche Einheit haben die Elemente der Steifigkeitsmatrix.
Die Knotenwärmeströme Qi sind dann von der Einheit W/m.
Das hier noch die Längeneinheit im Nenner erscheint ist Folge der ebenen Betrachtungsweise.
Die Knotenströme sind dann eigentlich die auf die Materialdicke bezogenen Ströme.

weitere JavaScript-Programme