CC BY
32
2016 Математика и механика № 3(41)
УДК 532.529:662.62 Б01 10.17223/19988621/41/8
О.В. Матвиенко, А.В. Андриасян, Н.А. Мамадраимова, А.О. Андропова
ИССЛЕДОВАНИЕ ДВИЖЕНИЯ ЧАСТИЦЫ В ФОРМЕ ВЫТЯНУТОГО ЭЛЛИПСОИДА ВРАЩЕНИЯ В ЗАКРУЧЕННОМ ПОТОКЕ
Проведено исследование движения частицы в форме вытянутого эллипсоида вращения в закрученном потоке. Установлено, что в зависимости от ориентации частицы в пространстве возможно как ее восходящее, так и нисходящее движение. С использованием модели дрейфа определены условия восходящего и нисходящего движения. Предложены зависимости, позволяющие определить компоненты вектора скорости частицы.
Ключевые слова: механика жидкости, частицы, дисперсная фаза, закрученный поток, сепарация.
В пылеочистительной технике большое распространение получили циклоны различных конструкций, однако принцип их работы одинаков и основан на использовании центробежной силы [1 - 3].
При центробежном разделении выбросу придается вращательное движение внутри циклонного аппарата, при этом твердые частицы отбрасываются центробежной силой на периферию аппарата к его стенке, так как центробежное ускорение в циклоне на несколько порядков больше ускорения силы тяжести, что позволяет удалить из выброса даже весьма мелкие частицы [4].
Для определения эффективности сепарации необходимо учитывать, помимо условий работы, свойства разделяемой смеси, дисперсность твердых частиц, вязкость дисперсионной среды, разницу плотностей разделяемых фаз, концентрацию вещества в жидкой фазе.
Отметим, что подавляющее большинство публикаций относится к изучению процессов разделения частиц сферической формы [5 - 9]. Однако в реальных ситуациях форма частицы может существенно отличаться от сферической. В результате распределения скорости жидкости и давления вблизи сферической и несферической частицы будет различным. И, как следствие этого, силы гидродинамического сопротивления, действующие на несферическую частицу, будут отличны от сил, действующих на частицу сферической формы. Это, в свою очередь, приведет к изменению траекторий движения частиц в аппарате и повлияет на се-парационные характеристики.
Целью работы является исследование движения в закрученном потоке неизометрических твердых частиц, имеющих форму вытянутого эллипсоида вращения.
Поле течения жидкости
Рассмотрим движение жидкости в зазоре между двумя коаксиальными цилиндрами разного размера. Оба цилиндра могут как вращаться с постоянными угловыми скоростями, так и покоиться. Поскольку рассматриваемое течение жидкости
можно считать плоским и осесимметричным, то уравнения гидродинамики в цилиндрических координатах примут вид [10]
ёг
г з ± | %
ёг 1 г
= 0, % = РД. (1)
ёг г
Здесь г - радиальная координата; у[(р - тангенциальная скорость жидкости; р1 -плотность жидкости; р - давление.
В качестве граничных условий используются условия прилипания на стенках:
г = гг : % = Ю , г = ге : % = Юеге . (2)
В уравнениях (2) гг, ге - радиусы внутреннего и внешнего цилиндров соответственно юг, юе - угловые скорости их вращения.
Проинтегрировав уравнение (1) с граничными условиями (2), получим радиальное распределение тангенциальной скорости жидкости:
у = аеге2 — Юг? г геЮе — Юг ())
% = 2 2 г 1 2 г • ()
г — г г — г г
е г е г
Если внешний и внутренний цилиндры вращаются с одной угловой скоростью юе = юг = ю , то распределение тангенциальной скорости соответствует закону вращения твердого тела Укр =юг .
Уравнения движения частицы
При движении частицы в потоке жидкости или газа на нее действует сила сопротивления со стороны несущей среды. Если вектор скорости частицы направлен в сторону большей полуоси эллипсоида или перпендикулярно ей, то сила сопротивления будет направлена в сторону, противоположную вектору скорости [)]. Если частица движется в направлении, перпендикулярном большей полуоси эллипсоида, то и в этом случае вектор скорости и вектор силы сопротивления будут направлены в противоположные стороны. При этом величина силы сопротивления в первом и втором случае будет различной. Если же вектор скорости центра масс частицы направлен под углом к главным осям эллипсоида, то векторы скорости и силы сопротивления перестают быть коллинеарными. Величина силы сопротивления будет определяться скалярным произведением тензора сопротивления К и вектора скорости частицы относительно несущего потока:
¥0 =— впцаК •(V, — Vр ), (4)
где ц - вязкость несущей среды, а - характерный размер частицы, V, - скорость несущего потока, V - скорость частицы.
В декартовой систем координат с использованием индексной формы записи и соглашения о суммировании Эйнштейна уравнение (4) может быть записано в виде
РВг = —впЦаКц (Уу — ур} ) . (5)
В случае сферической симметрии все ) диагональные компоненты тензора сопротивления будут равны между собой и закон сопротивления естественным образом переходит в хорошо известный закон Стокса.
В том случае, если частица является ортотропной, т.е. имеет 3 взаимно перпендикулярные плоскости симметрии, то тензор сопротивления в главных осях характеризуется главными значениями Кь К2 и К3. В случае осевой симметрии две из трех диагональных компонент тензора будут равны между собой, что имеет место для эллипсоида вращения [3]:
К =
Re
(ß2 - 1
3/2
18
К 2 = Кз =
(2ß2 -l)ln(ß + V(1 )-ß#^
(ß2 -1)
C
Re
(2ß2 - 3)ln(ß^(ß2 -1) ) + ß^/(ß2 -1)
Cr
(6)
(7)
где CD - коэффициент сопротивления твердой сферической частицы, ß = b/a характеризует отношение полярного и экваториального диаметра частицы. Формулы (6), (7) справедливы для вытянутых эллипсоидов вращения ß> 1.
Коэффициент сопротивления одиночной твердой частицы CD в простейшем случае является однозначной функцией относительного числа Рейнольдса Re = р\vl -vp\dI^ . В качестве характерного размера, необходимого для построения числа Рейнольдса, можно использовать эффективный диаметр частицы, равный диаметру сферической частицы того же самого объема. Для вытянутого эллипсоида вращения эффективный диаметр частицы равен d = 2a # . При низких числах Рейнольдса Re < 1 коэффициент сопротивления сферической частицы определяется формулой Стокса: CD = 24/Re. В переходной области кривая сопротивления описывается различными формулами. В частности, стандартную кривую сопротивления можно аппроксимировать степенными зависимостями Бабухи -Шрайбера [11]: CD = A Re-" . Значения параметров в диапазоне изменения числа Рейнольдса 1 < Re < 10 равны: A = 26.3, " = 0.8 . Для чисел Рейнольдса, лежащих I в диапазоне 10 < Re < 1000 рекомендуется ис-
пользовать следующие значения [12]: A = 12.3, " = 0.5 . При Re > 103 картина обтекания стабилизируется, что в первом приближении приводит к независимости CD от Re : CD = 0.44 . Последнее соотношение известно как закон сопротивления Ньютона.
Движение частицы в закрученном потоке удобно рассматривать в цилиндрической системе координат с базисными векторами er, еф, ez. Построим единичный вектор n, направление которого совпадает с направлением большей полуоси частицы. Тогда ориентацию частицы в пространстве можно задать через Рис. 1. Ориентация частицы в потоке углы прецессии у и нутации 6 (рис. 1). Fig. 1. Particle orientation in the flow
Предполагается, что в процессе движения частицы углы у и 6 остаются неизменными. Последнее возможно при стабилизации заряженной частицы электромагнитным полем [13].
При таком выборе координат компоненты тензора сопротивления определяются следующим образом:
Krr = ( - K2)cos2 у • sin2 6 + K2, Kpp = (K1 - K2) sin2 у • sin2 6 + K2, Kzz = (Ki - K2 )cos2 6 + K2, Knp = 1(Ki - K2 )sin(2у) • sin2 6 ,
Krz = 1 (Ki - K2 )cosy sin(26), Kz(p= 2(Ki - K2 )siny sin(26).
С использованием теоремы об изменении количества движения скорость центра масс частицы может быть найдена из решения системы дифференциальных уравнений:
% = 2 П [ Krz (ir - Vpr ) + Kp (p - Vpp ) + Kzz ( - Vpz )]-(p F~Pl) g $ (8) at 2 Ppa pp
^ = 2П-^^[Krp(vir -Vpr) + Kpp(vip-vj + Kzp(viz -Vpz)]-]$ (9) at 2 ppa r
^ = 2 П^ [ Krr ( - Vpr ) + Krp(vip-vj + Krz ( - Vpz ^.(10)
dt 2 p pa Pp r
Для определения траектории центра масс частицы используются кинематические соотношения [14]
drP dzP d Р P (11)
V =-— V =-— V = r -— (11)
pr dt' pz dt ' ppp p dt ' K '
Система уравнений (7 - 10) замыкается начальными условиями
t = 0 : Vpz = Vpz0 , Vpr = Vpr0 , Vpp = Vpp0 , zp = zp0 , rp = rp0 , Рp = Рp0.
Анализ результатов
Рассмотрим движение эллипсоидальных вытянутых частиц в закрученном потоке. Движение мелких частиц относительно несущей жидкости достаточно мало. В результате этого самые мелкие частицы движутся по винтовой траектории практически по цилиндрической поверхности. Крупные частицы, оттесняемые центробежной силой, движутся к стенкам внешнего цилиндра по конической поверхности. С увеличением скорости вращения внешнего цилиндра число витков, увеличивается, шаг винтовой линии уменьшается. С увеличением угловой скорости вращения цилиндров происходит рост значений центробежной силы, которая интенсифицирует радиальное движение частицы к внешнему цилиндру. В результате происходит увеличение угла конусности.
На рис. 2 показаны траектории частицы для различных значений угла нутации. Из рисунка видно, что в зависимости от ориентации частицы в пространстве возможно как ее восходящее, так и нисходящее движение.
2 N
0,3
0,110 0,105 0,100 0,095 0,090 0,085
-0,3
Рис. 2. Траектории движения вытянутых эллипсоидальных частиц: d = 10-4 м, в = 10 , ю, = юе = 10 рад/с; у = 0, 1 - 0 = 60° , 2 - 75°, 3 - 90°, 4 - 105°;
Fig. 2. Motion trajectories of prolate ellipsoidal particles:
d = 10-4m, p = 10 , ю, = юе = 10rad/s; y = 0, 1 - 0 = 60° , 2 - 75°, 3 - 90°, 4 - 105°
Такое движение объясняется особенностями действия на частицу силы сопротивления. Сила сопротивления, действующая на эллипсоидальную частицу, со стороны несущей жидкости характеризуется наличием горизонтальной и вертикальной составляющих. Горизонтальная составляющая силы сопротивления для тяжелых частиц, плотность которых превосходит плотность несущей среды, направлена к оси симметрии. Направление вертикальной составляющей силы сопротивления зависит от ориентации частицы в пространстве.
Влияние угла нутации 0 на движение частицы представлено на рис. 3, где показаны меридиональные проекции траектории, определенные для различных значений угла нутации. Для значений угла 0 <0<15° и 15°<0< 85° движение частиц происходит по нисходящей траектории, для 15°<0<85° - по восходящей. Максимальная высота подъема достигается при 0 = 43°, наибольшая глубина -при 0 = 142° .
Рис. 3. Меридиональные проекции траектории движения вытянутых эллипсоидальных частиц: d=10-3 м, у=0, Р=10, =<ве =10 рад/с; а: 1 - 6 = 0, 2 - 15°, 3 - 30°, 4 - 45°; б: 1 - 6 = 60° , 2 - 75°, 3 - 90°, 4 - 105°; в: 1 - 6 = 120°, 2 - 135°, 3 - 150°, 4 - 165° Fig. 3. Meridional projections of the motion trajectory of prolate ellipsoidal particles: d = 10-3 m, у = 0 , P = 10 , m¡ = = 10 rad/s; а: 1 - 6 = 0, 2 - 15°, 3 - 30°, 4 - 45°; б: 1 - 6 = 60° , 2 - 75°, 3 - 90°, 4 - 105°; в: 1 - 6=120°, 2 - 135°, 3 - 150°, 4 - 165°
Определим условия, при которых движение частицы является восходящим, а при каких - нисходящим.
Из анализа уравнений движения (8 - 10) следует, что на начальном участке происходит достаточно резкое изменение скорости частицы от начального значения до некоторого равновесного значения, определяемого балансом сил, действующих на частицу и не зависящего от начальных условий. Время выхода на равновесный режим движения может быть найдено как
г - Р^2 (12)
64 пцК2
Для мелких частиц время достижения динамического равновесия намного меньше характерного гидродинамического времени /¿уп — , что позволяет
для анализа движения использовать модель дрейфа частиц [9]. При этом скорость частиц определяется в предположении малости инерционных членов. Таким образом, компоненты вектора скорости частицы могут быть найдены из решения системы алгебраических уравнений:
2 п^г [ Кгг ( - УрГ) + Кгф(ф - Урф) + К22 ( - Ург)]—(р р-Р1) g ; (13) 2 Рра Рр
9п^[КгА( - V)+Крф(( -^Ф) +Кгф( -Vpг)] — 0 ; (14)
2 Рра
V2,
2nVaF^Кгг ( -Vpr) + Кгф ( -Vpv) + Krz ( -vр2)] = -(( р^ - С15)
Система уравнений (13) - (15) является системой нелинейных алгебраических уравнений. Преобразование этой системы приводит к кубическому уравнению, решение которого можно осуществить, например, методом Кардано. Однако полученное решение является достаточно громоздким. Поэтому остановимся на анализе нескольких частных случаев.
В первом случае будем полагать, что угол прецессии у = 0 . Для квазитвердого
вращения жидкости: vlz = 0, vlr = 0, vl<f = rar , система уравнений (13) - (15) примет вид
Krzvpr + Kzzvpz = - 2¿eqK (P,-Pl ) g ; (16)
9 P p
vpT=0; (17)
Krrvpr + KrZvpz = 2¿eqK ( Pp -P l ) ra2r . (18)
Pp
Решение системы (16) - (17) позволяет определить компоненты скорости частицы:
vpr = v* (Fr • Kzz - Krz ) , vpz = v* (Krr - Fr • Krz ) vpT = rar , (19)
где v* =— t -—— g - стационарная скорость осаждения частицы, большая
2 (Рр -Р1) 9 64 Р К
ось которой является вертикальной, Бг = ю2г^ - число Фруда.
На рис. 4 показаны зависимости радиальной и осевой составляющих скорости от угла нутации 6 , рассчитанные как с использованием дифференциальной модели (8) - (10), так и модели дрейфа (19). Сравнение полученных результатов позволяет сделать вывод о возможности применения модели дрейфа для анализа движения частиц. Как видно из рисунков, максимальная скорость радиального
Рис. 4. Зависимость радиальной (а) и осевой (б) составляющих скорости от угла нутации 0 : d = 10-3 м, у = 0, Р = 10, << = << = 10 рад/с; 1, 3 - расчет с использованием модели (10) - (13), кр. 2, 4 - с использованием модели дрейфа, кр. 1, 2 - Fr = 5 , кр. 3, 4 - Fr = 2 Fig. 4. Dependence of the (а) radial and (б) axial velocity components on nutation angle 0 : d = 10-3m, y = 0, P = 10, << = << = 10rad/s; 1, 3 -calculations according to models (10) - (13), 2, 4 - drift model is used, 1, 2 - Fr = 5, 3, 4 - Fr = 2
движения частицы к стенке достигается при значении угла 6 = 90°. В этом случае большая ось частицы ориентирована вдоль радиуса камеры, поэтому при движении в радиальном направлении частица испытывает наименьшее сопротивление. Напротив, при 6 = 0 сопротивление радиальному движению со стороны несущего потока становится максимальным. Это приводит к уменьшению значений vpr. С увеличением скорости вращения несущего потока эффект центробежной сепарации становится более значимым, поэтому с ростом значений числа Фруда Fr возрастает радиальная составляющая скорости движения частицы vpr.
Скорость гравитационного осаждения уменьшается при значениях угла нутации 0 <6< 45°, 135° <6< 180° и возрастает при 45° <6< 135°. При этом максимальная скорость осаждения частицы достигается при 6 = 135°, минимальная -6 = 45°. С увеличением значений числа Фруда Fr скорость гравитационного осаждения частицы уменьшается. При достаточно сильной закрутке потока движение частицы станет восходящим. С помощью уравнения (19) определим критическое значение числа Фруда Fr*, разделяющее режим восходящего и нисходящего движений: Fr* = Krr/Krz.
Если угол прецессии составляет у = 90°, то компоненты скорости частицы будут
KK
^ = -K-2 ' . 2 6 v* , vw = % , vpr = KLFr ■v* . (20)
Kj cos 6 + K2 sin 6 K2
При такой ориентации частицы возможно только нисходящее движение. Скорость гравитационного осаждения зависит от значений угла 6 . При 6 = 0 достигается максимальная скорость осаждения, при 6 = 90° - минимальная.
Рассмотрим теперь случай, когда угол нутации 6 = 0. В этом случае равновесные значения составляющих скорости частицы будут равны:
2 K2 (рp-Pi)
V = - 9 'eq J--P-g = - 2v* , vPT = % ,
2 . (Pp -Pi) 2 Fr ^ = 9 'eq-p^ »r = Kv*. (21)
Как следует из (21)? движение частицы для заданной ориентации осей всегда будет нисходящим. При угле нутации 6 = 0 скорость движения частицы не зависит от значения угла ф .
При значении угла 6 = 90° компоненты вектора скорости частицы будут
K 1
vpz =- ^ф =-—--
f w 2Fr v*
2
1 - 11 + 4Fr—
V %;
KФ
фф
v = -тг
Pr ъ^
агф
(ф- Vpф). (15)
Скорость гравитационного осаждения не зависит от величины угла прецессии и определяется только размером и формой частицы. Скорость движения частицы в тангенциальном направлении оказывается отличной от скорости движения жидкости и зависит от размеров, формы частицы и угловой скорости вращения жидкости. При этом, чем больше число Фруда, тем больше скорость движения частицы в тангенциальном направлении отличается от скорости несущего потока. Ско-
рость радиального движения частицы определяется не только размером и формой частицы, а также ее ориентацией в пространстве.
Таким образом, проведенные исследования показали, что движение частиц существенно определяется их положением в потоке, и для создания расчетных методик сепарации и пылеулавливания необходимо учитывать этот факт.
ЛИТЕРАТУРА
1. Crowe C., Sommerfeld M., Tsuji Ya. Multiphase Flows with Droplets and Particles. CRC Press, 1998. 472 p.
2. Островский Г.М. Прикладная механика неоднородных сред. СПб.: Наука, 2000. 359 с.
3. Кутепов A.M., Полянин Л.Д.. Запрянов З.Д. и др. Химическая гидродинамнка: справочное пособие. М.: Бюро Квантум, 1996.
4. Матвиенко О.В., Евтюшкин Е.В. Математическое исследование сепарации дисперсной фазы в гидроциклоне при очистке вязкопластических буровых растворов // Инженерно-физический журнал. 2011. Т. 84. № 2. С. 230-238.
5. Clift R., Grace J.R., Weber M.E. Bubbles, drops and particles. NY: Academ Press, 1978. 380 p.
6. Rubinow S.I., Keller J.B. The transverse force on spinning sphere moving in a viscous fluid //J. of Fluid Mech. 1965. V. 22. P. 385-400.
7. Матвиенко О.В., Данейко А.М. Исследование ударного взаимодействия частиц в потоке // Изв. вузов. Физика. 2013. Т. 56. № 9-3. С. 190-192.
8. Матвиенко О.В., Андропова А.О. Исследование движения частицы в потоке жидкости вблизи подвижной стенки // Вестник Томского государственного университета. Математика и механика. 2015. № 4. С. 85-92.
9. Матвиенко О.В., Евтюшкин Е.В., Андропова А.О. Применимости модели дрейфа частиц для моделирования переноса дисперсной фазы в потоке // Вестник Томского государственного университета. Математика и механика. 2015. № 5 (37). С. 76-83.
10. Шлихтинг Г. Теория пограничного слоя. М.: Наука, 1974.
11. Матвиенко О.В., Базуев В.П., Агафонцева М.В. Исследование динамики пузырька в закрученном потоке нелинейно-вязкой жидкости // Вестник Томского государственного архитектурно-строительного университета. 2012. № 4 (37). С. 144-156.
12. Хаппель Дж., Бреннер Г. Гидродинамика при малых числах Рейнольдса. М.: Мир, 1976.
13. Bottner C.U., Sommerfeld M. Numerical calculation of powder painting using the Euler/Lagrange approach // Powder Technology. 2002. V. 125. P. 206-216.
14. Матвиенко О.В., Агафонцева М.В. Численное исследование процесса дегазации в гидроциклонах // Вестник Томского государственного университета. Математика и механика. 2012. № 4(20). С. 107-118.
Статья поступила 11.02.2016 г.
Matvienko O.V., Andropova A.O., Andriasyan A.V., Mamadraimova N.A. (2016) INVESTIGATION OF THE PROLATE ELLIPSOIDAL PARTICLE MOTION IN A SWIRLING FLOW. Tomsk State University Journal of Mathematics and Mechanics. 2(41). pp. 74-85.
DOI 10.17223/19988621/41/8
In this paper, the investigations of the prolate ellipsoidal particles movement in a swirling flow were carried out.
The motion of small particles relative to the carrier fluid flow is insignificant. As a result, the smallest particles move along a helical path at an approximately cylindrical surface. Large particles affected by centrifugal force move to the walls of the outer cylinder along a conical surface. Increasing of the angular velocity of the cylinders leads to an increase in the centrifugal force which intensifies the radial motion of the particles to the outer cylinder. It is found that, depending on the particle orientation in the space, both upward and downward motions are possible.
The movement is caused by particular drag force influence on the ellipsoidal particle. The drag force is characterized by the horizontal and vertical components. The horizontal component for heavy particles with density value greater than that of the carrier medium is directed towards the axis of symmetry. The direction of the vertical component of the drag force depends on the orientation of the particles in the space.
Utilizing the particle slip velocity model, the regimes of upward and downward movement are defined. Dependencies for the determination of particle velocity components are proposed. It is found that the gravitational settling velocity of the particle decreases with an increase in the value of the Froude number. In the case of a strong swirling flow, the motion of the particle becomes upward.
Keywords: mechanics of fluid, particles, dispersed phase, swirling flow, separation.
MATVIENKO Oleg Viktorovich (Doctor of Physics and Mathematics, Tomsk State University,
Tomsk, Russia)
E-mail: matvolegv@mail.ru
ANDROPOVA Antonina Olegovna (Tomsk State University) E-mail: a.o.andropova@gmail.com
ANDRIASYANAlina Valerrievna (Tomsk State University of Architecture and Building) E-mail: alina-andriasyan@yandex.ru
MAMADRAIMOVA Nazira Abdiraimovna (Tomsk State University of Architecture and Building) E-mail: alina-andriasyan@yandex.ru
REFERENCES
1. Crowe C., Sommerfeld M., Tsuji Ya. (1998) Multiphase Flows with Droplets and Particles. CRC Press.
2. Ostrovsky G.M. (2000) Prikladnaya mekhanika neodnorodnykh sred [Applied mechanics of heterogeneous media]. Saint Petersburg: Nauka.
3. Kutepov A.M., Polyanin L.D.. Zapryanov Z.D. (1996) Khimicheskaya gidrodinamika [Chemical Hydrodynamics]. Moskow: Byuro Kvantum.
4. Matvienko O.V., Evtyushkin E.V. (2011) Mathematical study of hydrocyclone dispersed phase separation in clearing viscoplastic drilling fluids. Journal of Engineering Physics and Thermophysics. 84(2). pp 241-250. DOI 10.1007/s10891-011-0466-z.
5. Clift R., Grace J.R., Weber M.E. (1978) Bubbles, drops and particles. NY: Academ Press.
6. Rubinow S.I., Keller J.B. (1965) The transverse force on spinning sphere moving in a viscous fluid. J. of FluidMech. 22. pp. 385-400.
7. Matvienko O.V. Daneiko A.M. (2013) Issledovanie udarnogo vzaimodeystviya chastits v potoke [Investigation of shock interaction of particles in a stream]. Izvestiya vysshikh uchebnykh zavedeniy. Fizika. - Russian Physical Journal. 56(9/3). pp. 190-192.
8. Matvienko O.V., Andropova A.O. (2015) Issledovanie dvizheniya chastitsy v potoke zhid-kosti vblizi podvizhnoy stenki [Studying the particle motion in a fluid flow in the vicinity of a movable wall]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika -Tomsk State University Journal of Mathematics and Mechanics. 4(36). pp. 85-92. DOI 10.17223/19988621/36/10.
9. Matvienko O.V., Evtyushkin E.V., Andropova A.O. (2015) Issledovanie primenimosti modeli dreyfa chastits dlya modelirovaniya perenosa dispersnoy fazy v potoke [Studying the applicability of the algebraic slip model for predicting the dispersed phase motion in the flow]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 5 (37). pp. 76-83. DOI 10.17223/19988621/37/7.
10. Shlikhting G. (1964) Grenzschicht-Theorie. Springer-Verlag. (In German).
11. Matvienko O.V., Bazuev V.P., Agafontseva M.V. (2012) Issledovanie dinamiki puzyr'ka v zakruchennom potoke nelineyno-vyazkoy zhidkosti [Investigation of the dynamics of a bub-
ble in a swirling flow of nonlinear-viscous liquid]. Vestnik Tomskogo gosudarstvennogo arkhitekturno-stroitel'nogo universiteta - Bulletin of the Tomsk State Architectural University. 4(37). pp. 144-156.
12. Happel J., Brenner H. (1965) Low Reynolds number hydrodynamics: with special applications to particulate media. Prentice-Hall.
13. Bottner C.U., Sommerfeld M. (2002) Numerical calculation of powder painting using the Euler/Lagrange approach. Powder Technology. 125. pp. 206-216.
14. Matvienko O.V., Agafontseva M.V. (2012) Chislennoe issledovanie protsessa degazatsii v gi-drotsiklonakh [Numerical study of the degassing process in hydrocyclones]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 4(20). pp. 107-118.
