Математическая модель критической продолжительности пожара в помещении, смежном с очагом пожара

УДК 614.841

Математическая модель критической продолжительности пожара в помещении, смежном с очагом пожара

С.С. Лапшин, преп.

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

Ключевые слова: моделирование, пожар, температура.

Mathematical Model of Critical Fire Duration in Adjacent Room with Fire Source

S.S. Lapshin, Teacher

The article presents generalized equation allowing to define critical duration of fire in adjacent room. The decision is based on the integral model of fire. There is validation of the model with experimental and theoretical data. The equation is obtained with regard to not simplified solution of mass equation in a fire room.

Key words: modeling, fire, temperature.

Для решения проблемы повышения уровня безопасности людей на пожарах требуется разработка конструктивных и объемнопланировочных решений в строительстве с учетом динамики опасных факторов пожара (ОФП) [1].Современные методы исследования динамики ОФП в помещении базируются на различных математических моделях пожара [2].

Время от возникновения пожара до наступления ОФП называется критической продолжительностью пожара (КПП). Актуальность исследования динамики ОФП в смежном с очагом пожара помещении обусловлена отсутствием в нормативных документах аналитических решений (формул) для определения критической продолжительности пожара. В настоящее время для расчета динамики ОФП в смежном с очагом пожара помещении по интегральной модели применяют численные методы. Этот подход связан с большим количеством вычислений и может быть реализован только с помощью компьютера. В программу подготовки инженера пожарной безопасности не входит изучение численных методов, в результате чего провести такой расчет самостоятельно ему будет не по силам. При использовании численных методов дифференциальные уравнения заменяются их дискретными аналогами, следовательно, нет возможности установить функциональные зависимости ОФП от условий протекания пожара в помещении. Результаты численного расчета заведомо трудно проверить, так как программы, позволяющие производить такие расчеты, в большинстве своем являются интеллектуальной собственностью их разработчиков и доступ к их исходному коду, а следовательно, к уравнениям и алгоритмам, которые они используют, закрыт [3, 4]. Специалистам, производящим расчеты по таким программам, остается полагаться на то, что руководства по данным продуктам полностью отражают все нюансы расчета, и на отсутствие ошибок в программном коде. Как правило,

для верификации таких расчетов прибегают к аналитическим методам.

Применение для таких задач полевого метода наталкивается на еще большее количество проблем. Необходимость высококвалифицированных специалистов, наличие у них большого опыта в этой области, огромное количество вычислений (машинного времени), ряд нерешенных задач газовой динамики (турбулентность), а также мнения некоторых исследователей, что полевая модель до сих пор остается «вещью в себе» [1], - все это позволяет сделать вывод, что использование полевых моделей для расчета динамики ОФП является своего рода искусством. Как правило, полевые модели применяются для расчета ОФП в помещениях сложной конфигурации, в которых невозможно использовать интегральные модели.

Ниже, используя интегральный термодинамический метод описания пожара в помещении, созданный Заслуженным деятелем науки РФ, доктором технических наук, профессором Ю.А. Кошмаровым, в общем виде получена формула для определения динамики ОФП в помещении, смежном с очагом пожара.

Рассмотрим интегральную математическую модель начальной стадии пожара (НСП). Руководствуясь работой [2], примем следующие допущения: газовая среда внутри помещения при пожаре является смесью идеальных газов; в каждой точке пространства внутри помещения в любой момент времени реализуется локальное равновесие; режим газообмена является односторонним: поступление воздуха в помещение из окружающей среды отсутствует, а продукты горения выходят через все имеющиеся проемы в смежное помещение; среднее давление среды в помещении остается практически постоянным, равным давлению наружного воздуха; отношение теплового потока в ограждения к тепловыделению есть величина постоянная,

равная своему среднему значению на рассматриваемом интервале времени.

Динамика ОФП в помещении, где расположен очаг пожара. Запишем систему уравнений интегральной математической модели НСП с учетом принятых допущений:

- уравнение материального баланса

VI

= \-гн - вг;

- уравнение энергетического баланса

п\0% 0 - Ф1) - сРт1твг = °;

- уравнение баланса кислорода

V1

^102 . ^102

1 —= -Л\-о2 - ~вг

d т

р1т

ег;

- уравнение баланса токсичных продук-

тов горения

^ ^1ПГ

d т

= Л\-пг -^ПСв,

р1т

Г >

- уравнение баланса оптической плотности дыма

- уравнение состояния газовой среды

Р1т = р1тЯТ1т ,

где V - объем помещения, м3;

р1т, ^ю2, ^1ПГ = {^100, ^1С02 }, Н-1 - среднеобъемные

значения плотности газовой среды, парциальных плотностей компонентов газовой смеси и оптической плотности дыма, кгм-3; у - скорость выгорания горючей нагрузки (ГН), кгс-1; —

