Моделирование газификации твердого топлива в фильтрационном режиме тема диссертации и автореферата по ВАК РФ 01.04.17, кандидат физико-математических наук Салганский, Евгений Александрович

• Салганский, Евгений Александрович
• кандидат физико-математических наук
• 2004

Оглавление диссертации кандидат физико-математических наук Салганский, Евгений Александрович

ВВЕДЕНИЕ.

ГЛАВА 1. ЛИТЕРАТУРНЫЙ ОБЗОР.

ГЛАВА 2. ОПИСАНИЕ СТАЦИОНАРНЫХ ТЕМПЕРАТУРНЫХ

РАСПРЕДЕЛЕНИЙ В ПРОТИВОТОЧНОМ РЕАКТОРЕ

ФИЛЬТРАЦИОННОГО ГОРЕНИЯ.

2.1 Воздушная газификация. Постановка задачи.

2.2 Результаты расчетов и обсуждение.

2.3 Переходная волна. Постановка задачи.

2.4 Результаты расчетов и обсуждение.

2.5 Паровоздушная газификация. Постановка задачи.

2.6 Результаты расчетов и обсуждение.

ГЛАВА 3. НЕСТАЦИОНАРНАЯ МОДЕЛЬ, УЧИТЫВАЮЩАЯ

КИНЕТИКУ ПРОТЕКАЮЩИХ РЕАКЦИЙ, ЗАВИСИМОСТЬ

ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК ФАЗ ОТ ТЕМПЕРАТУРЫ И

СОСТАВА, КОНЕЧНЫЙ РАЗМЕР РЕАКТОРА.

3.1 Постановка задачи.

3.2 Результаты расчетов и обсуждение. 61 ЗАКЛЮЧЕНИЕ. 81 ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ. 84 ОБОЗНАЧЕНИЯ. 86 СПИСОК ЛИТЕРАТУРЫ.

Введение диссертации (часть автореферата) на тему «Моделирование газификации твердого топлива в фильтрационном режиме»

Под фильтрационным горением (ФГ) понимается распространение волн экзотермического превращения в пористой среде при фильтрации газа [1-4]. Специфическим элементом, определяющим особенность горения этого класса систем, является фильтрация газа, выступающего не только в роли участника химической реакции, но и теплоносителя, формирующего структуру волны горения. Другим важным фактором является возможность сильного разбавления реагентов инертной примесью, которая не накапливается в зоне реакции, а уносится с потоком.

При фильтрационном горении может наблюдаться явление «сверхадиабатических разогревов» (температура в зоне горения становится выше адиабатической), обусловленное концентрацией выделяющегося тепла в зоне горения за счет внутреннего теплообмена между газовой и конденсированной фазами [5-7]. Именно благодаря этому явлению становятся практически осуществимы различные технологические процессы в низкокалорийных системах, например, газификация низкосортных углей, переработка горючих отходов машиностроительного производства и др.

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

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

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

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

Предметом диссертационной работы является дальнейшее исследование фильтрационного горения твердого топлива с газообразным окислителем методом математического моделирования [9,10].

С учетом вышесказанного, для исследования ФГ предлагается двухтемпературная математическая модель стационарного процесса, учитывающая зависимость теплофизических характеристик фаз от температуры и состава, конечную длину реактора [11,12]. Можно ожидать, что учет зависимости теплоемкостей фаз от температуры и состава не даст точного равенства теплоемкостей потоков фаз, от соотношения которых зависят максимальные разогревы в зоне горения. Модель построена на принципах установления термодинамического равновесия в зоне реакции. Полагалось, что доминирующим процессом при формировании структуры волны горения является процесс теплообмена между фазами. Поэтому, не умаляя ценности детального рассмотрения кинетики химической реакции, для упрощения модели заложили простую химическую схему. Т.к. имеются работы, в которых рассматривалось влияние боковых теплопотерь на характеристики газификации, основное внимание было уделено такому виду теплопотерь как потери тепла через торцы реактора (вынос тепла из реактора продуктами реакции). Расчет позволяет оценивать состав газообразных продуктов сгорания.

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

Заключение диссертации по теме «Химическая физика, в том числе физика горения и взрыва», Салганский, Евгений Александрович

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ.

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

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

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

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

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

ОБОЗНАЧЕНИЯ. а - доля топлива в исходной смеси с - мольный расход топлива, моль/(с*м ) С - теплоемкость,

Дж/(моль*°С) л d - мольный расход диоксида углерода, моль/(с*м ) D - коэффициент диффузии, м2/с Еа - энергия активации, Дж/моль F - удельная поверхность частиц, 1/м о л

G - расход окислителя, м /(час*м )

