Теоретические основы химической технологии, 2023, T. 57, № 4, стр. 408-411

Полуэмпирический метод определения скоростей фаз при псевдоожижении и осаждении в автомодельной области

А. М. Трушин a*, М. А. Носырев a, Л. В. Равичев a, С. И. Фролова a, Л. С. Сальникова a, В. Е. Яшин a

a Российский химико-технологический университет им. Д.И. Менделеева
Москва, Россия

* E-mail: vnissok@list.ru

Поступила в редакцию 30.03.2023
После доработки 10.04.2023
Принята к публикации 20.04.2023

Полный текст (PDF)

Аннотация

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

Ключевые слова: безразмерная скорость осаждения, псевдоожижения, порозность, сферические, и несферические частицы, автомодельная область движения

ВВЕДЕНИЕ

В расчетах аппаратуры процессов со взвешенным слоем безразмерная приведенная скорость жидкости (отношение приведенной скорости жидкости к скорости витания) в однородном псевдоожиженном слое сферических и несферических частиц в автомодельной области движения (Re > 500) определяется по эмпирическим степенным зависимостям от порозности слоя полученным экспериментальным путем без теоретических предпосылок [1]. Безразмерная скорость стесненного осаждения сферических частиц (отношение скорости стесненного движения к скорости одиночной частицы)в поле силы тяжести в силу аналогии осаждения и псевдоожижения может определяться по тем же степенным зависимостям от порозности.

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

СКОРОСТЬ ОСАЖДЕНИЯ СФЕРИЧЕСКИХ ЧАСТИЦ

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

(1)
$\frac{3}{4}{{\xi }_{0}}\frac{{\rho w_{0}^{2}}}{d} = \left( {{{\rho }_{{\text{д}}}} - {{\rho }_{{\text{с}}}}} \right)g.$

При составлении баланса сил для стесненного осаждения примем следующие допущения:

1) порозность равномерно распределена по слою дисперсных частиц;

2) сопротивление, испытываемое группой частиц равно сумме сопротивлений, рассчитанных для каждой частицы подобным методом, как и для одиночной частицы. При этом коэффициент сопротивления ${{\xi }_{0}}$ заменяется на ${{\xi }_{{{\text{ст}}}}}$, учитывающий изменение поля скоростей и давлений, а также столкновение частиц при стесненном движении. Скорость движения частицы относительно стенок аппарата w0 заменяется на относительную скорость фаз wот;

3) выталкивающая сила Архимеда становится пропорциональной плотности слоя, а не плотности сплошной фазы как в случае одиночной частицы.

Из второго и третьего допущений получим баланс сил для каждой частицы слоя:

(2)
${{\xi }_{{{\text{ст}}}}}\frac{{\pi {{d}^{2}}}}{4}\frac{{w_{{{\text{от}}}}^{2}\rho }}{2} = \frac{{\pi {{d}^{3}}}}{6}\left( {{{\rho }_{{\text{д}}}} - {{\rho }_{{{\text{сл}}}}}} \right)g.$

Относительная скорость фаз может быть выражена через скорость движения частиц относительно стенок аппарата и порозность:

(3)
${{w}_{{{\text{от}}}}} = \frac{{{{w}_{{{\text{ст}}}}}}}{\varepsilon }.$

Разницу плотностей частиц и слоя можно представить следующим образом:

(4)
$\left( {{{\rho }_{{\text{д}}}} - {{\rho }_{{{\text{сл}}}}}} \right) = \varepsilon \left( {{{\rho }_{{\text{д}}}} - {{\rho }_{{\text{с}}}}} \right).$

Из уравнений (2)–(4) получим:

(5)
$\frac{3}{4}{{\xi }_{{{\text{ст}}}}}\frac{{w_{{{\text{ст}}}}^{2}\rho }}{{{{\varepsilon }^{2}}d}} = \varepsilon \left( {{{\rho }_{{\text{д}}}} - {{\rho }_{{\text{с}}}}} \right)g.$

Почленное деление уравнений (1) и (5) приводит к следующему выражению для безразмерной скорости стесненного осаждения:

(6)
$\frac{{{{w}_{{{\text{ст}}}}}}}{{{{w}_{0}}}} = {{\left( {{{\varepsilon }^{3}}\frac{{{{\xi }_{0}}}}{{{{\xi }_{{{\text{ст}}}}}}}} \right)}^{{0.5}}}.$

Эмпирическая зависимость, связывающая безразмерную приведенную скорость жидкости в однородном псевдоожиженном слое с порозностью в автомодельной области, имеет следующий вид [1]:

(7)
$\frac{{{{w}_{{{\text{пр}}}}}}}{{{{w}_{{{\text{вит}}}}}}} = {{\varepsilon }^{{2.4}}}.$

В силу аналогии однородного псевдожижения и осаждения эта зависимость сохраняется и для безразмерной скорости осаждения:

(8)
$\frac{{{{w}_{{{\text{ст}}}}}}}{{{{w}_{0}}}} = {{\varepsilon }^{{2.4}}}.$

Ввиду того, что коэффициент сопротивления для одиночной частицы (${{\xi }_{0}}$) применяется в расчетах постоянной величиной (в интервале 0.43–0.5), коэффициент сопротивления при осаждении в стесненных условиях согласно уравнениям (6) и (8) должен выражаться следующим образом:

(9)
${{\xi }_{{{\text{ст}}}}} = \frac{{{{\xi }_{0}}}}{{{{\varepsilon }^{{1.8}}}}}.$

Полученное выражение для коэффициента сопротивления при движении в стесненных условиях очевидно справедливо лишь при равномерном распределении порозности в слое осаждающихся частиц (первое допущение). Поскольку есть многочисленные данные о том, что при осаждении, всплытии пузырьков и псевдоожижении наблюдается групповое движение частиц дисперсной фазы [16], причем локальная объемная доля дисперсных частиц в группах (φлок) превышает среднюю по слою долю (φ), а локальная порозность соответственно меньше средней, уравнения (6) и (9) для группового движения неприменимы.

По данным работы [3] среднее соотношение локальной и средней доли твердых частиц (В) при осаждении в интервале значений φ (0–0.6) составляет 1.14. Принимая приближенно среднее соотношение локальной и средней доли твердой фазы постоянным для всех величин порозности, получим выражение для локальной порозности в осаждающихся группах частиц.

(10)
${{\varepsilon }_{{{\text{лок}}}}} = \left[ {1 - \left( {1 - \varepsilon } \right)1.14} \right].$

Для группового осаждения частиц уравнение (6) сохраняет свой вид при замене средней порозности на локальную. В связи с этим возникает необходимость найти зависимость ξст от локальной порозности. Поскольку общие закономерности осаждения в слое немного более концентрированным по твердой фазе не претерпят сильных изменений, можно предположить, что зависимость ξст от локальной порозности также имеет вид степенной функции:

(11)
${{\xi }_{{{\text{ст}}}}} = \frac{{{{\xi }_{0}}}}{{\varepsilon _{{{\text{лок}}}}^{n}}}.$

Уравнение (11) также как уравнение (9) соответствуют граничному условию при переходе от стесненного осаждения к свободному: ξст = ξ0 при ε = εлок = 1.

Показатель степени n в уравнении (11) можно найти по методу наименьших квадратов:

(12)
$\frac{d}{{dB}}\left\{ {{{{\int\limits_{0.44}^1 {\left[ {{{\varepsilon }^{{2.4}}} - {{{\left[ {1 - \left( {1 - \varepsilon } \right)B} \right]}}^{{\left( {3 + n} \right) \times 0.5}}}} \right]} }}^{2}}d\varepsilon } \right\} = 0.$

Решением уравнения (12) относительно В при различных значениях показателя степени n позволило найти величину n равной единице при условии равенства значения В (с небольшим приближением) ранее найденному соотношению локальной и средней доли твердой фазы (В = 1.14).

Таким образом, из уравнения (6), где средняя порозность заменяется на локальную, а также уравнений (8), (12) получен следующий результат: безразмерная скорость осаждения сферических частиц в автомодельной области равна локальной порозности в квадрате:

(13)
$\frac{{{{w}_{{{\text{ст}}}}}}}{{{{w}_{0}}}} = {{\left[ {1 - \left( {1 - \varepsilon } \right) \times 1.14} \right]}^{2}}.$

Ввиду того, что уравнение (13) получено при приближенном равенстве отношений локальной доли твердой фазы к средней, полученных из уравнения (12), а также с учетом допущения о выражении коэффициента сопротивления ξст в виде степенной функции от локальной порозности, было проведено сравнение величин безразмерных скоростей осаждения рассчитанных по уравнению (13) и уравнению (8) полученному обработкой экспериментальных данных. Результаты сравнительных расчетов приведены в табл. 1. Коэффициенты детерминации, найденные при сравнении расчетов по уравнениям (8) и (13) составляют 0.9997.

Таблица 1.  

Зависимость безразмерной скорости осаждения от порозности слоя

Номер смеси Форма частиц Расчетная формула Порозность, ε Коэффициент детерминации
0.44 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95
1 Сфера 8 0.1394 0.1895 0.2382 0.2930 0.3556 0.4250 0.5013 0.5853 0.6770 0.7765 0.8840 0.9997
Сфера 13 0.1308 0.1849 0.2372 0.2960 0.3612 0.4330 0.5112 0.5960 0.6872 0.7850 0.8890
2 Сфера 13 0.1308 0.1849 0.2372 0.2960 0.3612 0.4330 0.5112 0.5960 0,6872 0.7850 0.8890 0.9994
Куб 14 0.1317 0.1806 0.2285 0.2883 0.3452 0.4145 0.4915 0.5764 0.6695 0.7709 0.8810
3 Сфера 13 0.1308 0.1849 0.2372 0.2960 0.3612 0.4330 0.5112 0.5960 0.6872 0.7850 0.8890 0.9996
Цилиндр L/d = 2 14 0.1367 0.1863 0.2348 0.2889 0.3520 0,4213 0.4979 0.5822 0,6744 0.7746 0.8830
4 Сфера 13 0.1308 0.1849 0.2372 0.2960 0.3612 0.4330 0.5112 0.5960 0.6872 0.7850 0.8890 0.9997
Цилиндр L/d = 2.5 14 0.1394 0.1895 0.2387 0.2930 0.3556 0.4250 0.5013 0.5863 0.6770 0.7765 0,8840
5 Сфера 13 0.1308 0.1849 0.2372 0.2960 0.3612 0.4330 0.5112 0.5960 0.6872 0.7850 0.8890 0.9997
Цилиндр L/d=3 14 0.1413 0.1915 0.2404 0.2959 0.3581 0.4273 0.5040 0.5874 0.6782 0.7779 0.8849
6 Сфера 13 0.1308 0.1849 0.2372 0.2960 0.3612 0.4330 0.5112 0.5960 0.6872 0.7850 0.8890 0.9998
Цилиндр L/d = 4 14 0.1437 0.1943 0.2435 0.2990 0.3613 0.4305 0.5067 0.5902 0.6810 0.7796 0.8860
7 Куб 14 0.1317 0.1806 0.2285 0.2883 0.3452 0.4145 0.4915 0.5764 0.6695 0.7709 0.8810 0.9996
Куб 16 0.1228 0.1764 0.2285 0.2873 0.3528 0.4251 0.5041 0.5898 0.6822 0.7814 0.8870

СКОРОСТЬ ОСАЖДЕНИЯ НЕСФЕРИЧЕСКИХ ЧАСТИЦ

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

(14)
$\frac{{{{w}_{{{\text{ст}}}}}}}{{{{w}_{0}}}} = {{\varepsilon }^{m}}.$

Показатель степени m определяется по следующей формуле:

(15)
$m = 2.7{{\left[ {\frac{\pi }{6}{{{\left( {\frac{{{{d}_{s}}}}{{{{d}_{p}}}}} \right)}}^{3}}} \right]}^{{0.16}}},$
где ds – диаметр сферы с поверхностью равной поверхности несферической частицы; dp – диаметр круга, равного по площади миделевому сечению частицы.

