Аэродинамика вентиляции. Механика Аэрозолей. Учебное пособие
.pdf
11
Первый член левой части этого уравнения представляет собой отнесенное к единичному объему изменение количества движения в этом объеме за единицу времени, второй член - отнесенное к единичному объему изменение количества движения за счет конвекции в этом объеме за единицу времени.
Первый член правой части есть отнесенная к единице объема массовая сила, второй член – отнесенные к единице объема поверхностные силы.
Используя уравнение неразрывности получаем следующее:
|
|
|
|
|
|
|
|
|
|
|
|
ρ |
¶V |
= ρ |
|
|
|
|
|
|
|
||
f |
-V |
×Ñ ×V + Ñ × Ï |
i, j |
||||||||
|
|
||||||||||
|
¶t |
(1.14) |
|||||||||
Для ньютоновских жидкостей напряжение на некоторой площадке пропорционально скорости деформации сплошной среды (жидкости). При этом связь между давлением, скоростью деформации и компонентами тензора напряжений имеет вид:
|
|
|
¶ui |
|
¶u |
|
|
|
¶uk |
|
|
Ï i , j |
= - pδi, j |
+ μ |
+ |
|
j |
|
+ δi, j μ ' |
, i , j ,k = 1, 2, 3 |
(1.15) |
||
|
|
|
|
||||||||
|
|
|
¶x j |
|
|
|
|
|
¶xk |
|
|
|
|
|
|
¶xi |
|
|
|||||
где δi,j – символ Кронекера (δi,j = 1, если i=j, и δi,j = 0, если i≠j), u1, u2, u3 –
компоненты вектора скорости V , x1, x2, x3 – координаты радиус-вектора точки, µ – коэффициент динамической вязкости, µ’ – второй коэффициент
вязкости, связанный с объемной вязкостью χ = 2 3 μ + μ '.
Обычно объемной вязкостью пренебрегаю (кроме случаев рассмотрения распространения ударных, акустических волн), тогда выражение для тензора напряжений можно записать в виде:
|
|
|
¶ui |
|
¶u |
|
|
|
2 |
|
¶uk |
|
|
||
Ï i , j |
= - pδi, j |
+ μ |
+ |
|
j |
|
- |
δi, j |
, i, j, k = 1, 2,3 |
(1.16) |
|||||
|
¶x |
|
|
¶x |
|||||||||||
|
|
|
|
¶x |
j |
|
|
3 |
|
|
|
||||
|
|
|
|
i |
|
|
|
|
k |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Тензор напряжений разделяют на две части: |
|
||||||||||||||
|
|
|
|
Ï i , j |
= - pδi. j +τ i, j |
, |
|
(1.17) |
|||||||
где первое слагаемой в правой части – компоненты нормальных напряжений, а второе – касательных или вязких:
|
|
¶ui |
|
τ i, j |
= μ |
||
¶x j |
|||
|
|
||
|
|
|
|
¶u |
|
|
|
2 |
δi, j |
¶uk |
|
|
+ |
|
j |
|
- |
, i, j, k = 1, 2,3 |
(1.18) |
|||
¶x |
|
|
¶x |
||||||
|
|
|
3 |
|
|
|
|||
|
i |
|
|
|
|
k |
|
||
|
|
|
|
|
|
|
|
|
|
12
Течение несжимаемой вязкой жидкости с постоянным коэффициентом вязкости описывается следующим уравнением:
|
|
|
|
|
|
|
|
|
|
|
|
|
ρ |
¶V |
= ρ |
|
|
|
|
|
|
|
|
(1.19) |
|
f |
- Ñp -V |
×Ñ ×V + μÑ2 V . |
||||||||||
|
¶t |
|
||||||||||
Или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
du |
|
|
|
|
|
¶uk |
|
|
||
|
¶V |
|
|
|
|
|
|
|
¶ |
|
|
|
|
2 |
|
|
|||||||||
ρ |
= ρ f - Ñp -V ×Ñ ×V + |
|
μ |
dui |
+ |
|
j |
|
- |
δi, j μ |
|
(1.20) |
|||||||||||||
¶t |
¶x |
|
¶x |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
¶x |
j |
|
|
|
3 |
|
¶x |
k |
|
|||||||
|
|
|
|
|
|
|
|
|
|
i |
|
|
i |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
В проекции на декартову систему координат имеем три скалярных уравнения:
|
ρ |
∂Vx |
|
|
|||
∂t |
|||
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂Vy |
|
|
ρ |
||
|
|
∂t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂Vz |
|
|
ρ |
||
|
|
∂t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
∂ |
|
|
|
2 |
|
|
|
|
|
|
|
∂ |
V |
|
|
∂V |
y |
|
V |
z |
|
|
|||||||||||||||||||
+ Vx |
∂ |
x |
|
+ Vy |
∂ |
|
|
x |
+ Vz |
|
|
∂ |
x |
|
= ρ fx − |
∂ |
|
+ |
|
|
|
|
|
|
μ 2 |
|
|
|
x |
− |
|
|
|
− |
∂ |
|
+ |
||||||||||||||||||||||||||||||||||||||||||||
∂x |
∂y |
|
|
|
|
|
∂z |
|
|
∂x |
∂x |
|
3 |
|
|
|
|
|
|
|
|
∂y |
|
|
∂z |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂Vy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
+ |
|
∂ |
|
|
|
μ |
|
∂Vx |
|
+ |
|
|
|
|
+ |
|
|
∂ |
|
μ |
∂Vz |
|
|
+ |
|
∂Vx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
∂y |
|
∂y |
|
|
∂x |
|
|
∂z |
|
∂x |
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
∂V |
y |
|
|
|
|
|
∂V |
y |
|
|
|
|
|
|
∂V |
y |
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
∂ |
|
|
|
|
|
|
∂V |
y |
|
|
|
|
|
|
V |
x |
|
|
|
|
|
|
|
|
||||||||||||||||||||
+ Vx |
|
|
|
+ Vy |
|
|
|
|
|
+ Vz |
|
|
|
|
|
|
= ρ f y − ∂ |
|
|
+ |
|
|
|
|
μ |
|
|
|
|
+ |
|
∂ |
|
+ |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
∂x |
∂y |
|
|
∂z |
|
|
|
|
|
∂x |
∂x |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.21) |
||||||
|
|
∂ |
|
|
|
2 |
|
|
|
|
∂V |
|
|
|
|
V |
|
|
|
|
|
|
∂ |
V |
|
|
|
∂ |
|
|
|
|
|
|
∂V |
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
+ |
|
|
|
|
μ 2 |
|
|
|
|
|
|
|
− |
∂ x |
|
− |
|
|
+ |
|
|
|
|
μ |
|
|
|
|
+ |
|
∂ |
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
∂y |
|
3 |
|
|
|
|
|
|
|
|
|
|
∂z |
∂z |
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
+ V |
∂Vz |
|
|
+ V |
|
|
∂Vz |
+ V |
|
∂Vz |
= ρ f |
|
|
− |
∂p + |
|
|
|
∂ |
μ |
|
∂Vz |
+ |
∂Vx |
|
|
+ |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
y |
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
x |
∂x |
|
|
|
∂y |
|
|
|
|
|
z |
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z ∂x |
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
∂ |
|
|
|
|
|
∂V |
y |
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
∂ |
|
2 |
|
|
|
|
∂ |
V |
z |
|
|
|
|
|
|
V |
x |
|
|
|
∂V |
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
+ |
|
|
|
|
μ |
|
|
|
|
+ |
|
∂ z |
|
|
|
+ |
|
|
|
|
|
|
μ |
2 |
|
|
|
− |
∂ |
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
∂y |
|
|
|
|
|
|
|
∂y |
|
|
|
|
∂z |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Выделив в уравнениях компоненты тензора вязких напряжений, получим:
ρ |
|
¶Vx |
+V |
|
|
¶Vx |
+V |
|
¶Vx |
+V |
¶Vx |
|
= ρ f - |
¶p + |
|
¶τ xx |
+ |
|
¶τ xy |
+ |
|
¶τ xz |
|
|
|||||||||||
|
¶t |
|
¶x |
|
¶y |
¶z |
|
¶x |
¶y |
|
¶z |
|
|
||||||||||||||||||||||
|
|
|
|
x |
|
|
|
y |
z |
|
|
|
x |
|
¶x |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
¶Vy |
|
|
¶Vy |
|
|
|
¶Vy |
|
|
¶Vy |
|
|
|
|
|
¶p |
|
|
|
¶τ yx |
|
|
|
¶τ yy |
|
|
|
¶τ yz |
|
|
|||
|
ρ |
+Vx |
+Vy |
+Vz |
|
|
= ρ f y - |
+ |
|
|
|
+ |
|
|
+ |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.22) |
|||||||||||||||
|
¶t |
|
¶x |
¶y |
|
¶z |
|
¶y |
¶x |
|
¶y |
|
|
¶z |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
ρ |
|
¶Vz |
|
|
|
¶Vz |
|
|
|
¶Vz |
|
|
¶Vz |
|
= ρ |
|
|
|
¶p |
|
|
¶τ zx |
|
|
|
¶τ zy |
|
|
¶τ zz |
|
|
|
||
|
|
V |
|
+ |
V |
|
V |
|
f |
|
- |
+ |
|
+ |
|
+ |
|
|
|
||||||||||||||||
|
¶t |
|
¶x |
|
¶y |
¶z |
|
|
¶z |
|
|
¶x |
|
¶y |
|
¶z |
|
|
|||||||||||||||||
|
|
+ |
x |
|
|
y |
+ z |
|
|
z |
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где компоненты тензора имеют вид:
τ xx = 2 μ 2 |
¶ x - |
¶V |
y - |
¶ |
z , |
||||||||||
|
|
|
|
|
V |
|
|
|
|
|
V |
|
|||
|
3 |
|
¶x |
|
¶y |
¶z |
|
||||||||
|
|
|
2 |
|
¶V |
y |
|
|
V |
V |
|
||||
τ yy |
|
= |
μ 2 |
|
- |
¶ |
x |
- |
¶ |
z , |
|||||
|
|
|
|
|
|
|
|
||||||||
|
3 |
|
¶y |
|
¶x |
¶z |
|
||||||||
|
|
2 |
|
V |
|
|
|
V |
|
|
|
¶V |
y |
|
|
τ zz |
= |
μ 2 |
¶ z - |
¶ x |
- |
|
, |
||||||||
|
|
¶y |
|||||||||||||
|
3 |
|
¶z |
|
|
|
¶x |
|
|||||||
13
¶V τ = μ x
xy ¶y
τ = μ ¶Vz
xz ¶
x
¶V
τ yz = μ ¶ y
z
+ |
|
¶Vy |
= τ yx |
|
|
¶x |
|
||||
|
|
|
|
|
|
+ |
|
¶Vx |
= τ zx |
|
|
|
|
|
(1.23) |
||
|
¶z |
||||
|
|
|
|
|
|
+ |
|
V |
|
= τ zy |
|
|
¶ z |
|
|||
|
|
¶y |
|
|
|
Применение первого закона термодинамики к жидкости, протекающей через бесконечно малый объем приводит к следующему уравнению энергии:
¶Et + Ñ × EtV |
= |
¶Q - Ñ × q + ρ |
|
|
+ Ñ ×( Ï |
|
) , |
(1.24) |
|||
|
|||||||||||
|
|
|
|
||||||||
f |
|||||||||||
×V |
i, j ×V |
||||||||||
¶t |
¶t |
|
|
|
|
||||||
где Et – полная энергия единицы объема.
Первый член в левой части уравнения есть изменение полной энергии контрольного объема в единицу времени, второй – изменение полной энергии за счет конвекции через поверхность, ограничивающую контрольный объем, в единицу времени. Первый член в правой части – скорость тепловыделения внешних источников, отнесенная к единице объема; второй член – теплопотери за счет теплопроводности через контрольную поверхность в единицу времени; третий член – отнесенная к единице объема работа массовых сил над контрольным объемом; четвертый член – отнесенная к единице объема работа поверхностных сил над контрольным объемом.
Последние два слагаемых правой части можно заменить диссипативной функцией Ф, являющейся тепловым эквивалентом механической
мощности, затрачиваемой на вязкую деформацию жидкости. |
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Ô |
= Ñ ×(τ i, j V× |
|
) - (Ñ ×τ i, j |
)V× |
|
|
|
|
|
|
|
|
|
|
(1.25) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
В декартовой системе координат диссипативная функция принима- |
||||||||||||||||||||||||||||||||||||||||||
ет вид: |
|
|
|
∂V |
|
|
|
|
∂Vy |
|
|
|
∂V |
|
|
|
|
|
∂Vy |
|
|
|
∂V |
|
|
∂V |
|
∂Vy |
|
|||||||||||||||
Ô = μ |
|
|
|
|
|
2 |
|
2 |
|
|
|
|
2 |
|
|
|
2 |
|
2 |
|||||||||||||||||||||||||
2 |
|
|
x |
+ 2 |
|
|
|
|
|
+ 2 |
|
|
|
z |
|
|
+ |
|
|
+ |
|
|
x |
|
|
+ |
z |
+ |
|
+ |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
∂x |
|
∂y |
|
|
|
|
|
|
|
∂x |
|
|
∂y |
|
∂y |
∂z |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.26) |
|||||||||
|
∂V ∂V |
|
2 |
|
2 ∂V |
|
|
∂Vy |
|
|
∂V |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
x |
|
|
|
|
|
z |
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
+ |
|
|
|
+ |
|
|
|
|
|
− |
|
|
|
|
|
|
+ |
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
∂z |
|
|
|
|
|
|
3 |
|
∂x |
∂y |
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Введем величину энтальпии h = ε + p
ρ и получим:
|
|
|
|
14 |
|
|
|
|
||
ρ |
¶h +V |
×Ñh = |
¶p +V |
×Ñp + |
¶Q - Ñ × |
|
+ Ô |
(1.27) |
||
q |
||||||||||
|
¶t |
¶t |
¶t |
|
||||||
Используя закон Фурье для переноса энергии за счет теплопровод- |
||||||||||
ности: |
|
|
|
q = -λÑT , |
|
|
|
|
||
|
|
|
|
|
|
|
(1.28) |
|||
получаем: |
|
|
|
|
|
|
|
|
|
|
ρ |
¶h +V |
×Ñh = |
¶p +V |
×Ñp + |
¶Q + λÑ2T + Ô |
(1.29) |
||||
|
¶t |
¶t |
¶t |
|
||||||
Таким образом полную термодинамическую систему массообмена в газе составляют три уравнения: неразрывности, Навье-Стокса и энергии.
1.7. Виды сплошной среды
Экспериментальные данные показывают, что большинство сред обладает специфическим свойством: отсутствием или малостью касательных
напряжений pSτ, т.е. вектор p S можно считать перпендикулярным любой
площадке взаимодействия dS и равным нормальному напряжению pSn. Среду, обладающую таким свойством называют идеальной жидкостью или идеальным газом. Близки к таковым обычные воздух и вода при малых скоростях.
Понятно, что идеальная жидкость не единственно возможная модель сплошной среды, позволяющая определить компоненты тензора внутренних напряжений. Можно, например, рассматривать его компоненты как функции от деформации частицы: в этом случае среда называется упругой. В частном случае линейности это соотношение приобретает вид закона Гука. Изучением таких сред занимается теория упругости.
Особое место в механике сплошной среды занимает модель вязкой жидкости, предполагающая связь тензора внутренних напряжений с частными производными скорости по координатам. Имеется в виду эффект "трения" слоев вязкой жидкости между собой при наличии разности их поступательных скоростей. В частном случае линейности связь представляется в виде закона Навье-Стокса (или обобщенного закона вязкости Ньютона).
В теории вязкой жидкости m называется коэффициентом внутреннего трения или динамическим коэффициентом вязкости, ν = μρ – кинемати-
ческим коэффициентом вязкости (коэффициентом линейной вязкости),
|
|
15 |
|
|
|
|
|
|
|
|
|
ζ = λ + |
2 |
μ – вторым коэффициентом вязкости (коэффициентом объ- |
|||||||||
|
|||||||||||
3 |
|
|
|
|
|
|
|
|
|
|
|
емной вязкости). Размерность m, l и z в СИ: |
Н× с |
, |
а ν : |
м2 |
|||||||
|
|
2 |
|
|
|
. |
|||||
|
|
|
|
м |
|
|
|
|
|
с |
|
|
|
|
|
|
|
|
|
|
|
|
|
Нетрудно видеть, что упомянутые модели для идеальной и вязкой жидкости вводят еще одну неизвестную – давление p. Т.е. для замыкания системы уравнений движения сплошной среды оказывается необходимым еще одно скалярное соотношение. В этом качестве чаще всего применяются уравнения, представляющие различные гипотезы связи плотности и давления:
ρ = F( p) .
Если такое соотношение можно ввести, то жидкость называется баротропной. Выделяются следующие частные случаи.
1. |
F( p) = const – случай несжимаемой жидкости, или |
dρ |
= 0 . |
|
|||
|
|
dt |
|
2. |
F( p) = Cp , где C – постоянная, – случай изотермического процес- |
||
са.
3.F( p) = Cp1 / n , где C и n – постоянные, – случай политропического
процесса, n называется показателем политропы.
4.F( p) = pμ – уравнение Клапейрона-Менделеева для совершенного
RT
газа, где |
R = 8314,32 |
Дж |
– универсальная газовая постоянная, |
кмоль× град |
μкг
– масса вещества в кг, численно равная молекулярному весу, Tкмоль
–абсолютная температура, которую необходимо задавать еще одним дополнительным соотношением.
1.8. Свойства воздушной среды
Химический состав воздуха: азот - 78% , кислород – 21%, углекислый газ, инертные газы, водяной пар.
Плотность ρ = 353.4 . Если движение воздушного потока v << a, где a =
T
300 м/с – скорость звука в воздухе, то зависимостью плотности от давления можно пренебречь. В системах вентиляции плотность от давления мало зависит. T – температура в кельвинах.
16
При t = 0 °C ρнорм.усл = 1.29 кг / м3
При t = 20 °C ρстандарт.усл = 1.2 кг / м3
Вязкость газовой среды характеризуется коэффициентами динамической и кинематической вязкости.
Коэффициент динамической вязкости, в отличие от плотности, зависит только от температуры газа.
μвозд = 1.71×10−5 + 4.93×10−8T (Па×с) – динамическая вязкость;
ν возд = μρ - кинематическая вязкость.
Характер движения газового потока определяется отношением инерционных и вязкостных сил, которое характеризуется числом Рейнольдса
Re = vd ρ
μ
где v - средняя скорость течения, d – характерный размер течения. При определенных (критических) значениях числа Рейнольдса течение газа становится турбулентным, что характерно для систем промышленной вентиляции и газоочистки.
2.Приближенные методы описания движения воздушных потоков
2.1.Воздушные потоки вблизи всасывающих отверстий
Впомещении или открытом пространстве, из которого принудительно отсасывается воздух, возникает своеобразный воздушный поток, называемый стоком.
Причиной образования стока является разность между атмосферным давлением и разрежением воздуха в плоскости всасывающего отверстия, создаваемым работой отсасывающего устройства. Под действием разности давлений окружающий воздух со всех сторон устремляется к всасывающему отверстию, при этом по мере движения давление воздуха уменьшается, а скорость его движения и ускорение увеличиваются.
Секундный объем отсасывамого воздуха служит количественной мерой стока и носит название расхода.
Если сток воздуха к всасывающему отверстию развивается от твердых стенок, то ничто, не тормозит его свободного течения, и каждая частица воздуха в своем движении сближается с соседними, но не обгоняет их и не отстает от них. В воздухе не возникает касательных напряжений, не проявляются трение и вязкость. Течение можно считать потенциальным.
17
Точечный сток.
|
|
Точечный сток – пространственный воздушный |
||
|
|
поток, устремленный к одной точке, где он поглощает- |
||
|
|
ся. Точку исчезновения потока называют полюсом. |
||
|
|
Линии тока точечного стока – прямые, сходящиеся в |
||
|
|
|||
|
|
полюсе, а точки с одинаковым значением скорости |
||
|
|
расположены на одинаковом расстоянии от него. В |
||
|
|
случае свободного стока, которому не препятствуют |
||
|
|
твердые поверхности, точки с одинаковым значением |
||
Рис. 2.1. Свобод- |
скорости принадлежат сферическим поверхностям, |
|||
ный точечный |
скорость движения воздуха в произвольной точке вы- |
|||
сток |
ражается уравнением неразрывности |
|||
|
|
v = |
L0 |
|
|
|
|
(2.1) |
|
|
|
4πr 2 |
||
L0 - секундный расход воздуха; r – расстояние между произвольной точ-
кой пространства и полюсом стока; 4πr2 - поверхность сферы радиусом r. Примером свободного точечного стока может служить поток воз-
духа в помещении, вызванный отсосом через тонкую трубку.
Увеличить скорость движения воздуха на заданном расстоянии от полюса стока можно без увеличения расхода, если ограничить область течения твердыми непроницаемыми поверхностями. Для точечного стока из полупространства, ограниченного плоской стенкой, скорость движения по сравнению со скоростью движения свободного точечного стока удваивается:
v = |
L0 |
(2.2) |
2πr 2 |
В знаменателе стоит площадь полусферы. Такой сток воздуха называют полуограниченным. Примером полуограниченного стока является течение, образованное отсосом воздуха через небольшое отверстие в стене помещения.
Дальнейшее увеличение скорости течения может быть достигнуто путем еще большего сужения области подтекания воздуха к полюсу. Обобщающая формула для определения скорости точечного стока
v = |
L0 |
(2.3) |
Ψr2 |
где Ψ - телесный угол, под которым из полюса видна часть открытого пространства, откуда подтекает воздух.
18
Таблица 2.1.
Значения телесного угла для некоторых условий ограничения точечного стока
№ |
Поверхности, |
ограничи- |
Ψ |
|
вающие точечный сток |
|
|
1 |
Отсутствуют |
|
4 π |
2 |
Плоская стенка |
|
2 π |
3 |
Грани прямого двугранного |
π |
|
|
угла |
|
|
4 |
Грани прямого трехгранного |
π 2 |
|
|
угла |
|
|
5 |
Боковая поверхность конуса |
2π (1 − cosψ 2) |
|
|
с углом Ψ при вершине |
|
|
Чтобы найти составляющие скорости подтекания воздуха, следует выбрать координатную систему. Воспользуемся пространственной декартовой системой координат, совместив ее начало с полюсом.
Тогда, составляющие скорости будут вычисляться по следующим формулам:
vx = v |
x |
; |
v y = v |
y |
; |
vz = v |
z |
(2.4) |
|
|
r |
||||||
|
r |
|
r |
|
|
|||
где r = 
x2 + y2 + z2 - расстояние между произвольной точкой и началом
координат.
Линейный сток
Свободный линейный сток – это пространственный воздушный поток, устремленный к бесконечно длинной прямой линии, в которой поток поглощается. Воздух со всех сторон устремляется к линии поглощения потока, а точки с одинаковым значением скорости образуют цилиндрические поверхности.
Скорость движения воздуха в любой точке свободного линейного стока связана с расстоянием до поглощающей поток линии уравнением неразрывности:
v = |
L |
(2.5) |
|
2πlr |
|
||
L/l - удельный секундный расход воздуха; r – |
расстояние от произвольной |
||
точки до линии поглощения потока. |
|
||
Если с линией стока совместить оcь y, то составляющие скорости свободного линейного стока будут равны:
19
vx = − Lx 2 2πlr
(2.6)
vy = − Lz 2 2πlr
где x, z – координаты точки, в которой нужно найти значение скорости
воздуха; r = 
x2 + z2 .
Для линейного стока воздуха, ограниченного с одной стороны плоской стенкой, скорость течения воздуха по сравнению со свободным стоком удваивается:
v = |
L |
(2.7) |
|
πlr |
|||
|
|
Если линия поглощения расположена в ребре двугранного угла, грани которого составляют плоский угол ϕ , то скорость течения воздуха
на заданном расстоянии определяется формулой
v = |
L |
(2.8) |
|
ϕlr |
|||
|
|
Рассмотрим взаимодействие двух ограниченных с одной стороны линейных стоков одинакового удельного расхода, расположенных в плоскости x=0 параллельно оси y на расстоянии 2b один от другого.
Начало координат поместим в точку, разделяющую расстояние между полюсами стоков пополам.
Тогда, для произвольной точки полупространства, составляющие скорости, образованного независимым действием одного линейного стока равны:
v1x
v1z
|
L1x |
|
= − |
πl(x2 + (z − b)2 ) |
(2.9) |
L(z − b)
=− πl(x21+ (z − b)2 )
Составляющие скорости потока, образованного независимым действием другого стока составляют:
v2x = − |
|
L2 x |
|
(2.10) |
||
πl(x2 |
+ (z + b)2 ) |
|||||
|
|
|
||||
v2 z = − |
|
L1 |
(z + b) |
|
||
2πl(x2 |
+ (z + b)2 ) |
|
||||
20
Результирующие составляющие скорости потока, образованного взаимодействием двух одинаковых линейных свободных стоков, является суммой найденных составляющих:
vx = v1x + v2x
vz = v1z + v2z
Вплоскости симметрии vz = 0 . Нулевая скорость свидетельствует
об отсутствии потока; воздуха не проводится сквозь плоскость, и она может быть заменена непроницаемой поверхностью.
Таким образом, формулы отражают результат взаимодействия двух одинаковых свободных линейных стоков или результат взаимодействия одного линейного стока с непроницаемой поверхностью.
Сток воздуха к прямоугольным отверстиям
Рассмотрим прежде всего сток воздуха к бесконечно длинному отверстию конечной ширины, устроенному в плоской стене. Поместим начало координат в центр отверстия так, чтобы ось x была направлена навстречу потку воздуха. Ширина 2B, а секундный удельный расход воздуха L/l. Сток воздуха к отверстию представим как результат взаимодействия бесконечно большого числа бесконечно малых линейных стоков воздуха. Разделим мысленно всасывающее отверстие на бесконечно тонкие полоски шириной db.
Удельный расход воздуха, поглощаемый каждой из таких полосок, будет элементарным и при условии равномерного всасывания определится сле-
дующим соотношением dL = L db l 2lB
Выделим одну из элементарных полосок на расстоянии b от начала координат и допустим, что воздух из полупространства стекает только к ней одной. Образованный таким образом элементарный сток подчиняется закономерностям линейного стока.
Составляющие скорости воздуха, образуемой таким элементарным стоком:
dvx = dLx πlr 2
( ) (2.12)
dvz = dL z − b πlr2
Подставим в эти соотношения значения dL и r :
