Анализ столкновения осаждающихся твёрдых частиц со стенкой в вязкой жидкости

Полный текст

ВВЕДЕНИЕ

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

Так, в работе [1] анализировалась возможность столкновения осаждающихся мелких твёрдых тел в вязкой жидкости с нижней горизонтальной стенкой (дном сосуда). Автор работы [1] в примечании указывает, что этот раздел статьи о столкновении осаждающейся твёрдой частицы тела с днищем сосуда выполнен в соавторстве с Забабахиным Н.Е.

В статье [2] практически полностью воспроизведена суть работы со ссылкой на оригинальный теоретический и расчётный материал, ранее изложенный в [1]. В графической части статьи, как и в [1], твёрдая частица представлена как цилиндрическая частица с верхним горизонтальным днищем заданного радиуса и нижним сферическим, радиус которого не указан. Высота вертикальной цилиндрической части частицы не наименована и размер её также не указан.

В статье [2] дополнительно приведены результаты первых экспериментальных исследований, качественно подтверждающих теоретический вывод о невозможности в определённых условиях соударения выпуклых тел с ниже расположенной плоской стенкой сосуда. Однако никаких результатов сравнения поведения сферических частиц с расчётными оценками не приводится.

МЕТОД ПРОВЕДЕНИЯ ИССЛЕДОВАНИЯ

На рис. 1 показана принятая в работах [1, 2] форма частицы в виде цилиндра радиуса $R$ с верхним плоским днищем и с выпуклым нижним днищем радиуса $a$, укажем дополнительно высоту цилиндрической части $H$. Частица приближается к плоскому днищу сосуда со скоростью $\upsilon$.

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

Fig. 1. Configuration of collocation of a precipitating solid particle and a flat surface of a vessel head.

Ширина зазора между нижней поверхностью частицы и плоской стенкой $h=Y+r^2/\left(2R\right)$ со временем уменьшается. Таким образом, предполагается, что частица имеет форму цилиндра с плоской верхней и сферической нижней поверхностями. Предполагается также, что режим течения в кольцевом сечении между перемещающейся вниз частицей и плоской стенкой ламинарный. Распределение скорости $\upsilon$ по толщине зазора принимается параболическим, а при данном $r$ средняя по высоте скорость $u_{ср}=(2/3)u_{\max }$, а градиент скорости у поверхности равен $4u_{\max }/h$ или $6u_{\text{ср}}/h$.

Далее рассматривается элемент жидкого диска и действующие на него силы давления $P$ и трения $F$, сумма которых должна быть равна нулю. В результате определяется распределение давления в окружающей жидкости. Далее с учётом давления жидкости на торцы выделенного цилиндрического элемента жидкости определяется суммарная замедляющая сила и записывается уравнение движения осаждающего тела [1, 2]:

$M\frac{d\upsilon }{dt}=-\frac{3}{2}\pi \eta \upsilon \frac{a^4}{Y\cdot {\left(Y+a^2/2R\right)}^2}$. (1)

В уравнении (1) величина $M$ обозначает эффективную массу тела, превосходящую его истинную массу на величину присоединённой массы жидкости. Эта масса прибавляется к массе тела, движущегося неравномерно в жидкой среде для учёта воздействия среды на это тело [3]. Величина присоединённой массы зависит от формы тела, направления движения и плотности среды.

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

Поскольку $\frac{d\upsilon }{dt}=\upsilon \cdot \frac{d\upsilon }{dY}$, и, полагая, что $Y\ll (a^2/2R)$, выражение (1) упрощаем, и оно принимает вид:

$\frac{d\upsilon }{dY}=-6\pi \eta \cdot {\left(\frac{R^2}{MY}\right)}^2$. (2)

Интегрируя последнее выражение с граничным условием, что при $Y=Y_0$, $\upsilon ={\upsilon }_0$, получим

$\upsilon -{\upsilon }_0=-6\pi \eta \cdot \frac{R^2}{M}\cdot \ln \left(\frac{Y}{Y_0}\right)$. (3)

Полагая, что частица остановилась $\upsilon =0$, выражение (3) преобразуем в соотношение

