Теоретические основы химической технологии, 2021, T. 55, № 2, стр. 261-264

К вопросу о расчете скорости начала псевдоожижения

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

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

* E-mail: vnissok@list.ru

Поступила в редакцию 09.10.2020
После доработки 20.10.2020
Принята к публикации 20.11.2020

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

Аннотация

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

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

ВВЕДЕНИЕ

В наиболее распространенном методе расчета скорости начала псевдоожижения (νкр) перепад давления в слое, рассчитанный по уравнению Эргуна, приравнивается весу частиц слоя, отнесенному к площади поперечного сечения слоя, за вычетом подъемной силы Архимеда, пропорциональной плотности сплошной фазы. В уравнении для перепада давления используется порозность, близкая к порозности неподвижного слоя. Такой метод расчета, который может считаться стандартным, приводится в большинстве литературных источников, касающихся исследования процесса псевдоожижения [14], а также в учебной литературе последних изданий [58]. К сожалению, стандартный метод расчета содержит некоторые противоречия. Прежде всего, противоречием является то, что в баланс сил в начале псевдоожижения для сферических частиц подставляется порозность, равная 0.4, тогда как полное псевдоожижение начинается при несколько расширенном слое, при этом порозность увеличивается приблизительно на 10% до величины ε = 0.44 [3]. Это противоречие исчезает, если принять, что скорость начала псевдоожижения соответствует порозности слоя сферических частиц (0.4), которая превышает среднюю порозность неподвижного слоя сфер (0.38). Это свидетельствует о том, что при таком подходе начало псевдоожижения соответствует несколько расширенному слою по сравнению с неподвижным и отличается от начала полного псевдоожижения. Такое определение скорости начала псевдоожижения приближенно соответствует методу его определения по пересечению касательных на графике ∆Р = f(ν). Существенным противоречием является принятие коэффициентов сопротивления при расширении слоя до полного псевдоожижения постоянными, тогда как при скоростях, значительно меньших скорости начала псевдоожижения, наблюдается движение частиц [1]. Увеличение порозности слоя, сопровождающееся движением частиц, влияет на величину коэффициентов сопротивления.

РАСЧЕТ СКОРОСТИ НАЧАЛА ПСЕВДООЖИЖЕНИЯ И ПОЛНОГО ПСЕВДООЖИЖЕНИЯ

Коэффициенты сопротивления λ в неподвижном слое, принятые в уравнении Эргуна, определяются по следующим формулам:

(1)
$\begin{gathered} {{\lambda }_{1}} = \frac{{133}}{{Re}};\,\,\,\,Re = \frac{{\nu {{d}_{{\text{э}}}}\rho }}{\mu } < 0.{\text{2}};\,\,\,\,{{\lambda }_{{\text{2}}}} = {\text{2}}.{\text{34}}; \\ {\text{Re}} > {\text{7}}000, \\ \end{gathered} $
где ν – скорость в каналах неподвижного слоя; ${{d}_{{\text{э}}}} = \frac{{2\varepsilon d}}{{3\left( {1 - \varepsilon } \right)}}$ – эквивалентный диаметр каналов; d – диаметр сферических частиц.

Коэффициент сопротивления входит в уравнение, определяющее перепад давления в слое:

(2)
$\Delta Р = \lambda \frac{Н}{{{{d}_{{\text{э}}}}}}\frac{{\rho {{\nu }^{2}}}}{2},$
где Н – высота слоя.

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

Таким образом, для ламинарного и турбулентного движения имеем

(3)
${{\lambda }_{1}} = \frac{{{{k}_{1}}\left( \varepsilon \right)}}{{Re}};\,\,\,\,{{\lambda }_{2}} = {{k}_{2}}\left( \varepsilon \right).$

Чтобы найти функцию k1(ε) и k2(ε), на основе уравнений (1), (2) получим выражения для чисел Рейнольдса, рассчитываемых по приведенной скорости в конце стадии предварительного расширения слоя.

В ламинарной и автомодельной области:

(4)
${{\operatorname{Re} }_{{\text{п}}}} = \frac{8}{{9{{k}_{1}}\left( \varepsilon \right)}}\frac{{{{\varepsilon }^{3}}}}{{\left( {1 - \varepsilon } \right)}}{\text{Ar}},$
(5)
${{\operatorname{Re} }_{{\text{п}}}} = {{\left( {\frac{{{{\varepsilon }^{3}}{\text{Ar}}}}{{0.75{{k}_{2}}\left( \varepsilon \right)}}} \right)}^{{0.5}}}.$

С учетом выражений для скоростей свободного движения получим отношение приведенных скоростей к скорости витания:

(6)
${{f}_{1}}\left( \varepsilon \right) = \frac{{{{\nu }_{{\text{п}}}}}}{{{{\nu }_{{\text{в}}}}}} = \frac{{16{{\varepsilon }^{3}}}}{{{{k}_{1}}\left( \varepsilon \right)\left( {1 - \varepsilon } \right)}},$
(7)
${{f}_{2}}\left( \varepsilon \right) = \frac{{{{\nu }_{{\text{п}}}}}}{{{{\nu }_{{\text{в}}}}}} = \frac{1}{{1.73}}{{\left( {\frac{{{{\varepsilon }^{3}}}}{{0.75{{k}_{1}}\left( \varepsilon \right)}}} \right)}^{{0.5}}}.$

