Проинтегрировав (2.116) с учетом (2.114), (2.115), (2.117),
находим
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
arccos |
|
|
|
|
|
|
|
0,303u '3 |
|
|
2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
q1 |
q2 2 |
|
||||||||||
f |
4 1 |
q 2 |
2 e |
|
||||||||||||||||
|
|
q2 |
|
l |
|
|
|
|
1 |
|
||||||||||
|
|
|
|
1 |
ж t |
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
2 |
|
|
4 |
|
|
2 |
(n 1) |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
1 |
e |
|
q1 |
e |
q1 |
... |
e |
q1 |
P . |
(2.118) |
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
Сумма экспонентных функций в круглых скобках выражения (2.118) есть сумма членов геометрической прогрессии со знаменателем, поэтому запишем
|
|
|
|
|
|
|
|
|
2 |
|
|
4 |
|
|
|
|
2 (n |
1) |
|
|
|
2 |
n |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
|
q1 |
1 |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
1 |
e |
q1 |
e |
q1 |
... e |
|
q1 |
|
|
|
. |
(2.119) |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
|
q1 |
1 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
Кроме того, отметим, |
поскольку |
|
по физическому |
смыслу |
||||||||||||||||||||||||||||||
|
It |
>> |
|
n |
, т. е. n |
>> 1, а колебания затухающие, величиной P можно |
|||||||||||||||||||||||||||||||
|
u |
|
q1 |
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
пренебречь, приняв одновременно n |
|
nn . |
В результате из уравне- |
||||||||||||||||||||||||||||||||||
ний (2.118) и (2.119) получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
1,21u'3 |
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
f |
1 |
|
|
|
|
exp |
arccos |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
q12 |
|
q12 2 |
|
q1 |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
жlt |
|
|
|
q12 |
2 |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
exp |
|
2 lt |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u' |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.120) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
exp |
2 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где q1 − частота собственных колебаний ПКФ.
86
Вывод уравнения (2.120) приведѐн в сокращѐнном виде. Желающие получить более подробную информацию о решении представленной задачи могут обратиться к работе [9].
Частота собственных колебаний ПКФ
|
2,4A 2 |
|
0,5 |
|
|
q1 |
2 |
. |
(2.121) |
||
R |
|
|
|||
|
|
|
|
||
Скорость распространения волн по поверхности пузыря A определяется ее упругими свойствами, зависящими от поверхностного натяжения и плотности жидкости. Тогда из условий размерности следует записать
|
|
|
c |
0,5 |
|
|||
|
A |
|
|
, |
(2.122) |
|||
|
|
|
|
|||||
|
|
R |
|
|||||
|
|
|
ж |
|
||||
где с − коэффициент пропорциональности. |
|
|||||||
Согласно рис. 2.10, R sin |
dп |
. Подставив значение R в ра- |
||||||
|
||||||||
|
|
|
|
2 |
|
|
|
|
венство (2.121), запишем |
|
|
|
|
|
|
|
|
|
|
2 |
0,5 |
|
||||
A |
|
|
|
|
, |
(2.123) |
||
|
|
|
|
|||||
|
|
|
|
|||||
|
|
c1dп ж |
|
|||||
где c1 sin / c − коэффициент, определяемый экспериментально. Строго говоря, эта величина переменная, так как угол
зави-
сит от силы воздействия турбулентных пульсаций на поверхность пузыря, поэтому при решении каждой конкретной задачи должен быть определен отдельно.
Согласно работам [12, 13], пульсационную скорость можно выразить в виде зависимости
|
ELt |
1/ 3 |
|
|
u |
|
|
, |
(2.124) |
|
ж |
|||
|
|
|
|
|
где Lt (в дисперсных системах) |
является размерным |
значением |
||
масштаба турбулентности. |
|
|
|
|
87
В вопросе о выборе масштаба турбулентности единого мнения нет. Нами при решении данной задачи, для расчѐта масштаба турбулентности, принято уравнение, представленное в работе [9] (там же дана ссылка на первоисточник), так как оно в конечном итоге обеспечило наилучшую сходимость с опытными данными:
L |
dп (1 |
г ) |
. |
(2.125) |
|
|
|||
t |
6 |
|
|
|
|
г |
|
||
Подставляя уравнение (2.125) в (2.124), получим зависимость для определения пульсационной скорости:
|
Edп (1 |
г ) |
1,3 |
(2.126) |
|
u |
|
|
. |
||
6 г |
ж |
||||
|
|
|
Силу сопротивления среды колебанию ПКФ и частоту δ определим по данным работы [13]:
|
|
u' |
ж dy |
ds . |
(2.127) |
|||
Fтр ~ |
ж |
|
ds |
|
|
|
||
d |
d dt |
|||||||
В уравнении (2.100) слагаемое, характеризующее силу сопротивления, представим в виде пропорциональности
2 |
y |
~ |
|
Fтр |
. |
(2.128) |
|
t |
R |
ж ds |
|||||
|
|
|
|
В таком случае из (2.127) и (2.128) следует
c2 |
ж |
, |
(2.129) |
|
2Rd |
||||
|
|
|
где с2 − коэффициент пропорциональности.
Подставив уравнения (2.121)−(2.129) в равенство (2.120) и пренебрегая для маловязких жидкостей влиянием сил трения, получим выражение для определения осреднѐнной величины межфазной турбулентности в виде зависимости
88
|
|
|
0,054c3/ 2 E 2 / 3d |
13/ 6 |
(1 |
|
2 / 3 |
|
||
|
|
|
|
|
|
|||||
|
|
п |
г |
|
|
|||||
f |
|
1 |
|
|
|
. |
(2.130) |
|||
1/ 6 |
1/ 2 |
|
|
|
г |
|||||
|
|
|
|
|
|
|
||||
|
|
|
ж |
ж |
|
|
|
|
||
2.6.6. Определение диаметра пузырей
Значение диаметра пузыря, образующегося после дробления, найдѐм из следующих соображений. Предположим, что дробление пузырей происходит под воздействием динамического напора турбулентных пульсаций. При этом величина у в уравнении (2.102) достигает максимального значения, равного диаметру на момент дробления пузыря.
Величину ymax определим, исследовав уравнение (2.102) на максимум. Приравняв к нулю выражение (2.105), получим значе-
ние tm , при котором y |
ymax : |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
q t |
arctg |
q1 |
. |
|
|
|
|
(2.131) |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
1 m |
|
|
|
|
|
|
|
|
|
Функция I0 |
2,4 r |
|
имеет максимум при r |
0, I0 (0) |
1. Под- |
|||||||||
|
|
|
||||||||||||
|
R |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ставив эти значения и q1tm из (2.131) в равенство (2.105), получим |
||||||||||||||
|
|
|
|
|
1, 6u |
' |
|
|
|
δ |
|
q |
|
|
y |
max |
d |
п |
|
|
|
exp |
|
arctg |
1 |
. |
(2.132) |
||
|
(q12 δ2 )0,5 |
|
|
|||||||||||
|
|
|
|
|
|
q1 |
|
δ |
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||
Анализ уравнения (2.132) показал, что изменение коэффициента с2 в широком интервале не влияет на значение dп . По данным ра-
боты [13] рекомендуется с2 = 30. Исходя из этого влиянием экспоненты можно пренебречь; равенство (2.132) с учетом (2.121), (2.126) и (2.129) примет вид
|
11,6 |
0,6 |
|
|
0,4 |
|
||
dп |
|
г |
. |
(2.133) |
||||
|
|
|
|
|
||||
1,8 |
E |
0,4 0,2 |
1 г |
|||||
|
c1 |
|
ж |
|
|
|||
В уравнении (2.133) значение коэффициента с1 определялось на основе экспериментальных данных по измерению среднего диа-
89
метра пузырей в аппаратах с различными системами аэрации [9]. В барботажных аппаратах колонного типа за основу определения с1 в уравнении (2.133) были взяты данные из работ [14, 15]. Этот выбор диктовался тем, что авторы K. Akita и F. Ioshida провели исследования по измерению размеров пузырей в барботажной колонне методом фотографирования в различных жидкостях в широком диапазоне изменения их физических свойств. Кроме того, одновременно измеря-
лось истинное объемное газосодержание |
г , для вычисления которо- |
|||||||
го предложено эмпирическое критериальное уравнение |
|
|||||||
|
г |
|
0,2Bo |
0,125 |
Ga |
0,0833 |
Fr . |
(2.134) |
(1 |
г )4 |
|
|
|||||
|
|
|
|
|
|
|||
Этими же авторами было предложено уравнение для расчѐта диаметра пузырей:
|
gd 2 |
0,5 gd 3 0,12 |
w |
|
0,12 |
||
dп 26 |
a ж |
|
a |
|
г |
|
. (2.135) |
|
2 |
|
(gda ) |
0,5 |
|||
|
|
|
|
||||
|
|
|
ж |
|
|
||
Вычисленные по уравнению (2.135) значения диаметра пузырей подставляли в уравнение (2.133), откуда затем определяли коэффициент c1 . Входящую в уравнение (2.133) величину Е найдем из
уравнения (5.43) – см. разд. 5.
В результате обработки экспериментальных данных указанных работ было установлено, что коэффициент c1 не является величиной
постоянной и несколько зависит от размеров аппарата и физических свойств жидкости, слабо – от газосодержания. Для согласования опытных данных и уравнения (2.133) потребовалась необходимость представления c1 в виде зависимости
|
da g1/ 3 |
0,18 |
|
c1 0,28 |
. |
||
2 / 3 |
|||
|
|
Отсутствие в уравнении (2.133) вязкости жидкости и появление еѐ в полученном равенстве вполне понятно, так как при его выводе мы ею пренебрегли. Однако причиной появления в нѐм значи-
90