CC BY
16Вестник Сыктывкарского университета. Сер. 1. Вып. 2.1996
УДК 536.46
К АСИМПТОТИЧЕСКОЙ ТЕОРИИ ВОЛНЫ ГОРЕНИЯ В ГАЗАХ В. М. Холопов, С. И. Худяев
Особенностью асимптотики волны горения газов даже в предположении равенства единице числа Лыоиса (Lewis) является наличие второго погранслоя в зоне горения.
Введение. Построение асимптотики стационарной волны горения газов по большой энергии активации было начато в работе [1], где рассматривалась одностадийная реакция первого порядка. Дальнейшие работы показали существенную роль в конструкции асимптотики таких параметров, как порядок реакции п и число Лыоиса L. По поводу особого случая L = 0 (отсутствие диффузии реагента, характерное для конденсированных сред) можно обратиться к [2, 3] и приводимой там литературе. В работе [4] предпринята попытка рассмотреть произвольные п > 0 и L = 0(1). Однако, как было отмечено в [2], составные решения неявно содержат ограничение п ^ 1. В работе [3] при L = 0 такое ограничение удалось снять путем построения второго погранслоя в зоне горения. В настоящей работе показано, что при L — 1 второй погранслой в зоне горения оказывается необходимым элементом конструкции асимптотики для п > 1.
1. Основное уравнение. Внешнее решение. Задача о распространении стационарной волны горения в газах в случае одностадийной необратимой экзотермической реакции n-го порядка в предположении аррениусовской зависимости скорости реакции от температуры и в условиях подобия полей концентраций и температур описывается дифференциальным уравнением первого порядка
И:
© Холопов В. М., Худяев С. И., 1996.
связывающим безразмерный поток тепла р и безразмерную относительную температуру Безразмерные переменные и параметры определяются формулами
, ___ Г, - Т АГ _ Д7?с _ Т, - ?р _ т2с
где Т* = То + — - температура горения (температура полностью с
прореагировавшего вещества), То - начальная температура, <5 - тепловой эффект реакции, А - коэффициент теплопроводности, с -теплоемкость, Е - энергия активации, К - универсальная газовая постоянная, т - искомая массовая скорость горения, аррениусова
/ Е \ ■
функция /с(Т) = А:оехр^—^ту - задает температурную зависимость
скорости химической реакции, &о - предэкспоненциальный фактор. Кроме неизвестной функции р(1) в уравнении (1) искомой величиной является также безразмерная скорость горения Наличие лишнего граничного условия позволяет, в принципе, ее определить.
метр 7, имеющий смысл относительной ширины зоны реакции,- предполагается малым, что согласуется с представлениями о фронте пламени. Именно по этому параметру будет строиться асимптотика решения (1). Параметр /г. как видно из (2), ограничен О < ¡1 < 1. Еще один существенный парамер, входящий в уравнение
(1), - порядок реакции п. Для 0 ^ п ^ 1 асимптотика построена в
[2]. В этой же работе показано, что-при п > 1 предложенная конструкция не является асимптотикой решения уравнения (1).
В настоящей работе строится асимптотика решения (1). Как всегда, при рассмотрении стационарной задачи горения требуется оговорка. Задача (1), строго говоря, не имеет решения (см. особенность при £ = 1), Стационарная волна имеет смысл временной "промежуточной" асимптотики [5], которая, строго говоря, должна строиться одновременно с намеченной асимптотикой по 7. Однако (см. [5]) роль малого параметра в асимптотике по времени играет величина порядка ехр(—7-1), и при построении степенных асимптотик по у она не сказывается.
Правая часть уравнения (1) сосредоточена в малой окрестности точки < = 0 и экспоненциально мала при 7 —> 0 вне этой окрестности. Поэтому для построения асимптотики решения можно применить метод сращивания асимптотических разложений [6], построенных вне вышеупомянутой окрестности точки г = 0 (внешнее реше-
ние) и внутри нее (внутреннее решение), Устремляя 7 к нулю в (1) и учитывая условие р( 1) = 0, находим внешнее решение
p(t) = l~t. (3)
2. Первый погранслой. Для построения внутреннего реше-
t
ния перейдем в окрестности точки t = 0 к новой переменной г — —.
7
При этом уравнение (1) примет вид:
'^ + 7 = -ехр(-т-1—), Р(0) = 0. (4)
ат pw V 1 — /гуг /
Однако построение одного "погранслоя" в окрестности t = 0 относительно г = —, р = ;;(г) и сращивание его с внешним решением
т ;
позволяет получить правильную асимптотику лишь при п ^ 1 [2]. При п > 1 такое построение становится невозможным. Уже в этой простейшей модели горения при п > 1 необходимым элементом асимптотики оказывается второй погранслой в окрестности т — О относительно переменных
е = ? = а' =
в которых уравнение (4) переходит в уравнение '
«<о»=°- <6)
Из уравнения (6) видно, что асимптотику для и и q, а следовательно и для р, следует искать в виде разложений по степеням у вида l + ká, где Í, к - целые неотрицательные числа:
л, о i i 9 л
ш = + OJ]0J + W017 +W207 + ^li7 +Ш027 .+ •••>
9 = ío + íio7 + 9oi7a + 92072 + <7п71+а + 2о272а + • • •, ' (7) р -ро +Рю7 +Poi7a +Р2072 +Рн71+а +Р0272а + • • • •
Подставляя разложения (7) в уравнение (4), для определения коэффициентов разложения внутреннего решения р(т) получим уравнения:
dpo т"е 7
dr uQpQ ?
Ро(0) = 0, (8)
&г+1„_ф-(а»+га»+/,4 „„«»-о. (?)
ат ат V и'о ро '
^ = + (0) = 0. (10)
ат ат \ (¿о ро /
Интегрируя эти уравнения, находим
г
/ 2 Г \ 5
/ "" I »> — с* \ *
РоИ= (-■/ „ (11)
О
г г
1 / / X / ^10 , х И
РЮ = .... ....... ~ "
Ро./ ^о" '^'(Шо
о о
г
РоШ* ~ ^-Мт) - -г'-- /' 2<ГЖ (12)
Ан = "-^-Ро(т-)- (13)
В дальнейшем, для сращивания первого внутреннего решения со вторым, нам потребуется асимптотика полученных выше коэффициентов ро,р\о и роь при г —* 0. Опуская несложные вычисления, выпишем результат.
2 н±1 / п+1
(,,+ 1)(п* + 4« + 5)тг+ ь (и)
8(гг 4- 2) (п + 3)
___2_ / п + 1
Рхй ~ п + 3Г V + 2(п + 2)(п + 5)Т+
{п + 1){п2 - 2п - 19) 2 . 2(п + 2)2(?г + 3)(п + 5)(п 4- 7)Т + ''
Шо / 2 тЧ1/ п + 1 , (п + 1)(п2 + 4п + 5) о у 2ш0]1 и0{п + 1) V 2(п + 2) 8(п + 2)2(п + 3)
2 V 2 п + 2 п + 4 к ;
2(п + 3) V 2 V 2(п + 2)(п + 4)
_ ^'ш / 2 п|±Л _ п +1 ~ 2(^о у ^о(гг + 1)Г V 2 (п + 2)Г+
(„+!)(»'+ 4Я + 5) т2 у
8(п + 2) (п + 3) / [ '
3. Второй погранслой. Для построения второго внутреннего решения подставим разложения (7) в уравнение (б). Раскладывая правую часть уравнения по 7 и приравнивая коэффициенты при одинаковых степенях 7, получим уравнения для определения коэффициентов goj<7i(h9oi разложения второго внутреннего решения
■ + ! = 9о(0) =0, (17)
Ч u0qQ
^ + = + 9lo(0) = 0, (18) «ч <7о 4 «ч у
^о-тг + — (701 = -<*>10 [-77 + 1)--, ?oi(0) = 0, .(19)
«£ Qo 4 «ч у tfo
Уравнение (17) не интегрируется в квадратурах, и для его решения необходимо использовать приближенные методы, а уравнения (18) и (19) являются линейными, и их интегрирование при известной </о не составляет труда.
Для согласования с первым внутренним решением р(т) необходимо найти асимптотику <7о, ?кь tfoi при £ —> оо. Подставив с/р = С£к
« + 1 , 1 I 2
в уравнение (17), найдем к — —-— = 1 Н—, С = *—:--г, т.е.
2 а у аi0{n + 1)
/ ^
главный член асимптотики до имеет вид Л —т-^—— . Дальней-
У u,'0(n + 1)
шие члены можно найти методом неопределенных коэффициентов q0 = С £1+(1 +• + aoi £ ~ " + ~2 + ап 1 ~ " + " + - • •)• (20)
Подставляя (20) в (17) и определяя коэффициенты Ою, «20,
«п, йог, получим асимптотику qü при £ оо
Аналогично находим асимптотики д10 и <y0i при £ —» оо.
о с* + 1 ., i / cvfa + l)2 „2 \
910 = f(i - 2¿(2a+\fC° +...), (22)
а + 1- I / . а
.2
= ^'¿Т^Ч1 + С(2а + 1)(За+1)е
+ ...)■ (23)
2а /
Здесь выписаны лишь те члены асимптотики', которые будут участвовать в сращивании.
Подставляя во вторую из формул (5) разложение для д (7) и возвращаясь к переменной г, получим функцию
(24)
Заменим здесь <7о,5нь(/01 на асимптЬтики-(21)-(23) и перегруппируем слагаемые, расположив их по соответствующим степеням у.
За+ 2 У '\2а + 1 V (За+ 1)(За+2)
.3» + 1 1+±\ , »„з . а + 1 1 + 1
2 а •/. 2а
Эта формула должна быть согласована с разложением (7) для р(г), коэффициенты которого заменяются их асимптотиками (14)-(16):
, . , 2 п±1Л п + 1
(п+ 1)(п2+4п + 5) 2
-7" +
8(п + 2)2(п + 3)
/ 2 /
•Пп + зЧ + (п + 2)(п + 5)Г+
(п + 1)(п2 - 2п - 19)
2(п + 2)2(п + 3)(п + 5)(п + 7)
^ П + 1 у
2щ]1 шо(п 1) Л 2(п + 2) т У
Ц + . \ _
2(п + 3) V 2 Г 2 7 230
Г 2 П + 1 Г + ..Л. (26)
- 7 2^о\/^о(" + 1) V 2(п + 2) / 1
/ 2 2
Вспоминая, что С == < /—;——г, а == -, видим, что в разло-
V о;о(п +1) п — 1
жениях (25) и (26) у коэффициентов при 7°,71,7и старшие члены
совпадают. Именно в этом смысле и понимается согласование, т.е.
возможность сращивания.
Асимптотика составного решения уравнения (4) будет иметь вид:
■ • 4 о ' .
т т
-7(- - ^РоМ - ^ /У'«^-
^о 7 2ро У
о о
2 / п+1 л I 2 '
п + 3Г V + (п + 2)(п + 5)Т) 2и>о у ио(п + 1)Т )
+0(72 + 71+а+72а). (27)
3. Асимптотика скорости. Поскольку первое и второе внутренние решения согласуются при любых с^о^ю и , то эти величины могут быть найдены из условия сращивания внешнего решения (3) с внутренним решением
р(т) ^р4г) ^УРФ) +'уароЛт), (28)
где Ро,Рю,Ро\ определяются формулами (Г1)-(13). Условие сращивания имеет вид р(0) =.р(оо), что для коэффициентов ра,р\0,р01 дает
Ро - 1,' рю ~ -т, Р1к = 0, {кфО) при т-* оо. (29)
Из формул (11)-(13) с учетом (29) находим коэффициенты разложения безразмерной скорости и :
шо = 2Г(п + 1), (30)
+оо
' w,0 = 2(o;0 j ( 1 " /'Г(л + 3)), (31)
о
wbi = 0,
+00
где Г(п) ~ ] $п~1е~*й8 - гамма-функция Эйлера. Так как Ш(ц = 0,
о
то из-(,13.) имеем р01 = 0.
Чтобы отметить особенность асимптотики скорости горения, выпишем уравнения для коэффициентов ри и ро2
■dpil _ Фо /2^K)u;oi - ип(j0 ^ 2рюро1 - pnpo
dr d.T V oîq . pi
Н---ь—-Н--/¿Г Н--~j.tr ], ^ц(О) = 0, (33)
ЩРо Л()Ро ^'о Ро '
Фо2 фо /^01 - "'02^0 , р1\- ртрц . .^'о\Ро\\ /т п
--7— =--т~ -5--Ь --—о-+ +'-Ь Р 1Ни) = и,
dт с1т \ ^ ргй ~ОРО у
интегрируя которые и учитывая, что ил — 0 и ро1 = 0, получим
Ри = (35),
¿и0
Ро-2 = ~~Ро{т). (36).
¿и о
Откуда, с учетом условий сращивания (29), получаем сип =: 0> Ри = 0, (¿02 — 0, Р02 — 0. Аналогично, нетрудно показать, что вообще и>1к ■— 0? Р1к = 0 при к ф 0. Поэтому асимптотика скорости будет содержать лишь целые степени 7.
Итак, асимптотика безразмерной скорости и имеет вид
+оо
ш == 2Г(п + 1) + 27(2Г(п + 1) J(l~p0(s))ds~Mr(n + 3))+O(72), (37)
о
где po(s) - определяется формулой (11). Используя последнюю из формул (2), получаем асимптотику массовой скорости горения
+OÛ ■ ' '
m2 = tllBEll(2Т(п + 1) + 2Т(2Г(п + 1) У (1 - po(s))ds-
-ЯГ(п + 3))+0(72)). Составное решение уравнения (1)
р(*) =?(})> (39)
где р задается формулой (27). _ '
Разложение (37) допускает предельный переход при п —> 1. С учетом численного значения интеграла имеем
и = 2(1 + 67(0,448 ... — р)) + 0(72). (40)
Поправка порядка 7 может принимать значения разных знаков в зависимости от 0 < /t < 1.
Асимптотика скорости (33), вообще говоря, ухудшается с ростом
п.
Литература
1. Bush W. В., Fendell F. Е. Asymptotic Analysis of Laminar Flame propagation for General Lewis N ambers// Comb. Sci. and Technology. 1970. V.1.P.421-428.
2. Худяев С. И. К асимптотической теории стационарной волны горения // Хим. физика. 1987.Т.6.№ 5.С.681-691.
3. Ильин А. М., Худяев С. И. Об асимптотике стационарной волны горения в конденсированной среде // Хим. физика. 1989. Т.8.№ 4.С. 525-532.
4. Берман В. С., Рязанцев Ю. С. К анализу задачи о тепловом распространении пламени методом сращиваемых асимптотических разложений // ПММ. 1972. Т.36.С.659-666.
5. Худяев С. И. К математической теории распространения пламени // Вестник Сыктывкар, ун-та. Сер.1. Вып.1.1995.С,234-243.
6. НаЙфе А.Х. Методы возмущений. М.: Мир, 1976. 455 с.
Summary
Kholopov V. M., Klmdyaev S. I, To Asymptotic Theory of Combustion Wave in Gases
Peculiarity of asymptotic of combustion wave in gases even if Lewis number equality to 1 is presented is the second boundary layer in zone of combustion.
Сыктывкарский университет Поступила 26.03.96
