Diplom

advertisement
Федеральное агентство по образованию
ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ (ТГУ)
Физико-Технический факультет
Кафедра прикладной газодинамики и горения
Допустить к защите в ГАК
Зав. каф. прикладной газовой
динамики и горения
д-р физ.-мат. наук, профессор
_______________Г.Р.Шрагер
‘‘_____’’_____________200_г.
Бакалаврская работа
ВЛИЯНИЕ НЕСТАЦИОНАРНЫХ ЭФФЕКТОВ НА ДИНАМИКУ
ВСПЛЫТИЯ ПУЗЫРЬКА
Дабаева Дарима Бимбаевна
Руководитель
д.ф.-м.н., профессор
_________ В.А.Архипов
Автор работы
__________ Д.Б.Дабаева
Томск 2012
ABSTRACT
In this work we investigate the effect of non-stationary and hereditary forces
on the dynamics of spherical bubble rising. Calculations were carried out to find
stationary rising rate, taking the present mass. Due to these results an experimental
setup was developed. The results have practical applications in several important
processes related to ozone treatment of rivers, bubbling and the solution of some
environmental problems.
2
ОГЛАВЛЕНИЕ
Список обозначений…………………………………………………..………..4
Введение………………………………………………………………..……….5
1 Динамика всплытия пузырька без учета присоединённой массы…..…….7
Постановка задачи…………………………………………..…..7
Оценка стоксовского режима………………………………..… 9
1.2.1 Расчет Dmax пузырька в различных растворах…………….10
1.2.2 Расчет u0 для пузырьков с различными диаметрами……..10
1.3. Определение динамики всплытия пузырька…………………..13
1.3.1 Зависимость скорости всплытия от времени…………… .13
1.3.2 Нахождение зависимостей скорости всплытия от времени в
различных растворах……………………………………………...16
1.3.3 Зависимость скорости всплытия от пройденного расстояния
……………………………………………………………………...18
1.3.4 Нахождение зависимостей скорости всплытия от расстояния
в различных растворах…………………………………………....20
2 Динамика всплытия пузырька с учетом присоединенной массы………….22
2.1. Постановка задачи …………………………………………………..22
2.2. Оценка для стоксовского режима…………………………………..25
2.3. Определение динамики всплытия пузырька……………………….26
2.3.1 Зависимость скорости всплытия от времени……………….26
2.3.2 Нахождение зависимостей скорости всплытия от времени в
различных растворах ……………………………………………...26
2.3.3 Зависимость скорости всплытия от пройденного расстояния
……………………………………………………………………….29
2.3.4Нахождение зависимостей скорости всплытия от расстояния
в различных растворах ……………………………………………29
3 Экспериментальная установка ля исследования движения пузырька……..32
3.1 Описание установки …………………………………………………32
3.2 Описание видеокамеры ………………………………………………33
Заключение……………………………………………………………………….36
Список использованной литературы…………………………………………...37
1.1.
1.2.
3
СПИСОК ОБОЗНАЧЕНИЙ
A
18 ж
переменная;
г D2
B(
ж
 1) g переменная;
