[Yox,0 + (Yf, f - Yox,0)f] - Yox,0
, 0 < f < fst
Ypr = [Yox,0 + (Yf, f - Yox,0)f] - Yf, f 
, fst < f < 1
где стехиометрическое fst значение находим:
где Yox,0 - массовая доля кислорода в потоке окислителя, a Yf, f - массовая доля топлива в потоке газообразных продуктов пиролиза.
Очевидно, что при этом не учитывается влияние турбулентных пульсаций на химическую реакцию. Они могут быть учтены с помощью диффузионно-вихревой модели [13]. В этой модели, кроме уравнения переноса для f решается уравнение для Yf.
В ней в случае открытого пожара скорость реакции будет определяться местной концентрацией горючего, за исключением области вблизи источника продуктов пиролиза. При регулируемых вентиляцией пожарах в помещениях наблюдается дефицит воздуха, и, следовательно, потребление топлива будет определяться концентрацией кислорода. Третий член вводится для ограничения скорости реакции в холодных смесях
где С = 4, а В полагают равным 2.
Предположение для замыкания источникового члена (формула (3.31)) позволяет, помимо уравнения переноса для f, решать уравнение для массовой доли топлива и рассчитывать массовую долю каждого компонента упрощенной химической реакции. Модели этого типа успешно использовались при решении различных задач пожарной безопасности и оптимизации процесса горения в промышленных установках. Достоинством модели является ее простота. Она позволяет рассчитывать распределенное по объему выделение энергии, определяемое геометрией помещения и доступом воздуха. Можно определить концентрации CO2 и Н2O, если предположить, что они являются единственными продуктами горения.
Однако с помощью такой схемы нельзя учесть влияние конечности скорости химических реакций. Для корректного расчета концентраций продуктов неполного окисления, таких, как СО и сажа, необходима более усложненная модель.
Довольно перспективной является модель ламинарных элементов пламени [14, 16]. В ней предполагается, что горение происходит только в тонких ламинарных элементах пламени, входящих в турбулентное поле потока. Соотношения между мгновенным химическим составом и функцией смешения в таких условиях могут быть определены вычислительным путем, для простых горючих, таких, как метан и пропан, с достаточно хорошо известной кинетикой химических реакций. Однако встречающаяся на практике горючая нагрузка обычно имеет сложный химический состав, поэтому, из-за отсутствия соответствующих соотношений, в настоящее время данная модель мало применима для практических задач.
Радиационный теплоперенос
Наиболее простым способом учета радиационных тепловых потерь является так называемая cR -модель. Она состоит в том, что мощность тепловыделения в очаге горения путем занижения теплоты сгорания уменьшается на долю тепла cR, теряемую за счет излучения. Эта доля задается на основе экспериментальных данных в зависимости от вида топлива. Несмотря на кажущуюся примитивность, такая модель на начальной стадии пожара часто дает хорошие результаты.
Однако часто возникают задачи, требующие более точного моделирования радиационного теплопереноса.
Влияние радиационного теплопереноса выражается через источниковый член в уравнении сохранения энергии. Кроме того, радиационные потоки сильно влияют на температуры поверхностей стен помещения, а следовательно, на распространение пламени.
Основное уравнение радиационного переноса можно записать в виде
где I - интенсивность радиационного излучения в направлении W; s - расстояние в направлении W; Eg = s
- энергия, излучаемая абсолютно черным газом при температуре газа Tg; ka и ks - коэффициенты поглощения и рассеяния; Р(W, W') - вероятность того, что излучение в направлении W' после рассеяния попадет в телесный угол dW в окрестности направления W. Это уравнение необходимо интегрировать по всем направлениям и длинам волн. Для большинства практических задач точное решение невозможно, вместо него разработано несколько приближенных методов, которые и используются для моделирования динамики пожаров в помещениях.
Потоковые методы
Если разделить пространственное и угловое распределение интенсивности излучения, задачу можно существенно упростить. Этот подход используется в "потоковых методах" [15]. Если предположить, что спектральная интенсивность постоянна в пределах заданных интервалов телесного угла, то уравнение радиационного переноса сводится к нескольким связанным между собой обыкновенным линейным дифференциальным уравнениям относительно осредненных по пространству интенсивностей или потоков излучения.
Если телесные углы совпадают с поверхностями контрольного объема в декартовом пространстве и если предположить, что поток излучения через каждую поверхность однороден, то, обозначив через Fi+ тепловой поток, проходящий через контрольный объем в положительном направлении i, и через Fi- - поток в отрицательном направлении i, имеем
где ka и ks - локальные коэффициенты поглощения и рассеяния, а Еb - количество тепла, излучаемого контрольным объемом, если он является абсолютно черным.
Объединяя эти уравнения и дифференцируя их по xi получаем
где 
Уравнение имеет тот же вид, что и обобщенное уравнение сохранения (3.26), и может быть решено с помощью того же численного алгоритма. Вклад излучения в источниковый член уравнения энергии для каждого контрольного объема
Эта модель очень привлекательна для использования в полевых моделях, поскольку в ней используется тот же численный метод, что и для решения уравнений гидродинамики. Однако этот метод имеет ряд недостатков, среди которых одним из главных, применительно к пожарам, является неточность метода при моделировании радиационного переноса под углом к декартовой сетке.
Потоковые методы годятся, например, при определении радиационного переноса от припотолочного слоя к полу помещения, но они неточны вблизи очага, где скорость распространения фронта пламени может зависеть от переноса тепла, направленного под углом к сетке.
Метод дискретного радиационного переноса
Эта модель, разработанная Локвудом и Шахом [17], преодолевает основной недостаток потоковых методов. Для нее характерны некоторые черты методов Монте-Карло, а именно прохождение "лучей" электромагнитного излучения через вычислительную область между границами. Однако в отличие от методов Монте-Карло, где направления лучей генерируются случайным образом, в этой модели они выбираются предварительно, таким же образом, как выбирается расположение гидродинамической сетки. Метод включает в себя решение уравнения радиационного переноса вдоль путей этих лучей, выбираемых обычно таким образом, чтобы они приходили в центры граничных поверхностей гидродинамических контрольных объемов.
Скачать реферат