Известны эмпирические зависимости соотношения скоростей ${{f}_{1}}\left( \varepsilon \right) = {{\varepsilon }^{n}},$ $n = 4.65 + 20\frac{d}{D}$ (Re < 0.2), где d, D – соответственно диаметр частиц и аппарата. Часто используют показатель степени n = = 4.75, соответствующий соотношению $\frac{d}{D} = 5 \times {{10}^{{ - 3}}}.$

Формула вида, подобного уравнению (6), получена также на основе решения уравнения баланса сил в дифференциальной форме для медленного движения [9].

Таким образом, для ламинарной области имеем

(8)
${{f}_{1}}\left( \varepsilon \right) = {{\varepsilon }^{{4.75}}}.$

Для автомодельной области n = 2.4:

(9)
${{f}_{2}}\left( \varepsilon \right) = {{\varepsilon }^{{2.4}}}.$

Приравняв формулы (6) и (7) к формулам (8) и (9), получим данные для получения зависимостей k(ε).

Обработка полученных данных в виде степенной зависимости в интервале величин порозности 0.4–0.44 приводит к следующим формулам:

ламинарная область

(10)
${{k}_{1}}\left( \varepsilon \right) = 53{{\varepsilon }^{{ - 1}}};$

автомодельная область

(11)
${{k}_{2}}\left( \varepsilon \right) = 0.9{{\varepsilon }^{{ - 1}}}.$

Учитывая незначительное отличие величины порозности 0.4 от средней порозности упаковки сферических частиц 0.38, формулы (10) и (11) могут экстраполироваться до этой величины.

Перепад давления, необходимый для определения скорости начала псевдоожижения, можно найти из баланса сил, действующих на фазы кипящего слоя при ряде допущений, представленных в работах [1, 4]. Эти уравнения элементарны, но необходимость приведения вызвана неоднозначной трактовкой их решения, прежде всего в учебной литературе.

Сумма внешних сил, приходящихся на единицу объема каждой из фаз:

сплошная фаза

(12)
$ - \rho g\vec {i} - \nabla P - {{F}_{{R,c}}}\vec {i} = 0;$

дисперсная (твердая) фаза

(13)
$ - {{\rho }_{{\text{Т}}}}g\vec {i} + \left( {{{\rho }_{{\text{Т}}}}\left( {1 - \varepsilon } \right) + \varepsilon \rho } \right)g\vec {i} + \frac{{{{F}_{{R,c}}}\varepsilon }}{{1 - \varepsilon }}\vec {i} = 0,$
где $\vec {i}$ – единичный вектор, направленный вертикально вверх; FR,c сила сопротивления, приходящаяся на единицу объема сплошной фазы.

Вторая составляющая уравнения (13) является подъемной силой Архимеда, рассчитанной по удельному весу кипящего слоя. Третья составляющая получена из равенства абсолютных величин сил сопротивления, действующих на объем каждой из фаз: ${{F}_{{R,c}}}{{V}_{c}} = {{F}_{{R,{\text{д}}}}}{{V}_{{\text{д}}}};$ $\frac{{{{V}_{{\text{д}}}}}}{{{{V}_{{\text{д}}}} + {{V}_{с}}}} = 1 - \varepsilon .$

Из уравнения (13) получим силу сопротивления:

(14)
${{F}_{{R,c}}} = \left( {1 - \varepsilon } \right)\left( {{{\rho }_{{\text{T}}}} - \rho } \right)g.$

Сила сопротивления, приходящаяся на единицу объема слоя:

(15)
${{F}_{{R,c,{\text{д}}}}} = \varepsilon \left( {1 - \varepsilon } \right)\left( {{{\rho }_{{\text{T}}}} - \rho } \right)g.$

Из уравнений (12) и (14) получим градиент давления:

(16)
$\nabla Р = - \vec {i}\left( {\left( {1 - \varepsilon } \right){{\rho }_{{\text{Т}}}}g + \varepsilon \rho g} \right).$

В случае газовой сплошной фазы, пренебрегая малой величиной второй составляющей правой части (16), найдем

(17)
$\nabla P = \vec {i}\left( {1 - \varepsilon } \right){{\rho }_{{\text{T}}}}g.$

После интегрирования получим перепад давления в газовой фазе:

(18)
$\Delta P = \left( {1 - \varepsilon } \right){{\rho }_{{\text{T}}}}gH.$

Ввиду того что наиболее часто в промышленности псевдоожижающим агентом является газ, для определения скорости начала псевдоожижения и начала полного псевдоожижения будет использовано уравнение (18).

Баланс сил в виде равенства ∆Р, определяемого по уравнениям (2) и (18), приводит к следующему уравнению в безразмерных переменных:

(19)
$\frac{{9{{k}_{1}}\left( {{{\varepsilon }_{{{\text{кр}}}}}} \right)\left( {1 - {{\varepsilon }_{{{\text{кр}}}}}} \right)}}{{8\varepsilon _{{{\text{кр}}}}^{3}}}{{\operatorname{Re} }_{{{\text{кр}}}}} + \frac{{3{{k}_{2}}\left( {{{\varepsilon }_{{{\text{кр}}}}}} \right)\operatorname{Re} _{{{\text{кр}}}}^{2}}}{{4\varepsilon _{{{\text{кр}}}}^{3}}} = {\text{Ar}}.$

С учетом формул (10) и (11) получим

(20)
$\frac{{59.6\left( {1 - {{\varepsilon }_{{{\text{кр}}}}}} \right){{{\operatorname{Re} }}_{{{\text{кр}}}}}}}{{\varepsilon _{{{\text{кр}}}}^{4}}} + \frac{{0.68\operatorname{Re} _{{{\text{кр}}}}^{2}}}{{\varepsilon _{{{\text{кр}}}}^{4}}} = {\text{Ar}}.$

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

(21)
$\operatorname{Re} _{{{\text{кр}}}}^{2} + 52.6{{\operatorname{Re} }_{{{\text{кр}}}}}\, - 0.0376{\text{Ar}} = 0.$

Поскольку уравнение (19) включает в себя зависимость коэффициентов сопротивления от порозности его можно использовать для определения скорости полного псевдоожижения при ориентировочной величине критической порозности, соответствующей концу расширения слоя перед полным псевдоожижением, εкр = 0.44.

В результате получим следующее уравнение для определения скорости начала полного псевдоожижения сферических частиц:

(22)
$\operatorname{Re} _{{{\text{кр}}}}^{2}\, + 49.1\operatorname{Re} {}_{{{\text{кр}}}}\, - 0.0551{\text{Ar}} = 0.$

Уравнение для определения скорости начала псевдоожижения слоя сферических частиц по стандартной методике имеет следующий вид:

(23)
$\operatorname{Re} _{{{\text{кр}}}}^{2}\, + 51.4\operatorname{Re} {}_{{{\text{кр}}}}\, - 0.0366{\text{Ar}} = 0.$

ЗАКЛЮЧЕНИЕ

Коэффициенты полученного уравнения (21) и уравнения (23) стандартной методики незначительно различаются, что свидетельствует о близких величинах скорости начала псевдоожижения, рассчитанных по предложенной и стандартной методике. Учет зависимости коэффициентов сопротивления от порозности позволяет использовать уравнение (20) для определения приближенной величины скорости начала полного псевдоожижения. Расчеты по уравнениям (21) и (22) показывают, что скорость начала полного псевдоожижения может на 54% превышать скорость начала псевдоожижения.

ОБОЗНАЧЕНИЯ

FR,c сила сопротивления слоя, приходящегося на единицу объема сплошной фазы
g ускорение свободного падения, м/с2
$\vec {i}$ единичный вектор
ε доля сплошной фазы (порозность)
εкр критическая величина порозности
λ коэффициент сопротивления
ρ плотность сплошной фазы, кг/м3
ρТ плотность частицы, кг/м3
Ar число Архимеда
Reкр число Рейнольдса в критической точке

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

  1. Davidson J., Harrison F. Fluidization. London: Academic, 1971.

  2. Гельперин Н.И., Айнштейн В.Г., Кваша В.Б. Основы техники псевдоожижения. М.: Химия, 1967.

  3. Мухленов И.П., Сажин Б.С., Фролов В.Ф. Расчеты аппаратов кипящего слоя. Л.: Химия, 1986.

  4. Протодьяконов И.О., Чесноков Ю.Г. Гидромеханические основы процессов химической технологии. Л.: Химия, 1987.

  5. Касаткин А.Г. Основные процессы и аппараты химической технологии. М.: ТИД Альянс, 2004.

  6. Разинов А.И., Клинов А.В., Дьяконов Г.С. Процессы и аппараты химической технологии. Казань: КНИТУ, 2017.

  7. Дытнерский Ю.И. Основные процессы и аппараты химической технологии. Пособие по проектированию. М.: Химия, 1991.

  8. Айнштейн В.Г., Захаров М.К., Носов Г.А., Захаренко В.В., Зиновкина Т.В., Таран А.Л., Костанин А.Е. Общий курс процессов и аппаратов химической технологии. Книга 1. М.: Химия, 1999.

  9. Trushin A.M., Ravichev L.V., Yashin V.E., Nosyreva M.A., Komlyashev R.B. Semiempirical method for determining the rate of slow hindered motion // Theor. Found. Chem. Eng. 2019. V. 53. № 6. Р. 967. [Трушин А.М., Равичев Л.В., Яшин В.Е., Носырев М.А., Комляшев Р.Б. Полуэмпирический метод определения скорости стесненного медленного движения // Теор. осн. хим. технол. 2019. Т. 53. № 6. С. 617.]

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