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

Расчет динамики опасных факторов пожара в помещении с использованием интегральной математической модели пожара

Определение критической продолжительности пожара и времени блокирования эвакуационных путей

Прогнозирование обстановки на пожаре к моменту прибытия первых

Подразделений на тушение

Расчет огнестойкости ограждающих строительных конструкций

С учетом параметров реального пожара

Расчет динамики опасных факторов пожара в помещении с использованием зонной математической модели пожара

Заключение

Литература

Введение

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

-при создании и совершенствовании систем сигнализации и автоматических систем пожаротушения;

-при разработке оперативных планов тушения пожаров;

-при оценке фактических пределов огнестойкости;

И для многих других целей.

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

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

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

Исходные данные

Краткая характеристика объекта

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

Размеры цеха в плане:

- ширина =36 м;

- длина = 18 м;

- высота = 6м.

План цеха показан на рис.п.1.1

Рис. п.1.1. План отделочного цеха мебельного комбината

В наружных стенах помещения цеха имеется 3 одинаковых оконных проема, один из которых открытый. Расстояние от пола до нижнего края каждого оконного проема = 0,8 м. Высота оконных проемов = 2,4 м. Ширина каждого оконного проема = 6,0 м. Остекление оконных проемов выполнено из обычного стекла. Остекление разрушается при среднеобъемной температуре газовой среды в помещении, равной 300 0 C.

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

Отделочный цех имеет два одинаковых дверных проема, соединяющих цех с наружной средой. Их ширина равна 0,9 м и высота 2 м. При пожаре дверные проемы открыты.

Полы цеха бетонные, с асфальтовым покрытием.

Горючий материал представляет собой деревянные детали мебели, покрытые лаком. Горючий материал расположен на полу. Размер площадки, занятой горючим материалом: длина – 20 м, ширина – 10 м. Количество горючего материала составляет 10 тонн.

Сбор исходных данных

Геометрические характеристики объекта.

Выбирается положение центра ортогональной системы координат в левом нижнем углу помещения на плане (рис. п.1.1). Координатная ось x направлена вдоль длины помещения, ось y - вдоль его ширины, ось z - вертикально вдоль высоты помещения.

Геометрические характеристики:

помещение: длина L =36 м; ширина В = 18 м; высота Н = 6 м.

двери(количество дверей N д o =2): высота h д1,2 = 2,0 м; ширина b д1,2 = 0,9 м; координаты левого нижнего угла двери: у д1 = 10 м; х д1 = 0,0 м; у д2 = 7 м; х д2 = 36,0 м;

открытые окна (количество открытых окон N о o = 1): высота h о o 1 = 2,4 м; ширина b о o 1 = 6,0 м; координаты одного нижнего угла окна: x о o 1 = 3,0 м; у о o 1 = 0 м; z о o 1 = 0,8 м;

закрытые окна (количество закрытых окон N з o =2): высота h з o 1,2 = 2,4 м; ширина b з o 1,2 = 6,0 м; координаты одного нижнего угла окна: x з o 1 = 15 м; y з o 1 = 0,0 м; z T кр = 300 о С; x з o 2 = 27 м; y з o 1 = 0,0 м; z зо1 = 0,8 м; температура разрушения остекления T кр = 300 о С;

технологический проем(количество проемов N п o =1): высота h п1 = 3,0 м; ширина b п1 = 3,0 м; координаты левого нижнего угла проема: у п1 = 18 м; х п1 = 20,0 м.

ЗОННАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПОЖАРА В ПОМЕЩЕНИИ. Конвективная колонка Припотолочный слой Выводы по лекции Цели лекции: Учебные В результате прослушивания материала слушатели должны знать: опасные факторы пожара воздействующие на людей на конструкции и оборудование предельно допустимые значения ОФП методы прогнозирования ОФП Уметь: прогнозировать обстановку на пожаре. Прогнозирование опасных факторов пожара в помещении.

Поделитесь работой в социальных сетях

Если эта работа Вам не подошла внизу страницы есть список похожих работ. Так же Вы можете воспользоваться кнопкой поиск

ЛЕКЦИЯ

по дисциплине "Прогнозирование опасных факторов пожара"

Тема №6. «ЗОННАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПОЖАРА В ПОМЕЩЕНИИ. ЧИСЛЕННАЯ РЕАЛИЗАЦИИ ЗОННОЙ МОДЕЛИ»

План лекции:

Введение

Конвективная колонка
Припотолочный слой

Выводы по лекции

Цели лекции:

Учебные

В результате прослушивания материала слушатели должны знать:

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

Уметь: прогнозировать обстановку на пожаре.

Развивающие:

выделять самое главное
самостоятельность и гибкости мышления
развитие познавательного мышления

Литература

Д.М. Рожков. Прогнозирование опасных факторов пожара в помещении. Иркутск 2007 С.89
Ю.А.Кошмаров, М.П. Башкирцев Термодинамика и теплопередача в пожарном деле. ВИПТШ МВД СССР, М., 1987 г.
Лабораторный практикум «Прогнозирование опасных факторов пожара». Ю.А.Кошмаров, Ю.С.Зотов. 1997 г.
Ю.А.Кошмаров, В.В. Рубцов, Процессы нарастания опасных факторов пожара в производственных помещениях и расчет критической продолжительности пожара. МИПБ МВД России, М., 1999 г.

