ПРИКЛАДНАЯ ХИМИЯ И ХИМИЧЕСКАЯ ТЕХНОЛОГИЯ
УДК 544.478
В. Б. Репин, А. С. Балыбердин, Ф. Ш. Шарафисламов,
А. Ф. Махоткин
МАКСИМАЛЬНАЯ РАБОЧАЯ ПОВЕРХНОСТЬ СОТОВОГО КАТАЛИЗАТОРА
Ключевые слова: сотовый катализатор, максимальная рабочая поверхность.
В настоящей работе для сотовых катализаторов определены геометрические параметры, обеспечивающие максимальную рабочую поверхность контакта фаз. Рассмотрены сотовые катализаторы, которые содержат продольные каналы в форме треугольника, квадрата и шестиугольника. Показано, что максимальная поверхность достигается при выполнении равенства толщины перегородки величине эквивалентного диаметра продольных каналов блока сотового катализатора.
Keywords: cell catalyst, the maximum working surface.
In this paper, for cell catalysts determined the geometric parameters for maximum work surface contact between the phases. Considered honeycomb catalysts which comprise longitudinal channels to form a triangle, square and hexagon. It is shown that the maximum surface equality is achieved when the thickness partition size equivalent diameter of the longitudinal channels unit cell of the catalyst.
Широкое применение сотовые катализаторы нашли для очистки газовых выбросов теплоэлектростанций. Например, в Японии полностью перешли на сотовые катализаторы, которые изготавливают на основе диоксида титана. Стоимость технологии газоочистки на зарубежных теплоэлектростанциях достигает 30% от стоимости самой теплоэлектростанции. Сокращение капитальных и эксплуатационных затрат на газоочистку от окислов азота, достигаются путем изготовления дешевых катализаторов, полученных из железооксидных и железохромооксидных отходов различных производств [1].
Испытания гранулированного железохромоксидного катализатора на ТЭЦ-11 (г. Москва) показали [2], что при концентрации окислов азота на входе до 1,0 г/м3 концентрация на выходе составила 0,03 г/м3, что в 4 раза меньше допустимой нормы выброса.
Остаточная концентрация аммиака в газовой фазе составила 0,02 г/м3. За шесть месяцев непрерывной работы опытно-промышленной установки отложений продуктов неполного сгорания топлива и сажи на катализаторе не наблюдалось. Железохромоксид-ный гранулированный катализатор не только эффективно восстанавливает оксид азота аммиаком, но эффективно окисляет до нейтральных соединений продукты неполного сгорания топлива [3, 4]. За весь период непрерывных испытаний снижения эффективности катализатора по очистке газовых выбросов от окислов азота не обнаружено. Были изготовлены образцы сотовых катализаторов, которые прошли испытания по очистке газовых выбросов теплоэлектростанций от окислов азота в опытных установках на ТЭЦ-11 (г. Москва), ТЭЦ-1 (г. Набережные Челны) и Заинской ГРЭС. Ресурс надежной работы сотовых катализаторов больше, чем для гранулированных катализаторов.
Анализ современного состояния технологии производства сотовых катализаторов показывает, что для обеспечения наибольшей эффективности работы сотового, как и любого иного, катализатора
необходимо обеспечить максимальную рабочую поверхность контакта фаз.
Ниже приводятся результаты максимилиза-ции рабочей поверхности блока сотового катализатора, в котором продольные каналы имеют в поперечном сечении треугольник, квадрат, либо правильный шестиугольник [5, 6]. Полагается, что блок представляет собой цилиндр наружным диаметром Б и длиной Ь.
В качестве характерного поперечного размера внутреннего канала сотового катализатора принимается гидравлический диаметр канала (ф, равный учетверенной площади поперечного сечения канала, деленной на его периметр, 5 - толщина перегородки между каналами.
Полная рабочая поверхность блока 8кЕ вычисляется как сумма внутренней поверхности, образованной продольными каналами, 8к, боковой поверхности блока Б и удвоенной поверхности торцов блока Тк, образованной поперечными сечениями перегородок между продольными каналами. На самом деле рабочая поверхность катализатора определяется брутто-реакцией процесса. Например, при производстве серной кислоты на катализаторах осуществляется окисление двуокиси серы в трех-окись серы. Эта реакция экзотермическая, что обуславливает повышение температуры процесса и выхода этого параметра из оптимального по температуре. В связи с этим процесс каталитического окисления двуокиси серы в целевой продукт (803) проводят в многоступенчатом режиме и на каждой ступени поддерживают оптимальную температуру процесса путем добавления к продуктам реакции кислорода воздуха. Следовательно, в каждом индивидуальном слое катализатора целевая реакция не завершается. Тогда рабочая поверхность сотового катализатора будет совпадать с геометрической (смоченной) поверхностью сотового блока. При непосредственном контакте элементов сотового катализатора друг с другом торцевые поверхности пере-
крываются. В дальнейшем эта ситуация не рассматривается.
Поэтому
SkZ=Sk+Б+2Tk (1)
Площадь боковой поверхности блока не зависит от геометрии продольного канала блока и вычисляется как
Б=^
Величина площади торца блока зависит от конфигурации каналов и вычисляется как
Tk= (пО2/4)^^ (2)
- площадь живого сечения канала, N -количество каналов в блоке.
Внутренняя поверхность блока, образованная продольными каналами, вычисляется из следующего соотношения
Sk=ПкNkL (3)
Обозначим конфигурацию формы каналов следующим образом. Нижний индекс в уравнениях (2) и (3) принимает значения к = 3, 4, 6 и идентифицирует конфигурацию канала в сотовом катализаторе. Если к = 3, то это соответствует каналу в форме правильного треугольника, к = 4 - квадрата, к = 6 -правильного шестиугольника.
s3=(3*3(1/2)*d2)/4, s4=d2, s6=(3(1/2)*d2)/2 (4)
Величина Пк - периметр поперечного сечения канала, который вычисляется как
П3=3*3(1/2Ы П4=4d, П6=2*3(1/2^ (5)
Общее количество каналов (Кк), заполняющих поперечное сечение блока вычисляется, исходя из следующих рассуждений.
Рассмотрим квадрат покрытый, например треугольными элементами. Для поперечного сечения блока, например с треугольными каналами, можно выделить две взаимно перпендикулярные оси - горизонтальную и вертикальную. Для горизонтальной оси период ячейки треугольного элемента равен (30,5)^ + 5) и общее количество ячеек, укладывающихся на этой оси, равно ^) = (2/30,5)и. Для вертикальной оси период треугольной ячейки равен 1,5^ + 5) и общее количество ячеек, укладывающихся на вертикальной оси, равно ^2) = (2/3)и. Тогда общее количество ячеек, заполняющих с зазором 5 квадрат размером Б, равно ^вадрат = №^(2). А количество ячеек заполняющих круг диаметром Б равно N = (л/4)Мквадрат.
Окончательно получим соотношения для вычисления числа продольных каналов для блока круглого сечения в виде
^=(п*3(1/2)/3)*и2, ^=(п/4)*и2,
^=(0,5п/3(1/2))*и2 (6)
и=(О-б)/@+б)
Воспользовавшись соотношениями (4) - (6), можно показать, что как внутренняя поверхность сотового катализатора, так и поверхность торцов не зависят от формы поперечного сечения продольных каналов и равна
S3= S4= S6=пU2Ld (7)
Т3= Т4= Тб=(п/4)*(О2^2и2) (8)
Поэтому полная рабочая поверхность блока сотового катализатора также не зависит от формы поперечного сечения продольных каналов и вычисляется из следующего соотношения
5к^=п^0[(1-0,5(й^))*(й/0)*и2+1+0,5(0^)] (9)
1.2
О 4 8 12
Гидравлический диаметр канала, мм
Рис. 1 - Зависимость величины рабочей поверхности сотового катализатора от размера квадратного канала. Толщина перегородки 5 = 1 мм, диаметр блока D = 87,5 мм, длина блока L = 177 мм
Анализ уравнения (9) показывает, что при неизменной толщине перегородки (5 - Const), уменьшение размера окна (поперечного сечения каналов сотового блока) приводит к увеличению поверхности сотового катализатора. Эта тенденция сохраняется при снижении размера окна вплоть до величины d (рис. 1). Дальнейшее уменьшение размера окна (d < d) приводит к снижению рабочей поверхности сотового катализатора.
Максимальная рабочая поверхность блока сотового катализатора реализуется при выполнении условия равенства нулю производной функции по формуле (9). Тогда максимальная рабочая поверхность сотового блока катализатора реализуется при выполнении условия
d=6/(1-(6/L)) (10)
Соотношение (10) показывает, что условие максимальности рабочей поверхности не зависит от
диаметра блока сотового катализатора, а, следовательно, и от площади торцевых поверхностей блока. Для ситуации, когда толщина перегородки существенно меньше длины блока, соотношение (10) сводится к уравнению
в трехмерном пространстве с последующим вычислением поверхности и массы изделия (рис. 3).
В табл. 1 и табл. 2 приведены результаты расчетов по предлагаемым формулам и по программе КОМПАС - 3Б VII для различных ситуаций.
d=б
(11)
т. е максимальная рабочая поверхность реализуется при выполнении условия равенства толщины перегородки гидравлическому диаметру канала и вычисляется по формуле
8тах1=п1-0[((0/1-)-1)2*(С/40)+1+0,5(0/1_)] (12)
Если диаметр блока сопоставим с его длиной и обе эти величины существенно больше поперечного размера канала, то соотношение (12) принимает простой вид
Бтах^^О2/^))
(13)
Масса блока сотового катализатора рассчитывается по формуле М = Vp, где р - удельная плотность материала из которого изготовлен сотовый катализатор, V - объем занимаемый катализаторной массой
У=Т^
Тогда, масса блока равна
М=(п/4^(О2^2и2)
(14)
(15)
А для крупных блоков при выполнении условия максимальности рабочей поверхности получим
М=(3/16)^О2
(16)
Следовательно, для крупных блоков сотового катализатора его масса не зависит от геометрической формы продольных каналов и их поперечных размеров. Масса сотового блока определяется его наружным диаметром и длиной.
Воспользовавшись соотношениями (13) и (14) вычислим удельную поверхность блока при выполнении условия максимальности рабочей поверхности (11)
^тех1/ У)=(4^), м2/м3
(17)
Графическая зависимость максимальной площади поверхности катализатора от гидравлического диаметра продольного канала, представлена на рис. 2.
Приведенные вычисления не учитывали тот факт, что цилиндрическая форма блока катализатора искажает истинную конфигурацию элементарных ячеек, расположенных на периферии. Анализ влияния этого фактора проводился с использованием программы моделирования КОМПАС - 3Б У11, позволяющей конструировать катализаторный блок
8000-1
АППЛ —
оиии
111(1(1 —
•♦иии
1ППЛ —
іЧКи
0 - 1 1 1
0 4 8
Гидравлический диаметр канала, мм
12
Рис. 2 - Зависимость максимальной рабочей поверхности сотового катализатора от гидравлического диаметра продольного канала
Расчет по программе КОМПАС - 3Б V!! показал, что анализируемые величины для каналов различной конфигурации изменяются в пределах 8е = 2,03712 + 0,0222 м2, М = 2,42936 + 0,00044 кг (табл. 1) в случае крупных блоков, а в случае мелких блоков 8е = 0,04180 + 0,00037 м2, М = 0,04148 + 0,0003 кг (табл. 2). Такие незначительные отличия, в пределах 1%, указывают на то, что форма поперечного сечения канала не влияет на основные показатели блока сотового катализатора.
Сравнение результатов расчета по точному алгоритму с результатами расчета по предлагаемым формулам показало, что теоретические зависимости систематически занижают как рабочую поверхность блока сотового катализатора, так и его массу. Но эти различия не существенны.
Рис. 3 - Поперечное сечение сотового катализатора с каналами в форме правильных треугольников
Таблица 1 - Расчет массы и рабочей поверхности блока сотового катализатора (Б = 0,12 м,Ь = 0,177 м, а = 5 = 0,001 м, р = 1600 кг/м3)
Форма канала Площадь рабочей поверхности блока (м2)
3Б ¥11 Ур-е (9) Ур-е (13)
Треугольник 2,03842 2,023752 2,001764
Квадрат 2,03490
Шестиугольник 2,03803
Форма канала Масса блока (кг)
3Б ¥11 Ур-е (15) Ур-е (16)
Треугольник 2,42899 2,426867 2,402117
Квадрат 2,42980
Шестиугольник 2,42928
Таблица 2 - Расчет массы и рабочей поверхности блока сотового катализатора (Б =0,0328 м, Ь =0, 05 м, а =0,002 м, 5 =0,001 м, р =1600 кг/м3)
Форма канала Площадь рабочей поверхности блока (м2)
3Б ¥11 Ур-е (9) Ур-е (13)
Треугольник 0,042171 0,038608 0,0211241
Квадрат 0,041682
Шестиугольник 0,041548
Форма канала Масса блока (кг)
3Б ¥11 Ур-е (15) Ур-е (16)
Треугольник 0,041780 0,041664 0,050698
Квадрат 0,041117
Шестиугольник 0,041548
Так, для крупных блоков расчет площади рабочей поверхности по формуле (9) отличается на
0,7% от точного расчета, а по формуле (13), которая справедлива при выполнении условия максимальности рабочей поверхности сотового катализатора - не более чем на 1,83%. Аналогичные отклонения отмечаются и при вычислении массы блока (0,12% и 1,15%, соответственно).
Наибольшие расхождения наблюдаются для мелких блоков. При диаметре блока 32,8 мм и длине 50 мм различие между двумя алгоритмами расчета увеличивается (табл. 2). При вычислении рабочей
поверхности сотового блока по уравнению (9) отличие составляет 7,7%, а по уравнению (13), которое справедливо при выполнении условия максимальности рабочей поверхности блока - 49,5%. Для массы блока аналогичные расхождения при использовании уравнений (15) и (16) по сравнению с точным расчетом завышают этот показатель на 0,44% и 22,2%, соответственно. Следовательно, при уменьшении внешних размеров блока вклад периферийных ячеек, форма которых отличается от заданной, возрастает.
Литература
1. Махоткин, А.Ф. Теоретические основы очистки газовых выбросов производства нитратов целлюлозы / А.Ф. Махоткин. - Казань: изд. КГУ, 2003. - 268 с.
2. Махоткин, И.А. Разработка и анализ результатов внедрения в производство новых катализаторов для очистки газовых выбросов от оксидов азота на химических предприятиях и теплоэлектростанциях / И.А. Махоткин, С.П. Дмитриев, В.В. Жиляков, В.Б. Репин // Вестник Белгородского технологического университета. - 2004. -№8. - С.43-45.
3. Пат.2211728 РФ, МПК 7 В01.123/86,В01Б52/62. Катализатор для очистки газообразных выбросов от оксида углерода / А.Н. Бусыгин, В.В. Красильников, И.А. Махоткин [и др.]; заявл.12.11.2001; опубл.10.09.2003. - 4с.
4. Махоткин, И.А. Разработка технологии изготовления высокоэффективного сотового катализатора из порошкообразных железооксидных отходов для очистки газовых выбросов от оксидов азота / И.А. Махоткин, И.М. Тухватуллин, С.П. Дмитриев [и др.] // Современные проблемы технической химии: Материалы докладов Международной научно-технической и методической конференции. - Казань: КГТУ, 2004. - С.754-762.
5. Репин В. Б. Площадь поверхности сотового катализатора с продольными каналами в форме правильного треугольника. /В. Б. Репин, А. С. Балыбердин, Ф. Ш. Шарафисламов [и др.]//Вестник Казан. технол. ун-та. -2011. - Т. 14, № 23 - С 90-93
6. Репин В. Б. Оптимизация площади поверхности сотового катализатора. /В. Б. Репин, А. С. Балыбердин, Ф. Ш. Шарафисламов [и др.]//Вестник Казан. технол. ун-та.
- 2011. - Т. 14, № 18 - С 249-253
© В. Б. Репин - к.ф.-м.н, доц. каф. физики КНИТУ; А. С. Балыбердин - к.т.н., доцент, каф. ОХЗ КНИТУ, [email protected]; Ф. Ш. Шарафисламов - зав. лаб. каф. ОХЗ КНИТУ; А. Ф. Махоткин - д.т.н., проф., зав. каф. ОХЗ КНИТУ.