Построение математической модели сгорания топлива в цилиндре

 

МИНИСТЕРСТВО ОБРАЗОВАНИЯ  И НАУКИ УКРАИНЫ

Одесская национальная морская  академия

 

 

Кафедра Высшей математики

 

 

 

Реферат

С дисциплины : Математическое моделирование математических процессов

На тему : Построение математической модели сгорания топлива в цилиндре

 

 

 

 

 

Выполнил :

К-т СМФ

Группы 2161

 

 

 

 

 

Одесса – 2013

Содержание 

Введение …………………………………………………………………………..3

  1. Физические основы самовоспламенения и сгорания топлива..................................................……………………………………4

2.  Процесс сгорания ………………........................................………….…12

3. Математическая модель модель сгорания топлива в цилиндре дизеля........................................... ………………………………………….……16

Выводы ………………………………………………………………………..…23

Список литературы ………………………………………………………..…….24 

Введение

Актуальность  темы. На судах морского флота Украины в качестве главных установок в основном используются дизельные двигатели, эксплуатация которых сопровождается существенными расходами на топливные и смазочные материалы. Доля расходов на топливо и масло в общей сумме эксплуатационных расходов составляет около 40%, при этом стоимость топлив всё возрастает.

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

 

 

 

 

 

 

 

 

 

 

1.Физические основы  процессов самовоспламенения и  сгорания топлива

Сгорание топлива в  дизелях представляет собой процесс  окисления  молекул  углеводородного  топлива  кислородом  воздуха  с выделением теплоты.  Химическим процессам  окисления топлива  предшествуют  физические:  прогрев  капель  топлива  от воздуха  в  камере  сгорания  и  их  испарение.  Вследствие  малых размеров  капель  распыленного  топлива  (средний  диаметр  капель dk = 20-40 мкм) в факеле и высокой температуры воздуха в камере  сгорания  (более 500°С),  эти процессы  для каждой  образующейся  при впрыске капли топлива составляют очень короткий временной интервал  (менее 1  мс). Установлено, что процессу самовоспламенения топлива предшествуют также подготовительные  химические  процессы  -  распад молекул и образование промежуточных продуктов предпламенного окисления. Физико-химическая  подготовка  топлива к сгоранию  приводит  к запаздыванию  начала  выделения тепла в цилиндре  относительно  начала подачи топлива форсункой. На рисунке 1 показана индикаторная диаграмма судового среднеоборотного  дизеля  (n  =  720  об/мин)  в координатах р-φ. Здесь же  показаны  моменты начала  и конца подачи  топлива форсункой. Момент резкого возрастания давления  в конце процесса  сжатия,  отмеченный на рисунке точкой НВС, принято называть  началом видимого  сгорания топлива.  Этот момент  соответствует началу  интенсивного  выделения тепла вследствие начала  сгорания  топлива.  Угловой промежуток  от  начала  подачи топлива форсункой до  момента НВС,  обозначенный  на  рис. 1

Φi составляет отмеченное  выше  запаздывание  самовоспламенения топлива,  поэтому его называют  периодом  задержки  самовоспламенения топлива. Так как  φі   = 6nt1, то задержка самовос пламенения по времени определяется формулой

Для судовых дизелей значения  ti   составляют 2-15 мс.Период  задержки  самовоспламенения короче  для тех топлив, которые термически менее устойчивы. Склонность к распаду и к образованию легко воспламеняющихся смесей  зависит от структуры молекул углеводородов, составляющих топливо.Топлива нефтяного происхождения представляют  собой смесь большого числа углеводородов. С точки зрения способности к самовоспламенению в условиях дизеля их можно разделить на две группы:  1  -  углеводороды, имеющие цепную  структуру молекул (парафины);  2  -   углеводороды  с кольцеобразной структурой  соединения  атомов  углерода  в молекуле  (ароматики).  Структуры простейших  молекул углеводородов представлены на рис. 2

