или
div
−
D
k
grad
(
n
k
) +
v
k
n
k
+
∂
∂t
n
k
=
ψ
k
.
(3)
Длярассматриваемой физической ситуации функция
ψ
k
определя-
етсяследующими процессами:
1) подтеканием производственного газа UF
6
;
2) воздухообменом;
3) гидролизом;
4) коагуляцией.
Пусть
F
1
— плотность мощности внешних источников вещества
UF
6
(газ),
K
— кратность воздухообмена,
τ
1
,
τ
2
— периоды полувыве-
дениясоответственно веществ UF
6
(газ), UOF
4
(газ) за счет гидроли-
за,
τ
3
,
τ
5
— периоды полувыведениявеществ
UO
2
F
2
(газ), HF (газ)
за счет коагуляции,
λ
1
= (ln 2)
/τ
1
,
λ
2
= (ln 2)
/τ
2
,
λ
3
= (ln 2)
/τ
3
,
λ
5
= (ln 2)
/τ
5
. Тогда плотность мощности источников и стоков для
каждого вещества можно записать в виде:
ψ
1
=
−
(
λ
1
+
K
)
n
1
+
F
1
, ψ
2
=
λ
1
n
1
−
(
λ
2
+
K
)
n
2
,
ψ
3
=
λ
2
n
2
−
(
λ
3
+
K
)
n
3
, ψ
4
=
λ
3
n
3
−
Kn
4
,
ψ
5
= 2
λ
1
n
1
+ 2
λ
2
n
2
−
(
λ
5
+
K
)
n
5
, ψ
6
=
λ
5
n
5
−
Kn
6
.
Здесь процессы гидролиза, коагуляции и воздухообмена обусловли-
вают появление в выражении для
ψ
k
слагаемых, пропорциональных
n
1
, . . . , n
k
, а подтекание дает слагаемое, не зависящее от
n
1
, . . . , n
k
.
Пусть
D
k
≡
const,
v
k
≡
const. Тогда система уравнений непрерыв-
ности может быть записана в виде
∂
∂t
n
1
=
D
1
∆
n
1
−
(
v
1
,
grad
(
n
1
))
−
(
λ
1
+
K
)
n
1
+
F
1
(
x, t
)
,
∂
∂t
n
2
=
D
2
∆
n
2
−
(
v
2
,
grad
(
n
2
)) +
λ
1
n
1
−
(
λ
2
+
K
)
n
2
,
∂
∂t
n
3
=
D
3
∆
n
3
−
(
v
3
,
grad
(
n
3
)) +
λ
2
n
2
−
(
λ
3
+
K
)
n
3
,
∂
∂t
n
4
=
D
4
∆
n
4
−
(
v
4
,
grad
(
n
4
)) +
λ
3
n
3
−
Kn
4
,
∂
∂t
n
5
=
D
5
∆
n
5
−
(
v
5
,
grad
(
n
5
)) + 2
λ
1
n
1
+ 2
λ
2
n
2
−
(
λ
5
+
K
)
n
5
,
∂
∂t
n
6
=
D
6
∆
n
6
−
(
v
6
,
grad
(
n
6
)) +
λ
5
n
5
−
Kn
6
, x
∈
Q, t
∈
(
t
0
, t
1
)
.
(4)
Здесь
∆
— оператор Лапласа.
98
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2005. № 4
