Общая физико-математическая модель зажигания и горения древесины

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2010 Математика и механика № 2

УДК 521.534:536.245.022

А.М. Гришин ОБЩАЯ ФИЗИКО-МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЗАЖИГАНИЯ И ГОРЕНИЯ ДРЕВЕСИНЫ

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

Ключевые слова: математическое моделирование, древесина, пористая среда, фильтрация, теплопроводность, нагревание, сушка, горение.

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

В работах [1 - 24] изложены основы механики реагирующих сред, теории природных пожаров и катастроф.

Согласно [1 - 6], древесина представляет собой пористую реакционноспособную среду, физико-химические свойства которой изменяются при нагревании и горении. Поэтому для математического моделирования процессов зажигания и горения древесины целесообразно использовать модели реагирующих сред, представленные в работах [7 - 19]. Теплофизические свойства древесины можно частично определить на основе справочных данных [3 - 6], а также путем решения обратных задач механики реагирующих сред [18] с использованием результатов экспериментов по зажиганию и горению древесины различных пород деревьев.

Надо сказать, что предлагаемый комплексный метод исследования базируется не только на теоретических (законы сохранения массы, количества движения и энергии), но и на экспериментальных данных, которые характеризуют свойства древесины как пористой реакционноспособной среды [1 - 10]. Поэтому целесообразно использовать модель пористой реагирующей среды, развитую в [7 - 18], для которой должна быть создана новая база данных для описания процессов нагрева, пиролиза и горения газообразных и конденсированных процессов пиролиза древесины с использованием методов решения обратных задач механики реагирующих сред [20]. На основе этой модели с использованием итерационноинтерполяционного метода [21] и результатов, представленных в [22 - 24], могут быть решены конкретные задачи о возникновении и распространении пожаров в селах, поселках и городах. В данной работе на основе результатов [10, 11 - 17, 22

- 24] представлена общая физико-математическая модель зажигания деревянной стены в результате действия лесных или степных пожаров.

Физическая модель зажигания и горения стенки из древесины

Анализ наблюдений за горением древесины позволяет выделить следующие физико-химические процессы (этапы):

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

2) сушка древесины;

3) пиролиз древесины и образование конденсированного (кокса), газообразных (СО, Н2, СН4) и дисперсного (частиц сажи) продуктов пиролиза;

4) горение газообразных горючих продуктов пиролиза;

5) горение дисперсных продуктов пиролиза (частиц сажи);

6) горение конденсированного продукта пиролиза (коксика).

При горении древесины имеют место следующие физико-химические процессы:

1) кинетический режим горения, при котором скорость процесса определяется скоростью протекания гетерогенных химических реакций окисления конденсированных продуктов пиролиза древесины (угля);

2) диффузионный режим горения, при котором скорость разрушения древесины определяется скоростью диффузии окислителя (кислорода);

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

На рис. 1 представлена схема изменения состояния и химического состава древесины при ее нагревании и последующем горении.

Рис. 1. Схема изменения агрегатного состояния и химического состава древесины при ее горении. Стрелки обозначают: 1 - подвод теплоты в результате излучения, конвекции и теплопроводности; 2 - сушка древесины; 3 - пиролиз древесины; 4 - образование коксика; 5 - окисление газообразных и дисперсных продуктов пиролиза; 6 - образование пепла; 7 -образование и исчезновение частиц сажи; 8 - образование частиц дыма и их выветривание

Основная система уравнений

При выводе системы уравнений, выражающих законы сохранения массы, количества движения и энергии будем считать, что:

1) конденсированная фаза древесины в процессе горения состоит из четырех компонентов: исходного сухого материала, воды, конденсированного продукта пиролиза (угля) и продукта горения угля - золы;

2) при горении древесины в результате нагревания имеет место необратимая гомогенная реакция разложения

Здесь А1 - символ исходного конденсированного вещества - влажной древесины; А2 - символ воды; А3 - символ конденсированного продукта реакции (угля); А4а -символы газообразных горючих продуктов реакции разложения (пиролиза) - СО, Н2, СН4, v1, v'2, v'3 и v'4а - соответствующие стехиометрические коэффициенты.

При выводе стехиометрического соотношения (1), выражающего закон сохранения массы, использовалась схема, представленная на рис. 1.

Процессы зажигания и горения древесины целесообразно исследовать, используя законы сохранения массы, количества движения и энергии. При этом на основе данных [1 - 6] были сделаны следующие допущения:

