Разработка метода сокращения механизмов реакций и анализ чувствительности для моделирования процессов горения в энергоустановках тема диссертации и автореферата по ВАК РФ 05.07.05, кандидат технических наук Никандрова, Марина Викторовна

• Никандрова, Марина Викторовна
• кандидат технических наук
• 2007

Оглавление диссертации кандидат технических наук Никандрова, Марина Викторовна

Введение.

1. Обзор литературы и постановка задачи исследования.

1.1. Актуальность проблемы моделирования процессов горения.

1.2. Уравнения формальной химической кинетики и методы их решения

1.3. Методы анализа чувствительности.

1.4. Методы сокращения механизмов реакций.

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

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

2.1. Модель реактора идеального смешения и метод решения.

2.2. Алгоритм определения коэффициентов чувствительности.

2.3. Метод "зацепления" для сокращения механизма реакций.

3. Программное обеспечение и формирование полного механизма реакций

3.1. Описание программного комплекса.

3.2. Апробация модели реактора и метода определения коэффициентов чувствительности.

3.3. Формирование полного механизма реакций для реагирующей среды

Н + О + S + (N)".

4. Генерация и исследование сокращенных механизмов для реагирующей среды "Н + О + S + (N)".

4.1. Анализ чувствительности для полного механизма реакций.

4.2. Формирование и исследование LS - механизмов.

4.3. Формирование GS - механизмов для реагирующей смеси "H2S + воздух"

4.4. Апробация GS - механизма на задаче о фронте пламени.

4.5. Анализ чувствительности для сокращенного механизма.

Введение диссертации (часть автореферата) на тему «Разработка метода сокращения механизмов реакций и анализ чувствительности для моделирования процессов горения в энергоустановках»

Уменьшение загрязнения атмосферного воздуха токсичными веществами, выделяемыми промышленными предприятиями, автомобильным транспортом и т.д. является одной из важнейших проблем на сегодняшний день. Кроме того, актуальной является проблема сокращения расхода топлива на транспорте и в энергетических установках. Чтобы отвечать предъявляемым требованиям по экономичности и экологичности, необходимо оптимизировать процессы горения в высокотемпературных агрегатах. Желательно это осуществлять на стадии проектирования, путем математического моделирования. Такой подход является оптимальным не только с точки зрения удобства, но и с финансовой стороны. Описание процессов протекающих в таких установках затрагивает кинетику химических реакций, а следовательно предполагает работу с большими и сложными реакционными механизмами, включающими сотни реакций и десятки веществ. Моделирование, а также анализ таких механизмов является весьма сложной проблемой, которая влечет за собой ряд других задач, таких как развитие методов по сокращению механизмов, анализ чувствительности относительно констант скоростей химических реакций и т.д. Особое значение придается возможности внедрения таких разработок в программные оболочки, отвечающие не только простоте и удобству использования, но и современным тенденциям в компьютерных технологиях.

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

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

Во второй главе на основе уравнений химической кинетики в экспоненциальной форме была сформирована модель реактора идеального смешения. Был развит аналитический метод определения коэффициентов чувствительности газофазных реакций. Для сокращения сложных механизмов реакций газофазного горения был разработан метод "зацепления", позволяющий формировать скелетные механизмы (S - механизмы) с различной степенью приближения к полному механизму. В отличие от ранее известных подходов, данный метод генерирует достаточно компактные S - механизмы как для одного режима горения (LS - механизм), так и для заранее заданной области режимов (GS - механизм).

В третьей главе описан используемый программный комплекс, состоящий из двух инвариантных программ RISNEW и RISTV. В программе RISNEW был реализован метод определения коэффициентов чувствительности. Данный метод был апробирован путем сравнения выполненных расчетов для реагирующей среды "Н + О" с ранее известными данными. В свою очередь разработанный метод "зацепления" был внедрен в инвариантную программу расчета реактора идеального смешения (RISTV). На основании обобщения данных из ряда известных работ, для смеси "Н + S + О + (N)" был сформирован полный реакционный механизм горения (С - механизм), содержащий 25 веществ и 118 обратимых реакций.

Четвертая глава посвящена численным исследованиям. В частности были проведены расчеты:

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

