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

Интегральная модель пожара
Зонная модель пожара

Общие сведения о расчете пожаров. Опасные факторы пожара.

Расчет пожара (прогнозирование опасных факторов) необходим для оценки своевременности эвакуации и разработке мероприятий по ее совершенствованию, при создании и совершенствовании систем сигнализации, оповещения и тушения пожаров, при разработке планов пожаротушения (планирования боевых действий пожарных подразделений при пожаре), для оценки фактических пределов огнестойкости, проведении пожарно-технических экспертиз и других целей.
В развитии пожара в помещении обычно выделяют три стадии:
- начальная стадия - от возникновения локального неконтролируемого очага горения до полного охвата помещения пламенем; при этом средняя температура среды в помещении имеет не высокие значения, но внутри и вокруг зоны горения температура такова, что скорость тепловыделения выше скорости отвода тепла из зоны горения, что обуславливает само ускорение процесса горения;
- стадия полного развития пожара - горят все горючие вещества и материалы, находящиеся в помещении; интенсивность тепловыделения от горящих объектов достигает максимума, что приводит и к быстрому нарастанию температуры среды помещения до максимальных значений;
- стадия затухания пожара - интенсивность процесса горения в помещении снижается из-за расходования находящейся в нём массы горючих материалов или воздействия средств тушения пожара.
Однако в любом случае, как показывает уравнение «стандартного пожара», температура в очаге пожара через 1,125 мин достигает значения 365оС. Поэтому очевидно, что возможное время эвакуации людей из помещений не может превосходить продолжительности начальной стадии пожара.
В начальной стадии развития пожара опасными для человека факторами являются: пламя, высокая температура, интенсивность теплового излучения, токсичные продукты горения, дым, снижение содержания кислорода в воздухе, поскольку при достижении определённых уровней они поражают его организм, особенно при синергическом воздействии.
Исследованиями отечественных и зарубежных учёных установлено, что максимальная температура, кратковременно переносимая человеком в сухой атмосфере, составляет 149 0С, во влажной атмосфере вторую степень ожога вызывало воздействие температуры 55 0С в течение 20с и 70 0С при воздействии в течение 1с; а плотность лучистых тепловых потоков 3500 вт/м2 вызывает практически мгновенно ожоги дыхательных путей и открытых участков кожи; концентрации токсичных веществ в воздухе приводят к летальному исходу: окиси углерода (СО) в 1,0% за 2-3 мин, двуокиси углерода (СО2) в 5% за 5 мин., цианистого водорода (HCN) в 0,005% практически мгновенно; при концентрации хлористого водорода (HCL) 0,01- 0,015% останавливается дыхание; при снижении концентрации кислорода в воздухе с 23% до 16% ухудшаются двигательные функции организма, и мускульная координация нарушается до такой степени, что самостоятельное движение людей становится невозможным, а снижение концентрации кислорода до 9% приводит к смерти через 5 минут.
Совместное действие некоторых факторов усиливает их воздействие на организм человека (синергический эффект). Так токсичность окиси углерода увеличивается при наличии дыма, влажности среды, снижении концентрации кислорода и повышении температуры. Синергетический эффект обнаруживается и при совместном действии двуокиси азота и понижении концентрации кислорода при повышенной температуре, а также при совместном воздействии цианистого водорода и окиси углерода.
Особое воздействие на людей оказывает дым. Дым представляет собой смесь несгоревших частиц углерода с размерами частиц от 0,05 до 5,0 мкм. На этих частицах конденсируются токсичные газы. Поэтому воздействие дыма на человека также имеет, по-видимому, синергический эффект.
В действительности при пожаре выделяется значительно больше токсинов, воздействие которых достаточно хорошо изучено (табл. 1,2). Максимально допустимый уровень опасных (основных) факторов пожара, воздействие которого не приносит вреда человеку (табл.3), нормирован. Вырываясь из помещения, опасные факторы пожара, прежде всего дым, стремительно распространяются по коммуникационным путям здания.

Источники. 1-4, 6 - ГОСТ 12.1.004-91; 5 - ГОСТ 12.3.047-98; 7 - Кошмаров Ю. А. Прогнозирование опасных факторов пожара в помещении: Учеб. пособие. - М.: Академия ГПС МВД РФ, 2000.

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

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

Интегральная модель пожара

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

Зонная модель пожара

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

Рисунок 1.

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

Полевой (дифференциальный) метод расчета

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

Рис. 2. Расчеты с помощью полевой модели.

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

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

Критерии выбора моделей пожара для расчетов

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

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

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

Характеристика типовой пожарной нагрузки (примеры)

Здания I-II ст. огнест.; мебель+бытовые изделия
Низшая теплота сгорания, кДж/кг 13800,0
Линейная скорость пламени, м/с / Плотность ГЖ,кг/м3 0,0108
Удельная скорость выгорания, кг/м2-с 0,01450
Дымообразующая способность, Нпм2/кг 270,00
Потребление кислорода (О2), кг/кг -1,0300
Выделение газа:
углекислого (СОг), кг/кг 0,20300
угарного (СО), кг/кг 0,00220
хлористого водорода (НС1), кг/кг 0,01400

Здание I-II ст. огнест.; мебель+ткани
Низшая теплота сгорания, кДж/кг 14700,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/м3. 0,0108
Удельная скорость выгорания, кг/м2с 0,01450
Дымообразуюшая способность, Нпм2/кг. ...82,00
Потребление кислорода (O2), кг/кг -1,4370
Выделение газа:
углекислого (СО2). кг/кг...... 1,28500
угарного (СО), кг/кг 0,00220
хлористого водорода (НС1), кг/кг. 0,00600

Обществ.здания; мебель+линолеум ПВХ (0,9+0,1)
Низшая теплота сгорания, кДж/кг 14000,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/м3 0,015
Удельная скорость выгорания, кг/м2с.-. 0,01370
Дымообразуюшая способность, Нпм2/кг 47,70
Потребление кислорода (Ог), кг/кг -1,3690
Выделение газа:
углекислого (СО2), кг/кг 1,47800
угарного (СО), кг/кг 0,03000
хлористого водорода (НС1), кг/кг.. 0,00580