Введение

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

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

1. Постановка задачи о зонном моделировании.

В соответствии с вышесказанным в объеме помещения можно выделить три характерные зоны: конвективную колонку над очагом пожара, припотолочный слой нагретого газа и воздушную зону с практически неизменными параметрами состояния, равными своим начальным значениям. Математическая модель пожара, базирующаяся на разбиении пространства на характерные области, получила название трехзонной модели. Схема этой модели показана на рис. 6.1. На этой схеме использованы следующие обозначения: у к - координата нижней границы припотолочного слоя, отсчитываемая от поверхности горения; у ДВ - высота дверного проема; d э - эквивалентный диаметр очага горения; 2 h - высота помещения; G K - поток газа, поступающего в припотолочный слой из конвективной колонки, кг·с -1 ; G B - поток воздуха, поступающий в колонку из зоны III , кг·с -1 ;. G Г - поток вытесняемого газа из помещения, кг·с -1 ; ψ - скорость выгорания, кг·с -1 ; δ - расстояние от пола до поверхности горения, м.

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

При второй фазе нижняя граница II зоны расположена ниже верхнего края дверного проема. С наступлением второй фазы начинается процесс истечения нагретых газов из помещения через дверной проем. До наступления этой фазы имеет место лишь вытеснение (через дверной проем) холодного воздуха из III зоны.

Рис. 6.1. Схема трехзонной модели пожара:

I зона конвективной струи (конвективная колонка);

II - зона припотолочного нагретого газа; III - зона холодного

воздуха; IV - зона наружного воздуха (наружная атмосфера)

2. Конвективная колонка

Рассмотрим прежде всего I зону. Теория свободной конвективной струи к настоящему времени весьма детально разработана. Эта теория является одним из разделов вязкой аэродинамики газов. Она позволяет рассчитывать поля температур, плотностей и скоростей в конвективной колонке. Для определения температур и массовых расходов в сечениях конвективной колонки можно использовать формулы:

(6.1)

(6.2)

где Q пож - скорость тепловыделения, Вт; Q p H низшая теплота сгорания, Дж·кг -1 ; ψ уд - удельная скорость выгорания, кг·м -2 ·с -1 ; g - ускорение свободного падения, м·с -2 ; Т о и ρ 0 - температура и плотность холодного (окружающего) воздуха; G - расход газов через сечение струи, отстоящее от поверхности горения на расстояние у, кг·с -1 ; с р - изобарная теплоемкость газа, Дж·кг -1 ·К -1 ; - доля, приходящаяся на поступающую в ограждение теплоту от выделившейся в очаге горения; у - координата сечения колонки, отсчитываемая от поверхности горения, м; у 0 - расстояние от фиктивного источника тепла до поверхности горения, м.

С помощью формул (6.1) и (6.2) можно рассчитать расход газа из I зоны, поступающего во II зону, и его температуру. Для этого нужно положить координату у в формулах (6.1) и (6.2) равной координате нижней границы припотолочного слоя у к .

Расстояние от фиктивного источника тепла до поверхности горения вычисляется по формуле:

(6.3)

где F Г - площадь пожара, м 2 .

3. Припотолочный слой

Рассмотрим теперь II зону (припотолочный слой нагретых газов). Объем этой зоны в момент времени τ равен

где F П0 T - площадь потолка; у к - координата нижнего края припотолочного слоя газов. Масса газа, заключенная во II зоне, составляет величину т 2 = р 2 V 2 Давление в зоне II практически не меняется и остается равным начальному значению, т.е. Р 0 . Внутренняя (тепловая) энергия II зоны составляет:

Запишем уравнения материального баланса и энергии для II зоны применительно к первой фазе начальной стадии пожара:

(6.4)

(6.5)

где ρ 2 - средняя плотность во II зоне; Т 2 - средняя температура во II зоне; Q w 2 - тепловой поток от припотолочного слоя газа в ограждения, кВт.

Параметры состояния Т 2 и ρ 2 связаны между собой следующим уравнением:

(6.6)

Уравнение (6.6) следует из условия равенства давлений во всех зонах. Это условие является приближенным, но применимым для реальных пожаров.

Преобразуем уравнение энергии (6.5), используя уравнение (6.6):

или

и окончательно (6.7)

Из уравнения (6.1) следует:

(6.8)

Подставляя формулу (6.8) в уравнение (6.7), получим:

Примем, что (для начальной стадии φ= 0,66 ).

После дальнейших преобразований получим следующее уравнение:

(6.8а)

Подставим в это уравнение выражение для G k (6.2):

(6.9)

Отметим, что в этом уравнении

Введем обозначения:

Функции β(τ) и γ(τ) при горении твердых ГМ в момент времени τ = 0 равны нулю, так как F Г → 0. Уравнение (6.9) принимает вид:

(6.10)

Начальное условие.

Решение уравнения (6.10) при заданном начальном условии будем искать для интервала времени от τ = 0 до τ * , где τ * - момент окончания первой фазы начальной стадии пожара. После того как найдена функция у к (τ), находим G k = f 1 (τ) ; V 2 = f 2 (τ).

