Гидродинамика жидкой пленки на проницаемой цилиндрической поверхности Текст научной статьи по специальности «Физика»

Трулы БГТУ, 2017, Серия 3, № 2, с. 15-19

15

УДК 532.517

А. М. Волк

Белорусский государственный технологический университет

ГИДРОДИНАМИКА ЖИДКОЙ ПЛЕНКИ НА ПРОНИЦАЕМОЙ ЦИЛИНДРИЧЕСКОЙ ПОВЕРХНОСТИ

Выполнен анализ исследований гидродинамики пленочных движений в сепарационных, фильтровальных, тепло- и массообменных, газожидкостных аппаратах и реакторах. Проведено теоретическое исследование пленочного движения жидкой фазы на внутренней поверхности проницаемого цилиндра под воздействием массовых сил поля тяжести и закрученного газового потока. Впервые получены дифференциальные уравнения движения, найдены точные решения для составляющих скорости при условии прилипания пленки на поверхности проницаемого цилиндра и равенстве касательных напряжений на границе раздела фаз, определена толщина пленки. Найдены гидродинамические характеристики пленочного течения и аналитические условия для возможных режимов движения при условии отвода жидкой фазы через проницаемую поверхность.

Ключевые слова: аппараты, пленочное движение, закрученный поток, проницаемая поверхность, точные решения, гидродинамические характеристики, режимы движения.

A. M. Volk

Belarusian State Technological University

HYDRODYNAMICS OF LIQUID FILM ON A PERMEABLE CYLINDRICAL SURFACE

The studies of the hydrodynamics of film movements in the separation, filtration, heat and mass transfer, gas-liquid reactors and apparatuses are analyzed. The theoretical study of film movement in the liquid phase on the inner and outer surfaces of the cylinder under the influence of the gravitational field of mass forces and swirling gas flow is done. For the first time, the differential equations of motion is obtained, the exact solutions for the velocity components subject to adherence of the film on the surface of a permeable cylinder and equality of shear stresses at the interface are found, the film thickness is determined. The hydrodynamic characteristics of film flow and analytical conditions for the possible modes of motion are found.

Key words: apparatus, film motion, swirling flow, permeable surface, exact solutions, hydrody-namic characteristics, motion regimes.

Введение. Гидродинамика пленочного течения на проницаемых поверхностях имеет существенное значение для процессов фильтрования суспензий, отвода жидкой фазы в процессе сепарации газожидкостных потоков, при тепломассообмене. В этих случаях пленочное течение осуществляется на проницаемых поверхностях.

Отсос используется также для управления пограничным слоем и влияния на устойчивость ламинарного режима движения [1, 2].

Воздействия газового потока на пленку жидкости передаются посредством сил трения, возникающих на границе взаимодействия фаз.

Касательные напряжения на стенке трубы диаметром В = 2Я определяются гидродинамическим напором [1]:

= PrWz 4 2

(1)

При степенном распределении осевой скорости газового потока с показателем 1/7 коэффициент трения для гладких труб вычисляется по формуле Блазиуса [1]:

X0 =

0,3164

Re

0,25

(2)

который зависит от числа Рейнольдса Re = ZD.

При наличии пленочного течения на внутренней поверхности трубы волновая поверхность рассматривается как нерегулярная шероховатость стенок канала. Обобщением экспериментальных данных в этом случае для коэффициента трения на границе раздела фаз получено соотношение [3, 4]

X = X01 1 +

3005 D

(3)

Для определения тензора касательных напряжений в исследуемых диапазонах изменения скорости газового потока могут быть использованы рассмотренные ранее зависимости, которые дают осевую составляющую

V

т

R

16

Гидродинамика жидкой пленки на пронинаемой цилиндрической поверхности

Т z _

0,0791 ргЖ?

Re0

2

1 +

3005 D

(4)

Касательная составляющая зависит от угла закрутки потока и при этом будет

T=Tztg (¿в).

(5)

Установившееся пленочное течение характеризуется равновесием сил тяжести и возникающих сил трения.

Теоретически исследуются только автомодельные течения с постоянным отсосом или вдувом по всей поверхности.

Скорость оттока сплошной среды определяется свойствами проницаемой поверхности и перепадом давления на ней [5, 6]. При турбулентном режиме движения среды через отверстия в стенке на основании уравнения Бернул-ли получаем квадратичный закон:

U0 =

2AP__ 2AP