1. Древесина является двухтемпературной горючей недеформируемой пористо-дисперсной средой, состоящей из следующих фаз: ф1 - сухое горючее органическое вещество1; ф2 - вода, связанная с этим веществом в порах; ф3 - коксик (твердого продукта пиролиза древесины); ф4 - зола; ф5 - газовая фаза; ф6 - частицы сажи; ф7 - частицы дыма и ф8 - капельки воды на границе раздела «древесина - внешняя среда» в пасмурный день.

2. Скорость и температура дисперсных частиц совпадают со скоростью и температурой несущей газовой фазы в соответствующей точке пространстве.

3. Газодисперсная фаза состоит из тех же компонентов, что и в общей физикоматематической модели торфяных пожаров.

4. Силы инерции в уравнении движения газовой фазы малы по сравнению с силой давления, и уравнения сохранения количества движения можно использовать в форме квадратичного закона Дарси [4].

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

6. Температура газовой фазы совпадает с температурой дисперсных частиц дыма.

С учетом сделанных допущений имеем следующую систему уравнений, которая получена на основе сведений, представленных в работах [6 - 19]:

У1Л ^ у2А2 + у3А3 +ХУ4аА4а •

1 Сухая древесина состоит из лигнина и целлюлозы [4 - 6].

с1Т Л д

7, / , р? ЪСп = ~

ш ~~5 дх,

Чфф і

д Т д х

і

дТ_ д х,

р Дэфф а Ср а

д Са

д х

і У

+ак (Т, - Т) + Ср, (Т, - Т)(1 -аС )О1, +

+ Ср 7(Т, Т)О2, + 4-3, О3 ) + ?52°52 + ?53°53 + ?54°54 + ?2, (О8-) °8+));

^аса о + д

= О а +-

Ж

д х,

д са

эфф а

д х

Л дТ д

Ёр ф'ср? ~ёт=*;

- Са Є, а = 1,2,..., N ;

дТ

+?2,°2, + Ъя О3я + ?4,°4, + аК (Т Т, );

Р _д^1 = -О Р ^ = -О • Р дФ3 = а О -МсО -а О -О

р1 ^ , р2 ~ = О25; р3 ^ = аСО1, , ^ О3, а4О3, О6

дґ

дґ

дФ3

дґ

Мс

М1

р4

дф4 = а4о3, - о7я) = Я4,; дґ 3я

Ч,>

Ё са=1 Ёф? =1;

Р = ргОТЁ , Nг = N - 3;

а=1 М а

е = (1 -ас щ, + О2я + М о + о6я) + о7я) + о8я).

Здесь / - время, г - радиус-вектор любой точки; х}- (г = 1,2,3) - декартовы координаты; — символ полной производной по времени; 0> - массовая скорость обра-

&

зования газодисперсной фазы; к* - коэффициент фильтрации; срг, ср1 , ср; - теплоемкости при постоянном давлении отдельных компонентов газовой фазы, водяного пара и газообразных продуктов пиролиза; V) , |у| - компоненты и модуль ос-редненной скорости газа и дисперсных частиц; р,- - истинная плотность ,-й фазы;

= ^8+ - ^8- ; ^8(-), ^8(+) - массовые скорости конденсации паров и испарения свободной воды в газодисперсной фазе; д1;, д*, д3* и д4; - теплоты пиролиза древесины, испарения связанной воды и горения коксика; д52, д53, д54 - теплоты горения СО, Н2, СН4 ; р - давление газа в порах древесины; са - массовая концентрация а-компонента в газодисперсной среде; N - количество компонентов в газодисперсной среде; N = N - 3 - количество компонентов газовой фазы; N - общее число компонентов; Яа - массовая скорость образования а-компонента газовой фазы в результате пиролиза ЛГМ, испарения воды, гетерогенных и гомогенных химических реакций; ;), ^1*), ^8*) - массовые скорости образования час-

7

тиц сажи, дыма и капелек воды при конденсации ее паров; рг = Ё раг - плот-

а=1

ность газовой фазы; ра г - парциальные плотности компонентов газовой фазы

(индекс 1 соответствует кислороду, 2 - СО; 3 - Н2; 4 - СН4; 5 - СО2; 6 - М2; 7 -парам воды); саг =раг / рг, а = 1,2,...,7 - массовые концентрации для смеси газов;

