Использованные обозначения
V -вектор скоростей, Ей - число Эйлера, Ке - число Рейнольдса, р - давление, \ - функция тока, в - сферическая координата, ф - потенциал, а, Ь, с - полуоси эллипса, п, С- криволинейные координаты, х, у, г - декартовы координаты, /л - динамическая вязкость, X', X, - силы сопротивления, - элемент поверхности.
СПИСОК ЛИТЕРАТУРЫ
1. Ламб Г. Гидродинамика. - М.: Гостехиздат, 1947. -928 с.
2. Дудин И.В., Нариманов Р.К. Сопротивление медленно движущегося в вязкой жидкости трехосного эллипсоида // Препринт № 37. - Томск: Изд-во ТГУ, 2000. -11 с.
3. Хаппель Д., Бреннер Г. Гидродинамика при малых числах Рейнольдса. - М.: Мир, 1976. - 630 с.
УДК 536.46
К ВОПРОСУ НАХОЖДЕНИЯ ПОСТОЯННОЙ МАРКШТЕЙНА
К.О. Сабденов
Томский политехнический университет E-mail: [email protected]
На примере простой модели показано, что уравнения горения допускают решения со смещенными профилями температуры и участвующей в реакции горения концентрации вещества. Скорость искривленного пламени определяется однозначно, а длина Маркштейна оказывается существенно больше, чем в ранних теориях. Это позволяет, по меньшей мере, на порядок расширить область гидродинамической устойчивости по числу Рейнольдса.
Пытаясь объяснить, почему ламинарное пламя в газе устойчиво к гидродинамическим возмущениям Дж. Маркштейн предположил, что скорость пламени ип зависит от кривизны К его фронта [1]. Когда фронт обращен выпуклой стороной к горючей смеси ип уменьшается, а если вогнутой - увеличивается. Это означает, что изогнутый под действием возмущений фронт пламени стремится принять плоскую форму. В предложенной формуле Маркштейна
ип = и0(1 - ¡мК), (1)
где ип0 - скорость пламени с плоским фронтом, присутствует неопределенная величина ¡м размерности длины (постоянная Маркштейна). Им же были сделаны первые попытки экспериментального определения этой величины. Дальнейшие теоретические поиски [2-5 и др.] приводили к выражению постоянной Маркштейна в виде отношения коэффициента температуропроводности газа к на скорость иП
к (2)
I =к
К'
Но в этом случае объяснение гидродинамической устойчивости пламени ограничено. Эксперименты показывают устойчивость пламени при числах Рейнольдса Ке, по меньшей мере, на порядок превышающих критическое значение КеС1, которое следует из формул (1) и (2).
Кроме указанного Маркштейном механизма, приводящего к устойчивости пламени, может иметь место случай ухода возмущений из рассматриваемой области горения [3]. Дело в том, что в экспериментах пламя, как правило, занимает ограниченное пространство (пламя горелки Бунзена или в трубе, свечи). Эти виды пламени имеют касательную составляю-
щую скорости газа к поверхности горения, благодаря чему гидродинамические возмущения или затухают на стенках, или уходят в свободное пространство из области горения прежде, чем они успевают заметно вырасти [3]. Но такие явления наблюдается при достаточно больших числах Рейнольдса.
Механизм ухода возмущений не объясняет полностью наблюдающуюся гидродинамическую устойчивость пламени. Если основываться только лишь на таком механизме и на формулах (1, 2), то мы уже при относительно небольших числах Рей-нольдса (Ке -10...100) должны были видеть пламя с непрерывно колеблющимся фронтом. В реальности это не имеет место: колебание фронта, свидетельствующее о наличии гидродинамической неустойчивости, возникает внезапно и при больших числах Рейнольдса (Ке ~103).
Дальнейшие пути объяснения гидродинамической устойчивости пламени можно искать в самой теории ламинарного горения. Предложенные в работах [2-5] способы нахождения ¡м, несмотря на их математическую строгость расчетов, не исключают существования других подходов. Дело в том, что если в теории Зельдовича-Франк-Каменецкого одномерного пламени с плоским фронтом его скорость движения определяется единственным образом [5], то скорость распространения искривленного пламени без привлечения дополнительного физического принципа не является однозначно определенной [6, 7]. Это обстоятельство приводит к неожиданному результату при исследовании неодномерной диффузионно-тепловой устойчивости пламени [8]. Использование известной [5] схемы анализа на основе формулы (1) и без привлечения
дополнительных предположении относительно процессов в зоне химической реакции приводит к абсолютной неустойчивости пламени.
В работе [9] было высказано предположение, что при искривлении фронта пламени меняется его структура: происходит своеобразное рассогласование полей температуры и концентрации (выгорания). Имеется в виду сдвиг распределений температуры и реагирующего вещества относительно друг друга и по направлению движения плоского исходного пламени. Величина смещения считается пропорциональной обратной величине энергии активации.
Настоящая работа ставит целью дальнейшее развитие изложенных представлений об искривленном пламени на основе более точного, чем в [9] теоретического анализа процессов во фронте неодномерного пламени.
1. Нахождение скорости пламени при слабо искривленном фронте. Бесконечно тонкая зона химической реакции
Для описания процессов медленного горения, когда давление меняется слабо по сравнению с температурой и плотностью газа, воспользуемся моделью с одностадийной химической брутто-реакци-ей А^В [5]. Анализ задачи об искривленном пламени будем проводить в линейном приближении. Это означает, что величина деформации фронта и вызываемые деформацией возмущения малы, и можно пренебречь их квадратичными и более степенями.
Рассмотрим модель ламинарного пламени в следующей формулировке, согласно которой основное движение пламени происходит отрицательному по направлению координаты X [5]:
W = q(1 - b )exp
dT dT (д2 T д2T Л
+ v^-r = к\-т +-т I + —W(N,T),
dt' д x' I д x'2 д у'2 I cp
5-N- + vд£ = D(Щ + дЩ |_w(N,T),
дt' д x' I д x'2 д у 1 V ;
(l.l)
где Ж(М,Т) - скорость химической реакции; V - скорость газа; Т - температура; N - концентрация реагента А; Б - коэффициент диффузии; 0 - тепловой эффект реакции; С - время; у' - координата, поперечная оси X; ср - теплоемкость при постоянном давлении.
Примем, что химическая реакция имеет первый порядок:
W (N ,T) = k0N exp( _ RT),
(l.2)
где к0 - константа, зависящая от свойств исходной горючей смеси; Е - энергия активации; Я - универсальная газовая постоянная.
Выражения (1.1) и (1.2) в безразмерной форме имеют вид:
дв дв д2в д2в
-т— + w= —т +—T- + W (в,Ь),
д t д x д x2 д у2
дь дь т (д2ь д2ьЛ и..в
-т- + w— = Le\ —т +—т I + W(в,Ь), дt д x I д x2 д у2 I V ;
(1.3)
4 \i _
l
ßl м+в(1_м)
, q = ККе
где
в =
T _T T _T
-L К -L Л
: Cp (T_T )• T = T + QN - t
QN 0 b 0T c 0'
N о _N • N -
D T0 0 RT. x'v, y'v, V v Le = К ц = ß=—b- x= —- у = --1 = — -wv. к Tb E к к к v.
Символом N0 обозначена начальная концентрация реагирующего вещества. Здесь v, _ масштаб скорости. Его детальное определение нам сейчас не требуется, но в дальнейших расчетах удобно считать v, равным скорости движения стационарного пламени с плоским фронтом. Такого определения масштаба скорости мы будем придерживаться в дальнейшем.
Граничные условия. Пусть пламя не ограничено в пространстве, а возмущения имеют локальный характер по направлению у. Тогда необходимые условия имеют следующий вид: x^-оо: в = Ь = 0, x^+x>: дв/дx = дЬ/дx = 0,
у ^±оо: дв/ду = дЬ/ду = 0. (1.4)
Приведем анализ искривленного фронта пламени в приближении бесконечно тонкой зоны химической реакции и произвольных числах Льюиса. В этом случае достаточно рассматривать вместо (1.3) уравнения
дв дв д2в д2в дb дb т (д2b д2bЛ п гч -т—+ w-г— = —т + —т, — + w-г— = Le\ —т +—т I. (1.5) д t д x д x2 д у2 дt д x \д x2 д у2)
Наличие химической реакции учитывается разрывом (при x = 0) первых производных от температуры в и выгорания b [5]. Стационарные решения
в0 = exp(x), b0 = exp(x), x < 0- в0 = b0 = 1, x > 0, w0 = 1,
описывают плоский фронт пламени. При малом искривлении фронта мы ищем решение (1.5) в виде [9]
в^, x, у) = exp( x -£1)« exp( x)(1 _£),
b(t,x,y) = exp^1 4). (1.6)
Здесь £1, £2 _ смещение профилей температуры и выгорания. Причем
£2 = £(1+s), £1 = £l(t> у). (1.7)
Числовой параметр s, который считаем малым, но пока не конкретизируем, задает величину рассогласования распределений температуры и выгорания.
При x^x малые возмущения температуры и концентрации должны обратиться в ноль. Физически это означает, что при малой деформации фронта пламени его температура в приближении бесконечно тонкой зоны химической реакции не меняется. Поэтому решения в=Ь=1 в области x > 0 остаются неизменными с точностью до второго порядка по £1 и его производных. Как мы увидим ниже, такая деформация приводит только к изменению нормальной скорости пламени.
Выражения (1.6) удовлетворяют линеаризованным граничным условиям на фронте химической реакции с точностью до высших порядков по 4 и его производным на границе х=4.
Используя выражения (1.6) в ур. (1.5), учитывая представление у=1+у\ где У - изменение скорости пламени, и удерживая только линейные члены, приходим к следующей паре уравнений 5|1 , д4 д42 , д 42 дг ду дг ду2
Принимая в расчет представления (1.7), получим
дк + у = % ^ + = (1.8)
дг ду дг 1 + е ду2
Если считать скорость пламени не зависящей от кривизны его фронта, то мы должны полагать У=0. Тогда из (1.8) следует равенство 4=0, что означает абсолютную устойчивость пламени, т.к. его фронт просто недеформируемый. Полученный результат не зависит от того, что Ье>1 или Ье<1. Следовательно, 4=0 и при равенстве Ье=1. Значит, фронт пламени с бесконечно тонкой зоной химической реакции недеформируемый, и пламя абсолютно устойчиво к гидродинамическим возмущениям. Этот вывод противоречит теории Ландау-Дарье [5, 10].
Пусть теперь У ^0. Преобразуя (1.8), находим
У , = - (Ье-1)(1 + е) д 24
е ду2 '
Тогда уравнение для поверхности фронта пламени запишется как
= В = 1 + (Ье -!)(! + е).
dt
ду2
1 + Ps(Le -1) д
ду2
(1.10)
Приведенное выше предположение постоянства температуры пламени строго не обосновано. Это видно уже из того, что не были вовлечены в рассмотрение все возможные решения уравнений (1.5). В линейном приближении они, кроме приведенных (1.6), имеют еще соответственно для первого и второго уравнений решения вида [5]
f (x)exp(Qt + iky), g(x)exp(Qt + iky), i = (1.12)
с комплексной частотой Q и волновым числом k, где fx), g(x) удовлетворяют обыкновенным дифференциальным уравнениям
^- f - (Q + k2) f = 0, Le^-g- dg■ - (Q + Lek2)g = 0.
dx2 dx v u ' dx2 dx v
Поэтому мы рассмотрим сейчас задачу о слабо искривленном фронте пламени в расширенной постановке.
2. Расширенный анализ искривленного
фронта пламени
Пусть теперь возмущенные профили температуры и выгорания имеют вид
в=в°(x + |1) + F(t,x,y), b = b0(x + ^2) + G(t,x,y),
где Fи G - новые неизвестные функции, имеющие смысл возмущений. Считая их малыми по абсолютной величине, представим в и b в линеаризованной форме
в=в0( x)-^de0/dx + F (t, x, y), b = b°(x)-^2db°/dx + G(t,x,y), ( .
где e°(x), b°(x) являются решениями уравнений
Уточним вид параметра е. Его появление может быть обусловлено только процессами в зоне химической реакции, других видимых причин нет. Тогда при малой (порядка Р<<1), но конечной толщине зоны химической реакции возмущения температуры и концентрации, вызванные деформацией фронта, будут отличаться. Величина отличия будет порядка в, так как это рассогласование распределений температуры и концентрации вызваны процессами в зоне химической реакции. Тогда е~р. Кроме того, в предельном случае Р=0, когда е=0 и 4=4, пламя должно быть, как мы уже видели, абсолютно устойчивым. Причем, вне зависимости от величины отличия числа Льюиса от единицы: при е=0 из (1.8) следует 4=0. Сказанное означает, что е = *(Ье -1)0, (1.9)
где ^ - неопределенная константа. Тогда формула для изменения скорости пламени приобретает форму
w° 4f = Л + W 0 Ь0),
dx dx2 w° db0 = Le+ W(в°, b0),
dx dx2
(2.2)
где индекс 1 у символа 4 опущен для простоты записи. Фронт пламени будет меняться во времени согласно уравнению [9]
д4= В в =1 +Р*Ье (1.11)
дг В ду2' В Р* '
с первыми четырьмя граничными условиями из (1.4).
Подставив (2.1) в уравнения (1.3) и учитывая равенство w=1+W и (2.2) в промежуточных выкладках, после несложных преобразований получим линеаризованные уравнения
dF + dF + (<!к + dt dx dx I dt dy2
* de0 db0 db0 dxxb2 Ъ1/'
dG + dG + dbL fel + w'_ Le % ] = dt dx dx I dt dy I
= LeД*G + Щ-F + WG-Щ-df^ _4), * de0 db de0 dx Vb2 Ъ1/'
где производные от скорости химической реакции вычисляются через стационарные распределения температуры e° и выгорания b0.
В операторе Лапласа Д дифференцирование производится по переменным x, y. Выражения, стоящие в скобках левых частей приведенных вы-
ражении, следуя рассуждениям предыдущего раздела, приравняем нулю. Эта процедура даст две пары уравнении, первая пара - это уравнения для нахождения смещении 4, 4:
+w =-:нг, -В- + w = Ье-^г.
дt ду дt ду
Вторая же пара - это уравнения для нахождения откликов температуры и выгорания на изменение формы пламени:
¥+М. = AF+дЩ.Р +Щ.0+ (2.3)
дt дx * дв° дЬ дЬ dx
М + д0 _д G + дЖ¥+Ш-а-в, дWdв0 4
дt + дх = Д*0 + дв° ^ + дЬ° 0 в* дв° dx 4
где учтена связь (1.7) между 4 и 4 и опущен индекс "1" у смещения 4 Эволюция 4(1, У) подчиняется уравнению (1.11).
Остается теперь показать, что возмущения температуры ¥ и концентрации продукта реакции О, во-первых, падают с течением времени при произвольных малых начальных деформациях ( плоскоИ поверхности, во-вторых, они локализованы в зоне химическои реакции, а температура пламени не меняется. Для этого численно решались уравнения <¥+ <¥_<Э2¥,52¥+дЖр, а,дWdЬ04 (2.4)
д1 дх "дх2 + ду2 + дв" ¥ дЬ0 0 + в* дЬ" dx 4,
.М. + <д£ _ 1е ш+щ ]+
дг дх ^ дх2 ду2 )
+ дж ¥+Ж О-в, дЖ в 4,
дв дЬ дв" dx ъ'
п _ 1 + в,Ье
д I *ду2' * в* .
Начальные условия:
г=0: ¥=0=0, (=4(у). (2.5)
Граничные условия для ¥ и О:
х^—со: ¥=0=0; х^+да: д¥/дх=д0/дх=0;
у^±ж: д¥/ду=д0/ду=0.
Граничные условия для деформации задавались в двух формах:
1) у^±о: д(/ду = 0; 2) у^±»: (=0. (2.6)
Кроме того, брались различные варианты начального отклонения 40(у).
Стационарные распределения в", Ь получены как асимптотические (^го) решения нестационарной задачи
+ W (в,Ь), дЬ _ Ье Щ- + W (в,Ь), (2.7)
дг дх2 дг дх2
Щ _ д(1 - Ь0)
дв"
1 -д
с начальными условиями
г _ 0:в _Ь _ Дх), /(х) _
1, х ^ +о 0, х ^ -оо'
(2.8)
в[д+в°(1 -д) ]
-ехр
411 --
1
еж
дЬ0
_ -дехр
41 1 --
1
д+в"(1 -д)
Д+в"(1 -д),
в0 _в\х), 6° _ Ь(х).
и граничными условиями
х^-оо: в=Ь=0, х^+о: д(/дх =дЬ/дх=0. (2.9) Производные от Ж в (2.4) вычисляются по стационарным значениям температуры и выгорания и имеют вид
Полная схема решения поставленнои задачи такова: сначала решением уравнений (2.7) с начальными (2.8) и граничными (2.9) условиями находятся установившиеся распределения в°(х), Ь°(х). После этого начинается решение уравнений (2.4), для которых начальные и граничные условия задаются выражениями (2.5, 2.6).
Для численного решения всех уравнений применялась неявная схема. Причем, система (2.4) решалась методом расщепления по пространственной переменной Писмена-Рэдфорда [11] в квадрате 40x40. Функции источников вычислялись по значениям ¥, 0 на предыдущем временном слое.
На рисунке приведены распределения ¥, 0 в различные моменты времени, после возрастания в начальные периоды. Начальное отклонение 40(у) задавалось в виде:
4(у)=А, если уе [19,75; 20,25];
4(у)=0, если уё [19,75; 20,25]; .
А=10 с граничными условиями 2) из (2.6). Значения ¥, 0 берутся в сечении по у максимального отклонения 4(г, у), т.е. при у=20,0, а постоянные параметры д, ,, д, вполагались равными д=20,0; ,=2,0; д=0,15; в=0,1.
Как видно из приведенных графиков, возмущения температуры ¥и выгорания 0 сосредоточены в зоне пламени и убывают со временем. Изменения ¥ носят колебательный характер, причем это наблюдается вне зависимости от того, больше или меньше единицы число Льюиса. В то время как 0 (или концентрация продуктов реакции) стремится к нулю монотонно. Все изменения температуры и выгорания сосредоточены преимущественно в узкой зоне химической реакции, но незначительные затухания колебания могут распространяться и в сторону пламени, где находятся только продукты горения. Это означает, что природа зависимости скорости пламени от кривизны его фронта кроется в процессах, протекающих в самой зоне химической реакции, и не связана с распределением концентрации реагирующего вещества и температуры за пределами этой зоны. Поэтому утверждение пропорциональности е~в вполне обоснованно.
Обращает также на себя внимание большее по амплитуде значение выгорания по сравнению с максимумом температуры. Причем это имеет место вне зависимости от значения числа Льюиса. Если число Льюиса больше единицы, то при деформации фронта пламени в сторону продуктов горения происходит уменьшение температуры и выгорания свежей смеси. Если же деформация происходит в сторону свежей смеси (рисунок), то прирост температуры и выгорания положителен.
0,012' 0,010 0,008 0,006 0,004 0,002 0,000
0,0060,0050,0040,0030,0020,001 ■ 0,000
10
20
30
40
0 10 20 30 х 40 б)
Рисунок. Распределения возмущений температуры F и выгорания G в различные моменты времени: а) t = 2,0; б) t = 10,0. Le = 0,8
Прежде чем делать выводы по результатам настоящей работы, заметим, что однородные части уравнений для F и G совпадают с однородными частями системы (2.3). Так как отклонение у) поверхности пламени от плоского состояния быстро убывает ввиду большого значения коэффициента "диффузии" D,, то условие стремления к нулю F, ^ в сущности, должно совпадать с условием диффузионно-тепловой устойчивости пламени, полученным в приближении неизменности скорости пламени. Расчеты показали, что если числовые значения параметров q, ¡л, в лежат за пределами облас-
ти устойчивого режима горения, возмущения F и G с течением времени неограниченно растут. Изогнутый фронт пламени при этом стремится к плоской форме, согласно (1.11).
Все рассуждения проводились в рамках линейного приближения по величине кривизны фронта пламени, поэтому в пределах допустимости такого приближения решение задачи о пространственно искривленном ламинарном пламени может быть сформулировано следующим образом:
1. форма искривленной поверхности горения явля-а) ется решением параболического уравнения (1.11)
с эффективным коэффициентом диффузии;
2. скорость движения искривленного пламени является суммой скорости плоского пламени и поправки по формуле (1.10);
3. условие существования такого неодномерного пламени определяется диффузионно-тепловой устойчивостью ламинарного горения с постоянной скоростью пламени.
Теперь можно привести численную оценку по порядку величин устойчивого по отношению к гидродинамическим возмущениям размера L, пламени. Для этого достаточно воспользоваться известной формулой Маркштейна [1, 5], в которую необходимо внести множитель 1/(дв):
4п к
1' П-1
Для наблюдаемого экспериментально устойчивого пламени коэффициент и теплового расширения находится в пределах 6...9, в~0,1, к~10-5 м2/с, и„°~0,1 м/с [5]. Считая n/s(n-1)»1, находим Х,~0,01 м. Отсюда получаем для критического числа Рейнольдса согласующуюся с экспериментом оценку Кес1~103.
Что касается определения константы д, то этот вопрос остается открытым и находится на стадии проработки.
0
F
СПИСОК ЛИТЕРАТУРЫ
1. Нестационарное распространение пламени / Под ред. Дж.Г. Маркштейна. -М.: Мир, 1968. -220 с.
2. Рогоза Б.Е. О маркштейновской поправке к нормальной скорости распространения пламени // Физика горения и взрыва. - 1985. -Т. 21. -№ 3. -С. 45-48.
3. Зельдович Я.Б., Истратов А.Г., Кидин Н.И., Либрович В.Б. Гидродинамические течения и устойчивость искривленного фронта при распространении пламени в каналах. - М.: 1980. (Препр. ИПМ АН СССР; № 143).
4. Игнатьев С.М., Петухов Ю.И. Нелинейный анализ ячеистой структуры фронта пламени с учетом гидродинамических и диффузионно-тепловых процессов // Физика горения и взрыва. -1989. -Т. 25. - № 5. - С. 58-62.
5. Зельдович Я.Б., Баренблатт Г.И., Либрович В.Б., Махвила-дзе Г.М. Математическая теория горения и взрыва. - М.: Наука, 1980. -422 с.
6. Karpov A.I., Galat A.A., Bulgakov V.K. // Combustion Theory Modeling. -1999. -V. 3. - P. 535-546.
7. Сабденов К.О. К линейной теории искривленного пламени // Инженерно-физический журнал. — 2001. — Т. 74. — № 5. — C. 81-86.
8. Сабденов К.О. О диффузионно-тепловой неустойчивости ламинарного пламени // Инженерно-физический журнал. -2002. - Т. 75. - № 4. - C. 73-79.
9. Сабденов К.О. К разрешению парадокса Л.Д. Ландау о гидродинамической неустойчивости пламени // Фундаментальные и прикладные проблемы современной механики. - Томск: Изд-во ТГУ, 1998. -С. 84-85.
10. Ландау Л.Д. К теории медленного горения // Журнал экспериментальной и теоретической физики. -1944. - Т. 14. - № 6. -С. 240-244.
11. Самарский А.А. Введение в численные методы. - М.: Наука, 1982. -202 с.