6.1 Динамика опасных факторов пожара в протяженных помещениях

Двухмерная плоская математическая модель с учетом процессов горения использована для воспроизведения температурных, скоростных и концентрационных параметров, проведенных в Японии огневых экспериментов по сжиганию гептана в модели туннеля, имеющей высоту потолка 1,6 м, площадь

речного сечения 2,01 м ² и длину 21 м. Горение 18 л гептана инициировалось в стальном поддоне площадью 0,25 м ², расположенного на полу в средней части туннеля, открытого с обеих сторон, и продолжалось около 900 с. Как показали проведенные экспериментальные исследования, в установившейся стадии температура на выходе из туннеля достигала 350 ^оС, а скорость потока уходящего горячего газа увеличивалась до 3 м/с при скорости втекающего извне воздуха 0,5 м/с. Оптическая плотность дыма на выходе газов из туннеля составляла 3,33 1/м.

Рисунок 6.1 – Расчетные поля температур (а), оптической плотности газов (б )и концентрации СО (в) при горении гептана в центре туннеля длиной 21 м

Расчетная температура и оптическая плотность продуктов горения (рисунки 6.1.а, 6.1.б) близки к экспериментальным. А на рисунке 6.1.в представлены поля концентраций окиси углерода для того же момента 480 с. Согласуются и расчетные данные скоростей на входе и выходе из туннеля с экспериментальными данными.

Приведенные данные по составу продуктов горения позволили при записи реакции горения получить массовые стехиометрические коэффициенты (таблица 6.1).

Таблица 6.1 - Значения массовых стехиометрических коэффициентов для различных веществ

Вещество	Массовый стехиометрический коэффициент
	3,205
	2,611
	1,440
	0,0442
	0,110

6.2 Моделирование пожаров на складах лесоматериалов

Пожары, возникающие на предприятиях деревообработки, приносят значительный материальный ущерб, нередко они переходят на окружающую жилую застройку. Наиболее опасны пожары на открытых складах лесоматериалов. В этом случае они имеют классический вид массового пожара. Для таких пожаров характерны: разброс искр и головней в радиусе до 300 м, а при штормовой скорости ветра – в радиусе более 1км, большая скорость распространения пламени по штабелям.

В нашей стране проблемой тушения пожаров на складах лесоматериалов начали заниматься еще в 40–50 годы прошлого века. В эти годы изучение пожаров происходило в основном экспериментально на натуральных объектах. Древесина укладывалась в штабеля высотой до 6 м.

Установлено, что скорость распространения огня по штабелям зависит от многих факторов, в первую очередь, от вида и способа укладки древесины, ее влажности, направления ветра.

Проводимые опыты были чрезвычайно трудоемкими, поэтому их количество было ограничено, а полученные результаты практически не воспроизводимы.

Учитывая это обстоятельство, была предпринята попытка изучить процессы горения штабелей на моделях. Одним из первых во ВНИИПО такие исследования провел П.П. Девлишев.

При изучении горения штабелей на моделях автор попытался установить влияние геометрических размеров штабелей и отдельных элементов (брусков, досок), плотности распределения горючего в объеме на физические параметры горения: длительность, скорость распространения огня, скорость выгорания древесины, размеры факела пламени и др. Хотя в результате проделанной работы и не удалось получить критериальных уравнений, тем не менее, были установлены некоторые важные закономерности горения штабелей.

По экспериментальным данным, полученным П.П. Девлишевым, Н.А.Иванов, с использованием метода анализа размерностей получил эмпирические зависимости для скорости распространения фронта пламени по штабелям:

- при отсутствии ветра , м/с:

; (6.2)

- при ветре , м/с:

; (6.3)

где - скорость выгорания древесины, кг/(м² с);

- плотность укладки древесины, кг/м³;

- влагосодержание древесины;

- величина поверхности древесины в единице объема штабеля;

- высота штабеля, м;

- длина фронта пламени, м;

- скорость ветра, м/с.

Эти зависимости непригодны для расчета скорости распространения фронта пламени по территории лесосклада при наличии противопожарных разрывов.

С конца 70-х годов ХХ в. во ВНИИПО начали развиваться методы математического моделирования развития пожара на складах лесоматериалов.

В моделях распространения пожара теплоизлучением принимается следующая схема лучистого теплообмена: излучающая поверхность (как правило, плоскость) – одна из боковых граней штабеля.

Эта поверхность увеличивается на высоту факела пламени, формирующегося над штабелем; в зависимости от ориентации облучаемой и излучающей поверхностей по известным соотношениям рассчитываются коэффициенты теплообмена; степень черноты пламени разными авторами принимается от 0,6 до 1,0; температура факела варьируется в пределах от 900 до 1200 ^оС; за критерий воспламенения принимается либо соответствующее значение падающего теплового потока и время его действия, либо достижение на поверхности древесины температуры ее воспламенения.

Как показали эксперименты, при моделировании распространения пожара теплоизлучением важно учитывать влияние ветра на положение факела в пространстве. Под действием ветра факел отклоняется от вертикального положения над крышей штабеля, и тепловой поток с подветренной стороны от очага пожара увеличивается.

На рисунке 6.2 приведен график изменения скорости распространения пожара в зависимости от скорости ветра и направления движения фронта огня. Аналогичные зависимости были получены для определения периметра и площади пожара. Представленные фрагменты математической модели дают возможность учесть факторы, влияющие на процесс развития пожара.

Рисунок 6.2 – Скорость распространения фронта пламени для склада круглого леса: 1 - 1 м·с ^–1 ; 2 - 4 м·с ^–1; 3 - 7 м·с ^–1; 4 - 10 м·с ^–1; 5 - 13 м·с ^–1; 6 – 16 м·с ^–1 ; 7- 19 м·с ^–1

При расчете тепловых потоков факел пламени разбивается на две плоские прямоугольные излучающие системы:

- - боковая поверхность штабеля;

- - факел над крышей штабеля, наклоненный ветром на угол .

Таким образом, были установлены зависимости величин потоков не только от характеристик горящей древесины, но и от параметров атмосферы.

Угол наклона конвективной колонки вычисляют по формуле:

, (6.4)

где - вертикальная координата оси конвективной колонки, формирующейся над горячими штабелями, м.

Вертикальная координата оси конвективной колонки вычисляют по формуле:

, (6.5)

где - параметр плавучести;

- горизонтальная координата оси конвективной колонки, м.

Параметр плавучести вычисляют по формуле:

, (6.6)

где - низшая теплота сгорания древесины, Дж/кг;

- скорость выгорания древесины, кг/(м² с);

- площадь пожара, м²;

- температура воздуха, ˚С;

- скорость ветра, м/с.