са = рОГ /р ,а = 1,2,...,10 - массовые концентрации компонентов газодисперсной смеси; рТ - парциальные плотности отдельных компонентов газодисперсной фазы (а = 1,2,_10; первые семь значений рТ соответствуют компонентам газовой

фазы, а последние 8 - частицам сажи; 9 - частицам дыма; 10 - капелькам воды); р = р5 ф5 + р6 ф6 + р7 ф7 + р8 ф8 - плотность газодисперсной смеси; р5 -р8 - истинные плотности газовой фазы, частиц сажи, частиц дыма и капелек воды в мик-ропорах древесины; ф5 -ф8 - объемные доли газовой фазы, частиц сажи, дыма и воды в порах древесины; X - коэффициент молекулярной теплопроводности газовой фазы; ц - коэффициент динамической вязкости газа в порах; р, - эмпирический коэффициент в законе фильтрации (3); к, - коэффициент проницаемости в уравнении (3); Дэфф а - эффективный коэффициент диффузии а-компонента, в

том числе сажи, золы и капелек воды; ау - коэффициент внутреннего теплообмена; К5а - массовые скорости изменения компонентов газовой фазы (индекс 1 соответствует скорости изменения О2, 2 - горению СО, 3 - Н2, , 4 - СН4); К,- массовые скорости образования (исчезновения ) вещества конденсированных фаз (индекс 1 соответствует скорости пиролиза древесины, 2 - испарению воды, связанной с древесиной, 3 - горению коксика); Т и Т, - температуры газовой и конденсированной фаз; X, - коэффициент теплопроводности пористой конденсированной фазы (древесины); дк - тепловые эффекты химических реакций; к =1, 2, 3, где 1 соответствует реакции пиролиза, 2 - массовой скорости испарения связанной воды, 3 - скорости гетерогенного коксика; ас - коксовое число древесины; Я - универсальная газовая постоянная; Мс и М\ - атомарная и молекулярная массы углерода и кислорода; Ма - молекулярная масса а-компонента газовой фазы; как правило, Е = g - ускорение силы тяготения; Хэфф и цэфф - эффективные коэффициенты теплопроводности и вязкости газа; ау = єа, - коэффициент объемного теплообмена в слое древесины; , - удельная поверхность древесины с учетом микропор [4]; а, - коэффициент теплообмена типичного элемента древесины с окружающей средой; в (10) нижний индекс ««,» соответствует реакциям с участием конденсированных веществ; верхний индекс «,» - характеристикам дисперсной фазы.

Уравнение

- закон сохранения энергии в газовой фазе; (5) - уравнение сохранения массы а-компонента газовой фазы; (6) - закон сохранения энергии в конденсированной фазе; уравнения (7), (8) - описывают поведение объемных долей различных фаз, а выражение (9) представляет собой уравнение газовой фазы в порах.

Коэффициенты переноса, теплофизические и термокинетические постоянные

Итоговые гомогенные химические реакции в проницаемом слое древесины имеют вид [10 - 12]

1) СО + 2о2= С 02, 2) С Н4 + 2 02 = С 02 + 2Н2 О.

Уравнения химической кинетики для этих реакций запишутся как [12]

&.-^«Г.-^ехр^А.)- г„

ш ЯТХ

Ту2 = —к2 ^х-0’5 Т2 ехр

сИ 7 Яц

где х1 и х2 - предэкспоненциальные множители; у1 и у2 - молярно-объемная концентрация окиси углерода СО и метана; х1, х3 х5 - молярные концентрации СО, О2, СН4; к1, к2 - предэкспоненциальные множители; Е1 и Е2 - энергии активации соответствующих реакций; г1 и г2 - молярно-объемные скорости реакции.

Для определения скорости испарения связанной воды Я2, в многофазной среде - древесине - используется аналог закона Герца - Кнудсена [12]

