УДК 536.46
А.Г. Макаренко
ОСОБЕННОСТИ ФИЛЬТРАЦИОННОГО ГОРЕНИЯ ПОРИСТЫХ МАТЕРИАЛОВ
Рассматривается модель газогенератора, представляющего собой пористую пиротехническую смесь, заключенную в полузамкнутую оболочку, и работающего в режиме спутной фильтрации. Определяющие соотношения получены на основе уравнений тепломассопереноса, фильтрации газов и химической кинетики, с помощью которых рассчитаны критические значения длины заряда газогенератора и давления газа. Результаты, полученные в данной работе, легли в основу рекомендаций безопасного производства и применения пиротехнических газогенераторов.
Пористые материалы широко используются в процессах горения. Так в средствах пожаротушения широко применяются различные газогенерирующие и аэрозолеобразующие составы, используемые как их компоненты, так и самостоятельно. Эти составы представляют собой порошковые смеси с высокой пористостью. Многие пожароопасные материалы также являются пористыми (зерно, уголь, различные материалы в горючей упаковке и др.), которые способны зажигаться и гореть по фильтрационному механизму. Причем замечено, что материалы с высокой пористостью наиболее склонны к самовозгоранию [1]. Этому способствует развитая поверхность и низкая теплопроводность пористого материала. В связи с этим изучение закономерностей фильтрационного горения пористых материалов представляется весьма актуальным.
Среди средств пожаротушения весьма эффективными и надежными показали себя средства, созданные на основе газогенерирующих устройств для наддува и вытеснения дисперсных и жидких сред. Основным преимуществом подобных устройств является быстрота приведения в действие и низкая температура газообразных продуктов горения. Газогенерирующие элементы ЭГ-2 для огнетушителя ОПУ-2 выпускаются серийно уже более 5 лет и отлично зарекомендовали себя. Кроме того сами продукты горения могут успешно тушить очаги загорания. В частности, например, были проведены испытания газогенератора азота, который в течение 5 с выделял 1,5 м3 газа с температурой менее 70°С. В результате очаг пожара в замкнутом объеме был потушен. Сейчас широкое распространение получают системы объемного тушения (СОТ). Они обладают высокой огнетушащей способностью, но вместе с тем имеют высокую температуру продуктов горения, в которых, кроме того, значительное количество окислов щелочных металлов. В настоящее время имеются опытные образцы газоаэрозольного огнетушителя, продукты горения которого имеют температуру менее 100°С, а генерируемый аэрозоль содержит в своем составе N2, СО2, КС1 и ЫаС1. Подобные огнетушители могут найти применение при защите дорогостоящей оргтехники (компьютеров), складов продуктов, медпрепаратов и т.п., где даже в экстремальных ситуациях допустим контакт только с инертными материалами.
В большинстве случаев основу газогенераторов и газогенерирующих элементов составляет пористая пиротехническая смесь, горение которой происходит с фильтрацией газов.
Фильтрационным называют такое горение, при котором имеет место фильтрация газа через зону горения. Вследствие этого одним из основных факторов теплопереноса, наряду с кон-дуктивным, становится конвективный теплоперенос. Его значимость возрастает с увеличением пористости, так как рост пористости материала снижает его теплопроводность.
Различают несколько режимов фильтрационного горения: спутный и встречный, с принудительной и естественной фильтрацией, с фильтрацией химически активного и химически инертного по отношению к материалу газа. Режимы спутной и встречной фильтрации различаются направлением фильтрации газа по отношению к направлению движения фронта горения. Если они направлены в одну сторону - спутный, в противоположные - встречный. Одним из примеров фильтрационного горения с химически активным газом является горение пористого горючего при фильтрации окислителя. Если расход газа задан искусственно, то имеет место принудительная фильтрация. В естественном режиме фильтрация реализуется за счет градиентов давлений, возникающих при горении пористого материала.
Механизмы некоторых из этих режимов хорошо изучены в работах [2-5]. Так было установлено, что при спутной фильтрации, возможен такой режим горения, при котором достигается температура выше адиабатической. Это необычное явление объясняется тем, что происходит конвективный перенос тепла из зоны горячих продуктов горения в зону горения. Этот факт имеет большое значение для взрывопожаробезопасности, так как позволяет применять вещества, не способные к самопроизвольному зажиганию и горению в компактном беспорис-том состоянии.
Изучение волнового механизма фильтрационного зажигания пористого материала показало, что фильтрация нагретого газа значительно облегчает зажигание и в десятки раз сокращает время задержки зажигания по сравнению с кондуктивным нагревом беспористого материала [6]. Изучение теории фильтрационных процессов применительно к зажиганию пористых материалов позволило разработать математическую модель и экспериментально исследовать фильтрационное зажигание некоторых материалов [7]. При этом было установлено, что многие материалы, считающиеся сравнительно безопасными в компактном беспористом состоянии, могут быть весьма опасными в пористом.
Для изучения фильтрационного горения пористого материала была использована модель газогенерирующего элемента, представляющая собой пористую пиротехническую смесь, помещенную в полузамкнутую оболочку (рис.1). Смесь способна гореть без доступа воздуха и с
выделением газообразных продуктов горения. Если ее зажечь с "глухого" торца, то фронт и газообразные продукты горения будут распространяться в одном направлении - к "открытому" торцу. Реализуется режим спутной фильтрации.
Для оценки максимального давления, которое развивается в камере газогенератора при горении состава, можно решить несложную задачу, описывающую процесс газовыделения в полузамкнутой камере , т.е. когда она имеет отверстие в крышке или стенке. Здесь устанавливается стационарный режим, когда скорость газовыделения в реакто-
Рис. 1. Газогенерирующий элемент
ре т + и газоотвода из него т. сравниваются:
т+ = и0рк5+аг, т_ = Угрг8_, где щ - линейная скорость горения; - скорость истечения газа из камеры; рк -плотность состава; рг - плотность газа; 5+- площадь сечения реактора; & - площадь сечения газоотводного канала; ят-доля газообразных продуктов.
А. Л Л -
Связь аг и V имеет вид аг = V рг, где [V ] - см /г; рг -плотность при Т0 = 300 К, р=1атм.; иг определяется перепадами давлений в газоотводном канале:
Л
сір
сЫ
Т
Р+ -Р-
где р+ - давление в камере; р. - давление на выходе из канала длиной I (р=1атм); г] - динамическая вязкость газа.
Полагаемр+»р. и р+ -р. ~Р+- Тогда, приравнивая выражения для скоростей газовыделения и газоотвода, имеем
Л Р+
ЩРк$УРг = Ч~£Рг$-
и, далее, получаем выражение для давления
и„ркП5+
которое не учитывает параметры фильтрации газа, а также влияние температуры газа на максимальное давление в камере. Поэтому необходимо решать более сложную тепловую задачу, учитывающую процесс фильтрации газа.
В соответствии с теорией фильтрационного горения [2] поток фильтрующегося газа приводит к изменению температуры в зоне реакции, а следовательно, и скорости горения. Для описания этого процесса должна быть рассмотрена система уравнений тепло-массопереноса в пористой среде с нелинейными источниками, отражающая изотермическое превращение реагентов:
дГ _ _д_. дТ_ дГ_ , ш ,п
С' а ~ дхХ ас ' дк ’ ( }
% = + (2) & = Ш-Л=к0схр(-^)Г(фр\ (3)
vf= ~кг%'->ре= РАТ>Р? = А(1- (4)
Л = (А - р)' Ро>'с* = С%Р% + ср + сррр.
Начало координат (х=0) связано с фронтом горения; Ь - длина газогенератора; / - длина несгоревшей части исходной смеси.
В системе: Т - температура; Р - давление газа; // -степень(полнота) превращения вещества; V} - удельное газовыделение единицы объема исходной шихты ; и р0 - плотности газа и исходной шихты соответственно; Су - средняя теплоемкость при У=сот1; Ср - теплоемкость газообразных и конденсированных продуктов в газогенераторе при постоянном давлении; Q - теплота, выделяющаяся при сгорании 1 г исходной шихты, не зависящая от температуры; А -коэффициент теплопроводности; к0 - предэкспоненциальный множитель; Е - энергия активации; Я - универсальная газовая постоянная; ц - молекулярная масса газообразных продуктов; IV- скорость химической реакции; ку - коэффициент фильтрации.
Анализ такой системы уравнений возможен только при использовании численных методов. Однако модель можно значительно упростить, приняв во внимание особенности распространения фронта горения в условиях спутной фильтрации.
Прежде всего упростим уравнение (1). Заметим, что размеры зоны прогрева 1пр = и/а, где
А р
а а и ” линеиная скорость горения. Так как зона прогрева намного меньше длины по-
ристого элемента, то теплообменом на торцах можно пренебречь. Следовательно, его можно рассматривать как бесконечный. Преобразовав уравнение (1), связав систему координат с фронтом горения, получим
^%=&X%-c^P^v,Ь+^P^lЬu■ (5)
Интегрируя уравнение (5), находим
Тг~Т0 = Ы ■ (6)
При этом учитывается, что давление газа в точке х>0 постоянно. Это объясняется тем, что массовый газовый поток через исходные реагенты также не зависит от координаты х, если х>0.
Для вывода выражения зависимости скорости горения от основных параметров процесса был использован метод узкой зоны Зельдовича и Франк-Каменецкого. В соответствии с ним в волне горения выделены две зоны. Первая - зона реакции, где скорость химических превращений весьма высока. Вторая - зона прогрева, где тепловыделением от химической реакции можно
пренебречь, и нагрев происходит за счет тепла из первой зоны. В этом случае скорость горения имеет вид
тг сКТ2 и2=^ак0схр(--^г-)(-^-),
1 Т Т] (£
где а - коэффициент, зависящий от вида функции /(г7) = Г .сот ч и приблизительно
о
равный 1.
Для определения газового потока, фильтрующегося через пористое вещество, оцепим давление в сгоревшей смеси без истечении газа из него (замкнутая оболочка). Из экспериментов известно, что при горении выделяется 50 - 70 мл газа на 1 г исходной смеси при нормальных условиях. Использовав уравнения состояния газа, получили Р = 30 - 40 МПа, что намного больше экспериментально полученных (Рмах= 2 МПа). Из этого заключили, что наибольшее количество газа удалено из смеси. Этот факт позволил пренебречь переменным членом в уравнении фильтрации. Здесь принимается, что процесс фильтрации является квазистационарным процессом и скорость фильтрации больше скорости горения:
= -тК • (7)
Изменение давления от Ргор до Ратм происходит на участке несгоревшей смеси, где скорость реакции можно принять равной нулю. Тогда уравнение (7) имеет следующий вид:
йрУ,
-£- = о. (8)
Объединив уравнение фильтрации с уравнением (8), получаем
к «2 2
0 у —_____L-.Pl_Ел. /п\
21 ■ (9)
Распределение давления в несгоревшей части смеси описывается уравнением
^?гтг = 7- о»)
Для завершения системы уравнений используется уравнение массового баланса газообразных продуктов горения
с1М „ „ /11ч
"Л~ = О1)
где М- масса газа, приходящегося на единицу площади поперечного сечения смеси; q+ - газовыделение, приходящееся на единицу площади поперечного сечения смеси,
Я. =МгР0\ (12)
q. - отток газа, приходящийся на единицу площади поперечного сечения смеси,
к Р2 — Р2
1-=Р.',=-ф~-гГл- <13>
Длина сгоревшей части смеси составляет Ь-1, а несгоревшей /. Количества газа в этих частях соответственно составляют
М,=Р,Таш{1.-1) = -$-т, (14)
I I ,--------
М,= Шх,Тш)Ас-,Мг = \-ЛгЛ*? + Г?, (15)
о о 1 ,и
где т - пористость.
Из эксперимента известно, что давление во фронте горения больше 1 МПа. Тогда, пренебрегая Ро, получаем
, 2 т1Р»
■ М2= -ш-. (16)
Очевидно, что
^М2 _ Д>. 2т/ , , р,,( 2т , т л
л + а ~ к1т^+киЧяхт(п + кхтсот}-
Из условия стационарности следует
^ / 2т/ . ш/ \ . Л / 2т , т л_ .п ~Л.л /17л
^(звд + адои) + Лм(зд1т;п + ,А ВД( 2/ >• (17)
В результате получаем систему уравнений:
Т - Т =--—; (18)
ЛШ и сд)и-р^/
Ь(Д2-Р02) 1
..д.— 7^ (19)
м2 = ехр(-)(еКщ-~); (20)
с11/Ж=и; (21)
ЛР* ( 2т1 , гп1 \4_( 2т , т лр (р* ~ ч ,
КХТ^ + ЧЯТЫ + К?сот)Р'и-иУ'Р» ад/ 2/ (22)
Начальные условия: Т=Тп, 1=Ь, Р=Ро-
Проанализируем систему уравнений (18) - (22). Из физики процесса понятно, что для формирования квазистационарного профиля давления некоторая часть исходной смеси должна уже сгореть.
Примем, что моменту и соответствуют значения давления Р* и длины /*. Для «1* длина зоны фильтрации меньше, чем длина несгоревшей части. С увеличением / длина зоны фильтрации растет, а несгоревшей части уменьшается.
По условию стационарности вся масса газов, образовавшихся при горении, фильтруется в несгоревшую смесь. Тогда можно оценить максимальное давление:
р2 _ п2 _ ЩУАЪп
■Г*шах Г0 ~ I- ■ \
Л/у
Использовав (18) и (22) и подставив в них максимальные значения, получаем
т
Значение /* соответствует длине несгоревшего участка в момент развития максимума давления и имеет значение
'• = „,, 1/,г'~■ (24>
Р0УI + * т а х 5 В Т
Л ! •* ш
По условию квазистационарности процесса можно рассчитать время, зная скорость горения
I
и = - . (25)
и 4 ’
При условии 0<ки давление растет из-за роста скорости горения (нестационарный режим). Длина зоны фильтрации при 1>и равна длине несгоревшего участка. На этом участке давление уменьшается, так как уменьшается длина зоны фильтрации.
Теоретические результаты были подтверждены экспериментально и использованы при разработке способа и пиротехнического устройства получения газов [8].
Экспериментальные исследования показали, что в данных условиях могут самостоятельно гореть низкокалорийные системы с коэффициентом избытка окислителя а=0,1-0,3, в том числе вещества, не способные к самостоятельному горению на открытом воздухе.
Так при исследовании состава, специально разработанного для снаряжения газогенерирующих элементов, было установлено, что в обычных условиях он не способен к самостоятельному горению. Но если данный состав поместить в полузамкнутую оболочку, он сгорает спокойно. Однако в таких условиях скорость горения и уровень внутрикамерного давления
возрастают по длине пористого тела, о чем свидетельствуют профили давления, которые имеют вид ступенек с довольно крутым передним фронтом (рис.2), и график динамики газообразования (рис.З).
Сравнение этих результатов с замерами температур по длине заряда показало, что резкое повышение давления в каждой точке соответствует моменту прохождения этой точки фронтом горения. До этого, несмотря на опережающее движение части газов перед фронтом горения, давление в несгоревшей части заряда практически равно атмосферному. Позади фронта горения давление внутри камеры равномерно распределено по длине, однако его уровень и давление во фронте растут с увеличением сгоревшей части заряда.
Таким образом, можно сделать выводы.
1.Разработана модель расчета основных параметров процесса фильтрационного горения.
2.Получено выражение для расчета давления и скорости горения в зависимости от длины несгоревшей части смеси.
3.Получены выражения для максимального давления и времени достижения этого давления в процессе фильтрационного горения.
4.Проведена сравнительная оценка полученных расчетным путем данных с данными эксперимента.
5.Рассмотрена область практического применения реализуемого режима горения. Доказана возможность создания на его основе газо- и аэрозолеобразующих составов, которые используются как самостоятельные (аэрозольные и газовые
огнетушители), так и вспомогательные (для наддува порошковых огнетушителей и вытеснения дисперсных и жидких сред) средства пожаротушения
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Королъченко А.Я. Пожаровзрывоопасность промышленной пыли. М.: Химия, 1986.
2. Алдугиин А.П., Мержанов АТ. Теория фильтрационного горения- общие представления состояние исследований // Распространение тепловых волн в гетерогенных средах. Новосибирск: Наука. Сиб.отд. АН СССР, 1988.
3. Алдугиин А.П., Ивлева Т.П., Мержанов А.Г и др. Распространение фронта горения в пористых металлических образцах при фильтрации окислителя // Процессы горения в химической технологии металлургии. Черноголовка: ОИХФ АН СССР, 1975. С. 245-252.
4. Алдугиин А.П. О механизме горения СВ С-систем с газифицирующимся окислителем .// Проблемы технологического горения. Т. 1.Кинетика, термодинамика, механизм и теория горения. Черноголовка: ОИХФ АН СССР. 1981. С.11-21.
5. Алдугиин А.П, Сеплярский Б. С. Распространение волны экзотермической реакции в пористой среде при продуве газа //ДАН СССР. 1978. 241. №1. С. 72-75.
6. Сеплярский Б.С. Воспламенение коцденсированных систем при фильтрации газа // Физика горения и взрыва. 1991. №1,.С.З-12.
МПа
1-3 , им
Р и с. 2. Профили давлений внутри камеры. Время перемещения фронта горения: 1 - 2,0; 2 - 2,7; 3 - 3,2 с
1.0 2,0 1г.с
3.0
Р и с. 3. Динамика газообразования. Длина заряда: 1 - 30; 2 - 60; 3 - 90 мм