Рис. 1. Фазы процесса сгорания топлива НПФ,  КПФ -  начало и конец подачи топлива форсункой; НВС -  начало видимого сгорания топлива;  I -  задержка самовоспламенения;  II -  быстрое сгорание топлива;  III -  диффузионное сгорание и догорание топлива.

Углеводороды  первой  группы  обладают наименее прочной  связью углеродных  атомов,  поэтому  легко  самовоспламеняются в  дизеле.  Ароматические  углеводороды  (2  группа)  отличаются большой  прочностью  молекулы,  поэтому  их  способность  к  самовоспламенению мала.

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

Рис. 2 Структуры молекул углеводородов

В  качестве  эталона  наилучшей  воспламеняемости  принят углеводород  цетан  (C16H32),  ему присвоено ЦЧ=100;  в качестве эталона наихудшей воспламеняемости  принят  ароматический углеводород а-метилнафталин  (C7H10CH3), для которого ЦЧ=0.

Цетановое число смеси  этих двух углеводородов определяется  долей  цетана  в  смеси.  Очевидно,  чем  больше  цетановое число  топлива, тем оно лучше воспламеняется в дизеле и наоборот. Цетановые  числа легких маловязких дизельных  топлив  определяются  на  специальной  лабораторной  дизельной  установке (ИТ9-ЗД)  путем  сравнения  их  воспламеняемости  с  эталонными смесями цетана и а-метилнафталина. 

Значения  ЦЧ  дизельных  топлив  составляют  50-60.  Способность  к самовоспламенению у тяжелых  высоковязких топлив (мазутов)  ниже  (их  ЦЧ=30-35)  из-за  большого содержания в них ароматических углеводородов.  Для этих  топлив  оценка ЦЧ проводится  косвенными  методами,  так как из-за  их  высокой вязкости лабораторная установка не может использоваться.

Период  задержки  самовоспламенения  принято  считать первой  фазой  процесса  сгорания  топлива.  В  обобщенном  виде зависимость  задержки  самовоспламенения  топлива  в  дизеле

можно выразить следующим образом:

где  р ,Т   -  давление и температура воздуха в  камере  сгорания  в момент начала подачи топлива форсункой.

Экспериментально  установлено,  что  φi  нелинейно зависит от давления, температуры и ЦЧ. Увеличение давления, температуры,  ЦЧ  приводит  к сокращению φi и наоборот.  С частотой вращения  φi  связан  прямо пропорционально,  поэтому повышение оборотов двигателя приводит к его возрастанию.

На  величину  φi   влияют  и другие  факторы,  например,  конструкция и температура стенок камеры  сгорания, движение  воздушного  заряда  и др.  Этим  объясняется различие  величин φi  в двигателях  различного  конструктивного исполнения  даже  при одинаковых  значениях параметров р,Т ,ЦЧ ,n. В зависимости от сочетания перечисленных факторов  период  задержки  самовоспламенения в судовых дизелях составляет:  МОД -   3-6°п.к.в.; СОД – 5-15  °п.к.в.  В высокооборотных дизелях φi может достигать 20 и более град, п.к.в.

Вторая фаза процесса сгорания - фаза быстрого сгорания топлива -  начинается в точке НВС и  завершается в момент достижения  максимального  давления  газов  в  цилиндре  (точка  z  на рис. 1)  К  началу  сгорания  в  цилиндр  подана  часть  цикловой подачи  топлива  gці,  которая зависит от  величины  φі (при однофазном  впрыске связь между ними  примерно пропорциональная). Большая часть gці к моменту самовоспламенения подготовлена  к сгоранию:  испарилась,  пары топлива смешались с воздухом,  в топливовоздушной  смеси образовались химически активные части молекул -  радикалы.

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

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

