и неотражающие граничные условия на выходе из расчетной области
x
=
L
:
∂u
∂x
=
∂v
∂x
=
∂T
∂x
=
∂ρ
∂x
= 0
.
(15
б
)
Заметим
,
что при использовании криволинейной системы координат
граничные условия второго рода формулируются для изменения функ
-
ций вдоль подходящей координаты
.
На поверхности обтекаемого тела задаются условия прилипания
y
= 0 :
u
=
v
= 0
, T
=
T
w
,
∂p
∂y
= 0
.
(15
в
)
Для решения системы уравнений
(1)–(4)
воспользуемся методом
расщепления по физическим процессам
:
расчеты на каждом времен
-
ном слое будем проводить в два этапа
.
На первом этапе интегрируются
уравнения неразрывности и движения
,
а на втором этапе
—
уравнение
сохранения энергии
.
При необходимости организуются дополнитель
-
ные внутренние итерации между указанными двумя этапами
.
Рассмотрим преобразование системы уравнений
(1)–(4)
для ее
интегрирования в произвольной криволинейной системе координат
ξ
=
ξ
(
x, y
)
,
η
=
η
(
x, y
)
.
Уравнения
(1)–(3)
представим в следующем
виде
:
∂
f
∂t
+
∂
E
∂ξ
+
∂
F
∂η
=
∂
E
v
∂ξ
+
∂
F
v
∂η
+ H;
(16)
f =
1
J
ρ
ρu
ρv
;
(17)
E =
1
J
ρU
ρuU
+
ξ
x
p
ρvU
+
ξ
y
p
;
(18)
F =
1
J
ρV
ρuV
+
η
x
p
ρvV
+
η
y
p
;
(19)
E
v
=
1
J
∙
Re
0
ξ
x
τ
xx
+
ξ
y
τ
xy
ξ
x
τ
xy
+
ξ
y
τ
yy
;
(20)
F
v
=
1
J
∙
Re
0
η
x
τ
xx
+
η
y
τ
xy
η
x
τ
xy
+
η
y
τ
yy
;
(21)
ISSN 0236-3941.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Машиностроение
”. 2005.
№
3 11