Математическое моделирование образования и выгорания сажи в цилиндре дизеля 2 ч10,5/12,0 при работе на рапсовом масле и этаноле Текст научной статьи по специальности «Химические технологии»

УДК 621.436

БОТ 10.18698/2308-6033-2017-8-1649

Математическое моделирование образования и выгорания сажи в цилиндре дизеля 2Ч 10,5/12,0 при работе на рапсовом масле и этаноле

© В. А. Лиханов, А.Н. Козлов, М.И. Арасланов

Вятская государственная сельскохозяйственная академия, г. Киров, 610017, Россия

Рассмотрен один из перспективных способов снижения дымности отработавших газов дизеля — использование альтернативных топлив, применение которых позволяет решить задачу по замещению нефтяного топлива. Представлена математическая модель образования сажи в цилиндре дизеля при сгорании рапсового масла и этанола, основанная на химической кинетике пиролиза углеводородов топлива и высокотемпературном ацетиленовом механизме. На базе современных представлений о механизме образования зародышей частиц сажи в условиях камеры сгорания дизеля предложен химизм сажеобразования при работе на рапсовом масле и этаноле. Результаты проведенного расчета уровня дымности отработавших газов дизеля 2Ч 10,5/12,0 на рапсовом масле и этаноле показывают, что при работе на альтернативных топливах на номинальном режиме концентрация сажи в отработавших газах снижается в 3,7раза по сравнению с использованием дизельного топлива. Полученные результаты теоретических расчетов удовлетворительно согласуются с экспериментальными данными.

Ключевые слова: дизель, рапсовое масло, этанол, камера сгорания, пиролиз, двойная система топливоподачи, сажа, дымность отработавших газов

Изучение процессов сажеобразования в камере сгорания (КС) дизелей связано прежде всего с постоянно ужесточающимися требованиями к экологическим характеристикам энергетических установок. Известно, что в результате высокой излучательной способности сажи значительно возрастают потери энергии в цилиндре дизеля в виде теплового излучения пламени. Часть тепла, выделяемая при сгорании дисперсного углерода, отводится конвективным путем. Таким образом, образование сажи снижает индикаторный КПД дизеля, в том числе за счет неполноты и несвоевременности окисления углеводородов топлива [1].

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

ность сажи [2]. Известно, что канцерогены, адсорбированные на поверхности твердых частиц, действуют на живые клетки активнее, чем канцерогены в чистом виде [3]. В образцах сажи из выхлопной системы или отобранных при низких температурах отработавших газов (ОГ) обнаруживаются тяжелые углеводороды, сконденсированные в агломераты и превращенные в смолистые вещества. Попадая в легкие человека, частицы сажи вызывают раздражение, способствуют появлению злокачественных новообразований. Кроме того, исследования показали, что и сельскохозяйственные растения способны накапливать сажу в значительном количестве [4].

Один из перспективных способов снижения дымности дизелей — это применение альтернативных топлив [5, 6]. Для оптимизации процессов сгорания и максимального снижения выбросов сажи при работе дизеля на альтернативном топливе необходимо иметь полное представление об основных стадиях быстропротекающих процессов в цилиндре. Для этого, в частности, необходимо разработать адекватную математическую модель образования и выгорания сажи в цилиндре дизеля, описывающую основные стадии сажеобразования от стадии термического разложения углеводородов исходного топлива и их производных до реакций газификации на поверхности частиц сажи. Такая модель должна не только с достаточной точностью позволять рассчитывать показатели сажесодержания в ОГ, но и давать представление о динамике процессов сажеобразования в цилиндре дизеля при работе на альтернативных видах топлива — этаноле и рапсовом масле (РМ).

Образование и выгорание частиц сажи в КС дизеля — сложные физико-химические процессы, интенсивно протекающие в короткие промежутки времени в нестационарных температурных условиях, при переменных значениях давления и концентрации компонентов газовой смеси.

По температурным условиям, наличию кислорода, концентрации углеводородов топлива объем КС дизеля неоднороден. При работе дизеля на РМ и этаноле с двойной системой топливоподачи в КС можно условно выделить несколько зон.

В зонах КС дизеля с дефицитом кислорода воздуха при высоких температурах углеводороды топлива подвергаются пиролизу с образованием молекул с меньшим количеством атомов углерода. По мере роста температуры устойчивыми остаются молекулы с более короткими цепями, а также циклические молекулы и радикалы бензольного ряда. При температурах свыше 1500 °С остаются только молекулы К2, О2, Н2, Н2О, СО2, СО, С2Н2 и некоторые углеводородные радикалы, устойчивые к распаду.

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

В зонах с недостаточной концентрацией углеводородного топлива, где температура на границах факелов впрыскиваемых РМ и этанола понижена вследствие испарения спирта, процесс термического разложения углеводородов замедляется. В них возможно развитие холоднопламенных реакций, ведущих к образованию и осаждению на стенках цилиндра кислородсодержащих молекул, в частности пер-оксидов и альдегидов [7].