$Y_k=Y_0\exp \left(\frac{M{\upsilon }_0}{6\pi \eta R^2}\right)$. (4)

В работах [1, 2] приводится один и тот же пример расчёта: $M$=1 мг, $R$=0,05 см, $Y_0$=1 мм, $\eta$=0,01 г/см· с при ${\upsilon }_0$=– 1 см/с. В обеих статьях расчёт по формуле (4) даёт значение $Y_k$=0,012 см, т. е. частица не касается днища сосуда.

Кроме того, в работе [1] дополнительно приводится второй пример расчёта с отличающимися исходными данными: $M$=1 г, $R$=0,5 см, $Y_0$=1 см, $\eta$=0,01 г/см·с при ${\upsilon }_0$= –1 см/с. Расчёт по формуле (4) даёт значение $Y_k$=10–8 см, что много меньше ожидаемой шероховатости поверхности, т. е. частица практически касается днища сосуда [1].

В этих двух примерах, несмотря на значительное различие в массах частиц, в тысячу раз, начальная скорость ${\upsilon }_0$= –1 см/с принята одинаковой, что традиционно требует специального обоснования, которое отсутствует.

УТОЧНЕНИЕ МЕТОДИКИ РАСЧЁТА ПРИДОННОГО ОСАЖДЕНИЯ ЧАСТИЦ

Вместе с тем использование допущения $Y\ll (a^2/2R)$ не является препятствием для проведения интегрирования уравнения (1). Проведя интегрирование дифференциального уравнения (1) в пределах от $Y_0$ до $Y_k$, получим выражение:

$\upsilon -{\upsilon }_0=6\pi \eta \cdot \frac{R^2}{M}\cdot \left[\ln \left(1+\frac{a^2}{2RY}\right)+\frac{1}{1+a^2/\left(2RY\right)}-\ln \left(1+\frac{a^2}{2R{Y}_o}\right)-\frac{1}{1+a^2/\left(2R{Y}_0\right)}\right]$ (5)

Полагая, что частица остановилась при $\upsilon =0$, точнее близком к нулю, выражение (5) преобразуем в соотношение:

${\upsilon }_0=6\pi \eta \cdot \frac{R^2}{M}\cdot \left[\ln \left(1+\frac{a^2}{2R{Y}_0}\right)+\frac{1}{1+a^2/\left(2R{Y}_0\right)}-\ln \left(1+\frac{a^2}{2R{Y}_k}\right)-\frac{1}{1+a^2/\left(2R{Y}_k\right)}\right]$ (6)

В сравнении с предшествующим вариантом, рассмотренном в работах [1, 2] (см. выражение (4)), возникла зависимость от формы частицы, которая в [2] была принята в виде цилиндра радиуса a c выпуклым дном радиуса $R$, обращенным в сторону днища сосуда.

В табл. 1 представлены результаты расчётов, привёденных в работах [1, 2] и полученных в настоящей работе. Во всех расчётах динамическая вязкость воды равна $\eta$=0,01 г/см·с

Таблица 1. Сравнение расчётных характеристик осаждения твёрдых частиц вблизи горизонтального днища

Table 1. Comparison of calculated properties of solid particles precipitation near a horizontal vessel head

Сравнение результатов расчётов, выполненных в работах [1, 2], с расчётом с использованием соотношений, полученных в настоящей работе, показывает, что весьма значительно уточняется значение конечного расстояния $Y_k$ от вершины тела до неподвижной стенки. Кроме того, появилась возможность определить локальное значение скорости частиц в ${\upsilon }_k$ см/с в этих точках.

Отметим также, что расчётные значения величины $Y_k$ находятся на уровне микронов и менее, что соответствует шероховатости поверхности, достигаемой при её шлифовании [4], и представляются физически недостижимыми в реальных условиях заполненного жидкостью сосуда.

С использованием полученных результатов оценим приемлемость допущений, использованных в работах [1, 2], а именно условие

$Y\ll (a^2/2R)$. (7)

