Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

HEAT AND EFFICIENCY MASS TRANSFER IN PACKED LAYERS

Laptev A.G. 2 Basharov M.M. 1 Lapteva E.A. 2
1 OAO «TANECO»
2 VPO «Kazan State Power Engineering University»
The article focuses on the processes of turbulent transfer of heat and mass in one – and two-phase media in a stationary Packed and granular layers with chaotic and regular structure. Considered a one-dimensional model of heat and mass transfer in a turbulent boundary layer and theoretical expressions for the average transfer coefficients. Use of the properties of conservativeness of the laws of friction in the boundary layer in the presence of perturbations that allowed to apply the one-dimensional model for calculation of heat and mass transfer in Packed layers. The calculations are made, the coefficients of heat and mass transfer are compared with experimental data of various authors. The calculated efficiency of mass transfer in the gas phase at film mode and shows the agreement with the experiment. A mathematical model was developed and theoretical expressions are used when designing and modernization of industrial units in a petrochemical complex.
heat transfer
mass transfer
nozzle
boundary layer
turbulence
1. Ajerov M.Je., Todes O.M., Narinaskij D.A. Apparaty so stacionarnym zernistym sloem: Gidravlicheskie i teplovye osnovy rascheta. Himija. L., 1979. 176 р.
2. Blinichev V.N., Komlev V.G., Zaharov V.M. i dr. Issledovanie kojefficientov soprotivlenija i teplootdachi sloja nasadki // Izv. vuzov. Himija i him. tehnol. 1987. T. 30, no. 2. рр. 124–126.
3. Djakonov S.G., Elizarov V.V., Elizarov V.I. Teoreticheskie osnovy proektirovanija promyshlennyh apparatov himicheskoj tehnologii na baze soprjazhennogo fizicheskogo i matematicheskogo modelirovanija: Feder. agentstvo po obrazovaniju, Kazan. gos. tehnol. un-t. Kazan: KGTU, 2009. 456 р.
4. Djakonov S.G. Modelirovanie massoteploperenosa v promyshlennyh apparatah na osnove issledovanija laboratornogo maketa / S.G. Djakonov, V.I. Elizarov, A.G. Laptev // TOHT. 1993. T. 27, no. 1. рр. 4–18.
5. Laptev A.G, Mineev N.G., Malkovskij P.A Proektirovanie i modernizacija apparatov razdelenija v nefte- i gazopererabotke. Kazan: 2002. 220 р.
6. Laptev A.G. Matematicheskoe modelirovanie processov masso- i teplootdachi v gazovoj faze nasadochnyh kolonn / A.G. Laptev, S.G. Djakonov // Himicheskaja promyshlennost. 1993. no. 6. рр. 48–51.
7. Laptev A.G. Modeli pogranichnogo sloja i raschet teplomassoobmennyh processov. Kazan: Izd-vo Kazansk. un-ta, 2007. 500 р.
8. Laptev A.G., Farahov T.M.,. Lapteva E.A Modeli javlenij perenosa v neuporjadochennyh nasadochnyh i zernistyh slojah // Teor. osnovy him. tehn. 2015. T.49. no. 4. рр. 407–414.
9. Laptev A.G. Modeli perenosa i jeffektivnost zhidkostnoj jekstrakcii. Kazan: Kazansk. gos. jenerg. un-t, 2005. 229 р.
10. Laptev A.G., Nikolaev N.A., Basharov M.M. Metody intensifikacii i modelirovanija teplomassoobmennyh processov. M.: Teplotehnik, 2011. 336 р.
11. Laptev A.G., Farahov T.M. Model massootdachi v zernistyh i nasadochnyh slojah // Izv. Vuzov. Himija i himicheskaja tehnologija. 2013. T. 56. рр. 92–96.
12. Laptev A.G., Lapteva E.A. Obobshhenie gidrodinamicheskoj analogii dlja razlichnyh uslovij obtekanija poverhnosti // Vestnik Kazan. tehnol. un-ta (KHTI). 2013. T. 16. no. 23. рр. 64–69.
13. Laptev A.G., Lapteva E.A. Opredelenie kojefficientov turbulentnogo peremeshivanija v odno- i dvuhfaznyh sredah po modeli Tejlora // Fundamentalnye issledovanija. 2015. no. 2 (chast 13). рр. 2810–2814.
14. Ramm V.M. Absorbcija gazov. Himija. M., 1976. 656 р.
15. Farahov M.I., Laptev A.G., Farahov T.M. Opredelenie kojefficientov massootdachi v gazovoj faze oroshaemyh nasadochnyh kolonn v protivotoke // Mezhdunarodnyj zhurnal prikladnyh i fundamentalnye issledovanija. 2014. no. 9. рр. 50–53.
16. Farahov M.I., Laptev A.G. Jenergojeffektivnoe oborudovanie razdelenija i ochistki veshhestv v himicheskoj tehnologii // Vestnik Kazans. tehnol. un-ta (KHTI). 2011. no. 9. рр. 152–158.