- по формированию сокращенного LS - механизма реакций;

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

- по определению коэффициентов чувствительности для GS - механизма;

- по апробации GS - механизма на модели фронта горения.

Был также проведен анализ полученных численных данных.

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

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

Заключение диссертации по теме «Тепловые, электроракетные двигатели и энергоустановки летательных аппаратов», Никандрова, Марина Викторовна

ВЫХОД

Рис. 3.5. Блок-схема программы расчета RISTV

3.2. Апробация модели реактора и метода определения коэффициентов чувствительности

С целью обоснования достоверности нашей программы RISNEW было проведено сравнение с результатами, полученными по инвариантной программе "Aurora", включенной в известный пакет CHEMKIN [70], позволяющей рассчитывать характеристики реактора идеального смешения и определять коэффициенты чувствительности. Следует отметить, что в "Aurora", в отличие от RISNEW:

- применяется традиционная форма уравнений химической кинетики;

- используется другая база данных по веществам (с другой формой полиномиальных зависимостей Hj =f(T); S° = f{T))\

- уравнение энергии применяется в дифференциальной форме;

- частные производные Якобиана и коэффициенты чувствительности определяются численно.

Сравнение было выполнено для реагирующей среды "Н + О", включающей 9 веществ (О, 02, Н, Н2, ОН, Н20, Н02, Н202, N2) и 19 реакций, приведенных в таблице 3.1. В ней также представлены коэффициенты аррениусо-вой зависимости констант скорости прямых направлений реакций: k+s = A*T"S ехр(-£* /R0T), с размерностями для [E+s] = кал/моль,

R0 =1,987 кал/(моль-К), [Т] = К. Соответствующие скорости реакций имели размерности [Ws+] = моль/см3-сек.

При выполнении сравнительных расчетов значения коэффициентов ,п*,Е* обеих программ были идентичны, но использовались разные базы о термодинамических свойствах веществ (в "Aurora" применялась база данных [72], а в RISNEW - база данных TTI [2]).

ЗАКЛЮЧЕНИЕ

По результатам выполненного диссертационного исследования можно сделать следующие выводы:

1. Для сокращения сложных механизмов реакций был разработан метод "зацепления". Данный метод создан для условий реактора идеального смешения и позволяет на базе некоторого полного механизма формировать компактные S - механизмы как для одной точки (LS - механизмы), так и для заданной области изменения параметров аок, Т, Р реагирующей среды (GS - механизмы).

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

3. На основе базового программного обеспечения были созданы инвариантные программы расчета химически неравновесных процессов в реакторе: RISNEW (с расчетом коэффициентов чувствительности ) и RISTV (с генерацией LS - и GS - механизмов). Апробация разработанных методов и программ подтвердила их достоверность, показала хорошее качественное и количественное совпадение с известными результатами, полученными другими авторами.

4. На базе работ [77, 83] был сформирован сложный реакционный механизм (С - механизм) для реагирующей среды "Н + О + S +(N)", включающий 118 реакций и 25 веществ. Для этой среды практически отсутствуют публикации по созданию сокращенных механизмов и анализу чувствительности.

5. Для этой реагирующей среды (реагенты "H2S + воздух") были выполнены следующие численные исследования: а) кинетический расчет по полному механизму реакций с получением коэффициентов чувствительности. На основе этих данных, при использовании традиционного подхода [36, 52, 70, 87, 93, 96] количество реакций удалось сократить почти вдвое, однако ни одно из участвующих веществ не было исключено. Кроме того, проводить собственно анализ чувствительности оказалось затруднительно, ввиду громоздкости С -механизма. Поэтому было предложено выполнять такой анализ по сокращенному механизму, полученному после применения метода зацепления. б) по формированию LS - механизмов для проверки разработанного метода "зацепления", оценки влияния показателя сокращения механизма и возможности включения заданных веществ. Была показана возможность генерации сокращенных механизмов с различной степенью приближения к С - механизму (в зависимости от показателя при этом исходный механизм может быть сокращен в несколько раз (при приемлемой погрешности). Была продемонстрирована способность включения в сокращенный механизм заданного вещества (что является важным при прогнозировании вредных выбросов); в) по генерации GS - механизмов для области изменения параметров аок= l.l,5; Т- 1200-2000 К; Р = 1.5атм., свойственной работе ряда энергоустановок. Показатель сокращения варьировался в интервале С, = 0.0,4. При этом достигалось сокращение исходного С - механизма в 3.4 раза по реакциям и примерно в два раза по веществам (при сокращении объема "кинетических" расчетов приблизительно в 15.25 раз). Погрешность GS - механизмов по сравнению с С - механизмом оставалась незначительной при £ =0,1; 0,2; 0,3, но становилась неприемлемой при C, = 0,4. Для практического применения рекомендуется механизм GS (£ = 0,3). г) в рассматриваемой области изменения аок, Т, Р для механизма GS(£ = 0,3) был проведен анализ чувствительности с предварительным определением коэффициентов Ris, RTs. Данные, полученные в результате выполнения этих расчетов позволили выделить набор "критических" реакций для различных веществ и зон области исследования, а также определить допустимые погрешности для констант скоростей реакций.

