Roland 5000s Manual de usuario Pagina 103

  • Descarga
  • Añadir a mis manuales
  • Imprimir
  • Pagina
    / 195
  • Tabla de contenidos
  • MARCADORES
  • Valorado. / 5. Basado en revisión del cliente
Vista de pagina 102
§3 Diskretisierung des gekoppelten Systems 75
tor abh
¨
angen, erh
¨
alt man nichtlineare FEM–FCT–Verfahren, die in [KMT04, KM05]
beschrieben sind. F
¨
ur den Fall, dass die Gewichte alle gleich Eins sind, erh
¨
alt man die
Galerkin–Finite–Element–Methode (III.17) zur
¨
uck [JS08].
Dem linearen FEM–FCT–Verfahren, dass in [Kuz09] beschrieben wird, liegt folgende
Idee zugrunde: Man ersetzt u
k
in Gleichung (III.22) durch eine Approximation, die
mit einem expliziten Schema berechnet wird. Zu diesem Zweck, definiert man einen
Zwischenwert
u
k1/2
:=
u
k
+ u
k1
2
und setzt diesen in (III.22) ein
r
ij
= 2m
ij
u
k1/2,i
u
k1,i
2m
ij
u
k1/2,j
u
k1,j
t
k
d
ij
u
k1/2,i
u
k1/2,j
. (III.23)
Die Approximation von u
k1/2
kann mittels des Vorw
¨
arts–Euler–Verfahrens berechnet
werden, indem man im diskreten Zeitpunkt t
k1
einen halben Zeitschritt t
k
/2 macht,
˜u
= u
k1
t
k
2
M
1
L
Lu
k1
f
k1
. (III.24)
Die numerischen Fl
¨
usse des linearen FEM–FCT–Verfahrens erh
¨
alt man dann durch
Einsetzen dieser Approximation in (III.23),
r
ij
= t
k
m
ij
v
k1/2,i
v
k1/2,j
d
ij
˜u
i
˜u
j

,
mit den abk
¨
urzenden Schreibweisen
v
k1/2,i
=
M
1
L
f
k1
Lu
k1

i
,
˜u
i
= u
k1
+
t
k
2
v
k1/2,i
.
Zur Berechnung der Gewichte wird der Algorithmus von Zalesak [Zal79] benutzt.
Eine Motivation f
¨
ur diesen, findet sich zusammen mit einer ausf
¨
uhrlichen Diskussion in
[KM05]. Analog zu [JS08] wird hier aus Gr
¨
unden der Vollst
¨
andigkeit der Algorithmus
nur kurz wiedergegeben:
1. Berechne die Summen der positiven und der negativen antidiffusiven Fl
¨
usse im
Knoten i
P
+
i
=
N
X
j=1,j6=i
max {0, r
ij
} , P
i
=
N
X
j=1,j6=i
min {0, r
ij
} ,
Vista de pagina 102
1 2 ... 98 99 100 101 102 103 104 105 106 107 108 ... 194 195

Comentarios a estos manuales

Sin comentarios