В работах [1, 2] не содержится никакой информации о величине радиуса частицы . Поэтому при проведении расчётов в рамках используемой модели принималось, что радиусы цилиндрической части частицы и радиус её выпуклой донной части равны, т. е. $a=R$. В этом случае условие (7) преобразуется в соотношение

$Y\ll R/2$. (8)

Расчётные данные, проведённые в табл. 1, показывают, что условие (8) выполняется во всех случаях при расположении частиц вблизи нижней неподвижной поверхности. Условие (8) не выполняется в исходном состоянии с использованием соотношений приведённых в [1, 2].

Ранее выполненные исследования [5] показали, что в случае твёрдых частиц произвольной формы целесообразно использовать три характерных диаметра частицы: $d_v^3=6a^2(H+2a/3)$ — куб диаметра шара, эквивалентного объёму частицы; $d_s^2=a(3a+2H)$ — квадрат диаметра шара, эквивалентного по площади боковой поверхности частицы; $d_m=2a$ — диаметр шара, эквивалентного по площади миделевому сечению частицы. В выражении для числа Рейнольдса используется эквивалентный диаметр частицы $d_e=(2d_s+d_m)/3$. С использованием таких представлений рассчитывались силы гидравлического сопротивления для твёрдых частиц, выполненных в виде двух сферических частиц, диске, иглоподобных эллипсоидах. Во всех случаях для двух различных ориентаций частиц относительно вертикального направления осаждения частиц получено удовлетворительное соответствие с результатами точных численных расчётов [5]. Этот подход был распространён и на осаждение полимодальных двухфазных смесей. В этом случае учитывалось также столкновение между частицами разных фракций твёрдой фазы. В дальнейшем такой подход использовался для расчёта движения полимодальных смесей твёрдых частиц в горизонтальных, вертикальных и наклонных трубопроводах [6]. При этом дополнительно учитывался процесс диффузии твёрдых частиц в вертикальном направлении.

Во всех случаях при исследовании осаждающихся твёрдых частиц и движении их в потоках жидкости необходима полная информация о фактической геометрической конфигурации твёрдых частиц. В работах [1, 2] фигурирует лишь один непосредственно геометрический параметр — радиус нижнего выпуклого днища твёрдой частицы $R$. Величина $M$ определяет массу частицы, но без указания плотности вещества частицы, что не позволяет определить её объём и далее определить, скорее даже оценить, величины $a$ и $H$. Поэтому, к сожалению, не представляется возможным корректно рассчитать начальное значение скорости ${\upsilon }_0$.

Скорость осаждения твёрдой частицы произвольной формы при использовании величины коэффициента гидравлического сопротивления по формуле Стокса для сферических частиц можно представить в виде зависимости [5]:

${\upsilon }_0=({\rho }_s-\rho )\cdot \frac{g{d}_v^3{d}_e}{18\eta d_m^2}$. (9)

Положим так же, как было принято в [1, 2], что при осаждении крупных части $a=H=R=0,5 см$, а в случае мелких частиц — $a_{*}=H_{*}=R_{*}=\mathrm{0,5}\text{ мм}$.

Используя выражение (9), определим отношение скоростей осаждения крупных частиц к скоростям осаждения мелких частиц.

$\frac{{\upsilon }_{01}}{{\upsilon }_{02}}={\left(\frac{d_{v1}}{d_{v2}}\right)}^3\cdot \left(\frac{{d_e}_1}{d_{e2}}\right)\cdot {\left(\frac{d_{m2}}{{d_m}_1}\right)}^2$, (10)

где ${d_{v1}}^3=6a^2(H+2a/3)=\mathrm{1,25}{\text{ см}}^{\text{3}}$; ${d_v}_2^3=6{a_{*}}^2(H_{*}+2a_{*}/3)=\mathrm{1,25}{\text{ мм}}^{\text{3}}$; $d_{m1}=2a=1\text{ см}$; $d_{s1}={\left[a\right(3a+2H\left)\right]}^{1/2}=\mathrm{1,118}\text{ см}$; $d_{e1}=(2d_{s1}+d_{m1})/3=\mathrm{1,0787}\text{ см}$; $d_{s2}={\left[a_{*}\right(3a_{*}+2H_{*}\left)\right]}^{1/2}=\mathrm{1,118}\text{ мм}$; $d_{m2}=2a_{*}=2\times \mathrm{0,5}=1\text{ мм}$; $d_{e2}=(2d_{s2}+d_{m2})/3=\mathrm{1,0787}\text{ мм}$.