P(1 + S" f2)

(6)

Для потоков, не имеющих относительной продольной скорости, коэффициент расхода у зависит только от свойств проницаемой поверхности. Идельчиком [7] при обработке экспериментальных данных получена формула зависимости коэффициента расхода у от относительной площади / отверстий:

Y-2 _

(0,707^1-7 +1 + f)2

f2

(7)

Для исследования процессов разделения многофазных потоков необходимо математическое моделирование движения закрученного газового потока и пленочного течения жидкости с учетом взаимодействия фаз как внутри непроницаемого цилиндра, так и при оттоке жидкой фазы на проницаемых поверхностях.

Экспериментальные исследования более полно характеризуют происходящие процессы, но они не всегда позволяют выделить и проанализировать доминирующие факторы.

Гидродинамика данного процесса в значительной степени зависит от оттока жидкой фазы. В настоящее время практически не имеется адекватной модели и метода расчета процесса разделения фаз при пленочном течении суспензии.

Наиболее полными и достоверными данные и расчеты будут при сочетании экспериментальных и теоретических методов анализа исследуемых процессов.

Математическая модель. Рассмотрим установившееся осесимметричное течение вязкой

несжимаемой жидкости по внутренней стенке проницаемого цилиндра под воздействием закрученного газового потока (рис. 1). Ось г цилиндрической системы координат направим вниз по оси цилиндров.

Рис. 1. Схема двухфазного пленочного течения на проницаемой цилиндрической поверхности

Скорость оттока жидкой фазы П0 на некотором элементарном цилиндре длиной Аг будем считать постоянной. Объемный расход несжимаемой жидкости через цилиндрические поверхности равной длины будет одинаков: 2шПгДг = 2п^иоДг . Отсюда находим ради-

п и оК

альную скорость в пленке жидкости и г =-.

Тогда из уравнения неразрывности

V _ 0

1 Э t л \ dUz

1 rUr )+"1Т

r dr dz

(8)

получим

U

dz

_ 0 и Uz _ Uz (r).

dP

Принимаем у _ — _ const. С учетом oce-d z

симметричности ^^ = 0 уравнения Навье -Эф

Стокса для составляющих скорости и давления приводим к виду:

r