Библиотеки, архивы; книги, журналы на стеллажах
Низшая теплота сгорания, кДж/кг 14500,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/м3 0,0103
Удельная скорость выгорания, кг/м2с 0,01100
Дымообразуюшая способность, Нпм2/кг 49,50
Потребление кислорода (О2), кг/кг -1,1540
Выделение газа:
углекислого (СО2), кг/кг 1,10870
угарного (СО), кг/кг 0,09740
хлористого водорода (НС1), кг/кг. .0,00000

Верхняя одежда; ворс, ткани (шерсть+нейлон)
Низшая теплота сгорания, кДж/кг 23300,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/м3 0,0835
Удельная скорость выгорания, кг/м2-с 0,01300
Дьшообразуюшая способность, Нпм2/кг 129,00
Потребление кислорода (О2), кг/кг -3,6980
Выделение газа:
углекислого (СО2), кг/кг 0,46700
угарного (СО), кг/кг 0,01450
хлористого водорода (HС1), кг/кг 0,00000

Резинотехн. изделия; резина, изделия из нее
Низшая теплота сгорания, кДж/кг 36000,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/м3.... 0,0184
Удельная скорость выгорания, кг/м2-с 0,01120
Дымообразуюшая способность, Нп м2/кг 850,00
Потребление кислорода (О2), кг/кг -2,9900
Выделение газа:
углекислого (СО2), кг/кг 0,41600
угарного (СО), кг/кг.. 0,01500
хлористого водорода (НС1), кг/кг 0,00000

Автомобиль; 0,3*(резина, бензин)+0,15*(ППУ, искожа ПВХ)+0,1* эмаль
Низшая теплота сгорания, кДж/кг 31700,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/м3 0,0068
Удельная скорость выгорания, кг/м2 с 0,02330
Дымообразуюшая способность, Нп м2/кг 487,00
Потребление кислорода (О2), кг/кг. -2,6400
Выделение газа:
углекислого (СО2), кг/кг 1,29500
угарного (СО), кг/кг 0,09700

Кабинет; мебель+бумага (0,75+0,25)
Низшая теплота сгорания, кДж/кг.14002,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/м3 0,042
Удельная скорость выгорания, кг/м2с.0,01290
Дымообразуюшая способность, Нпм2/кг.. 53,00
Потребление кислорода (О2), кг/кг. .-1,1610
Выделение газа:
углекислого (СО2), кг/кг...0,64200
угарного (СО), кг/кг....... 0,03170
хлористого водорода (НС1), кг/кг. , 0,00000

Помещение, облицованное панелями; панели ДВП
Низшая теплота сгорания, кДж/кг 18100,0
Линейная скорость пламени, м/с / Плотность ГЖ, кг/мЗ 0,0405
Удельная скорость выгорания, кг/м2с 0,01430
Дымообразуюшая способность, Нпм2/кг 130,00
Потребление кислорода (О2), кг/кг -1,1500
Выделение газа:
углекислого (СО2), кг/кг 0,68600
угарного (СО), кг/кг 0,02150
хлористого водорода (НС1), кг/кг.... г.. 0,00000

Литература

Федеральный закон РФ от 22 июля 2008 г. № 123-ФЗ «Технический регламент о требованиях пожарной безопасности».
ГОСТ 12.1.004-91* Пожарная безопасность. Общие требования.
ГОСТ 12.1.033-81* Пожарная безопасность. Термины и определения.
СП 118.13330.2012 Общественные здания и сооружения.
СНиП 21-01-97* Пожарная безопасность зданий и сооружений.
Холщевников В.В., Самошин Д.А. Парфененко А.П., Кудрин И.С., Истратов Р.Н., Белосхов И.Р.Эвакуация и поведение людей при пожарах: Учеб. пособие. - М.: Академия ГПС МЧС России, 2015. - 262 с.

МИНИСТЕРСТВО РОССИЙСКОЙ ФЕДЕРАЦИИ ПО ДЕЛАМ ГРАЖДАНСКОЙ ОБОРОНЫ, ЧРЕЗВЫЧАЙНЫМ СИТУАЦИЯМ И ЛИКВИДАЦИИ ПОСЛЕДСТВИЙ СТИХИЙНЫХ БЕДСТВИЙ

Федеральное государственное учреждение "Всероссийский ордена "Знак Почета" научно-исследовательский институт противопожарной обороны"

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

Представлено описание основных уравнений полевого метода моделирования пожаров, известного в зарубежной литературе под наименованием CFD (computational fluid dynamics). Указана рекомендуемая область применения метода. Изложен порядок проведения расчетной оценки пожарной опасности конкретных объектов.

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

Рекомендации разработаны сотрудниками ФГУ ВНИИПО МЧС России канд. техн. наук A.M. Рыжовым, д-ром техн. наук И.Р. Хасановым, канд. техн. наук А.В. Карповым, А.В. Волковым, В.В. Лицкевичем, канд. техн. наук А.А. Дектеревым.

СПИСОК ОБОЗНАЧЕНИЙ

С m , С 1 , С 2 - константы в модели турбулентности;

с Р - удельная массовая изобарная теплоемкость, Дж/(кг×К);

f - функция смешения;

G k - генерация турбулентности за счет вынужденной конвекции, Па/с;

G B - генерация турбулентности за счет естественной конвекции, Па/с;

g - ускорение свободного падения, м/с 2 ;

H k - теплота образования k -гo компонента смеси, Дж/кг;

Удельная массовая энтальпия смеси, Дж/кг;

k - кинетическая энергия турбулентных пульсаций, м 2 /с 2 ;

m - масса, кг;

р - динамическое давление, Па;

R - приведенная газовая постоянная, Дж/(кг×К);

s - стехиометрическое отношение;

S Ф - источниковый член;

t - время, с;

Т - термодинамическая (абсолютная) температура, К;

u , v , w - проекции вектора скорости соответственно на оси х , у , z в декартовых и х , r , j в цилиндрических координатах, м/с;

Y k - массовая концентрация k -го компонента смеси, кг/кг;

b - коэффициент объемного расширения, 1/К;