Интенсивность  сгорания  во  второй  фазе  оценивают  средней  скоростью  нарастания  давления  при  сгорании — wcp ,  которая определяется  как отношение приращения  давления  от  момента НВС до  точки z  к продолжительности второй  фазы  в° п.к.в.,  а также степенью  повышения давления  газов в цилиндре λ=pz/pc, которые называют показателями динамичности рабочего  цикла.  Величину wcp определяют  по  индикаторной  диаграмме (см. рис. 1). Величина wcp  показывает,  на  сколько бар повышается  давление  при повороте  коленчатого вала  на  1  градус  в среднем за фазу быстрого  сгорания. Чем ее величина больше, тем «жестче»

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

Для МОД Wcp — 2 - 3  бар/°п.к.в для СОД wcp — 3 - 4  бар/°п.к.в. и для ВОД Wcp -  6 - 8  бар/°п.к.в.  При больших значениях wcp  работа двигателя становится жесткой,  появляются  стуки в цилиндре, увеличивается износ деталей.

Очевидно, что динамические показатели  связаны  с  величиной  φi чем он  больше,  тем больше  будет gці,   тем большая доля цикловой  подачи  сгорит  во  второй фазе, тем выше  будут значения  λ и wcp,  и наоборот.

Третья  фаза  процесса  сгорания  длится  от  момента  достижения pz  в цилиндре  до  полного выгорания топлива.  В отличие от точек НВС и z, момент окончания сгорания топлива по индикаторной диаграмме определить невозможно. Расчетными методами  с использованием  экспериментальных диаграмм  установ лено,  что в судовых МОД и СОД продолжительность третьей фазы составляет 30—40° п.к.в.  Подача  топлива,  как правило,  в начале  третьей фазы  процесса  сгорания  еще продолжается,  в камере  сгорания  располагаются горящие топливные факелы, достигшие своего максимального размера.

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

Значительная часть воздушного заряда в силу неравномерного  распределения  топлива  по  объему  камеры  сгорания  и  малого  времени,  отводимого  на  смесеобразование  и  сгорание,  не участвует  в  сгорании.  Температура воздуха вне топливных факелов существенно ниже (порядка  800-1000°С), чем в зонах горения. Таким образом, в реальном рабочем  процессе дизеля имеет место высокая степень неоднородности  концентраций топливовоздушной смеси и температур  по  камере  сгорания. В связи с этим температура газов и коэффициент избытка воздуха,  используемые при расчете рабочих процессов дизелей (в отличие от  давления газов), являются условными параметрами, усредненными по камере сгорания.

Характер  сгорания  капель  топлива,  не  успевших  сгореть  к началу  третьей фазы процесса, а также поступивших в камеру сгорания  в процессе  продолжающегося впрыска, существенно отличается от  сгорания  во  второй  фазе.  Так  как  в  горящем факеле очень высокие температуры -   в среднем по факелу 2500-2700°С, то процессы  прогрева и воспламенения паров  многократно сокращаются. Скорость сгорания топлива лимитируется  физическими процессами: испарением капли и взаимодиффузией паров топлива и воздуха вокруг капли. Принято считать, что сгорание топлива в 3  фазе определяется процессами турбулентной диффузии. В связи с  сокращением содержания воздуха в факеле и медленным его притоком в  зону горения диффузионное  сгорание к  концу фазы  переходит в  догорание топлива с низкими скоростями.

Сокращению процесса догорания  топлива способствует организованное  (вихревое) движение воздушного заряда в камере сгорания.  Продолжительность 3 фазы также зависит от продолжительности  впрыскивания топлива форсункой  φвпр= Фнпф  + фкпф Сокращение φвпр  приводит к сокращению длительности 3 фазы и в целом процесса  сгорания  топлива и наоборот.

  1. Процесс сгорания