Преобразуем уравнение материального баланса (6.4). Интегрируя его, получаем:

(6.11)

После преобразований из формулы (6.11) получаем:

(6.12)

После вычислений плотности ρ 2 определяется средняя температура в припотолочном слое газа:

(6.13)

Уравнение баланса для токсичного газа (продукт горения) во II зоне имеет вид:

(6.14)

где ρ n - парциальная плотность токсичного газа; L - количество (масса) токсичного газа, образующаяся при сгорании 1 кг горючего материала. Из формулы (6.14) следует формула:

(6.15)

где М τ - количество (масса) ГМ, выгоревшего к моменту времени τ.

Уравнение дыма для припотолочного слоя имеет вид:

и, следовательно:

Исходя из выше изложенного, имеем уравнение с разделяющимися переменными, с помощью которого рассчитывается изменение координаты границы припотолочного слоя в течение времени:

где:

при условии: y 0 = const ;

Выводы по лекции : зонная модель представляет собой опять же частный случай интегральной модели для припотолочного слоя, и с применением известных теорий, в частности - теории конвективной колонки.

PAGE 6

Другие похожие работы, которые могут вас заинтересовать.вшм>
10172.		Основные понятия и уравнения интегральной математической модели пожара в помещении	53.24 KB
	Основные понятия и уравнения интегральной математической модели пожара в помещении. Основные понятия математической модели пожара в помещении. Допущения интегрального метода термодинамического анализа пожара.
10170.		ГАЗООБМЕН ПОМЕЩЕНИИ И ТЕПЛОФИЗИЧЕСКИЕ ФУНКЦИИ, НЕОБХОДИМЫЕ ДЛЯ ОПИСАНИЯ ЗАМКНУТОГО ПОЖАРА	576.18 KB
	Распределение давлений по высоте помещения. Плоскость равных давлений и режимы работы проема. Распределение перепадов давлений по высоте помещения. Побудителем движения газа через проемы является перепад давлений т.
10173.		МАТЕМАТИЧЕСКАЯ ПОСТАНОВКА ЗАДАЧИ О ДИНАМИКИ ОПАСНЫХ ФАКТОРОВ НАЧАЛЬНОЙ СТАДИИ ПОЖАРА	101.99 KB
	Вопросы обеспечения безопасности людей, зданий и сооружений сегодня являются приоритетными. При этом, наиболее актуальными остаются вопросы, связанные с обеспечением пожарной безопасности. Наряду с огромным материальным ущербом, пожары продолжают уносить жизни людей.
7866.		Экономико-математическая модель создания МТЛЦ	16.16 KB
	Следует иметь в виду что отправитель продукции не всегда отдает предпочтение наиболее дешевому варианту по тарифам и прочим платежам перевозчику и экспедиторам. Обобщая вышеприведенные рассуждения можно сделать вывод что в процессе выбора транспортнотехнологических систем доставки продукции должны учитываться разносторонние интересы клиентов и различных видов транспорта. Товарооборот между продавцом и покупателем рассматриваемой продукции О будет уменьшаться а объем национального продукта также сократится Н. В такой ситуации...
1538.		Математическая модель диска с изгибающими нагрузками	1.12 MB
	Множество алгоритмов математического программирования, решающих задачи оптимального проектирования, реализовано в виде программных библиотек или в качестве части пакетов универсальных программных комплексов. Общим недостатком этих алгоритмов является низкая скорость сходимости и высокая вероятность получить неоптимальный результат.
16733.		МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЭКОНОМИЧЕСКОГО ЦИКЛА ЖЮГЛЯРА	726.28 KB
	В частности сокращение спроса ведет к сокращению производства а сокращение производства в свою очередь ведет к дальнейшему сокращению спроса; наличие определенной инерционности запаздывания реакции экономики на изменение условий например запаздывание в изменениях уровня инвестиций по отношению к изменению спроса; усиление финансовой системой обратных положительных связей и временных лагов в экономике за счет влияния на процессы кредитов спекулятивных операций и т. Положительная обратная связь между инвестициями и изменением...
5810.		Экономико-математическая модель по оптимизации производственной структуры в ООО «Пшеница»	77.63 KB
	Экономикоматематическая модель даёт возможность определить основные параметры развития производства для текущего и перспективного планирования может использоваться для анализа сложившейся структуры производства позволяющего выявить более целесообразные пути использования ресурсов и возможности увеличения объёмов производства продукции опираясь на фактические данные за предшествующие годы. Под оптимальной производственной структурой сельскохозяйственного предприятия следует понимать такие количественные соотношения между отдельными...
21763.		Математическая модель системы автоматического регулирования высоты жидкости в герметизированной емкости	3.32 MB
	Но магистральная линия создания принципиально новых и совершенствования существующих технических устройств - это реализация возможностей, открывающихся при использовании результатов фундаментальных исследований. Этим, в частности, объясняется и современный акцент в инженерном образовании на фундаментальную научную подготовку. Решающую роль при реализации результатов таких исследований играет математическое моделирование.