г
CD- безразмерный коэффициент сопротивления;
D-диаметр пузырька;
Dmax-максимальный диаметр пузырька для которого выполняется условие
стоксовского режима;
Fa- сила Архимеда;
Fc-сила сопротивления;
Fc1- сила вязкого сопротивления;
Fc2- сила, связанная с присоединенной массой;
Fc3- сила Бассэ;
Fg-сила тяжести;
Fi- силы, действующие на пузырек;
g-ускорение свободного падения;
m- масса пузырька;
mпр- приведенная масса;
R-радиус пузырька;
S-площадь миделева сечения;
t-время;
to-безразмерный параметр, равный u0/g;
u-скорость всплытия пузырька;
u0-скорость стационарного всплытия;
V-объем пузырька;
x-пространственная координата;
x0=u0·t0- масштаб расстояния;
z=Au-B –переменная;
y=u/u0-переменная;
µ-динамическая вязкость жидкости;
ξ=x/x0- переменная для замены;
π-математическая константа;
ρг- плотность газа(пузырька );
ρж-плотность жидкости, в которой находиться пузырек;
τ=t/t0-переменная для замены;
Критерии подобия:
Re-число Рейнольдса.
4
Введение
В
данной
работе
исследуется
влияние
нестационарных
и
наследственных сил на динамику всплытия пузырька.
Процессы движения пузырьков в жидкости в поле силы тяжести
играют важную роль в целом ряде технологических процессов и в
технических
устройствах.
В
качестве
примеров
можно
привести
барботажные процессы, которые используются при озонировании воды,
получении газированных напитков, а так же в процессах
обогащения
полезных ископаемых методом флотации.
Большую роль пузырьковые течения играют при возникновении
кавитации, которое приводит к разрушению грибных винтов теплоходов, и
других водных транспортных средств.
Все указанные процессы зависят от концентрации, размеров и скорости
всплытия пузырьков.
В
большинстве
литературы
динамику
всплытия
пузырька
рассматривают без учета нестационарных связей, однако данный эффект
играет важную роль в ряде многих процессов, отмеченных ранее. В связи с
этим целью работы являлось теоретическая оценка режимов всплытия с
учетом нестационарных сил и разработана экспериментальная установка для
изучения данного вопроса.
Данная работа состоит из трех разделов. В первом рассмотрена
динамика
всплытия
пузырька
в
стационарном
теоретические расчеты скорости пузырьков
режиме.
Приведены
в различных растворах. При
движении пузырьков в режиме ускорения на них действуют дополнительные
силы: сила, приведенной массы, связанная с присоединенной массой и сила
Бассэ.
Во втором разделе рассматривается с учетом присоединенной массы.
5
Результаты, полученные в первых двух разделах, используются при
разработке специальной установки. В заключительном третьем разделе
описывается экспериментальная установка.
6
1 Динамика Всплытия пузырька без учета присоединенной массы
1.1Постановка задачи
Для
определения
динамики
всплытия
сферического
пузырька
диаметром D в жидкости, то есть для нахождения зависимости скорости
всплытия пузырька от времени u(t), рассмотрим уравнение движения
пузырька в форме второго закона Ньютона:
m
du
  Fi
dt
i
(1.1)
Где m=(πD3/6)ρг-масса пузырька диаметром Dи плотностью газа ρг
u - вектор скорости пузырька
Fi -силы, действующие на пузырек.
В проекции на ось z, направленную вертикально вверх (противоположно
направлению вектору ускорения свободного падения g),действующие на
пузырек силы можно представить в виде:
сила тяжести:
сила Архимеда :
Fg  mg  
Fa  Vg 
Где ρж- плотность жидкости
V-объем пузырька
г -плотность газа
7
 D3
6
 D3
6
г g
ж g
(1.2)
(1.3)
сила сопротивления:

Fc  CD S u 2
2
(1.4)
2
Где S=πD /4- площадь миделевого сечения
СD-безразмерный коэффициент сопротивления.
Для
стоксовского
режима
движения
пузырька(Re<=1)
коэффициент
сопротивления равен
СD=24/Re
(1.5)
Где число Рейнольдса определяется зависимостью Re  жuD /  , где  динамическая вязкость жидкости
подставляя (1.5) в (1.4) получим :
 D 2  жu 2 24ж
Fc  
.
.
 3 uDж
4
2  ж uD
Подставляя действующие силы Fg (1.2), Fa
движения (1.1):
(1.6)
(1.3) и Fc (1.6) в уравнение
 D3
du
 D3
 D3
ж .  
. г g 
ж g  3 uDж (1.7)
6
dt
6
6
Преобразуем уравнение (1.7), разделив обе части на
 D3
6
ж
du
18
 g(
 1) 
