Numerische Modellierung und Bewertung der Umsetzbarkeit ... - GSI

Simulation der vereinfachten thermischen Abschirmung
− ϕ(uMP , vj)∆tλ
2∆v2
T
ρc
l1:=
n+1/2 n+1/2
j+1,k − 2Tj,k −l1T n+1/2
j+1,k
= −T n+1/2
j,k
+(1+2l1)T n+1/2
j,k
Lineares Gleichungssystem:
+ T n+1/2
j−1,k
− ∆tλ
2∆z2 n
Tj,k+1 − 2T
ρc
l2:=
n j,k + T n
j,k−1
+ T n j,k − ∆t
ɛAσ(T
dρc
4 j,k − T 4 A) + ɛIσ(T 4 j,k − T 4 I )
S2(uMP ,vj):=
n+1/2
−l1T
j−1,k = l2T n j,k+1+(1−2l2)T n j,k+l2T n j,k−1−S2(uMP , vj)
A11 ∗ T n+1/2 = A12 ∗ T n − b1
Diskretisierung der Punkte auf dem Rand δΩ :
T (vj,z3)+4T (vj,z2)−3T (vj,z1)
∆z = 0
T (vj, zOP ) = TK
−T (v3,zk)+4T (v2,zk)−3T (v1,zk)
∆v = 0
3T (vNP ,zk)−4T (vNP −1,zk)+T (vNP −2,zk)
∆v = 0
Zweiter Teilschritt
Diskretisierung der Punkte innerhalb des Gebiets Ω :
23
für alle j = 1, . . . , NP
für alle k = 2, . . . , OP − 1.
dλ ϕ(uMP , vj) ∂2T ∂v2 + ∂2T ∂z2
= dρc ∂T
∂t + ɛAσ T 4 j,k − T 4 A + ɛIσ T 4 j,k − T 4 A
⎡
⎢
dλ ⎣ϕ(uMP , vj)
T n+1/2 n+1/2
j+1,k −Tj,k vj+1−vj
−
n+1/2
Tj,k −T n+1/2
j−1,k
vj−vj−1
1/2(vj+1 − vj−1)
+
T n+1 n+1
j,k+1−Tj,k zk+1−zk
n+1 n+1
Tj,k −Tj−1,k−1 −
zk−zk−1
1/2(zk+1 − zk−1)
(3.12)
⎤
⎥
⎦