Уравнение (13) полученное для сферических частиц, может быть использовано и для несферических (кубических и цилиндрических) при условии равенства коэффициентов формы (отношений скорости осаждения несферических и сферических частиц для свободного и стесненного осаждения). В этом случае при почленном делении уравнений (1) и (5) получится тоже уравнение (6). Кроме того, необходимо, чтобы коэффициент В оставался прежним (1.14). Эти условия достаточно хорошо выполняются для цилиндрических частиц при соотношениях длины и диаметра больше двух. Результаты безразмерных скоростей, рассчитанных по уравнениям (13) и (14) при этом условии достаточно близки (коэффициенты детерминации составляют 0.9996–0.9998).

В случае кубических частиц величина В, характеризующая объединение частиц в группы должна быть больше 1.14, так как при столкновении частиц кубической формы возможно их объединение по плоскостям граней, поэтому коэффициент 1.14 для этого случая в уравнении (13) заменен на 1.16. Таким образом, для кубических частиц предлагается следующее выражение:

(16)
$\frac{{{{w}_{{{\text{ст}}}}}}}{{{{w}_{0}}}} = {{\left[ {\left( {1 - \left( {1 - \varepsilon } \right) \times 1.16} \right)} \right]}^{2}}.$

Сравнению безразмерных скоростей, определенных по уравнениям (14) и (16) соответствует коэффициент детерминации 0.9996.

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

(17)
$\frac{{{{w}_{{{\text{ст}}}}}}}{{{{w}_{0}}}} = \frac{{{{w}_{{{\text{пр}}}}}}}{{{{w}_{{{\text{вит}}}}}}} = {{\left[ {\left( {1 - \left( {1 - \varepsilon } \right)B} \right)} \right]}^{2}},$
где В = 1.14 для сферических и цилиндрических частиц при указанных выше соотношениях L/d; В = 1.16 для кубических частиц.

ЗАКЛЮЧЕНИЕ

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

ОБОЗНАЧЕНИЯ

ξ0 коэффициент сопротивления при свободном осаждении
ξст коэффициент сопротивления при стесненном осаждении
ε средняя доля сплошной фазы (порозность), м33
εлок локальная порозность, м33
g ускорение свободного падения, м/с2
w0 скорость движения одиночной частицы, м/с
wст скорость стесненного движения частиц дисперсной фазы, м/с
wот относительная скорость фаз, м/с
ρд плотность частицы, кг/м3
ρс плотность сплошной фазы, кг/м3
ρсл плотность слоя, кг/м3
n показатель степени
В отношение локальной доли дисперсной фазы к средней доле, м33

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

  1. Davidson J., Harrison F. Fluidization. Academic Press. London and New York. 1971. P. 847.

  2. Couper J.R., Penney W.R., Fair J.R., Walas S.M. // Gulf professional publishing. 2005. P. 801

  3. Trushin A.M., Dmitriev E.A., Nosyrev M.A. The Use of the Variation Method to Define the Velocity of Constrained Turbulent Motion of Disperse Particles in Liquid // Theor. Found. Chem. Eng. 2017. V. 51. № 2. P. 155–158. [Трушин А.М., Дмитриев Е.А., Носырев М.А. Использование вариационного метода для определения скорости стесненного турбулентного движения дисперсных частиц в жидкости // Теорет. основы хим. технологии. 2017. Т. 51. № 2. С. 86–89].

  4. Giusti A., Lucci F., Soldati A. Influence of the lift force in direct numerical simulation of upward/downward turbulent channel flow laden with surfactant contaminated microbubbles // Chemical Engineering Science V. 60. № 6. P. 6176–6187

  5. Bakker P.J., Heerties P.M. Porosity distributions in a fluidized bed // Chemical Engineering Science V. 12. № 4. 1960. P. 260–271

  6. Равичев Л.В., Беспалов А.В. К расчету порозности псевдоожиженного зернистого катализаторного слоя // Химическая промышленность сегодня. 2007. № 1. С. 4–9.

Дополнительные материалы отсутствуют.