u
2
dt
г
г D
Распишем Для случая стационарного всплытия, то есть
8
. г
(1.8)
du
0
dt
ж
 1) D 2  г
г
( ж  г ) D2
u0 
g
g
18ж
18ж
(
Из
(1.9)
видно,
что
скорость
стационарного
всплытия
(1.9)
пузырька
пропорциональна квадрату от диаметра.
1.2 Оценка стоксовского режима
Найдем максимальный диаметр пузырька Dmax, для которого выполняется
условие
стоксовского
режима
(Re<<1).
Подставим
уравнение
(1.9)
выражение для числа Рейнольдса, получим:
 жuD  ж D (  ж   г ) D 2
 ж (  ж   г ) D3
Re 

.
g
g
2
ж
ж
18ж
18ж
(1.10)
Из уравнения (1.10), при Re=1 определим Dmax:
Dmax 
3
18ж 2
ж ( ж  г ) g
9
(1.11)
1.2.1 Рассчет Dmax пузырька в различных растворах
Имеем следующие параметры расчёта:
g  9.80665 м / с 2
. ж  103 кг / м3
 г  1.205кг / м3
А) Для раствора воды с вязкостью µ=10-3 Па·с
Подставив условия в (10) получим выражение для Dmax.
Dmax  1.22 103 м  1.22 мм
Найдем значение скорости для соответствующего максимального диаметра
Dmax. Для этого подставим в уравнение (8) значение D.
Подставив, мы получаем численное значение u0=0.77м/с.
Б) Для раствора с вязкостью µ= 10-2 Па·с
Dmax=0.005685м=5,685мм.
Для этого значения находим u0 =1.7м/с
В) Для раствора с вязкостью µ= 10-1 Па·с
Dmax=0.026м=26мм
Для этого значения находим u0=3.67м/с
1.2.2 Расчет u0 для пузырьков с различными диаметрами
А) Для раствора воды с вязкостью µ=10-3 Па·с
Найдем численные значения u0 для различных D. Составим таблицу 1.1.
10
Таблица 1.1
D,
мм
0
0.4
0.8
1
1.2
U0,
м/с
0
0.08705 0.3482 0.54406 0.77
1.5
2
3
4
1.22
2.176
4.8
8
Представим наглядно, на рис.1.1
Мы видим что график зависимости представляет собой, ветвь
параболы, направленную верх . То есть с увеличением диаметра скорость
увеличивается
квадратично.
Рисунок 1.1
Зависимость скорости стационарного всплытия u0 от диаметра пузырька D
для раствора с вязкостью µ= 10-3 Па·с
Б) Для раствора с вязкостью µ= 10-2 Па·с.(Таблица 1.2 и рисунок 1.2)
11
Таблица 1.2
D, мм 0
0.8
1
2
3
4
5.6
7
8
U0,
м/с
0.034
0.054
0.217
0.489
0.87
1.706
2.66
3.48
0
Рисунок 1.2
Зависимость скорости стационарного всплытия u0 от диаметра пузырька D
для раствора с вязкостью µ= 10-2 Па·с
В) Для раствора с вязкостью µ= 10-1 Па·с (Таблица 1.3 и рисунок 1.3)
Таблица1.3
D, мм 0
2
4
5
10
15
26
30
U0, м/с
0.021
0.087
0.136
0.54
1.224
3.67
4.89
0
12
Рисунок 1.3
Зависимость скорости стационарного всплытия u0 от диаметра пузырька D
для раствора с вязкостью µ= 10-1 Па·с
1.3.Определение динамики всплытия пузырька
1.3.1определение скорости всплытия пузырька от времени
Для этого решим уравнение (1.8)

du
18
 g ( ж  1) 
u
2
dt
г
г D
B(
Примем
ж
 1) g
г
и
A
18 ж
г D2
Тогда уравнение можно переписать в виде:
13
du
 B  Au
dt
du
 ( Au  B )
dt
Au  B  z
dz  Adu
du 
dz
A
1 dz
 z
A dt
dz
  Adt