Я = ,2М2 А, Ф2 [к2, ехр

2" (2пЯТ1М2 )0’5 .

Здесь s2 - удельная поверхность воды; М2 - молекулярная масса воды; А, - эмпирический множитель; р2 - парциальное давление паров воды; п =3,14; е2, -энергия активации процесса испарения воды; к2, - предэкспоненциальный множитель.

Для нахождения парциального давления паров воды в древесине р2 в (12) используется закон Дальтона [10 - 13]. Тогда для р2 имеем выражение

М

Р2 = РC2 — ,

М 2

где М - молярная масса смеси; р - давление в смеси газов; с2 - массовая концентрация паров воды.

Эффективные коэффициенты диффузии определяются по формуле Фристрома

- Вестенберга [12, 19]

N

X

А* = (1 - са )(X )-1

1=1 <1 1^а

7 [

,; Ра)а (В 1,а I кТО0,17

Формула для коэффициентов теплопроводности компонентов газовой фазы X1,1 = 1, 2, ..., Ы, взята из [19]

С N

X, = X0(0,115 + 0,354-Я-), Х5 =]ГХгсг , Х° = 8,3210-2

N

где с , - теплоемкости ,-го компонента при постоянном давлении; ^(2,2) - эм-

пирические постоянные [19].

4

Коэффициенты теплоемкости компонентов газовой фазы X = Х. в Ф. +Х5Ф5 =

а, + ^Т + с-/Т2 (а,,6,,с, - эмпирические постоянные) брались из [8], а значе-

ния величин ср5 и Х.5 в конденсированной фазе взяты из [10 - 12]:

N

4

С

'р5 = X СР ]С] , Х = Xх. вФ. +Х5Ф5 .

,=1

. =1

Здесь с, - массовые концентрации компонентов газовой фазы; N - число компонентов газовой фазы; ср, - теплоемкости отдельных компонентов; X - коэффициент теплопроводности многофазной среды; К и Х5 - коэффициент тепло -проводности компонентов конденсированной и газовой фаз.

Массовые скорости химических реакций Я^ Яи, Я2 в уравнениях (2), (3)-(7) имеют вид

где Я1, Яи, Я2 - размерные величины, упомянутых выше скоростей реакций, а величины

Над очагом горения древесины возникает конвективная колонка, содержащая дисперсные (частицы кокса и сажи) и газообразные продукты горения древесины). Математическое описание параметров состояния этой колонки (скорости, давления, температуры, концентраций компонентов) представляет сложную задачу математической физики. Эта система уравнений может быть получена как частный случай из общей математической модели лесных пожаров [12 - 16]. Для упрощения этой системы были сделаны следующие допущения:

1. В результате горения древесины возникает конвективная колонка из нагретого очагом пожара воздуха, а также продуктов пиролиза и горения древесины.

Я1 П1 Я1х М1 Г1 , Я2 П2 Я1х Я2х + 2 М2 Г2,

Я3 = - п3 Я3х - М3 г1 /2 - 2 М3 г2, Я4 = М4 (г1 + г2),

Я5 = П4 Я1х - М5 Г2 ;

а

О математическом моделировании процессов переноса в приземном слое атмосферы при горении древесины

2. Для упрощения постановки задачи конвективную колонку можно считать осесимметричной.

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

4. На границе раздела атмосферы и очага горения древесины должны выполняться соотношения, выражающие законы сохранения массы, количества движения и энергии, примеры которых даны в [12].

Кроме того, для математического моделирования такой нестационарной задачи, как зажигание и горение древесины, необходимы начальные условия, аналогичные условиям, представленным в [14].

База данных общей математической модели горения древесины и методики численного решения задач

В основную систему уравнений (2) - (10) входят теплофизические коэффициенты, скорости сушки и реакций пиролиза и горения и множество эмпирических постоянных, совокупность которых представляет базу данных общей математической модели. Компоненты базы данных могут быть определены, с одной стороны, из литературных источников, а с другой - на основе экспериментальных исследований с помощью решения обратных задач механики реагирующих сред [20]. Очень важно, чтобы все компоненты базы данных были согласованы друг с другом. Это осуществлялось методом их варьирования в пределах их точности и выбора из условия наилучшего согласования соответствующих расчетных и экспериментальных данных. Очень важно при этом, чтобы решение прямых и обратных задач математической теории лесных пожаров при численном решении определялось с контролируемой точностью, что представляет собой сложную математическую задачу, так как основная система уравнений (2) - (10) нелинейна. Для численного решения прямых задач теории лесных пожаров использовались разработанные в ТГУ итерационно-интерполяционный метод [18] и специальные методики численного счета, основанные на применении метода Патанкара -Сполдинга.