Gj - приток — сток j-го компонента вещества, моль/(м3*с) h - мольный расход водорода, моль/(с*м ) Н - энтальпия, Дж/моль i - мольный расход инертного компонента, моль/(с*м2) ко - предэкспонент, 1/с

L - длина реактора, м

Lif- координата левого торца реактора, м

Lpr - координата правого торца реактора, м

L* - критическая длина реактора, м ш - мольный расход монооксида углерода, моль/(с*м2)

М - молярная масса, кг/моль п — мольный расход азота, моль/(с*м2)

Nu - число Нуссельта

Р - давление, Па р(с) - давление паров углерода

Рг - число Прандтля

Q - теплоприход от химических реакций, Дж/моль г - характерный размер частиц, м R - универсальная газовая постоянная,

Дж/(моль*°С)

Re - число Рейнольдса t - время, с Т - температура, °С

Tif - температура выгружаемого инертного компонента, °С

Трг - температура выходящих газообразных продуктов сгорания, °С

TgG - температура газовой фазы перед зоной горения, °С р л

Ts - температура конденсированной фазы перед зоной горения, С Т° - температура горения, °С и - линейная скорость горения, м/с

0 - скорость фильтрации воздуха, м/с л w - мольный расход паров воды, моль/(с*м ) л х - мольный расход кислорода, моль/(с*м ) х - текущая координата, м у - доля газообразного компонента в газовой фазе ОС - коэффициент теплообмена, Дж/(м2*с*°С)

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

Г] - степень выгорания топлива

Я - коэффициент теплопроводности, Вт/(м*°С) - динамическая вязкость, Па* с

9 - пористость л -з

Р - плотность, кг/м (моль/м )

ИНДЕКСЫ g - газовая фаза s - конденсированная фаза

If - значение на левом торце реактора рг - значение на правом торце реактора

1 - инертный компонент х - кислород п- азот с - углерод ш - монооксид углерода d - диоксид углерода h - водород w - пары воды j, k, 1 - порядковый номер pd - продукты re - реагенты О - начальное условие

ЗАКЛЮЧЕНИЕ.

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

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

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

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

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

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

Учет кинетики протекающих реакций позволяет описывать режимы с неполным расходованием реагентов, что актуально в случаях: а) низкой реакционной способности твердого топлива (например, графит); б) высоких скоростей фильтрации газообразного окислителя; в) высокой концентрации водяного пара в составе газообразного окислителя и т.д. Зону горения при этом нельзя считать бесконечно узкой из-за конечной скорости протекающих химических реакций. Потребление реагентов может происходить не полностью вследствие больших градиентов температур и "замораживания" соответствующих реакций.

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

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

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

Список литературы диссертационного исследования кандидат физико-математических наук Салганский, Евгений Александрович, 2004 год

1. Алдушин А.П., Сеплярский Б.С. Распространение волны экзотермической реакции в пористой среде при продуве газа // Докл. АН СССР, 1978, Т. 241, № 1, с. 72-75.

2. Алдушин А.П., Сеплярский Б.С. Инверсия структуры волны горения в пористой среде при продуве газа // Докл. АН СССР, 1979, Т. 249, № 3, С. 585-589.

3. Рабинович О.С., Гуревич И.Г. // ИФЖ, 1983, Т. 44, № 1, С. 75-80.

4. Бабкин B.C., Вежба И., Карим Г.А. // Физика горения и взрыва, 2002, Т. 38, № 1,С. 3-11.

5. Алдушин А.П., Мержанов А.Г. Теория фильтрационного горения: общие представления и состояние исследований / Распространение тепловых волн в гетерогенных средах. Новосибирск: Наука, сиб. отд-ие, 1988, С. 952.

6. Манелис Г.Б. Сверхадиабатика // Природа. 1996. № 3-4. С. 43.

7. Выжол Ю.А. Сверхадиабатический режим фильтрационного горения гетерогенных систем: Дис. . канд. физ.-мат. наук. Черноголовка, ОИХФ АН СССР, 1990 г.

8. А.Р. Aldushin, I.E. Rumanov, B.J. Matkowsky. Maximal Energy Accumulation in a Superadiabatic Filtration Combustion Wave // Comb, and flame, 1999, V.l 18, pp. 76-90.

9. Математическое моделирование высокотемпературных процессов в энергосиловых установках/ В.Е. Алимасов, А.Ф. Дрегалин, В.Г. Крюков, В.И. Наумов.-М.:Наука, 1989, 256 с.

10. Ю.Правда. 1986. 31 января.

11. Манелис Г.Б., Салганский Е.А., Фурсов В.П., Глазов С.В., Салганская М.В. Численное исследование газификации углерода в волне фильтрационного горения // Наука-Производству, 2001, №8, С. 28-32.