z
Проинтегрируем обе части, в пределах от z до z0 , и по времени от 0 до t:
ln z |zz0   At |t0
ln z  ln z0   At
Воспользуемся свойством логарифмов для преобразования уравнения.
z
ln( )   At
z0
z
 exp( At )
z0
Окончательно, мы получаем уравнение для зависимости z(t):
z  z0 exp( At )
Из начальных условий при t=0 u=0, найдем значение z0.
14
(1.12)
du
 Au  B   z
dt
z0 : 0  0  B
z0  B
Вернемся к замене, подставив значение z.
Au  B   B exp( At )
(1.13)
Преобразуем,
u
B
1  exp   At  
A
(1.14)
Подставим значения A и B в уравнение (1.12) получаем,
ж
 1) g

 18ж  
г
u
1  exp  t

2 
18ж 

D
г


2
г D
(
(1.15)
Покажем что при t→∞ u=u0.При t→∞ exp→0, значит u=B\A. Преобразуем
В\А.
ж
 1) g
г
(  ж   г )  гг D 2 (  ж   г ) gD 2
B


g

18ж
A
г
18ж
18ж
г D2
(
А выражение
(  ж   г ) gD 2
18ж
и есть u0.
15
Значит
формулу
(1.13)

 18 ж
u  u0 1  exp  t
2

D
г


можно
представить
в
виде

 .

1.3.2 Нахождение зависимостей скорости всплытия от времени в
различных растворах
А) Для раствора воды с вязкостью µ=10-3 Па·с
С помощью приведенных ранее таблицы 1.1. и рис.1.1 найдем u=u(t)
Значит, u=u0(1-exp(t*(-18*10-3)/1.275*D2))
Графики зависимостей u(t) с различными диаметрами, представлен на рис.1.4
Рисунок 1.4
Зависимость скорости всплытия u от времени t
для раствора с вязкостью µ= 10-3 Па·с
16
Б) Для раствора с вязкостью µ= 10-2 Па·с
С помощью таблицы 1.2 и рисунка 1.2 u(t) будет
u=u0(1-exp(t*(-18*10-2)/1.275*D2))
График зависимостей u(t) представлен на рисунке 1.5
Рисунок 1.5
Зависимость скорости всплытия u от времени t
для раствора с вязкостью µ= 10-2 Па·с
В) Для раствора с вязкостью µ= 10-1 Па·с
С помощью таблицы 1.3 и рисунка 1.3 u(t) будет
u=u0(1-exp(t*(-18*10-1)/1.275*D2))
График зависимостей u(t) представлен на рисунке 1.6
17
Рисунок 1.6
Зависимость скорости всплытия u от времени t
для раствора с вязкостью µ= 10-1 Па·с
1.3.3 1определение скорости всплытия пузырька от пройденного
расстояния
Найдем теперь зависимость скорости от пройденного пузырьком расстояния.
Для этого решим уравнение (1.8) для х. Рассмотрим в неявном виде. Для
этого введем безразмерные параметры:
t0 
u0
g
  tt
0
x0  u0  t0
  xx
0
y u
u0
(1.14)
Тогда исходное наше уравнение примет вид
18
t0 du
u
u
 1
u0 dx
u0
(1.15)
Перейдя к безразмерным переменным τ, ξ, y получим
y
dy
 1 y
d
(1.16)
Проинтегрировав дифференциальное уравнение (3’) с нулевыми начальными
условиями получим
  ( y  ln(1  y))
(1.17)
Теперь рассмотрим решение, так что бы зависимость получилась в явном
виде.
du
 B  Au
dt
du dx
du
u
 B  Au
dx dt
dx
udu  ( B  Au )dx
Проинтегрировав (1.17) с учетом начальных условий x=0, y=0 получим:
x
B
A2
A 
A
u

Ln
(1

u)
 B
B 
ж
 1) g
