CC BY
19
МЕХАНИКА
Н. Л. Великанов, В. А. Наумов, С. И. Корягин
ДИНАМИКА ДИСПЕРСНОЙ ЧАСТИЦЫ, СДУВАЕМОЙ С ПЛОСКОЙ ГОРИЗОНТАЛЬНОЙ ПОВЕРХНОСТИ ЛАМИНАРНЫМ СЛОЕМ КУЭТТА
Определен характер влияния силы Саффмэна на выход частицы из слоя. Продольная скорость частицы может оказаться больше скорости жидкости, тогда сила Саффмэна будет направлена к поверхности. В этих условиях нужна большая начальная вертикальная скорость, чтобы частица достигла внешней границы слоя.
Character of influence Saffmena force on the particles out of the layer is detected. The longitudinal velocity of the particle can be more fluid speed, then the power Saffmena will be directed to the surface. In these circumstances we need a big initial vertical velocity of the particle to reach the outer layer of the border.
Ключевые слова: дисперсные частицы, динамика, линейный профиль скорости.
Key words: dispersed particles, dynamics, linear velocity profile.
В расчетах дисперсных течений в водотоках, в системах гидро- и пневмотранспорта при математическом моделировании образования ветропесчаных потоков необходимо исследовать движение частиц вблизи твердой поверхности [1 — 3]. Однако при этом нередко без всякого обоснования не учитываются подъемные силы.
При движении дисперсной частицы в сдвиговом потоке возникает поперечная сила, обусловленная разностью давлений в направлении градиента скорости (чем меньше скорость, тем больше давление) [4; 5]. В классической работе Саффмэна [6] для случая линейного профиля скорости несущей среды и сферической частицы получена формула
где р/ — плотность жидкости; б — диаметр частицы; Ц/, Щ — продольная скорость жидкости и частицы соответственно; У — поперечная координата; V — коэффициент кинематической вязкости жидкости. Причем коэффициент С а = 6,46. Анализ Саффмэна был ограничен условиями малости чисел Рейнольдса [7]
90
g2 ( FSaf = CSaf Р f (Uf - Up ) V
( л U Y/2
8
Re = |VK| 5 / v < 1, Rem=|QR| 52 v < 1,
© Великанов Н. Л., Наумов В. А., Корягин С. И., 2016
Вестник Балтийского федерального университета им. И. Канта. Сер.: Физико-математические и технические науки. 2016. № 4. С. 90 —
Ие А=-
52 д и
V д У
«1, а=^
_ Рк
Ие
1/2
д и
д У
-«1,
где У„ , Ю„
модуль относительной линейной и угловой скорости
частицы соответственно.
Рассмотрим поступательное движение сферической частицы, сдуваемой с плоской горизонтальной поверхности ламинарным слоем Куэтта. Ось Ох направим по поверхности вдоль течения, а ось Ог — вертикально вверх. Кроме сил тяжести, Архимеда, гидродинамического сопротивления, эффекта присоединенных масс, необходимо учесть силу Саффмэна. Для малых чисел Рейнольдса проекции уравнения движения на продольную и поперечную (вертикальную) оси имеют следующий вид [4; 8]:
1+Х] ^ = Р и - и), Р=^
2] й^ 1 р 52
Х = ^,
1+2
йТ
= -р ^ - g + а(иг - и )
Гд Щ У/2 3,08Х^
д У
а
(1)
(2)
91
Учтем, что в слое Куэтта продольная скорость жидкость растет линейно с расстоянием от стенки и^ = 2 /Ь. Перейдем в (1) — (2) к безразмерным переменным, выбрав в качестве характерных параметров и„ — скорость жидкости на внешней границе слоя и Ь — его толщину:
«■=Ы1+1
а = 3,08
йи й1
йш -ш й1 Я
БИ(м.
г - и
-| + а(г - и),
1+х
18X1 2
Я = ^ Л '
Ие5 =
и 5
4Я
1+ Х/2
1 =
Я
1 -X
18ХРг 1+ Х/2
Рг = ^
gL
и=ш=^ х = Х ¿ =2 t = и*Т
и и , ш и„, х Ь, г Ь, Ь ■
(3)
(4)
л=Ь,
5
Начальные условия:
х(0) = ¿(0) = 0, хх(0) = и0, ¿(0) = ш0 > 0.
(5)
Задача Коши (3), (4), (5) имеет аналитическое решении, впервые полученное в [8] для определения условий осаждения частицы в пристеночном слое. Применим метод [8] для решения задачи динамики частицы, сдуваемой с плоской горизонтальной поверхности.
а
V
р
р
V
92
Выразим из уравнения (3) вертикальную координату частицы и ее производные:
2 = и • StkM + и, т = 2 = и • StkM + и, т = и • StkM + и.
Подставим эти выражения в уравнение (4):
и • StkM + и = -(и • StkM + й) / StkM -Е, + а(и • StkM + и - и).
После преобразования отсюда получаем линейное неоднородное дифференциальное уравнение третьего порядка с постоянными коэффициентами:
и +-
Stkм
-и +
V
-а
и = -
Stkм
Начальные условия с учетом (3):
и(0) = и0, и(0) = -u0/StkM, и(0) = ( т0 + u0/StkM) / StkM.
(6)
(7)
Стандартная процедура решения задачи Коши (6) — (7) после подстановки в (4) дает при т ^ 0, т0 > 0 и а^ 1/ Stk2M :
= с1ехр (X11) + с2ехр (X 2t)-В, В =
l•StkN
1 - а • StkM
^) = £- (ехр (X ^)- 1)+"Х2- (ехр (X 21)-1)-Bt,
X1 X 2
с =_ 1 2
1
с, =-
т0 + и0
л/а-
1/ StkM + 7а
(8) (9) (10)
тг
- иа4а-
1/StkM -Та
X х =--1--у[а, X 2 = —
ш 1
StkM StkM
(11)
Исследуем, при каких условиях частица достигнет внешней границы слоя. Заметим, что йт(0)/^ = -т>0/Stkм - аио < 0, т. е. в начале движения т(£) — убывающая функция. Рассмотрим сначала случай, когда функция т(£) имеет минимум. Из условия йт/ ^ = 0 следует
йехр^^) + C2exp(X2t) = 0, так что функция т(£) достигает минимума в момент времени
^ =
2>/а
1п
V X2C2 J
(12)
Согласно (10), (11) < 0, а > 0, тогда Ь > 0 при
-Х1С1 > Х2С2 > 0. Последнее условие эквивалентно системе неравенств
Stk,
-^¡a
w0 + u0Va+-
Stk,
-^¡a
-1 StkM
-sfä
w,
- u0Va-
StkM
--Va
StkM
y[a I w0 - un\fa -
1/ StkM --\/a
> 0.
(13)
Из (13) последнего неравенства следует ад > -(амо + БЬкш, что выполняется автоматически при м0 ^ 0. Условие (13) может быть удовлетворено двумя способами:
a < -
StkM
a > -
StkM
w0 < uQy[a-
w0 > uQ-Ja-
>/a -1/ StkM '
4 .
л/a -1/ StkM
(14)
(15)
Наиболее интересный для практики случай, когда либо м0 = 0, либо м0 — малая положительная величина. Тогда из (14) т0 < 0, но это значит, что частица не покинет поверхность. Поэтому условия (14) далее не рассматриваются. Подстановка условий (15) в (8) с учетом (12) показывает, что в этом случае и>(Ь) > 0. Следовательно, при выполнении (15) частица обязательно покинет слой через верхнюю границу.
Однако найденные условия (15) являются достаточными для ухода частицы из слоя, но не необходимыми. Действительно, частица может покинуть слой через верхнюю границу при монотонно убывающей ш(Ь), но достаточно большой начальной скорости т0. Примером служат линии 3 на рисунках 5 — 6 и на рисунках 9 — 10. Найдем момент времени ^ в который т(Ь) = 0 (после чего частица начинает двигаться вниз):
C1exp(Xtt2) + C2exp(^2t2) = B.
(16)
Из (16) численным методом определяется Далее необходимо найти условия, при которых > 1, для чего должно быть численно решено уравнение
[exp(Xtt2) - 1]ct/X1 + [exp(^2t2) - 1]c2/- Bt2 = 1.
(17)
Поиск условий выхода частицы из слоя через верхнюю границу путем решения (16) — (17) является весьма трудоемкой задачей. Целесообразно ее решать построением зависимостей (8), (9) в среде Mathcad при различных значениях критериев подобия, что показано на рисунках 1—10.
93
\
1
1
94
"IV 0.08 0.06 0.04 0.02
О
• • * • J у /
V • • у / *
* * * V * У
и Л*-.-, чч
1 2
10
15 20
25
30
35 40
Рис. 1. Изменение безразмерной вертикальной скорости частиц по времени при Я = 0,033, X = 0,0012, Рт = 400, и0 = 0 и различных значениях Ш0: 1 - Ш0 = 0,03, 2 - Ш0 = 0,04, 3 - Ш0 = 0,05, 4 - Ш0 = 0,08
Рис. 2. Изменение безразмерной вертикальной координаты частиц по времени при Я = 0,033, X = 0,0012, Рт = 400, и0 = 0 и различных значениях Ш0: 1 - Ш0 = 0,03, 2 - Ш0 = 0,04, 3 - Ш0 = 0,05, 4 - Ш0 = 0,08
0.06 0.04 0.02
0
* * * * *
* * * У р*
К 1
1 2 ......
• _
10
15
20
25
Рис. 3. Изменение безразмерной вертикальной скорости частиц по времени при Я = 0,033, X = 0,0012, Рт = 400, ш0 = 0,08 и различных значениях и0: 1 - и0 = 0, 2 - и0 = 0,04, 3 - и0 = 0,06, 4 - щ = 0,08
Рис. 4. Изменение безразмерной вертикальной координаты частиц по времени при Я = 0,033, X = 0,0012, Рг = 400, Wo = 0,08 и различных значениях М0: 1 - М0 = 0, 2 - М0 = 0,04, 3 - М0 = 0,06, 4 - М0 = 0,08
95
0.08 0.06 0.04 0.02
0
■ * * •
4„- • • •
W" *
V" "-Л 2
\ 1 ч %
10
15
20
25
30
35
Рис. 5. Изменение безразмерной вертикальной скорости частиц по времени при R = 0,05, X = 0,002, w0 = 0,08, u0 = 0 и различных числах Фруда: 1 - Fr = 100, 2 - Fr = 150, 3 - Fr = 155, 4 - Fr = 250
9 • y y
• • * • — ta V *
• * * * y > 4 t V
Vy 1 V 1
/ --- X l i l i t
0 5 10 15 20 25 30 35 40 t
Рис. 6. Изменение безразмерной вертикальной координаты частиц по времени при R = 0,05, X = 0,002, w0 = 0,08, u0 = 0 и различных числах Фруда: 1 - Fr = 100, 2 - Fr = 150, 3 - Fr = 155, 4 - Fr = 250
, У г
г
п »2
*.Ч3 • * •
• • ■С:-
1 0 *
Рис. 7. Изменение безразмерной вертикальной скорости частиц по времени при К = 0,05, Рт = 600, ш0 = 0,08, и0 = 0 и различных отношениях плотностей: 1 - X = 0,001, 2 - X = 0,002, 3 - X = 0,003, 4 - X = 0,004
0 2 4 б 8 10 12 14 16 15 *
Рис. 8. Изменение безразмерной вертикальной координаты частиц по времени при К = 0,05, Рт = 600, ш0 = 0,08, и0 = 0 и различных отношениях плотностей: 1 - X = 0,001, 2 - X = 0,002, 3 - X = 0,003, 4 - X = 0,004
• • 4.*
* * • *
1 • 1 * с * Ц *
3
XI * * I- ___
0 10 20 30 40 50 60 t
Рис. 9. Изменение безразмерной вертикальной скорости частиц по времени при X = 0,0012, Рт = 600, ш0 = 0,08, и0 = 0 и различных значениях К: 1 - К = 0,029, 2 - К = 0,03, 3 - К = 0,0304, 4 - К = 0,032
0.Î
0.6
0.4
0.2
m m * V f * t'
V * * //
• • • * 4 'if' T + i
• *V> 1
0
10 15 20 25 30 35 40 45
97
Рис. 10. Изменение безразмерной вертикальной координаты частиц по времени при X = 0,0012, Рт = 600, ш0 = 0,08, и0 = 0 и различных значениях К: 1 - К = 0,029, 2 - К = 0,03, 3 - К = 0,0304, 4 - К = 0,032
Таким образом, можно считать установленным характер влияния силы Саффмэна на выход частицы из слоя. Параметр а характеризует отношение силы Саффмэна к силе гидродинамического сопротивления. При Mo = 0 чем больше а, тем меньше должна быть начальная скорость частицы для ее выхода из слоя. При м0 > 0 зависимость «выходной» скорости от а более сложная: функция имеет минимум, после которой указанная зависимость становится возрастающей. Действительно, продольная скорость частицы M может оказаться больше скорости жидкости, тогда сила Саффмэна будет направлена к поверхности. В этих условиях нужна большая начальная вертикальная скорость, чтобы частица достигла внешней границы слоя.
Полученные результаты могут служить основой для совершенствования расчетных схем гидросистем [1; 2; 9].
Список литературы
1. Великанов Н. Л., Наумов В. А., Космодамианский А. С. Моделирование осаждения твёрдых частиц в пульпопроводе // Наука и техника транспорта. 2011. № 2. С. 69-78.
2. Великанов Н.Л., Наумов В. А., Примак Л. В. Осаждение частиц взвесей в воде // Механизация строительства. 2013. № 7. С. 44—48.
3. Великанов Н. Л., Наумов В. А. Взвешенные частицы в водотоке и их осаждение // Вода: химия, экология. 2014. № 2. С. 114 — 119.
4. Naumov V.A. Influence of Saffman's lift force on the motion of a particle in a Couette layer // J. of Engineering Physics and Thermophysics. 1995. Vol. 68, № 5. P. 840-844.
5. Наумов В. А. Динамика дисперсной частицы в вязкой среде // Математическое моделирование. 2006. Т. 18, № 5. С. 27 — 36.
6. Saffman P. G. The lift on a small sphere in a slow shear flow // Journal Fluid Mechanics. 1965. Vol. 22. P. 385-400.
7. Асмолов Е. С. О динамике сферической частицы в ламинарном пограничном слое // Известия АН СССР. МЖГ. 1990. № 6. C. 91—96.
98
8. Брич М. А. Выпадение капель из дисперсного потока на нагретую поверхность // Тепло- и массоперенос: итоги и перспективы. Минск, 1985. С. 50—55.
9. Великанов Н. Л., Корягин С. И. Проблемы эксплуатации и ремонта металлических и железобетонных конструкций. Калининград, 2015.
Об авторах
Николай Леонидович Великанов — д-р техн. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград.
Е-mail: monolit8@yandex.ru
Владимир Аркадьевич Наумов — д-р техн. наук, проф., Калининградский государственный технический университет.
Е-mail: van-old@rambler.ru
Сергей Иванович Корягин — д-р техн. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград. Е-mail: SKoryagin@kantiana.ru
About authors
Prof. Nikolay Velikanov — I. Kant Baltic Federal University, Kaliningrad. E-mail: monolit8@yandex.ru
Prof. Vladimir Naumov — Kaliningrad State Technical University. E-mail: van-old@rambler.ru
Prof. Sergey Koryagin — I. Kant Baltic Federal University, Kaliningrad. E-mail: SKoryagin@kantiana.ru