На участках факела спирта с достаточной концентрацией кислорода при распространении фронта пламени происходит выгорание спирта с образованием углекислого газа и воды по реакции

С2Н50Н + 302 ^ 2С02 +3Н20. (1)

Процесс пиролиза молекул РМ происходит в две стадии: декар-боксилирование и разложение олефина. Декарбоксилирование — начальная стадия термического распада молекул карбоновых кислот РМ и происходит при более низких температурах, поэтому скорость брутто-реакции распада молекулы РМ в зоне высокотемпературного пиролиза лимитирована скоростью разложения олефина, образовавшегося в результате декарбоксилирования молекулы карбоновой кислоты РМ.

Для определения скорости разложения углеводородов РМ необходимо провести химическую аппроксимацию процессов декар-боксилирования и пиролиза одной мономолекулярной реакцией термической конверсии карбоновой кислоты РМ с произвольно ограниченным числом стабильных компонентов:

СхН7С00Н ^ С02 + У1С3Н6 + у2С2Н4 + У3С2Н2 + у4СН4, (2)

где СхН7 С00Н — условная молекула карбоновой кислоты РМ; у I — стехиометрический коэффициент.

Стехиометрические коэффициенты в уравнении (2) были определены с учетом индуктивного и мезомерного эффектов, наблюдаемых при термическом распаде олефинов.

Заключительный этап глубокого пиролиза — взрывное разложение ацетилена с образованием химических зародышей частиц сажи, достаточно подробно рассмотренное в монографии [8]. Общий вид

уравнения реакции разложения ацетилена может быть записан следующим образом [9]:

С2Н2 ^ С2Н + Н. (3)

Науглероживание радикала зародыша сажи происходит при взаимодействии с ацетиленом за счет присоединения атомов углерода и усиления их связей. При достижении числа атомов углерода 80.. .100 и более появляется физическая поверхность сажистой частицы [10]. Таким образом, происходит фазовый переход к твердой конденсированной фазе по следующей брутто-реакции:

СтН2 + С2Н2 ^ С(т+2)Н2 + Н2, (4)

где т — количество атомов углерода в зародыше сажи.

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

При работе дизеля на РМ и этаноле, даже в зонах с богатой топ-ливовоздушной смесью, параллельно с образованием и ростом сажевых частиц происходит их выгорание из-за наличия атомов кислорода, поступающего в КС вместе с топливом. Принято, что процесс газификации твердого углерода сажевых частиц может происходить за счет взаимодействия с радикалами, содержащими кислород, молекулами воды, монооксидом и диоксидом углерода, а также за счет реакции со свободным кислородом. Химический механизм включает в себя гетерогенные и гомогенные реакции окисления, идущие на поверхности и в пограничном объеме углеродной частицы [9]:

С + О2 ^ СО2; (5)

С + 0,5О2 ^ СО; (6)

С + СО2 ^ 2СО; (7)

С + 2Н2О ^ СО + 2Н2; (8)

2СО + О2 ^ 2СО2; (9)

2Н2 + 02 ^ 2Н20. (10)

Полное сажевыделение дизеля при работе на этаноле и РМ с двойной системой топливоподачи (ДСТ) складывается из сажи, образовавшейся по различным механизмам в пламени богатой топливно-воздушной смеси (ТВС) этанола и запального РМ (рис. 1).

Химические зародыши сажевых частиц

Физические зародыши сажевых частиц С^о_1оо^Н2

+ С2Н2 + С6Н5

+ со2 + н20

Процессы коагуляции и поверхностного роста

+ 02 + СО + ОН

Процессы газификации и выгорания

Процессы агломерации и агрегации

Содержание сажи в цилиндре дизеля при работе на РМ и этаноле с ДСТ

Рис. 1. Химизм сажеобразования при работе дизеля на этаноле с использованием ДСТ с воспламенением от запальной струи РМ

На сегодняшний день в научной среде существует несколько гипотез о ключевом механизме процесса сажеобразования. В условиях подачи в КС двух топлив, различных по своим физико-химическим свойствам, в цилиндре дизеля характерно наличие значительной температурной и концентрационной неоднородности топлив и кислорода в объеме КС, следовательно, образование сажи в нем потенциально может идти с помощью различных механизмов.

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

образования сажи (через полициклические ароматические углеводороды), предполагаем, что эти реакции не вносят существенного вклада в конечный выход сажи. С некоторой долей вероятности можно говорить о частичном фракционировании негорящих капель в топливном факеле РМ и на поверхностях КС. Однако допускаем, что этот процесс не вносит существенного вклада в суммарное сажеобразование. Предполагаем, что для этанола не характерен процесс сажеобразования путем высокотемпературной термической полимеризации и дегидрогенизации парожидкостного ядра испаряющихся капель, предложенный Н.Ф. Раз-лейцевым [11], ввиду высокой скорости испарения.

