4.3.3. Расчет коэффициента теплопоглощения (коэффициента
Как уже указывалось, коэффициент теплопотерь представляет собой отношение суммарного теплового потока в ограждения Qw к скорости тепловыделения в очаге горения Qпож, т.е.
(4.70)
Скорость тепловыделения в каждый момент процесса развития пожара вычисляется по формуле
(4.71)
где FT - площадь пожара, м² .
Суммарный тепловой поток в ограждения (выполненные из кирпича или из материалов, близких ему по своим теплофизическим свойствам) при значениях среднеобъемной температуры среды в помещении T0 < Тт < Ткр (где Т0290 К и Гкр343 К) можно рассчитать с помощью эмпирической формулы [2]:
(4.72)
где q0, a, b1 - размерные эмпирические константы (q0=4,07 Вт•м2; а = 0,8 К-1; b1=0,00065 К-2); Fw - суммарная площадь поверхности ограждений, м2.
Для помещений, представляющих собой прямоугольный параллелепипед, суммарная площадь поверхностей ограждений вычисляется по формуле
(4.73)
где l1, l2 - ширина и длина помещения, м; h - половина высоты помещения, м. При этом предполагается, что площадь проемов пренебрежимо мала по сравнению с величиной Fw. Кроме того, при использовании этой формулы предполагается возможным не учитывать наличие предметов и оборудования, находящихся внутри помещения .
Следует сделать замечание по поводу формулы (4.72). При Тт → Т0 (т.е. в первый момент процесса развития пожара при t→ 0) из этой формулы следует, что Q*w → 0. Однако в первый момент времени поступление тепла в ограждения происходит главным образом за счет лучеиспускания от пламени. Поток лучистой энергии от пламени к ограждениям по мере задымления помещения постепенно уменьшается. При сильном задымлении излучение от пламени рассеивается в задымленной среде, заполняющей помещение. С учетом сказанного радиационный поток от очага горения к ограждениям при FRAD << Fw можно оценить по формуле [2]
(4.74)
где С = 5,7 Вт•м-2 К-4 - коэффициент излучения; г- степень черноты пламени; FRAD - площадь поверхности излучения, м2; Гпл - температура пламени, К; со - коэффициент, учитывающий ослабление радиационного потока из-за задымления.
Коэффициент ω в начальные моменты времени, когда среда еще достаточно прозрачна, равен единице. При сильном задымлении, которое может иметь место в конце начальной стадии пожара, этот коэффициент равен нулю. С учетом сказанного величину этого коэффициента можно приближенно оценить с помощью формулы
(4.75)
Эта формула представляет собой линейную интерполяцию зависимости ω(Тт) в интервале температур Т0 < Тт < Ткр, т.е. в интервале времени, равном критической продолжительности пожара.
Из вышесказанного следует, что для расчета суммарного теплового потока в ограждения при начальной стадии пожара нужно использовать формулу
(4.76)
I4
С учетом вышеизложенного значение коэффициента теплопотерь следует вычислять для каждого момента времени в начальной стадии пожара (т.е. при 0 <<ткр) по формуле
(4.77)
Из этой формулы следует, что величину ср, лишь условно можно назвать коэффициентом, так как она является функцией Тт() и изменяется во времени.
В интервале времени, равном критической продолжительности пожара, среднеобъемная температура среды в помещении не сильно отличается от "предельного" значения Тпред приблизительно равного 345 К. Если температура среды перед пожаром Т0 = 293 К (что является типичным начальным условием), то в начальной стадии пожара выполняется условие
(4.78)
С учетом этой оценки правую часть уравнения (4.77) можно упростить, отбрасывая второй член первого слагаемого в квадратных скобках. После этой операции получим уравнение
А7»и
(4.79)
Среднее значение коэффициента теплопотерь в интервале температур T0<Tт<Tкр , т.е. в интервале времени, равном критической продолжительности пожара, определяется путем операции осреднения значения коэффициента * , т.е.
(4.80)
где ΔТт = Тт -Т0; ΔТкр = Ткр-Т0; - средний коэффициент теплопотерь.
Далее отдельно рассмотрим пожары в помещении при горении ТГМ и ГЖ. Вначале дается анализ пожаров в помещениях при горении ТГМ. При оценке радиационной составляющей теплового потока от пламени можно принять, что площадь поверхности излучения FRAD равна площади пожара Fг, т.е. FRAD = Fг.
4/
еС\
Прежде всего, рассмотрим пожары при круговом распространении пламени по слою ТГМ. Для таких пожаров в помещении уравнение (4.79) преобразуется с помощью формулы (4.43):
(4.81)
Из формулы (4.52) следует, что
(4.82)
где - среднее значение коэффициента теплопотерь.
«1
При начальной стадии развития пожара выполняется условие
Из этого условия следует, что
(4.83)
С учетом формул (4.82) и (4.83) уравнение (4.81) преобразуется в следующую формулу:
(4.84)
где *, - коэффициент теплопотерь при температуре среды Tm(); -средний коэффициент теплопотерь в интервале времени, равном критической продолжительности пожара.
Поставим полученное выражение для * в формулу (4.80). После интегрирования получается уравнение для расчета среднего коэффициента теплопотерь, которое можно представить в виде
(4.85)
uде ; ;
; .
Безразмерный комплекс Г характеризует макрокинетику горения ТГМ. Безразмерный комплекс Ф есть обобщенная геометрическая характеристика помещения. Этот комплекс можно назвать критерием формы помещения. Если помещение имеет формулу куба, то Ф = 6. Для помещений, форма которых отличается от куба, Ф > 6.
При вычислении суммарной площади поверхности ограждений Fw (и следовательно, критерия формы Ф) можно не учитывать наличие проемов, потому что здесь рассматриваются пожары в помещениях с
F
относительно малыми проемами (т.е. помещения, у которых 0).
Безразмерный комплекс Δrad является параметром влияния радиационного теплообмена и представляет собой отношение тепла,
теряемого из-за излучения единицей площади поверхности ТГМ, охваченной пламенем, к теплу, выделяющемуся на этой единичной площадке вследствие горения. Для большинства ТГМ значение параметра влияния радиационного теплообмена составляет малую величину. Например, если горючим материалом является древесина, у которой
Tпл103 К, уд Qрн 4•105 Дж•м-2•с-1, то параметр влияния радиационного
теплообмена ΔRAD 0,06.
Решение уравнения (4.85) можно получить численным методом. Этим методом были рассчитаны средние коэффициенты теплопотерь для пожаров в различных по форме и размерам помещениях при круговом распространении пламени по слою ТГМ. Результаты расчетов представлены на рис. 4.1.
Полученные результаты при условиях, когда 6 < Ф < 24 и 0,4 < ГФ < 2, с
достаточной для практики точностью аппроксимируются формулой
(4.86)
При линейном распространении пламени по поверхности ТГМ уравнение (4.79) преобразуется с помощью формулы (4.45) в следующее уравнение:
(4.87)
Из ранее полученной формулы (4.52) при п=2 следует
(4.88)
где - среднее значение коэффициента теплопотерь.
Рис 4 1 Зависимость среднего коэффициента теплопотерь
от параметра ГФ для пожаров при круговом
распространении пламени по слою древесины
Уравнение (4.87) преобразуется с помощью формулы (4.88) и соотношения (4.83) в формулу
Ф.
(4.89)
Подставим полученное выражение для Ф. в формулу (4.80). После интегрирования получается уравнение для расчета среднего коэффициента теплопотерь в случае линейного распространения пламени
(4.90)
где ; ΔRAD - то же, что и в формуле (4.85).
Уравнение (4.90) можно преобразовать в следующее:
(4.91)
где
Уравнение для расчета среднего коэффициента теплопотерь при условиях, когда в помещении горит ГЖ, выводится так же, как уравнения (4.85) и (4.91).
Следует сказать, что эффективная поверхность радиации FRAD, которая содержится в формуле (4.74), при неустановившемся режиме выгорания ГЖ во времени изменяется. Чтобы учесть это, можно принять,
что
Формула для коэффициента теплопотерь ф при нестационарном горении жидкости с учетом ранее сделанной оценки (4.78) и зависимости (4.48) принимает следующий вид:
(4.92)
Преобразуем уравнение (4.92) с помощью зависимости (4.52). В результате получим выражение
(4.93)
После несложных преобразований уравнение (4.93) принимает следующий вид:
(4.94)
Значение среднего коэффициента теплопотерь для интервала времени, равного критической продолжительности пожара, вычисляется путем подстановки уравнения (4.94) в уравнение (4.80). После интегрирования получается уравнение
(4.95)
где (4.96)
Введем следующие обозначения:
;
С учетом введенных обозначений уравнение (4.95) принимает вид
(4.97)
В заключение анализа, посвященного среднему коэффициенту теплопотерь, рассмотрим гипотетический процесс выгорания ГЖ, когда ст → 0. Формулы для расчета критической продолжительности пожара и среднего коэффициента теплопотерь существенно упрощаются, если сделать допущение, что процесс стабилизации горения жидкости происходит мгновенно. В этом случае уравнение, описывающее изменение средней температуры во времени, имеет следующий вид:
(4.98)
Уравнение для определения коэффициента теплопотерь * в этом случае имеет вид
(4.99)
Формула для вычисления среднего коэффициента теплопотерь получается из уравнения (4.80) после подстановки в него выражения (4.99):
(4.100)
Формула для расчета критической продолжительности по температуре в этом случае имеет вид
(4.101)
где - средний коэффициент теплопотерь, который вычисляется по формуле (4.100).
Эта формула совпадает с рекомендуемой государственным стандартом формулой для расчета критической продолжительности пожара при горении ГЖ. Дополнением к рекомендации ГОСТ является полученная здесь формула (4.100) для определения коэффициента теплопотерь.
Из анализа следует, что формулу, рекомендуемую государственным стандартом, можно принять, строго говоря, лишь при условиях, когда
Ткр >:> ст (т.е. когда продолжительность начальной стадии пожара много больше времени стабилизации горения ГЖ).
- Глава 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