CC BY
88
СПИСОК ЛИТЕРАТУРЫ
1. Валеев К.Г, Карелова О.Л., Горелов В.И. Оптимизация линейных систем со случайными коэффициентами. - М.: Изд-во РУДН, 1996.
2. Гихман И.И., Скороход А.В. Стохастические дифференциальные уравнения. - Киев: Наукова думка, 1968.
3. Кац И.Я. Метод функций Ляпунова в задачах устойчивости и стабилизации систем случайной структуры. - Екатеринбург, 1998.
4. Мильштейн Г.Н., Репин Ю.М. О воздействии марковского процесса на систему дифференциальных уравнений // Дифференциальные уравнения. - 1969. - Т. 5. - № 8. - С. 1371-1384.
5. Тихонов В.И. Миронов В.А. Марковские процессы. - М.: Советское радио, 1977.
УДК 536.46
ТЕОРИЯ СПОНТАННОЙ ДЕТОНАЦИИ В ГАЗАХ. Ч. 2. МОДЕЛИРОВАНИЕ ВЗРЫВНЫХ ПРОЦЕССОВ
К.О. Сабденов
Томский политехнический университет E-mail: [email protected]
Рассматривается система уравнений газовой динамики горения на основе исходных понятий «нормальная скорость» и «поверхность горения». Турбулентное пламя опиывается одним нелинейны/м параболическим уравнением. Проведеныi расчеты взрывных процессов в трубах для ряда горючих смесей. Сравнение теоретических данных по длине и времени перехода медленного горения в детонацию с экспериментальны/ми результатами показы/вает их удовлетворительное согласие.
Изложенные в [1] рассуждения привели к представлению о турбулентном пламени в трубе радиуса а как о хаотически блуждающей поверхности с фрактальной размерностью (, где ее элементарный участок движется относительно газа с нормальной скоростью ып ламинарного пламени. Поверхность горения топологически может быть и многосвязной. Весь газ состоит из двух компонент - свежей смеси с массовой долей С (концентрацией) и продукта горения. Как и С, все термодинамические и гидродинамические параметры потока газа по отношению к гидродинамическому хаосу имеют такой же смысл средних величин, который они имеют в не турбулентных средах по отношению к хаосу молекулярному. Т.е. остающиеся макроскопическими временные и пространственные масштабы турбулентности много меньше аналогичных масштабов рассматриваемых ниже процессов. Для пламени в таком макроскопическом описании можно ввести коэффициент диффузии Б [2]:
D = 2aun (1 + B—)4 un
s = df - 2,
(1)
содержащей турбулентную пульсацию ы' скорости потока газа. Превращение исходной свежей смеси в продукт горения описывается выражением для скорости химической реакций Ф
ф(С) = (i + в iL-)s C(1 - C)2 a Un
(2)
приходящим на смену закону Аррениуса. Свободный параметр В зависит от свойств газовой смеси.
Система уравнений газовой динамики горения
Пусть р, р, Т и ы - средние по сечению трубы давление, плотность, температура и скорость газа. В общих случаях без ограничений на скорость дви-
жения турбулентного пламени, газодинамические параметры и концентрация свежей смеси С, могут быть найдены решением уравнений [2]
^ + и^ = 1^1 Т Гдр + и+ ЯФ(С),
дt дх у р \ о1 дх) ср
рдС + ри = 4- ЪрЩ--рФ(С), н д^ н дх дх н дх н у ь
f PU +1 (PU2
p) = -—pu \u\,
-dp- + p = 0, p = pRgT, D = 2una ■ K(u', Un),
Ф(С) = —к(u', un) ■ C(1 - C)2, K(u' , un) = (1 + B u.)s,
a un
(3)
где 0 - тепловой эффект горения горючей газовой смеси с теплоемкостью ср при постоянном давлении и газовой постоянной Д,; у - показатель адиабаты; ?, х - время и координата вдоль оси трубы.
Уравнения (3) не замкнуты: не хватает связи между пульсацией скорости и скоростью ы потока газа. Также система (3) должна быть дополнена связью фрактальной размерности с параметром, характеризующим турбулентное пламя. Таких параметров здесь две: число Рейнольдса Re и отношение ы/ып. В дальнейшем мы будем пользоваться простейшей формулой замыкания [3] ы=^Я|ы|. Что касается фрактальной размерности поверхности горения, то функциональную зависимость (==((Ке,ы /ып) пока не удается установить.
На ряде примеров покажем непротиворечивость ур. (3) ранее установленным теоретическим и экспериментальным фактам. Известно [4], что с точностью до ошибок эксперимента концентрационные и тепловые пределы ламинарного горения и детонации или совпадают, или мало различаются. В качестве примера в табл. 1 приведены экс-
периментальные данные по концентрационным пределам детонации и нормального распространения пламени ряда газов в смеси с воздухом [4]. Из (3) следует, что эти пределы совпадают. Действительно, если и„=0, что имеет место за концентрационными пределами, то ур. (3) будут описывать химически инертную среду, где заведомо не может быть ни детонации, ни медленного горения. Таким образом, сформулированная модель горения газов по концентрационным пределам двух предельных режимов - дефлаграции и детонации не противоречит эксперименту (возможность описания детонационного горения показана ниже).
Таблица 1. Концентрационные пределы детонации и медленного горения для ряда горючих газов в смеси с воздухом
Горючее Концентрационный предел детонации (и нормального распространения пламени), %
Нижний Верхний
Этан 2,87 (3,0) 2,57 (2,20)
Пропан 1,98 (1,90) 3,32 (3,10)
Бутан 3,55 (2,25) 12,2 (12,5)
Этилен 7,37 (9,50) 6,18 (8,50)
Пропилен 14,7(32,0) 10,4 (10,8)
В [5] анализом ур. (3) показано, что в однородной турбулентной среде (и'=еош11) пламя распространяется со скоростью
uf = un (1 + Bu' / un)1
(4)
и что она является собственным значением нелинейной краевой задачи. Последнее обстоятельство является одним из важнейших математических характеристик волн горения [6]. При отсутствии турбулентности скорость пламени в точности равна экспериментально наблюдаемому значению: «=«„. Если интенсивность турбулентности высока, то выражение (4) переходит в формулу Вилюнова-Кли-мова. Т.о., используемые представления горения, основанные на понятиях «поверхность горения» и «нормальная скорость пламени» не противоречат известным полуэмпирическим и статистическим подходам в теории турбулентного горения. Более того, если ранее число 5 в формуле Вилюнова-Кли-мова присутствовало как формальный параметр, то здесь оно имеет явную геометрическую трактовку, указывая на меру сложности поверхности горения.
Расчет перехода медленного горения в детонацию
Переход медленного горения в детонацию бы-строгорящих смесей изучался на основе системы (3), в которой первое уравнение переписано в удобной для численного решения форме и всем уравнениям придана дивергентная форма. Переход к безразмерным величинам проводился по следующим масштабным параметрам
и, = ипМ, р„=р0, р,= р0ми2, ( , = а /ип, х,, = аМ, М = ^р0 / р0ип2.
Здесь р0, р0 - начальные значения давления и плотности газовой смеси. Численно решались уравнения в безразмерных переменных:
д ( 1 и2
Р^
д t \у-1
Y
p + р — 1= M1KpC (1 - С)2
д
дх 1 y -1 2
д „ д „ д (2Kр дС
— рС + — риС =—I-2--
д t дх дх I M2 дх
- KрС(1 - С)2
д д
— ри +--(ри2 + p) = -\M ри\и\
д t дх
K = (1 + BM-J\\u |)1
д д
— р +--ри = О,
д t дх
(5)
\= 0,0032 -
О,221
(Re U |)°Д
, Re =
2aunM
m1=м.
Р 0
Выражение для коэффициента сопротивления Я (формула Никурадзе) взято из [7]. Коэффициент кинематической вязкости V считался постоянным.
Сформулируем начальные условия. При выводе выражений для Ф(С) (1) и Б (2) заведомо предполагается наличие фронта горения. Это означает, что ур. (3) или (5) принципиально не могут описывать стадию воспламенения горючей смеси. Начальной момент горения должен моделироваться методами классической теории горения, изложенной, например, в [6]. Поэтому при решении (5) начальные условия необходимо задавать в виде сформировавшегося фронта пламени.
Потери тепла в стенку трубы в расчетах не учитываются. Начало взрывного горения происходит в центре трубы неограниченной длины. Это эквивалентно зажиганию у закрытой части при адиабатических условиях на жестком торце.
Начальные условия имеют вид: р = 1,5, С = 0, и = 0, р = 0,3, для — £< х < +е; (6) р = 1, С = 1, и = 0, р = 1 - для других х.
Константа е подбиралась таким образом, чтобы размер области «начального» горения для всех смесей было равным 3а. Показатель адиабаты и коэффициент кинематической вязкости брались равными 1,4 и 1,5 10-5 м2/с. Радиус трубы а=1,3-10-2 м.
Система ур. (5, 6) решалась численно по явной 2-го порядка точности по пространственной координате схеме Годунова-Колгана.
Моделирование перехода медленного горения в детонацию проводилось для смесей 4Н2+О2,
2Н2+О2, Н2+О2,
Н2+2О2,
CH4+2O2, С2Н2+воздух,
2СО+О2, по которым имеются наиболее полные экспериментальные данные, включая условия проведения опытов. Все начальные данные соответствуют нормальным условиям.
Сравнение расчетной длины Яв преддетона-ционного участка и времени Б перехода медленного горения в детонацию с экспериментальными
данными [8-10] показывает их удовлетворительное согласие (табл. 2). Основные размерные параметры (начальные давление и плотность топливной смеси, тепловой эффект реакции, нормальная скорость пламени и пр.), определяющие Яв и ¡в, отвечают условиям проведения экспериментов.
Если параметр В подобран для согласования расчета с экспериментом, то близость значения показателя 5 к единице легко обосновать. В экспериментальных исследованиях факельного турбулентного горения газовых смесей значение 5«0,7...0,8 [10]. При сильной турбулентности, имеющей место в период ускорения пламени, фрактальная размерность поверхности горения 5^1. Но предельно допустимое значение 5=1 (т.е. 4=3) реализуется лишь при Re=<». Поэтому здесь и далее 5=0,8.0,9.
Таблица 2 Сравнение теории с экспериментом по определению протяженности преддетонационного участка Нв/2а и времени перехода %
Смесь ВЫ Яв/2а tD, мс
теор. эксп. теор. эксп.
2Н2+02 25 (0,8) 45 38 1,4 1,2
Н2+2О2 15 (0,9) 86 92 9,5 8,3
4Н2+02 30 (0,8) 54 79 5,9 3,1
Н2+О2 30 (0,9) 44 41 1,6 1,8
СН4+202 20 (0,9) 57 56 2,1 2,1
С2Н2+воздух 22 (0,9) 67 92 5,7 6,1
2С0+02 40 (0,9) 39,6 37,5 3,8 3,0
1,8 1,6 1,4 1,2 1,0 0,8 0,6 0,4 0,2 0,0
а)
0 0,5 1,0 1,5 2,0 2,5 3,0
б)
0,0 0,5 ' 1,0 ' 1,5 ' 2,0 ' 2,5 ' 3,0
Эволюция распределения давления и скорости для смеси Н2+202 приведена на рис. 1-3, для смеси СН4+202 - на рис. 4-6. Центру симметрии соответствует х=0. На остальных рисунках индексы у символов указывают на приведенные в подрису-ночных подписях моменты времени.
При сверхзвуковой скорости движения турбулентного пламени в его области расположения наблюдаются быстрые пространственные изменения всех физических параметров (рис. 1, б), характеризующих горение. Так как значение приведенной скорости химической реакции Ф(С) практически целиком определяется турбулентной составляющей, то положение максимумов на профиле скорости и совпадает с максимумами Ф(С). Окрестности точек максимумов функции Ф(С) можно считать, как и в классической теории ламинарного горения, областями локализации турбулентного пламени. Они характеризуются высокой температурой и быстрым изменением концентрации исходного продукта. На рис. 1, б, видны два следующих друг за другом турбулентных пламени. Такая двойственность объясняется тем, что за первым фронтом турбулентного пламени формируется волна разрежения, где скорость газа (как и давление и плотность) быстро падает. Это влечет за собой снижение скорости химической реакции Ф(С). В результате исходная смесь не успевает сгореть. Во втором фронте пламени снова активизируется скорость горения, что вновь сопровождается интенсификацией горения, ростом давления, скорости и плотности газа.
Рис. 1. Стадия ускорения фронта пламени в смеси Н-2+20-2 для I: а) 7,2, б) 9,0 мс: 1 - р, 2 - и, 3 - р, 4 - С
Скорость движения второго фронта турбулентного горения превышает скорость первого. Это обусловлено тем, что на процессы горения, согласно принятой здесь модели, существенным образом сказывается влияние газодинамических факторов. Скорости движения фронтов горения зависят от характера эволюции (следующих друг за другом) ударных волн, которые сами же их генерируют.
Если скорость пламени меньше скорости звука в исходной смеси, то пространственное распределение газодинамических параметров и концентрации С гладкое, без ярко выраженных изломов (рис. 1, а). Здесь максимальная скорость газа составляет несколько десятков м/с, а скорость турбулентного пламени относительно газа немногим превышает скорость ламинарного горения ип.
На рис. 1, а, отчетливо видна в точке х«2,6 слабая волна сжатия впереди фронта пламени, наличие которой обусловлено начальными данными (6).
На рис. 2 приведены этапы формирования детонационного режима горения. Вторая волна турбулентного горения здесь уже догнала переднюю волну, в результате чего образовалась одна область турбулентного горения с более высокой скоростью химической реакции, большими скачками газодинамических параметров, чем первые две. После такого слияния волн турбулентного горения переход медленного горения в детонацию еще больше ускоряется.
х
6
3
2
0
х
40'
30'
20'
10'
а)
0,0
i ■ i ■ i ■ i ■ i ■ i
0,5 1,0 1,5 2,0 2,5 3,0
-1,5
Рис. 2.
40'
30'
20'
10'
б)
-I-1-1-1-1-1-1-1-1-1-1
0,5 1,0 1,5 2,0 2,5 3,0
а)
0,0 0,5 1,0 1,5 2,0 2,5
X
3,0
б)
1
3,0
Переход пересжатой детонационной волны в стационарную в смеси Н2 +20-1
Волнообразный вид распределений давления, скорости и плотности за фронтом пламени вызван блуждающими возмущениями, возникшими при догорании остатков свежей смеси, и частично - при слиянии двух первоначальных фронтов турбулентного горения. Эти волны, взаимодействуя с ударной волной и «стенкой» (с центром симметрии), с течением времени исчезают. На достаточно больших временах, когда завершено формирование стационарной детонационной волны, а последняя значительно удалена от стенки, блуждающие волны практически не наблюдаются. Максимум давления находится всегда немного впереди пламени. Положение последней определяется по максимуму скорости химической реакции. Образовавшаяся впоследствии пересжатая детонационная волна (рис. 3) выходит с течением времени (с точностью до ошибок аппроксимации) к режиму, соответствующему на детонационной адиабате точке Чепмена-Жуге. Экспериментальное значение скорости детонационной волны в смеси Н2+202«1927 м/с. Расчетная же скорость в используемой модели - 1857 м/с. Также несколько заниженной оказалось отношение плотности газа в ударной волне рув к р0: 5,6 вместо 6,0.
3,0'
Возникновение детонационной волны в смеси Н-+20- для 1:1) 9,4, 2) 9,7, 3) 9,9 мс
0,0
3,5' 3,0' 2,5' 2,0' 1,5'
1
a)
—I—
0,5
—I—
1,0
—I—
1,5
—I
2,0
Л
:: 0,0
б)
0,0
—i—
0,5
—i—
1,0
—i—
1,5
—I
2,0
Рис. 4. Ускорение слабо турбулентного пламени в смеси СН4+202 для I: а) 0,79, б) 1,18 мс; 1 -р, 2 - и, 3 -р, 4 - С
Качественная картина развития процесса перехода медленного горения в детонацию для смеси СН4+202, не отличается от той, что имеет место для смеси Н2+202. Здесь тоже могут наблюдаться два фронта (волны) турбулентного горения (рис. 4, 5).
0
X
X
3
4
2
0
X
1
0
3
4
2
X
На рис. 4 показано последовательное развитие второго фронта, который движется и усиливается быстрее первого. В нем скорость турбулентного горения пока еще относительно низка. Если пик давления в первой волне за время от 0,79 (а) до 1,18 (б) мс увеличился всего на ед. %, то во второй волне давление возросло примерно в 1,4 раза от первоначального значения. Эволюционное развитие второго фронта турбулентного горения идет столь быстро, что через еще 0,24 мс (рис. 5, а) амплитуда давления в нем становится в 1,68 раза больше амплитуды давления в первом фронте горения. 8т 7 6 5 4 3 2
а)
0,0
12 10 8 6 4 2
0,5
1,0
X
1,5
2,0
б)
1,0
X
Рис. 5. Фронт турбулентного пламени с двумя очагами горения в смеси СН4+202 для 1: а) 1,42, б) 1,71 мс; 1 - р, 2 - и, 3 - р, 4 - С
Также сильно различаются и скорости движения указанных волн турбулентного горения. Скорость первой волны меняется со временем слабо, на рис. 4, 5 она «1100 м/с. Скорость движения второй волны быстро растет: если в интервале времени на рис. 4 она имеет среднюю скорость «300 м/с, то в промежутке от 1,18 (рис. 4) до 1,58 мс (рис. 5) составляет «2000 м/с.
Несмотря на незначительное отличие этой величины от скорости стационарной детонационной волны 2146 м/с пики давления и плотности в ударной волне впереди фронта пламени (рис. 5) в 1,5 и более раз ниже, чем при детонационном горении. Слияние двух волн турбулентного горения в одно целое происходит ко времени 1,86 мс. В этот момент область максимального давления не имеет остроконечного пика (рис. 6, а), как при других временах.
50
40
30
20
10
0,0
а)
0,0 0*5 ' 10 ' 1,5 ' 2,0
б)
0,0 0,5 1,0 1,5 2,0
X
р
-г
0,5
-I-
1,0
т-
1,5
т
2,0
в)
Рис. 6. Завершающая стадия возникновения детонационной волны в смеси СН4 +202 для £ 1) 1,86, 2) 1,97, 3) 2,1 мс
В смеси CH4+2O2 наблюдается переход медленного горения к стационарной детонационной волне через пересжатое состояние. Переход из пересжатой в обычную детонационную волну происходит с отрывом от фронта детонации небольшого возмущения. Это возмущение, слабо заметное на рис. 3, легко обнаруживается на рис. 6, а, как выпуклость на профиле давления позади детонационного пика. Оно распространяется в противоположную сторону по отношению к направлению движения детонационной волны. Ретонационная волна, возникающая в смеси 4H2+O2, имеет ярко выраженную форму ударной волны (рис. 7). При проведении расчетов для этой смеси обнаружено, что ам-
X
4
3
2
X
плитуда ретонационной волны растет по мере увеличения амплитуды ударной волны, идущей впереди ускоряющегося фронта турбулентного пламени.
X
Рис. 7. Ретонационная волна (х = 0,5) и волна начального возмущения (х = 6,4) в смеси 4Н2+02. 1 - р, 2 - и, 3 - р, 4 - С
Физическая природа появления обратной волны заключается в необходимости «уноса» некоторого избытка энергии из пересжатой детонационной волны для того, чтобы последняя могла выйти к состоянию стационарной детонационной волны.
СПИСОК ЛИТЕРАТУРЫ
1. Сабденов К.О. Теория спонтанной детонации в газах. Ч. 1. Формулировка основных положений // Известия Томского политехнического университета. - 2005. - Т. 308. - № 3. - С. 16-22.
2. Сабденов К.О., Миньков Л.Л. К фрактальной теории перехода медленного горения в детонацию в газах // Физика горения и взрыва. - 1998. - Т. 34. - № 1. - С. 70-78.
3. Зельдович Я.Б. К теории возникновения детонации в газах // Журнал технической физики. - 1947. - Т. 17. - № 3. - С. 3-26.
4. Нетлетон М. Детонация в газах. - М.: Мир, 1989. - 280 с.
5. Сабденов К.О. Нестационарные задачи горения газовых смесей, жидких и твердых взрывчатых веществ и ракетных топлив: Дис. ... канд. физ.-мат. наук. - Томск, 1999. - 115 с.
Отрицательное направление движения волны возмущения, которую называют ретонационной волной [8, 10], следует из закона сохранения импульса.
Выводы
- Уравнения газовой динамики горения газовых смесей допускают формулировку на основе исходных понятий «поверхность горения» и «нормальная скорость» распространения пламени.
- В такой формулировке уравнений газовой динамики горения при моделировании турбулентного пламени отсутствует проблема замыканий уравнений типа Фридмана-Келлера, но возникает другая - проблема нахождения численных значений меры Хаусдорфа и фрактальной размерности поверхности горения.
- Разработанная теория дает корректную физическую картину перехода медленного горения в детонацию, а рассчитанные для газовых смесей теоретическое время и протяженность предде-тонационного участка удовлетворительно согласуются с экспериментальными данными. Также согласие теоретических результатов с экспериментом имеет место по возникновению пересжатой детонационной и ретонационной волн.
6. Зельдович Я.Б., Баренблатт Г.И., Либрович В.Б., Махвила-дзе Г.М. Математическая теория горения и взрыва. - М.: Наука, 1980. - 478 с.
7. Лойцянский Л.Г. Механика жидкости и газа. - М.: Наука, 1987. - 840 с.
8. Бартльме Ф. Газодинамика горения. - М.: Энергоиздат, 1981. -312 с.
9. Щелкин К.И. К теории возникновения детонации в газовых смесях в трубах // Доклады АН СССР. - 1939. - Т. 23. - № 3. -С. 636-640.
10. Льюис Б., Эльбе Г. Горение, пламя и взрывы в газах. - М.: Мир, 1968. - 592 с.
