Моделирование процессов в горючих спреях Текст научной статьи по специальности «Математика»

УДК 517.9+ 621.436

МОДЕЛИРОВАНИЕ ПРОЦЕССОВ В ГОРЮЧИХ СПРЕЯХ

© 2011 С. С. Сажин1, В. А. Соболев2, Е. А. Щепакина2

'Университет Брайтона, Великобритания 2Самарский государственный аэрокосмический университет имени академика С. П. Королева (национальный исследовательский университет)

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

Инвариантные многообразия, негладкие нелинейности, редукция, дизельный двигатель, горючий

спрей.

Модель воспламенения в дизельном двигателе

Процесс воспламенения рассматривается в предположении об одностадийной химической реакции согласно закону Аррениуса с постоянным предэкспонен-циальным фактором. Это допустимо, если концентрация горючего настолько мала, что изменение процентной концентрации расходуемого кислорода можно не учитывать. При таком подходе не принимаются во внимание некоторые важные эффекты, свойственные процессу самовоспламенения, такие как формирование холодных пламен, но в то же время можно адекватно описывать общий баланс энергии. Следуя работам [1, 2], в модель включаем уравнения баланса энергии, радиуса капли и концентрации горючего газа, которые в безразмерной форме имеют вид -в

7 (1 + Рв Г -Г = н„ (М) - НІ (М); (1)

ат

-г3

-т

-е 2 Нь (в ,ц);

~Т = ~ Н (в ,ц) + у Н (в ,ц), ат

(2)

(3)

где Нк (в,ц) = ц ехр Н (в) = Єї г [в (1 + Рв)

в

1 + рв

1/2

+ Є3 Г

((1 + рв )4 -1)1

в - безразмерная температура газа; ц -безразмерная концентрация горючего газа; г - безразмерный радиус капель; р, 7

- традиционные параметры в теории теплового взрыва Семёнова; Р - приведённая начальная температура по отношению к так называемой температуре активации; 7 представляет конечную безразмерную адиабатическую температуру термически изолированной системы после взрыва. Характерные значения этих параметров малы по отношению к единице для типичных газовых смесей благодаря высокой экзотермичности химической реакции и высокой энергии активации [1]. Параметры е1, е2 характеризуют взаимодействие между газовой и жидкой фазами, а параметр е3 характеризует влияние тепловой радиации и представляет собой отношение коэффициентов лучевой и конвективной теплопередачи. Характерные значения у лежат в диапазоне 10-100 и выше [1].

Начальные условия для уравнений (1)-(3) имеют вид

в(0) = в0, т = % г (0) = 1. (4)

Интеграл энергии для уравнений (1)-(3) может быть представлен следующим образом:

Л = ^0 - Ь ln

/ 1 + рв Л 1 + Рв0

+

У

e 2

- (1 - г3).

(5)

/ exp р У -1 ho + \ + ln(1 + рв0) \ - 1

V _ e 2 _ ) 0

Уравнения (1)-(3) позволяют найти условия равновесия. Естественно ожидать, что завершение процессов испарения и горения приводит к равенствам г = 0,h = 0. Это позволило найти выражение для температуры в конце процессов:

ввfinal ь

V

Редукция системы

Метод интегральных многообразий обычно применяется для исследования систем дифференциальных уравнений с достаточно гладкими правыми частями, хотя есть примеры успешного применения этого метода даже для редукции (понижения размерности) систем с разрывными правыми частями [3]. В данной работе обсуждается вопрос о возможности редукции систем с негладкими, более того, с нелипшицевыми нелинейностями. Напомним, что векторная функция F(z) называется липшицевой или, более точно, удовлетворяющей условию Липшица с константой L в некоторой области изменения векторного аргумента z,

если имеет место неравенство | F(Zj) - F(z2)| < L | Zj - z2 | для любых Zj z2 из этой области.

В последние годы появились статьи, посвященные исследованию моделей химической кинетики с нелипшицевыми нелинейностями (см., например, [1]). К сожалению, в этих работах сначала излагаются общеизвестные результаты теории интегральных многообразий сингулярно возмущённых систем с гладкими правыми частями, а затем рассматриваются физические задачи, к которым эти результаты

неприменимы. Следует отметить, что можно ввести понятие положительно (отрицательно) инвариантного множества, т.е. множества, обладающего следующим свойством: если точка траектории принадлежит этому множеству в некоторый момент времени =о, то соответствующая положительная (I > (отрицательная

(I < 1йУ) полутраектория принадлежит этому множеству целиком. Если множество является и положительно и отрицательно инвариантным, то оно инвариантно. Для систем с липшицевыми правыми частями эти три понятия эквивалентны. Ясно, что притягивающие положительно инвариантные множества могут использоваться для понижения размерности моделей, но обоснование такой возможности выходит за рамки теории метода интегральных многообразий и требует применения другого математического аппарата.

Введение новых переменных

3

уг

+ Ц при-

e2

p = ln(1 + рв), q = г3, водит систему (1)-(3) к виду

b^ = G(s, p, q) - e, q"1 F(p, q); р dt

dq

dt

—г^=- e1 q F(p, q);

d- =- G(s, p, q), dt

(6)

(7)

(8)

где G(s, p, q)

F (p, q)

1

/ \ У

s-------q

e2

2

\

exp

ep/2 + e 3 q1/3

ep -1

pep

к у

(e4 p -1).

V ' /

Типичные значения

параметров в случае дизельного двигателя подчиняются неравенствам 1/ е2 << у / /3 < 1 (для значений параметров, которые будут обсуждаться ниже у/3 = 0,268, 1/е2 = 0,0058 при начальной температуре на поверхности капель 600 К и 1/е2 = 0,0084 при 300 К). Это даёт возможность различать характерные значения величин

1/ е2, у / 3 в системе (6)-(8), вводя ие-

2

p

e

рархию переменных: переменная д является быстрой, переменная р - более медленная и переменная 5 - самая медленная. Таким образом, система (6)-(8) может рассматриваться как сингулярно возмущённая с малым параметром V = 1/ е2 .

Полагая V = 0, из уравнения (7) для быстрой переменной получаем уравнение медленной поверхности 0 = - е1 д1/3 Р(р, q) . Заметим, что функция

F(p,q) принимает положительные значения, следовательно, данное уравнение даёт единственное решение q=0, описывающее медленную поверхность системы (6)-(8). Плоскость q=0 представляет собой притягивающее положительно инвариантное множество дифференциальной системы (6)-(8).

Нетрудно заметить, что функция Р(р,д) в окрестности нуля не удовлетворяет условию Липшица по переменной д. Нарушение этого условия ведёт к тому, что не выполняется теорема о существовании интегрального многообразия [4, 5] и сам метод нахождения многообразия в виде разложения по степеням V неприменим. Именно на данном этапе наблюдаются разногласия в теории и её применении в статье [1].

Для обоснования допустимости перехода от полной системы (1)-(3) к системе меньшей размерности на плоскости д=0 покажем, что эта плоскость является притягивающим положительно инвариантным множеством. Для этого достаточно проверить, что положение равновесия д=0 так называемого присоединённого уравнения

^ =-Лдш -Вд2/3

асимптотически устойчиво [4, 5]. Для проверки рассмотрим функцию Ляпунова —(д)=1/2 ц2, производная которой в силу присоединённого уравнения

-— = - Л д4/3 - Вд53, (д > 0) отрица-

ds

тельно определена, что и означает выпол-

нение требования асимптотической устойчивости.

Можно выделить два этапа в динамике химической системы. Первый отвечает процессу испарения и траектории системы (6)-(8), соответствующие этому этапу химической реакции, за очень короткий отрезок времени выходят на плоскость д=0. При этом значения переменной д изменяются очень быстро на фоне почти неизменных значений р и 5. После окончания процесса испарения (в данном случае при д=0, или в исходных переменных при г=0) порядок системы (1)-(3) может быть понижен и можно переходить к исследованию зависимости между переменными в и Ц .

Испарение капель

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

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

считается, что теплообмен между пламенем и каплей осуществляется за счёт конвекции.

Используя замену переменной т =пт в (6)-(8), получаем следующую систему:

[Оp, ч) - е 1 диз р(р, ч)];

у ёр -—= V

Р ёт

= - е 1 чш р(p, ч);

ёт

^ =- V О(s, p, ч). ат

Полагая малый параметр V равным нулю, из последней системы имеем ёр &

= ^=- = 0 и, в соответствии с началь-ат ат

ными условиями (4), получаем р ° 1п(1 + рв0); ^ = пу + Ч о.

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

V 44 =-Лч1,ъ - В ч2/3 ат

где Л = е1 в0 (1 + рв0 У

(9)

41/2

В=ее ((1+рв0)4-і); ф)=1.

Уравнение (9) можно представить в

виде ёг3

V

ёт

= - Аг - В г

Непосредственное интегрирование этого уравнения с учетом (4) даёт

3

е 2 В

А А + В

1 - г (т)-----1п-----------

В А + Вг(т)

3

' еуар

Є 2 В

А+В

(10)

системы, можно изменять параметр теуар,

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

Воспламенение газовой смеси

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

Для значений времени т >теуар безразмерный радиус капель равен нулю и система (1)-(3) принимает вид

в

у(1 + рв)1 — = л ехР ёт

1 + рв

ёц

ёт

= - л ехр

в

1 + рв

ч х

Исключая переменную т из этой системы, приходим к уравнению

у(1 + рв)'

-1 ёв

Исходя из полученной взаимосвязи между безразмерными временем и радиусом капель, можно вычислить время испарения т =теуар , полагая Г(теуар ) = 0 :

ёч

Непосредственное интегрирование приводит к соотношению р

в (т )

(1 + рв0 ) ехр Гр (ц о- л) ^ -1

_ 1 у 0 _

Этот параметр играет важную роль. Имея аналитическое выражение для его подсчёта и варьируя начальные данные

где в 0 = в (теуар X Ч 0 = Ч (теуар ). Следует отметить, что зависимость между в и ч при т >теуар может быть получена из (5) при г=0. Эта зависимость показывает, как бу-

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

Работа выполнена при поддержке гранта РФФИ 10-08-00154 а.

Библиографический список

1. Sazhin, S.S. Thermal ignition analysis of a monodisperse spray with radiation [Text] / Sazhin S.S., Feng G., Heikal M.R. and others // Combustion and Flame. - 2001.

- V. 124. - Issue 4. - P. 684-701.

A NEW MATHEMATICAL TOOL FOR MODELLING THE PROCESSES IN FUEL SPRAYS

© 2011 S.S. Sazhin1, V.A. Sobolev2, E.A. Shchepakina2

'University of Brighton, UK 2Samara State Aerospace University named after academician S. P. Korolyov

(National Research University)

New concepts of positively (negatively) invariant manifolds for non-linear systems with non-Lipschitzian nonlinearities have been introduced. It is pointed out that the order reduction of these systems can be performed, using these new concepts. It is shown that the singularly perturbed system of three equations (droplet radii, gas temperature and fuel vapour concentration) describing evaporation and ignition of monodispersed sprays possesses the attractive positively invariant manifold. This allows us to reduce this system to the system of two equations describing gas temperature and fuel concentration when the droplet evaporation process has been completed. An additional equation describes the reduction of droplet radii when the variations of gas temperature and fuel vapour concentration can be considered small.

Invariant manifolds, non-smooth nonlinearities, reduction, diesel engine, fuel spray.

Информация об авторах

Сажин Сергей Степанович, кандидат физико-математических наук, профессор Школы информатики, техники и математики, Университет Брайтона, Великобритания. Область научных интересов: моделирование процессов тепло- и массопередач в спреях, моделирование двигателей внутреннего сгорания. E-mail: S.Sazhin@brighton.ac.uk.

2. Sazhin, S.S. Order reduction of a non-Lipschitzian model of monodisperse spray ignition [Text] / Sazhin S.S., Shchepakina E., Sobolev V. // Mathematical and Computer Modelling. - 2010. - V. 52 -Issues 3-4. - P.529-537.

3. Соболев, В. А. Декомпозиция разнотемповых систем с разрывными управлениями [Текст] / В.А. Соболев, Л.М. Фридман // Автоматика и телемеханика. -1988. - № 3. - С. 29-34.

4. Singular perturbations and hysteresis [Text] / Eds. M.P. Mortell, R.E. O'Malley, A. Pokrovskii [et al.] - Philadelphia: SIAM, 2005. - 344 p.

5. Соболев, В. А. Редукция моделей и критические явления в макрокинетике [Текст] / В. А. Соболев, Е.А. Щепакина. -М.: Физматлит, 2010. - 320 с.

Соболев Владимир Андреевич, доктор физико-математических наук, профессор кафедры технической кибернетики, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). Область научных интересов: динамика многотемповых управляемых систем, математическое моделирование. E-mail: hsablem@yahoo.com.

Щепакина Елена Анатольевна, доктор физико-математических наук, профессор кафедры технической кибернетики, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). Область научных интересов: динамические системы, математическое моделирование. E-mail: shchepakina@yahoo.com.

Sazhin Sergei, candidate of physical and mathematical sciences, professor of School of Computing, Engineering and Mathematics, University of Brighton, UK. Area of research: mathematical modelling of fluid dynamics, heat transfer and combustion processes in internal combustion engines. Е-mail: S.Sazhin@brighton.ac.uk.

Sobolev Vladimir Andreevich, doctor of physical and mathematical sciences, professor of the department of technical cybernetics, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Area of research: dynamics of multi-scale control systems, mathematical modeling. Е-mail: hsablem@yahoo.com.

Shchepakina Elena Anatolyevna, doctor of physical and mathematical sciences, professor of the department of technical cybernetics, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Area of research: dynamical systems, mathematical modeling. Е-mail: shchepakina@yahoo.com.