В процессе сгорания топлива  происходит выделение теплоты, преобразуемой  в двигателе в полезную механическую работу. Самовоспламенение и  сгорание топлива  в дизеле отличаются чрезвычайно  сложными физико-химическими процессами, точное математическое описание которых  не получено до настоящего  времени. Даже  при современных методах  расчета  с  помощью ЭВМ сгорание описывается полуэмпирическими  уравнениями. В  начале XX века профессор  МВТУ В.И. Гриневецкий предложил  метод расчета сгорания, который  был впоследствии развит его учеником Е.К. Мазингом и с некоторыми усовершенствованиями  используется  до  настоящего  времени.  Это  объясняется  исключительной  простотой метода при  вполне  приемлемой точности определения показателей  двигателя.Суть  метода  заключается  втом,  что  реальное  изменениедавления  в  процессе  сгорания заменяется  условно  эквивалентным  подводом  тепла  по изохоре  с — z’  и  изобаре z’ — z ,  как показано на рисунке  4.  В действительном процессе  участок  сгорания  топлива  характеризуется  плавной кривой  с - m - у .

Максимальное  давление  сгорания  достигается  в  точке  m после ВМТ поршня  и является  функцией  процесса  сгорания. Завершается процесс сгорания  в начале  процесса  расширения  в точке у. Предложенная  схематизация  соответствует  идеальному термодинамическому циклу  со  смешанным подводом тепла,  однако  в  предложенном  методе,  в  отличие  от  идеального  цикла, учитывается  изменение  массы  рабочего  тела  вследствие  сгорания  топлива,  изменение  его  термодинамических  параметров  в зависимости  от температуры и состава, теплообмен  со стенками цилиндра.

В  начале  процесса  сгорания  (точка  с)  известны  давление, температура, масса и  объем  pабочего  тела,  соответственно  -  рс, Тc Мс  и Vc. Напомним, что при отсутствии утечек рабочего тела в процессе  сжатия Мс = Ма.   Задачей расчета процесса  сгорания является  определение этих  же  параметров  в точках  z '  и z  . Промежуточные состояния рабочего  тела  в рассматриваемом методе не рассчитываются.

В точке z’ известен только объем, так как Vz  — Vc. В точке z при допущении о завершении химических реакций окисления топлива может быть  определена  масса  рабочего  тела Мz= Мс+∆ M , где ∆M, кмоль - приращение массы рабочего тела вследствие сгорания топлива. В данной постановке задача  не  разрешима из-за  слишком большого  числа неизвестных параметров.  В.И.  Гриневецкий предложил считать известной величину  максимального давления  сгорания,  принимая  его при проектировании  нового  двигателя на  основании данных  двигателя -  прототипа.  Это позволяет определить  координаты  точки z   на p -V   диаграмме и уменьшить число неизвестных термодинамических параметров в точке z  до двух -  Тz   и Vz.

Процесс сгорания рассчитывается исходя из  1  кг сжигаемого  топлива,  массы  компонентов  рабочего  тела  выражены  в  киломолях.  Температура  рабочего  тела  в  точке  z  определяется  из уравнения первого закона  термодинамики,  которое применительно к процессу  сгорания  может быть  выражено  в конечных разностях в виде

где  ∆Ucz=Uz-Uc   -  приращение  внутренней  энергии рабочего тела на участке сгорания ∆Lz’zвнешняя работа,  совершаемая  рабочим  телом  на  изобарном  участке  процесса  сгорания;

∆Q  -  количество теплоты, расходуемое при  сгорании  1  кг топлива  на  повышение  внутренней  энергии  рабочего  тела  и  совершение механической работы на участке  с - z’-z .

При  полном  сгорании  1  кг  топлива  выделится  количество тепла,  равное  QH.   В pассматриваемом  методе  расчета  сгорания принимается,  что  ∆QCZ  = QH-  Qhc  -  Qохл ,  где Qhc-   тепло,  не выделившееся  в точке z вследствие того,  что не  все топливо успело сгореть к этому моменту;  Qохл  -  тепло, отведенное от рабочего  тела  в стенки  камеры  сгорания  на  участке сгорания c-z. Отношение ∆Qcz к QH  называют коэффициентом использования тепла в точке z.

