где
q
rad
—
вектор плотности интегрального радиационного теплового
потока
.
В данном случае при расчете функций
U
i,j
, V
i,j
использовалась про
-
цедура сглаживания локальных пульсаций скорости
.
Например
,
сгла
-
женная функция
U
i,j
рассчитывалась по следующему алгоритму
:
а
)
Δ
+
= ˜
U
i,j
+1
−
˜
U
i,j
,
Δ
−
= ˜
U
i,j
−
˜
U
i,j
−
1
;
б
)
если
(Δ
+
Δ
−
)
>
0
,
то
U
i,j
—
не модифицируется
;
в
)
если
(Δ
+
Δ
−
)
<
0
,
то
U
i,j
= (1
−
2
ω
) ˜
U
i,j
+
ω
˜
U
i,j
−
1
+ ˜
U
i,j
+1
,
где
ω
≤
0
,
5
,
˜
U
i,j
—
значение сеточной функции
,
получаемой в резуль
-
тате численного интегрирования уравнений движения
.
Схема Т
-2.
Это простейшая противотоковая схема
1-
го порядка
,
при получении которой интегро
-
интерполяционным методом
,
функции
ρc
p
U/J
и
ρc
p
V /J
полагаются постоянными в пределах элементарного
объема интегрирования
,
а их знак определяет способ расчета произ
-
водных температуры
.
Пятиточечная конечно
-
разностная схема
(67)
в
этом случае имеет следующие коэффициенты
:
A
i,j
=
a
+
c
p
−
i
+
a
−
h
p
−
i
p
i
;
B
i,j
=
−
a
−
c
p
+
i
+
a
+
h
p
+
i
p
i
;
(88)
ˉ
A
i,j
=
b
+
c
q
−
j
+
b
−
h
q
−
j
q
j
;
ˉ
B
i,j
=
−
b
−
c
q
+
j
+
b
+
h
q
+
j
q
j
;
(89)
C
i,j
=
c
p,i,j
ρ
i,j
J
i,j
τ
+
A
i,j
+
B
i,j
+ ˉ
A
i,j
+ ˉ
B
i,j
;
(90)
F
i,j
=
c
p,i,j
ρ
i,j
T
n
i,j
J
i,j
τ
+
Q
dis,i,j
+
Q
p,i,j
+
Q
α,i,j
+
Q
xy,i,j
+
Q
rad,i,j
,
(91)
где
a
±
c
=
1
2
(
a
i,j
± |
a
i,j
|
) ;
a
i,j
=
ρ
i,j
c
p,i,j
U
i,j
J
i,j
;
b
±
c
=
1
2
(
b
i,j
± |
b
i,j
|
) ;
b
i,j
=
ρ
i,j
c
p,i,j
V
i,j
J
i,j
.
Остальные функции такие же
,
как в схеме
Т
-1
.
Схема Т
-3.
В этой схеме для расчета температурного поля исполь
-
зуется
AUSM-
алгоритм
.
Здесь так же
,
как и в схеме
Т
-1
,
уравнение со
-
хранения энергии преобразуется к виду
ρc
p
∂T
∂t
+
∂
∂ξ
ρc
p
U
J
T
+
∂
∂η
ρc
p
V
J
T
−
−
T
∂
∂ξ
ρc
p
U
J
−
T
∂
∂η
ρc
p
V
J
=
R
T
.
(92)
ISSN 0236-3941.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Машиностроение
”. 2005.
№
3 21