6. Апробация = 0,3) - механизма проводилась на задаче расчета горения во фронте пламени, путем сравнения данных, полученных по С - и GS - механизмам. Результаты сравнения позволили констатировать, что в исследуемой области изменения аок, Т, Р GS(^ = 0,3) - механизм можно применять для прогнозирования характеристик процессов горения в высокотемпературных энергоустановках.

Список литературы диссертационного исследования кандидат технических наук Никандрова, Марина Викторовна, 2007 год

1. Абдуллин A.JL Математическое моделирование процессов во фронте пламени с использованием методологии НРИС/ Авиационная техника, №4, 2003, стр. 41-44.

2. Алемасов В.Е., Дрегалин А.Ф., Тишин А.П., Термодинамические и теплофизические свойства продуктов сгорания. Справочник Т.1, ВИНИТИ, Москва, 1971.

3. Алемасов В.Е., Дрегалин А.Ф., Черенков А.С. Основы теории физико-химических процессов в тепловых двигателях и энергоустановках. М.:Химия, 2000, 520 с.

4. Аствацатуров А.Е., Булыгин 10.И., Жигулин И.Н., Ладоша Е.Н., Сакаев Э.К., Яцепко О.В. Компьютерная модель термогазодинамики и химических превращений в поршневом двигателе. Известия вузов. Сев.-Кавк. Регион. Технические науки, 2001. №1, с. 75-82.

5. Бсзменова Н.В., Шустов С.А., Гидаспов В.Ю., Иванов И.Э. Моделирование течения продуктов сгорания в соплах ЖРДМТ. // Математическое моделирование, 1999, Т. 11, №6, стр.45-51.

6. Белов Г.В., Иориш B.C., Юнгман B.C. Моделирование равновесных состояний термодинамических систем с использованием IVTANTHERMO для Windows. // Теплофизика высоких температур. -2000, т. 38, No. 2, с. 191 196.

7. Блинов В.И., Худяков Г.И. Диффузионное горение жидкостей. М.: Изд-во аи СССР, 1961. - 207 с.

8. Вильяме Ф.А. Теория горения. М.: Наука, 1971.-616 с.

9. Волков В.А., Гидаспов В.Ю., Пирумов У.Г. и др. Численное моделирование течений реагирующих газокапельных и газовых смесей в экспериментах по воспламенению метонола. // Теплофизика высоких температур, Том 36, №3, 1988, стр. 424-434.

10. Воронецкий А.В., Макаров Д.В., Скибин А.П., Филимонов JI.A., Югов В.П. Расчет процессов горения и охлаждения в малогабаритной камере сгорания//Математическое моделирование, 1999, т.11, №4, стр.29-36.

11. Давлетшин Р.Ф., Яценко О.В. КИНКАТ автоматизированная система разработки сложных радиационно-кииетических моделей. Мат. моделирование процессов управления и обработки информации: Межвед. сб. -М.: МФТИ, 1993.-С. 113-123.

12. И.Двигатели внутреннего сгорания: Устройство и работа поршневых и комбинированных двигателей./ Учебник для студентов вузов/ под ред. А. С. Орлина, М. Г. Круглова.З-е изд., перераб. и доп. М.: Машиностроение, 1980.