В статье рассматриваются процессы переноса в одно- и двухфазных средах в стационарных насадочных и зернистых слоях с хаотичной засыпкой при турбулентном режиме движения газов и жидкостей сквозь слой. Такие слои широко используются для проведения как гидромеханических процессов (например, проточные статические смесители, насыпные фильтры), так и тепломассообменных и реакционных процессов. Сюда относятся некоторые процессы гетерогенного катализа, сушка в слое, нагревание материала в шахматных и доменных печах, адсорбция, абсорбция, ректификация и т.п. [1, 5, 14, 16].

Ниже показано применение моделей пограничного слоя для получения расчетных выражений средних коэффициентов тепло- и массоотдачи и их использование для расчета эффективности тепло- и массопереноса в каналах и аппаратах с насадками при турбулентном режиме (Reэ > 50) [14].

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

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

Среднее касательное напряжение на стенке τст находят, используя уравнение баланса сил, а если это затруднительно, то выражают через среднюю объемную или массовую скорость диссипации энергии в аппарате [4, 6, 7]. Например, для неупорядоченного насадочного или зернистого слоя такое выражение имеет вид [7, 8]

laptev01.wmf (1)

где κ – коэффициент пропорциональности; ρ – плотность фазы кг/м3; laptev02.wmf – средняя скорость диссипации энергии, Вт/м3; ν – коэффициент молекулярной вязкости, м2/с. Среднюю скорость диссипации энергии laptev03.wmf можно выразить через фиктивную скорость среды w0 в полном сечении аппарата без насадки

laptev04.wmf (2)

где ξ – коэффициент гидравлического сопротивления; aν – удельная поверхность насадки, м23; εсв – удельный свободный объем насадки; w0 – средняя скорость, м/с.

Или через среднюю истинную скорость среды laptev05.wmf

laptev06.wmf (3)

Используя эквивалентный диаметр насадки laptev07.wmf из (3), можно получить

laptev08.wmf (4)

Выражение (2) следует из зависимости для средней диссипации энергии

laptev09.wmf (5)

где ΔP – перепад давления, Па; S – площадь сечения аппарата, м2.

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

В теории пограничного слоя широко используется понятие динамической скорости

laptev10.wmf.

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

laptev11.wmf (6)

где Re – число Рейнольдса laptev12.wmf.

В работах [7, 8] при турбулентном режиме движения среды (Reэ > 40) для насадок установлен коэффициент пропорциональности в интервале 1,7–2,0. Возьмем среднее значение κ = 1,85.

Модели тепло- и массоотдачи

Для определения средних коэффициентов тепло- и массоотдачи используются модели пограничного слоя в одномерной постановке с применением различных функций турбулентного обмена [6–9]. Записывается сопротивление переносу тепла и массы компонента поперек пограничного слоя в виде

laptev13.wmf

