+
∂
∂x
(
σ
xx
u
+
σ
xz
v
)+
∂
∂z
(
σ
zx
u
+
σ
zz
v
)+
∂
∂x
κ
(
T
)
∂T
∂x
+
∂
∂z
κ
(
T
)
∂T
∂z
+
X
k
h
k
m
k
∂
∂x
ρD
k
(
T
)
∂c
k
∂x
+
∂
∂z
ρD
k
(
T
)
∂c
k
∂z
;
p
=
RTn
=
X
i
R
m
i
c
i
!
ρT
=
ρT
X
i
R
i
c
i
;
ε
=
c
v
T
+
X
k
h
k
n
k
ρ
=
c
v
T
+
X
k
h
k
c
k
;
σ
xx
= 2
μ
∂u
∂x
−
2
3
μ
∂u
∂x
+
∂v
∂z
;
σ
zz
= 2
μ
∂v
∂z
−
2
3
μ
∂u
∂x
+
∂v
∂z
;
σ
zx
=
σ
xz
=
μ
∂v
∂x
+
∂u
∂z
.
Здесь
ρ
— плотность;
p
— давление;
u, v
— компоненты скорости;
c
i
=
ρ
i
/ρ
— массовые концентрации отдельных компонентов, участву-
ющих в реакции;
E
=
ε
+ (
u
2
+
v
2
)
/
2
— плотность полной энергии
смеси;
ε
— плотность внутренней энергии;
R
— универсальная газовая
постоянная;
R
i
=
R/m
i
,
m
i
— масса моля
i
-го компонента;
n
— моляр-
ная плотность смеси;
σ
ij
— компоненты тензора вязких напряжений;
c
v
=
P
i
c
vi
c
i
— удельная теплоемкость смеси при постоянном объеме;
h
i
— удельная энтальпия образования
i
-го компонента;
κ
(
T
)
,
μ
(
T
)
—
коэффициенты теплопроводности и вязкости смеси соответственно;
D
i
(
T
)
— коэффициент диффузии
i
-го компонента в многокомпонент-
ной среде;
∂c
i
∂t
хим
— изменение концентрации
i
-го компонента за
счет химических превращений.
Для описания химической кинетики использовалась редуцирован-
ная схема окисления водорода [7]
Н
2
+
О
2
→
2ОН
;
Н
+
О
2
→
ОН
+
О
;
2Н
+
М
→
Н
2
+
М
;
2ОН
+
М
→
Н
2
О
2
+
М
;
Н
+
Н
2
О
2
→
Н
2
+
НО
2
;
ОН
+
Н
2
→
Н
2
О
+
H
;
Н
2
+
О
→
ОН
+
Н
;
40
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2008. № 4