Подставляя приведённые численные значения в выражение (10), получим, что отношение ${\upsilon }_{01}/{\upsilon }_{02}=100$, т. е. скорости осаждения частиц значительно различаются, а не равны, как предполагается в работах [1, 2].

Применительно к сферическим осаждающимся частицам такой же вывод следует и из известной формулы Стокса для осаждающихся сферических частиц: ${\upsilon }_0=({\rho }_s-\rho )\frac{g{d}^2}{18\eta }$. В этом случае при различии диаметра частиц в 10 раз скорости осаждения также различаются в 100 раз.

При определении скорости свободного осаждения твёрдых частиц, как правило, не учитывается влиянием Броуновского движения на скорость свободного осаждения частиц. Как отмечается в [7], Броуновское движение оказывает важное влияние на двухфазные потоки, включающие мельчайшие твёрдые частицы. Если исходить из условия, что кинетическая энергия осаждающихся частиц много больше энергии теплового движения молекул, то условие, что Броуновским движением можно пренебречь, записывается в виде [5].

$\frac{1}{2}\cdot \frac{\pi d^3}{6}\cdot \rho {\upsilon }^2>>\frac{3}{2}\cdot kT$, (11)

где $k=\mathrm{1,38}\times {10}^{-23}\text{ Дж/К}$ — постоянная Больцмана, $T$ — температура в градусах Кельвина.

Характерное значение диаметра частиц, при котором необходимо учитывать влияние Броуновского движения на скорость осаждения твёрдых частиц, определяется выражением [5].

$d\gg {\left[\frac{5832kT{\nu }^2}{\pi {\left(\frac{{\rho }_s}{\rho }-1\right)}^2{g}^2{\rho }_s}\right]}^{\frac{1}{7}}$. (12)

Расчёты величины правой части в формуле (10) при различных значениях плотности материала твёрдых частичек в воде при $T=293°K$ представлены в табл. 2.

Таблица 2. Значения критического диаметра частиц твёрдой фазы

Table 2. Values of critical diameter of solid particles

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

ЗАКЛЮЧЕНИЕ

Проведённый анализ показал, что использование в работах [1, 2] допущения, что $Y\ll \left(a^2/2R\right)$, является физически некорректным при анализе осаждения твёрдых частиц вблизи плоского днища. Показано, что возможно получение аналитического решения без привлечения указанного допущения, что позволяет более точно определить координаты частицы вблизи плоского днища.

ДОПОЛНИТЕЛЬНО

Вклад авторов. А.С. Кондратьев ― постановка задачи исследования, анализ публикаций по теме статьи, разработка обобщённой расчётной модели, анализ результатов расчётов, редактирование текста рукописи; П.П. Швыдько — библиографический поиск, проведение численных исследований, подготовка текста рукописи, создание изображений. Авторы подтверждают соответствие своего авторства международным критериям ICMJE (все авторы внесли существенный вклад в разработку концепции, проведение исследования и подготовку статьи, прочли и одобрили финальную версию перед публикацией).

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

Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования.

Об авторах

Александр Сергеевич Кондратьев

Автор, ответственный за переписку.
Email: ask41@mail.ru
ORCID iD: 0009-0001-2030-5165
SPIN-код: 6739-0496

д-р техн. наук, профессор кафедры «Промышленная теплоэнергетика»

Павел Петрович Швыдько

Email: Shvydko.P@gmail.com
ORCID iD: 0009-0007-8030-6403
SPIN-код: 8347-3886

руководитель проекта

Список литературы

Дополнительные файлы

Доп. файлы

Действие

1. JATS XML

Скачать

2. Рис. 1. Конфигурация взаиморасположения осаждающейся твёрдой частицы и плоской поверхности днища сосуда.

Скачать (52KB)

Метаданные