CC BY
92
Features of the use of interactive programs like GAS2 to simulate the injection process in reciprocating internal combustion engines are considered. The results for the velocity field in the combustion chamber and the engine cylinder space are obtained.
Key words: internal combustion engine, the method of large particles, the interactive software solutions for the computer.
Lukin Sergei Eduardovich, master, manager of software department, lu-kinse82@yandex. ru, Russian, Tula, Tula State University,
Potapov Sergei Alexandrovich, candidate of technical sciences, docent, la-ma0 71 @mail. ru, Russia, Tula, Tula State University.
Chesnokov Sergei Alexandrovich, doctor of technical sciences, professor, sa-chesn45@,mail. ru, Russia, Tula, Tula State University
УДК 621.436
ОБРАЗОВАНИЕ ТОКСИЧНЫХ ОКСИДОВ АЗОТА И УГЛЕРОДА В АВТОМОБИЛЬНЫХ ДВИГАТЕЛЯХ
С.Э. Лукин, С.А. Потапов, С.А. Чесноков
Рассмотрены особенности химических реакций горения и образования токсичных оксидов азота и углерода в ДВС. Получены результаты по содержанию токсичных компонентов.
Ключевые слова: дизельный двигатель, коэффициент избытка воздуха, образование токсичных оксидов в двигателях.
Двигатели внутреннего сгорания (ДВС) вносят существенный вклад в загрязнение окружающей среды. Они являются одним из главных источников токсичных веществ, выбрасываемых в атмосферу, особенно в крупных городах.
Повышение экологической чистоты рабочих процессов в ДВС требует всестороннего изучения, но особое внимание следует уделять химической составляющей процесса горения.
Рассмотрим особенности химических реакций горения и образования токсичных оксидов в ДВС [1].
1. Внутри фронта горения текущий состав смеси является существенно неравновесным. Это требует решения задач химической кинетики горения для сотен реакций.
2. На выходе фронта при высокой температуре горения (> 2000 К) состав продуктов сгорания (кроме N0) устанавливается практически равновесным.
3. За фронтом в объеме камеры сгорания (КС) при турбулентном перемешивании струй оттока различного состава (горючих и окисляющих) и более низких температурах инициируются лишь медленные реакции догорания. Их дополняют очень медленные и поэтому сильно не уравновешенные реакции образования токсичного оксида N0 из инертного газа -азота воздуха.
Механизмы образования оксидов N0 и СО. Анализ результатов кинетических расчетов в монографии [1] показал, что основными активными компонентами в реакциях N0 и СО являются радикалы О и ОН.
В настоящее время для токсичных компонентов общепринятыми являются:
1) кинетический механизм образования оксида N0 в объеме КС N + 0 « N0 + N N + 02 « N0 + 0; N + 0H « N0 + ^
2) механизм догорания в объеме КС образовавшегося во фронте оксида СО
га + oн « га2 + 0; га + 0 « га2; га + 02« га2+0.
В первом приближении содержание оксида СО, как и других углеводородов, при расширении в цилиндре будем считать равновесным.
Равновесные концентрации (мольные доли) любого компонента определяются тремя параметрами: температурой, давлением и коэффициентом избытка воздуха. Равновесное содержание компонентов легко вычисляется по типовым программам. Температура и давление при горении и расширении в цилиндре предварительно рассчитываются в рамках двух-зонной модели горения. Коэффициент избытка воздуха а задается соотношением массы воздушного заряда в цилиндре и цикловой массы впрыскиваемого топлива.
При выборе топлива стремятся обеспечить высокую температуру горения, что обуславливает большое значение термического К.П.Д. цикла двигателя. Однако при больших температурах возрастает роль диссоциации (распада) крупных молекул СО2 и Н2О с образованием токсичного оксида углерода СО и водорода Н2.
Выбор коэффициента избытка воздуха а также достаточно противоречивый.
1. При а > 1 и избытке воздуха в горючей смеси наблюдается полное сгорание топлива, токсичный оксид СО выгорает и превращается в безвредный углекислый газ СО2, но этот же избыток воздуха (и кислорода) ускоряет окисление азота воздуха с образованием более токсичного (~ в 40 раз) оксида азота N0.
2. При а < 1 и избытке топлива для заметного окисления азота кислорода недостаточно, доля N0 стремится к нулю, но наблюдается неполное сгорание топлива с образованием значительного количества СО.
В работах В.А. Звонова [2], Р.З. Кавтарадзе [3] и др. рассмотрены три механизма образования оксидов азота.
1. «Термический», предложенный Я.Б. Зельдовичем, в котором ведущую роль играют свободные атомы кислорода и азота:
N + N0 = N + 0; N + 02 = N0 + 0; N + ОН = N0 + Н. (1)
2. «Быстрый» при активном участии углеводородных радикалов, предложенный Фенимором. Ведущей реакцией в нем является реакция связывания молекул азота радикалом СН, содержание которого внутри фронта горения быстро увеличивается (на несколько порядков), а затем резко падает [1]:
СН + N = HCN + N. (2)
Образовавшиеся азотные соединения быстро окисляются до N0 по механизму, указанному Гаррисом:
HCN + Н = CN + Н2; HCN + 0Н = CN + Н20;
CN + 02 = га + N0; CN + 0Н ® га + КН; (3)
КН + 0Н ® N0 + Н2; МН + N0 ® N + 0Н.
Константы скорости последних трех реакций определены Баульхом.
3. «Быстрый» через образование закиси азота ^0, предложенный Мальте:
N + 0 + М = ^0 + М; ^0 + 0 = N0 + N0; ^0 + 0 = N + 02.
Параметры механизма реакций (1) - (3) приведены в таблице.
Детальный кинетический механизм образования оксида N0
№ п/п Реакции (... + 0, кДж/моль) Прямые Обратные
А, моль, л, с п Е, КДж / моль А, моль, л, с п Е, кДж/ моль
1 N + N0 = N + 0 + 314 1,32+10 0 0 0,69+11 -0,026 314,5
2 N + 02 = N0 + 0 + 134 1,81+5 1,5 25,0 0,24+05 1,56 158,4
3 N + 0Н = N0 + Н + 165 4,2+10 0 0 0,14+14 -0,57 207,9
4 Ш + N = HCN + N - 8,4 8,0+8 0 23,3 1,5+11 -0,43 16,3
5 HCN + Н = CN + Н2 - 87,6 1,2+12 -0,6 51,9 3,2+9 0 10,5
6 HCN + 0Н = CN + Н20 - 26,5 2,0+8 0,6 10,5 2,9+6 1,15 -1,2
7 CN + 02 = ТО + N0 + 458 3,2+8 0 0 1,3+8 0,11 229,5
8 CN + 0Н ® ТО + Ж + 229 6,0+10 0 0 - - -
9 ЫН + 0Н ® N0 + Н2 + 303 5,0+10 0 0 - - -
10 ЫН + N0 ® N + 0Н + 411 1,2+10 0 0 - - -
Коэффициенты А, п, Е определяют константы скорости реакций: к = АТпехр( Е / ЯТ).
2. Тестовый расчет содержания N0 (и СО) в цилиндре ДВС.
Содержание любых углеводородов, в том числе оксида СО, считаем равновесными. Их мольные доли однозначно определяются коэффициентом избытка воздуха, температурой и давлением.
Для оксида N0 примем, что ведущими являются три реакции Зельдовича (1). Теодорчук и Рихтер [4] показали, что систему из трех уравнений химической кинетики можно свести к уравнениям, в которых константы скорости обратных реакций заменены константами равновесия, а затем равновесными долями компонентов:
сСхш /Л = [-Р(851 + 52 + S3) + S1 + 5(52 + 53)]ЯТ/р; (4)
сХн / Ж = [-8(Р51 + 52 + S3) + S1 + Р( 52 + 53 )]ЯТ / р, (5)
где в = х^о [N0]; 5 = хN /Г^]; х^о и хN - неравновесные доли, [N0] и [Ы] - равновесные доли. Скорости прямых реакций (1) - (3) из таблицы использованы в виде
51 = ^[N0]^]; ^2 = к2[02] №; 53 = к3[0Н] [N]. Имеющихся данных достаточно, чтобы проинтегрировать систему уравнений (4) - (5). В то же время их можно свести к одному уравнению [4], предположив квазистационарность процесса образования азота N поскольку этот компонент смеси образуется и расходуется только в рассматриваемых реакциях (11.1). Из условия
Жхн / Ж = 0,
следует выражение
5 = 51 + в( 52 + 53) /(Р5 + 52 + 53) = (5 + в)/(1 + РБ), (6)
где 5 = 51/(52 + 53).
Дальнейшее упрощение связано с преобразованием правой части уравнения (4). Вынесем из выражения в квадратных скобках сумму (5 2 + 53) и подставим в (6):
5 + Р
1 + Р5
(1-Р5) + 5 - Р
(52 + 53).
Откуда имеем
25 - 25Р 2 1 + Р5
(52 + 53) = 251(1 - Р 2)
V 2 3 1 + Р5
Подставим последнее выражение вместо квадратной скобки в уравнение (4), получим
^ = 2 (1 - Р2)^. (7)
Л р Ч + Р5
Здесь Р = XNо /[N0], 5 = 51 /(52 + 53). При решении уравнения (7) для условий ДВС необходимо ввести в программу зависимости для температуры Т и давления р, полученные при расчете рабочего процесса ДВС, а также аппроксимацию равновесного состава углеводородов.
Общие результаты. Расчеты проводились для одноцилиндрового дизельного двигателя ТМЗ-450Д-03 производства ОАО «Тульский маши-
3 3
ностроительный завод»: рабочий объем цилиндра 0,454-10" м , диаметр цилиндра 0,085 м, ход поршня 0,08 м, номинальная мощность 8,1 кВт, степень сжатия 20. Доля остаточных газов в цилиндре 5 %, коэффициент избытка воздуха аср = 1 и 3.
Расчеты на ЭВМ общей кинетической задачи для двухзонной модели горения в ДВС показали, что основная часть N0 образуется в зоне продуктов сгорания по «термическому» механизму (1). Генерация «быстрых» N0 происходит только во фронте горения, когда в смеси появляются радикалы СН. Доля «быстрых» N0, образующихся через закись азота ^0, пренебрежимо мала. При обеднении топливовоздушной смеси (до а = 3) уменьшение образования оксида азота по «термическому» механизму происходит из-за снижения температуры в зоне продуктов сгорания, а «быстрых» - из-за уменьшения концентрации и времени существования радикалов СН. Выход «быстрых» N0 при снижении уровня температур уменьшается не так сильно, как «термических», поэтому их относительная доля несколько увеличивается.
Приближенные расчеты по формуле (7) качественно подтверждают данные ЭВМ, однако погрешность в определении конечного содержания оксида N0 в отработавших газах достигает 30 %.
Отметим, что для любых режимов работы ДВС, при температурах < 2200 К, реакции образования и распада оксида N0 замедляются («замораживаются») настолько, что его доля остается неизменной - («закаливается»). Аналогичный эффект наблюдается при догорании СО - его содержание стабилизируется, если температура падает ниже 1200 К [1].
Список литературы
1. Чесноков С.А. Химический турбулентный тепломассообмен в двигателях внутреннего сгорания. Тула: Изд-во ТулГУ, 2005. 466 с.
2. Звонов В.А. Токсичность двигателей внутреннего сгорания. М.: Машиностроение, 1981. 160 с.
3. Кавтарадзе З.Р., Кавтарадзе Р.З. Анализ механизмов образования и методов расчета концентрации оксидов азота в поршневых двигателях // Транспорт на альтернативном топливе. 2011. № 6. С. 65-71.
4. Teodorczyk A., Rychter T.J. Mathematical-model of nitric oxide formation in an Si piston engine // Journal of Technical Physics. 1987. Vol 28. P. 47-66.
Лукин Сергей Эдуардович, магистр, нач. отдела программного обеспечения, [email protected], Россия, Тула, Тульский государственный университет,
Потапов Сергей Александрович, канд. техн. наук, доц., lama071 @mail.ru, Россия, Тула, Тульский государственный университет,
Чесноков Сергей Александрович, д-р техн. наук, проф., sachesn45@,mail.ru, Россия, Тула, Тульский государственный университет
FORMATION OF TOXIC NITROGEN AND CARBON OXIDES
IN THE ENGINE
S.E. Lukin, S.A. Potapov, S.A. Chesnokov
The features of chemical reactions of combustion and formation of toxic nitrogen oxides and carbon in the engine are considered. Results on the content of toxic components are obtained.
Key words: diesel engine, the excess air ratio, the formation of toxic oxides engines.
Lukin Sergei Eduardovich, master, manager of software department, lu-kinse82@yandex. ru, Russian, Tula, Tula State University,
Potapov Sergei Alexandrovich, candidate of technical sciences, docent, la-ma0 71 @mail. ru, Russia, Tula, Tula State University,
Chesnokov Sergei Alexandrovich, doctor of technical sciences, professor, sa-chesn45@,mail. ru, Russia, Tula, Tula State University