3211.		Математическая модель непропорционального (эксцедентного) перестрахования. Общая схема. Численный пример	67.57 KB
	Непропорциональное страхование – или Страхование эксцедента убытка (stop-loss, передается то, что выше опред. суммы, котор.зависит от r). Перестрахование редко вступает в действие, но в этих случаях не несёт рисков – распределение убытка несимметрично.
12153.		Математическая модель межрайонных корреспонденций (передвижений) на индивидуальном транспорте в условиях высокого уровня автомобилизации	17.78 KB
	Это требует учета при моделировании корреспонденций ограниченных возможностей районов по размещению прибывающих автомобилей. Величина этих затрат должна зависеть от соотношения объема автомобильных прибытий в район и его возможностей принять такой объем и увеличиваться по мере нарастания объема прибытий. в районах прибытия для автомобильного транспорта возникают такие дополнительные затраты при которых в ходе расслоения корреспонденций формируются объемы прибытий на автомобилях порождающие именно эти величины дополнительных затрат. Для...

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

Рис. 1.2. Схема зонной модели пожара в помещении

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

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

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

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

Ниже сформулированы основные уравнения зонной математической модели пожара в помещении.

Уравнение баланса массы . При отсутствии проемов в верхней части помещения и без учета механической вентиляции уравнение баланса массы припотолочного слоя дыма записывается в виде

M - масса слоя дыма, кг;

τ - время с момента возникновения пожара, с;

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

Если свободная граница находится ниже основания очага, будет справедливым очевидное равенство G = Ψ (где Ψ - массовая скорость газификации горючей нагрузки, кг/с). При τ = 0 уравнению баланса массы отвечает начальное условие M (0) = 0.

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

U - внутренняя энергия слоя дыма, Дж;

Q - тепловой поток, подводимый из конвективной колонки или непосредственно из очага горения, кг/с;

Q - тепловой поток, отводимый в ограждающие конструкции, Вт;

P - статическое давление газов в задымленном слое, Па;

V - объем задымленного слоя, м 3 .

Если свободная граница находится ниже основания очага, то

Q = ( Q - I )ψ,

- массовая полнота сгорания;

Q - низшая теплота сгорания ГН, Дж/кг;

I - энтальпия продуктов газификации ГН, Дж/кг.

Если же свободная граница находится выше основания очага, то

Q = C T G ,

где C и T изобарная теплоемкость и температура газов в конвективной колонке на высоте свободной границы, Дж/(кг·К) и К соответственно.

Используя соотношения термодинамики, уравнение возможно преобразовать к конечному виду

(С Р /R ) (dV / d )= Q - Q ,

где C и T - изобарная теплоемкость и приведенная газовая постоянная задымленного слоя, Дж/(кг·К). При τ = 0 этому уравнению отвечает начальное условие V (0) = 0. Как показывают численные оценки, значения С Р и R в данном уравнении допустимо принять постоянными и равными значениям этих параметров для нормальной атмосферы.

Дополнительные соотношения . Уравнения позволяют рассчитать изменение во времени массы M и объема V задымленного слоя, если определить соотношения для входящих в эти уравнения неизвестных переменных G , T , Ψ и (так как значения , Q , и C могут считаться постоянными, а величиной I можно пренебречь). Кроме того, необходимо задать соотношения для расчета основных параметров - высоты свободной границы Y и температуры слоя дыма T .

Из теории стационарной свободной конвективной струи имеем

G =Ψ + 0,21(Y - Y ) ((1 – χ ) g Q / (C T )) ,

T = ((1 – χ ) g Q / (C G )) + T ,

Для интегрирования системы уравнений пожара с заданными начальными условиями можно использовать стандартную программу (метод Рунге - Кутта) с автоматическим выбором шага интегрирования. Шаг интегрирования выбирается в соответствии с погрешностью интегрирования. Как правило, следует задавать очень невысокую погрешность.

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

2. Расчет динамики опасных факторов пожара в помещении

Александренко М.В. 1 , Акулова М.В. 2 , Ибрагимов А.М. 3

1 Студент,

Ивановский государственный политехнический университет

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПОЖАРА

Аннотация

В статье рассмотрено – виды математических моделей пожара и их область применения. Математическое моделирование позволяет спрогнозировать динамику пожара в помещениях зданий различного назначения, а следовательно позволяет вывести исследование пожарной опасности объектов на качественно новый этап развития, обеспечить переход от сравнительных методов к прогнозным, учитывающим условия эксплуатации объекта.

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

Alexandrenko M.V. 1 , Akulova M.V. 2 , Ibragimov A.M. 3

Ivanovo State Polytechnic University

MATHEMATICAL MODELLING OF THE FIRE

Abstract

The article considers types of mathematical models of the fire and their scope. Mathematical modeling allows to predict dynamics of the fire in rooms of buildings of different function and consequently allows to bring research of fire danger of objects to qualitatively new stage of development, to provide transition from comparative methods to expected, considering object service conditions.

Keywords : mathematical model, fire.

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

Моделирование пожара в помещениях основано на представлении пожара как физического явления передачи тепла и массы в соответствующих условиях его развития. Условия развития пожара характеризуются видом пожарной нагрузки и конструктивно-планировочными характеристиками здания (помещения).

По типу математического аппарата различают следующие модели: детерминированные; вероятностные; смешанные (детерминированные – вероятностные); имитационные.