г
18 ж
A
г D2
B(
Где
19
(1.18)
1.3.4 Нахождение зависимостей скорости всплытия от расстояния в
различных растворах
А) Для раствора воды с вязкостью µ=10-3 Па·с
Подставляя в (1.18) значения u от 0 до значения B/A и значения диаметров
пузырьков из таблицы 1.1.получим графики зависимости u(x).
Рисунок 1.7
Зависимость скорости всплытия u от расстояния х
для раствора с вязкостью µ= 10-3 Па·с
20
Б) Для раствора с вязкостью µ= 10-2
Па·с
Рисунок 1.8
Зависимость скорости всплытия u от расстояния х
для раствора с вязкостью µ= 10-2 Па·с
В) Для раствора с вязкостью µ= 10-1 Па·с
Рисунок 1.9
Зависимость скорости всплытия u от расстояния х
для раствора с вязкостью µ= 10-1 Па·с
21
2 Динамика всплытия пузырька с учетом присоединенной массы.
2.1 Постановка задачи
Процессы движения частиц дисперсной фазы (твердых частиц, капель
и пузырьков) в поле силы тяжести играют важную роль в природе, в
различных технологических процессах и в проблемах экологии. В качестве
примеров можно привести закономерности образования атмосферных
осадков, процессы седиментации и барботажа, распространения облака
жидко-капельных аэрозолей при аварийном сбросе авиационного топлива и
токсичных компонентов жидких ракетных топлив при отделении ступеней
ракет-носителей [8].
В большинстве работ по гравитационному осаждению частиц
рассматривается уравнение движения, учитывающее силу тяжести, силу
Архимеда и силу сопротивления [9]:
du
u 2
m
 mg  Vg  C D S
,
dt
2
где
g
(2.1)
– масса и объем частицы; u – скорость гравитационного осаждения;
m, V
– ускорение свободного падения;

безразмерный коэффициент сопротивления;
– плотность жидкости;
S
CD
–
– площадь миделева сечения
частицы.
Для
стоксовского
режима
движения
(Re<0.1)
коэффициент
сопротивления сферической частицы равен
CD 
где
D
24 24 

,
Re uD
– диаметр частицы;

(2.2)
– коэффициент динамической вязкости
жидкости.
С учетом выражения для коэффициента сопротивления (2.2) уравнение
движения частицы (2.1) примет вид:
22


du
 g 1 
 p
dt

где
 18u

,
  D2

p
(2.3)
– плотность материала частицы.
Аналогичное уравнение можно записать и для случая всплытия
p
пузырька в жидкости.
В ряде практически важных задач необходимо учитывать влияние
нестационарных и «наследственных» сил, действующих на частицу [3-4].
При нестационарном движении тела как в идеальной, так и в вязкой
жидкости на него со стороны жидкости действует сила, связанная с
присоединенной массой [5]. Для шара сила сопротивления направлена
противоположно вектору ускорения шара du / dt и равен
Fпр  
где
V du
,
2 dt
– объем частицы.
Уравнение движения шара с учетом
(2.4)
V
m
Fпр
(4) примет вид
du
V du
F
dt
2 dt
или
  du

 F,
  p  V
2
dt


где
F
(2.5)
– действующие на шар внешние силы (сила тяжести, сила Архимеда и
сила сопротивления).
Коэффициент при du / dt в (5) можно рассматривать как некоторую
эффективную массу шара, которая складывается из массы самого шара
m   pV
и присоединенной массы
mпр  V / 2 , которая равна
половине массы жидкости, вытесненной шаром. Особенно значителен вклад
Fпр при движении пузырька.
23
На частицу, движущуюся с переменной скоростью в вязкой жидкости
действует еще одна дополнительная сила, зависящая от предыстории
движения («наследственная» сила Бассэ [5]). Для жидкости большой
вязкости она представляется через линейный интегральный оператор или
определяется численным решением интегро-дифференциального уравнения,
полученного Бассэ [3].
Л.Д. Ландау получил аналитическое решение для суммарной силы
сопротивления [5]:
 1 3
