CC BY
39УДК 621.43:681.5
DOI: 10.30987/article_5cda64cebb6ba2.01228392
А.А. Обозов, М.А. Новиков, Е.В. Дмитриевский
РАСЧЕТ ХАРАКТЕРИСТИК ТЕПЛОВЫДЕЛЕНИЯ ПО ЗАДАННОМУ ЗАКОНУ ПОДАЧИ ТОПЛИВА И ПРОВЕРКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ НА ДИЗЕЛЬНОМ ДВИГАТЕЛЕ
Приведены результаты расчета характеристик тепловыделения по заданному закону подачи топлива по разработанной программе в среде Microsoft Excel. Получено экспериментальное подтверждение математической модели взаимосвязи
подачи топлива и тепловыделения в цилиндре дизельного двигателя ,3243E.
Ключевые слова: дизельный двигатель, программа расчета, закон подачи топлива, характеристики тепловыделения, опытная проверка программы.
A.A. Obozov, M.A. Novikov, E.V. Dmitrievsky
CALCULATION OF HEAT GENERATION ACCORDING TO SPECIFIED LAW OF FUEL FEEDING AND SIMULATOR VERIFICATION BASED ON DIESEL ENGINE
The results of a heat generation characteristic calculation according to a specified law of fuel feeding on the program developed in Microsoft Excel environment are shown. There is obtained an experimental confirmation of the model of the interconnection be-
Характеристики процесса сгорания в дизельных двигателях внутреннего сгорания с объемным и объемно-пленочным смесеобразованием зависят от динамики и структуры факелов распыленного топлива, их числа и ориентации в пространстве, взаимодействия факелов с воздушным зарядом цилиндра и стенками камеры сгорания. Это необходимо учитывать при разработке математической модели смесеобразования и сгорания в дизельных двигателях, а также при решении практических вопросов по интенсификации, оптимизации и управлению характеристиками процесса сгорания [1-6]. При этом достоверность математической модели и эффективность практических решений зависят от того, насколько точно учтены особенности динамики и структуры факелов распыленного топлива в дизельном двигателе.
Методы экспериментального исследования характеристик распыливания топлива дизельными форсунками не дают полной информации о динамике и структуре нестационарной импульсной струи топлива, впрыскиваемой в высокотемпературную среду с переменной плотностью,
tween fuel feeding and heat generation in the cylinder of D243E diesel engine.
Key words: diesel engine, calculation program, law of fuel feeding, heat generation characteristics, program test verification.
особенно при начавшемся горении. Такие исследования проводятся, как правило, на специальных безмоторных моделирующих установках (дизельных бомбах). Так как измерить параметры топливного факела непосредственно в цилиндре работающего дизельного двигателя технически трудно из-за высоких давлений и температур цикла, излучения пламени, нестабильности впрыска топлива от цикла к циклу и др., известны только отдельные попытки изучения динамики и структуры топливных факелов непосредственно на работающем дизельном двигателе или его модели путем киносъемки факела и стробоскопического отбора проб газового топлива в различных точках поперечного сечения факела.
Для исследования процессов распы-ливания топлива в дизельных бомбах применяют фотографирование и киносъемку топливного факела, отбор топлива из различных его зон с помощью сотовых уловителей или трубчатых зондов, улавливание капель распыленного топлива на законченную пластинку. Фото- и кинорегистрация дает представление о геометрии топливного факела и динамике ее изменения, но почти не раскрывает внутреннюю
структуру в связи с его высокой оптической плотностью. Вследствие этого практически не дают дополнительной информации о структуре топливного факела дизельных двигателей и более сложные методы: шлирен-метод и голография. Отбор топлива зондами и впрыск на пластинки позволяют определить распределение потоков жидкого топлива на различных расстояниях от оси факела или сопла и дисперсность распыливания, суммарную за период впрыска. Таким образом, измерения также не дают представления о внутренней структуре факела в различных его сечениях и динамике ее изменения. Полезную информацию о дисперсности распыленного топлива в определенном сечении факела можно получить с помощью установки со стробоскопическим приспособлением в виде вращающегося диска со щелью, через которую в определенный момент впрыска только часть факела попадает на улавливающую пластинку, расположенную за диском. Однако получаемые таким методом данные недостаточно достоверно отображают действительную структуру факела в его поперечном сечении, так как взаимодействие топливной струи с диском стробоскопа нарушает нормальное развитие факела, а в связи с различием скоростей капель в его поперечном сечении диск стробоскопа вырезает из факела объем в форме эпюры скоростей капель. Вследствие этого основную часть массы топлива, улавливаемого на пластинку, составляют капли, летящие с большой скоростью в плотном осевом потоке струи.
В многочисленных исследовательских работах отечественных и зарубежных авторов, в том числе в [1; 3-5], рассмотрены кинетические уравнения испарения и выгорания распыленного топлива в цилиндре дизельного двигателя, учитывающие характеристики впрыска и распыли-вания топлива, параметры заряда цилиндра, физико-химические характеристики топлива, особенности его испарения в ди-
зельном двигателе, кинетику цепных реакций в неоднородной топливовоздушной смеси. На основании аналитических и эмпирических зависимостей протекания указанных процессов разработан метод расчета процесса сгорания топлива и характеристик тепловыделения в дизельных двигателях [5].
По рассмотренной методике была разработана программа расчета характеристик тепловыделения по заданному закону подачи топлива в среде Microsoft Excel. Выполненный расчет для дизельного двигателя Д243Е подтвердил хорошее согласование расчетных характеристик тепловыделения и индикаторных диаграмм с экспериментальными данными.
Исходные данные для расчета (дизельный двигатель Д243Е):
Молекулярная масса воздуха m=28,9 кг/кмоль.
Диаметр цилиндра D=0,110 м.
Ход поршня £=0,125 м.
Постоянная КШМ À=R/L=0,2717.
Степень сжатия е=16.
Количество сопловых отверстий форсунки ic= 6.
Диаметр сопловых отверстий dc=0,2
мм.
Длина свободного пролета факела Lr=0,038 м.
Характеристики топлива:
Теоретически необходимое количество воздуха для сгорания 1 кг топлива M0=0,5 кмоль.
Плотность топлива (при 323 К) рт = 825 кг/м3.
Динамический коэффициент вязкости топлива (при 323 К) Цт = 3 10-3 Пас.
Коэффициент поверхностного натяжения топлива (при 323 К) = 28 10-3 Н/м.
Коэффициенты расчетных формул :
Константа времени испарения крупных капель Az =4,6 c-1.
Коэффициенты уравнения безразмерной характеристики впрыска топлива
о = ХФ): a=3,9; è=1,73; с=3,5; mi=1,5; m2=2; даз=2; m4=4; Ф1=0,333. Коэффициенты уравнения характеристики тепловыделения на участке топливоподачи:
ТУ=1.
a0= 410 м3/(кг с); a1=10-3 с; Коэффициент турбулизации заряда
«2=15 м3/(кгс); ¿0=0,1.
Минимальное значение степени использования воздушного заряда цилиндра ^во=0,42; её абсцисса Ф20=0,33. 66
Условная энергия активации пред-пламенных реакций £а=23,8-103 кДж/кмоль.
Коэффициент в формуле расчета среднего диаметра капель Ез2=1,7.
Коэффициент поправочной функции константы испарения топлива >=1.
Коэффициент в формуле определения дальнобойности топливного факела
Аф=2,9.
Коэффициент в формуле угла конуса топливного факела ^ст=0,0055.
Степень уменьшения скорости испарения топлива на стенке (минимальная) /0 =0,7.
Параметры режима работы двигателя:
Частота вращения вала двигателя п =1500 мин-1.
Давление продувочного воздуха (абс.) Р, =0,09 МПа.
Температура наддувочного воздуха Т, =315 К.
Коэффициент наполнения Пу=0,96.
Коэффициент остаточных газов 7=0,05.
Коэффициент избытка воздуха при сгорании ор=1,86.
Продолжительность впрыска топлива ^впр=29°пкв.
Момент (угол поворота коленчатого вала до ВМТ) начала впрыска топлива ^н=17°пкв.
Цикловая подача топлива
Сц=1,247 10-4 кг.
Недожог топлива (в долях от цикловой подачи) ^т=0,01.
Момент (угол поворота коленчатого вала) окончания расчета тепловыделения по уравнениям участка топливоподачи ^к=90°пкв.
Расчет характеристик тепловыделения 1. Расчет параметров заряда:
Рабочий объем цилиндра V =жБгБ /4
м .
м .
Объем камеры сжатия у = у /(s -1),
Количество свежего заряда M =п -103рУ /(8,312Т), кмоль.
см IV .Г S П \ ~ s /
Количество газов в цилиндре M = M (1 + v), кмоль.
V см \ ' /
Цикловая подача топлива (если в исходных данных Gц=0).
Коэффициент избытка воздуха при сгорании (если в исходных данных Gц=0). 2. Расчет констант испарения топлива:
Средняя скорость истечения топлива из распылителя форсунки
и = 24Чп , м/с. 0 0,75р 7Ч т
' ' т ее' впр
Критерий М, М = £ /рёс Ст ). Критерий Вебера = ир а /С ■
0 ' т е ~ т
Плотность заряда в конце условного сжатия р = тм / V, кг/м3.
Симплекс р = р /р .
Средний поверхностный диаметр капель
^ = 106 ЕъгёМ °'0733 /(рЖв )0266, мкм. Давление в цилиндре в конце услов-
ного сжатия
Рс = PsS
МПа.
Константа испарения топлива К = 1/(106 р ), м2/с.
Относительная константа испарения топлива (теоретическая) Ь = 1012 К / <, с-1.
и.т. 32
Поправка, учитывающая фактические условия испарения, у — у/ = уП • 10 3 =1.
Относительная константа испарения (действительная) Ь = УЬ , с-1.
и и.т.
Расчет временных и угловых параметров диаграммы
Объем цилиндра в момент начала впрыска топлива,
Ун = Ус (1 + (s- 1)k/2),
м
где
ая= 1 - cos ( + — (1 - cos(arcsin (Л sin ())).
Л
Давление газов в момент фн
Р» = Р, ((У + У,)/ Ун Г > МПа. Температура газов в момент фн Т = 103 рУ /(8,312M ), К.
н J. н н \ 7 V /
Период задержки воспламенения
I T Ea
т = 3,8 • 10 6 (1 -1,6 • 104 n) Te8 312Тн, с;
' \Рн
А( = 6пт ,°пкв.
Момент воспламенения
( =( +А( ,°пкв
Момент окончания впрыска топлива р = р + р ,°пкв.
f к.в. f н f впр
Продолжительность испарения и сгорания крупных капель z = A ¡(b а 0,6), с.
и z \ и р '
Продолжительность сгорания
р = р —Ар.+ 6nz ,°пкв.
~ я i впр i i и
Момент окончания горения
Р«, =Ре ,°пкВ.
Момент достижения факелом стенки камеры сгорания
z = dMl2 p¡(d UWe 0'21M016), с;
ст ф ст* \ c U /
р = р + 6nz ,°пкв.
i ст ' н ст
Продолжительность взаимодействия фронта факела со стенкой
р= 6nA 2z JpWe0 32 ¡M0 07 ,°пкв.
' ф ст ст V '
Коэффициент в формуле, определяющей характер испарения топлива на стенке,
ст =(1 — Хо )¡ 0,485. Расчет характеристик впрыска и испарения топлива
Текущий ход поршня
S = L + Я — (RCos р + 4 L — Я2 (Sin р)2). Относительное время впрыска (рис.
1)
Ф = (Р — Рн )/(пр .
Доля впрыснутого топлива:
(
если 0 < Ф < Ф ^ — = a
v mi
1 фт1 —b фm2 m
Л
иначе — = — Ф).
1
1
— фт3--фт 4
v т3 т4
Интегральная характеристика впрыснутого топлива (рис. 2)
а = ах +а2
Текущее время от начала впрыска
Т = {Р-РН)/6п •
Квадрат относительного диаметра испаряющейся капли
В=1-Ьт при В < 0 ^ В = 0 •
Степень уменьшения скорости испарения топлива на стенке:
если р<рст ^ Х = 1-
Относительное время взаимодействия фронта факела со стенкой:
если р>р ^ Ф = (р-р )/р•
Т Т ст ст V т т ст ' ' ф
Степень уменьшения скорости испарения топлива на стенке:
если р>р ^ г = 1 - с •
' т ст т
Скорость испарения топлива (рис. 3):
da
если ря <р<ркв ^ — ~ = -(1 — Въп)Х, . dz z
если р <р<р ^ — = — ((1 — К(z — zmpГ — B3 2)z
dz z
впр
Доля испарившегося топлива (рис. 4)
р
— = ¡(d- ¡ dz)d^¡6n.
Рис. 1. Дифференциальные (dol/dФ, do2/dФ, doi/dФ) и интегральные (a¡, а2, a¿) характеристики процесса впрыска
Рис. 2. Характеристика впрыска топлива
-2 = с
т, 1-1 О^С 01 23456789 10
^ ~ ^ л ¡л Рис. 4. Доля впрыснутого топлива от, доля
Рис. 3. Скорость испарения топлива аси /ах ^ ^ :
испарившегося топлива ои, характеристика
тепловыделения Х
Доля паров, образовавшихся за период задержки воспламенения,
1 = °и к ] при т = т.
Расчет характеристик тепловыделения на участке топливоподачи
Коэффициенты расчетных формул:
А = а0; А = «1; А = а2п°'5.
Объем цилиндра в момент воспламенения (если т = к) :
ag = 1 - cos ( +1 (1 - cos(arcsin (Л sin ( ))) ;
Л
У. = Ус (1 + (s-1)a./2) .
Тепловыделение при начальной вспышке топлива (если (>( и x0 < a):
- функция выгорания паров топлива, образовавшихся за т (рис. 5),
-р
о CL,
-2000—1
X
-10
-5 0 5
Угол поворота коленчатого вала , град
10
Рис. 5. Функция скорости выгорания паров топлива (образовавшихся за период т) Р0 и функция Р2
P0 = AG (аш-Х0Jb0*,+ x)/У.; скорость тепловыделения
f ^ da"
dx l „ da — = l p +—и dT v dT
1 + AiP +
dT
- доля выгоревших паров топлива, образовавшихся за период задержки воспламенения,
*=1 л 'к;
к 6п
- общая доля выгоревшего топлива
к ах йт
х = I---.
к йт 6п
Тепловыделение после вспышки топлива (если х0 >аи. и к <к^к):
- скорость догорания топлива
Р2 = АО„ (ар — х\аи — х) / V;
- относительная скорость тепловыделения
йх йт
da „I .( . da
—и+p I /11+A—-
dT ) v dT
доля выгоревшего топлива
= í
dx d(. dT 6n
v
Расчет характеристик тепловыделения на участках развитого горения и догорания
Коэффициент h = 1 /(л/2<Ф0 )•
Коэффициент с =(1 -£0)/0,485 •
Параметры точки фк.
хк = х\рк ] - значение х при р = рк ;
Фк = (Рк -Рв ) /Р, ;
^ = 1 - 2СвИФк /(Ж ехр(ИФк )2)
а =Е а / х ;
тк ^ ек р к'
А3 =^х / ¿А / (атк (1 - А т - Хк )хк )
Показатели процесса сгорания на участке рк<р<ркг:
относительная продолжительность
горения
степень использования воздушного
заряда
£ = 1 - 2с ИФ /(4Ж ехр((НФ )2)); - текущее значение коэффициента избытка воздуха в зоне горения
а =£а /х;
т ^ е р
- относительная тепловыделения (рис. 6)
dx
скорость
dz
= A«„(l -Am - x)x '
6)
доля выгоревшего топлива (рис. 4, dx dp
x = x, + j
dz 6n
Рис. 6. Интегральная характеристика тепловыделения (доля выгоревшего топлива)
Таким образом, разработанная программа расчета характеристик тепловыделения по заданному закону подачи топлива в среде Microsoft Excel позволяет получить исчерпывающие количественные значения характеристик тепловыделения в цилиндре
двигателя. Выполненный расчет для дизельного двигателя Д243Е подтвердил хорошее согласование расчетных характеристик тепловыделения и индикаторных диаграмм с экспериментальными данными
СПИСОК ЛИТЕРАТУРЫ
1. Алемасов, В.Е. Основы теории физико-химических процессов в тепловых двигателях и энергетических установках / В.Е. Алемасов [и др.]. - М.: Химия, 2000. - 520 с.
2. Двигатели внутреннего сгорания: учебник: в 3 кн. / В.Н. Луканин, И.В. Алексеев, М.Г. Шатров [и др.]; под ред. В.Н. Луканина. - М.: Высш. шк., 1995.
3. Дмитриевский, Е.В. Электронное управление двигателями внутреннего сгорания: учеб. пособие / Е.В. Дмитриевский, А.А. Обозов, М.А. Новиков. - Брянск: БГТУ, 2018. - 195 с.
4. Иванченко, Н.Н. Рабочий процесс дизелей с камерой в поршне / Н.Н. Иванченко, Б.Н. Семенов, В.С. Соколов. - Л.: Машиностроение, 1972. - 228 с.
5. Разлейцев, Н.Ф. Моделирование и оптимизация процесса сгорания в дизелях / Н.Ф. Разлейцев. -Харьков: Изд-во Харьков. ун-та, 1980. - 169 с.
6. Рогалев, В.В. Курсовое проектирование двигателей внутреннего сгорания: Тепловой расчет: учеб. пособие /В.В. Рогалев, Е.В. Дмитриевский. - Брянск: БГТУ, 2017. - 131 с.
1. Alemasov, V.E. Fundamentals of Theory of Physical-Chemical Processes in Heat-Engines and Power Plants / V.E. Alemasov [et al.]. - M.: Chemistry, 2000. - pp. 520.
2. Internal Combustion Engines: manual in 3 Vol. / V.N. Lukanin, I.V. Alexeev, M.G. Shatrov [et al.]; under the editorship of V.N. Lukanin. -M.: Higher School, 1995.
3. Dmitrievsky, E.V. Electronic Control of Internal Combustion Engines: manual / E.V. Dmitrievsky, A.A. Obozov, M.A. Novikov. - Bryansk: BSTU, 2018. - pp. 195.
4. Ivanchenko, N.N. Working Process of Diesel Engines with Chamber in Piston / N.N. Ivanchenko, B.N. Semyonov, V.S. Sokolov. - L.: Mechanical Engineering, 1972. - pp. 228.
5. Razleitsev, N.F. Modeling and Optimization of Combustion Procedure in Diesel Engines / N.F. Razleitsev. - Kharkov: Kharkov University Publishers, 1980. - pp. 169.
6. Rogalev, V.V. Course Design of Internal Combustion Engines: Heat Calculation: manual / V.V. Rogalev, E.V. Dmitrievsky. - Bryansk: BSTU, 2017. - pp. 131.
Статья поступила в редакцию 1.03.19. Рецензент: д.т.н., профессор Брянского государственного инженерно-технологического университета
Заикин А.Н.
Статья принята к публикации 23. 04. 19.
Сведения об авторах:
Обозов Александр Алексеевич, д.т.н., доцент, профессор кафедры «Тепловые двигатели» Брянского государственного технического университета, e-mail: [email protected] .
Новиков Максим Александрович, аспирант кафедры «Тепловые двигатели» Брянского государ-
ственного технического университета, e-mail: novikovmax 14@yandex. r .
Дмитриевский Евгений Владимирович, к.т.н., доцент кафедры «Тепловые двигатели» Брянского государственного технического университета, email: Dmitr-bs tu@ yandex.r .
Obozov Alexander Alexeevich, Dr. Sc. Tech., Assistant Prof. of the Dep. "Heat Engines", Bryansk State Technical University, e-mail: [email protected]. Novikov Maxim Alexandrovich, Post graduate student of the Dep. "Heat Engines", Bryansk State Tech-
nical University, e-mail: novikovmax [email protected]. Dmitrievsky Evgeny Vladimirovich, Can. Sc. Tech., Assistant Prof. of the Dep. "Heat Engines", Bryansk State Technical University, e-mail: [email protected].
