Все вышесказанное обуславливает необходимость разработки и практического применения научного подхода к проблеме предупреждения и ликвидации лесных пожаров. Особое внимание необходимо уделять опасности перехода лесных пожаров на близлежащие промышленные объекты.
Помимо инструментальных средств контроля промышленного объекта и сопредельной территории важное значение имеет разработка новых подходов к прогнозированию пожарной безопасности такого объекта на базе технологий математического моделирования и геоинформационного мониторинга [2].
Цель настоящей работы – разработка детерминированной математической модели для оценки параметров воздействия лесного пожара на промышленный объект.
Физико-математическая постановка
Основное допущение в настоящей работе – промышленный объект подвергается воздействию низового лесного пожара (интенсивность пожара может быть сильная, средняя и слабая в зависимости от предполагаемого сценария). Передняя кромка лесного пожара распространяется с малой скорости в условиях отсутствия ветра. Форма кромки лесного пожара в окрестности промышленного объекта может быть описана прямой в силу ее большой протяженности. Передача тепла от кромки пожара осуществляется за счет излучения. Рассматривается сценарий, когда ограждающие конструкции промышленного объекта выполнены из древесины (для определенности из сосновой древесины). В слое ограждающей конструкции преобладает кондуктивный теплоперенос. Влага в приповерхностных слоях ограждающих конструкций отсутствует – рассматривается сценарий катастрофических лесопожарных погодных условий. основным параметров безопасности промышленного объекта в указанных условиях является установление факта воспламенения (или не воспламенения) ограждающих конструкций промышленного объекта. В качестве критериев воспламенения может быть использован теплофизический подход, описанный в работе [4]. Использование такого подхода апробировано при решении задач прогностического моделирования лесной пожарной опасности [5].
Геометрия области решения с учетом предложенных допущений представлена на рисунке 1. В таблице 1 приведены данные [4], которые могут быть использованы в качестве критериев воспламенения ограждающих конструкций промышленного объекта.
Рис. 1. Геометрия области решения: 1 - кромка лесного пожара, 2 - ограждающие конструкции объекта, 3 – внутреннее пространство объекта, qff – тепловой поток от кромки лесного пожара
Таблица 1
Экспериментальные данные по условиям зажигания древесины сосны [4]
Время задержки зажигания, с |
Тепловой поток, кВт/м2 |
Температура поверхности, К |
63.5 |
12.5 |
658 |
45.0 |
21 |
700 |
11.1 |
42 |
726 |
2.6 |
84 |
773 |
0.4 |
210 |
867 |
Сформулированная математическая модель (1) – (2) с краевыми и начальными условиями (4) – (8) решена методом конечных разностей. Для решения разностных аналогов одномерных дифференциальных уравнений использовался метод прогонки.
Уравнение теплопроводности для ограждающей конструкции:
(1)
Уравнение теплопроводности для внутреннего пространства:
(2)
Граничные условия:
На границе воздействия лесного пожара:
(3)
На остальных внешних границах объекта:
(4)
(5)
(6)
На границах двух сред:
(7)
(8)
Начальные условия:
(9)
Где Ti - температура (1 – ограждающая конструкция, 2 – внутренне пространство объекта), Te - температура окружающей среды; ri, ci, li - плотность, теплоемкость и теплопроводность (1 – ограждающая конструкция, 2 – внутренне пространство объекта); t – время; x, y – пространственные координаты; qff – тепловой поток от кромки лесного пожара.
Результаты и обсуждение
В результате математического моделирования получены поля температур в структуре промышленного объекта. Было исследовано несколько базовых сценариев по влиянию излучения от лесного пожара на промышленную постройку. Во-первых, варьировалась величина плотности теплового потока от фронта низового пожара. Исходными данными послужили результаты экспериментальных выжиганий на специальных площадках Красноярского края [3]. Во-вторых, варьировалось время непрерывного воздействия излучения лесного пожара на промышленную постройку. На рисунках 2 и 3 представлены типичные поля температур в структуре промышленного объекта. Анализ полученных результатов позволяет сделать вывод, что при низких значениях теплового потока от фронта лесного пожара поверхность ограждающих конструкций не нагревается до критических температур при подходе фронта пожара на расстояние 2-3 метра от промышленного объекта. Следует отметить, что есть корреляция между высотой пламени и тепловым потоком из зоны горения [3]. Чем выше факел пламени, тем больше интенсивность тепловыделения из зоны горения [3]. Для сценария низового лесного пожара высокой интенсивности достижение критических температур [4] наблюдается для вариантов и кратковременного и более длительного воздействия.
a)
б)
Рис. 2. Распределение температуры в промышленном объекте или жилой постройке в различные моменты времени (тепловой поток от фронта пожара 4300 Вт/м2):
a - t=30 с; b – t=120 с
Анализ распределений температуры по глубине ограждающей конструкции показывает, что до критических температур нагревается поверхностный слой конструкционного материала толщиной примерно 1 – 1,5 сантиметра. Результаты экспериментов с деревянными конструкционными материалами показывает [7], что такой глубины прогрева достаточно для образования необходимого объема газообразных продуктов пиролиза. Эти газы воспламеняются и приводят к возгоранию промышленной постройки. Кроме того, газообразные горючие продукты пиролиза могут воспламениться в результате попаданию нагретой частицы из фронта лесного пожара в зону их смешивания с окислителем.
Были проведены вычислительные эксперименты по оценке влияния скорости распространения фронта лесного пожара. Рассматривался сценарий низового лесного пожара высокой интенсивности, который приближается к промышленному объекту или жилой постройке с постоянной скоростью.
Рис. 3.Распределение температуры в строении при скорости фронта лесного пожара u=0,5 м/мин
Рис. 4. Распределение температуры в строении при скорости фронта лесного пожара u=2 м/мин
Рис. 5. Распределение температуры в строении при скорости фронта лесного пожара u=4 м/мин
Анализ результатов численных расчетов показал слабую зависимость условий зажигания от скорости распространения низового лесного пожара в указанном диапазоне значений (u=0,5 – 4 м/мин). Это обстоятельство позволяет использовать приближение неподвижной кромки пожара при сценарных расчетах условий зажигания деревянного строения.
Выводы
В процессе настоящей работы разработана детерминированная математическая модель для оценки параметров воспламенения конструкций промышленного объекта или жилой постройки при воздействии лесного пожара. Математическая модель разработана на основе предложенных допущений. Использование такого математического аппарата может внести значительный вклад в развитие нового поколения геоинформационных систем, включающих программные компоненты для оценки параметров воспламенения построек на территории промышленного объекта или населенного пункта.
Рецензенты:Барахнин В.Б., д.т.н., с.н.с., Институт вычислительных технологий СО РАН, г. Новосибирск;
Немова Т,Н., д.т.н., профессор, Томский государственный архитектурно-строительный университет , г. Томск.
Библиографическая ссылка
Барановский Н.В., Сагалаков А.В. УПРОЩЕННАЯ ДЕТЕРМИНИРОВАННАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЗАЖИГАНИЯ ДЕРЕВЯННОГО СТРОЕНИЯ ОТ ФРОНТА ЛЕСНОГО ПОЖАРА // Современные проблемы науки и образования. – 2015. – № 1-1. ;URL: https://science-education.ru/ru/article/view?id=18584 (дата обращения: 11.10.2024).