13. М.Демидов Л.Г., Шандыба В.А., Щеглов ГШ. Горение и свойства горючих веществ. М.: Химия, 1981, 272 с.

14. Дмитров В.И. Простая кинетика. Новосибирск: 11аука, 1982. -381 с.

15. Евлампиев А.В., Фролов С.М., Басевич В.Я., Беляев А.А. Глобальные кинетические механизмы на основе ламинарного пламени для моделирования турбулентных реагирующих течений. 4.IV. Диффузионное горение, Химическая физика, 2001. Т.20. №11. С.21.

16. Егоров Б.В., Комаров В.Н., Маркачев Ю.Е. и др. Кластеризация продуктов горения в следе сверхзвукового самолета//Математическое моделирование, 1999, т. 11, №4, стр. 70-82.

17. Зельдович Я.Б., Баренблатт Г.И., Либрович В.Б., Махвиладзе Г.М. Математическая теория горения и взрыва. М.: Паука, 1980, 478 с.

18. Калиткин Н.Н. Численные методы. Гл. ред. Физико-математической литературы изд-ва «Наука», М., 1978, 512 с.

19. Камзолов В.И., Пирумов У.Г., Расчет неравновесных течений в соплах. //Изв. АН СССР, Механика жидкости и газа, 1966, No 6, стр. 25-33.

20. Кнорре Г.Ф. Топочные процессы. М.: Наука, 1951, 328 с.

21. Козлов Г.И., Кнорре В.Г. Использование одноимпульсной ударной трубы для исследования суммарной кинетики термического разложения метана. Инженерно физический журнал, №10, т.4., 1961, стр.4.

22. Кондратьев В.Н., Никитин Е.Е. Кинетика и механизм газофазных реакций.-М.: Наука, 1974. 512 с.

23. Кондратьев В.Н., Константы скорости газофазных реакций. Справочник, М.-Наука, 1974, 512 стр.

24. Кондратьев В.Н. Определение констант скоростей газофазных реакций. -М.: Наука, 1971.-96 с.

25. Короленок Е.В., Нагорнов О.В. Моделирование приземных концентраций озона городского региона//Математическое моделирование, 2002, т. 14, №4, стр.80-94.

26. Крюков В.Г., Абдуллин A.JL, Дуригон А. Применение метода сплайн-функции для решения задач химической кинетики.//Вестпик КГТУ, Казань, 2004. Том 3.

27. Крюков В.Г., Наумов В.И., Абдуллин A.J1., Демин А.В., Исхакова P.J1. Математическое моделирование реагирующих течений на базе реакторного подхода//Физико-химическая кинетика в газовой динамике. 2004. Том 2.

28. Крюков В.Г., Наумов В.И., Абдуллин А.Л., Демин А.В., Тринос Т.В. Горение и течение в агрегатах энергоустановок. Москва, Янус-К, 1997. -304 с.

29. Маликов К.Ю., Романюха Н.Ю., Четвертушкин Б.Н. Образование окислов азота в промышленных печах струйно-факельного нагре-ва.//Математическое моделирование, 1998, т. 10, №9, стр.41-52.

30. Мальцев В.М., Мальцев М.И., Кашиоров Л.Я. Основные характеристики горения. М.: Химия, 1977, 320 с.

31. Математическое моделирование горения распыленного угля./ Труды IV школы-семинара "Проблемы тепломассопереноса в энергетических установках", Казанский филиал РАН, 2004. -стр.43-55.

32. Мейтис J1. Введение в курс химического равновесия и кинетики: Пер. с англ.-М.: Мир, 1984.-484 с.

33. Никандрова М.В., Крюков В.Г., Исхакова P.JI. Аналитическое определение коэффициентов чувствительности реакций для условий реактора идеального смешения.//Физико-химическая кинетика в газовой динамике. 2006. Том 4.

34. Никандрова М.В., Крюков В.Г., Исхакова P.JI. Сокращение механизмов реакций методом "зацеплсния".//Физико-химическая кинетика в газовой динамике. 2007. Том 5.

