Analysis of collisions of precipitating solid particles with a wall in a viscous liquid
- Authors: Kondratiev A.S.1, Shvydko P.P.2
-
Affiliations:
- Moscow Polytechnic University
- Sistemotekhnika
- Issue: Vol 17, No 3 (2023)
- Pages: 233-239
- Section: Hydraulic and pneumatic systems
- URL: https://journals.eco-vector.com/2074-0530/article/view/501750
- DOI: https://doi.org/10.17816/2074-0530-501750
- ID: 501750
Cite item
Abstract
BACKGROUND: When calculating velocities of free precipitation of solid particle in a viscous Newtonian liquid, it is assumed that influence of incoming vessel head can be neglected. Necessity of considering this factor grows with a vessel head getting closer and with decrease of particles’ geometrical size.
AIMS: Development of the method of calculating solid particles’ precipitation rate in a Newtonian fluid considering the incoming vessel head and definition the limits of its application.
METHODS: The known analytical model of motion of precipitating disperses solid particles when approaching a vessel head is analyzed. It is shown that the known model of calculation of solid particles’ precipitation rate near the vessel head needs to be refined.
RESULTS: The proposed modification of the analytical model of definition the solid particles’ precipitation rate considering the existence of a vessel head helps to ensure more credible results from physical standpoint. Recommendations regarding the influence of the Brownian motion of fluid’s molecules on motion of precipitating solid particles are given.
CONCLUSION: Practical value of the study lies in ability to calculate solid particles’ precipitation rate near to a vessel head.
Full Text
ВВЕДЕНИЕ
При исследовании осаждения твёрдых частиц в жидкой среде предполагается, что среда состоит из частиц дисперсной твёрдой фазы и вязкой ньютоновской жидкости. В отдельных работах анализируется процесс осаждения мелких частиц твёрдой фазы в жидкой фазе в нижней придонной части двухфазной смеси.
Так, в работе [1] анализировалась возможность столкновения осаждающихся мелких твёрдых тел в вязкой жидкости с нижней горизонтальной стенкой (дном сосуда). Автор работы [1] в примечании указывает, что этот раздел статьи о столкновении осаждающейся твёрдой частицы тела с днищем сосуда выполнен в соавторстве с Забабахиным Н.Е.
В статье [2] практически полностью воспроизведена суть работы со ссылкой на оригинальный теоретический и расчётный материал, ранее изложенный в [1]. В графической части статьи, как и в [1], твёрдая частица представлена как цилиндрическая частица с верхним горизонтальным днищем заданного радиуса и нижним сферическим, радиус которого не указан. Высота вертикальной цилиндрической части частицы не наименована и размер её также не указан.
В статье [2] дополнительно приведены результаты первых экспериментальных исследований, качественно подтверждающих теоретический вывод о невозможности в определённых условиях соударения выпуклых тел с ниже расположенной плоской стенкой сосуда. Однако никаких результатов сравнения поведения сферических частиц с расчётными оценками не приводится.
МЕТОД ПРОВЕДЕНИЯ ИССЛЕДОВАНИЯ
На рис. 1 показана принятая в работах [1, 2] форма частицы в виде цилиндра радиуса с верхним плоским днищем и с выпуклым нижним днищем радиуса , укажем дополнительно высоту цилиндрической части . Частица приближается к плоскому днищу сосуда со скоростью .
Рис. 1. Конфигурация взаиморасположения осаждающейся твёрдой частицы и плоской поверхности днища сосуда.
Fig. 1. Configuration of collocation of a precipitating solid particle and a flat surface of a vessel head.
Ширина зазора между нижней поверхностью частицы и плоской стенкой со временем уменьшается. Таким образом, предполагается, что частица имеет форму цилиндра с плоской верхней и сферической нижней поверхностями. Предполагается также, что режим течения в кольцевом сечении между перемещающейся вниз частицей и плоской стенкой ламинарный. Распределение скорости по толщине зазора принимается параболическим, а при данном средняя по высоте скорость , а градиент скорости у поверхности равен или .
Далее рассматривается элемент жидкого диска и действующие на него силы давления и трения , сумма которых должна быть равна нулю. В результате определяется распределение давления в окружающей жидкости. Далее с учётом давления жидкости на торцы выделенного цилиндрического элемента жидкости определяется суммарная замедляющая сила и записывается уравнение движения осаждающего тела [1, 2]:
. (1)
В уравнении (1) величина обозначает эффективную массу тела, превосходящую его истинную массу на величину присоединённой массы жидкости. Эта масса прибавляется к массе тела, движущегося неравномерно в жидкой среде для учёта воздействия среды на это тело [3]. Величина присоединённой массы зависит от формы тела, направления движения и плотности среды.
Далее для упрощения конечных расчётных выражений полагается, что расстояние от нижнего, придонного контура частицы до неподвижного днища сосуда мало, что позволяет значительно упростить сформулированные расчётные выражения.
Поскольку , и, полагая, что , выражение (1) упрощаем, и оно принимает вид:
. (2)
Интегрируя последнее выражение с граничным условием, что при , , получим
. (3)
Полагая, что частица остановилась , выражение (3) преобразуем в соотношение
. (4)
В работах [1, 2] приводится один и тот же пример расчёта: =1 мг, =0,05 см, =1 мм, =0,01 г/см· с при =– 1 см/с. В обеих статьях расчёт по формуле (4) даёт значение =0,012 см, т. е. частица не касается днища сосуда.
Кроме того, в работе [1] дополнительно приводится второй пример расчёта с отличающимися исходными данными: =1 г, =0,5 см, =1 см, =0,01 г/см·с при = –1 см/с. Расчёт по формуле (4) даёт значение =10–8 см, что много меньше ожидаемой шероховатости поверхности, т. е. частица практически касается днища сосуда [1].
В этих двух примерах, несмотря на значительное различие в массах частиц, в тысячу раз, начальная скорость = –1 см/с принята одинаковой, что традиционно требует специального обоснования, которое отсутствует.
УТОЧНЕНИЕ МЕТОДИКИ РАСЧЁТА ПРИДОННОГО ОСАЖДЕНИЯ ЧАСТИЦ
Вместе с тем использование допущения не является препятствием для проведения интегрирования уравнения (1). Проведя интегрирование дифференциального уравнения (1) в пределах от до , получим выражение:
(5)
Полагая, что частица остановилась при , точнее близком к нулю, выражение (5) преобразуем в соотношение:
(6)
В сравнении с предшествующим вариантом, рассмотренном в работах [1, 2] (см. выражение (4)), возникла зависимость от формы частицы, которая в [2] была принята в виде цилиндра радиуса a c выпуклым дном радиуса , обращенным в сторону днища сосуда.
В табл. 1 представлены результаты расчётов, привёденных в работах [1, 2] и полученных в настоящей работе. Во всех расчётах динамическая вязкость воды равна =0,01 г/см·с
Таблица 1. Сравнение расчётных характеристик осаждения твёрдых частиц вблизи горизонтального днища
Table 1. Comparison of calculated properties of solid particles precipitation near a horizontal vessel head
Обозначения, размерность | Источник | |||
[1] | (6) | (6) | ||
, г | 1 | 1 | 0,001 | 0,001 |
, см | 0,5 | 0,5 | 0,05 | 0,05 |
, см/с | -1 | -1 | -1 | -1 |
, см | 1 | 1 | 0,1 | 0,1 |
, см | 10-8 | 5,4×10-11 | 0,01 | 1,177×10-3 |
, см/с | 5,54×10-5 | 2,777×10-4 |
Сравнение результатов расчётов, выполненных в работах [1, 2], с расчётом с использованием соотношений, полученных в настоящей работе, показывает, что весьма значительно уточняется значение конечного расстояния от вершины тела до неподвижной стенки. Кроме того, появилась возможность определить локальное значение скорости частиц в см/с в этих точках.
Отметим также, что расчётные значения величины находятся на уровне микронов и менее, что соответствует шероховатости поверхности, достигаемой при её шлифовании [4], и представляются физически недостижимыми в реальных условиях заполненного жидкостью сосуда.
С использованием полученных результатов оценим приемлемость допущений, использованных в работах [1, 2], а именно условие
. (7)
В работах [1, 2] не содержится никакой информации о величине радиуса частицы . Поэтому при проведении расчётов в рамках используемой модели принималось, что радиусы цилиндрической части частицы и радиус её выпуклой донной части равны, т. е. . В этом случае условие (7) преобразуется в соотношение
. (8)
Расчётные данные, проведённые в табл. 1, показывают, что условие (8) выполняется во всех случаях при расположении частиц вблизи нижней неподвижной поверхности. Условие (8) не выполняется в исходном состоянии с использованием соотношений приведённых в [1, 2].
Ранее выполненные исследования [5] показали, что в случае твёрдых частиц произвольной формы целесообразно использовать три характерных диаметра частицы: — куб диаметра шара, эквивалентного объёму частицы; — квадрат диаметра шара, эквивалентного по площади боковой поверхности частицы; — диаметр шара, эквивалентного по площади миделевому сечению частицы. В выражении для числа Рейнольдса используется эквивалентный диаметр частицы . С использованием таких представлений рассчитывались силы гидравлического сопротивления для твёрдых частиц, выполненных в виде двух сферических частиц, диске, иглоподобных эллипсоидах. Во всех случаях для двух различных ориентаций частиц относительно вертикального направления осаждения частиц получено удовлетворительное соответствие с результатами точных численных расчётов [5]. Этот подход был распространён и на осаждение полимодальных двухфазных смесей. В этом случае учитывалось также столкновение между частицами разных фракций твёрдой фазы. В дальнейшем такой подход использовался для расчёта движения полимодальных смесей твёрдых частиц в горизонтальных, вертикальных и наклонных трубопроводах [6]. При этом дополнительно учитывался процесс диффузии твёрдых частиц в вертикальном направлении.
Во всех случаях при исследовании осаждающихся твёрдых частиц и движении их в потоках жидкости необходима полная информация о фактической геометрической конфигурации твёрдых частиц. В работах [1, 2] фигурирует лишь один непосредственно геометрический параметр — радиус нижнего выпуклого днища твёрдой частицы . Величина определяет массу частицы, но без указания плотности вещества частицы, что не позволяет определить её объём и далее определить, скорее даже оценить, величины и . Поэтому, к сожалению, не представляется возможным корректно рассчитать начальное значение скорости .
Скорость осаждения твёрдой частицы произвольной формы при использовании величины коэффициента гидравлического сопротивления по формуле Стокса для сферических частиц можно представить в виде зависимости [5]:
. (9)
Положим так же, как было принято в [1, 2], что при осаждении крупных части , а в случае мелких частиц — .
Используя выражение (9), определим отношение скоростей осаждения крупных частиц к скоростям осаждения мелких частиц.
, (10)
где ; ; ; ; ; ; ; .
Подставляя приведённые численные значения в выражение (10), получим, что отношение , т. е. скорости осаждения частиц значительно различаются, а не равны, как предполагается в работах [1, 2].
Применительно к сферическим осаждающимся частицам такой же вывод следует и из известной формулы Стокса для осаждающихся сферических частиц: . В этом случае при различии диаметра частиц в 10 раз скорости осаждения также различаются в 100 раз.
При определении скорости свободного осаждения твёрдых частиц, как правило, не учитывается влиянием Броуновского движения на скорость свободного осаждения частиц. Как отмечается в [7], Броуновское движение оказывает важное влияние на двухфазные потоки, включающие мельчайшие твёрдые частицы. Если исходить из условия, что кинетическая энергия осаждающихся частиц много больше энергии теплового движения молекул, то условие, что Броуновским движением можно пренебречь, записывается в виде [5].
, (11)
где — постоянная Больцмана, — температура в градусах Кельвина.
Характерное значение диаметра частиц, при котором необходимо учитывать влияние Броуновского движения на скорость осаждения твёрдых частиц, определяется выражением [5].
. (12)
Расчёты величины правой части в формуле (10) при различных значениях плотности материала твёрдых частичек в воде при представлены в табл. 2.
Таблица 2. Значения критического диаметра частиц твёрдой фазы
Table 2. Values of critical diameter of solid particles
, кг/м3 | 1250 | 1500 | 2000 | 2500 | 3000 | 3500 | 4000 | 4500 |
, мм | 0,116 | 0,095 | 0,078 | 0,069 | 0,064 | 0,060 | 0,057 | 0,054 |
Из приведённых расчётных данных следует, что пренебречь влиянием Броуновского движения молекул несущей жидкости на движение осаждающихся твёрдых частиц можно лишь для частиц высокой плотности.
ЗАКЛЮЧЕНИЕ
Проведённый анализ показал, что использование в работах [1, 2] допущения, что , является физически некорректным при анализе осаждения твёрдых частиц вблизи плоского днища. Показано, что возможно получение аналитического решения без привлечения указанного допущения, что позволяет более точно определить координаты частицы вблизи плоского днища.
ДОПОЛНИТЕЛЬНО
Вклад авторов. А.С. Кондратьев ― постановка задачи исследования, анализ публикаций по теме статьи, разработка обобщённой расчётной модели, анализ результатов расчётов, редактирование текста рукописи; П.П. Швыдько — библиографический поиск, проведение численных исследований, подготовка текста рукописи, создание изображений. Авторы подтверждают соответствие своего авторства международным критериям ICMJE (все авторы внесли существенный вклад в разработку концепции, проведение исследования и подготовку статьи, прочли и одобрили финальную версию перед публикацией).
Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи.
Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования.
ADDITIONAL INFORMATION
Authors’ contribution. A.S. Kondratiev ― statement of the research problem, search for publications on the topic of the article, writing the text of the manuscript, development of a generalized analytical model, analysis of calculation results, editing the text of the manuscript; P.P. Shvydko ― bibliographic search, numerical research, preparation of the text of the manuscript, creating images. All authors made a substantial contribution to the conception of the work, acquisition, analysis, interpretation of data for the work, drafting and revising the work, final approval of the version to be published and agree to be accountable for all aspects of the work.
Competing interests. The authors declare that they have no competing interests.
Funding source. This study was not supported by any external sources of funding.
About the authors
Alexander S. Kondratiev
Moscow Polytechnic University
Author for correspondence.
Email: ask41@mail.ru
ORCID iD: 0009-0001-2030-5165
SPIN-code: 6739-0496
Dr. Sci. (Tech.), Professor of the Industrial Heat Power Engineering Department
Russian Federation, MoscowPavel P. Shvydko
Sistemotekhnika
Email: Shvydko.P@gmail.com
ORCID iD: 0009-0007-8030-6403
SPIN-code: 8347-3886
Project Manager
Russian Federation, MoscowReferences
- Zababakhin EI. Some cases of viscous fluid motion. In: Questions of modern experimental and theoretical physics. Leningrad: Nauka; 1984:58–68. (in Russ).
- Zdeshchits VM, Sidorenko VD. Collision of solids in a viscous fluid. Visnik KTU. 2010;26:59–61. (in Russ).
- Sedov LI. Plane problems of hydrodynamics and aerodynamics. Moscow: Nauka; 1980. (in Russ).
- GOST 9378-93. Obraztsy sherokhovatosti poverkhnosti (sravneniya). Mezhgosudarstvennyy standart. Minsk: Mezhgosudarstvennyy sovet po standartizatsii, metrologii i sertifi-katsii; 2002. (in Russ). [дата обращения: 17.05.2023] Режим доступа: https://meganorm.ru/Data2/1/4294847/4294847416.pdf
- Kondratyev AS. Settlement of polymodal solids in Newtonian fluids. Moscow: Sputnik+; 2014. (in Russ).
- Kondratyev AS, Nya TL, Shvydko PP. Methodology and methods for calculating the movement of dispersed solid particles in pipelines. Moscow: Mos. Politekh; 2020. (in Russ).
- Batchelor J. The influence of Brownian motion on the average stress in a suspension of spherical particles. In: Hydrodynamic interaction of particles in suspensions. Mechanics. New in foreign science. Moscow: Mir;22:124–153. (in Russ).