12. Е.А. Салганский, В.П. Фурсов, С.В. Глазов, М.В. Салганская, Г.Б. Манелис Модель паровоздушной газификации твердого горючего в фильтрационном режиме // Физика горения и взрыва, 2003, Т. 39, № 1, С. 44-50.

13. Алдушин А.П. Теория фильтрационного горения: Дис. . д-ра физ.-мат. наук. Черноголовка, ОИХФ АН СССР, 1981 г.

14. Aldushin А.Р., Seplyarsky. B.S. Sov. Phys. Dokl. 23:483-486 (1992).

15. Aldushin, A. P., in Progress in Aeronautics and Astronautics, vol. 173 (W. A. Sirignano, A. G. Merzhanov, and L. DeLuca, Eds.), American Institute of Aeronautics and Astronautics, 1997. p. 95-115.

16. Schult, D. A., Matkowsky, B. J., Volpert, V. A., Fernandez-Pello, A. C.,

17. Combust. Flame 104:1-26 (1996). 17.Schult, D. A., Matkowsky, B. J., Volpert, V. A., Fernandez-Pello, A. C., Combust. Flame 101:471^90 (1995).

18. Aldushin, A. P., Matkowsky, B. J., Schult, D. A., Combust. Sci. Technol. 00:1-67(1997).

19. Aldushin, A. P., Matkowsky, B. J., Schult, D. A., J. Eng. Math. 31:205-2341997).

20. Aldushin, A. P., Matkowsky, B. J., Schult, D. A., Combust. Flame 107:51-175 (1996).

21. Aldushin, A. P., Matkowsky, B. J., Shkadinsky, K. G., Shkadinskaya, G. V., Volpert, V. A., Combust. Sci. Technol. 99:313-343 (1994).

22. Aldushin, A. P., Matkowsky, B. J., Combust. Sci. Technol. 140:259-2931998).

23. Wahle C.W., Matkowsky B.J., Combust. Flame 124:14-34 (2001).

24. Babkin V.S. Filtration combustion of gases. Present state of affairs and prospects // Pure and Appl. Chem. 1993. V. 65. P. 335-344.

25. К.В. Добрего, С.А. Жданок Физика фильтрационного горения газов / Минск: Ин-т тепло- и массообмена имени А.В. Лыкова НАНБ, 2002, 203 с.

26. Футько С.И., Жданок С.А. Химия фильтрационного горения газов / Минск: Бел. навука, 2004. 319 с.

27. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике / М.: Наука, 1967.

28. Алдушин А.П., Хайкин Б.И., Мержанов А .Г. // Докл. АН СССР, 1974, Т. 215, №3.

29. Алдушин А.П., Хайкин Б.И., Мержанов А.Г. // Докл. АН СССР, 1972, Т. 205, № 2.

30. Баренблатт Г.И., Зельдович Я.Б. // УМН, 1971, Т. 26, № 2.

31. А.П. Алдушин. Теплопроводностный и конвективный режимы горения пористых систем при фильтрации теплоносителя // Физика горения и взрыва, 1990, Т. 26, № 2, С. 60-68.

32. Новожилов Б.В. Нестационарное горение твердых ракетных топлив / М.: Наука, 1973.

33. Алдушин А.П., Каспарян С.Г. // Физика горения и взрыва, 1981, Т. 17, № 6.

34. Е.М. Тонкопий, Г.Б. Манелис, С.В. Куликов Численная модель гетерогенного горения в пористой среде // Химическая физика, 1992, Т. 11, № 12, С. 1649-1654.

35. Е.М. Тонкопий, Г.Б. Манелис, С.В. Куликов Численная модель горения углерода в пористой среде с учетом кинетики химических превращений // Физика горения и взрыва, 1993, Т. 29, № 3, С. 136-139.

36. Беккер А.В., Полианчик Е.В., Глазов С.В. Неединственность стационарных режимов при фильтрационном горении углерода // Тез. докл. XXI Всероссийский симпозиум молодых ученых по химической кинетики. МО, Клязьма. 2003. С. 76

37. Беккер А.В. Макрокинетическое моделирование сверхадиабатического фильтрационного горения углеродосодержащих материалов: Дис. . канд. физ.-мат. наук. Черноголовка, ИПХФ РАН, 2004 г.

38. Howard J.B., Williams G.C. and Fine D.H. Kinetics of Carbon Monoxide Oxidation in Postflame Gases // International Fourteenth Symposium on Combustion. The Combustion Institute. Pittsburgh. P.A. 1973. P. 975

39. Wakao N., Kaguil S. Heat and mass transfer in porous media. Gordon and Breach Science publ., 1982.41. «Термодинамические свойства индивидуальных веществ» под ред. В.П.Глушко и др., изд. третье, Москва, «Наука», 1978.

