Теоретические основы химической технологии, 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
- EDN: VJTNVR
- DOI: 10.31857/S0040357123040140
Аннотация
Получено выражение для определения безразмерной скорости в процессах осаждения и псевдоожижения для сферических, цилиндрических и кубических частиц в автомодельной области.
ВВЕДЕНИЕ
В расчетах аппаратуры процессов со взвешенным слоем безразмерная приведенная скорость жидкости (отношение приведенной скорости жидкости к скорости витания) в однородном псевдоожиженном слое сферических и несферических частиц в автомодельной области движения (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.$Относительная скорость фаз может быть выражена через скорость движения частиц относительно стенок аппарата и порозность:
Разницу плотностей частиц и слоя можно представить следующим образом:
(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]:
В силу аналогии однородного псевдожижения и осаждения эта зависимость сохраняется и для безразмерной скорости осаждения:
Ввиду того, что коэффициент сопротивления для одиночной частицы (${{\xi }_{0}}$) применяется в расчетах постоянной величиной (в интервале 0.43–0.5), коэффициент сопротивления при осаждении в стесненных условиях согласно уравнениям (6) и (8) должен выражаться следующим образом:
Полученное выражение для коэффициента сопротивления при движении в стесненных условиях очевидно справедливо лишь при равномерном распределении порозности в слое осаждающихся частиц (первое допущение). Поскольку есть многочисленные данные о том, что при осаждении, всплытии пузырьков и псевдоожижении наблюдается групповое движение частиц дисперсной фазы [1–6], причем локальная объемная доля дисперсных частиц в группах (φлок) превышает среднюю по слою долю (φ), а локальная порозность соответственно меньше средней, уравнения (6) и (9) для группового движения неприменимы.
По данным работы [3] среднее соотношение локальной и средней доли твердых частиц (В) при осаждении в интервале значений φ (0–0.6) составляет 1.14. Принимая приближенно среднее соотношение локальной и средней доли твердой фазы постоянным для всех величин порозности, получим выражение для локальной порозности в осаждающихся группах частиц.
(10)
${{\varepsilon }_{{{\text{лок}}}}} = \left[ {1 - \left( {1 - \varepsilon } \right)1.14} \right].$Для группового осаждения частиц уравнение (6) сохраняет свой вид при замене средней порозности на локальную. В связи с этим возникает необходимость найти зависимость ξст от локальной порозности. Поскольку общие закономерности осаждения в слое немного более концентрированным по твердой фазе не претерпят сильных изменений, можно предположить, что зависимость ξст от локальной порозности также имеет вид степенной функции:
Уравнение (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]:
Показатель степени m определяется по следующей формуле:
(15)
$m = 2.7{{\left[ {\frac{\pi }{6}{{{\left( {\frac{{{{d}_{s}}}}{{{{d}_{p}}}}} \right)}}^{3}}} \right]}^{{0.16}}},$Уравнение (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.
Исходя из вышеизложенного, безразмерные скорости частиц при осаждении и жидкости при псевдоожижении могут быть определены по уравнению:
где В = 1.14 для сферических и цилиндрических частиц при указанных выше соотношениях L/d; В = 1.16 для кубических частиц.ЗАКЛЮЧЕНИЕ
С учетом неравномерного распределения порозности в слое осаждающихся и псевдоожиженных частиц, безразмерная скорость стесненного осаждения и безразмерная приведенная скорость жидкости при псевдоожижении могут быть приравнены локальной порозности в квадрате. Такая полуэмпирическая зависимость найдена для сферических, кубических и цилиндрических частиц.
ОБОЗНАЧЕНИЯ
ξ0 | коэффициент сопротивления при свободном осаждении |
ξст | коэффициент сопротивления при стесненном осаждении |
ε | средняя доля сплошной фазы (порозность), м3/м3 |
εлок | локальная порозность, м3/м3 |
g | ускорение свободного падения, м/с2 |
w0 | скорость движения одиночной частицы, м/с |
wст | скорость стесненного движения частиц дисперсной фазы, м/с |
wот | относительная скорость фаз, м/с |
ρд | плотность частицы, кг/м3 |
ρс | плотность сплошной фазы, кг/м3 |
ρсл | плотность слоя, кг/м3 |
n | показатель степени |
В | отношение локальной доли дисперсной фазы к средней доле, м3/м3 |
Список литературы
Davidson J., Harrison F. Fluidization. Academic Press. London and New York. 1971. P. 847.
Couper J.R., Penney W.R., Fair J.R., Walas S.M. // Gulf professional publishing. 2005. P. 801
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].
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
Bakker P.J., Heerties P.M. Porosity distributions in a fluidized bed // Chemical Engineering Science V. 12. № 4. 1960. P. 260–271
Равичев Л.В., Беспалов А.В. К расчету порозности псевдоожиженного зернистого катализаторного слоя // Химическая промышленность сегодня. 2007. № 1. С. 4–9.
Дополнительные материалы отсутствуют.
Инструменты
Теоретические основы химической технологии