В систему уравнений (2) - (10), представленную выше, входит ряд параметров, которые характерны не только для торфяных пожаров, но и для горения древесины (Xж, Ба - коэффициенты теплопроводности и диффузии, теплофизические и термокинетические коэффициенты). К сожалению, из-за ограниченности полноценных экспериментальных данных по возникновению и распространению пожаров конкретной базы данных для математической модели (2) - (10) нет. Для ее создания необходимы результаты соответствующих экспериментальных исследований и использование методов решения обратных задач механики реагирующих сред [20]. После определения параметров, характеризующих теплофизические и реакционные свойства древесины, необходимо расчеты проводить с использованием правдоподобных допущений о порядках упомянутых выше величин, а в дальнейшем в качестве основного метода решения задач теории горения древесины целесообразно использовать итерационно-интерполяционный метод и некоторые приближенные аналитические методы [21].

Выводы и предложения по развитию теории зажигания и горения древесины при пожарах в деревнях, поселках и городах

По своей структуре (физическим, химическим и прочностным) свойствам древесина и торф значительно отличаются друг от друга. В иллюстрированном энциклопедическом словаре [1] дано следующее определение понятия древесины: «древесина - это ткань высших растений, служащая для проведения воды и растворов солей от корней к листьям и другим органам. Древесиной называют также срубленные и реликтовые стволы деревьев (деловая и строительная древесина). Древесину различных пород, имеющую красивый цвет, блеск и оригинальную текстуру, называют декоративной и используют для отделки мебели, интерьеров и художественных промыслов. Самая тяжелая древесина у железного дерева (1490 кг/м3), а самая легкая - у бальзового дерева (40 - 384 кг/м3)».

В этом же словаре говориться о том, что «торф - осадочная порода, продукт неполного разложения растений в условиях болот» [1]. Кроме растительных остатков торф содержит аморфное органическое вещество (гумус), минеральные примеси и воду. В естественном состоянии - однородная по составу и окраске плотная масса черного или коричневого цвета. Он используется как удобрение, теплоизоляционный материал и реже - как топливо.

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

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

1. Создание экспериментальных установок и методик для определения теплофизических характеристик для определения коэффициентов теплопроводности и диффузии в древесине в обычных условиях.

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

3. Создание детерминированно-вероятностной модели для прогноза возникновения горения древесины в конкретных метеорологических и техногенных условиях с использованием геоинформационных технологий.

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

5. Исследование предельных условий распространения горения древесины и разработка новых способов и устройств для их тушения.

Решение этих вопросов позволит в первом приближении завершить создание математической теории горения древесины и использовать ее для разработки как способов, так и средств для борьбы с горением древесины, а также прогнозов экологических последствий этих пожаров. На основании имеющегося опыта исследования задач теории лесных пожаров [9, 12 - 16] можно утверждать, что для разработки теории горения древесины необходимо сочетать экспериментальные и теоретические методы исследования. В качестве методической основы теории горения древесины целесообразно использовать механику многофазных реагирующих сплошных сред и метод математического моделирования с применением технологии распараллеливания вычислительных операций на суперкомпьютерах,

что позволит существенно уменьшить затраты машинного времени и получать поля температур и другие характеристики горения древесины в режиме, опережающем реальное время развития этого пожара. Для сокращения числа эмпирических входных данных, например эмпирических коэффициентов тепло- и массо-обмена, целесообразно использовать теорию сопряженных задач и численные методы механики реагирующих сред [7, 10, 11, 14, 20], а также методы решения некорректных (обратных) задач математической физики [21] и теории катастроф [22 - 24].

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

Кроме того, в рамках предложенной физико-математической модели можно определить поля плотности излучения над очагом горения древесины в различные моменты времени, что в принципе позволяет, используя метод решения обратных задач [20] и математическую теорию катастроф [22 - 24], создать новую методику обнаружения и диагностики состояния среды в очаге горения древесины по данным аэрокосмического мониторинга.

Для математического моделирования конкретных задач о горении древесины (например, при пожарах в городах, поселках и деревнях) можно использовать итерационно-интерполяционный метод [21]. Кроме того, необходимо задать конкретные граничные и начальные условия. Используя экспериментальные данные об упомянутых выше пожарах и математическую технологию решения обратных задач математической физики, необходимо определить коэффициенты переноса, тепломассообмена, энергии активации и предэкспоненциальные множители различных химических реакций [21]. Знание упомянутых выше характеристик позволит осуществлять достаточно точное математическое моделирование зажигания и горения древесины с использованием вычислительных технологий, развитых в [21 - 24].

