УДК 66.796.2
МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ОЧИСТКИ ОТРАБОТАВШИХ ГАЗОВ ХИМИЧЕСКИХ ПРОИЗВОДСТВ И ДИЗЕЛЬНЫХ АГРЕГАТОВ ОТ ТВЕРДЫХ ЧАСТИЦ СВС-ФИЛЬТРАМИ
В.В. Евстигнеев, А.Л. Новоселов, В.И. Пролубников, Н.П. Тубалов
Алтайский государственный технический университет им. И.И. Ползунова. г. Барнаул E-mail: kef@agtu.secna.ru
Разработана математическая модель фильтрации и выгорания сажи, отработавших газов дизелей через пористую поверхность СВС-материалов, полученных на основе оксидов железа и алюминия.
Вводная часть
Очистка отработавших газов в том числе твердых частиц (ТЧ) химических производств и дизельных агрегатов является одной из важнейших экологических проблем и сложной технической задачей решение которой зависит от различных факторов. В числе которых каталитическая очистка различных газов, физические свойства пористого материала. Нами предложена модель процессов фильтрации с использованием материалов, полученных на основе самораспространяющегося высокотемпературного синтеза (СВС), которая позволяет производить фильтрующие блоки для дизелей в соответствии со стандартом ЕВРО-4.
1. Состояние отработавших газов на входе в пористый проницаемый фильтр для их очистки от твердых частиц
Объемный расход отработавших газов через пористый проницаемый фильтр твердых частиц связан с параметрами источника в химическом производстве или дизеля и зависит от следующих из них: коэффициента избытка воздуха а; теоретически необходимого количества воздуха для сгорания 1 кг , кг воздуха
топлива /0, ----------- —; часового расхода топлива
кг топлива
дизелем От, кг/ч; давления наддувочного воздуха Рк, МПа (для дизелей без наддува РК=Р0); температуры воздуха после компрессора Тк, К (для дизелей без наддува ТК=Т0); газовой постоянной для возду-
К К 8,314 „ ,,
ха Кв =^~ = 28 96; Рв - действительного коэффициента молекулярного изменения. Тогда для четырехтактных двигателей объемный расход отработавших газов определяется по выражению, приведенному в работах В.И. Смайлиса, например [1]: а • 10 • От • Тк - р.
V„, =-
Pk -106
(1)
или = 4,148-10 “\ОТТк / Рк )(а+ 0,0675). (2)
Учитывая, что удельная масса отработавших газов дизелей при 273 К составляет /„=1,3 кг/м3, для другой любой температуры она определяется по известному выражению:
(/ог • Т)
Динамическая вязкость отработавших газов дизеля при Т0=273 К составляет рог=1,5-106 кгс/м2. Для вычисления динамической вязкости отработавших газов при любой температуре Тог удобно использовать эмпирическую формулу Сутерлянда:
Мог = Мо
(1 + )/Т(
(1+сг)/Т
где Сг - постоянная отработавших газов, для дизелей Сг=173.
Величину кинематической вязкости газов можно определить по выражению:
V =ЬГК
ог у ’
I ог
где #=9,81 - земное ускорение свободного падения, м/с2. Средняя скорость отработавших газов при фильтровании связана с площадью ¥ф поверхности фильтра и пористостью (Пфм):
и =
F. • n.
Режим течения определяется числом Рейнольдса. Зная среднюю скорость газов, их плотность р0„ газовую постоянную Вох, динамическую вязкость газов, можно определить число Рейнольдса из выражения:
Re =
и р
ГО*
Рог
Известно, что в капиллярных порах металлокерамических фильтров число Рейнольдса может иметь величины в пределах Re«0,001...0,1 или Re=-•2Wфм•-/3«ог(1-Ифм).
¡t Я 5
. вн н- - _Ен 1 Рн . 1 Fh
Утог = '
Т
Рис. 1. Схема трассы выпуска дизеля с фильтром твердых частиц: 1) выпускной канал; 2) выпускной коллектор; 3) соединительный трубопровод; 4) завихритель потока газов; 5) фильтр твердых частиц; 6) выпускной патрубок; 7) турбокомпрессор. Здесь участок ЕН - первая ступень очистки; Fн - вторая ступень очистки от твердых частиц
ог
Систему выпуска дизеля с сажевым фильтром или вернее фильтром твердых частиц можно представить в виде схемы, приведенной на рис. 1.
Система впуска нагревательной печи в химическом производстве масел представлена на рис. 2.
Рис. 2. Схема процесса выпуска нагревательной печи в химическом производстве масел: 1) бак с топливом; 2) воздуходувка; 3) насос; 4) форсуночное устройство; 5) камера сгорания; 6) перегородка; 7) печь; 8) трубы; 9) экран; 10) фильтр
Технологическая схема производства резинотехнических изделий представлена на рис. 2, 3, где отмечены этапы вредных выбросов в атмосферу.
Рис. 3. Технологическая схема производства резинотехнических изделий
2. Моделирование процесса фильтрования отработавших газов в пористых проницаемых СВС-блоках
Пористые проницаемые фильтрующие СВС-материалы имеют ряд особенностей, отличающий их от других материалов. Средний диаметр пор в них определяется составом шихты и технологией
изготовления, существует некоторая упорядоченность пор, поперечных стенке фильтров, направления пор и различная извилистость последних [2].
Течение отработавших газов в стенке пористого проницаемого фильтра имеет очень сложный характер, так как поток постоянно изменяет направление, огибая соединительные узлы в интерметаллическом каркасе.
Несмотря на существование многослойности в структуре, для СВС-фильтров невозможно использовать модель, предложенную Н.А. Фуксом [3]. А движение отработавших газов не может быть ламинарным по причине структуры и характера потока и структуры пористых фильтров, где присутствуют как поперечные, так и продольные течения у стенки пористого блока.
Процесс фильтрации отработавших газов в пористых материалах можно разделить на несколько этапов, следующих один за другим. Сначала твердые частицы, которые должны быть отделены от потока отработавших газов, подводятся к поверхности фильтра, затем они примыкают к поверхности фильтра, а на следующем этапе проникают в поверхность фильтра и задерживаются в порах. Чем больше размер частицы, тем больше вероятность ее улавливания в пористой стенке.
Движение крупных частиц, например, кокса подчиняется закону Стокса, и для частиц размерами 1 106...1 104 м скорость их осаждения зависит от плотностей ртч и рог соответственно твердых частиц и сажи, отработавших газов; dm - среднего диаметра твердой частицы; цог - динамической вязкости отработавших газов:
.. = (Ртч Рог)g • dL
т 18 •Мог ■
Для очень мелких частиц в силу вступают законы броуновского движения. Эти частицы постоянно диффундируют поперек линий течения газов, поэтому степень очистки газов от мелких частиц тем выше, чем меньшая скорость прохождения отработавших газов через пористый фильтр.
Следующим фактором, который влияет на процессы улавливания твердых частиц в фильтрующем элементе, является седиментация. Она проявляется лишь тогда, когда частицы пребывают внутри достаточно времени, чтобы осесть.
Относительно отработавших газов дизелей важно, в зависимости от конструктивных особенностей, скоростного и нагрузочного режимов знать закон распределения твердых частиц по отдельным фракциям. Функция логарифмически нормального распределения массы твердых частиц по диаметрам имеет вид:
0(8) = П J exp(- (lg°;-У5о) '1 dlgS, (3)
lg с
2lg2 ст
где 5 - размер (диаметр) частицы, мкм; 550 - медиана распределения массы частиц, медианый диа-
метр, мкм, Это размер частиц, при которых масса всех т частиц мельче или крупнее; 550 составляет 50 %; Б - проход, %, доля массы т и прошедшие через СВС-блок с заданной пористостью; ^ст - стандартное среднеквадратичное отклонение величины логарифмов диаметров от их среднего значения; \ё2ст - дисперсия величины ^ст.
Интеграл, входящий в выражение (3) не может быть выражен через элемент функции. Для вычисления искомой функции ее преобразуют в функцию нового аргумента / - нормального распределения величины. Среднее значение ее равно 0.
1М 5- 1я550
t = -
lgCT
После замены:
D(8) = F (t) = F
lg ст - lg 8 5 lg ст
100
\[2ж
Jexp |- - Idt,
Эффективность захвата твердых частиц фильтром в результате диффузии описывается выражением:
Л,
(R -и-Гтч)2
Однако в этом выражении неопределенным остается коэффициент Дтч. Получение его величины связано с экспериментом для конкретных фильтров и конкретных дизелей. С другой стороны эксперимент пока остается единственным путем его определения.
Эффективность захвата частиц за счет сил инерции определяется критерием Стокса [4]:
Stk =
2r р и
тч гтч 9 ЛогЯог ’
Инерционный захват увеличивается с увеличением среднего диаметра пор d„, плотности отработавших газов, скорости фильтрации газов:
Stk2
п =-------------------------
“ Stk3 + 0,77 Stk2 + 0,22
Если представить пористый СВ С-фильтр как многослойную структуру с толщиной 8ст и толщиной слоя 8М, можно представить забивание верхнего слоя фильтра следующим выражением:
(
Ф2 =9i • exp
-и-8 (1 + r /8 )
т_________ст V_____тч сл-'
2>/3(г +8 )8
т сл сл
где тт - время прохождения газом стенки фильтра; ф - забивание фильтра начальное.
Коэффициент поглощения твердых частиц пористым СВС-фильтром вследствие зацепления можно описывать [2]:
8^ =
282 -л/3
•и + -
3r2
8ст +8„,
В порах с извилистостью, какие и есть в СВС-материалах, существуют силы:
- центробежная ^ = ттч 1вр ютч;
-трения ^ = Кст ир;
где штч - масса твердой частицы; 1ер - расстояние от частицы до оси вращения; атч - угловая скорость твердой частицы; ир - радиальная скорость.
Масса твердой частицы может быть условно представлена как:
тту = Уштч (Ртч - Рог )>
где Уштч - объем частицы при ее условно шаровой форме, УШт„ = 3п т, г,, - радиус условной шаровой частицы.
Учитывая, что ш^ =^~;Рцбс = 4лг3(д„ -рог)/1
и в случае равновесного движения ¥,^=Встур, радиальная скорость частицы, движущейся в поле центробежной силы, будет:
2Г 2 и (Ртч -Рог )
и =-------------------,
Р 91 /Л0г
а время падения твердых частиц составит:
4,5 Лог К
Т , =-------------—
пад / \ 2 ’
(Ртч -Рог ) §Г
где ку - высота установки фильтра твердых частиц. Теоретические скорости движения твердых частиц можно вычислить на основании выше приведенных зависимостей для условно принятых частиц шарообразной формы, а при испытании фильтров для очистки газов целесообразно отдельно определить скорости осаждения частиц различного происхождения.
3. Противодавления, создаваемые пористыми проницаемыми фильтрами твердых частиц на выпуске дизеля
Выполняя свои функции по очистке отработавших газов от твердых частиц, фильтры создают на выпуске значительные противодавления, которые, в свою очередь, влияют на пропускную способность, а также мощность и часовой расход топлива печью или дизелем. По мере осаждения на поверхности пористого фильтра и проникновения в поры твердых частиц возрастает сопротивление движению отработавших газов. Для пористых бумажных и тканевых фильтров очистки воздуха от пыли [5] предложено выражение для определения гидравлического сопротивления, которое по аналогии можно в первом приближении использовать и для пористых СВС-материалов:
др = 817 Лог иог (1 - Птч ) х
ё2 • п
тч тч
х(0,82-106ёОЧ\м(1 -Птч)2/3По +о0гГрКД /Рч),
где иог - средняя скорость отработавших газов, м/с; йтч - средний диаметр твердых частиц, итч - порис-
тость твердых частиц; щм - пористость фильтрующего материала; ртч - плотность твердых частиц, кг/м3;
- время между регенерациями фильтра, ч; Лог - динамическая вязкость отработавших газов; п0 - число твердых частиц в м3 газа; КД - показатель дымности отработавших газов по оптическому дымомеру.
Суммарное противодавление, создаваемое фильтром твердых частиц можно представить как сумму противодавлений, создаваемых отдельными ступенями очистки:
др _ др+дрп+ДРШ + • • •+д р..
При расчете фильтра твердых частиц с блоками из СВС-материалов предлагается использовать показатель удельной газовой нагрузки. Исходя из того, что объемный расход газов с учетом особенностей нагрузки и условий сгорания в дизеле определяется по выражениям (!) и (2), показатель удельной газовой нагрузки можно представить для первого фильтрующего элемента:
Чог _
V,
т.е. одновременно речь идет о скорости фильтрации через пористую стенку. Это нам дает возможность исходить из скоростей фильтрации, рекомендуемых в предшествующих исследованиях [6].
Предварительные экспериментальные исследования, проведенные авторами, показали, что ухудшение показателей работы дизеля для фильтров с различной пористостью происходит вначале плавно до некоторого предела дРпр, а затем показатели резко ухудшаются. Величину дРпр необходимо выбирать таким образом, чтобы превышение ее не наступало раньше времени между техническими обслуживаниями ТО-2. В противном случае сажевый фильтр оборудуется системами регенерации в процессе эксплуатации. Время между регенерациями иногда называют сажеемкостью фильтров.
Исходя из того, что для более полного поглощения фильтрами твердых частиц необходимо время пребывания газов в пористой стенке не менее т„да4,6-10“2 с, определяют объемную скорость, имеющую связь с временем пребывания в пористой стенке как ю0=3б00/г„. Из условия знания объемной скорости определяется объем пористой массы: V -т
V ___ ог п
3600 ■
При создании высокоэффективных фильтров по очистке отработавших газов пыли использованы положения, изложенные В.И. Волковым [7]. Для режима турбулентной фильтрации в существенно нелинейном режиме, когда коэффициент сопротивления фильтра квадратично зависит от числа Рейнольдса. Своими предварительными исследованиями мы обнаружили такую зависимость для фильтров из СВС-материалов.
Если СВС-материалы в некотором приближении рассматривать как зернистые среды, то за основу при описании турбулентной фильтрации можно положить модель Эргана.
Существование турбулентных пульсаций газов заложено своим дискретным принципом работы механизмов газораспределения дизелей и здесь в основу влияния может быть использована система уравнений Навье-Стокса. Решение задач, поставленных В.И. Волковым по определению критических чисел Рейнольдса, начиная с которых движение газов в порах СВС-структур становится турбулентным независимо от проходного сечения, выявление составляющих пульсаций скорости, приводящих к квадратической зависимости сопротивления на выпуске от числа Рейнольдса, приводит к очень важному положению, позволяющему решать задачу относительно полного перепада давления по толщине пористой стенки. При этом будем полагать, что рассмотрение газодинамики внутри пористой стенки базируется на двух моделях: законе Дарси при низких числах Рейнольдса и двухчленном уравнении Эргана - при высоких числах Рейнольдса, связывающих дР со среднерасходной скоростью внутри пористого слоя [5]:
др/5 _
150ЛогОог (1 - Пфм )2 , 1,75 Рог Ч2 (1 - Пфм )
пфп1„ -
где дР/5т - полный перепад давления на толщине стенки 5Ш; ик, рк, лог - соответственно расходная скорость газа в поре, плотность и динамическая вязкость газов; пфм - пористость; йш - диаметр шарика.
Уравнение Навье-Стокса и уравнение неразрывности для зернистого слоя, где координата г направлена по направлению расходной скорости, представляется в виде:
дих
~дГ
+и.
дих
дх
+иу
дП
1 дР
_ —-^- + v
Р дх
д2Пх
дх2
дП,.
дг
-+их
дП,.
дх
-+и,.
ду д2Пх + ду2 дП
+ и.
1 др
Р дУ
2
+ v
д и,
дх2
ду
д 2П,
+и.
ду2
ди,
дг
+и.
ди,
дх
1 дР
_ ——+ V
Р дг
+и,
ди
дх2
ди,
ду
д2и ,
+ и.
ду2
дих
д
д2их д,,2
и
д
д2иу д,,2
дЦ
д
д2и, д,2
(4)
Уравнение неразрывности запишется следующим образом:
дих ди ди,
-+ ^г^ + ^г,_ о.
(5)
дх ду д,
Из симметрии физической картины течения следует, что
ТТ тт диг диг
их _ иу; 17 _-^- (6)
С учетом соотношения (6) уравнения (4) можно переписать в другом виде:
dU „ r, dU гт dU х- + 2-U,-^ + U.
dt
1 др
= -—^- + v р дх
дх
dz
„ d2Ux d2Ux
2—7х + - х
дх
dz2
дt
дх
1 др ( д 2Uz
= -—^- + v\ 2—^Z+
р д.
дх
дz
д 2Uz д.2
Решение системы найдено в виде стоячих волн, а окончательно уравнение Эргана получило вид:
др ! 5 _ 178Лог иог (1 ~ Пфм )2 + 1,5 Рог Ч2 (1 " Пфм )
пф - d2
фм ш
- d„.
Последнее выражение, как наиболее верно описывающее явление оказалось удобным при расчете толщины стенки пористого СВС-фильтра. Нами при условии выполнения ограничений по противодавлениям получено выражение:
ДР. = ■
КАт^ (1 - Пфм )(2,4 -lnRe.)
Р6ПфмГ (2 - lnRej )
+HCH (C. - CCÍ)) -
dt
ср.- Мг
т =.
’ pj V
Q
dt
= OnFpJ (TJ - L )■
Общее количество тепла, теряемое 1 м2 поверхности корпуса фильтра путем лучеиспускания и конвекции можно рассчитать по выражению:
( ( j \4 / пл \4
т 1 0
100
+ Кф т - т0)1,
Mf=-
V -дс
f
1000
Теплота, получаемая вследствие сгорания СО, СН, ТЧ определяется как:
^Н _ МСННСН _ 50034 • МСН > кДж/ч;
0ТЧ _ МтчЯтч _ 34034 • Мтч, кДж/ч;
Qco _ мсоЯсо _10091 -Мсо> кДж/ч.
Учет действительного теплового состояния фильтрующих блоков дает возможность рассчитывать скорости выгорания сажи на поверхности фильтров, используя реакции:
C+O2
2C+O2
C+H2O-
-> CO2
-> 2CO,
-> CO+H2.
Скорости этих реакций зависят от температуры, а константы скоростей имеют вид:
К+1 _ 1,085-106р0 -Т4’^ехр' 164500
Это выражение можно использовать с условием создания допустимых противодавлений выпуску фильтрами нескольких ступеней j очистки.
4. Температурное состояние пористых СВС-блоков
Температурное состояние фильтров можно определить с учетом отвода тепла через стенки и подвода - с отработавшими газами и за счет доокисле-ния продуктов неполного сгорания:
T = T . + (M HCO (C1 - с2! ) + H C - с21) +
pj 1 j v ог CO v CO CO ■' тч v тч тч ■'
K+2 = 14,3 -10-4 PO T~0,5 exp
K+3 = 1,0125
P P
1 H2O -1 CO2
P P
y H2 CO y
exp
RT
93660
RT
58200
RT
где Тпк - температура поверхности стального корпуса; Т0 - температура окружающей среды; Кш - константа лучеиспускания; для листовой стали Км=2,8; Кф - конвективный фактор; для средних атмосферных условий Кф=3,5.
Часовые выбросы с отработавшими газами продуктов неполного сгорания рассчитывают по выражению:
На основе выражений, представленных в разделах 1 —3 настоящей работы, был разработан программный комплекс для моделирования течения газов через пористую проницаемую стенку.
Результаты работы позволяют определить как площадь поверхности фильтрующего элемента, так и его геометрические характеристики.
При расчетах фильтрующих элементов для химических производств вводились условия с учетом специфики, производительности источников выбросов, температуры выбросов и их химический состав.
5. Результаты моделирования
Для получения основных конструктивных размеров пористых проницаемых фильтров для очистки отработавших газов дизелей от твердых частиц была использована модель, в основу которой были положены выражения, приведенные в предыдущих разделах. Расчеты выполнены для дизелей Д-461 (6ЧН13/14) и Д-442 (4ЧН13/14). Конкретизация двигателей необходима в виду того, что двигатели имеют различный рабочий объем, различные частоты вращения и различные степени форсировки по среднему эффективному давлению, а, следовательно, различные расходы отработавших газов, различные выбросы твердых частиц, различный дискретный состав твердых частиц.
В результате моделирования и расчетов были предложены варианты основных конструктивных
т
параметров. Основные результаты расчетов приведены в таблице, где показаны варианты конструктивных элементов фильтров твердых частиц.
Таблица. Параметры дизелей и нейтрализаторов отработавших газов
Параметры дизеля и нейтрализатора Обозна- чения Размер- ность Типы дизеля
Д-461 Д-442
Полный рабочий объем цилиндров i V л 11,6 7,44
Количество цилиндров I шт. 6 4
Мощность номинальная N кВт 125,6 87,7
Часовой расход топлива Gт кг/ч 29,3 19,1
Часовой расход воздуха Q м3/ч 714 421
Частота вращения коленвала n мин-1 1700 1900
Объем пористой массы VM м3 1,22-10-2 0,73-10-2
Толщина стенок блоков sa м 12-10-2 12-10-2
Количество блоков n шт. 1 1
Противодавление выпуску ДР кПа 5,5 5,1
Извилистость пор £ - 1,1 1,1
При определении основных параметров фильтров с СВС-блоками учитывались особенности, связанные с тем, что при увеличении толщины стенки блока удельная поверхность возрастает непропорционально. Это, в частности, объясняется тем, что число сквозных пор уменьшается, появля-
г м м
ется большое количество пор с тупиками , извилистость пор увеличивается.
СПИСОК ЛИТЕРАТУРЫ
1. Смайлис В.И. Малотоксичные дизели. — Л.: Машиностроение, 1972. —128 с.
2. Евстигнеев В.В., Вольпе Б.М., Горколь ДА. Создание принципов обобщений СВС-технологий дисперсного композиционного материала // Труды АлтГТУ. Вып. 1. —Барнаул, 1993. —С. 3-18.
3. Фукс Н.А. Механика аэрозолей. — М.: АН СССР, 1956. — 352 с.
4. Коузов П.А. Основы анализа дисперсного состава промышленных пылей и измельченных материалов. — Л., 1974. — 280 с.
Учитывались положения об унификации блоков по отношению рядам двигателей для обеспечения определенной номенклатуры фильтров.
Расчеты показали, что при допустимых противодавлениях на выпуске до 5,9 кПа при одной ступени обеспечивается очистка 70 % от твердых частиц.
Конструктивные размеры блоков были положены в основу создания модификаций фильтров для очистки отработавших газов дизелей от твердых частиц.
В результате математического моделирования удалось определить объем пористой массы со средним диаметром пор 120 мкм, количество блоков, толщину их стенок, а также противодавление выпуску, создаваемое нейтрализаторами. Эти данные для дизелей связаны в таблице.
6. Выводы
В результате изучения процессов фильтрации отработавших газов дизелей через пористую стенку СВС-блока была получена математическая модель процессов фильтрации и выгорания сажи, включая блоки: определения состояния отработавших газов на входе в фильтр; процессы фильтрации газов в пористых проницаемых СВС-блоках; определение противодавлений, создаваемых пористыми блоками на выпуске дизеля и температурного состояния пористых СВС-блоков. Модель дала возможность учесть выгорание сажи на поверхности фильтров для расчета действительного теплового состояния блоков. Математическая модель реализована в виде программного для ПЭВМ и использована при оценке конструктивных параметров фильтров.
5. Фукс НА. Механика аэрозолей. — М.: АН СССР, 1956. — 352 с.
6. Груданов В.Я. Физико-химические теплообменные процессы в каталитических нейтрализаторах с утилизацией теплоты отработавших газов // Двигателестроение. —1991. —№ 1. —С. 47—49.
7. Волков В.И. Режим турбулентной фильтрации // Повышение экологической безопасности автотракторной техники: Сб. статей. — Барнаул, 1999. — С. 72—77.