й2иг 1 ( и0Я _ ^ йи„

йг2

г V V

йг

Ц

'Ч 1Г ЦЯЛ^ 1( М+ = 0,

йг

2

йг г

V

ф

Цф + Ц2Я2

г V V

дР дг

В результате получаем систему обыкновенных дифференциальных уравнений. Это означает, что решение и = и (г) будет автомодельным. Выполним переход к безразмерной коорди-

и0 Я

нате г = г / Я , обозначим а =- и получим:

йи, (а-1) йи2 = Р£-уЯ2

йг

Ц

(9)

й2иф (а-1) йиф (а +1)

йг2

дР дг

=Р

йг

(и2+и02 Я2 ^

ч

г

„3

г

иф= 0, (10)

(11)

/

Для граничных условий используем значения составляющих скорости жидкости на цилиндрической поверхности и компонентов тензора касательных напряжений:

ди,

г ц _ Цф ^

дг г

(12)

на границе раздела фаз.

Частные решения уравнений (9-11) ищем в

виде гк и получаем общие решения:

и, = с + с2г ^ г2; (13)

2ц(а _ 1)

иф="Г + с4г г

а+1

(14)

За граничные условия принимаем условие прилипания на стенке и равенство касательных напряжений (12) на границе раздела фаз:

Ш-= =иф| г=1 =0; т, =_

Ц йи,

Я йг

г =1_8

тф=_Ц

Г д_и±_ Цф ^ Ядг Яг

(15)

(16)

'=1_8

Из условия равновесия сил, действующих на газовый поток,

п(я _ 8^АР = 2п(я _ 8)1;/

получим

АР 2< у = — =

2т,

/ я (1 _8) я (1 _8)'

Учитывая граничные условия (15, 16), получаем распределение скорости в пленке жидкости:

и =

тЯ

+

Р&Я2

ц(а_2)(1 _8) " ца(а_2)(1 _8) х(1 _ г а)_

х

РёЯ2

- + -

Т, Я

2ца(а_ 2) ц(а_ 2)(1 _8) ЯТф(1 _8 )2 (1

(1 _ г 2), (17)

и,

ча+2

1 _8) + 2ц

- _ г~а+Ч. (18)

Найдем объемный расход жидкой фазы на единицу периметра:

д = Я | и гМг =

1_8

pgЯ3

ц(а_2)(1 _8) " ца(а_2)(1 _8) 1 _(1 _8)2 1 _(1 _8)а+

Р&Я

- + -

а+2

т Я2

2ц(а_ 2) ц(а_ 2)(1 _8) 1 _(1 _8)2 1 _(1 _8)

2

4

(19)

Изменение данного объемного расхода по длине описывается уравнением

, = _и о.

Разложение правой части в ряд до четвертой степени включительно имеет вид

Т;Я2

Ц 2

+

Р&Я

Ц

283 3а - 8

-+-

3!

4!

4!

84+...

- +...

+

Р

г

г

2

+

х

2

х

2

4

3

д

18

Гидродинамика жидкой пленки на проницаемой цилиндрической поверхности

Тг б2 1

Ц 2

1 иоб

1

+ -8

( П02б2 б2 ^ —0--1--

ч V2 Я2 ,

ч

+...

+

+

р^б3

ц

I+-1 Г зПоб

3 241 v

8 Я'+-

(20)

Данное разложение показывает, что изменение гидродинамических характеристик вследствие оттока жидкой фазы происходит, когда

и 0 б

безразмерный комплекс Реинольдса - со-

V

измерим с единицей.

Среднее значение тангенциальной составляющей скорости пленки и перепад давления в радиальном направлении будет:

- 1 1

Пф = б | П^г = ф

(1 -б )2

' 1-б

х

бЦ а(1 -б) (1 -б )

а+2

-х

+ 2

а+2 а+2

1п (1 -б)

ДР =Р |

1 < и2

1-8

ф+ ио Я ^

г

г

3

ёг .

(21)

(22)

Исследуем пленочное течение воды в проницаемом цилиндре диаметром 0,15 м, под воздействием газового потока со средней расходной скоростью 15 м/с с углом закрутки в 45°. Расход жидкости примем равным 0,5 м3/ч, а скорость оттока жидкой фазы - 0,001 м/с. В данном случае из уравнения (20) получена толщина пленки 1,1 мм, а на рис. 2, 3 представлены графики составляющих скорости в зависимости от переменной у = Я - г .

Графические зависимости показывают, что отток жидкой фазы через проницаемую поверхность существенно изменяет профили скоростей жидкой пленки. Распределение скоростей в пленке жидкости на проницаемой поверхности является линейным.

Вывод. Создание моделей процессов разделения, происходящих в конструируемых аппаратах, является важной и достаточно сложной задачей. Математическое описание позволяет выполнить расчет основных характеристик,

исследовать влияние на процесс различных факторов, сочетать теоретические и экспериментальные методы исследования, определять оптимальные конструкционные и технологические параметры.

П г, м/с

0,20

0,15

0,10 0,05

0 0,0002 0,0004 0,0006 0,0008 У, м

Рис. 2. Осевая составляющая скорости в пленке жидкости

0 0,0002 0,0004 0,0006 0,0008 У, м

Рис. 3. Касательная составляющая скорости в пленке жидкости

Полученная математическая модель дает возможность определить гидродинамические характеристики пленочного течения под воздействием закрученного газового потока с учетом оттока жидкой фазы и учитывать их при исследовании процессов фильтрования, а также найти оптимальные соотношения между геометрическими и режимными параметрами исследуемых процессов.

Литература

1. Шлихтинг Г. Теория пограничного слоя. М.: Наука, 1969. 742 с.

2. Ерощенко В. М., Зайчик Л. И. Гидродинамика и тепломассообмен на проницаемых поверхностях. М.: Наука, 1984. 274 с.

3. Уоллис Г. Б. Одномерные двухфазные течения. М.: Мир, 1972. 440 с.

4. Соколов В. И., Доманский И. В. Газожидкостные реакторы. Л.: Машиностроение, 1976. 216 с.

5. Берд Р., Стьюарт В., Лайтфут Е. Явления переноса. М.: Химия, 1974. 688 с.

1

6. Гельперин Н. И. Основные процессы и аппараты химической технологии: в 2 кн. М.: Химия, 1981. 812 с.

7. Идельчик И. Е. Справочник по гидравлическим сопротивлениям. М.: Машиностроение, 1975. 560 с.

References

1. Schlikhting G. Teoriya pogranichnogo sloya [The theory of the boundary layer]. Moscow, Nauka Publ., 1969. 742 p.

2. Yeroshchenko V. M., Zaychik L. I. Gidrodinamika i teplomassoobmen na pronitsayemykh poverkhnostyakh [Hydrodynamics and heat and mass transfer on permeable surfaces]. Moscow, Nauka Publ., 1984. 274 p.

3. Uolis G. B. Odnomernye dvukhfaznye techeniya [One-dimensional two-phase flows]. Moscow, Mir Publ., 1972. 440 p.

4. Sokolov V. I., Domanskiy I. V. Gazozhidkostnye reaktory [Gas-liquid reactors]. Leningrad, Mashinostroenie Publ., 1976. 216 p.

5. Berd R., St'yuart V., Laytfut Ye. Yavleniya perenosa [Transport phenomena]. Moscow, Khimiya Publ., 1974. 688 p.

6. Gel'perin N. I. Osnovnye protsessy i apparaty khimicheskoy tekhnologii [The main processes and apparatus of chemical technology]: in 2 books. Moscow, Khimiya Publ., 1981. 812 p.

7. Idel'chik I. Ye. Spravochnik po gidravlicheskim soprotivleniyam [Reference book on hydraulic resistance]. Moscow, Mashinostroenie Publ., 1975. 560 p.

Информация об авторе

Волк Анатолий Матвеевич — кандидат технических наук, доцент кафедры высшей математики. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: [email protected]

Information about the author

Volk Anatoliy Matveevich - PhD (Engineering), Assistant Professor, the Department of Higher Mathematics. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]

Поступила 05.05.2017