Окончательно процессы образования и газификации сажи в цилиндре дизеля 2Ч 10,5/12,0, работающем на этаноле и запальном РМ, можно представить в упрощенном виде как систему брутто-реакций:

С17Н33СООН К1 > СО2 + 2С3Н6 + 4С2Н4 + С2Н2 +СН4, С2Н5ОН К2 > С2Н4 + Н2О, С3Н6 ——^ С2Н4 + 0,5С2Н2 + 0,5Н2, СН4 К4 > 0,5С2Н4 +Н2,

С2Н4 —К5 > С2Н2 +

С2Н2-——> С2Н+ Н

СтН2 +С2Н2-С(т+2)Н2 + Н2

Сп + С2Н2-^^ Сп+2 + Н2,

С + О2 ——^ СО2, С + О2 К10 > СО, С + СО2 К11 > СО, С + Н2О К12 > СО + Н2, Н2 + О2 —^ Н2О,

(11)

где К1 — константа скорости 1-й реакции; Сп — количество атомов углерода на поверхности частицы.

Для аналитического описания процесса сажевыделения используем методы решения прямой задачи химической кинетики. Обобщенную химическую модель процесса сажевыделения описывает система уравнений (11). Для определения скорости образования химических зародышей частиц сажи необходимо провести кинетический

расчет первых пяти уравнений, проходящих в газовой фазе. Для этого составим систему уравнений для расчета мгновенных концентраций основных компонентов газовой смеси, участвующих в образовании ацетилена и влияющих на скорость сажеобразования в зонах богатой ТВС:

ёМ1 _ к1М1 ёф 6п йМ2 _ к2М2 ёф 6п

ёМ3 1

ёф 6п ёМ4 1

_—(2кМ - къМъ),

ёф 6п ёМ5 1

_—(к1М1 - к4М4),

ёф 6п ёМ6 1

_ —(4к1М1 + к2М2 + к3М3 + 0,5к4М4 - к5М5),

(12)

ёф 6п ёМ7 1

_ —(0,5км3 + к4М4 + к5М5),

_ — (0,5к3М3 + к1М1 + к5М5),

ёф 6п йМ8 _ к2М2 ёф 6п ёМ9 _ к1М1 ёф 6п

где М1 — текущая концентрация /-го компонента газовой смеси; ф — угол поворота коленчатого вала двигателя относительно верхней мертвой точки (ВМТ); к1 — константа /-й реакции системы (11); п — частота вращения коленчатого вала двигателя, мин-1. Здесь индексы текущей концентрации системы (12) соответствуют следующим компонентам газовой смеси: 1 — РМ; 2 — этанол; 3 — пропилен; 4 — метан; 5 — этилен; 6 — водород; 7 — ацетилен; 8 — вода; 9 — углекислый газ. Константы скоростей реакций могут быть вычислены с использованием справочных данных.

Для связывания процессов образования и горения компонентов газовой смеси в условиях КС дизеля используем уравнение [12], учитывающее подвод компонентов с испаренным топливом, химические преобразования и выгорание в процессе основного горения:

_(м) (м) (м)

ёф I dф I ёф ) I ёф ^

(13)

где 1М1 / ёф — суммарная скорость изменения концентрации /-го компонента газовой смеси; (ёМ/ /ёф) — скорость изменения концентрации /-го компонента, связанного с подводом испаренного топлива; (ёМ / ёф)к — скорость изменения концентрации /-го компонента вследствие предпламенных химических превращений топлива; (ёМ/ / ёф)х — скорость изменения концентрации /-го компонента в

процессе основного горения.

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