6
Fc  2 R 3 g  
t
R
 3 R 2
где

t 
,
 
– коэффициент кинематической вязкости жидкости;
(2.6)
R
– радиус
шара.
Подставляя в (2.6) диаметр шара
динамической вязкости
D  2R и коэффициент
   , получим для силы сопротивления
следующую формулу:
Fc  Fc1  Fc 2  Fc3 ,
где
(2.7)
Fc1  3Dgt – сила вязкого сопротивления;
Fc 2  Vg / 2 – сила, связанная с присоединенной массой;
Fc3  3gD 2 t – сила Бассэ.
Для определения динамики всплытия сферического пузырька
диаметром D в жидкости, то есть для нахождения зависимости скорости
всплытия пузырька от времени u(t) с учетом присоединенной массы
рассмотрим уравнение движения пузырька в форме второго закона Ньютона:
m
du
V du
F
dt
2 dt
Преобразуем данное выражение.
24
(2.8)
du
V
(m 
)F
dt
2
 жV
du  D3
 D3
 D3
(
ж 
)
г g 
 ж g  3 ж Du
dt 6
2
6
6
г
18 D ж
du
 D3
 g(

(1

))

u
3
3
dt
 D V
ж
( D  V )  ж
(2.9)
Преобразуем уравнение (2.9), расписав V:
 18

du 6 
  g (1  г )   2 ж
dt 7 
ж  D  ж
 
 u
 
(2.10)
2.2.Оценка для стоксовского режима
Распишем Для случая стационарного всплытия, то есть
получаем
Найдем
du
0
dt
(  ж   г ) gD 2
u0 
18ж
максимальный
диаметр
пузырька
(2.11)
Dmax,
для
которого
выполняется условие стоксовского режима (Re<<1). Подставим уравнение
(2.8) в выражение для числа Рейнольдса получим:
.
 жuD  ж D (  ж   г ) D 2
 ж (  ж   г ) D3
Re 

.
g
g
2
ж
ж
18ж
18ж
(2.12)
Из уравнения (2.12), при Re=1 определим Dmax:
25
Dmax 
3
18ж 2
ж ( ж  г ) g
максимальный диаметр пузырька Dmax, для которого выполняется условие
стоксовского режима, полностью совпадает с аналогичным расчетом для
всплытия пузырька с учетом присоединенной массы. Таким образом, мы
можем воспользоваться расчетами Dmax пузырька в различных растворах,
приведенными в п.1.2.1.
2.3 Определение динамики всплытия пузырька
2.3.1 1определение скорости всплытия пузырька от времени
Решение уравнения (2.8) будет в виде:
u  u0 (1  exp(
Или
18 D ж
 t ))
3
( D  V )  ж
u  u0 (1  exp(
108   ж
 t ))
2
7  D ж
(2.11)
Метод получения данного соотношения аналогичен приведенному
решению уравнения (2.7) в п.2.3.1
2.3.2 Нахождение зависимостей скорости всплытия от времени в
различных растворах
Имеем следующие параметры расчёта:
g  9.80665 м / с 2
. ж  103 кг / м3
 г  1.205кг / м3
26
А) Для раствора воды с вязкостью µ=10-3 Па·с
С помощью приведенных ранее таблицы 1.1 и рис.1.1 найдем u=u(t)
Значит
108 103
u  u0 (1  exp(
 t ))
2
3
7  D 10
Графики зависимостей u(t) с различными диаметрами, представлен на
рисунке 2.1
Рисунок 2.1
Зависимость скорости всплытия u от времени t
для раствора с вязкостью µ= 10-3 Па·с
Б) Для раствора с вязкостью µ= 10-2 Па·с
С помощью таблицы 1.2 и рисунок 1.2. u(t) будет
108 102
u  u0 (1  exp(
 t ))