laptev14.wmf (7)

где α – коэффициент теплоотдачи, Вт/(м2∙К); β – коэффициент массоотдачи, м/с; δ – толщина пограничного слоя, м; a, D – коэффициенты молекулярного переноса тепла и массы (диффузии), м2/с; aт, Dт – коэффициенты турбулентного переноса тепла и массы, м2/с; y – поперечная координата в пограничном слое, м; ρ – плотность среды, кг/м3; cp – удельная теплоемкость среды, Дж/(кг/К).

Например, из выражения (7) с использованием функций турбулентной вязкости Г.А. Дрейцера, Deissler и др. получены формулы [10]:

laptev15.wmf (8)

laptev16.wmf (9)

где значение R′ для плоского пограничного слоя без возмущения R′ = 5,303; Rδ – безразмерная толщина пограничного слоя; Sc – число Шмидта Sc = ν/D; Pr – число Прандтля, laptev17.wmf.

Значение R′ фактически является одним из пределов интегрирования в формулах (7) и близко к значению безразмерной толщины вязкого подслоя в трехслойной модели (~5). При наличии возмущений это значение будет отличаться от пяти. Такими возмущениями являются кривизна и шероховатость поверхности, сужение и расширения потока в насадочном слое. Для коррекции параметров пограничного слоя, который отличается от плоского, в работах [6–11] используется консервативность его законов к возмущениям (подход С.С. Кутателадзе, А.И. Леонтьева и др.).

С использованием данного подхода, выполняя аналогичный вывод [8] для насадки, получим

laptev18.wmf (10)

Формула для расчета среднего значения Rδ для насадки имеет вид [8]

laptev19.wmf (11)

Записывая динамическую скорость в виде (6) из формул (8) и (9) с параметрами (10), (11) получим

laptev20.wmf (12)

laptev21.wmf (13)

где laptev22.wmf laptev23.wmf – числа Нуссельта и Шервуда.

Расчеты по выражениям (12), (13) дают удовлетворительное согласование с экспериментальными данными и с уравнениями, полученными с применением моделей Прандтля и Кармана [8].

Результаты расчетов

Для проверки полученных выражений (12)–(13) используем как известные критериальные выражения, так и экспериментальные данные различных авторов в виде графиков.

Массоотдачу в газовой фазе от насадочного слоя изучали путем возгонки твердых тел (нафталина) и испарения воды с поверхности пористых тел [1, 14]. Экспериментальные данные обобщены критериальным выражением

laptev24.wmf (14)

где по данным Гильденблата А = 0,407, m = 0,655 (Reэ = 40–104); по данным Аэрова и Умника А = 0,395, m = 0,64 (Reэ = 30–2000) и по данным Шулмена А = 0,45, m = 0,64 (при Reэ = 40–3000).

В работе [11] получена приближенная обобщающая зависимость

laptev25.wmf (15)

а в работе [12] в результате обобщения гидродинамической аналогии Чилтона – Кольборна на градиентные потоки получено

laptev26.wmf (16)

Результаты расчетов по формуле (13) удовлетворительно (±15 %) согласуются с результатами, полученными по уравнениям (14)–(16) для различных насадок при 50 < Reэ < 104.

На рис. 1 представлены результаты расчета по формуле (13) и экспериментальные данные, обобщенные в работе [14]. Кольца 25 и 50 мм. Возгонка нафталина.

pic_27.tif

Рис. 1. Корреляция данных по массоотдаче в газовой фазе насадочных колонн: 1 – расчет по уравнению (13); 2 – эксперимент [14]

На рис. 2 дано сравнение по числу Нуссельта при теплоотдаче в нерегулярном насадочном слое из цилиндров и параллелепипедов (ξ = 11,6/Re0,25) с опытными данными для воздуха [2].

pic_28.tif

Рис. 2. Теплоотдача в слоях насадочных элементов: 1 – расчет по уравнению (12); 2 – опытные данные [2]

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

