5.1. Схема трехзонной модели пожара:
I - зона конвективной струи (конвективная колонка);
II - зона припотолочного нагретого газа; III - зона холодного
воздуха; IV - зона наружного воздуха (наружная атмосфера)
Рассмотрим прежде всего I зону. Теория свободной конвективной струи к настоящему времени весьма детально разработана. Эта теория является одним из разделов вязкой аэродинамики газов. Она позволяет
Рассчитывать поля температур, плотностей и скоростей в конвективной колонке. Для определения температур и массовых расходов в сечениях конвективной колонки можно использовать формулы [4]
(5.1)
(5.2)
где ; Qпож - скорость тепловыделения, Вт; Qрн -теплота сгорания, Дж/кг; уд - удельная скорость выгорания, кг/м2•с; g -ускорение свободного падения, м/с2; Т0 и 0 - температура и плотность холодного (окружающего) воздуха; G - расход газов через сечение струи, отстоящее от поверхности горения на расстояние у, кг/с; ср - изобарная теплоемкость газа, Дж/кг-К; - доля, приходящаяся на поступающую в ограждение теплоту от выделившейся в очаге горения; у - координата сечения колонки, отсчитываемая от поверхности горения, м; у0 - расстояние от фиктивного источника тепла до поверхности горения, м.
С помощью формул (5.1) и (5.2) можно рассчитать расход газа из I зоны, поступающего во II зону, и «То температуру. Для этого нужно положить координату у в формулах (5.1) и (5.2) равной координате нижней границы припотолочного слоя ук
Расстояние от фиктивного источника тепла до поверхности горения вычисляется по формуле [4]
(5.3)
где Fr - площадь пожара, м2.
Рассмотрим теперь II зону (припотолочный слой нагретых газов). Объем этой зоны в момент времени т равен
,
где Fпот - площадь потолка; уk - координата нижнего края припотолочного слоя газов. Масса газа, заключенная во II зоне, составляет величину m2=2V2. Давление в зоне II практически не меняется и остается равным начальному значению, т.е. р0. Внутренняя (тепловая) энергия II зоны составляет
Запишем уравнения материального баланса и энергии для II зоны применительно к первой фазе начальной стадии пожара:
(5.4)
(5.5)
где 2 - средняя плотность во II зоне; Т2 - средняя температура во II зоне;
Qw2 - тепловой поток от припотолочного слоя газа в ограждения, кВт.
Параметры состояния T0 и р2 связаны между собой следующим уравнением:
(5.6)
Уравнение (5.6) следует из условия равенства давлений во всех зонах. Это условие является приближенным, но применимым для реальных пожаров.
Преобразуем уравнение энергии (5.5), используя уравнение (5.6),
или
и окончательно (5.7)
Из уравнения (5.1) следует
(5.8)
где Qwl = χQпож
Подставляя формулу (5.8) в уравнение (5.7), получим
Примем, что const (для начальной стадии 0,55).
После дальнейших преобразований получим следующее уравнение:
(5.8а)
Подставим в это уравнение выражение для Gk (5.2):
(5.9)
Отметим, что в этом уравнении
Введем обозначения
Функции () и () при горении твердых ГМ в момент времени = 0 равны нулю, так как Fr → 0.
Уравнение (5.9) принимает вид
(5.10)
Начальное условие yk(=0)=2h-δ.
Решение уравнения (5.10) при заданном начальном условии будем искать для интервала времени от т = 0 до * , где *, - момент окончания первой фазы начальной стадии пожара. После того как найдена функция ук(), находим
Обратимся к уравнению материального баланса (5.4). Интегрируя его, получаем
(5.11)
После преобразований из формулы (5.11) получаем
(5.12)
После вычислений плотности р2 определяется средняя температура в при потолочном слое газа
(5.13)
Уравнение материального баланса для токсичного газа (продукт горения) во II зоне имеет вид
(5.14)
где п - парциальная плотность токсичного газа; L - количество (масса) токсичного газа, образующаяся при сгорании 1 кг горючего материала.
Из формулы (5.14) следует формула
(5.15)
где М - количество (масса) ГМ, выгоревшего к моменту времени .
Уравнение дыма имеет вид
(5.16)
и, следовательно,
(5.17)
Рассмотрим частный случай, когда
Qпож=удQрнFг=const; Fг=const; δ=0
В этом случае будем иметь
где
Далее следует отметить, что y0=const;
Следовательно, имеем уравнение с разделяющимися переменными
С помощью этого уравнения рассчитывается изменение координаты границы припотолочного слоя в течение времени.
- Глава 1. Интегральная математическая модель пожара в помещении................14
- Глава 6. Дифференциальные (полевые) математические модели пожара в помещении..................................................................................................................88
- Введение. Общие сведения о методах прогнозирования опасных факторов пожара в помещении
- Глава 1. Интегральная математическая модель пожара в помещении
- Исходные положения и основные понятия интегрального метода термодинамического анализа пожара
- 1.2. Дифференциальные уравнения пожара
- Глава 2. Дополнительные уравнения интегральной математической модели пожара для расчета расходов уходящих газов и поступающего через проемы воздуха
- 2.1. Исходные положения
- 2.2. Распределение давлений по высоте помещения
- 2.3. Плоскость равных давлений и режимы работы проема
- 2.4. Распределение перепадов давлений по высоте помещения
- 2.5. Формулы для расчета расхода газа, выбрасываемого через прямоугольный проем
- 2.6. Формулы для расчета расхода воздуха, поступающего через прямоугольный проем
- 2.7. Влияние ветра на газообмен
- Глава 3. Дополнительные уравнения интегральной модели пожара для расчета теплового потока в ограждения и скорости выгорания горючих материалов
- 3.1. Приближенная оценка величины теплового потока в ограждения
- 3.2. Эмпирические методы расчета теплового потока в ограждения
- 3.3. Полуэмпирические методы расчета теплового потока в ограждения
- 3.4. Методы расчета скорости выгорания горючих материалов и скорости тепловыделения
- Глава 4. Математическая постановка и методы решения задачи о прогнозировании офп на основе интегральной математической модели пожара в помещении
- 4.1. Классификация интегральных моделей пожара
- 4.2. Интегральная математическая модель пожара для исследования динамики офп и ее численная реализация
- 4.3. Интегральная математическая модель начальной стадии пожара и расчет критической продолжительности пожара
- 4.3.1. Постановка задачи и ее решение
- 4.3.2. Расчет критических значений средних параметров состояния среды в помещении
- 4.3.3. Расчет коэффициента теплопоглощения (коэффициента
- Глава 5. Зонная математическая модель пожара в помещении
- 5.1. Схема трехзонной модели пожара:
- Глава 6. Дифференциальные (полевые) математические модели пожара в помещении математическая модель расчета тепломассообмена при пожаре в помещении
- 6.1. Особенности и упрощения термогазодинамической картины пожара
- 6.2.Структура полевой модели расчета тепломассообмена
- Основные уравнения
- 6.3. Основные уравнения полевой модели
- 6.4. Уравнения для расчета процесса прогрева строительных конструкций
- 6.5. Расчет турбулентного тепломассообмена
- 6.5.6. Уравнения (6.17)-(6.23) позволяют определить коэффициенты турбулентной вязкости, теплопроводности и диффузии, входящие в уравнения полевой модели (6.2)-(6.6).
- 6.6. Моделирование радиационного теплообмена
- 6.7. Расчет процесса выгорания горючей нагрузки
- 6.8. Моделирование горения
- 6.9. Условия однозначности
- 6.10. Моделирование действий систем пожаротушения
- 6.11. Моделирование действий систем механической вентиляции и дымоудаления
- 6.12. Метод численного решения дифференциальных уравнений
- Заключение
- Литература
- 129366, Москва, ул. Б. Галушкина, 4