Наиболее эффективным инструментом прогноза и изучения пожаров являются детерминированные математические модели.

Наряду с детерминированным моделированием следует отметить и вероятностные оценки распространения пожара на основе статистической обработки данных по реальным пожарам.

Приведём краткую характеристику каждой из моделей.

Детерминированные математические модели

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

–интегральные (модели первого поколения);

–зонные (модели второго поколения);

–полевые (CFD) (модели третьего поколения).

1.1. Интегральные математические модели

Интегральный (однозонный) метод является наиболее простым методом моделирования пожаров. Суть интегрального метода заключается в том, что состояние газовой среды оценивается через осредненные по всему объему помещения термодинамические параметры. Соответственно температура ограждающих конструкций и другие подобные параметры оцениваются как осредненные по поверхности. На основе интегрального метода были разработаны, в частности, рекомендации .

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

1.2. Зонные математические модели

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

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

1.3. Полевые математические модели

Полевые модели, обозначаемые в зарубежной литературе аббревиатурой CFD (computational fluid dynamics), являются более мощным и универсальным инструментом, чем зональные; они основываются на совершенно ином принципе. Вместо одной или нескольких больших зон в полевых моделях выделяется большое количество (обычно тысячи или десятки тысяч) маленьких контрольных объемов, никак не связанных с предполагаемой структурой потока. Для каждого из этих объемов с помощью численных методов решается система уравнений в частных производных, выражающих принципы локального сохранения массы, импульса, энергии и масс компонентов. Таким образом, динамика развития процессов определяется не априорными предположениями, а исключительно результатами расчета.

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

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

Вероятностные математические модели

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

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

Смешанные (детерминированные – вероятностные) математические модели

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

Имитационные математические модели

Имитационное моделирование представляет интерес в исследовании сложных систем при априорной неопределенности. В модели может быть задано вероятное протекание пожара, вероятные законы распределения и распространения тепловых потоков, имитируется процесс работы конструкций.

Моделирование пожара в помещении и оценка его воздействия на строительные конструкции состоит из следующих основных этапов:

Анализ конструктивно-планировочных характеристик помещения;

Определение вида, количества и размещения пожарной нагрузки;

Определение вида возможного пожара и его базовых параметров;

Выбор метода расчета и проведение расчета, оценка вероятностных характеристик пожара;

Анализ огнестойкости конструкций, определение эквивалентной продолжительности стандартного испытания.

Заключение

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

Литература

Клуб студентов «Технарь». Конспекты по математическим моделям [Электронный курс] URL: http://www.c-stud.ru (дата обращения 10.03.2015)
Расчет необходимого времени эвакуации людей из помещений при пожаре: Рекомендации. – М.: ВНИИПО МВД СССР, 1989. – 22 с.
Методические рекомендации «Применение полевого метода математического моделирования пожара в помещениях.
ГОСТ 12.1.004-91* Пожарная безопасность. Общие требования.
СНиП 21-01-97* Пожарная безопасность зданий и сооружений.

References

Club of students “Technician”. Abstracts on mathematical models of URL: http://www.c-stud.ru (date of the address 10.03.2015)
Calculation of necessary time of evacuation of people from rooms at the fire: Recommendations. – M.: VNIIPO MVD USSR, 1989. – 22 s.
Methodical recommendations “Application of a field method of mathematical modeling of the fire in rooms.
GOST 12.1.004-91 * Fire safety. General requirements.
SNiP 21-01-97 * Fire safety of buildings and constructions.

ЛЕКЦИЯ

по дисциплине "Прогнозирование опасных факторов пожара"

Тема №2. «Основные понятия и уравнения интегральной математической модели пожара в помещении»

План лекции:

1.2 Среднеобъемная плотность газовой среды

Лекция 2.

Цели лекции:

Учебные

В результате прослушивания материала слушатели должны знать:

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

Уметь: прогнозировать обстановку на пожаре.

Развивающие:

выделять самое главное
самостоятельность и гибкости мышления
развитие познавательного мышления

Литература

Д.М. Рожков Прогнозирование опасных факторов пожара в помещении. Иркутск 2007. С.89
Ю.А.Кошмаров, М.П. Башкирцев Термодинамика и теплопередача в пожарном деле. ВИПТШ МВД СССР, М., 1987 г.
Ю.А.Кошмаров Прогнозирование опасных факторов пожара в помещении. Москва 2000. С.118
Ю.А.Кошмаров, В.В. Рубцов, Процессы нарастания опасных факторов пожара в производственных помещениях и расчет критической продолжительности пожара. МИПБ МВД России, М., 1999 г.

Лекция 1. Основные понятия математической модели пожара в помещении

1.1 Допущения интегрального метода термодинамического анализа пожара

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

1) С позиций термодинамики газовая среда, заполняющая помещение с проемами (окна, двери и т.п.), как объект исследования есть открытая термодинамическая система (рис. 1.1).

2)Ограждающие конструкции (пол, потолок, стены) и наружный воздух (атмосфера) являются внешней средой по отношению к этой термодинамической системе. Граница между термодинамической системой и внешней средой (контрольная поверхность) показана условно на рис. 1.1 пунктирной линией. Эта система взаимодействует с внешней средой путем тепло- и массообмена. В процессе развития пожара через одни проемы выталкиваются из помещения нагретые газы, а через другие поступает холодный воздух.

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