35. Оран Э., Борис Дж. Численное моделирование реагирующих потоков: Пер. с англ. М.: Мир, 1990. - 660 е., ил.37.0сновы практической теории горения. Под ред. Померанцева В.В. -Ленинград, Энергоатомиздат, 1986, 312 с.

36. Пирумов У.Г. Вычислительная механика и проблемы охраны окружающей среды.//Математическое моделирование, 2002, т. 14, №4, стр.80-94.

37. Пирумов У.Г., Росляков Г.С. Газовая динамика сопел. М.: Наука, Гл. изд. физ.мат. литературы, 1990, 368 с.

38. Рэмсден Э.Н. Начала современной химии: Справ, изд.: Пер. с англ./Под ред. В.И. Барановского, А.А. Белюстина, А.И. Ефимова, А.А. Потехина-Л.: Химия, 1989.-784 с.

39. Самарский А.А., Михайлов А.П. Математическое моделирование. Идеи. Методы. Примеры. -М.: Наука. Физматлит, 1997. 320 с.

40. Скребков О.В., Каркач С.П. Колебательная неравновесность и электронное возбуждение в реакции водорода с кислородом за ударной вол-ной.//Физико-химическая кинетика в газовой динамике. 2004. Том 2.

41. Талантов А.В. Основы теории горения. Казань, КАИ, 1975, 252 с.

42. Турчак Л.И. Основы численных методов. М.: Наука, 1987. - 320 с.

43. Уилкинсон Д.Х., Алгебраическая проблема собственных значений, М.: Наука, 1970.

44. Физико-химические процессы в газовой динамике. Компьютеризированный справочник в 3-х томах, Т.1: Динамика физико-химических процессов в газе и плазме. Под ред. Черного Г.Г., Лосева С.А., Москва, Изд. МГУ, 1995,350 с.

45. Физико-химические процессы в газовой динамике. Компьютеризированный справочник в 3-х томах, Т.2: Физико-химическая кинетика и термодинамика. Под ред. Черного Г.Г., Лосева С.А., Москва, Изд. МГУ, 2002, 368 с.

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

47. Фролов Ю.В. Теория горения и взрыва. М.: Наука, 1981,412 с.

48. Химия горения: Пер. с англ./ Под ред. У. Гардинера, мл. М.: Мир, 1988. -464с., ил.

49. Эмануель Н. М., Кнорре Д.Г., Курс химической кинетики. 4-е изд., пере-раб. и доп. М.: Высшая школа, 1984.-463 с.

50. Alzueta M.U., Bilbao R., Glarborg P., Inhibition and Sensitization of Fuel Oxxidation by S02. Combustion and Flame, Vol. 127, pp. 2234-2251, 2001.

51. Baulch D.L., Drysdale D.D., Home D.G., Lloyd A.C. Evaluated kinetic data for high temperature reactions, Butterworth, London, 1972.

52. Baulch D.L., Duxbury J., Grant S J., Montague D.C. Evaluated kinetic data for high temperature reactions, Journal Phys.Chem. Ref. Data, Vol.10, Supplement No. 1., 1981.

53. Blom J.G. and Verwer J.G., A Comparison of Integration Methods for Atmospheric Transport-Chemistry Problems. Journal of Computational and Applied Mathematics, No 126, pp. 381-396, 2000.

54. Brown R.C., Kolb C.E., Yetter R.A., Dryer F.L., Rabitz H. Kinetic modeling and sensitivity analysis for B/H/O/C/F combination systems. Combustion and Flame, Vol.' 101, pp. 221-238, 1995.

55. Burcat A., Third Millennium Ideal Gas and Condensed Phase Thermochemi-cal Database for Combustion. Technion Aerospace Engineering (TAE) Report # 867 January, Faculty of Aerospace Eng. Technion Israel Institute of Technology, 2001.

56. Burks T.L., Oran E.S. A Computational Study of the Chemical Kinetics of Hydrogen Combustion, NRL Memorandum Report 4446, Naval Research Laboratory, Washington, DCAD A094384., 1981.

57. Chen J.Y. A General procedure for constructing reduced reaction mechanisms with given independent reactions. Combustion Science and Technology, Vol. 57, pp.89-94, 1988.