расход газов, уходящих из помещения через проемы, кгс-1; п - коэффициент полноты сгорания ГН; 0% - теплота сгорания ГН, Джкг-1; ф -

коэффициент теплопоглощения; Ср - теплоемкость среды при постоянном давлении, Джкг-1К'1; Т1т - среднеобъемная температура газовой среды в помещении, К; -0 - количество кислорода, необходимое для сгорания единицы массы ГН; 1~ПГ - количество токсичных газов, образующихся при сгорании единицы массы ГН; й -дымообразующая способность ГН, Нпм2кг-1; Р1т - среднеобъемное давление газовой среды в помещении, Па; Я - газовая постоянная, 8,31 кДжкг-1К-1 .

Уравнения (3)-(5) можно обобщить:

dт р1т

где

?1/ = {^102, ^1ПГ, ^1ДП, Н-1},

-/ = {-гн ; -п^02; п-С02; п-С0; п-ДП; й}.

Тогда, чтобы решить систему уравнений (1)-(5), необходимо и достаточно решить уравнение (7).

Исходя из принятых допущений и используя уравнение состояния (6), можно получить следующее соотношение [2]:

р0Т0 = р1тТ1т .

Решая уравнение (7) совместно с (2), получим

где С0/ = {р1т,^°02,^1002,^100,^0ДП,^°} значения

плотности газовой среды, парциальных плотностей компонентов газовой смеси и оптической плотности дыма в момент времени т = 0;

?1/ = {р1т,^102,^1С02,^100,^1ДП,^1} - среднеобъем-

ные значения плотности газовой среды, парциальных плотностей компонентов газовой смеси и оптической плотности дыма на отрезке времени

т = [т0, тКР ] ; -/ ={-гН , -п-02, л-С02, л-С0 , п-ДП,й} -

количество вещества, выделившегося в результате термической деструкции ГН -гн и реакции горе-

Лор

ния; 7=

СрТ0р0

Для определения КПП с учетом уровня рабочей зоны используем формулу Меркушки-ной-Зотова-Тимошенко:

