∂ρ
g
∂t
+
r ∙
ρ
g
~v
g
= 0
,
∂ρ
g
~v
g
∂t
+
r ∙
(
ρ
g
~v
g
~v
g
+
p
E) =
μ
Δ
~v
g
,
∂ρ
g
E
g
∂t
+
r ∙
ρ
g
E
g
+
p
ρ
g
~v
g
=
λ
ΔΘ
g
+
ω
(1)
при следующих граничных условиях
:
а
)
на поверхности горения топлива
~v
∙
~n
=
−
1
ρ
(
ρ
S
−
ρ
)
D, ~v
∙
~τ
I
= 0
, I
= 1
,
2
,
Θ = Θ
g
;
(2)
б
)
на остальных твердых поверхностях
~v
= 0
, ~n
∙
λ
r
Θ
g
= 0;
(3)
в
)
в выходном сечении
ρ
g
=
ρ
ge
;
и начальных условиях при
t
= 0
:
θ
=
θ
0
, ρ
=
ρ
0
, ~v
= 0
.
(4)
Здесь
ρ
g
—
плотность газа
;
~v
g
—
вектор скорости
;
μ
—
коэффициент
вязкости
;
E
g
—
полная энергия
,
E
g
=
C
V
Θ
g
+
v
2
2
,
C
V
—
теплоемкость
при постоянном объеме
,
Θ
g
—
температура газа
,
v
2
—
квадрат моду
-
ля скорости
,
определяемый следующим образом
:
v
2
=
v
2
r
+
v
2
z
+
v
2
θ
;
Рис
. 1.
Расчетная схема задачи течения продуктов сгорания в РДТТ
:
1
—
твердотопливный заряд
;
2
—
корпус
;
3
—
сопло
;
4
—
мембрана
;
5
—
камера
сгорания
46
ISSN 1812-3368.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Естественные науки
”. 2005.
№
3