58. Constanza V., Seinfeld J.H. Stochastic sensitivity analysis in chemical kinetics, Journal Chem. Phys., Vol.74, pp.3852-3858, 1981.

59. Dickinson R.P., Gelinas R.J. Sensitivity analysis of ordinary differential equation systems. Journal Comput. Phys., Vol.21, pp. 123-143, 1976.

60. Eaton A.M., Smooth L.D., Hill S.C., Eatough C.N. Components, Formulations, Solutions, Evaluation and Application of Comprehensive Combustion models. Progress in Energy and Combustion science 25, 1999, pp. 387 436.

61. Eggels R.L., Goey L.P.H. Mathematically reduced reaction mechanisms applied to adiabatic flat hydrogen/air flames. Combustion and Flame, Vol. 100, pp. 559-570, 1995.

62. Forsth M. Sensitivity analysis of the reaction mechanism for gas-phase chemistry of H2 = 02 mixtures induced by a Hot Pt surface. Combustion and Flame, Vol. 130, pp.241-260, 2002.

63. Frenklach M., Bornside D.E. Shock-Initiated Ignition in methane-propane mixtures, Combustion and Flame, Vol. 56, pp. 1-27, 1984.

64. Gear C.W. DIFSUB for solutions of ordinary differential equations. Communications of ACM, Vol.14, Nu3, pp. 185-190, 1971.

65. Gear C.W. Numerical Initial Value Problems in Ordinary Differential Equations, Prentice-Hall, Englewood Cliffs, NJ, 1971.

66. Gelinas R.J., Vajk J.P. Systematic sensitivity analysis of air quality simulation models, Science Applications Inc., Pleasanton, CA. PB80-112162., 1978.

67. Glarborg P., Miller J.A., Kee R.J. Kinetic Modeling and Sensitivity Analysis of Nitrogen Oxide Formation in Well-Stirred Reactors. Combustion and Flame, Vol.65, pp. 177-202, 1986.

68. Glaude P. A., Battin-Leclerc F., Fournet R., Warth V., Come G.M., Scacchi G. Construction and simplification of a model for the oxidation of alkanes. Combustion and Flame, Vol. 122, pp.451 -462, 2000.

69. Gordon S, McBridc B.J., NASA SP-273 Computer Program for Calculation of Complex Chemical Equilibrium Compositions, Rocket Performance, Incident and Reflected Shocks and Chapman-Jouguet Detonations, NASA, Washington, 245p., 1971.

70. Hampson R.F. Chemical Kinetic and Photochemical data sheets for atmospheric reactions, Report No. FAA/EE-80-17, Department of Transportation, Washington, DC AD H091631 ., 1980.

71. Jones W.P., Lindstedt R.P., Global Reaction Schemes for Hydrocarbon Combustion. Combustion and Flame, Vol. 73, pp.233-242, 1988.

72. Ju Y., Niioka T. Reduced kinetic mechanism of ignition for non-premixed hydrogen/air in a supersonic mixing layer. Combustion and Flame, Vol.99, Issue 2, pp.240-246, 1994.

73. KaIamatianos S., Park Y.K., Vlachos D.G. Two-parameter continuation algorithms for Sensitivity analysis, parametric dependence, reduced mechanisms, and stability criteria of ignition and extinction. Combustion and Flame, Vol. 112, pp.45-61, 1998.

74. Kee R.J, Rupley F.M., Meeks E., Miller J.A., "CHEMKIN: A Software package for the analysis of gas-phase chemical and plasma kinetics", Sandia National Laboratories Report SAND96-8216, CA, 2000.

75. Kramer M.A, Kee R.J, Rabitz H. CHEMSEN: A computer code for sensitivity analysis of elementary chemical reaction models, Sandia Report SAND-82-8230, Sandia National Laboratories, Albuquerque, NM, 1984.

76. Kramer M.A, Rabitz H, Calo J.M, Kee R.J. Sensitivity analysis in chemical kinetics: Recent developments and computational comparisons, International Journal of Chemical Kinetics, Vol.16, pp.559-578, 1984.

77. Kuo, K.K, Principles of Combustion, J, Wiley & Sons, Singapore, 1986.