2
3
7  D 10
График зависимостей u(t) представлен на рисунке 2.2
27
Рисунок 2.2
Зависимость скорости всплытия u от времени t
для раствора с вязкостью µ= 10-2 Па·с
В) Для раствора с вязкостью µ= 10-1 Па·с
С помощью таблицы 1.3 и рисунка1.3 u(t) будет
108 101
u  u0 (1  exp(
 t ))
2
3
7  D 10
График зависимостей u(t) представлен на рисунке 2.3
Рисунок 2.3
Зависимость скорости всплытия u от времени t
для раствора с вязкостью µ= 10-1 Па·с
28
2.3.3 От пройденного расстояния
Решение уравнения (2.8) будет в виде:
x
7B
6  A2
A
Где
A 
A
u

Ln
(1

u)
 B
B 
18   ж
D2  ж
B  (1 
(2.12)
г
) g
ж
Метод получения данного соотношения аналогичен приведенному
решению уравнения (1.8) в п.1.3.2.
2.3.4 Нахождение зависимостей скорости всплытия от расстояния в
различных растворах
А) Для раствора воды с вязкостью µ=10-3 Па·с
Подставляя в (2.12 ) значения u от 0 до значения B/A и значения
диаметров пузырьков получим графики зависимости u(x) (рисунок 2.4)
29
Рисунок 2.4
Зависимость скорости всплытия u от расстояния х
для раствора с вязкостью µ= 10-3 Па·с
Б) Для раствора с вязкостью µ= 10-2 Па·с
С помощью значений диаметров пузырьков из таблицы 1.2получим графики
зависимости u(x) (рисунок 2.5)
Рисунок 2.5
Зависимость скорости всплытия u от расстояния х
для раствора с вязкостью µ= 10-2 Па·с
30
В) Для раствора с вязкостью µ= 10-1 Па·с
С помощью таблицы 1.3.
График зависимостей u(x) представлен на рисунке 2.6
Рисунок 2.6
Зависимость скорости всплытия u от расстояния х
для раствора с вязкостью µ= 10-1 Па·с
Итак мы видим что всплытие пузырька с учетом присоединенной
массы идет медленнее, чем без учета присоединенной массы.
31
3 Экспериментальная установка для исследования движения
пузырьков
3.1 Описание установки
а
б
Рисунок 3.1 – Экспериментальное устройство для исследования
всплытия пузырьков: а – Схема экспериментальной установки: 1 – кювета; 2
– устройство для получения пузырьков;
3 – видеокамера; 4 – масштабная линейка;
б – Фотография экспериментальной установки
Исследование проводилось на установке, состоящей из кюветы с
рабочей жидкостью, устройства для генерации пузырьков и системы
32
визуализации
процесса
всплытия
одиночного
пузырька.
Блок-схема
экспериментальной установки представлена на рисунке 3.1. Кювета 1
призматической
формы
размером
150  150  600
мм
изготовлена
из
оптического стекла толщиной 5 мм. В кювету, установленную строго
вертикально, заливался водно-глицериновый раствор объемом 13.5 л.
Пузырьки воздуха получали с помощью медицинского шприца или
спринцовки 2 объемом 0.15 л, установленных в нижней части кюветы. Для
получения пузырьков разного размера использовались сменные иглы
диаметром (0.63) мм. Система визуализации включала источник света
(электрическую лампу мощностью 100 Вт) и цифровую видеокамеру 3 типа
«Видеоскан-285/п-usb Время экспозиции варьировалось в пределах от 1/350
до 1/1500 с. Для измерения размеров и скорости всплытия пузырька снаружи
кюветы в плоскости движения пузырька установлена масштабная линейка с
ценой деления 1 мм.
С
учетом
геометрические
заданной
размеры
разрешающей
пузырьков
способности
определялись
с
видеокамеры,
относительной
погрешностью 1%.
3.2 Описание видеокамеры
Камера выполнена в цилиндрическом корпусе и имеет интерфейс USB
2.0. Аналого-цифровой преобразователь установлен на камере, что позволяет
получить высококачественное изображение.
Камера имеет следующие технические параметры (Таблица 3.1).
33
Таблица 3.1
Тип ПЗС матрицы
ICX285Al
Формат изображения
1392*10406
Размер пиксела (мкм)
.45*6.45
Размер СCD матрицы (мм)
8.8*6.6
Кадровая частота (Гц)
7.7
Тип затвора
Электронный затвор
Время накопления
3.5 мкс-10 мин
Охлаждение матрицы
Термоэлектрическое, -25 С
Мощность
термоэлектрического 5
холодильника Пельтье (Вт)
Разрядность оцифровки
12
Режимы синхронизации
Внутренняя, программная и внешняя
Посадочное место для крепления С-mount или CS-mount
объектива
Интерфейс камера-компьютер
USB 2.0
Длина кабеля камера-компьютер
5 метров
Потребляемый ток
0.45
Диапазон рабочих температур (град.)
-10-50
Порядок регулировки заднего отрезка объектива:
 Ослабить 3 стопорных винта регулировки заднего отрезка (рисунок 3.2)
 Вращая латунное кольцо с резьбой установить необходимый задний
