УДК 697.922.2/697.942.2
Гидродинамический расчет орошаемой колонны с пористыми насадочными телами
Д-р техн. наук А. В. ЦЫГАНКОВ1, д-р техн. наук В. А. ПРОНИН2, Д. И. ШПИЛИН3, А. Е. АЛЕШИН
1ра1Ма-Ш@шЮрго^рЬ^и, 2таюг.Ршшп@таД.:га, 3shpilinspb@gmail.com
Университет ИТМО Институт холода и биотехнологий 191002, Санкт-Петербург, ул. Ломоносова, 9
В данной статье представлена методика расчета гидродинамических характеристик вертикальной колонны с насадочными телами различной формы. Такие колонны используются для очистки и дезодорации атмосферного воздуха, содержащего в себе примеси вредные для человека и окружающей среды. В основу методики заложен принцип математического моделирования, а так же приводятся основные формулы, для оценки гидродинамических процессов, с использованием критерия Рейнольдса, закона Дарси, а так же формулы Блазиуса и Никурадзе. Авторами предложена методика гидродинамического расчета орошаемой колонны, позволяющий оценить процессы, протекающие в колоннах, такие как: расход воздуха, пленочный и фильтрационный расход жидкости, средняя скорость воздуха и жидкости, а так же толщина пленки жидкости. Предложенный метод позволяет качественно оценить гидродинамические процессы, происходящие в насадочной колонне орошаемого типа. Дальнейшие исследования позволят получить количественные характеристики и разработать инженерную методику расчета подобных систем. Ключевые слова: орошаемые колонны, пористые насадочные тела, дезодорация, очистка воздуха, гидродинамический расчет.
Calculation of hydrodynamics for spray tower with porous packing
D. Sc. A. V. TSYGANKOV1, D. Sc. V. A. PRONIN2, D. I. SHPILIN3, A. E. ALESHIN
1pallada-ltd@mfopra.spb.su, 2maior.Pronin@mail.ru, 3shpilinspb@gmail.com University ITMO Institute of Refrigeration and Biotechnologies 191002, Russia, St. Petersburg, Lomonosov str, 9
The article deals with calculation methods of hydrodynamics for spray tower with porous packing. The towers of this type are used for cleaning and deodoriza-tion of contaminated outdoor air. Method is based on the mathematical modeling principle. Main formulae for hydrodynamics estimation are given. Reynolds criteria, Darcy's law, Blasius and Nikuradze formulae are used. The calculation methods of hydrodynamics for spray tower have been proposed. It allows evaluating tower processes e. g. air consumption, film and filtration liquid consumption, air and liquid average velocity and the thickness of liquid film. Further research will allow to obtain quality characteristics and to develop engineering calculation methods for the systems of this type.
Keywords: spray tower, porous packing, deodorization, air cleaning, calculation of hydrodynamics.
Очистка и дезодорация больших объемов атмосферного воздуха является актуальной задачей для систем жизнеобеспечения современных мегаполисов. Формирование городских агломерации и одновременно дефицит
земель пригодных для жилищного строительства приводит к тому, что жилые районы приближаются к очистным сооружениям систем канализации, водоотведения и полигонам твердых бытовых отходов. Для дезодорации больших объемов воздуха вблизи этих объектов широко используются орошаемые насадочные колонны. Одним из перспективных направлений повышения эффективности очистки является изготовление насадочных тел (кольца Рашига, кольца Палля, седла Берля и пр.) из пористых материалов [1].
Рассмотрим метод расчета вертикальной колонны с насадкой такого типа (рис) [2]. Исходные геометрические данные для расчета: 5", Ж — площадь поперечного сечения и объем засыпки, соответственно; М — масса
I
засыпки телами /-го типоразмера; 5., м т. — площадь поверхности, объем и масса отдельного насадочного тела /-го типоразмера, соответственно.
Вычислим удельную поверхность контакта засыпки — а и долю свободного объема засыпки — е:
е
n Mts,
a = -
m,.
W
е
e = -
n M Wi
i=0 m,.
W
где п — количество типоразмеров насадочных тел.
В качестве расчетной модели засыпки примем изотропное пористое тело со сквозными вертикальными плоскими каналами (рисунок). Ширину канала — И и расстояние между каналами — t определим из условия равенства
i=0
ХОЛОД: ТЕХНИКА И ТЕХНОЛОГИИ
35
Пористый материал
Уф =-k grad(-+ у),
рических характеристик пор
к =
Рж
, тогда
=
m я
L
- Уж
Полагая, что в пористой среде и каналах давление линейно изменяется по высоте засыпки, градиент давления вычислим по формуле:
5ф = (1 - е)
Тотда Ga = ^ Уф
(3)
Рассмотрим ламинарное стационарное течение жидкой пленки при противотоке газа. Совместим поверхность канала, по которой движется пленка с координатной плоскостью у — г, а ось у с направлением движения газа. Запишем уравнение движения для рассматриваемого случая:
Схема гидродинамического расчета орошаемой колонны
а и е в реальной модельной засыпки. Так как количество каналов за единицу длины равно Ш, то t = 2/а.
Жидкий реагент подается в верхнюю часть колонны и стекает вниз по стенкам каналов и фильтруется через пористое тело:
Об = О + Оф, (1)
общ пл ф^ у '
где Оо6щ — общий объемный расход жидкости; Опл — пленочный расход; Оф — фильтрационный расход.
Движение вязкой жидкости через пористую среду описывается законом Дарси [3]
Иж ~t + 1--Рж g = 0
dx dy
(4)
где V — скорость движения жидкости вдоль оси у.
Интегрируя уравнение (4) с граничными условиями vж = 0 при х = 0 (условие прилипания) и vж = vгр при х = 5 получим
Цж v +
dP dy '
Vе
^х2
1
/ V /
dx ' 2
¥=о,
о
% - ** К*2 "&).
(2)
где Уф — скорость фильтрации; k — коэффициент фильтрации; р — плотность жидкости; g — ускорение свободного падения; Р — давление жидкости; y — невели-ровочная высота относительно нижней границы пористой среды.
Примем, что в любом горизонтальном сечении насадки давление в пористом теле и каналах величина постоянная, тогда просачивание жидкости из пористой среды в канал не происходит, и граничные условия для уравнения (2) можно сформулировать в виде P = Рп при y = 0; P = 0 при y = L,
здесь Р — давление подачи воздуха в колонну, L = W/s— высота засыпки.
После интегрирования получим
Рж gL
Коэффициент фильтрации k выразим через коэффициент пористости с, который зависит только от геомет-
ж 5 2ЦжЧ>
Здесь у — скорость жидкости на границе раздела пленки и воздуха, 5 — толщина жидкой пленки на поверхности канала.
Определим среднюю по сечению пленки скорость жидкости у, и пленочный расход жидкости, рассчитанный по этой скорости
v = — Т v dx =
ж 5 J о ж
12ц
dP
dy '
G 255 _
G =-_
пл t ж
(5)
(6)
где д — динамический коэффициент вязкости; у — удельный вес жидкости.
Для того, чтобы связать толщину пленки и скорость на границе раздела сред, рассмотрим взаимодействие воздушного потока жидкой с пленкой.
Ширину канала, свободного для прохода газовой фазы обозначим 4:
4 = t - 25.
Переместим начало координатной системы (х, у, г) на ось симметрии канала. Режим движения воздуха в канале определяется по критерию Рейнольдса:
= Кекр,
т р
где V — скорость воздуха вдоль оси у; д — динамический коэффициент вязкости воздуха; р — плотность воздуха; d = 24 — гидравлический диаметр плоского канала; Reкр — критическое число Рейнольдса.
Пренебрегая действием объемных сил, для стационарного ламинарного движения воздуха в канале можно записать в виде:
= —. (7)
dx2
dy
Граничные условия для уравнения (7) зададим в виде:
m
ж
c
V
У
v = Угр при x = ± ^/2.
После интегрирования уравнения движения и вычисления постоянных коэффициентов интегрирования получим
х2 &
v = v —
гр 8m dy
1 — 4
/ л2 x
(8)
При турбулентном режиме движения воздуха вычислим осредненную по ширине канала скорость воздуха
v =1 fX vdx
X J-X'2
=Vrp - V0>
(9)
здесь vQ — средняя скорость воздуха относительно поверхности жидкой пленки
_ X2 ОР
у° = ^ Т-12ц ах
Выразим градиент давления с помощью коэффициента гидравлического трения X:
ар
ах 2Х'
В общем случае коэффициент гидравлического трения выражается через число Рейнольдса степенной
к
функцией вида 1 = К + • Коэффициенты К12 и т
зависят от режима движения жидкости и формы канала [1]. В частности, для рассмотренного выше случая
96
ламинарного режима — 1 = —. При турбулентном
Re
режиме могут быть использованы формулы Блазиуса, Никурадзе [4] и др.
Выразим скорость V через X с учетом уравнения (3)
V = V - Щ * \£1 р
Общий расход воздуха через насадку вычисляется по формуле
G =Х vS = 3l - — I vS.
(10)
С другой стороны, расход определяется разницей напора и расходной характеристикой вентилятора, которая, в общем случае, может быть задана в виде полинома:
=е :=0 . (и)
Таким образом, для определения искомых величин: расхода воздуха, пленочного и фильтрационного расхода жидкости, средних скоростей воздуха и жидкости, а так же толщины пленки, необходимо решить систему уравнений (2), (3), (5), (6), (9)-(11).
В заключение необходимо отметить, что предложенный метод позволяет получить только качественные оценки гидродинамических процессов, проходящих
в насадочной колонне. Количественные оценки могут быть получены при использовании имитационной модели построенной на основе методов вычислительной гидродинамики, при верификации которой предложенный метод расчета может оказаться весьма полезным.
Список литературы
1. Кутателадзе С. С. Теплопередача и гидродинамическое сопротивление. — Энергопромиздат, 1990. 367 с.
2. Пронин В. А., Молодое М. А., Щпилин Д. И. Газовоздушные выбросы пищевых предприятий и способы их устранения. // Научный журнал НИУ ИТМО. Серия «Процессы и аппараты пищевых производств». 2013. № 2.
3. Леонтьев Н. Е. Основы теории фильтрации. — М.: Изд-во ЦПИ при механико-математическом факультете МГУ, 2009. 88 с.
4. Гейер В. Г., Дулин В. С., Заря А. Н. Гидравлика и гидропривод: Учеб для вузов. — 3-е изд., перераб. и доп. — М.: Недра, 1991.
5. Цыганков А. В. Моделирование движения гидродинамической смазки в узлах трения приборов с учетом шероховатости рабочих поверхностей // Известия высших учебных заведений. Приборостроение. 2004. Т. 47. № 2. С. 43-48.
6. Зуланке Р., Виттен П. Дифференциальная геометрия и расслоения. — М.: Мир, 1975.
7. Хоеалыг Д. М., Синицина М., Бараненко А. В., Цой А. П. Энергоэффективность и экологическая безопасность техники низких температур // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование» 2014. № 1.
8. Полянин А. Д., Зайцев В. Ф. Нелинейные уравнения математической физики. Точные решения. Справочник. — М.: ФМЛ, 2002. 432 с.
9. Федяевский К. К., Войткунский Я. И., Фаддеев Ю. И. Гидромеханика. — Л.: Судостроение, 1968. 568 с.
References
1. Kutateladze S. S. Teploperedacha i gidrodinamicheskoe soprotivlenie. — Jenergopromizdat, 1990. 367p.
2. Pronin V. A., Molodov M. A., Shhpilin D. I. Nauchnyj zhurnal NIUITMO. Serija «Processy i apparaty pishhevyh proiz-vodstv». 2013. No 2.
3. Leont'ev N. E. Osnovy teorii fil'tracii. — M.: Izd-vo CPI pri mehaniko-matematicheskom fakul'tete MGU, 2009. 88 p.
4. Gejer V. G., Dulin V. S., Zarja A. N. Gidravlika i gidrop-rivod: Ucheb dlja vuzov. — 3-e izd., pererab. i dop. — M.: Nedra, 1991.
5. Tsygankov A. V. Izvestija vysshih uchebnyh zavedenij. Priborostroenie. 2004. T. 47. № 2. S. 43-48.
6. Zulanke R., Vitten P. Differencial'naja geometrija i rass-loenija. — M.: Mir, 1975.
7. Khovalyg D. M., Sinityna M., Baranenko A. V., Tsoy A. P. Nauchnyj zhurnal NIU ITMO. Serija «Holodil'naja teh-nika i kondicionirovanie» 2014. No 1.
8. Poljanin A. D., Zajcev V. F. Nelinejnye uravnenija matematicheskoj fiziki. Tochnye reshenija. Spravochnik. — M.: FML, 2002. 432 p.
9. Fedjaevskij K. K., Vojtkunskij Ja. I., Faddeev Ju. I. Gidromehanika. — L.: Sudostroenie, 1968. 568 p.