78. Lam S.H. Using CSP to understand complex chemical kinetics. Combustion Science and Technology, Vol.89, Issue 5-6, pp.375-404, 1993.

79. Lam S.H, Goussis D.A. The CSP method for simplifying kinetics. International Journal of Chemical Kinetics, Vol.26, pp.461, 1994.

80. LEEDS Reaction Kinetics Database. School of Chemistry, University of Leeds, 2004.

81. Lindstedt R.P, Maurice L.Q, "Detailed Chemical Kinetic Model for Aviation Fuels", Journal of Propulsion and Power, Vol. 16, No 2, pp. 187-195, 2000.

82. Maas U, Pope S.B. Simplifying chemical kinetics: Intrinsic low-dimensional manifolds in composition space. Combustion and Flame, Vol. 88, Issue 3-4, pp. 239-264, 1992.

83. Massias A, Diamantis D, Mastorakos E, Goussis D.A. An algorithm for the construction of global reduced mechanisms with CSP data. Combustion and Flame, Vol. 117, pp.685-708, 1999.

84. MilIer D, Frenklach M, Sensitivity analysis and parameter estimation in dynamic modeling of chemical kinetics, Int. J. Chem. Kin, Vol.15, pp. 677-696, 1983.

85. Naumov V.I, Krioukov V.G, Abdullin A.L, 2003. Chemical kinetics software system for the propulsion and power engineering In. 41-st Aerospace

86. Sciences Meeting and Exhibit. Reston. VA. American Institute of Aeronautics and Astronautics. Vol 1, pp. 1-12.

87. Nieckele A.O., Naccache M.F., Gomes M.S.P., Carneiro J.E., Serfaty R. Evaluation of models for combustion processes in a cylindrical furnace. ASME-IMECE, International Conference of Mechanical Engineering, New York, 2001.

88. Niemann H., Schmidt D., Maas U. An efficient storage scheme for reduced chemical kinetics on orthogonal polynomials. Eng. Math., Vol. 31, pp. 131, 1996.

89. Norris A.T., Pope S.B. Modeling of extinction in turbulent diffusion flames by the velocity-dissipation-composition PDF method. Combustion and Flame, Vol. 100, pp. 211-220, 1995.

90. Peters N., Kee R.J. The computation of stretched laminar methane-air diffusion flames using a reduced four-step mechanism. Combustion and Flame, Vol. 68, pp. 17-29, 1987.

91. Rabitz H., Kramer M., Dacol D. Sensitivity Analysis in Chemical Kinetics, Ann. Rev. Phys. Chem., 34, pp. 419-461, 1983.

92. Stoer J., Bulirsch R. Introduction to Numerical analysis, Chpt.7, Springer-Verlag, New-York, 1980.

93. Tianfeng L., Yiguang J., Lam C.K. Complex CSP for chemistry reduction and analysis. Combustion and Flame, Vol. 126, pp. 1445-1455, 2001.

94. Tilden J.W., Constanza V., McRae G.J., Seinfeld J.H., Sensitivity analysis of chemically reacting systems in: Modeling of chemical reaction systems, Ebert K.H., Deuflhard P., Jager W., Springer Verlag, New York, 69-91,1981.

95. Verwer J.G., Hundsdorfer W.H., Blom J.G. Numerical time integration for air pollution models, Mas-r-9825, Centrum Voor Wis Kunde en Informatica, Amsterdam, 1998.

96. VIachos D.G. Mathematically reduced reaction mechanisms applied to adia-batic flat hydrogen/air flames. Chem. Eng. Sci., Vol. 51, pp. 3979, 1996.

97. Warnatz J. Rate coefficients in the C/H/O system, in: Combustion Chemistry, W.C. Gardiner (ed.), Springer-Verlag, New York, pp.173-197, 1984.

98. Westbrook C.K., Dryer F.L. Chemical Kinetic Modeling of Hydrocarbon Combustion, Prog./Energy Combust. Sci., Vol.10, pp. 1-57, 1984.

99. Westley F. Tables of experimental rate constants for chemical reactions occurring in combustion (1971-1977), Report NBSIR-81-2254, National Bureau of Standards, Washington, DC.PB81-205429., 1981.