Физический смысл этого коэффициента заключается  в том, что  его  величина  (всегда <1)  отражает  потери  тепла  вследствие неполноты  сгорания топлива и теплообмена. Согласно опытным данным  его  значения  для  малооборотных  и  среднеоборотных дизелей составляют 0,75-0,9; для высокооборотных -  0,7-0,85.С  учетом  отмеченного  выше  ∆QCZ= ĘZQH.   Дальнейшиепреобразования  уравнения сводятся  к подстановке в него развернутых выражений для внутренней  энергии и работы  применительно к начальной и конечной  точкам  процесса  (с и z). Выражения для расчета внутренних  энергий приведены в таблице  1.

Наименование

Точка с 

Точка z

Срстав смеси газов

Свежий заряд воздуха  Мв и остаточные газы Мг

Продукты сгорания топлива М и остаточные газы Мг

Масса смеси, кмоль

Мсвг

Мz=М+Мг

Температура, К

Тс

Tz

Средняя молярная изохорная  теплоемкость кдж/(кмоль*К)

Пренебрегая влиянием остаточных газов сvc=20+0.0024Tc

Выражаем по формулам c’vc=20+0.0024Tz; c‘‘vz=21.5+0.0035Tz;

 

Расчетная формула для  внутреней энергии, кдж

Uc=Mc*cvc*Tc

Uz=Mz*cvz*Tz


Формула   может  быть  записана  в  виде ĘzQh=Uz-Uc+pzVz-pzVc . После подстановки в нее формул для внутренних  энергий и масс рабочего  тела в начальной и конечной  точках  процесса  сгорания  и несложных алгебраических  преобразований  получим уравнение  сгорания  топлива  в окончательном виде

 

Где λ=pz/pc

В  результате  подстановки  в  левую  часть  уравнения уже  известных  из  расчетов  процессов  наполнения  и  сжатия  и принятых из рекомендованных пределов значений параметров, в каждом  конкретном  случае  расчета  будет  получено  некоторое числовое  значение  С.  В  правой  части  уравнения  теплоемкость смеси  газов является функцией  искомой  температуры  Tz, поэтому в итоге уравнение сгорания  приводится  к виду ATZ2  + BTZ — С = 0,  где А,  В и С - числа.  Уравнение можно привести к виду и  решать  его методом  последовательных  приближений,  задавая в  качестве  первого  приближения  значение  температуры  в  знаменателе формулы из диапазона  1700-1900 К.

Из формулы  очевидно, что  с увеличением  λ, Ęzс  и уменьшением а   температура в конце сгорания увеличивается и наоборот. В зависимости от значений этих величин для судовых дизелей характерен диапазон  Tz   = 1700 -1900 К. Как уже отмечалось ранее, значение максимального давления  сгорания должно быть принято по данным двигателя-прототипа.

Для  определения  объема  в  точке  z  запишем  уравнения  состояния рабочего для начальной  и конечной точек процесса сгорания  pzVz = RMzTz   и pcVc = RMcTc. Разделив  почленно первое на второе, получим

Определив  из  этого уравнения  значение  степени предварительного  расширения  рабочего  тела ρ,  получим искомое значение  Vz =  Vcρ. Таким образом, все термодинамические параметры состояния рабочего тела в точке z определены.

3. Математическая  модель сгорания топлива в  цилиндре дизеля

Первое  слагаемое  в  правой  части  уравнения  QH=(dgx/dφ)  представляет  собой скорость  выделения тепла при сгорании  топлива,  кДж/(°п.к.в.)  В свою  очередь,  dgx/dφ  называют  скоростью  сгорания  топлива,  кг/(°п.к.в.).  Текущую  массу сгоревшего  топлива gx  можно представить как gx = gцx,  где gц - цикловая  подача  топлива,  кг/цикл;  х =  gx/gц   -   относительное

количество  сгоревшего  топлива.  Очевидно,  что х изменяется  от нуля в начале процесса сгорания до  1  в его конце  (полагаем, что топливо  сгорает  полностью).  Первая  производная  от  х  по  углу поворота коленчатого вала dx/dφ называется относительной скоростью  сгорания топлива и имеет размерность  1/(°п.к.в.).

