Моделирование мезомасштабных процеcсов переноса примеси при наличии нестационарного источника тепла (лесного пожара)

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Моделирование мезомасштабных процеcсов переноса примеси при наличии нестационарного источника тепла (лесного пожара)

А. А. Вязников

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

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

Введение© Вязников А. А., Шварц К. Г., 2013

Работа выполнена в рамках тематического плана НИР ПГНИУ №5.5061.2011 "Математическое моделирование мезомасштабных атмосферных процессов переноса вредных примесей от нестационарных тепловых источников".

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

В работе Гришина [5] дается краткий обзор математических моделей природных пожаров. Одной из известных лесопирологических моделей является модель Доррера [4], которая позволяет достаточно просто определить контур лесного пожара. В то же время этот подход не позволяет учитывать взаимное влияние приземного слоя атмосферы и фронта лесного пожара друг на друга [5], так как скорость распространения фронта считается эмпирически заданной.

Уравнения модели

Для решения задачи использована квазидвумерная модель распространения многокомпонентной примеси от источника с учетом термической неоднородности подстилающей поверхности [1-3]:

(1)

, (2)

(3)

(4)

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

Для учета нестационарного источника пожара вводится уравнение фронта

(5)

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

Задача решается со следующими начальными условиями:

; (6)

граничные условия:

Численные расчеты

Расчеты проводились на основе системы уравнений (1)-(7).

Рассматривается лесная площадка размера L L при L = 40 км, фоновая температура в лесу 15оС, толщина слоя D = 1000 м.

Для расчетов использована явная конечно-разностная схема [1] на сетке 200 200 узлов. На площадке в узле с номером (70, 155) находится стационарный источник торфяного пожара. Температура источника 450оС. Ветер северо-западный - азимут 330о. Скорость ветра на границе слоя 2 м/c.

Через некоторое время начинается низовой лесной пожар. Лесной пожар перемещается с течением времени и является нестационарным источником примеси и тепла. Температура пожара также бралась равной 450оС.

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

Параллельный алгоритм расчета реализован на языке Java, расчеты проводились на базе процессора Intel Core i3.

Численный эксперимент включал два этапа: расчеты со стационарным источником и расчеты с нестационарным источником.

Рис. 1. Поле концентрации примеси при t=1 (стационарный источник)

Рис. 2. Поле средней в слое температуры воздуха при t=1 (стационарный источник)

На рис. 1 и 2 представлены результаты расчета со стационарным источником: поле концентрации примеси и средней в слое температуры имеют форму факела, направленного вдоль вектора скорости ветра.

Рис. 3 а. Фронт пожара при t=1.5 (нестационарный источник)

Рис. 3 б. Поле концентрации примеси при t=1.5 (нестационарный источник)

Рис. 3 в. Фронт пожара при t=3 (нестационарный источник)

Рис. 3 г. Поле концентрации примеси при t=3 (нестационарный источник)

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

Над фронтом пожара возникает слабый вихрь (рис. 4, а), который не влияет на искривление траектории примеси и поле средней температуры воздуха в слое (рис. 4, б).

В расчетах единица модельного времени соответствует 20000с, что примерно составляет 5,5 часов. Рисункам соответствует безразмерное значение времени t.

Скорость пожара равна 0,2 м/c, максимальная скорость ветра меняется менее чем на 0,1 процента. Максимальная концентрация примеси 1,91 в долях ПДК, максимальная угловая скорость вихря составляет в безразмерном виде.

Рис. 4 а. Поле среднего вихря в слое при t=1 (нестационарный источник)

Рис. 4 б. Поле средней в слое температуры воздуха при t=1 (нестационарный источник)

Выводы

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

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

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

температура пожар лесной двумерный

Список литературы

Шварц К.Г., Шкляев В.А. Численное моделирование атмосферных мезомасштабных процессов переноса многокомпонентной примеси при торфяном пожаре // Вычислительная механика сплошных сред. 2012. Т. 5, № 3. С. 274-283.

Шатров А.В., Шварц К.Г. Численное моделирование атмосферных мезомасштабных процессов переноса примесей в окрестности города Кирова // Вычислительная механика сплошных сред. 2010. Т.3, №3. С. 117-125.

Аристов С.Н., Шварц К.Г. Вихревые течения в тонких слоях жидкости. Киров: ВятГУ, 2011. 207 с.

Доррер Г.А. Динамика лесных пожаров. Новосибирск: Изд-во СО РАН, 2008. 404 с.

Гришин А.М. О математическом моделировании природных пожаров и катастроф // Вестник Томского государственного университета. 2008. №2. С. 105-114.

Размещено на Allbest.ru