Г Ф - коэффициент переноса;

e - скорость диссипации кинетической энергии турбулентности, м 2 /с 3 ;

F - обобщенная переменная;

l - коэффициент теплопроводности, Вт/(м×К);

m - ламинарная динамическая вязкость, Па×с;

m t - турбулентная динамическая вязкость, Па×с;

m эфф - эффективная динамическая вязкость, Па×с;

v - кинематическая вязкость, м 2 /с;

r - плотность, кг/м 3 ;

s k , s e - аналоги критерия Прандтля для уравнений кинетической энергии турбулентных пульсаций и скорости ее диссипации;

c R - доля тепла, теряемая за счет излучения.

ВВЕДЕНИЕ

В последние годы во многих странах мира (Англия, США, Япония, Австралия и др.) наметился переход к гибкому (объектно-ориентированному) нормированию, которое позволяет наиболее оптимальным образом обеспечить пожарную безопасность объекта с учетом его индивидуальных особенностей, в отличие от "жесткого" нормирования, предписывающего соблюдение определенных положений для любого объекта, относящегося к данному классу.

В ряде отечественных норм также реализуются элементы гибкого нормирования, например в ГОСТ 12.1.004-91* и СНиП 21-01-97* .

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

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

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

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

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

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

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

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

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

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

1. ОБЩИЕ ПОЛОЖЕНИЯ

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

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

2. ОБЛАСТЬ ПРИМЕНЕНИЯ

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

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

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

выбора рационального варианта противопожарной защиты конкретных объектов.

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

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

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

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

помещений, в которых один из геометрических размеров гораздо больше остальных;

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

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

3. ОСНОВЫ ПОЛЕВОГО МЕТОДА МОДЕЛИРОВАНИЯ ПОЖАРОВ

3.1. Основные уравнения

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

Уравнение сохранения массы:

Уравнение сохранения импульса:

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

где - статическая энтальпия смеси;

Н k - теплота образования k -го компонента; - теплоемкость смеси при постоянном давлении; - радиационный поток энергии в направлении x j .

Уравнение сохранения химического компонента k :

Для замыкания системы уравнений (3.1)-(3.5) используется уравнение состояния идеального газа. Для смеси газов оно имеет следующий вид:

где R о - универсальная газовая постоянная; M k - молярная масса k -гo компонента.

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

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

Это уравнение очень похоже на исходное уравнение (3.1). Отличие состоит в появившемся дополнительном члене , который представляет собой турбулентный перенос массы из-за флуктуации плотности и скорости.

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

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

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

При этом уравнения сохранения принимают вид, аналогичный системе (3.9)-(3.12), однако они учитывают флуктуации плотности, что существенно при рассмотрении областей, где происходит горение.

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

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

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

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

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

где k 1/2 - характеристическая скорость; k = /2 - турбулентная кинетическая энергия; l - характеристическая длина смешения; - константа.

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

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

Величина v t определяется с помощью модели турбулентности. Наибольшее распространение при моделировании пожаров получила k -e модель. В ней решаются два уравнения переноса, аналогичные уравнениям (3.9)-(3.12): одно для турбулентной кинетической энергии k и второе для вязкостной диссипации этой энергии e во внутреннюю энергию жидкости. Уравнение переноса для k можно вывести из осредненных по времени уравнений сохранения импульса:

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

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

Используя концепцию вихревой вязкости, уравнение (3.18) можно записать в виде

где С 1 , С 2 , С 3 и s e - эмпирические константы. Источниковые члены, обусловленные вязкостными напряжениями и плавучестью, определяются выражениями:

Систему уравнений (3.9)-(3.12), (3.18), (3.23) часто записывают в форме обобщенного уравнения переноса:

где Ф - консервативная величина (скаляр), Г Ф - соответствующий ей коэффициент переноса; S Ф - источникоый член.

Уравнение (3.26) описывает сохранение импульса при Ф = h , сохранение энергии при Ф = u i , сохранение массы при Ф = 1, сохранение массы компонентов при Ф = Y k , перенос кинетической энергии турбулентности при Ф = k и скорости ее диссипации при Ф = e.

3.3. Модели горения

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

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

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

При моделировании пожаров часто бывает достаточно представить процесс горения в виде одной одноступенчатой реакции:

F + sO ®(1 + s )P , (3.27)

где F , О и Р обозначают массы горючего, окислителя и продукта соответственно.

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

где b = Y f - (Y 0 /s ) - консервативная переменная Шваба-Зельдовича, а индексы f и 0 относятся к горючему и окислителю соответственно. Если предположить, что коэффициенты диффузии компонентов равны, становится возможным избавиться от источникового члена при определении степени смешения топлива и окислителя. Если реакция необратима и можно предположить, что она протекает бесконечно быстро, то локальные массовые доли можно определить непосредственно через среднее по времени значение функции смешения f :

где Y ox ,0 - массовая доля кислорода в потоке окислителя, a Y f , f - массовая доля топлива в потоке газообразных продуктов пиролиза.

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

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

где С = 4, а В полагают равным 2.

Предположение для замыкания источникового члена (формула (3.31)) позволяет, помимо уравнения переноса для f , решать уравнение для массовой доли топлива и рассчитывать массовую долю каждого компонента упрощенной химической реакции. Модели этого типа успешно использовались при решении различных задач пожарной безопасности и оптимизации процесса горения в промышленных установках. Достоинством модели является ее простота. Она позволяет рассчитывать распределенное по объему выделение энергии, определяемое геометрией помещения и доступом воздуха. Можно определить концентрации CO 2 и Н 2 O, если предположить, что они являются единственными продуктами горения.

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

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

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

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

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

Основное уравнение радиационного переноса можно записать в виде

