Abbildung 1: Kupferplatte mit thermischen und elektrischen

Praktikum 4
zur Vorlesung „Multi-Physik-Simulationen“
Prof. Dr. Christian Schröder
Multi-Physik-Simulation I – Thermo-Elektrische Kopplung – Joulesche Erwärmung eines
Kupferblechs
Jedes elektrisch leitfähige Material erwärmt sich aufgrund des Ohmschen Widerstands, sobald ein
Strom hindurch fließt. Dies nennt man Joulesche Erwärmung. In den vorangegangenen
Praktikumsaufgaben lagen jeweils Systeme vor, die durch eine einzige Differentialgleichung,
nämlich die Wärmeleitungsgleichung, beschrieben werden können.
Für die Simulation der Jouleschen Erwärmung ist es allerdings nötig, den Effekt der Wärmeleitung
mit dem Effekt der elektrischen Leitfähigkeit zu koppeln. Hierfür gibt es in COMSOL bereits einen
vordefinierten Modus „Joule Heating“.
Untersuchen Sie diese Kopplung mit Hilfe von COMSOL anhand einer 1 Meter x 1 Meter dünnen
(wie dünn?) Kupferplatte, an die eine Spannung von 0,1 V angelegt wird. Die Platte wird von einem
Luftstrom gekühlt, so dass die Temperatur am Rand 300 K beträgt (siehe Abb. 1).
Abbildung 1: Kupferplatte mit thermischen und elektrischen Randbedingungen
Aufgabe 1:
Führen Sie mit Hilfe von COMSOL eine stationäre, zweidimensionale Simulation des oben
beschriebenen Problems durch. Nutzen Sie dafür den in COMSOL vordefinierten Modus Joule
Heating.
•
•
•
•
Erzeugen Sie ein neues Projekt.
Wählen Sie die Optionen 2D, Heat Transfer → Joule Heating (jh), Stationary.
Zeichnen Sie die Geometrie: Square, Side length = 1m.
Wählen Sie das Material Kupfer (Copper) aus der Materialbibliothek und weisen Sie es
dem Gebiet der Kupferplatte zu.
•
•
Übertragen Sie die folgende Zeile in die Tabelle unter Parameters.
V0
0.1[V]
0.1[V]
Elektrisches Potential
Übertragen Sie die Randbedingungen aus Abb. 1 und der Aufgabenbeschreibung in
COMSOL.
• Für die elektrischen Randbedingungen müssen Sie zwei Knoten Electric Potential
im Model Builder erzeugen. Sie finden diesen Bedingungstyp unter Electric
Currents → Electric Potential.
• Tragen Sie unter Voltage für V0 den jeweiligen Wert aus der Abbildung ein.
• Für die thermischen Randbedingungen müssen Sie einen Knoten Temperature
(Heat Transfer → Temperature) erzeugen und für T0 den in Abb. 1 gezeigten Wert
eintragen.
• Passen Sie die Gebietsbedingung Joule Heating Model 1 an. Der entsprechende
Knoten sollte bereits automatisch von COMSOL erzeugt worden sein. Wählen Sie
unter Conduction Current für Electric Conductivity Linearized resistivity aus.
Sonst müssen Sie nichts verändern.
• Vernetzen Sie die Geometrie sinnvoll und führen Sie die Berechnung durch.
• Stellen Sie als Ergebnis den Temperaturverlauf über der Platte in x-Richtung
grafisch dar.
Aufgabe 2:
Mit COMSOL ist es nicht nur möglich vordefinierte Kopplungen zu nutzen. Sie können auch selbst
Kopplungen erzeugen. Führen Sie die gleiche Simulation nochmals in einem neuen COMSOLProjekt durch und greifen Sie dabei nicht auf den vordefinierten Modus Joule Heating zurück! Auf
diese Weise bekommen Sie einen ersten Eindruck von den Möglichkeiten, die COMSOL im
Hinblick auf Multi-Physik-Kopplungen bietet.
Bei diesem Beispiel handelt es sich um eine sogenannte Multi-Feld-Kopplung (in einem
gemeinsamen Gebiet). Die Formulierung als System partieller Differentialgleichungen lautet in
diesem (vereinfachten) Fall:
−∇ k ∇ T −Q V =0
∇ T  ∇ V =0
Die Kopplung erfolgt über die Gleichungen
QV = T ∣∇ V∣2 Joulesche Erwärmung
und
T =
0
Leitfähigkeit.
1T −T 0
Die letzte Gleichung beschreibt die Temperaturabhängigkeit der Leitfähigkeit wie sie für Metalle bei
nicht zu niedrigen Temperaturen zu beobachten ist. Hierbei ist der Referenzwiderstand
0 =1.72e-8 m an der Referenztemperatur T 0=298 K . Der Temperaturkoeffizient ist für
Kupfer
=0.0039
1
.
K
Statt die vordefinierte Kopplung aus der ersten Aufgabe zu benutzen, sollen Sie die Kopplung nun
selbst implementieren. Dazu müssen Sie zwei Ausdrücke Q_selfmade und sigma_selfmade in
Ihr Modell einpflegen. Formeln für beide Ausdrücke sind oben gegeben. Sie müssen die Ausdrücke
für Q(V) (das ist Ihr Q_selfmade) und T  (das ist Ihr sigma_selfmade) jetzt so
formulieren, dass nur Variablen benutzt werden, die COMSOL auch bekannt sind!
Sie können dazu selbstdefinierte Größen (Parameters) und die abhängigen Variablen Ihrer
Physikmodi (das sind hier die Temperatur T und das elektrische Potential V) benutzen.
Überlegen Sie sich, wie die Ausdrücke für Q_selfmade und sigma_selfmade in COMSOL
aussehen müssen!
Hinweis: Sie benötigen dafür die Komponenten des elektrischen Potentials. Die beiden benötigten
Komponenten sind einfach über Vx und Vy zu verwenden.
Vorgehensweise in COMSOL:
•
•
•
•
•
•
•
•
•
•
•
•
•
•
•
Starten Sie ein neues Projekt in COMSOL.
Wählen Sie die Option 2D.
Fügen Sie zwei Physikmodi zum Modell hinzu.
• Heat Transfer
• Electric Currents
Wählen Sie die Option Stationary.
Zeichnen Sie die Geometrie (Square mit Seitenlänge 1m).
Legen Sie die folgenden Konstanten unter Global Definitions → Parameters an.
• v0: 0.1[V]
• sigma0: 5.814e7[S/m] (Das ist der Kehrwert von der gegebenen Größe 0 )
• T0: 298[K]
• beta: 0.0039[1/K]
Fügen Sie das Material Kupfer (Copper) zum Modell hinzu.
Geben Sie die selbstdefinierten Größen Q_selfmade und sigma_selfmade unter Model
1 → Definitions → Variables 1 ein. Diesen Knoten müssen Sie selbst erzeugen. Denken
Sie daran, dass Sie COMSOL bekannte Größen verwenden müssen (s.o.).
Nachdem Sie die Kopplungsgrößen definiert haben, können Sie sie nun in den einzelnen
Physikmodi verwenden. Legen Sie dazu zunächst einen Knoten Heat Source 1 unter Heat
Transfer (ht) an.
Tragen Sie unter General Heat Source für Q die zuvor definierte Größe Q_selfmade ein.
In einer Richtung haben Sie somit die Kopplung zwischen den beiden Modi definiert.
COMSOL bestimmt die Größe Q des Modus Heat Transfer nun in Abhängigkeit der Größe
 des Modus Electric Currents.
Es gibt allerdings auch eine Kopplung in die andere Richtung. Dadurch, dass  unter
Electric Currents temperaturabhängig ist, hängt diese Größe wiederum vom Modus Heat
Transfer ab, da dort die Temperatur berechnet wird.
Tragen Sie unter Electric Currents → Current Conservation 1 für  die selbst
definierte Größe sigma_selfmade ein. Die Kopplung zwischen den beiden Modi ist somit
vervollständigt!
Ergänzen Sie nun noch die Randbedingungen aus Aufgabenteil 1!
Vernetzen Sie das Modell und berechnen Sie die Lösung.
Vergleichen Sie die Ergebnisse mit denen aus Aufgabenteil 1.

Download Report