- •Прогнозирование опасных факторов пожара Курс лекций
- •С. А. Сазонова
- •Введение
- •Исходные понятия и общие сведения о методах прогнозирования опасных факторов пожара в помещениях
- •§1. Основные понятия и уравнения интегральной математической модели пожара в помещении
- •1.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. Основные положения зонного моделирования пожара, численная реализация зонной математической модели
- •§6. Основы дифференциального метода прогнозирования опасных факторов пожара, численная реализация дифференциальной математической модели
- •Заключение
- •Библиографический список
- •Приложение № 6 к пункту12 Методики [6].
- •I. Порядок проведения расчета
- •II. Классификация и область применения методов математического моделирования пожара
- •III. Интегральная математическая модель расчета газообмена в здании, при пожаре
- •IV. Математическая двухзонная модель пожара в здании
- •Оглавление
- •Прогнозирование опасных факторов пожара
- •3 94006 Воронеж, ул. 20-летия Октября,84
2.4. Распределение перепадов давлений по высоте помещения
Разность между наружным и внутренним давлениями на разной высоте может быть различна. Другими словами, перепад давлений изменяется с высотой [4].
Рассмотрим сначала область, лежащую выше ПРД. На высоте, соответствующей координате у (см. рис. 2.2), разность давлений определим с помощью уравнений (2.3) и (2.8). Из этих уравнений следует
Δp= pвн - pнар =(pm -pa) + g(ρm –ρa) (h-y), (2.11)
где Δp - перепад давлений в области у > у*.
Изменение этой разности давлений при увеличении высоты на малую величину dy получим, дифференцируя уравнение (2.11)
d(Δp)= - g(ρm –ρa) dy, (2.12)
Чтобы получить распределение перепадов давлений по высоте помещения над ПРД. проинтегрируем правую чисть уравнения (2.12) по у в пределах от у = у* до текущего значения у, а левую - в пределах от Δр = 0 (так как при у = у* перепад Δр = 0) до текущего значения Δр.
Δp= = g(ρa –ρm) (y-y*), (2.13)
В области, лежащей ниже ПРД, распределение перепадов давлений определяется аналогичным образом. Перепад давлений в этой области определяется как разность наружного и внутреннего давлений, т.е.
δp= pнар - pвн, (2.14)
где δр - перепад давлений в области у<у*.
После математических операций, аналогичных использованным выше, получается следующая формула, описывающая распределение перепадов давлений по высоте помещения под ПРД (т.е. в области у < у*):
δp= - g(ρa – ρm) (y*-y). (2.15)
Знак "минус" перед правой частью в уравнении (2.15) означает, что перепад давлений δp действует в направлении, противоположном перепаду давлений Δp.
Абсолютное значение перепада давлений ниже ПРД определяется по формуле
|δp|= g(ρa – ρm) (y*-y). (2.15a)
2.5. Формулы для расчета расхода газа, выбрасываемого через прямоугольный проем
Если проем целиком расположен выше ПРД, то через все участки этого проема происходит истечение газа из помещения. Скорость газа зависит от перепада давлений (уравнение Бернулли) [4] :
. (2.16)
Подстановка в это уравнение вместо Δр выражения (2.13) дает формулу, позволяющую рассчитать распределение скоростей по высоте проема:
(2 16а)
Разобьем проем по высоте на малые участки, площадь которых составляет величину b∙dy, где b - ширина проема. Расход через такой малый участок проема равен произведению площади этого участка проема на скорость истечения и на плотность истекающего газа, т.е.
. (2.17)
Расход газа через весь проем получим, интегрируя правую часть уравнения (2.17) по высоте проема, т.е.
, (2.18)
где уН - координата нижнего края проема, м; уВ - координата верхнего края проема, м.
После подстановки в уравнение (2.18) выражения (2.16а) получим следующий результат:
. (2.19)
В результате после интегрирования получается следующая формула для расчета расхода газа через проем при условии, когда уН ≥ у*:
. (2.20)
В том случае, когда плоскость равных давлений пересекает проем и распределяет его на две части, количество газа, выбрасываемого из проема, вычисляется также путем интегрирования выражения (2.17). Однако в этом случае интегрирование правой части уравнения (2.17) производится в пределах от у* до уа, т.е.
, (2.21)
Формула для расчета расхода уходящих газов через проем при смешанном режиме его работы, которая получается после интегрирования, имеет следующий вид:
. (2.22)