Метод расчета сверхзвукового обтекания сферы на основе AUSM конечно-разностных схем - page 10

M
U
i,j
1/2
=
M
U
i,j
1/2
+
L
+
M
U
i,j
1/2
R
;
(40)
M
U
i,j
1/2
+
L
=
 
1
4
M
U
i,j
1
+ 1
2
при
M
U
i,j
1
6
1
,
1
2
M
U
i,j
1
+
M
U
i,j
1
при
M
U
i,j
1
>
1;
(41)
M
U
i,j
1/2
R
=
 
1
4
M
U
i,j
1
2
при
M
U
i,j
6
1
,
1
2
M
U
i,j
M
U
i,j
при
M
U
i,j
>
1;
(42)
M
U
i,j
=
U
i,j
a
i,j
,
a
i,j
=
r
γ
p
i,j
ρ
i,j
.
(43)
Конечно
-
разностные формулы для
S
η
i
+1/2
,j
,
S
η
i
1/2
,j
, M
V
i,j
записыва
-
ются аналогично
.
Аппроксимация производных потоков импульса
,
обусловленного
действием давления
,
используется в виде
E
p
∂ξ
=
S
p,ξ
i,j
+1/2
S
p,ξ
i,j
1/2
ξ
j
+1/2
ξ
j
1/2
;
(44)
F
p
∂η
=
S
p,η
i
+1/2
,j
S
p,η
i
1/2
,j
η
i
+1/2
η
i
1/2
;
(45)
S
p,ξ
i,j
+1/2
=
p
+
L
 
0
ξ
x
/
J
ξ
y
/
J
 
i,j
+
p
R
 
0
ξ
x
/
J
ξ
y
/
J
 
i,j
+1
;
(46)
p
+
L
=
 
1
2
p
i,j
1 +
M
U
i,j
при
M
U
i,j
6
1
,
1
2
p
i,j
M
U
i,j
+
|
M
U
i,j
|
M
U
i,j
при
M
U
i,j
>
1;
(47)
p
R
=
 
1
2
p
i,j
+1
1
M
U
i,j
+1
при
M
U
i,j
+1
6
1
,
1
2
p
i,j
+1
M
U
i,j
+1
|
M
U
i,j
+1
|
M
U
i,j
+1
при
M
U
i,j
+1
>
1;
(48)
S
p,ξ
i,j
1/2
=
p
+
L
 
0
ξ
x
/
J
ξ
y
/
J
 
i
1
,j
+
p
R
 
0
ξ
x
/
J
ξ
y
/
J
 
i,j
;
(49)
p
+
L
=
 
1
2
p
i,j
1
1 +
M
U
i,j
1
при
M
U
i,j
1
6
1
,
1
2
p
i,j
1
M
U
i,j
1
+
M
U
i,j
1
M
U
i,j
1
при
M
U
i,j
1
>
1;
(50)
16 ISSN 0236-3941.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Машиностроение
”. 2005.
3
1,2,3,4,5,6,7,8,9 11,12,13,14,15,16,17,18,19,20,...27
Powered by FlippingBook