Зависимости  dx/dφ = fi(φ)  и х=f2(φ) часто называют законом сгорания топлива, записанным в дифференциальной и интегральной формах  соответственно. Сгорание  топлива в значительной мере зависит от «закона»  его подачи  в цилиндр gт=f(φ), где gт -  масса впрыснутого в цилиндр топлива. По аналогии со сгоранием топлива введем понятие  дифференциальной  и  интегральной характеристик впрыска топлива      dơ/dφ = f3(φ)  и ơ = f4(φ).

Здесь ơ = gт/gц -  относительная масса впрыснутого топлива. Характеристики впрыска определяются, главным образом,  формой кривой  давления  топлива перед распылителем форсунки. Протекание  процессов впрыска и сгорания  топлива существенно зависят от закона подачи топлива.

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

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

В  инженерной  и  исследовательской  практике  чаще  всего используются наиболее простые модели сгорания, позволяющие экономить время  счета ЭВМ. Хорошее  соответствие  опытных  и расчетных  данных  достигается  путем  выбора коэффициентов в уравнениях,  описывающих характеристики  сгорания.  Большое распространение получила формула И.И. Вибе, которая имеет вид

где  φz -  условная продолжительность  сгорания топлива; m - показатель характера процесса сгорания.

Используются также формулы  других авторов, именуемые «эмпирическими».  Большим недостатком этих формул является то,  что они не  учитывают характеристику  впрыска топлива. Кроме того,  эмпирические  формулы носят частный характер  и не  могут применяться при моделировании рабочих процессов дизелей других типов без предварительных исследований.

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

В  этом методе между  ơ  и х  введена промежуточная  функция w =f5(φ), отражающая динамику подготовки  топлива к сгоранию. Параметр w= gw/gц  представляет собой текущее относительное количество подготовленного к сгоранию топлива. Полагаем,  что подготовленное  к сгоранию  топливо связывает воздух в  объеме  факела  стехиометрическим  соотношением.  Положив также,  что  скорость  протекания  физических  процессов  пропорциональна  общей  текущей  поверхности  капель  распыленного топлива,  запишем  уравнение  для  относительной  скорости  подготовки топлива к сгоранию в следующем виде:

 

Где величина,  пропорциональная  суммарной поверхности капель  топлива в камере  сгорания;  у = f6(φ)-корректирующая функция.

В  уравнении  использованы  также  следующие  обозначения:  х-эмпирический показатель, учитывающий  степень неравномерности  распыливания топлива;  при  абсолютно равномерном  распыливании  (все  капли  имеют  одинаковый  диаметр) X = 2/3, при неравномерном  распыливании  х  = 0.8-1,0.

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

 

 

структурно формула состоит из трех комплексов. Первый из них

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

и  отражает  снижение  скорости  подготовки  топлива  к  сгоранию по мере роста  концентрации  подготовленного  топлива и  снижение  массы  «чистого»  воздуха  в  среднем  по  объему  камеры  сгорания.

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

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

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

Схема  расчета  характеристик  сгорания  иллюстрируется рис.  5 Момент  начала  подачи  топлива  φнпф  и характеристика впрыска должны  быть  заданы  в исходных  данных.  Относительная скорость подготовки к сгоранию рассчитывается по уравнению с учетом  зависимостей  на протяжении всех  периодов  процесса  сгорания.  Расчет  dw/dφ ведется совместно с уравнениями ,  отражающими изменение параметров рабочего тела в цилиндре. Текущее количество подготовленного к сгоранию  топлива определяется  интегрированием уравнения на каждом шаге счета.

Период  задержки  самовоспламенения  может  быть  рассчитан  по  одной  из  эмпирических  формул.  Для  судовых  дизелей используется формула

Построение математической модели сгорания топлива в цилиндре