^ ехр(-1,4У I,

рз -50 у Ч 2Л 1

У

где 5рз - предельно допустимое значение

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

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

тРз =

1кр

А|п (1 -в

аА

где А и п - параметры, зависящие от условий распространения горения по твердой ГН [2, с. 62];

Р/ =

^ (рКт -р0) ?1* - р0

для расчета температуры;

гРЗ г0

р ?1у -?1'

р/ =------- для расчета парциальных плотно-

?1* -?1 ]

стей газовой среды.

Динамика ОФП в помещении, смежном с очагом пожара. Запишем систему уравнений, при этом, как и в работе [5], сделаем допущение, что режим газообмена является односторонним, т.е. продукты горения из помещения очага пожара через общий проем выходят в смежное помещение:

- уравнение материального баланса

дР2т _

дт

_ пЄі — ©2;

- уравнение энергетического баланса

СрТ1тп©1

- уравнение баланса кислорода

V,

дЯ20 _ п ^10

©1 —

^20

д т р1т р2т

- уравнение баланса токсичных продуктов горения

'___ Я________

Ц2 ^2ТС _ п кгвв-! — ^02;

д т р1т р2т

- уравнение баланса оптической плотности дыма

р1т р2т

- уравнение состояния газовой среды

Р2т = р2тЯТ2т .

Здесь в1 - расход газов из помещения очага пожара в смежное с ним, кгс-1; в2 - расход газов из смежного с очагом пожара помещения в атмосферу, кгс-1; р/т, |/0, |/ге, ц - среднеобъемные значения плотности газовой среды, парциальных плотностей кислорода и токсичных продуктов горения, оптической плотности

дыма, кг м-3; Ср - теплоемкость среды при по-1 -1

стоянном давлении, Джкг "К-1; Т1т - среднеобъемная температура газовой среды в помещении, К; ф2 - коэффициент теплопоглощения; V/ -объем помещения, м3; п = Р^+ ^Рщ) - коэффициент проемности, характеризующий газообмен помещения очага пожара и смежного с ним [6]; Рй - площадь проема между смежными помещениями, м2; ^Рщ - сумма площадей от-

2

крытых проемов в помещении очага пожара, м ; т - время, с; / = 1, 2 (индекс 1 обозначает помещение очага пожара, индекс 2 - помещение, смежное с очагом пожара).

Уравнения (14)-(18) можно записать в обобщенном виде:

2 _ п в1 — в2,

2 д т - 1 - 2

р1т р2т

где С2/ ={р2т, ^20, %2Т0, ^2} - среднеобъемные значения плотности газовой среды, парциальных плотностей токсичных компонентов газовой смеси, оптической плотности дыма на отрезке времени т = [т0, тКР] .

Решая уравнение (20), получим

С2/ _

С*

1 — Ф2

с0—с*

1 —ф2 — к +С°ехр {—І2 (1 — Ф2 )),

(1 — ехр ( (1 — ф2 ))) +

(ехр (—72к) — ехр ( (1 — ф2 ))) -

V,

.Гудт, к _ -Ц2-

Г , пЦ1

где

г = СрТ0-тй г0 7 = па

^ = (1 -ф1) Ор , 2= V2

Для определения КПП рассчитываем критические значения среднеобъемных параметров пожара с учетом высоты рабочей зоны в смежном помещении:

сКР _ с0

Ь/ч Ъ/1

с Р3 — с0

ч ч

У-ехр ( 1,4-У-2Ь Ч 2Ь

Критическая масса ГН есть масса ГН, при сгорании которой в рабочей зоне смежного помещения параметры пожара достигнут своих опасных для человека значений

г кр

^ 2 Ч

С1*

1-

Ф2

1 — ехр

™МКР

с0, — С1 (

1 — ф2 — к

па

V.

+Сі0,ехр ^Цамкр

Уравнение (23) решается методом последовательных приближений. ОФП, наступающий первым, характеризуется наименьшей, по сравнению с другими, величиной Мкр.

Проверка адекватности модели проведена с помощью экспериментальных данных, полученных в работе [7], и теоретических расчетов работы [6]. Эксперимент был проведен на модели клуба-театра на 600 мест, выполненной в 1/6 натуральной величины и состоящей из двух помещений: сценического и зрительного зала. Для проверки теории использованы данные, полученные в двух сериях экспериментов. Первая проводилась при закрытых проемах в сценическом помещении (опыт 7). Вторая - при открытых (опыт 1). Проем между смежными помещениями был открыт во всех опытах. Более подробно описание эксперимента приведено в работе [7]. На рис. 1-3 приведены графики изменения среднеобъемной температуры, где модель 1 - значения среднеобъемной температуры, полученные с помощью уравнений, представленных выше; модель 2 - с помощью уравнений в работе [6]. Максимальные расхождения составили: модели 2 и эксперимента - 4,9 %; модели 1 и эксперимента - 3 %; модели 1 и модели 2 - 4,5 %.

360 350

340 --

330

К 320

,а р 310

у т 3 о о _ _

а р290

е п 280

м е 2 -ч о

н 260 --

250

30

60

Эксперимент - Модель 1 Модель 2

90 120

Время, с

Рис. 1. Изменение среднеобъемной температуры

360

а

р

у

т

а

р

е

30 60

-0 —Эксперимент

------Модель 1

-А-----Модель 2

90

120

Время, с

Рис. 2. Изменение среднеобъемной температуры

0 30 60

--- — Эксперимент

--------Модель 1

--------Модель 2

90 120

Время, с

Заключение

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

Полученное обобщенное решение позволяет находить среднеобъемные значения ОФП в смежном с очагом пожара помещении. Отличие температурных кривых модели 1 и модели 2 объясняется тем, что при решении уравнений модели 1 использовалось обобщенное решение уравнений интегральной модели для очага пожара, которое получено в работе [4] без упрощения балансового уравнения массы, принятого в работе [2].

Для оценки влияния этого упрощения на точность расчета ОФП, уточнения коэффициентов проемности п и теплопоглощения ф требуется проведение дальнейших теоретических и экспериментальных исследований.

Список литературы

1.Кошмаров М.Ю. Моделирование динамики начальной стадии пожара в помещениях, зданиях и сооружениях при воспламенении горючей жидкости: дис... канд. техн. наук. - М.: Академия ГПС МЧС России, 2004.

2.Кошмаров Ю.А. Прогнозирование опасных факторов пожара в помещении: Учеб. пособие. - М.: Академия ГПС МВД России, 2000.

3.Deibjerg T., Husted B.P., Bygbjerg H., Westerman

D. Argos User's Guide / Danish Institute of Fire and Security Technology, 2003.

4.Wade C.A. BRANZFIRE Technical Reference Guide BRANZ Study Report No 92. Building Research Association of New Zealand, Judgeford, 2004.

5.Лапшин С.С., Тараканов Д.В. Обобщенное решение системы уравнений начальной стадии пожара в помещении // Вестник Ивановского института ГПС МЧС России. -2008. - №1. - С. 25-28.

6. Коршунов И.В., Кошмаров М.Ю. Математическая модель начальной стадии пожара в театре с колосниковой сценой. Ч. II: Экспериментальная проверка математической модели // Пожаровзрывобезопасность. - 2006. - Т.15. - № 2. -С. 17-23.

7.Кривошеев И.Н. Исследование начальной стадии развития пожара в зрелищных предприятиях (с целью обоснования допустимого времени эвакуации): дис... канд. техн. наук. - М.: ВИПТШ МВД СССР, 1973.

Рис. 3. Изменение среднеобъемной температуры

0

Лапшин Сергей Сергеевич,

Ивановский институт ГПС МЧС России, преподаватель кафедры гражданской защиты, e-mail: wfxdfx@gmail.com