4) Термодинамическая система в свою очередь совершает работу, выталкивая нагретые газы во внешнюю атмосферу. Эта термодинамическая система взаимодействует также с ограждающими конструкциями путем теплообмена. Кроме того, в эту систему с поверхности горящего материала (т.е. из пламенной зоны) поступает вещество в виде газообразных продуктов горения.

Рис. 1.1. Схема пожара в помещении:

Контрольная поверхность;1 - ограждения; 2 - проемы (окна, двери); 3 горящий материал; G г - расход уходящих газов; G в - расход поступающего холодного воздуха; ψ- скорость выгорания материала

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

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

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

(2.1)

где р - локальное давление, Н·м -2 ; ρ - локальная плотность, кг·м -3 ; R - газовая постоянная, Дж·кг -1 К -1 ; Т - локальная температура, К.

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

7) Особенностью рассматриваемой термодинамической системы (т.е. газовой среды в помещении) является то, что ее объем (т.е. пространственная конфигурация) в процессе развития пожара практически не изменяется. В связи с этим в интегральном методе описания состояния термодинамической системы, коей является газовая среда в помещении, используются "интегральные" параметры состояния термодинамической системы среднеобъемные параметры - среднеобъемную плотность газовой среды и среднеобъемную (удельную) внутреннюю энергию.

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

1.2 Среднеобъемная плотность газовой среды в помещении представляет собой отношение массы газа, заполняющего помещение, к объему помещения, т.е.

(2.2)

где М - масса газа, заполняющего помещение, кг; V - свободный объем помещения, м 3 . Нижний индекс т , используемый здесь и далее, представляет собой первую букву в немецком слове mittel (средний). Следует отметить, что

(2.3)

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

(2.4)

Газовая среда в помещении представляет собой смесь кислорода, азота и продуктов горения. В процессе развития пожара количественное соотношение между компонентами смеси изменяется. В интегральном методе описания процесса изменения массы i -го компонента смеси в течение времени используется параметр, называемый среднеобъемной парциальной плотностью i -го компонента смеси.

1.3 Среднеобъемная парциальная плотность i -го компонента представляет собой отношение массы i -го компонента смеси (например О 2 ), содержащейся в объеме помещения, к объему помещения, т.е.

(1.5)

(2.5)

где М, - масса i -го компонента, находящегося в помещении, кг. Отметим, что с формальной точки зрения среднеобъемная парциальная плотность i -го компонента есть результат осреднения по объему помещения всех значений локальной парциальной плотности этого компонента, т.е.

(2.6)

где ρ i , - локальное значение парциальной плотности i -го компонента, кг·м -3 .

1.4 Среднеобъемная (удельная) внутренняя энергия представляет собой отношение внутренней тепловой энергии всего газа, заполняющего помещение, к объему помещения, т.е.

(2.7)

где и - внутренняя энергия всей газовой среды, заполняющей помещение, Дж. С формальных позиций среднеобъемная внутренняя энергия газовой среды есть результат осреднения по объему всех значений локальной удельной (объемной) внутренней энергии, т.е.

(2.8)

где U V - локальное значение удельной (объемной) внутренней энергии, Дж·м -3 . Локальные значения удельной объемной внутренней энергии и удельной массовой внутренней энергии связаны между собой простым соотношением, которое имеет следующий вид:

(2.9)

где и - локальное значение удельной массовой внутренней энергии газа, Дж·кг. Отметим здесь, что между локальным значением удельной массовой внутренней энергии и локальной температурой идеального газа существует простая взаимосвязь, а именно

(2.10)

где c v - изохорная теплоемкость газа, Дж·кг·К.

В интегральном методе описания процесса изменения состояния термодинамической системы (т.е. газовой среды в помещении) вместо среднеобъемной внутренней энергии используется параметр состояния, называемый среднеобъемным давлением. Эти два параметра в формальном отношении являются взаимозаменяемыми. Покажем это. Формулу (2.8) можно преобразовать с помощью выражений (2.9) и (2.10)

(2.11)

Если теперь воспользоваться уравнением Клапейрона (2.1), то формулу (2.11) можно преобразовать и получить следующее выражение:

(2 . 12)

где p - локальное давление, Н·м -2 ;

к = C p / C V - отношение изобарной и изо хорной теплоемкостей идеального газа (показатель адиабаты). С достаточной для практики точностью можно считать, что показатель адиабаты во всех точках внутри помещения есть одна и та же постоянная величина. С учетом этого замечания формулу (2.12) можно преобразовать:

(2.13)

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

(2.14)

где р т - среднеобъемное давление, Н·м -2

Сравнивая выражения (2.13) и (2.14), получим следующее соотношение между среднеобъемной внутренней энергией и среднеобъемным давлением:

(2.15)

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

Степень нагретости газовой среды характеризуется в среднем отношением внутренней энергии этой среды к ее массе. Отношение этих физических величин можно представить с помощью формул (2.2), (2.7) и (2.15) в следующем виде:

(2.16)

Если правую и левую части равенства (2.16) поделить на изохорную теплоемкость, то получится следующее выражение:

