УДК 532.517
А. М. Волк, кандидат технических наук, доцент (БГТУ)
ДВИЖЕНИЕ ВЯЗКОЙ ЖИДКОСТИ В ПРОСТРАНСТВЕ МЕЖДУ ПРОНИЦАЕМЫМИ ПЛОСКОСТЯМИ
Исследование механики движения вязкой жидкости в пространстве, ограниченном проницаемыми поверхностями, имеет научное и практическое значение. Структура потока в пространстве между поверхностями и в пограничном слое играет важную роль при решении задач тепло- и массообмена и разделении гетерогенных систем. В данной работе выполнено математическое моделирование и исследованы гидродинамические характеристики движения сплошных потоков между движущимися горизонтальными проницаемыми плоскостями. Получено точное решение для уравнения Навье - Стокса, описывающего рассматриваемое движение при постоянной скорости в нормальном направлении, а также зависимости для оценки пограничного слоя и гидродинамические характеристики движения в основном потоке. Найдены интегральные характеристики течения, показано, что направление и интенсивность движения жидкости через проницаемую поверхность имеют существенное значение.
The research of mechanics of viscous fluid motion in the space limited to permeable surfaces has scientific and practical value. The stream structure in the space between surfaces and in the interface is of high importance for the problems solving of heat-and-mass transfer and divisions of heterogeneous systems. The mathematical modeling has been carried out and hydrodynamic characteristics of movement of solid streams between moving horizontal permeable planes have been investigated in this work. The exact solution for the Navier - Stoks equation has been received. This equation is to describe motion under study at constant speed in a normal direction. The relations for estimation of interface and hydrodynamic characteristics of movement in the main flow have been received. The integral characteristics of stream have been found, it has been shown that the direction and intensity movement of fluid through permeable surfaces is of an high value.
Введение. Исследование механики движения вязкой жидкости в пространстве, ограниченном проницаемыми поверхностями, имеет научное и практическое значение. Структуру потока в пространстве между поверхностями и в пограничном слое важно учитывать при решении задач тепло- и массообмена и разделении гетерогенных систем.
При математическом моделировании процессов переноса необходимо принимать во внимание гидродинамические факторы, взаимодействие потоков с ограничивающими поверхностями, отток среды через проницаемую поверхность.
Теоретически исследуются только автомодельные течения с постоянным отсосом или вдувом по всей поверхности.
В работе [1] исследование течения на плоской проницаемой поверхности и в трубе рассматривается в соответствии с теорией пограничного слоя, изложенной Шлихтингом в его монографии [2].
Аналитические решения имеют предпочтение по сравнению с экспериментальными результатами. Они позволяют получить интегральные характеристики течения [3]: - толщину вытеснения
толщину потери импульса
8 2 = J
hW f
0 Wh
1 --
W„
\
W
dy;
h у
толщину потери энергии
hWx f
8 = J^-
83 J w
0 Wh
W 1--^
Y
W
dy.
h У
Точное аналитическое решение для продольного обтекания плоской пластины с равномерным оттоком исследовано в работах [1, 2]. Движение потока с осевой скоростью V в зависимости от поперечной координаты у описано уравнением
U,
dU
■ = v-
d 2U
dy dy
(1)
Данное уравнение имеет аналитическое решение при известных скоростях оттока ио и на бесконечности V/ :
f
U = U„
Uo y Л
1 - e
U0 < 0.
(2)
81=Ц1 - W dy;
0 V h
\
Для толщины вытеснения и толщины потери импульса соответственно получены значения:
8,= У
- U0
48
ISSN 1683-0377. Труды БГТУ. 2011. № 6. Физико-математические науки и информатика
82 =
1
2-U
Касательное напряжение т0 на стенке не зависит от вязкости:
(
т0 = ц
dU
= p(-Uo)U .
j=0
Отмечается, что отток значительно увеличивает силу трения с поверхностью, оказывает стабилизирующее действие на ламинарный пограничный слой, вызывает перераспределение профиля скорости, уменьшает степень турбулентности потока. При сильном оттоке силы трения на поверхности не зависят от вязкости, профили скорости становятся близкими к плоским, может наблюдаться полный отток массы, вблизи оси канала появляется зона обратных токов [1, 2].
Математическая модель. Рассмотрим стационарное движение вязкой жидкости в пространстве между проницаемыми плоскостями. Декартовую систему координат х, у свяжем со второй плоскостью, через которую происходит отток жидкости. Пусть первая плоскость находится на расстоянии к от второй и движется с относительной скоростью Жк = Ж1 - Ж2. Ось х направим параллельно вектору Жк, а ось у - по нормали к поверхности. Считаем, что движение жидкости в нормальном направлении происходит с постоянной скоростью Ж0 (рис. 1).
У h
wh
ZZ3ZZZZZZ3Z
Wo
^0
Рис. 1. Схема движения жидкости между проницаемыми плоскостями
Из граничных условий
, 4=0 =0;
=к = ^к определяем произвольные постоянные и получаем
Wx = -
W„
1 - exp
Wo h
1 - exp
Wo у"
(5)
Безразмерный комплекс Яеу = |Ж0| к / V является числом Рейнольдса и связывает зависимость продольного течения жидкости с ее движением в нормальном к плоскостям направлении.
Интегрируя распределение скорости, получаем среднее значение продольной составляющей скорости жидкости в пространстве между проницаемыми плоскостями:
Wx = Wh
1
1
1 - exp(- Re у) Re>
(6)
Находим производную продольной скоро -сти и касательные напряжения сил трения:
dWx dy
WoWh
v[1 - exp(- Re у)]
exP
Wo у"
(7)
т = n-
dWx
pWoWh
dy
1 - exp(- Re у)
На поверхности с оттоком
т = —
PWoWh
1 - exp(- Re у)
expl —°у |. (8)
(9)
Если Re у > 5, то можем принять
Wx = Wh
1 - exp
Wo у
и Wx = Wh (1 -1„у
i1 - Re -).
При больших значениях Яе у продольное течение на поверхности с оттоком характеризуется большим градиентом скорости:
dWx
dy
PWoWh
у =o
В данной системе координат Ж0 < 0. Уравнение Навье - Стокса для рассматриваемого движения принимает вид [3, 4]
Wo
dWx
d 2WX
■ = v-
dy йу2 Находим его общее решение:
Жх = С1 + С2 ехр I Ж(} у
(3)
(4)
а касательные напряжения сил трения незначительно зависят от вязкости.
В этом случае т| у=0 = -рЖ0Жк.
Принимая
dWx
dy
= 0,01, получим оценку
7=8
для толщины пограничного слоя:
WoWh
exp
Wo81 = o,o1; 8 = —— In Wh
1ooWoWh
v
V
V
V
x
V
V
V
Построим график скорости (рис. 2) воды в пространстве между движущимися с относительной скоростью 0,5 м/с плоскостями при расстоянии между ними 0,5 м. Нормальную скорость примем равной 0,01 м/с.
0,05.
0,04-
0,03
0,02'
0,01
0,00
0,0 0,1 0,2 0,3 0,4 0,5 Скорость жидкости, м/с
Рис. 2. Скорость жидкости в пространстве между проницаемыми плоскостями
Анализируя полученные результаты можем сделать вывод, что распределение продольной составляющей скорости жидкости в пограничном слое близко к линейному:
Жх =-
«у
а_в основном потоке является равномерным: Жх = Жк.
Выполнив интегрирование, находим: 1) толщину вытеснения
81 =-
V
к
Ж0 Яе у
2) толщину потери импульса
8 2 =-
V
к
2Ж0 2Яе у
3) толщину потери энергии
8 - 5 V =
5к
6 Ж0 6Яе у
Применяя теорему энергии для плоского ламинарного пограничного слоя, получаем изменения потока энергии на единицу длины в пограничном слое:
±(83)- ^№ ,
ёх
ду
2Ж02Жк
Такое же значение получается и при использовании теоремы энергии для турбулентного течения в пограничном слое:
I ^ )= 21Р
(дЖх РI ду
\
ёу =
Это означает, что возникающие при оттоке жидкости силы трения вызывают большие затраты энергии на нагревание. Данные затраты характеризуют потери, вызванные разделением потоков при относительном движении проницаемых плоскостей и движением среды в нормальном направлении.
Заключение. Полученное точное решение уравнения движения вязкой жидкости в пространстве между проницаемыми плоскостями и найденные гидродинамические характеристики течения показывают, что направление и интенсивность движения жидкости через проницаемую поверхность имеют существенное значение.
Характер течения у поверхности с притоком и в основном потоке определяется вязкостью жидкости, зависит от линейной скорости этой поверхности и интенсивности движения в нормальном направлении. Течение у поверхности с оттоком характеризуется большим градиентом скорости, наличием пограничного слоя. Силы трения на этой поверхности, вообще говоря, не зависят от вязкости жидкости.
Полученная математическая модель с учетом реологических факторов может применяться для анализа процессов разделения многофазных потоков в пространстве, ограниченном проницаемыми поверхностями.
Литература
1. Козлов, Л. Ф. Ламинарный пограничный слой при наличии отсасывания / Л. Ф. Козлов. -Киев: Наукова думка, 1968. - 195 с.
2. Шлихтинг, Г. Теория пограничного слоя / Г. Шлихтинг. - М.: Наука, 1969. - 742 с.
3. Лойтянский, Л. Г. Механика жидкости и газа: учеб. для вузов / Л. Г. Лойтянский. - 6-е изд., перераб. и доп. - М.: Наука, 1987. - 840 с.
4. Ландау, Л. Д. Механика сплошных сред / Л. Д. Ландау, Е. М. Лифшиц. - М.: Гостехиздат, 1953. - 123 с.
Поступила 28.02.2011
2
V
V
V