Двухфазные потоки в насадках

Рассмотрим применение приведенных выражений для расчета коэффициентов массоотдачи в орошаемых насадках при противоточном пленочном режиме с газовым потоком. До начала подвисания жидкости есть рекомендации производить расчет Shг по уравнениям для сухих насадок. Действительно, поправка Шульмана на орошение при пленочном режиме близка к единице [15].

В работе [3] получены опытные данные по гидравлическому сопротивлению и массоотдаче в газовой фазе для регулярной металлической рулонной насадки. Насадка образована сдвоением металлическими лентами, одна из которых имеет гофры треугольной формы. На боковых поверхностях гофр выполнены лепестки в виде круговых сегментов, при этом хорды сегментов гофр расположены под углом друг к другу. Свободный объем насадки εсв = 0,95, удельная поверхность av = 480 м23, эквивалентный диаметр dэ = 0,0079 м. Опыты проводились при скорости воздуха 0,5 до 3,0 м/с и плотности орошения от 5 до 30 м3/(м2∙ч).

В результате обобщения опытных данных получены коэффициенты сопротивления сухой ξо и орошаемой насадки ξор в виде

laptev27.wmf (17)

laptev28.wmf (18)

где laptev29.wmf – число Рейнольдса жидкой фазы; laptev30.wmf – плотность орошения, м3/(м2∙с).

Коэффициент массоотдачи в газовой фазе находился экспериментально при увлажнении воздуха водой (S = 0,7). На рис. 3 представлены результаты расчетов коэффициента массоотдачи laptev31.wmf по формуле (13) и опытные данные [3].

Согласование с расхождением в пределах погрешности эксперимента

Если при расчете насадочных колонн используется объемный коэффициент массоотдачи, то запишем

laptev32.wmf (19)

где ψw – коэффициент смачиваемости поверхности насадки.

Структуру потока газовой фазы можно оценить по числу Пекле по перемешиванию [13]

laptev33.wmf (20)

Для рассмотренной насадки при wг = 0,5 м/с; Reэ = 263,3; ξ = 0,19 и H = 1,0 м по формуле (20) получаем Рег = 401,6, что указывает на практически идеальное вытеснение потока.

pic_29.tif

Рис. 3. Зависимость коэффициентов массоотдачи в газовой фазе от скорости газа при плотности орошения 10 (м32∙ч), dэ = 0,0079, м. 1.wmf – эксперимент [3], линия – расчет по (13)

Тогда эффективность массоотдачи в газовой фазе запишем по известному выражению

laptev34.wmf (21)

где N – число единиц переноса; F = avψvSкH – поверхность контакта фаз, м2; Vг = wгSк – объемный расход, м2/с; Sк – площадь поперечного сечения колонны, м2.

В результате имеем

laptev35.wmf (22)

Для рассмотренного выше примера при wг = 0,5 м/с получаем η = 0,999, т.е. данная насадка обеспечивает высокую эффективность разделения исследованной смеси.

pic_30.tif

Рис. 4. Зависимость эффективности массоотдачи насадки из колец Рашига от высоты слоя: 1 – Reэ = 100; 2 – Reэ = 500 (Scг = 1)

На рис. 4. даны результаты расчетов эффективности массоотдачи в насадочном слое из колец Рашига размером 50×50 мм (av = 110 м23; εсв = 0,95).

Заключение

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

Статья выполнена в рамках проектной части государственного задания в сфере научной деятельности. Задание № 13.405.2014/К.

Рецензенты:

Зиятдинов Н.Н., д.т.н., профессор, заведующий кафедрой «Системотехника», ФГБОУ ВПО «Казанский национальный исследовательский технологический университет», г. Казань;

Ваньков Ю.В., д.т.н., профессор, заведующий кафедрой «Промышленная теплоэнергетика и системы теплоснабжения», ФГБОУ ВПО «Казанский государственный энергетический университет», г. Казань.