Автор благодарит профессора, д.т.н. А.С. Якимова за обсуждение работы.

ЛИТЕРАТУРА

1. Бородулин В.И., Горкин А.П., Гусев А.А. и др. Новый иллюстрированный энциклопедический словарь. М.: БСЭ, 1999. 911 с.

2. Пожарная безопасность. Энциклопедия. М.: ФГУ ВНИИ ПО МЧС России, 2007. 416 с.

3. Мищенко К.П., Равдель А.А. Краткий справочник физико-химических величин. Л.: Химия, 1972. 200 с.

4. Харук Е.В. Проницаемость древесины газами и жидкостями. Новосибирск: Наука, 1976. 187 с.

5. НовосельцевВ.Д. Справочник лесничего. М.: Агропромиздат, 1986. 352 с.

6. Чудинов С.В., Трофимов А.Н. и др. Справочник лесохимика. М.: Лесная промышленность. 1987. 271 с.

7. Гришин А.М. Математическое моделирование некоторых нестационарных аэротермо-химических явлений: дис. ... докт. физ.-мат. наук. Томск: Том. гос. ун-т, 1974. 432 с.

8. Гришин А.М., Якимов А.С. Исследование переходных режимов горения пористых реагирующих веществ // Численные методы механики сплошных сред. 1976. Т. 7. № 6. С. 39 - 51.

9. Гришин А.М. Математические модели лесных пожаров. Томск: Изд-во Том. ун-та, 1981. 277 с.

10. Гришин А.М., Фомин В.М. Сопряжённые и нестационарные задачи механики реагирующих сред. Новосибирск: Наука, 1984. 319 с.

11. Алексеев Б.В., Гришин А.М. Физическая газодинамика реагирующих сред. М.: Высш. школа, 1985. 464 с.

12. Гришин А . М. Математическое моделирование лесных пожаров и новые способы борьбы с ними. Новосибирск: Наука, 1992. 408 с.

13. Бурасов Д.М., Гришин А.М. Математическое моделирование низовых лесных и степных пожаров. Кемерово: Практика, 2006. 134 с.

14. Linn R.R. A transport Model for Prediction of wild fire Behavior Los Alamos National Laboratory. 1997. 194 p.

15. Гришин А.М. О математическом моделировании торфяных пожаров // Вестник ТГУ. Математика и механика. 2008. № 3(4). С. 85 - 95.

16. Гришин А.М. Общие математические модели лесных и торфяных пожаров и их приложения // Успехи механики. 2008. Т. 2. № 4. С. 41 - 89.

17. Гришин А.М., Якимов А.С., Рейн Г., Симеони А. О математическом и физическом моделировании возникновения торфяных пожаров // ИФЖ. 2009. Т. 82. № 5. С. 1210 - 1217.

18. Гришин А.М., Якимов А.С. Математическое моделирование возникновения торфяных пожаров // Изв. вузов. Физика. 2009. Т. 52. № 2/2. С. 112 - 121.

19. Гиршфельдер Дж., Кертис Ч., Берд Р. Молекулярная теория газов и жидкостей. М.: ИИЛ, 1961. 929 с.

20. Гришин А.М., Зинченко В.И., Кузин А.Я. и др. Решение некоторых обратных задач механики реагирующих сред. Томск: Изд-во Том. ун-та, 2006. 418 с.

21. Гришин А.М., Зинченко В.И., Ефимов К.Н. и др. Итерационно-интерполяционный метод и его приложения. Томск: Изд-во Том. ун-та, 2004. 320 с.

22. Гришин А.М. Моделирование и прогноз катастроф. Ч. 1. Томск: Изд-во Том. ун-та, 2003. 526 с.

23. Гришин А.М. Моделирование и прогноз катастроф. Ч. 2. Кемерово: Практика, 2005. 560 с.

24. Гришин А.М., Петрин С.В., Петрина Л.С. Моделирование и прогноз катастроф. Ч. 3. Томск: Изд-во Том. ун-та, 2006. 576 с.

СВЕДЕНИЯ ОБ АВТОРЕ:

ГРИШИН Анатолий Михайлович - доктор физико-математических наук, профессор, заведующий кафедрой физической и вычислительной механики Томского государственного

университета. E-mail: fire@mail.tsu.ru

Статья принята в печать 29.04.2010 г.