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

Sξ
i,j+1/2 =
1
2
M U
i,j+1/2 + M U
i,j+1/2  ρa/J
ρua/J
ρva/J  
i,j
+
+
1
2
M U
i,j+1/2 − M U
i,j +1/2  ρa/J
ρua/J
ρva/J  
i,j+1
, (33)
Sξ
i,j−1/2 =
1
2
M U
i,j−1/2 + M U
i,j−1/2  ρa/J
ρua/J
ρva/J  
i,j−1
+
+
1
2
M U
i,j−1/2 − M U
i,j−1/2  ρa/J
ρua/J
ρva/J  
i,j
, (34)
Sη
i+1/2,j =
1
2
M V
i+1/2 ,j + M V
i+1/2,j  ρa/J
ρua/J
ρva/J  
i,j
+
+
1
2
M V
i+1/2 ,j − M V
i+1/2 ,j  ρa/J
ρua/J
ρva/J  
i+1,j
, (35)
Sη
i−1/2,j =
1
2
M V
i−1/2,j + M V
i−1/2,j  ρa/J
ρua/J
ρva/J  
i−1 ,j
+
+
1
2
M V
i−1/2,j − M V
i−1/2,j  ρa/J
ρua/J
ρva/J
 
i,j
, (36)
M U
i,j+1/2 = M U
i,j+1/2
+
L + M U
i,j+1/2
−
R , (37)
M U
i,j+1/2
+
L =  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;
(38)
M U
i,j+1/2
−
R =  −
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;
(39)
ISSN 0236-3941. Вестник МГТУ им. Н .Э. Баумана . Сер . “Машиностроение”. 2005. № 3 15