40. Калиткин Н.Н. Численные методы. М.: Наука, 1978.

41. Канторович Б.В. Введение в теорию горения и газификации твердого топлива. М.: Металлургиздат, 1960.

42. Глазов С.В. и др. Экспериментальное и теоретическое исследование процессов газификации низкокалорийных горючих материалов в сверхадиабатическом режиме.// Научный отчет. ИПХФ РАН, Черноголовка, 1999.

43. Вещунов И.С. Экспериментальное исследование фильтрационного горения составов с активированным углем при различном соотношении компонентов. Дипломная работа. ИПХФ РАН, Черноголовка, 2002.

44. Haslan R.T. // Industrial Engineering Chemistry, V. 15, P. 679, 1983

45. Лайров H.B., Розенфельд E.H., Хаустович Г.Н. Процессы горения топлива и защита окружающей среды / Москва: Металлургия, 1981

46. Howard I.B., Williams G.S., Fine D.H. Kinetics of carbon monoxide oxidation in post flame gases // Fourteenth Symposium on Combustion, The Comb. Inst., Pittsburg PA, 1973, P. 975

47. Baulch P.L. // Combustion and Flame, 98:59-79, 1994

48. Govind R., Shan J. Modeling and simulation of an entrained flow coal gasifire // AIlhE Journal, 1984,V. 30, № 1, P. 79-91

49. Bews I.M., Hayhurst A.N., Richardson S.M., Taylor S.G. // Combustion and Flame, 124:231-245 (2001)

50. Field M.A., Fill D.V., Morgan B.B., Hawksley P.G.W. Combustion of pulverized coal / BCURA, leaher head, England, 1967

51. Dean A.J., Davidson D.F., Hanson R.R. // J. Phys. Chem. 95:183 (1983)

52. Нечаев B.H., Рогайлин M.H., Новые методы сжигания топлива и вопросы горения / М.: Наука, 1972

53. Головина Е.С., Климов А.А., Об истинной кинетической константе гетерогенной газификации С+С02 // ФГВ, 1999, Т. 35, № 4, С. 48-51

54. Lee S., Angus J.C., Edwards R.V., Gardrer N.C. // AIlhE Journal 30:583 (1981)

55. Предводителев A.C. и др. Горение углерода / Изд. АН СССР, M.-JL, 1949

56. Fritz O.W., Huttinger K.J., Active sites and intrinsic rates of carbon-gas reaction // Carbon, 1993, V. 31, № 6, P. 923-930

57. Young B.C. // Proceeding of the International conference on coal science, Verlag Glue Kauf Essen, 1981, P. 260

58. Glarborg P., Alzueta M.V. // Combustion and Flame, 115:1 (1998)

59. Libby P.A., Blake T.R. // Combustion and Flame, 36:139 (1979)

60. Головина E.C. // Физика горения и взрыва, 2002, Т. 38, № 4, С. 25-34

61. Lizzo A., Hong Jiang, Radovic R., On the kinetics of carbon (char) gasification: reconciling model with experiments // Carbon, 1990, V. 28, № 1, P. 7-19

62. Хитрин JI.H., О нестационарных диффузионно-кинетических процессах в гетерогенных системах // ХТТ, 1967, № 6, С. 110-117

63. Бабий В.Н., Буланов Д.В., Артемьев В.П., Макрокинетические константы процесса разложения NO, С02 и Н20 в восстановленной пылевидной среде // XII Симпозиум по горению и взрыву, Черноголовка, 2000, Технология, С. 120

64. Самсонов В.Г., Головина Е.С., Арабаджиев Б.Г., Якушев А.А., Газификация углерода водяным паром при температурах 1600-2400К // VI Всесоюзный Симпозиум по горению и взрыву, (1980, Алма-Ата) ОПХФ АН СССР, Черноголовка, 1980, С. 66-68

65. Cho S.V., Yetter R.A., Dryer F.L. // J. Comp.Phys. 102:160 (1992)

66. Yetter R.A., Dryer F.L., Rabitz H. // Comb. Science Tech., 79:97-128 (1991) 69.Smith I.W. The intrinsic reactivity of carbons to oxygen // Fuel, 1978, V. 57,1. P. 409-414

67. Головина E.C. Высокотемпературное горение и газификация углерода / М.: Энергоатомпродукт, 1983

68. Li С., Brown Т.С. Carbon oxidation kinetics from evolved carbon oxide analysis during temperature- programmed oxidation // Carbon, 2001, V. 39, P. 725-732.

69. Головина E.C. Об окислении некоторых углей // Известия АН СССР. ОТН. 1949, № 9, С. 1343-1351

70. Спицин В.Б. Паровоздушная газификация графита в волне фильтрационного горения. Дипломная работа. ИПХФ РАН, Черноголовка, 2001.