(2.17)

Комплекс в левой части выражения (2.17) имеет размерность "Кельвин". Этот комплекс представляет собой параметр состояния рассматриваемой термодинамической системы, который называется среднемассовой температурой газовой среды, т.е.

(2.18)

С помощью выражения (2.18) можно преобразовать формулу (2. ] 7) и в результате получить следующее уравнение:

(2.19)

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

1.5 Дым и его влияние на термодинамические параметры среды

Газовая среда, заполняющая помещение при пожаре, содержит в себе мельчайшие твердые частицы. Следует отметить, что доля тепловой энергии, приходящейся на эти частицы, пренебрежимо мала по сравнению с внутренней энергией газовой среды, находящейся в помещении. Не существенным является также вклад этих частиц в суммарную массу среды, заполняющей помещение при пожаре. Поэтому можно не учитывать присутствие этих частиц при вычислениях таких параметров состояния среды, как среднеобъемная плотность, среднеобъемное давление и среднемассовая температура. Однако присутствие этих частиц сильно изменяет оптические свойства среды в помещении. В результате рассеяния энергии световых волн из-за многократного диффузного отражения от этих мельчайших частиц (их диаметр приблизительно равен 0,2-4 мкм) ухудшается видимость. Оптические свойства среды, находящейся в помещении, характеризуются среднеобъемной оптической плотностью дыма.

Среднеобъемная плотность (концентрация) дыма представляет собой отношение оптического количества дыма, находящегося в помещении, к объему помещения, т.е.

(2.20)

где S - оптическое количество дыма, Нп·м 2 ; µ m - среднеобъемная оптическая плотность дыма, Нп·м -1 . Здесь сокращением «Нп» обозначено слово "Непер". Оптическое количество дыма в помещении есть произведение средней концентрации твердых частиц на объем помещения и эффективное сечение экстинкции, т.е.

S = NVx , (2.21)

где N - средняя концентрация частиц, т.е. число частиц, приходящееся на единицу объема, м -3 ; χ - эффективное сечение экстинкции, м 2 . Чем выше оптическая плотность (концентрация) дыма, тем хуже видимость в помещении. Оптическая плотность дыма и дальность видимости связаны между собой следующим приближенным соотношением:

(2.22)

где l вид - дальность видимости, м.

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

Дымообразующая способность горючего материала есть оптическое количество дыма, образующегося при сгорании единицы массы горючего материала, т.е.

D = J χ , (2.23)

где D - дымообразующая способность ГМ, Нп·м 2 ·кг -1 ; J - число частиц, образующихся при сгорании единицы массы горючего материала, кг -1 ; χ - эффективное сечение экстинкции частиц, м 2 .

Лекция 2. Дифференциальные уравнения пожара

Уравнения пожара описывают в самом общем виде изменение среднеобъемных параметров состояния газовой среды в помещении в течение времени (в процессе развития пожара). Эти уравнения были сформулированы в 1976г. проф. Ю.А. Кошмаровым (статья "Развитие пожара в помещении" в научном сборнике ВНИИПО МВД СССР "Горение и проблемы тушения пожаров". М.: ВНИИПО МВД СССР, 1977).

Уравнения пожара являются обыкновенными дифференциальными уравнениями. Они вытекают, как и большинство уравнений математической физики, из фундаментальных законов природы - первого закона термодинамики для открытой термодинамической системы и закона сохранения массы. Подробный вывод этих уравнений приведен в учебнике Ю.А. Кошмарова и М.П. Башкирцева "Термодинамика и теплопередача в пожарном деле" (М., ВИПТШ МВД СССР, 1987). Ограничимся здесь кратким изложением рассуждений, используемых при выводе уравнений пожара.

Первое уравнение - уравнение материального баланса пожара в помещении - вытекает из закона сохранения массы. Применительно к газовой среде, заполняющей помещение, этот закон можно сформулировать так: изменение массы газовой среды в помещении за единицу времени равно алгебраической сумме потоков массы через границы рассматриваемой термодинамической системы. Под границей системы здесь подразумевается воображаемая контрольная поверхность, ограничивающая пространство, внутри которого заключена рассматриваемая газовая среда. На рис. 1.1 эта поверхность условно показана пунктирной линией. Часть этой поверхности совпадает с поверхностью ограждений (стены, пол, потолок). Там, где находятся проемы, эта поверхность является воображаемой. Объем пространства, заключенный внутри этой поверхности, называется свободным объемом помещения и обозначается буквой V . Введем следующие обозначения:

а) G B - расход поступающего воздуха из окружающей атмосферы в помещение, который имеет место в рассматриваемый момент времени процесса развития пожара, кг∙с -1 ;

б) G Г - расход газов, покидающих помещение через проемы в рассматриваемый момент времени, кг∙с -1 ;

в) ψ - скорость выгорания (скорость газификации) горючего материала в рассматриваемый момент времени, кг∙с -1 ;

г) ρ m V - масса газовой среды, заполняющей помещение в рассматриваемый момент времени, кг.

За малый промежуток времени, равный dx , будет иметь место малое изменение массы газовой среды. В то же время можно считать, что значения G Г , G B и ψ в течение этого малого промежутка времени остаются практически неизменными. С учетом вышесказанного уравнение материального баланса для газовой среды в помещении записывается следующим образом:

