О моделировании систем очистки отработавших газов ДВС с использованием нейтрализаторов и тепловых аккумуляторов фазового перехода

В связи с неблагополучным состоянием окружающей природной среды в городах Российской Федерации органы государственной власти проводят политику улучшения экологических показателей автотранспортных средств. В нашей стране с 2006 года Постановлением Правительства РФ № 609 от 12.10. 2005 г. введены в действие международные экологические требования к колесным транспортным средствам: Правила № 83 ЕЭК ООН, Правила № 49 ЕЭК ООН и Правила № 24 ЕЭК ООН, значительно ужесточившие нормы эмиссии вредных веществ с отработавшими газами (ОГ) автомобильных двигателей [1].

Известно, что одним из основных методов обезвреживания ОГ двигателей внутреннего сгорания (ДВС) автомобилей является их каталитическая нейтрализация, осуществляемая путем пропус- кания ОГ через каталитический нейтрализатор (КН). Однако на определенных режимах работы дизельного двигателя эффективность КН может быть незначительной. В табл. 1 представлены экспериментальные данные, полученные В.Н. Ложкиным при испытаниях дизеля КамАЗ-740 по 13-тиступенчатому циклу согласно Правилам ЕЭК ООН № 49 [2]. Как следует из анализа опытных данных, наиболее низкая степень очистки углеводородов С Н и оксида углерода CO наблюдается в режиме холостого хода и на малых нагрузках. Это объясняется низкой температурой ОГ, недостаточной для начала реакции нейтрализации. Поэтому возникает проблема тепловой стабилизации матрицы каталитического конвертора на эксплуатационных режимах работы двигателя. Одним из альтернативных вариантов решения этой про- блемы является использование теплового аккумулятора фазового перехода (ТАФП) [3,4]. Бортовой ТАФП обеспечивает температурную стабилизацию тепловых процессов путем «заимствования» недостающей или «сброса» излишней теплоты.

КН, совмещенный с ТАФП, состоит (см. рис.) из двух основных частей: собственно каталитического конвертора и ТАФП, объединенных в едином цилиндрическом корпусе 1. ТАФП включает в себя входной патрубок 2, диффузор 3 и теплообменник, состоящий из центральной капсулы, выполненной в виде цилин- дра, и нескольких коаксиально расположенных заполненных теплоаккумулирующим материалом (ТАМом) цилиндрических капсул 4 с образованием между ними кольцевых зазоров 5. Цилиндрические капсулы 4 и кольцевые зазоры 5 имеют одинаковые радиальные размеры (толщину). Между цилиндрическим корпусом 1 и теплообменником находится слой тепловой изоляции 6. Каталитический конвертор включает в себя блочный реактор 7, «глухую» полость расширения 8 и выходной патрубок 9.

Таблица 1 – Выборочные технические и экологические показатели дизеля КамАЗ-740 на режимах нагрузочных характеристик по Правилам ЕЭК ООН № 49 с каталитическим нейтрализатором

cd * Ри ^	N e , кВт	n дв , мин -1	т г , К	G Т , кг/ч	G B , кг/ч	G г , кг/ч	с сн ,	млн -1	л сн , %	с со ,	млн -1	п со , %
							О Щ	о R НН о С			55 о У с
1	0	800	373	1,8	159	161	80	80	0	270	210	22
2	10,2	1800	433	5,8	531	537	136	130	4	1230	730	41
3	25,4	1800	493	9,7	520	530	120	115	4	820	620	24
4	50,8	1800	573	14,7	523	538	135	125	7	840	500	41
5	76,3	1800	673	19,6	513	533	165	90	45	930	380	59
6	101,7	1800	813	25,4	510	535	380	100	74	3080	420	86
7	0	800	393	1,8	165	167	105	95	9,5	420	330	21
8	122,5	2600	833	34,2	764	798	386	90	76	1730	290	83
9	91,9	2600	738	26,9	762	789	150	75	50	780	170	78
10	61,3	2600	653	19,6	756	776	114	65	43	720	210	71
11	30,6	2600	573	14,0	747	761	105	75	29	870	180	45
12	12,3	2600	523	10,9	745	756	125	120	4	1110	1050	5
13	0	800	413	1,8	162	164	130	125	4	450	400	11

