CC BY
76
2013
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Математика и механика
№ 2(22)
МЕХАНИКА
УДК 536.46
А.А. Дементьев, А.Ю. Крайнов ИССЛЕДОВАНИЕ ВЛИЯНИЯ ОТНОСИТЕЛЬНОГО ДВИЖЕНИЯ
ВЗВЕСИ ИНЕРТНЫХ ЧАСТИЦ НА СКОРОСТЬ ФРОНТА ГОРЕНИЯ
ГАЗОВОЙ СМЕСИ
Исследуется зависимость скорости горения газовой смеси с инертными частицами от концентрации горючего, взятого в стехиометрическом соотношении с окислителем и в недостатке. Приведено сравнение результатов расчетов, полученных на основе двух математических моделей: с учетом теплового расширения, без его учета. Получены зависимости скорости распространения фронта горения от размера и массовой концентрации инертных частиц при различных начальных концентрациях горючего.
Ключевые слова: газовзвесь, ламинарное пламя, инертные частицы, скорость горения.
Газодисперсные системы, к которым относятся различные газовзвеси инертных частиц, а также частиц органических горючих, часто используются в химической, пищевой промышленности, инертные порошки используются для пожаротушения [1-5]. Исследование и моделирование процессов, происходящих при воспламенении и горении газовзвесей, изучение качественных и количественных характеристик - скоростей горения, температур, плотностей и т.д., является необходимым для эффективного развития промышленных технологий.
В данной работе представлены результаты теоретического исследования горения газовой смеси со взвешенными в ней мелкими инертными частицами с учетом теплового расширения газа. Проведено сравнение результатов с результатами вычислений по модели, в которой не учитывается тепловое расширение газовой фазы и движение частиц.
Физико-математическая постановка задачи
Горючая газовая смесь с взвесью инертной пыли занимает область 0 < х < ®. Считается, что реакция в газовой фазе проходит по второму порядку. Частицы одного размера, сферической формы, с радиусом гк, в начальный момент времени распределены равномерно в газовой смеси. В координате х = 0 размещена горячая непроницаемая стенка с температурой равной адиабатической температуре сгорания газа в изобарических условиях. Такое граничное условие моделирует зажигание у «закрытого торца» [6]. Горение протекает в ламинарном режиме при постоянном давлении [7-9]. Предполагается, что число Льюиса и коэффициент теплопроводности не зависят от температуры; теплообмен между газовой фазой и частицами описывается законом Ньютона; при движении на частицы со стороны газа
действует сила сопротивления, описываемая законом Стокса, из-за малой объемной концентрации частиц пренебрегаем их обратным влиянием на движение газа.
Система уравнений, описывающая процессы в горючем газе с инертными частицами, в безразмерной форме запишется в виде
<90 _ д0g і д2Є,, pk B / \ _
—- + u—- =---------2---к— ( - 0к Кр^^ехр
дт д| р- д|2 р- ^ - к} ^ 1 2
1 + в
. дС1 = Ье д2С1 _
Р- д|2
дС _ _
+и = ^ :Ът - -С1С2 ехР
дт д£
1 + в
дС2
■ + и
дт д£,
дС2 _ Ье д С2 ГУ2^2 ЇІ0 \ —
р- д|2
^іМі
20
УР-С1С2 еХР
1 +(
др8 дР-и п —- +—— = 0 ;
дт д£,
дйі+£М.=0.
дт д£,
д®к -д0к
—- + п—-
дт д£,
дРкп + дРкп
2
(08 -0к ).
3 РтС_ _ч рк(и - п);
дт д| 2 уЫи
р- =(1 -^У)/(і + Р0- ) .
Граничные условия:
д0- (°°,т)
(1)
(2)
(3)
(4)
(5)
(6) () (8)
0- (0, т) = 0 дС2 (0, т)
=0,
= 0,
д|
дС1 (да, т)
д0к (0,т)=п дРк (0,т)
=0,
= 0
д|
дС2 (да, т) д^
=0
дС1 (0, т)
д|
=0,
= 0, и (0, т) = 0, п (0, т) = 0. (9)
д£, д£,
Начальные условия:
0я (5,0) =-1/у, ек (1,0) = -1 у, С[ (|,0) = Сш, с2 (|,0) = 1,
и (|,0) = о, * (|,0) = 0, Рк (5,0) = 1. (10)
Уравнение (1) описывает перенос тепла в газовой смеси за счет теплопроводности, конвекции, с учетом теплообмена между газовой фазой и частицами, тепловыделения за счет химической реакции; (2) - уравнение концентрации горючего; (3) - уравнение концентрации окислителя; (4), (5) - уравнения сохранения массы газа и частиц; (6) - уравнение энергии частицы, с учетом ее движения и теплообмена с газом; (7) - уравнение движения частиц. (8) - уравнение состояния идеального газа, записанное при постоянном давлении.
При записи системы уравнений выбраны следующие масштабы: времени и = КТаС& / ^Ерё0¥10У2окоехр(-Е/(ЯТа))); координаты х* = (^*/(сгрг0))1/2; скорости
V* = х*И*; плотности газа р^ - начальная плотность газа; приведенной плотности частиц рк0 - начальная приведенная плотность частиц; температуры Т* = ЯТа2/Е, где Та - адиабатическая температура продуктов сгорания газа в изобарических условиях. Безразмерные переменные: т = t/t* - время; £ = х/х* - координата; 0&к = (Е/ЯТа2)/(Т^к - Та) - температура; С1 = 71 / 710, С2 = 72 / У20 - относительные концентрации горючего и окислителя в газовой фазе, где 71 =р1 / , 72 =р2 / ря,
ря = р^ / р^0 - плотность газа; рк = рк / рк0 - приведенная плотность частиц;
м = и / V* - скорость газа; V = V / V* - скорость частиц. Индексы: g - газ, к - частицы, 0 - начальные условия, 1 - горючее, 2 - окислитель, * - масштаб, черта над обозначением - безразмерная величина.
Безразмерные параметры (критерии): Рг = \с^Х - число Прандтля, Ми = агк/Х -число Нуссельта, Ье = Бср^Х - число Льюиса, С = ск /с& - отношение теплоемкостей, В = 4пrk3ckpk0n/(3cgpg0) - начальная массовая концентрация частиц,
X = гкскрк0/2Ыики - размер частиц, у = сЯТа2/(<2ЕУ10), в = ЯТа/Е - параметры, характеризующие тепловой эффект и энергию активации химической реакции.
Обозначения: t - время, х - координата, Т - температура, 7 - отношение плотности окислителя/горючего к плотности газовой смеси, р - плотность, и - скорость газа, V - скорость частиц, с - удельная теплоемкость, X - коэффициент теплопроводности, О - коэффициент диффузии, п - коэффициент вязкости, д - молекулярная масса, V - стехиометрический коэффициент, а - коэффициент теплоотдачи, гк - радиус частиц, п - количество частиц в единице объема, Я - универсальная газовая постоянная, Q - тепловой эффект реакции, к0 - предэкспоненци-альный множитель, Е - энергия активации химической реакции.
Математическая модель (1) - (10) описывает распространение пламени в газовой смеси с инертными частицами с учетом теплового расширения. При и = 0, V = 0, рё = 1, модель переходит в теплодиффузионную модель горения газодисперсных сред [10-13].
Система уравнений (1) - (10) решалась численно. Уравнения (1) - (3) решались по неявной разностной схеме, конвективное слагаемое аппроксимировалось разностями против потока. Значения температуры газа на новом временном слое использовались для нахождения плотности газа из уравнения (8). Затем из уравнения сохранения массы (4) определяли скорость газа. Уравнения (5) - (7) решались по явной разностной схеме.
Расчеты проводились при значениях параметров: у = 0,1, в = 0,07, Ье = 1, Ыи = 1, VI д = 16, V2Д2 = 64, 710 = 0,05, 720 = 0,2, Рг = 1, С = 1, В = 0,15. В расчетах варьировались значения С в интервале от 1 до 0,7, параметр % в интервале от 2000 до 500000, В в интервале от 0 до 0,25.
В численных экспериментах горение газовзвеси инициировали передачей тепла от горячей стенки. В момент времени, когда фронт пламени достигал безразмерной координаты £ = 500, граничное условие 0g (0, т) = 0 заменялось граничным условием д0g (0, т)1д| = 0. Скорость фронта определялась скоростью перемещения координаты с концентрацией горючего С1 = 0,5.
На рис. 1 представлены распределения параметров среды в последовательные моменты времени с интервалом безразмерного времени Ат = 120000 для модели без учета теплового расширения, на рис. 2 - для модели с учетом расширения при значениях % = 50000, С = 0,8, В = 0,15.
Температура частиц во фронте горения ниже, чем у газовой смеси, за фронтом их температуры выравниваются (рис. 1, а, 2, а). Концентрация горючего, взятого в недостатке, за фронтом пламени обращается в ноль. Окислитель, взятый в избытке, расходуется до концентрации С2 = 0,2 (рис. 1, б, 2, б).
Сі,С2
Рис. 1. Распределения температуры (а) (сплошная линия - газ, пунктир -инертные частицы), концентраций (б) (сплошная линия - горючее, пунктир - окислитель). Модель без учета расширения
0
Сі,С2
и, Ж
Р, Рк
Рис. 2. Распределения температуры (а), концентрации (б), скорости (в), плотности (г); (а, в, г: сплошная линия - газ, пунктир - инертные частицы), (б: сплошная линия - горючее, пунктир - окислитель). Модель с учетом расширения
0
Перед установившимся фронтом горения скорость газа и частиц имеет одинаковое значение, позади фронта скорость газа и частиц равна нулю, во фронте прогрева газовой фазы частицы имеют скорость выше скорости газовой фазы из-за их механической инерции (рис. 2, в). Частицы в области химических превращений имеют увеличенную приведенную плотность (рис. 2, г). В силу этого эффекта в зоне интенсивных химических реакций локально присутствует большее количество частиц, которые локально понижают температуру в зоне реакции сильнее, чем это получается в расчетах на основе модели без учета расширения газа и относительного движения частиц. Подогретые в зоне химических реакций частицы и имеющие скорость большую, чем газ, переносят тепло в зону прогрева перед фронтом горения. Это способствует растяжению волны горения в газе.
На рис. 3 представлены зависимости скорости горения от начальной концентрации горючего, полученные на основе модели без учета теплового расширения (рис. 3, а) и с его учетом (рис. 3, б) при массовой концентрации частиц В = 0,15 и трех значениях параметра х, характеризующего размер частиц. Видно, что тепловое расширение увеличивает скорость распространения фронта пламени в чистом газе на величину от 20 до 30 % (кривые 1 на рис. 3, а, б). Присутствие крупных инертных частиц в газе не сказывается на скорости распространения пламени в модели без учета теплового расширения газа (кривая 2, рис. 3, б). Расчет на основе модели с учетом теплового расширения показывает, что скорость пламени уменьшается, несмотря на то, что частицы во фронте пламени не успевают нагреться, забрав тепло из зоны химической реакции. Этот эффект обусловлен меньшей температурой газа позади фронта реакции и в связи с этим меньшим объемным расширением газа. С уменьшением размера частицы начинают влиять на скорость химической реакции во фронте горения, уменьшая в нем температуру (кривые 3, 4, рис. 3, а, б). Существует интервал значений х, при которых частицы слабо влияют на скорость пламени при стехиометрическом и близком к нему соотношении окислителя и горючего и сильно влияют на скорость при значительном отклонении состава газовой смеси от стехиометрической (кривые 3, рис. 3, а, б).
Рис. 3. Зависимости значений скоростей фронта пламени от концентрации горючего: а - модель с учетом расширения (1 - незапыленная газовая смесь, 2 - х = 500000, 3 - х = 50000, 4 - х = 2000); б - модель без учета расширения (1 - незапыленная газовая смесь, 2 -х = 300000, 3 - х = 50000, 4 - х = 2000)
На рис. 4 представлены зависимости относительной безразмерной скорости распространения фронта пламени от параметров дисперсной фазы - В и х при двух значениях начальной концентрации горючего в газовой смеси (С1 = 0,85 и
С1 = 1), рассчитанные по модели, учитывающей тепловое расширение. У0 - значение скорости фронта горения для чистой газовой смеси, которое имеет значение при С1 = 1 - У0 = 0,02, при С1 = 0,85 - У0 = 0,01. Параметр х изменялся в следующем интервале: 2000 < х < 500 000, параметр В : 0 < В < 0,25. Из рис. 4 видно, что при фиксированном х с увеличением массовой концентрации скорость пламени уменьшается. Чем мельче частицы, тем их влияние на скорость пламени значительнее.
Рис. 4. Зависимости значений относительных скоростей фронта пламени от параметров дисперсной инертной фазы, пунктир - концентрация C1 = 0,85, сплошная линия - C1 = 1: а - зависимость от величины, характеризующей размер частиц % (1 - B = 0,05, 2 - B = 0,1, 3 -B = 0,15, 4 - B = 0,2, 5 - B = 0,25,); б - зависимость от массовой концентрации B (1 - % = 500 000, 2 - % = 120 000, 3 - % = 50 000, 4 -X=2000)
Сформулирована математическая модель горения газовой смеси с взвешенными в ней инертными частицами, учитывающая тепловое расширение газа и относительное движение фаз. Получены зависимости скорости распространения пламени от массовой концентрации, размера инертных частиц и концентрации горючего в газовой смеси.
ЛИТЕРАТУРА
1. Бесчастнов М.В., Соколов В.М. Предупреждение аварий в химических производствах. М.: Химия, 1979. 392 с.
2. Таубкин С.И., Таубкин И.С. Пожаро- и взрывоопасность пылевидных материалов и технологических процессов их переработки. М.: Химия, 1976. 264 с.
3. Абдурагимов И. М., Говоров В. Ю., Макаров В. Е. Физико-химические основы развития и тушения пожаров - М.: Высшая инженерная пожарно-техническая школа МВД СССР, 1980.
4. Баратов А.Н., Вогман Л.П. Огнетушащие порошковые составы. М.: Стройиздат, 1982.
5. Корольченко А.Я. Пожаровзрывоопасность промышленной пыли. М.: Химия, 1986. 216 с.
6. Krainov A.Yu., Dementiev A.A. Effect of inert particles on the rate of flame propagation in igniting it from a «closed end» and an «open end» // Proc. 7th International Seminar on Flame Structure 7ISFS, July 11 - July 19, 2011, Novosibirsk, Russia. P. 43.
7. Мержанов А.Г., Хайкин Б.И., Шкадинский К.Г. Установление стационарного распространения пламени при зажигании газа накаленной поверхностью // ПМТФ. 1969. № 5. С. 42-48.
8. Шкадинский К.Г., Барзыкин В.В. Закономерности зажигания газов накаленной поверхностью с учетом диффузии и гидродинамики // ФГВ. 1968. Т. 4. № 2. С. 176-181.
9. Дементьев А.А., Крайнов А.Ю. К вопросу о распространении ламинарного пламени в газе с инертной пылью // Физика горения и взрыва. 2011. Т. 47. № 4. С. 70-75.
10. Лисицын В.И., Руманов Э.Н.,Хайкин Б.И. О периоде индукции при воспламенении совокупности частиц // Физика горения и взрыва. 1971. Т. 7. № 1. С. 3-9.
11. Зельдович Я.Б., Баренблатт Г.И., Либрович В.Б., Махвиладзе Г.М. Математическая теория горения и взрыва. М.: Наука, 1980. 478 с.
12. Krainov A.Yu. Flame Propagation in the suspension of combustible particles and gases // Archivum Combustionis. 1996. V.16. No. 3-4. P. 189-198.
13. Крайнов А.Ю. Распространение пламени в смеси горючих газов и частиц // ФГВ. 2000. Т. 36. № 2. С. 3-9.
Статья поступила 12.09.2012 г.
Dementiev A.A., Krainov A.Yu. STUDYING THE INFLUENCE OF RELATIVE MOTION OF SUSPENDED INERT PARTICLES ON THE RATE OF THE GAS MIXTURE COMBUSTION FRONT. The dependence of the burning rate of a gas mixture with inert particles on concentration of fuel taken in the stoichiometric ratio with the oxidant and under the lack of fuel is studied. The compared calculation results wer obtained based on two mathematical models: with allowace for the thermal expansion and without regard to it. The dependences of the propagation velocity of the combustion front on the size and mass concentration of the inert particles were obtained at different initial concentrations of fuel.
Keywords: gas mixture, laminar flame, inert particles, the burning rate.
DEMENTIEV Aleksandr Aleksandrovich (Tomsk State University)
E-mail: daa024@ftf.tsu.ru
KRAINOV Alexey Yurievitch (Tomsk State University)
E-mail: akrainov@ftf.tsu.ru