(2.24)

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

Уравнение (2.24) называется уравнением материального баланса пожара.

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

(2.25)

Уравнение баланса токсичного продукта горения:

(2.26)

Уравнение баланса оптического количества дыма:

(2.27)

В этих уравнениях использованы следующие обозначения: ρ 1 , - среднеобъемная парциальная плотность кислорода, кг · м -3 ; ρ 2 - среднеобъемная парциальная плотность токсичного продукта горения, кг · м -3 ; μ м - объемная оптическая концентрация дыма, Нп · м -1 .

В правой части уравнения (2.25) - уравнения баланса массы кислорода - использованы, кроме ранее указанных, следующие обозначения: х 1в - массовая доля кислорода в поступающем воздухе; средняя массовая доля кислорода в помещении; L 1 - стехиометрический коэффициент для кислорода (количество кислорода, необходимое для сгорания единицы массы горючего материала), кг∙кг -1 ; η - коэффициент полноты сгорания; n 3 , - коэффициент, учитывающий отличие концентрации кислорода в уходящих газах от среднеобъемной концентрации кислорода.

В правой части уравнения (2.26) - уравнения баланса токсичного продукта горения - использованы, кроме ранее указанных, следующие обозначения: L 2 - стехиометрический коэффициент для продукта горения (количество продукта горения, образующегося при сгорании единицы массы горючего материала), кг∙кг -1 ; средняя массовая доля токсичного газа в помещении; п 2 - коэффициент, учитывающий отличие концентрации токсичного газа в уходящих газах от среднеобъемной концентрации этого газа.

В правой части уравнения (1.36) - уравнения баланса оптического количества дыма - использованы, кроме ранее указанных, следующие обозначения: n 3 - коэффициент, учитывающий отличие оптической концентрации дыма в уходящих газах от среднеобъемного значения оптической концентрации дыма; F w - площадь поверхности ограждений (потолка, пола, стен), м 2 ; к с - коэффициент седиментации частиц дыма на поверхностях ограждающих конструкций, Нп · с -1 . Коэффициент седиментации по физическому смыслу есть скорость осаждения частиц дыма.

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

С учетом всего сказанного получается следующее уравнение энергии пожара:

(2.28)

Левая часть этого уравнения есть скорость изменения внутренней тепловой энергии газовой среды в помещении за единицу времени в рассматриваемый малый промежуток времени dτ , т.е.

(2.29)

В правой части уравнения (2.28) первый член представляет собой количество тепла, поступающего за единицу времени в газовую среду в результате горения (скорость тепловыделения). Второй член есть поток энергии в помещение, поступающий вместе с продуктами газификации (пиролиз, испарение) горючего материала. Здесь величина i r - энтальпия этих продуктов. Третий член представляет собой сумму внутренней тепловой энергии поступающего за единицу времени воздуха и работы проталкивания, которую совершает внешняя атмосфера. Четвертый член есть сумма внутренней тепловой энергии, которую уносят за единицу времени уходящие газы, и работы выталкивания, которую совершает рассматриваемая термодинамическая система. Пятый член представляет собой тепловой поток, поглощаемый ограничивающими конструкциями и излучаемый через проемы.

Представленные выше пять дифференциальных уравнений содержат шесть неизвестных функций p m (τ), p m (τ), Т m (τ), р 1 (τ), р 2 (τ) и  m (τ) . Эту систему уравнений дополняет алгебраическое уравнение - усредненное уравнение состояния (2.19).

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

(2.30)

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

(2.31)

где G o в - массовый расход подачи огнетушащего вещества, кг∙с -1 . Соответствующим образом изменяются в этом случае и остальные дифференциальные уравнения пожара.

Как уже говорилось, в уравнениях пожара искомыми (неизвестными) функциями являются среднеобъемные параметры газовой среды, а независимой переменной является время. Кроме этих переменных величин, уравнения содержат целый ряд других физических величин, которые можно разделить на две группы. К первой группе относятся величины, заданные условиями однозначности, которые представляют собой сведения о размерах помещения (объем V и поверхность ограждений F w ) и свойствах горючего материала (теплота сгорания Q р н , стехиометрические коэффициенты L 1 , L 2 , дымообразующая способность D , энтальпия продуктов горения i n . Ко второй группе относятся те величины, которые зависят, помимо всего прочего, от параметров состояния среды в помещении. К этим величинам относятся массовые расходы поступающего через проемы воздуха G B и уходящих через проемы газов G Г , тепловой поток, поглощаемый ограждающими конструкциями и излучаемый через проемы Q w , коэффициент полноты сгорания η , скорость тепловыделения ηQ p н ψ . Для вычисления значений физических величин, относящихся ко второй группе, необходимо располагать дополнительными уравнениями.

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

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

В интегральной математической модели мы оперируем с интегральными характеристиками термодинамической системы. Этот подход не требует каких-либо допущений и оговорок о том, как распределены локальные значения термодинамических параметров состояния по объему помещения. Здесь не уместны оговорки такого, например, типа: "предположим, что температурное поле является однородным", или часто используемое выражение о "размазанности" того или иного параметра состояния газовой среды.

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

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