Развитие метода ленточно-адаптивных сеток на основе схем TVD для решения задач газовой динамики - page 4

вид:
R
i
k
=
 
1
1
0
0
1
ˉ
v
1
as
1
ˉ
v
1
s
2
s
3
ˉ
v
1
+
as
1
ˉ
v
2
as
2
ˉ
v
2
s
1
0
ˉ
v
2
+
as
2
ˉ
v
3
as
3
ˉ
v
3
0
s
1
ˉ
v
3
+
as
3
H
a
v
i
g
ii
|
v
|
2
2
ˉ
v
1
s
2
ˉ
v
2
s
1
ˉ
v
3
s
1
ˉ
v
1
s
3
H
+
a
v
i
p
g
ii
 
,
(7)
если
s
1
6
= 0
;
R
i
k
=
 
1
1
0
0
1
ˉ
v
1
as
1
ˉ
v
1
s
2
0
ˉ
v
1
+
as
1
ˉ
v
2
as
2
ˉ
v
2
s
1
s
3
ˉ
v
2
+
as
2
ˉ
v
3
as
3
ˉ
v
3
0
s
2
ˉ
v
3
+
as
3
H
a
v
i
p
g
ii
|
v
|
2
2
ˉ
v
1
s
2
ˉ
v
2
s
1
ˉ
v
2
s
3
ˉ
v
3
s
2
H
+
a
v
i
p
g
ii
 
,
(8)
если
s
2
6
= 0
;
R
i
k
=
 
1
1
0
0
1
ˉ
v
1
as
1
ˉ
v
1
s
3
0
ˉ
v
1
+
as
1
ˉ
v
2
as
2
ˉ
v
2
0
s
3
ˉ
v
2
+
as
2
ˉ
v
3
as
3
ˉ
v
3
s
1
s
2
ˉ
v
3
+
as
3
H
a
v
i
p
g
ii
|
v
|
2
2
ˉ
v
3
s
1
ˉ
v
1
s
3
ˉ
v
2
s
3
ˉ
v
3
s
2
H
+
a
v
i
p
g
ii
 
,
(9)
если
s
3
6
= 0
. Здесь
s
j
=
P
i
j
/
p
g
ii
.
Столбец
Φ
i,n
k
+1
/
2
численной диссипации имеет следующие компо-
ненты:
Φ
i,n
k
+1
/
2
= (
φ
i,n
1
,k
+1
/
2
, . . . , φ
i,n
5
,k
+1
/
2
)
T
;
(10)
φ
i,n
j,k
+1
/
2
=
g
i,n
j,k
+
g
i,n
j,k
+
ψ
(
λ
i
a
i,n
j,k
+1
/
2
+
γ
i,n
j,k
+1
/
2
)
α
i,n
j,k
+1
/
2
;
g
i,n
j,k
=
s
i,n
j,k
+1
/
2
max[0
,
min(
|
˜
g
i,n
j,k
+1
/
2
|
,
˜
g
i,n
j,k
+1
/
2
s
i,n
j,k
+1
/
2
)];
s
i,n
j,k
+1
/
2
= sign(˜
g
i,n
j,k
+1
/
2
);
90
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2011. № 2
1,2,3 5,6,7,8,9,10,11
Powered by FlippingBook