При разработке математических моделей функционирования системы КН-ТАФП использованы следующие основные допущения:

• 1.    В начальный момент времени т = 0 в период зарядки ТАФП ТАМ находится в твердом состоянии, а в период разрядки – в жидком, при этом его температура выровнена по объему и равна температуре фазового перехода Т ф .

• 2.    В обратимых процессах фазового перехода ТАМа плавление-кристаллизация при т > 0 границы раздела фаз сформированы, температурное поле ТАМа в растущей фазе линейно, а

• температура исчезающей фазы равна температуре фазового перехода ТФ .

• 3.    Теплопроводность ТАМа в продольном направлении отсутствует.

• 4.    Процесс фазового превращения ТАМа принимается одномерным и моделируется чистой теплопроводностью. При этом границы раздела фаз неизменны по форме и в каждый момент времени представляют собой цилиндрические поверхности, расположенные концентрично по отношению к цилиндрическим стенкам капсулы (вариант осесимметричной задачи).

• 5.    Коэффициенты переноса не зависят от температуры.

• 6.    Движение потока ОГ по каналам реактора рассматривается применительно к стационарному режиму работы ДВС как установившееся одномерное адиабатное течение газа.

• 7.    Для стационарного режима работы дизеля в каналах каталитического реактора устанавливается равновесный процесс между явлениями тепломассообмена и кинетики химических преобразований.

• 8.    Процесс теплообмена между стенкой канала и потоком ОГ осуществляется только по внешней поверхности активного слоя, т.е. теплообменом внутри пор пренебрегаем.

• 9.    Поток ОГ равномерно распределяется по всей совокупности каналов блочного каталитического реактора.

  
  
  

Рисунок – Принципиальная схема КН, совмещенного с ТАФП : цилиндрический корпус; 2,9 – входной и выходной патрубки; 3 – диффузор; 4 – цилиндрические капсулы с ТАМом; 5 – кольцевые зазоры; 6 – слой тепловой изоляции; 7 – блочный реактор; 8 – «глухая» полость расширения

Ниже представлена математическая модель функционирования системы КН-ТАФП на примере режима отдачи теплоты (разрядки ТАФП).

В процессе разрядки ТАФП часть отдаваемой им теплоты рассеивается в окружающей среде. Уравнение теплового баланса для ТАФП в этом случае имеет вид:

q Г <т) = Праз -q ТАФП (т) ,       (1)

где: qГ (г)  - тепловая мощность, полу чаемая потоком ОГ, Вт; qТАФП (т) - тепловая мощность, отдаваемая при разрядке ТАФП, Вт; зар – энергетический КПД процесса разрядки.

Процесс конвективного теплообмена ОГ со стенками капсул при т > 0 описывается следующим уравнением: Ттафп (О = «г F • [Тст (т) - ТСр (т)], (2) где:   г - коэффициент теплоотдачи от стенок капсул к ОГ, Вт/(м2 К); Тст(т) -температура стенок капсул в момент времени , К; F - суммарная площадь цилиндрических поверхностей всех капсул, м2; Тс(т) - средняя температура ОГ в полости ТАФП в момент времени , вычисляемая по формуле

Т_вх +Т _ВЬix(0 вх вых

Т ср ⁽ ) -         2        ,

где: Т вх , Т _ВЬ1Х (т) - температуры ОГ на входе и выходе из ТАФП, К.

Кроме того, при т > 0 справедливы следующие уравнения:

T ТАФП(т) = Х^ТЛ Т ^ф Т с ст ^(Т)

F ц ;  (4)

У СО

^т тафп СО - рТ^В • rr • У (^Т) • F u ,    ⁽⁵⁾

d где: у(т) - толщина закристаллизовавшегося слоя ТАМа в момент времени т , м;

^Т _Т ^В –  коэффициент теплопроводности

К); ^Т _Т ^В - плотность твердого ТАМа, кг/м³; r _Т - удельная теплота фазового перехода плавление-кристаллизация, Дж/кг.

