Метод расчета сверхзвукового обтекания сферы на основе AUSM конечно-разностных схем - page 14
Q
(2)
dis,i,j
= 2
ξ
y,i,j
+1
v
i,j
+1
J
i,j
+1
−
ξ
y,i,j
−
1
v
i,j
−
1
J
i,j
−
1
1
2
q
j
+
+
η
y,i
+1
,j
v
i
+1
,j
J
i
+1
,j
−
η
y,i
−
1
,j
v
i
−
1
,j
J
i
−
1
,j
1
2
p
i
2
;
(76)
Q
(3)
dis,i,j
= 2
α
v
i,j
y
i,j
J
i,j
2
;
(77)
Q
(4)
dis,i,j
=
U
i,j
+1
J
i,j
+1
−
U
i,j
−
1
J
i,j
−
1
1
2
q
j
+
V
i
+1
,j
J
i
+1
,j
−
V
i
−
1
,j
J
i
−
1
,j
1
2
p
i
2
;
(78)
Q
(5)
dis,i,j
=
−
2
3
U
i,j
+1
J
i,j
+1
−
U
i,j
−
1
J
i,j
−
1
1
2
q
j
+
+
V
i
+1
,j
J
i
+1
,j
−
V
i
−
1
,j
J
i
−
1
,j
1
2
p
i
+
α
v
i,j
y
i,j
J
i,j
2
;
(79)
Q
p,i,j
=
U
i,j
p
i,j
+1
−
p
i,j
−
1
2
q
j
+
V
i,j
p
i
+1
,j
−
p
i
−
1
,j
2
p
i
1
J
i,j
;
(80)
Q
xy,i,j
= 2
k
i,j
T
s
+1
i
+1
,j
+1
−
T
s
+1
i
−
1
,j
−
1
+
T
s
+1
i
−
1
,j
+1
+
T
s
+1
i
+1
,j
−
1
4
p
i
q
j
+
+
(
k
i,j
+1
−
k
i,j
−
1
)
T
s
+1
i
+1
,j
−
T
s
+1
i
−
1
,j
4
p
i
q
j
+
+
(
k
i
+1
,j
−
k
i
−
1
,j
)
T
s
+1
i,j
+1
−
T
s
+1
i,j
−
1
4
p
i
q
j
;
(81)
Q
α,i,j
=
α
λ
i,j
Re
Pr
1
y
i,j
J
i,j
"
ξ
y,i,j
T
s
+1
i,j
+1
−
T
s
+1
i,j
−
1
2
q
j
+
η
y,i,j
T
s
+1
i
+1
,j
−
T
s
+1
i
−
1
,j
2
p
i
#
;
(82)
k
i,j
=
λ
i,j
Re
Pr
(
ξ
x
η
x
+
ξ
y
η
y
)
i,j
J
i,j
;
(83)
a
±
h
=
1
2
(
a
h,i,j
+
a
h,i
±
1
,j
) ;
b
±
h
=
1
2
(
b
h,i,j
+
b
h,i,j
±
1
) ;
(84)
a
h
=
λ
i,j
Re
Pr
η
2
x
+
η
2
y i,j
J
i,j
;
b
h
=
λ
i,j
Re
Pr
ξ
2
x
+
ξ
2
y i,j
J
i,j
;
(85)
p
+
i
=
η
i
+1
−
η
i
;
p
−
i
=
η
i
−
η
i
−
1
;
q
+
j
=
ξ
j
+1
−
ξ
j
;
q
−
j
=
ξ
j
−
ξ
j
−
1
;
(86)
Q
rad,i,j
= (
div
q
rad
)
i,j
L
0
ρ
0
u
0
c
p,
0
J
i,j
,
(87)
20 ISSN 0236-3941.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Машиностроение
”. 2005.
№
3
1...,4,5,6,7,8,9,10,11,12,13
15,16,17,18,19,20,21,22,23,24,...27