((1x1 (ск)обр (^выг , (14)

где (с1Г / ёх) — суммарная скорость изменения сажесодержания в цилиндре дизеля; ((г / 1х)о6р — скорость образования сажи;

((г / ёт)выг — скорость выгорания сажи; N — текущее значение

сажесодержания в цилиндре, кг.

Скорость сажеобразования в предпламенной зоне пропорциональна концентрации ацетилена [12] и определяется по выражению

(ёЫ! _ _ 107 крРГС2Н2Гобр^уд (15) — _ В1-;-, (15)

( ёф)обр 6п

где В1 — эмпирический коэффициент; кр — константа скорости поверхностного роста; р — текущее давление в цилиндре, МПа; гС2Н2 —

объемная концентрация ацетилена, определяемая по уравнению материального баланса; Гобр — текущее значение сажесодержания в

зоне образования, кг; £уд — удельная поверхность частицы сажи.

Скорость выгорания частиц сажи в условиях дизельного пламени, согласно работе [7], определяем по выражению

'(г! _ ^У^Гехр

ёф) Уа

V т /выг V

Е2 Т

КГт ( Г

+2к ехр

~КГт ( Г*

(16)

где В2 — эмпирический коэффициент; qв — цикловая подача воздуха, кг/цикл; а — суммарный коэффициент избытка воздуха; % — доля

выгоревшего топлива; V — текущий объем цилиндра, м3; к *, Т — координаты полюса химических реакций по зависимости Аррениуса; Е1, Е2 — энергия активации, кал/моль, соответственно для реакций (5) и (6); Я — универсальная газовая постоянная; Тт — среднемассовая температура сажистых частиц.

Константа скорости поверхностного роста частиц сажи вычисляется по данным П. А. Теснера [12]:

кр _ 17,1-10-4 ехр(-138 000/ЯТ). (17)

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

Т + Т

-1 г 1 -1 т

Тобр

Тобр _ г ^ т . (18)

2

Под среднемассовой температурой Тт понимается среднеинте-гральная температура частиц, распределенных в объеме, осредняю-щая действительную температуру частиц в зависимости от их местонахождения от Тт;п в пристеночной зоне до Ттах в пламенной зоне. За определяющую температуру процесса выгорания принимаем среднемассовую температуру Твыг _ Тт .

По данным исследований [12], температура частиц сажи в пламени находится в диапазоне 2600... 2800 К, что значительно превышает осредненную температуру газов. Очевидно, что такие высокие температуры частиц нельзя распространять на весь излучающий объем. Для определения Тт необходимы данные о распределении температуры частиц в объеме КС. Выгорание частиц сажи происходит преимущественно в зоне продуктов сгорания, поэтому характер изменения Тт в цикле соответствует изменению температуры пламени Тп . Тогда упрощенно Тт можно определить как

Твыг _ Тт _ Тп -ДТ, (19)

где Тп — температура пламени; ДТ — величина, характеризующая снижение среднемассовой температуры частиц сажи относительно Тп .

Величина ДТ на основании экспериментальных данных составляет 300.400 К. Отношение температуры пламени к температуре газов [7] в общем виде можно определить функцией

Тп _ 1,2 а0,3

Тг

( р ( с

V рм У

с

V тн У

х ехр-

2,4

ф/ - 5,5ф2тах (ц / qцМ )

0,333

5,5ф2тах (£7ц / qцМ ) ,

0,333

(20)

где р — среднее индикаторное давление; р М — среднее индикаторное давление на режиме максимального крутящего момента; ст — средняя скорость движения поршня; Стн — средняя скорость движения поршня на номинальном режиме работы двигателя; ф ^ — текущий угол поворота коленчатого вала (п.к.в.), отсчитываемый от начала видимого горения; ф 2 тах — угол 2-го максимума тепловыделения; qц — цикловая подача топлива; qцМ — цикловая подача топлива на режиме максимального крутящего момента.

Начальное содержание реагирующих веществ примем пропорционально их объемной концентрации в исходной смеси г и количеству топлива а(фн), поданного в цилиндр к моменту ф_фн. Величина а(ф н) определяется законом топливоподачи и скоростью испарения топлив, а значения начальных концентраций РМ и этанола соответственно равны

1н _ 1/М; (21)

Г2н _и / М2 ,

(22)

где М2 — сумма молей газовой смеси; и — соотношение количества молей этанола, поданного в цилиндр дизеля, к одному молю запального РМ.

Основываясь на принятых предпосылках, уравнениях (12) и (14), математическую модель сажевыделения, в зависимости от угла п.к.в. дизеля при работе на РМ и этаноле с использованием ДСТ, с учетом подвода испаренных топлив в газовую смесь и выгорания компонентов смеси, можно записать в виде системы дифференциальных уравнений:

М Л А ¿о к1М1

1 = агМормГ1Н —1 --^ -М1

¿ф ам2

ёф дМ,

ёф 6п

, ✓ до2 к2М2 = агМоэГ2н —- -М2

1 ¿X °!-Х ¿ф' 1 ¿X

ёф

6п

1

3 =-—(2кМ1 - кзМз)-Мз

°!-Х ¿ф 1 ¿X

ёф 6п дМ4 1

4 =—(к1М1 - к4М4) - М4

ёф 6п

°!-Х ¿ф 1 ¿х °!-Х ¿ф'

= — (4к1М1 + к2М2 + к3М3 + 0,5к4М4 -кМ5)-М5

ёф 6п

1 ¿х О-Жр'

= — (0, ^к3М3 + к4М4 + кМ5) - М6 ёф 6п

=—(0,5кМ3 + к1М1 + кМ5) - М7 ёф 6п

¿М8 к2М2

1 ¿х

°!-Х ¿ф' 1 ¿X

°1-Х ¿ф'

(23)

ёф 6п ¿М9 = к1М1 ёф 6п

Г Л аг

= К г ъ

ёф

N ёф

N ёф

обр Л

= N

У выг

Л Г <т

Л 1 ¿ф

а ёф

+В

10 крРгс2н2

уд

6п

- N

обр

1 ¿X

°!-Х ¿ф'

а-аг ¿о^ а ёф

р I ¿ф-^

+В2 NSуд А+#обр

1 ¿х °!-Х ¿ф'

где аг — локальный коэффициент избытка воздуха в зоне горения; М0РМ, М0э — стехиометрическое количество воздуха, необходимого для полного сгорания одного моля условной молекулы РМ и эта-

та; о1, о 2, — относительная доля испарившегося РМ и этанола в цилиндре дизеля и приведенный по теплоте сгорания суммарный подвод испаренных топлив соответственно; Ын — начальное саже-содержание в цилиндре; А — функция, введенная для упрощения записи уравнения выгорания сажи в системе.

нола соответственно; гн — начальная концентрация 1-го компонен

Учитывая известное процентное содержание атомов кислорода, водорода и углерода по массе, средний состав РМ можно приближенно выразить в виде формулы молекулы олеиновой кислоты С18Н34О2. В соответствии с принятой аппроксимированной формулой молекулы уравнение реакции полного выгорания запального топлива примет вид

C18H34O2 +25,5O2 ^ 18CO2 + 17H2O. (24)

На основании уравнений (1) и (24) определим необходимое количество молей воздуха для сгорания одного моля топлива для этанола и РМ:

М0э = 3/0,21 = 14,28; (25)

М0РМ = 25,5/0,21 = 121,42. (26)

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

N = У Г ЧвСд, (27)

где уг — коэффициент остаточных газов; Сд — действительная

концентрация сажи в ОГ.

Численно интегрируя систему дифференциальных уравнений (20) с заданным вектором начальных условий, получаем мгновенные концентрации основных компонентов газовой смеси в реакционной зоне. Значение концентрации сажи в момент открытия выпускного клапана принимаем за результирующую концентрацию сажи в ОГ дизеля, т. е. не учитываем газификацию сажи на участке выпускной коллектор — пробоотборник. Для сравнения с результатами экспериментального измерения дымности ОГ полученная расчетная концентрация может быть переведена в шкалу Bosch [13] по эмпирической формуле

Г 10 ^1206

[С1 = 5651 ln-I , (28)

L J l 10 - Bosch/

где [С ] — расчетная концентрация сажи в ОГ.

На основе представленной модели были проведены теоретические расчеты показателей сажесодержания в цилиндре дизеля 2Ч 10,5/12,0 при работе на этаноле и запальном РМ.

В представленном расчете при работе дизеля на альтернативных топливах на номинальном режиме:

Математическое моделирование образования и выгорания сажи в цилиндре дизеля... Установочные углы опережения впрыскивания

топлив п.к.в. до ВМТ © , град........................................34

Частота вращения п, мин-1..................................................1800

Среднее эффективное давление ре , МПа.........................0,588

Цикловые подачи топлива, мг/ цикл:

для этанола дэ..................................................................52

для запального РМ дРМ ..................................................13

По результатам расчетов видно, что процесс пиролиза топлив инициируется в момент начала тепловыделения (рис. 2). До начала тепловыделения происходит постепенное увеличение концентрации паров топлива, связанное с продолжающимся испарением.

Рис. 2. Концентрации углеводородов топлива и продуктов пиролиза в реакционной зоне КС при работе дизеля 2Ч 10,5/12,0 на этаноле и РМ на номинальном режиме работы дизеля:

1 — Н2; 2 — С2Н2; 3 — СН4; 4 — С3Н6; 5 — С2Н4; 6 — этанол; 7 — РМ

С развитием кинетического сгорания и прогревания топливовоз-душной среды увеличивается скорость мономолекулярных реакций распада топлива, в связи с этим увеличивается концентрация продуктов пиролиза. Вместе с тем интенсивно растет расход топлива, связанный с выгоранием. После 5° п.к.в. после ВМТ концентрация газифицированного топлива (РМ и этанола) в реакционной зоне КС минимальна. Несмотря на то, что процесс испарения продолжается, испарившееся топливо быстро подвергается термическому разложению и выгоранию.

С ростом температуры в КС дизеля увеличивается скорость разложения продуктов пиролиза углеводородов топлива с образованием ацетилена. С этим связано смещение пика концентрации ацетилена на более поздний угол п.к.в. Ацетилен образуется в ходе пиролиза в КС до конца процесса испарения топлива, т. е. до 25° п.к.в. после ВМТ на номинальном режиме работы дизеля. Постепенное снижение концентрации ацетилена на более поздних углах связано с окислением и сажеобразованием. В ходе этих гетерогенных реакций выделяется водород, максимальная концентрация которого смещена на более поздний угол п. к.в. относительно пика концентрации ацетилена. Последующее снижение концентрации Н2 объясняется постепенным преобладанием процесса окисления в связи с диффузией кислорода из других зон КС и снижением интенсивности образования водорода в ходе пиролиза углеводородов топлива и сажеобразования.

С некоторой задержкой относительно начала тепловыделения (4...6° п.к.в.) в КС дизеля при сгорании этанола и РМ начинается процесс интенсивного сажеобразования (рис. 3). Скорость газификации частиц сажи на начальном этапе минимальна преимущественно по причине относительно низкой температуры в зоне продуктов сгорания. По мере увеличения осредненной температуры в КС скорость газификации сажи увеличивается. После 20° п.к.в. при работе дизеля на рассматриваемом режиме процесс газификации сажи начинает преобладать над процессом образования, и массовое содержание сажи —расч в цилиндре дизеля уменьшается.

Ы, мг 0,2 0,1 0

ВМТ 20 40 60 80 100 <р°, п.к.в

Рис. 3. Изменение массового содержания сажи -расч в цилиндре дизеля 2Ч 10,5/12,0 на номинальном режиме работы на этаноле и РМ:

1 —расч; 2 —вых. расч

Снижается и скорость сажеобразования, поскольку уменьшается концентрация ацетилена в предпламенной зоне. К моменту открытия выпускного клапана (при ф = 126° п.к.в.) массовое содержание сажи

в цилиндре дизеля стабилизируется и становится равным -вых. расч.

Предложенная математическая модель позволяет рассчитать концентрации сажи и основных компонентов газовой смеси в реак-

ционной зоне КС на протяжении всего процесса сгорания, а также содержание сажи в ОГ на различных скоростных и нагрузочных режимах работы дизеля. Вместе с тем данная модель позволяет получить ценную информацию о динамике протекания основных этапов образования и выгорания сажи в цилиндре при работе дизеля на этаноле и РМ.

Расчетное значение концентрации сажи в цилиндре дизеля при открытии выпускного клапана составляет Свых. расч = 0,114 г/м , что соответствует уровню дымности ОГ, равному 2,2 ед. по шкале Bosch, — это в 3,7 раза ниже, чем при работе двигателя на дизельном топливе.

Результаты расчетов дымности ОГ дизеля 2Ч 10,5/12,0 при работе на РМ и этаноле по приведенной математической модели удовлетворительно соответствуют полученным экспериментальным данным.

ЛИТЕРАТУРА

[1] Матиевский Д. Д., Свистула А.Е. Влияние переменности состава топливо-воздушной смеси в зоне горения на сажевыделение, параметры рабочего цикла и индикаторный КПД цикла дизеля. Ползуновский вестник, 2002, № 1, с. 10-16.

[2] Матвеев С.Г., Орлов М.Ю., Чечет И.В., Семенов А.В. Взаимосвязь между выбросами сажи и бенз(а)пирена тепловыми двигателями на углеводородном топливе. Вестник Самарского государственного аэрокосмического университета, 2009, № 3, с. 210-214.

[3] Шаромов И.М., Зонов А.В., Чупраков А.И. Образование бенз(а)пирена и зародышей сажи. Улучшение эксплуатационных показателей двигателей внутреннего сгорания. Материалы IV Междунар. науч.-практ. конф. «Наука — Технология — Ресурсосбережение». Сб. науч. трудов. Киров, Вятская ГСХА, 2011, с. 154-159.

[4] Денисова Е.С. Аккумуляция некоторыми сельскохозяйственными растениями техногенной пыли сажевых заводов. Омский научный вестник, 2014, № 2 (134), с. 196-199.

[5] Торопов А.Е. Улучшение экологических показателей дизеля 4Ч 11,0/12,5 при работе на метаноло-топливной эмульсии путем снижения дымности отработавших газов. Дис. ... канд. техн. наук. Киров, 2010, 172 с.

[6] Россохин А.В. Результаты исследований влияния применения этаноло-топливной эмульсии на показатели процесса сгорания и сажесодержания в цилиндре и отработавших газах дизеля 4Ч 11,0/12,5 в зависимости от нагрузки. Молодой ученый, 2015, № 15 (95), с. 171-174.

[7] Ложкин В.Н. Исследование динамики и термических условий сажевыделения при сгорании распыленного топлива в цилиндре дизелей. Дис. ... канд. техн. наук. Ленинград, 1978, 228 с.

[8] Лиханов В.А., Шаромов И.М. Снижение дымности отработавших газов дизеля 4Ч 11,0/12,5 путем применения этаноло-топливных эмульсий. Киров, Вятская ГСХА, 2012, 144 с.

[9] Петриченко Р.М., Батурин С.А., Исаков Ю.Н., Пугачев Б.П. Элементы системы автоматизированного проектирования ДВС. Алгоритмы прикладных программ. Ленинград, Машиностроение, 1990, 328 с.

[10] Россохин А.В. Особенности расчета сажевыделения в дизелях. Материалы VI Междунар. науч.-практ. конф. «Наука — Технология — Ресурсосбережение». Сб. науч. трудов. Киров, Вятская ГСХА, 2013, с. 196-199.

[11] Разлейцев Н.Ф. Моделирование и оптимизация процесса сгорания в дизелях. Харьков, Высшая школа, 1980, 169 с.

[12] Батурин С.А. Физические основы и математическое моделирование процессов сажевыделения и теплового излучения в дизелях. Дис. ... д-ра техн. наук. Ленинград, 1982, 441 с.

[13] Кулешов А.С. Программа расчета и оптимизации двигателей внутреннего сгорания ДИЗЕЛЬ-РК. Описание математических моделей, решение оптимизационных задач. Москва, Изд-во МГТУ Н.Э. им. Баумана, 2004, 123 с.

Статья поступила в редакцию 03.04.2017

Ссылку на эту статью просим оформлять следующим образом: Лиханов В.А., Козлов А.Н., Арасланов М.И. Математическое моделирование образования и выгорания сажи в цилиндре дизеля 2Ч 10,5/12,0 при работе на рапсовом масле и этаноле. Инженерный журнал: наука и инновации, 2017, вып. 8. http://dx.doi.org/10.18698/2308-6033-2017-8-1649

Лиханов Виталий Анатольевич — д-р техн. наук, профессор, заведующий кафедрой тепловых двигателей, автомобилей и тракторов Вятской государственной сельскохозяйственной академии (Вятской ГСХА). Область научных интересов: улучшение эксплуатационных показателей дизелей путем применения альтернативных видов топлива, в первую очередь природного газа и спиртов; теоретические и экспериментальные исследования рабочих процессов, протекающих в цилиндрах автотракторных дизелей при работе на альтернативных видах топлива, и образования основных токсичных компонентов; разработка систем дозирования и регулирования; создание новых смесевых топлив на основе спиртов. Автор более 600 научных работ. e-mail: lihanov.va@mail.ru

Козлов Андрей Николаевич — ассистент кафедры тепловых двигателей автомобилей и тракторов Вятской ГСХА. Область научных интересов: улучшение показателей работы дизеля путем применения альтернативных топлив. e-mail: dnka59@mail.ru

Арасланов Марат Ильдарович — ассистент кафедры тепловых двигателей автомобилей и тракторов Вятской ГСХА. Область научных интересов: улучшение показателей работы дизеля путем применения альтернативных топлив. e-mail: araslanov.89@mail.ru

Mathematical modeling of soot generation and burnout in a 2Ch 10,5/12,0 diesel cylinder operating on rapeseed oil and ethanol

© V.A. Likhanov, A.N. Kozlov, M.I. Araslanov Vyatka State Agricultural Academy, Kirov, 610017, Russia

We consider one of the promising methods of decreasing diesel exhaust gas opacity, namely using alternative fuels, employing which makes it possible to replace petroleum fuel. We present a mathematical model for soot generation in a diesel cylinder during rapeseed oil and ethanol combustion. The model is based on the chemical kinetics for pyrolysis of fuel hydrocarbons and a high-temperature acetylene mechanism. We suggest a soot generation chemism for diesels operating on rapeseed oil and ethanol, based on contemporary understanding of the soot particle precursor formation mechanism in diesel combustion chambers. Calculation results for the exhaust gas opacity level of the 2Ch 10,5/12,0 diesel using rapeseed oil and ethanol show that while operating on alternative fuels under nominal conditions, soot concentration in exhaust gases decreases by 3.7 times as compared to using diesel fuel. The results of theoretical computations obtained are in satisfactory agreement with the experimental data.

Keywords: diesel, rapeseed oil, ethanol, combustion chamber, pyrolysis, dual fuel supply system, soot, exhaust gas opacity

REFERENCES

[1] Matievskiy D.D., Svistula A.E. Polzunovsky vestnik — Polzunovsky vestnik (Herald of I.I. Polzunov Altai State Technical University), 2002, no. 1, pp. 10-16.

[2] Matveev S.G., Orlov M.Yu., Chechet I.V., Semenov A.V. Vestnik Samarskogo gosudarstvennogo aerokosmicheskogo universiteta — Vestnik (Herald) of the Samara State Aerospace University. Aerospace and Mechanical Engineering, 2009, no. 3, pp. 210-214.

[3] Sharomov I.M., Zonov A.V., Chuprakov A.I. Obrazovanie benz(a)pirena i zarodyshey sazhi. Uluchshenie ekspluatatsionnykh pokazateley dvigateley vnutrennego sgoraniya [Benzo[a]pyrene and soot precursor formation. Improving performance characteristics of internal combustion engines]. Materialy IV Mezhdunar. nauch.-prakt. konf. Nauka — Tekhnologiya — Resursosberezhenie. Sb. nauch. trudov [Proc. of the 4th International scientific and applied conference, Science — Technology — Resource Conservation]. Kirov, Vyatka State Agricultural Academy Publ., 2011, pp. 154-159.

[4] Denisova E.S. Omskiy nauchnyy vestnik — Omsk Scientific Bulletin, 2014, no. 2 (134), pp. 196-199.

[5] Toropov A.E. Uluchshenie ekologicheskikh pokazateley dizelya 4Ch 11,0/12,5 pri rabote na metanolo-toplivnoy emulsii putem snizheniya dymnosti otrabotavshikh gazov. Diss .... kand. tekhn. nauk [Improving environmental peformance of the 4ChN 11.0/12.5 diesel running on methanol/fuel emulsion by means of decreasing the exhaust gas opacity. Cand. eng. sc. diss.]. Kirov, 2010, 172 p.

[6] Rossokhin A.V. Molodoy uchenyy — Young Scientist, 2015, no. 15 (95), pp. 171-174.

[7] Lozhkin V.N. Issledovanie dinamiki i termicheskikh usloviy sazhevydeleniya pri sgoranii raspylennogo topliva v tsilindre dizeley. Diss. ... kand. tekhn. nauk [Investigating the dynamics and thermal conditions of soot generation during combustion of fuel sprayed in diesel cylinders. Cand. eng. sc. diss.]. Leningrad, 1978, 228 p.

[8] Likhanov V.A., Sharomov I.M. Snizhenie dymnosti otrabotavshikh gazov traktornogo dizelya 4Ch 11,0/12,5 putem primeneniya etanolo-toplivnykh emulsiy [Reducing exhaust gas opacity of the 4ChN 11.0/12.5 tractor diesel by means of using ethanol/fuel emulsion]. Kirov, Vyatka State Agricultural Academy Publ., 2012, 144 p.

[9] Petrichenko R.M., Baturin S.A., Isakov Yu.N., Pugachev B.P. Elementy sistemy avtomatizirovannogo proektirovaniya DVS. Algoritmy prikladnykh programm [Elements of computer aided design systems for internal combustion engines. Software application algorithms]. Leningrad, Mashinostroenie Publ., 1990, 328 p.

[10] Rossokhin A.V. Osobennosti rascheta sazhevydeleniya v dizelyakh [Specifics of soot formation computations in diesels]. Materialy VIMezhdunar. nauch.-prakt. konf. Nauka — Tekhnologiya — Resursosberezhenie. Sb. nauch. trudov [Proc. of 6th International scientific and applied conference, Science — Technology — Resource Conservation]. Kirov, Vyatka State Agricultural Academy Publ., 2013, pp. 196-199.

[11] Razleytsev N.F. Modelirovanie i optimizatsiya protsessa sgoraniya v dizelyakh [Modeling and optimising the combustion process in diesels]. Kharkiv, Vysshaya Shkola Publ., 1980, 169 p.

[12] Baturin S.A. Fizicheskie osnovy i matematicheskoe modelirovanie protsessov sazhevydeleniya i teplovogo izlucheniya v dizelyakh. Diss. . dokt. tekhn. nauk [Physical foundations and mathematical modeling of soot generation and heat emission processes in diesels. Dr. eng. sc. diss.]. Leningrad, 1982, 441 p.

[13] Kuleshov A.S. Programma rascheta i optimizatsii dvigateley vnutrennego sgoraniya DIZEL-RK. Opisanie matematicheskikh modeley, reshenie optimizatsionnykh zadach [DIZEL-RK software for computing and optimising internal combustion engine parameters. Description of mathematical models, solutions to optimisation problems]. Moscow, BMSTU Publ., 2004, 123 p.

Likhanov V.A., Dr. Sc. (Eng.), Professor, Head of Department of Heat Engines, Automobiles and Tractors, Vyatka State Agricultural Academy. Specialises in improving performance characteristics of diesel engines by using alternative fuels, primarily natural gas and alcohols; theoretical and experimental studies of work cycles of automotive diesel engine cylinders running on alternative fuels; theoretical and experimental studies of primary toxic component formation; developing metering and regulating systems; formulating new alcohol-based mixed fuels. Author of over 600 scientific publications. e-mail: lihanov.va@mail.ru

Kozlov A.N., Assist. Lecturer, Department of Heat Engines, Automobiles and Tractors, Vyatka State Agricultural Academy. Specialises in improving performance characteristics of diesel engines by using alternative fuels. e-mail: dnka59@mail.ru

Araslanov M.I., Assist. Lecturer, Department of Heat Engines, Automobiles and Tractors, Vyatka State Agricultural Academy. Specialises in improving performance characteristics of diesel engines by using alternative fuels. e-mail: araslanov.89@mail.ru