Тепловая мощность q Г (т) при т > 0 определяется по балансовому уравнению q г (т) = GT се г [Твых (О" Твх ],   (6)

где: G _г – массовый расход ОГ, кг/с; с _г – удельная массовая теплоемкость ОГ, Дж/( кг∙К).

Система уравнений (1 – 6), в которых неизвестными функциями являются q Г ⁽О , q ТАФП СО , У (^, ^Т вых ⁽О , Т ст (т) , описывают процесс кристаллизации ТАМа в системе КН-ТАФП при следующих начальных и граничных условиях:

q тафп ^{( 0)} - q г ^{( 0)} - ^{0 ;}

0< у (т)< §Т; У ( 0 ) = 0 ;

Т ст ( 0 ) ~Т ф;           ^>

_           _*

Т (ft\-T

Твых(0) Т ,          J где Т – полутолщина капсулы, м.

Для аналитического решения вышеупомянутой системы уравнений введем следующие безразмерные комплек-

•< [( N- Т]р_аз + 2 ) 2+8- 1 ] ² -( N -^раз + 2 )

;

сы:

- безразмерная температура нагреваемой поверхности капсулы

Т ф - Т _ст (т) ст                      ;

Т ф Т _вх

2 N раз ст ^{ср     2 N} раз ^;

- средняя безразмерная температура ОГ в полости ТАФП

g вых

2- N" Праз' ( 1-2-%т )

2 N _раз

Т    Тс о(^т)

^{ф    ср} ;         (9)

ср

Т ф   Т _вх

- безразмерная температура ОГ на выходе из ТАФП

Они показывают, что процесс кристаллизации ТАМа в режиме чистой теплопроводности определяется тремя безразмерными параметрами – t , N и _зар .

^{Т ф   Т вых} ;     (10)

вых

Т ф   Т _вх

- безразмерная толщина закристаллизовавшегося ТАМа в режиме чистой

теплопроводности ос г- у (т)

ТВ ;

Т

- безразмерное время т-(аг)2 (Тф-Твх)

t       ТВ ТВ         ;

_{Т Т}   r _T

Рассматривая внешнедиффузионную задачу гетерогенного катализа ОГ с учетом принятых допущений, можно получить функциональную связь между температурой стенки канала каталитического конвертора, разностью концентраций

реагирующего вещества в ядре потока и у активной каталитической поверхности, а также на входе и выходе из блочного каталитического конвертора. Эта зависимость теоретически обосновывает сте-

пень конвертирования вредных веществ от теплового состояния конвертора [1, 2]

- число теплопередачи

F

N ^{г ц} .                (13)

G _г с _г

~ =~ qo^_.

ст    вх cгV

( С нач - С вых) + ( С 0 - С 1 )•

2 D 3

a _г

После тривиальных алгебраических преобразований система уравнений (1)-(6) с учетом введенных соотношений (8)-(13) и краевых условий (7) сводится к решению следующих уравнений:

ст      d ;(14)

dt ст ср ст;(15)

N раз 1 вых ;(16)

ср ст

0 ср =0,5-(0 ВыХ +1)

с начальным условием

0 при t 0 .(18)

Решением системы безразмерных уравнений (14 – 17) являются следующие функции:

где: Тст , Твх – температуры соответственно стенки канала блочного конвертора и ОГ на входе в реактор (на выходе из ТАФП), К; q – тепловой эффект реак- ции, Дж/кг; сгV – удельная объемная теплоемкость ОГ, Дж/(м3∙К); Снач, Свых, С0, С1 — концентрации реаги- рующего компонента соответственно на входе в реактор, на выходе из реактора, в ядре потока и у активной поверхности каналов реактора, кг/м3; D – истинный коэффициент диффузии, м2/с; аг – коэффициент температуропроводности реагирующего компонента, м2/с.

Рассматривая изображенный на рис. КН как теплообменный аппарат, получаем уравнение теплового баланса системы КН-ТАФП:

Q Г Q ТАФП Q КАТ Q ПОТ , (24)

где: Q _Г – общее количество теплоты, получаемое ОГ при прохождении через систему КН-ТАФП за время    , Дж;

QТАФП – количество теплоты, отдаваемое ТАФП потоку ОГ при зарядке или получаемое потоком ОГ от ТАФП при его разрядке за время , Дж;  QКАТ – количество теплоты, выделяющееся в блочном реакторе в процессе гететеро-генной каталитической реакции за время , Дж;  QПОТ – суммарные тепловые потери КН за время , Дж.

Приведенная математическая модель функционирования системы КН-ТАФП позволяет исследовать и рассчитывать процессы отдачи теплоты и гетерогенного катализа. Она является основой для инженерной методики расчета подобных систем.

Рассмотрим пример вариантного расчета ТАФП в системе КН-ТАФП для городского автобуса ЛиАЗ-5256 с механической коробкой передач и дизельным двигателем КамАЗ-7408.10. Путем численного расчета городского цикла на дороге для городских автобусов, заданного согласно ГОСТ 20306-90 графической зависимостью V (скорость) – S (пройденный путь), можно получить среднеинтегральные значения температуры и расхода ОГ, характеризующие реальный режим работы ДВС в режиме городского автобуса. Они являются наиболее вероятными величинами и представляют собой основные исходные параметры утилизации отходящей теплоты ДВC, предназначенные для теплового расчета ТАФП и некоторых других теплоутилизационных агрегатов, а также для выбора ТАМа. Так, при температуре окружающего воздуха Т 253 К они составляют: массовый расход отработавших газов - G _г ^расч = 345 кг/ч, их температура – Т _г ^расч = 510 К. Приведем величину Т _г ^расч к температуре Т 296К (при этой температуре были получены экспериментальные данные, представленные в табл.). Для этого используем приближенную формулу, полученную А.А. Сорокиным в работе [5]:

зим t _г

лет tг    кt     лет зим

I        -(Т о ТО 0   ), лет

T ₀

зим

T ₀

где: t _г ^зим , t _г ^лет – температуры ОГ соответственно в зимний и летний периоды, ºС;

Т ₀ ^зим , Т ₀ ^лет – температуры окружающего воздуха соответственно в зимний и летний периоды, К; к – коэффициент, учитывающий увеличение удельного расхода топлива при зимней эксплуатации.

Полагая к 1 , по формуле (25)

получаем: Т € _г^расч  576К при Т ^лет 296

К. Исходя из значения Т € _г^расч 576 К, выбираем ТАМ – гидроксид натрия NaOH , основные теплофизические свойства которого представлены в табл. 2.

Зададимся геометрией теплообменника, изображенного на рис.: r 20 мм, r 25 мм, r 35мм, r 40 мм, r 50 мм, r 55мм, r 65мм, r 70 мм, r 80мм, r 85мм, r ₁₁ 95 мм, r 100 мм, где r - расстояние от оси

центральной капсулы до цилиндрической поверхности «кольцевой» капсулы.

Примем, что при работе дизеля КамАЗ-7408.10 на холостом ходу ТАФП должен подогревать его ОГ, поступающие в ТАФП, от Т вх 373К до Т _вых 553 К. Отдельным расчетом опре-

деляем, что коэффициент теплоотдачи от стенок капсул к ОГ аг =62 , 4 Вт/(м²^К).

Тогда необходимая безразмерная температура ОГ на выходе из ТАФП

т 1 -т

Тф Твых вых — ^   ^   _ •••

Т ф Т_вх

572-553

... =---------- = 0 , 0955 ;

572 373

,

средняя безразмерная температура ОГ в полости ТАФП

,  ₌ ° ф -Т _С р (т)

ср

Т ф   Т _вх

572  463

572  373

= 0 , 547 ,

а максимальное значение безразмерной

толщины закристаллизовавшегося

а г- у С0_

1 ТВ     ’”

Т

_ 62 , 4-0 , 005

ее                   - 0 , 17 3 »

1 , 8

Приняв г|_раз = 0 , 9 и стно уравнения (19) и (22) _ст и N , получаем:

решая совме-относительно 0 ст = 0 , 0803,

N = 2,15. Зная 9_ст =0 , 0803, по формуле (8) вычисляем Т _ст 556 К. Из уравнения (13) вычисляем значение площади поверхности теплообмена:

_ G_rс-N _ 161-1024-2 , 15        2

F тт —            —                  — 1,58 ^^ .

ц      ос_г         3600-62 , 4

Таблица 2 – Основные теплофизические свойства гидроксида натрия NaOH

Наименование теплофизической величины и ее размерность	Численное значение
Температура фазового перехода плавление-кристаллизация, К	572
Удельная теплота фазового перехода, кДж/кг	393
Плотность в твердой фазе, кг/м3	2130
Плотность в жидкой фазе, кг/м3	1780
Коэффициент теплопроводности в твердой фазе, Вт/(м∙К)	1,8

Наконец, разрешая уравнение (20) относительно t , получаем, что безразмерное время полной кристаллизации ТАМа при выбранных параметрах составит t 0,314. Из уравнения (12) вычисляем физическое время полной кристаллизации ТАМа :

ТВ ТВ т= < • ЛТ • РТ • rT  =

(а г ) ² Т ф- Т _В х )

0 , 314-1 , 8-2130-393000  _Г11

еее = -------------------------- = 611 С е

62 , 4 2 -( 572-373 )

Исходя из принятой геометрии теплообменника и значения F 1,58 м2, оконча- тельно определяем, что его длина должна составлять 410 мм.

Таким образом, размеры ТАФП практически сопоставимы с размерами штатного глушителя шума.

Л и т ература

• 1.    Саватеев А.И. Модификация систем выпуска отработавших газов пожарных автомобилей разогреваемыми каталитическими конверторами: Автореф. … канд. техн. наук / СПб. ун-т МВД России. – СПб., 2002. – 25 с.

• 2.    Ложкин В.Н. Теория и практика безразборной диагностики и каталитической ней-трализации отработавших газов дизелей: Дисс.... докт. техн. наук /СПб. гос. техн. ун-т. – СПб., 1994. – 444 с.

• 3.    Ложкин В.Н., Шульгин В.В. Теория и практика применения тепловых аккумуляторов фазового перехода для улучшения экологических и топливно-экономических показателей автотранспортных средств // Вопросы охраны атмосферы от загрязнения: Информационный бюллетень № 2 (22). – СПб. : НПК «Атмосфера» при ГГО им. А.И. Воейкова, 2000. – С. 40-53

• 4.    Патент РФ 2204027 С1 МКИ 7F 01 N 3/00. Каталитический нейтрализатор / В.Н. Ложкин, В.В. Шульгин, С.Д. Гулин, Г.М. Золотарев (РФ). -№ 2001129630/06; Заявлено 01.11.2001; Опубл. 10.05.2003, Бюл. № 13

• 5.    Сорокин А.А. Система предпусковой тепловой подготовки двигателей строительных и дорожных машин в условиях отрицательных температур окружающей среды при возведении военностроительных комплексов: Дисс.... канд . техн. наук / Воен. инж.-строит. ин-т. – СПб., 1994. – 158 с.