где I - интенсивность радиационного излучения в направлении W; s - расстояние в направлении W; E g = s - энергия, излучаемая абсолютно черным газом при температуре газа T g ; k a и k s - коэффициенты поглощения и рассеяния; Р (W, W") - вероятность того, что излучение в направлении W" после рассеяния попадет в телесный угол d W в окрестности направления W. Это уравнение необходимо интегрировать по всем направлениям и длинам волн. Для большинства практических задач точное решение невозможно, вместо него разработано несколько приближенных методов, которые и используются для моделирования динамики пожаров в помещениях.

3.4.1. Потоковые методы

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

Если телесные углы совпадают с поверхностями контрольного объема в декартовом пространстве и если предположить, что поток излучения через каждую поверхность однороден, то, обозначив через F i + тепловой поток, проходящий через контрольный объем в положительном направлении i , и через F i - - поток в отрицательном направлении i , имеем:

где k a и k s - локальные коэффициенты поглощения и рассеяния, а Е b - количество тепла, излучаемого контрольным объемом, если он является абсолютно черным.

Объединяя эти уравнения и дифференцируя их по x i получаем:

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

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

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

Эта модель, разработанная Локвудом и Шахом , преодолевает основной недостаток потоковых методов. Для нее характерны некоторые черты методов Монте-Карло, а именно прохождение "лучей" электромагнитного излучения через вычислительную область между границами. Однако в отличие от методов Монте-Карло, где направления лучей генерируются случайным образом, в этой модели они выбираются предварительно, таким же образом, как выбирается расположение гидродинамической сетки. Метод включает в себя решение уравнения радиационного переноса вдоль путей этих лучей, выбираемых обычно таким образом, чтобы они приходили в центры граничных поверхностей гидродинамических контрольных объемов.

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

Для каждого луча при его прохождении от одной границы до другой решается уравнение радиационного переноса (3.32). Если для краткости ввести: коэффициент ослабления k e = k a + k s , оптическую глубину элемента ds * = k e ds и модифицированную энергию излучения

то уравнение переноса можно переписать в виде

Для элементарного контрольного объема, в котором температуру можно считать постоянной, уравнение можно проинтегрировать и привести к виду

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

где I n и I n +1 - соответственно значения интенсивности излучения, входящего и выходящего из n -го контрольного объема;

ds* - оптическая длина контрольного объема.

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

где N - общее количество лучей, dА - площадь поверхности ячейки.

4. ЗАМЫКАНИЕ ОСНОВНОЙ СИСТЕМЫ УРАВНЕНИЙ.

УСЛОВИЯ ОДНОЗНАЧНОСТИ

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

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

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

условия на твердых негорючих поверхностях;

условия на плоскости (оси) симметрии;

условия, характеризующие работу приточно-вытяжной вентиляции;

условия на свободной границе;

условия на поверхности горючего.

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

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

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

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

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

Рассчитывается безразмерное расстояние у + до ближайшего пристеночного узла:

где k p - значение кинетической энергии турбулентности, рассчитанное при решении соответствующего уравнения переноса с использованием граничного условия на стенке k = 0; у р - размерное расстояние от ближайшего пристеночного узла до стенки, м.

Рассчитывается значение безразмерной скорости и + :

Определяется значение безразмерной энтальпии h + :

h + = Pr t (u + +П),

где Pr t - турбулентное число Прандтля; П - сопротивление ламинарного подслоя переносу энергии:

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

где h w - энтальпия ближайшего узла внутри стенки; h p - энтальпия ближайшего пристеночного узла.

Значение скорости диссипации турбулентной кинетической энергии определяется из соотношения

На плоскости (оси) симметрии традиционно используется условие v n = 0 для нормальной компоненты скорости и условие d Ф/dn = 0 - для остальных переменных.

Для описания вентиляционного потока, подаваемого (удаляемого) через границу расчетной области, как правило, задается значение скорости потока. При этом в случае входящего потока задаются также значения для остальных консервативных величин, в случае выходящего потока для них используется условие d Ф/dn = 0.

При моделировании пожаров часто встречаются участки границы, через которые возможно течение газовой среды как внутрь расчетной области, так и из нее (дверные и оконные проемы, люки дымоудаления и т.п.). Используемые на таких границах граничные условия можно разделить на два типа: условия с заданной нормальной скоростью и условия с заданным давлением. В условиях первого типа значение скорости задается не явно, а, в виде условий типа dv n /dn = 0 или d 2 v n /dn 2 = 0. Значение давления на границе при этом определяется из решаемых уравнений. В условиях второго типа давление может задаваться как в явном виде, так и в форме dp /dn = 0. При этом величина нормальной скорости вычисляется с использованием значения давления. Для касательных компонент скорости и в том и в другом случае обычно используются условия dv /dn = 0.

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

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

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

Стандартная k -e модель турбулентности не позволяет смоделировать этот эффект, поэтому использование в качестве граничных условий значений k и e, соответствующих параметрам потока горючего, приводит к занижению значений турбулентной вязкости в области пламени и, как следствие, к завышению значений скоростей и температур в области пламени и восходящей свободно-конвективной струи . Строгого решения задачи о постановке этих граничных условий на данный момент не существует. Для практических расчетов в качестве граничных условий используют искусственные значения k и e , обеспечивающие разумную величину турбулентной вязкости в области пламени без рассмотрения процессов, протекающих в тонком слое вблизи поверхности горючего. Так, проведенные исследования показали, что хорошие результаты при использовании k -e модели в сочетании с диффузионно-вихревой моделью горения дает использование значений k = 0,3 м 2 /с 2 и e = 1×10 -6 м 2 /c 3 .

5. ПОРЯДОК ПРОВЕДЕНИЯ РАСЧЕТНОЙ ОЦЕНКИ ПОЖАРНОЙ ОПАСНОСТИ КОНКРЕТНОГО ОБЪЕКТА

Порядок проведения расчетной оценки пожарной опасности конкретного объекта в виде блок-схемы представлен на рис. 1.

Сбор исходных данных включает в себя изучение:

объемно-планировочных решений объекта;

теплофизических характеристик ограждающих конструкций и размещенного на объекте оборудования;

вида, количества и расположения горючих материалов;

количества и вероятного расположения людей в здании;

материальной и социальной значимости объекта;

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

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

вероятность возникновения пожара;

возможная динамика развития пожара;

наличие и характеристики систем противопожарной защиты (СППЗ);

вероятность и возможные последствия воздействия пожара на людей, конструкцию здания и материальные ценности;

соответствие объекта и его СППЗ требованиям противопожарных норм.

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

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

Рис. 1. Порядок проведения расчетной оценки пожарной опасности объекта

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

Приложение

ПРИМЕР РАСЧЕТА

Характеристика объекта

Рассматриваемое пятиэтажное здание II степени огнестойкости является многофункциональным комплексом и включает в себя спальную зону с номерами, административно-бытовую часть и учебные помещения. Пожарная нагрузка представлена офисной и бытовой мебелью, оргтехникой, горючими материалами отделки помещений. В здании одновременно могут находиться 255 человек, которые распределены по этажам следующим образом: на 1-м этаже 34 человека; на 2-м - 48; на 3-м - 96; на 4-м - 59; на 5-м - 18 человек.

Система противопожарной защиты представлена:

тепловыми пожарными извещателями;

незадымляемыми лестничными клетками;

системой оповещения о пожаре 2-го типа;

внутренним пожарным водопроводом и первичными средствами пожаротушения.

Качественный анализ пожарной опасности объекта

С точки зрения пожарной опасности особенностями рассматриваемого объекта являются:

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

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

наличие эвакуационных путей через галереи и помещения, открытые в объем атриума;

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

возможность массового присутствия людей в одном помещении.

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

Выбор критерия пожарной опасности

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

Выбор сценария пожара

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

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

Очаг пожара возникает в бильярдной на первом этаже. Пламя распространяется по мебели (бильярдный стол, кресло, открытый шкаф). Максимальная площадь горящей поверхности 5,2 м 2 , максимальная мощность пожара 2 МВт. Динамика развития очага пожара определяется характерной скоростью распространения фронта пламени по горизонтали 3 см/с и по вертикальным поверхностям - 0,1 см/с и охватывает всю поверхность горючих материалов за 120 с.

Рис. 2. Схема системы помещений

Формулировка математической модели

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

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

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

Результаты моделирования

Первоначально развитие пожара происходит в пределах помещения очага пожара (бильярдной). К моменту времени 30 с происходит заполнение дымом верхней части помещения очага и начинается выход продуктов горения через открытый дверной проем (двустворчатая дверь 2´1,7 м), а через нижнюю часть проема в помещение поступает воздух, поддерживающий горение. Далее происходит выход продуктов горения в объем атриума (рис. 3) и их растекание под галереей 2-го этажа.

Рис. 3. Поля температур (К) в вертикальном сечении атриума в момент времени 90 с

Образуется плоская конвективная колонка, поднимающаяся к потолку атриума. К моменту времени 90 с струя продуктов горения поднимается до уровня 4-го этажа. Задымление галерей 2-го и 3-го этажей при этом не происходит. В то же самое время продолжается растекание продуктов горения под галереей 2-го этажа. К моменту времени 120 с конвективная колонка достигает потолка атриума и начинается радиальное растекание продуктов горения (рис. 4, а ). При этом происходит задымление ближайшей к колонке части галереи 5-го этажа и блокирование одного из эвакуационных выходов (рис. 4, в ).

Рис. 4. Поля температур (К) в вертикальном сечении атриума (а), горизонтальном сечении под потолком 1-го этажа (б) и сечении на уровне 1,7 м от пола 5-го этажа в момент времени 120 с

К моменту времени 180 с продукты горения в объеме атриума опускаются до уровня 2-го этажа (рис. 5). При этом происходит полное задымление галереи 5-го этажа и блокирование обоих эвакуационных выходов на 4-м этаже. На третьем этаже (рис. 6, а ) большая часть галереи остается свободной от дыма и блокируется только один эвакуационный выход. Задымление на 2-м этаже (рис. 6, б ) на уровне 1,7 м незначительно, и все эвакуационные выходы свободны. Эвакуационные выходы на первом этаже остаются свободными. К моменту времени 240 с происходит опускание дымовых газов до пола первого этажа и полное блокирование эвакуационных выходов на всех этажах (рис. 7).

5-й этаж - t 5,1 = 120 с; t 5,2 = 180 с;

4-й этаж - t 4,1 = 180 с; t 4,2 = 180 с;

3-й этаж - t 3,1 = 180 с; t 3,2 = 240 с;

2-й этаж - t 2,1 = 240 с; t 2,2 = 240 с; t 2,3 = 240 с;

1-й этаж - t 1,1 = 240 с; t 1,2 = 240 с.

Сравнение расчетных значений критерия пожарной опасности с критическими значениями

Таким образом, в результате проведенного расчета получены количественные значения критерия оценки пожарной опасности. Данные значения необходимо сравнить с критическими, а именно со значениями времени эвакуации людей, полученными согласно методике ГОСТ 12.1.004-91*, приложение 2, п. 2.4. Значения расчетного времени эвакуации и времени блокирования эвакуационных путей для каждого этажа здания приведены в табл. 1.

Таблица 1

Наименование участка эвакуации	Количество людей, чел.	Расчетное время эвакуации t p , с	Время блокирования эвакуационных путей t бл , с	Выполнение условия t p £ t бл
Первый этаж				Выполняется
Второй этаж				Выполняется
Третий этаж				Выполняется
Четвертый этаж				Выполняется
Пятый этаж				Выполняется

Сравнение значений, приведенных в таблице, показывает, что условия безопасной эвакуации людей выполняются.

Анализ правильности выбора сценария

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

Заключение о пожарной опасности объекта

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

ЛИТЕРАТУРА

1. ГОСТ 12.1.004-91* Пожарная безопасность. Общие требования.

2. СНиП 21-01-97* Пожарная безопасность зданий и сооружений.

3. Расчет необходимого времени эвакуации людей из помещений при пожаре: Рекомендации. - М.: ВНИИПО МВД СССР, 1989. - 22 с.

4. Рыжов A . M . Моделирование пожаров в помещениях с учетом горения в условиях естественной конвекции // Физика горения и взрыва. - 1991. - Т. 27, № 3. - С. 40-47.

5. Computer modelling of aerodynamics and movement of aerosol in volumes of complex geometry / L.P. Kamenshchikov, V.I. Bykov, S.P. Amel"chugov, A.A. Dekterev // Proc. of the 2 nd Int. Seminar on Fire and Explosion Hazard of Substances and Venting of Deflagrations. Moscow, 1997. - P. 512-521.

6. Cox G., Kumar S. Field Modelling of Fire in Forced Ventilated Enclosures // Comb. Science and Tech. - 1987. - Vol. 52. - P. 7-23.

7. Lewis M.J., Moss M.B . and Rubini P.A. (1997) CFD modelling of combustion and heat transfer in compartment fires // Proc. of V Int. Symp. On Fire Safety Science. - P. 463-474.

8. Патанкар С. Численные методы решения задач теплообмена и динамика жидкостей. - М.: Энергоатомиздат, 1984. -150 с.

9. Провести исследования и разработать методические рекомендации применения фундаментального полевого метода моделирования динамики развития пожаров и распространения их опасных факторов в помещениях зданий различного назначения: Отчет о НИР (аннот.) // ВНИИПО МВД России. -П.3.4.Д.002.2001; Код "Фундамент". - Этап 1. - М., 2001. - 51 с.

10. Провести фундаментальные исследования процесса развития пожара внутри и вне помещений и зданий различного назначения с использованием методов вычислительной гидродинамики, изучить закономерности процесса и сформулировать предложения в НПБ: Отчет о НИР (заключ.) // ВНИИПО МВД России. - П.3.4.Д.001.98, Код "Закономерности". - М., 2000. - 144 с.

11. Сох G. Combustion Fundamentals of Fire. - London: Academic Press, 1995. - 476 p.

12. Baum H.R., McGrattan K.B., Rehm R.G. Three dimensional simulations of fire plume dynamics // Proc. of V Int. Sump. "Fire Safety Science", 1997. - P. 511-522.

13. Magnussen B.F. and Hjertager B.H. (1976) On mathematical modelling of turbulent combustion with special emphasis on soot formation and combustion. 16 th Sump. (Int.) Combust. The Combustion Institute. - Pittsburgh, PA. - P. 719-729.

14. Peters N. (1986) Laminar flamelet concept in turbulent combustion. 21 th Symp. (Int.) Combust. The Combustion Institute. - Pittsburgh, PA. - P. 1231-1250/

15. Patankar S.V. and Spalding D.B. (1973) A computer model for three-dimensional flow in furnaces. 14 th Symp. (Int.) Combust. The Combustion Institute. - Pittsburgh, PA. - P. 605-614.

16. Tuovinen H. (1994) Modelling of laminar diffusion flames in vitiated environment, Proc. of IV Int. Symp. on Fire Safety Science. - P. 113-124.

17. Lockwood F.C. and Shah N.G. (1981) A new radiation solution method for incorporation in general combustion prediction procedures. 18 th Symp. (Int.) Combust. The Combustion Institute. -Pittsburgh, PA. - P. 1405-1414.

18. Методы расчета температурного режима пожара в помещениях зданий различного назначения: Рекомендации. - М.: ВНИИПО МВД СССР, 1988. - 56 с.

19. Термогазодинамика пожаров в помещениях / В.М. Астапенко, Ю.А. Кошмаров, И.С. Молчадский, А.Н. Шевляков . - М.: Стройиздат, 1988. - 448 с.

20. Белов И.А., Исаев С.А., Коробков В.А. Задачи и методы расчета отрывных течений несжимаемой жидкости. - Л.: Судостроение, 1989. - 150 с.

21. Jayatillake C.L.V. The influence of Prandtl number and surface roughness on the resistance of laminar sublayer to momentum and heat transfer // Progress in Heat and Mass Transfer. - 1969. - № 1. - P. 193-329.

22. Tuovinen H. (1997) CFD modelling of underventilated fires // Proc. of the 2 nd Int. Seminar on Fire and Explosion Hazard of Substances and Venting of Deflagrations, Moscow, 1997. - P. 113-124.

23. Weckman E.J. and Strong A.B. Experimental Investigation of the Turbulence Structure of Medium Scale Methanol Pool Fires // Combustion and Flame. - 1996. - Vol. 105, № 3. - P. 245-266.

24. Карпов А.В., Крюков А.П., Рыжов A . M . Полевое моделирование процессов тепло- и массопереноса в пламени и восходящей свободноконвективной струе //Пожаровзрывобезопасность. - 2001. - Т. 10, № 2. - С. 35-41.

25. Modelling thermal radiation in open liquid pool fires /K.C. Adiga, D.E. Ramaker, PA. Tatem, F.W. Williams // Proc. of III Int Symp. on Fire Safety Scince. - 1989. - P. 241-250.

26. Turbulent diffusion flames with large buoyancy effects E. Gengembre, P. Cambray, D. Karmed and J.C. Bellet // Combustion Science and Technology. - 1984. - Vol. 41. - P. 55-67.

27. Modelling Buoyant Turbulent Diffusion Flames in Coherent Flame-sheet model / С .A. Blunsdon, Z. Beeri, W.M.G. Malalesekera, J.C. Dent // Symposium on Fire and Combustion, ASME Winter Annual Meeting Chicago: ASME. - 1994. - P. 79-88.

28. Welch S., Rubini P . SOFIE, Simulations of Fires in Enclosures, User Guide. - Cranfield University, 1996.

Список обозначений

Введение

1. Общие положения

2. Область применения

3. Основы полевого метода моделирования пожаров

3.1. Основные уравнения

3.2. Моделирование турбулентности

3.3. Модели горения

3.4. Радиационный теплоперенос

3.4.1. Потоковые методы

3.4.2. Метод дискретного радиационного переноса

4. Замыкание основной системы уравнений.

Условия однозначности

4.1. Граничные условия на твердых негорючих поверхностях

4.2. Граничные условия на плоскости (оси) симметрии

4.3. Граничные условия, характеризующие работу приточно-вытяжной вентиляции

4.4. Граничные условия на свободной границе

4.5. Граничные условия на поверхности горючего

5. Порядок проведения расчета при оценке пожарной опасности конкретного объекта

Приложение. Пример расчета

Интегральная математическая модель пожара представляет собой систему обыкновенных дифференциальных уравнений, описывающих изменение среднеобъёмных параметров состояния газовой среды в помещении в процессе развития пожара. Они следуют из фундаментальных законов природы? первого закона термодинамики для открытой термодинамической системы и закона сохранения массы. Впервые интегральная модель была сформулирована профессором Ю.А. Кошмаровым в 1976 году.

Более подробно интегральная модель пожара описана в приложении 6 к приказу МЧС России от 30.06.2009 №382.

Ограничения интегральной модели

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

достаточно большой источник пожара;

относительно небольшой объем помещений;

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

Таким образом, интегральную модель можно применять при следующих условиях:

для зданий, содержащих развитую систему помещений малого объема простой геометрической конфигурации;

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

для предварительных расчетов с целью выявления наиболее опасного сценария пожара.

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

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

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

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

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

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

Полевые модели основываются на системе дифференциальных уравнений в частных производных. Результатами решения данной системы уравнений являются поля распределения температур, скоростей, концентраций компонентов газовой среды в каждый момент времени. Программа FDS (Fire Dynamics Simulator) реализует вычислительную гидродинамическую модель (CFD) тепломассопереноса при горении. FDS решает уравнения Навье-Стокса для низкоскоростных температурно-зависимых потоков. Базовым алгоритмом является определенная схема использования метода предиктора-корректора второго порядка точности по координатам и времени.

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

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

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

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

опасный пожар прогнозирование моделирование

САНКТ-ПЕТЕРБУРГСКИЙ УНИВЕРСИТЕТ ГПС МЧС РОССИИ ______________________________ _________________________ Кафедра Правового и Кадрового обеспечения КОНТРОЛЬНАЯ РАБОТА по курсу: «Физико-химические основы развития и тушения пожаров» Тема: Исходные понятия и общие сведения об опасных факторах пожара и методах их прогнозирования. Выполнил: студент института заочного и дистанционного обучения Гр. № 508 зачетная книжка № в-0876 специальность 280104.65 Габдуллин Динар Дамирович Санкт-Петербург 2011г.

Содержание

Введение………………………………………………………… ………………	3стр.
Опасные факторы пожара ………………………………………… ………...	4стр.
Пламя как опасный фактор пожара……………………………………………	4стр.
Искры как опасный фактор пожара……………………………………………	4стр.
Повышенная температура как опасный фактор пожара…………………….	5стр.
Дым как опасный фактор пожара……………………………………………...	5стр.
Пониженная концентрация кислорода как опасный фактор пожара………..	5стр.
Концентрация токсичных веществ как опасный фактор пожара…………...	5стр.
Разрушение конструкций как опасный фактор пожара……………………..	6стр.
Отравление угарным газом как опасный фактор пожара……………………	6стр.
Методы прогнозирования пожара …………………………………………..	7стр.
Классификация интегральных математических моделей пожара…………...	7стр.
Интегральная модель пожара…………………………………………………..	9стр.
Зонная модель пожара………………………………………………………….	9стр.
Полевой (дифференциальный) метод расчета………………………………..	11стр
Критерии выбора моделей пожара для расчетов……………………………..	12стр
Заключение…………………………………………………… …………………	13стр
Список использованной литературы…………………………………………..	14стр

Введение

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

Опасные факторы пожара

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

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

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

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

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

Пламя как опасный фактор пожара

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

Искры как опасный фактор пожара

Самое частое и, вместе с тем банальное - это когда «из искры возгорится пламя»: здесь враг виден, если можно так выразиться - в лицо. Маленькая искра, перерастающая в открытое пламя - и, как следствие, большие неприятности: лесные и степные пожары, пожары в сельскохозяйственных и промышленных постройках, административных зданиях, жилых помещениях, движимом имуществе. Как правило, огромные материальные убытки. Однако что касается людей, то открытый огонь на них редко воздействует, людей поражают преимущественно испускаемые пламенем лучистые потоки, поражающие открытые участки тела. Весьма опасны ожоги от горящей одежды, особенно из синтетических тканей, которая трудно тушится и так же трудно сбрасывается.

Повышенная температура как опасный фактор пожара

Следующий фактор пожара - повышенная температура окружающей среды - может как усугубить действие предыдущего, так и выступить самостоятельным источником материальных убытков и физических страданий людей, вызванных пожаром от самовозгорающихся предметов и материалов. Наибольшая опасность для людей исходит от нагретого воздуха, который при вдыхании, обжигает верхние дыхательные пути и приводит к удушью и смерти. К летальному исходу приводит и вызванный этим фактором пожара перегрев, из-за чего из организма интенсивно выводятся соли, нарушается деятельность сосудов и сердца. Достаточно побыть несколько минут в среде с температурой в 100 °С - как сразу же теряется сознание и наступает смерть. Вместе с тем, губительное влияние на человека оказывает и продолжительное облучение инфракрасными лучами с интенсивностью около 540 Вт/м. Также при повышенной температуре окружающей среды часты ожоги кожи.

Дым как опасный фактор пожара

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

Пониженная концентрация кислорода как опасный фактор пожара

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

Концентрация токсичных веществ как опасный фактор пожара

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

Разрушение конструкций как опасный фактор пожара

Разрушение конструкций это еще один из опасных факторов пожара приводящих к травмам увечьям и гибели людей находящихся в зоне разрушения.
В первые 10-20 минут пожар распространяется вдоль горючего материала и в это время помещение заполняется дымом. Температура воздуха поднимается в помещении до 250-300 градусов. Через 20 минут начинается объёмное распространение пожара.
Спустя ещё 10 минут наступает разрушение остекления. Увеличивается приток свежего воздуха, резко прогрессирует развитие пожара и температура достигает 900 градусов.
После того, как выгорают основные вещества, конструкция здания теряет свою несущую способность и в это время происходит обрушение выгоревших конструкций.

Отравление угарным газом как опасный фактор пожара

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

Методы прогнозирования пожара

Классификация интегральных математических моделей пожара

Современные научные методы прогнозирования Опасных Факторов Пожара основываются на математическом моделировании, т.е. на математических моделях пожара. Математическая модель пожара описывает в самом общем виде изменение параметров состояния среды в помещении в течение времени, а также изменение параметров состояния ограждающих конструкций этого помещения и различных элементов технологического оборудования.
Методы прогнозирования ОФП различают в зависимости от вида математической модели пожара. Математические модели пожара в помещении условно делятся на три класса (три вида): интегральные, зонные, полевые (дифференциальные).
1. Интегральная модель пожара позволяет получить информацию, т.е. сделать прогноз, о средних значениях параметров состояния среды в помещении для любого момента развития пожара. При этом для того, чтобы сопоставлять (соотносить) средние (т. е. среднеобъемные) параметры среды с их предельными значениями в рабочей зоне, используются формулы, полученные на основе экспериментальных исследований пространственного распределения температур, концентраций продуктов горения, оптической плотности дыма и т. д.
2. Зонная модель позволяет получить информацию о размерах характерных пространственных зон, возникающих при пожаре в помещении, и средних параметров состояния среды в этих зонах. В качестве характерных пространственных зон можно выделить, например, припотолочную область пространства, в начальной стадии пожара, область восходящего над очагом горения потока нагретых газов и область незадымленной холодной части пространства.
3. Полевая дифференциальная модель позволяет рассчитать для любого момента развития пожара значения всех локальных параметров состояния во всех точках пространства внутри помещения.
Перечисленные модели отличаются друг от друга объемом той информации, которую они могут дать о состоянии газовой среды в помещении и взаимодействующих с нею конструкций на разных этапах (стадиях) пожара. В этом отношении наиболее детальные сведения можно получить с помощью полевой модели.
В математическом отношении три вышеназванных вида моделей пожара характеризуются разным уровнем сложности.
Интегральная модель пожара в своей основе представлена системой обыкновенных дифференциальных уравнений. Искомыми функциями выступают среднеобъемные параметры состояния среды, независимым аргументом является время.
Основу зонной модели пожара в общем случае составляет совокупность нескольких систем обыкновенных дифференциальных уравнений. Параметры состояния среды в каждой зоне являются искомыми функциями, а независимым аргументом является время. Искомыми функциями являются также координаты, определяющие положение границ характерных зон.
Наиболее сложной в математическом отношении является полевая модель. Ее основу составляет система уравнений в частных производных, описывающих пространственно-временное распределение температур и скоростей газовой среды в помещении, концентраций компонентов этой среды (кислород, оксид и диоксид углерода и т.д.), давлений и плотностей. Эти уравнения включают реологический закон Стокса, закон теплопроводности Фурье, закон диффузии, закон радиационного переноса и т.п. В более общем случае к этой системе уравнений добавляется дифференциальное уравнение теплопроводности, описывающее процесс нагревания ограждающих конструкций. Искомыми функциями в этой модели являются плотность и температура среды, скорость движения газа, концентрации компонентов газовой среды, оптическая плотность дыма (натуральный показатель ослабления света в дисперсной среде) и т.д. Независимыми аргументами являются координаты х, у, z и время т.

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

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

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

Интегральная модель пожара

Интегральная математическая модель пожара описывает в самом общем виде процесс изменения во времени состояния газовой среды в помещении.
С позиций термодинамики газовая среда, заполняющая помещение с проемами (окна, двери и т.п.), как объект исследования есть открытая термодинамическая система. Ограждающие конструкции (пол, потолок, стены) и наружный воздух (атмосфера) является внешней средой по отношению в этой термодинамической системе. Эта система взаимодействует с внешней средой путем тепло- и массообмена. В процессе развития пожара через одни проемы выталкивается из помещения нагретые газы, а через другие поступает холодных воздух. Количество вещества, т.е. масса газа в рассматриваемой термодинамической системе, в течении времени изменяется. Поступление холодного воздуха обусловлено работой проталкивания, которую совершает внешняя среда. Термогазодинамическая система в свою очередь совершает работу, выталкивая нагретые газы во внешнюю атмосферу. Эта термодинамическая система взаимодействует также с ограждающими конструкциями путем теплообмена. Кроме того, в эту систему с поверхности горящего материала (т.е. из пламенной зоны) поступает вещество в виде газообразных продуктов горения.
Состояние рассматриваемой термодинамической системы изменяется в результате взаимодействия с окружающей средой. В интегральном методе описания состояния термодинамической системы, коей является газовая среда в помещении, используются «интегральные» параметры состояния – такие, как масса всей газовой среды и ее внутренняя тепловая энергия. Отношение этих двух интегральных параметров позволяет оценивать в среднем степень нагретости газовой среды. В процесс развития пожара, значения указанных интегральных параметров состояния изменяются.

Зонная модель пожара

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

Рисунок 1

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

Полевой (дифференциальный) метод расчета

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

Рисунок 2

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

Критерии выбора моделей пожара для расчетов

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

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

зональный метод:

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

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

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

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

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

ЛЕКЦИЯ

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

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

План лекции:

Введение

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

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

Цели лекции:

1. Учебные

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

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

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

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

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

Литература

1. Д.М. Рожков. Прогнозирование опасных факторов пожара в помещении. – Иркутск 2007 С.89
2. Ю.А.Кошмаров, М.П. Башкирцев Термодинамика и теплопередача в пожарном деле. ВИПТШ МВД СССР, М., 1987 г.
3. Лабораторный практикум «Прогнозирование опасных факторов пожара». Ю.А.Кошмаров, Ю.С.Зотов. 1997 г.
4. Ю.А.Кошмаров, В.В. Рубцов, Процессы нарастания опасных факторов пожара в производственных помещениях и расчет критической продолжительности пожара. МИПБ МВД России, М., 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
	Это требует учета при моделировании корреспонденций ограниченных возможностей районов по размещению прибывающих автомобилей. Величина этих затрат должна зависеть от соотношения объема автомобильных прибытий в район и его возможностей принять такой объем и увеличиваться по мере нарастания объема прибытий. в районах прибытия для автомобильного транспорта возникают такие дополнительные затраты при которых в ходе расслоения корреспонденций формируются объемы прибытий на автомобилях порождающие именно эти величины дополнительных затрат. Для...