отрезок.
 Затянуть 3 стопорных винта.
34
Рисунок 3.2
Стопорные винты камеры
Программа Viewer, поставляемая вместе с камерой «ВИДЕОСКАН»,
позволяет
управлять
основными
параметрами
ввода
усилением, экспозицией, выбором режима синхрониации.
35
изображения:
Заключение
В данной работе было исследовано влияние нестационарных и
наследственных сил на динамику всплытия пузырька. Проведена
теоретическая оценка режимов всплытия пузырька в различных растворах, с
учетом приведенной массы.
Для оценки динамики всплытия были получены расчетные данные,
скорости всплытия от времени, представленных на рисунках 1.4-1.9, а так же
от пройденного расстояния (рисунки 2.1-2.6)
Данные результаты использовались при разработке экспериментальной
установки.
36
Список использованной литературы
1.Архипов
В. А. Курс лекций по теории и практике закрученных
потоков./ В. А. Архипов. — Томск : Изд-во ТГУ, 1999. — 60 с.
2.Архипов
В.А.
Образование
вторичных
капель
при
ударном
взаимодействии капли с поверхностью жидкости//прикладная — механика и
техническая физика.—2005.—№1.—С.55-62.
3.Висицкий Е.Г. Движение частиц в вязкой жидкости под действием
силы тяжести и вибрации при наличии силы Бассе // Прикладная математика
и механика-2009.- Т. 73. Вып. 5.- С. 763-775.
4.Водопьянов И.С. О нестационарном осаждении сферической твердой
частицы в вязкой жидкости // Известия РАН. Механика жидкости и газа.2010.- № 2.- С. 98-107.
5.Ландау Л.Д. Теоретическая физика. Т. 6. Гидродинамика/ Ландау
Л.Д., Лифшиц Е.М. - М.: Наука, 1986. 736 с.
6.Майков И.Л. Численная модель динамики капли вязкой жидкости //
Вычислительные методы и программирование.—2009.—№1.—С.148-157.
7.Невский Ю.А. О роли нестационарных и «наследственных» сил в
задачах гравитационной конвекции суспензий // Вест. МГУ. Сер. 1.
Математика, механика.- 2008.- № 4.-С. 37-40.
8.Усанина А.С. Динамика и устойчивость формы капель и пузырьков
при течении вязкой жидкости: дис. … канд. физ.-мат. наук : 01.02.05 / А.С.
Усанина ; Томский государственный университет. – Томск, 2011. – 158 с.
8.Шиляев М.И. Аэродинамика и тепломассообмен газодисперсных
потоков/ Шиляев М.И., Шиляев А.М.- Томск.; Изд-во Том. гос. арх.-строит.
ун-та, 2003. 272 с.
10. O.A. Druzinin Influence of basset force on particle dynamics in twodimensional flows // Physica-D76(1994).- 34-43
37
Download