CC BY
123
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ АДИАБАТИЧЕСКОГО ТЕПЛОВОГО ВЗРЫВА ДЛЯ РЕАКЦИИ ОКИСЛЕНИЯ ВОДОРОДА
Е.Е. Пигасов, В.К. Рябинин, Ю.М. Ковалев
В данной работе представлены результаты математического моделирования адиабатического реактора, заполненного газовой смесью водорода с воздухом. Несмотря на то, что реакция окисления (горения) водорода очень широко используется в практических целях, и существует достаточно большое количество публикаций, разработка кинетических схем горения водорода является актуальной и в настоящее время. В работе были апробированы различные кинетические схемы горения водорода и было показано, что наиболее точное описание экспериментальных данных удается получить, используя полную кинетику окисления. Показано, что использование сокращенной кинетики окисления допустимо лишь для описания высокотемпературной асимптоты, а упрощение механизма до одной брутто-реакции пригодно только для приблизительных оценок.
Ключевые слова: водород, химическая кинетика, горение и взрыв, период индукции.
Введение
Несмотря на то, что реакции окисления (горения) очень широко используются в практических целях, и существует достаточно большое количество публикаций, разработка кинетических схем горения водорода является актуальной и в настоящее время.
На протяжении длительного времени водород рассматривается как эффективное рабочее тело для двигательных установок вследствие его исключительной химической активности в сочетании с воздухом и кислородом в качестве окислителя, а также благодаря отсутствию токсических продуктов сгорания. Имеется ряд прямых свидетельств использования чистого водорода в качестве горючего для прямоточных воздушно-реактивных двигателей (ПВРД) и для камер сгорания жидкостных ракетных двигателей (ЖРД). Весьма часто [1] можно встретить упоминание об использовании водорода, как присадки к углеводородным горючим, особенно применительно к процессам горения в двигателях внутреннего сгорания (ДВС) с искровым зажиганием и воспламенением от сжатия (дизельные моторы).
Столь широкое использование водорода как энергоносителя в сочетании с очевидными проблемами безопасности технических устройств, где предполагается применение или появление водорода, требует ясного представления об его взрывных характеристиках при горении в специально организованном процессе (как в камере сгорания) или при аварийной ситуации (как при тяжелой аварии с потерей охлаждающего компонента на АЭС). При рассмотрении рабочего процесса в ПВРД, ЖРД, ДВС следует заботиться о высокой эффективности и устойчивости преобразования горючей смеси в камерах сгорания (ПВРД, ЖРД, ДВС) и газогенераторах (ЖРД). Особое внимание здесь следует обращать на переходные режимы горения при изменении рабочих параметров. Дополнительные вопросы появляются в связи с проблемой использования в двигательных установках комбинированных типов горючего, представляющих собой смесь водорода с углеводородами [2].
Одним из основных параметров, характеризующих реагирующую смесь, является адиабатический период индукции, расчет которого и рассмотрен в данной работе для различных кинетических механизмов окисления водорода.
Постановка задачи
Рассмотрим адиабатический реактор [3] - реактор, в котором отсутствует всякий тепло-и массообмен с окружающей средой, следовательно все тепло, получаемое в ходе экзотермических реакций окисления идет на разогрев системы. В таком реакторе тепловой взрыв происходит при любых начальных условиях. Отрезок времени от начала процесса до момента резкого увеличения температуры носит название адиабатического периода индукции (задержка воспламенения).
Расчет адиабатического периода индукции сводится к решению системы обыкновенных дифференциальных уравнений для описания поведения температуры смеси и концентраций реагентов:
где Т - температура, £ - время, г - число реакций, т - число реагентов, Н - энтальпия, Ср - теплоемкость, С - концентрация, М - молярная масса, Е - энергия активации, ко,Ь - кинетические параметры, г' - стехиометрический коэффициент для ]-компоненты в г-реакции [1].
Теплоемкость Ср и энтальпия Н рассчитываются при помощи аппроксимации их зависимостей от температуры.
От использования тепловых эффектов реакций в явном виде при расчете температуры лучше отказаться вследствие их высокой температурной зависимости. В данной модели тепловой эффект учитывается в соответствии с законом Гесса, что делает модель более корректной.
При решении поставленной задачи необходимо решить систему дифференциальных уравнений, состоящую из одного уравнения (1) и т уравнений (2), где т - число реагентов. Стоит отметить, что, несмотря на кажущуюся простоту, данная система является весьма сложной для численного решения, поскольку в ходе процесса, описываемого данными уравнениями, скорости роста величин концентрации и температуры очень сильно изменяются, что делает систему не решаемой простыми явными методами (коэффициент жесткости достигает 1012).
Проблемы неустойчивости решения можно избежать, если использовать так называемые неявные методы решения. Другая возможность избежать этой проблемы - устранение быстрых стадий при помощи разделения временных масштабов [4].
Для решения данной системы уравнений разработана программа в среде Ма1ЬаЬ, которая использует встроенный солвер оёе15в, реализующий неявный метод Гира переменного порядка точности.
При расчетах использовались разные кинетические механизмы. Для простой модели весь процесс представлен одностадийной брутто-реакцией окисления молекулярного водорода:
Кинетические данные реакции: ¿о = 2,14 ■ 10141; Е = 129000 Дж/моль [5].
Начальные концентрации Н2 - 0,6667, О2 - 0,3333 (в мольных долях) при начальном давлении Р = 1 атм.
Для более точной модели применяются более сложные кинетики (табл. 1).
т
г
т
(1)
г
т
(2)
г= 1
3=1
2Н2 + О2 = 2Н2О.
В данном механизме реакции участвуют 6 компонентов:
н2, 02, Н20, н, о, он.
Таблица 1
Короткая кинетика окисления водорода [6]
Схема реакции К+ (моль,см,с) К- (моль,см,с)
0 е1 ^ І ^ ^ НН 1 ^ ^ _|_ о к к ^ + +м ^ ++++им+и ммк0+ог,+ ООп 11 м II ? 0 11 11 к о ^ ^ 4 + 5 + + 7 + + с к с % к с дококкоо 3ж1014ежр(-8, 81/Г) 3ж1014ежр(-4, 03/Т) 3ж1014ежр(-3, 02/Т) 3ж1014ежр(-3, 02/Т) 1, 85ж1017ежр(-54/Т)/Т 9, 66ж1018ежр(-62, 2/Т)/Т 8ж1016ежр(-52, 2/Т)/Т 5, 8ж1016ежр(-60, 6/Т)/Т 2, 48ж1013ежр(-0, 66/Т) 1, 30ж1013ежр(-2, 49/Т) 1, 33ж1015ежр(-10, 95/Т) 3,12ж1015ежр(-12, 51/^) 1016 1017 1016 6х1014
Здесь Т - абсолютная температура, деленная на 1000.
Наиболее полная кинекика содержит 44 реакции - 22 прямых и 22 обратных (табл. 2).
Таблица 2
Полная кинетика окисления водорода [7]
№ Реакция Прямая реакция Обратная реакция
ко Ь Е/Ко ко Ь Ко О
1,2 Н + 02 = ОН + 0 2,2є7 0 2590 9,3є7 0 10250
3,4 0 + Н2 = ОН + Н 1,8є4 1 4480 8,3є3 1 3500
5,6 Н + Н02 = 20Н 2,5є8 0 950 1,2є7 0 20200
7,8 2Н + М = Н2 + М 6,4є5 1 0 2,2є8 0 48300
9,10 Н + 0Н + М = Н20 + М 1, 4є11 2 0 2,2є10 0 52900
11,12 Н + Н202 = Н02 + Н2 6,3є6 0 550 6,8є7 0 9240
13,14 20Н + М = Н202 + М 9,1є2 0 -2550 1,2є11 0 22900
15,16 Н + Н202 = Н20 + 0Н 1,7є6 0 1900 7,3є5 0 9400
17,18 Н02 + 0Н = Н2 0 + 02 1,0є7 0 0 1,35є8 0 36480
19,20 Н202 + 0 = Н02 + 0Н 2,8є7 0 3200 6,25є6 0 9990
21,22 Н02 + Н02 = Н202 + 02 8,5є6 0 500 5,14є7 0 21670
23,24 Н + 02 = 0Н + 0 2,2є8 0 8450 1,3є7 0 0
25,26 Н + 02 + М = Н02 + М 5,0є3 0 -500 7,0є9 0 23000
27,28 20Н = Н20 + 0 6,3є6 0 550 6,8є7 0 9240
29,30 20 + М = 02 + М 1,1є3 0 0 4,51є11 0,35 59960
31,32 Н + Н02 = Н2 + 02 2,5є7 0 350 5.5є7 0 29100
33,34 Н02 + Н2 = Н20 + 0Н 1,2є5 0 12080 1,25є3 0,44 38820
35,36 0Н + Н202 = Н20 + Н02 1,0є7 0 910 2,8є7 0 16500
37,38 20Н = Н2 + 02 6,17є5 0 12830 1,91є8 0 21890
39,40 Н + Н02 = Н20 + 0 1,0є7 0 0 8,69є6 0 28330
41,42 0Н + М = 0 + Н + М 1,57є10 0 51460 1,0є4 0 0
43,44 Н02 + 0 = 0Н + 02 1,0є7 0 0 1,33є6 0,27 27600
В данном механизме реакции участвуют 9 компонентов:
Н2, 02, N2, Н2О, Н, О, ОН, НО2, Н2О2.
Анализ результатов
Проведение расчетов связано с выявлением соотношений точности и скорости счета для различных кинетик. Так для реакции окисления водорода был получен ряд результатов, отображенный на рисунке.
0.1
О 01
1 ю'3
1 ю'4
1 ю"5
5104 6 -10 7 10 8 -10 9 10 4 0.001 0.0011 0.0012
1.Т
--- Высокотемпературная асимптота [8)
--- Низкотемпературная асимптота [8]
+ + 3 ксперименгальнье даи-ье [9]
О О Экслерименгальньеда1-№е[10]
• • • • • В естбрук, кинетика воспламенешя метана
---Б абушэк. Короткая №иеп«акжеп«а
---Б абушок, Полная кинетика
--- Брутто-реакция
Зависимость адиабатического периода индукции от температуры
Как видно из графика, экспериментальные данные распределены вдоль двух асимптот, рассмотренных в [8]. Наиболее точное описание экспериментальных данных удалось получить, используя полную кинетику [7]. Использование сокращенной кинетики [6] допустимо лишь для описания высокотемпературной асимптоты, а упрощение механизма до одной брутто-реакции пригодно лишь для приблизительных оценок, хотя во многих случаях достаточно и этого, особенно если принимать во внимание скорость счета при использовании такой кинетики.
Работа выполнена при поддержке РФФИ. Грант № 13-01-00072.
Литература
1. Bode Shrestha, S.O. Hydrogen as an Additive to Methane for Spark Ignition Engine Applications / S.O. Bode Shrestha, G.A Karim // Intern. J. of Hydrogen Energy. - 1994. -V. 24, № 5. - P. 521-586.
2. Гельфанд, Б.Е. Водород: параметры горения и взрыва / Б.Е. Гельфанд, О.Е. Попов. -М.: Физматлит, 2008.
3. Рябинин, В.К. Математическое моделирование адиабатического периода индукции для метан-кислородных смесей в широком диапазоне начальных давлений и температур / В.К. Рябинин, Ю.М. Ковалев // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2013. - Т. 6, № 1. - С. 56-71.
4. Варнатц, Ю. Горение. Физические и химические аспекты, моделирование, эксперименты, образование загрязняющих веществ / Ю. Варнатц, У. Маас, Р. Диббл. - М.: Физматлит, 2003.
5. Бабушок, В.И. Глобальные кинетические параметры для высокотемпературных газофазных реакций / В.И. Бабушок, А.Н. Дакданча // ФГВ. - 1993. - № 4. - С. 48-80.
6. Бабушок, В.И. Тестовые примеры моделирования кинетики сложных реакций / В.И. Бабушок, А.Н. Дакданча. - Красноярск, 1990.
7. Максимов, Э.И. Физическое моделирование химических процессов в реакторе идеального перемешивания / Э.И. Максимов, А.Г. Мержанов // ДАН СССР. - 1975. - Т. 6, № 2. - С. 402-405.
8. Димитров, В.И. Простая кинетика / В.И. Димитров. - Новосибирск: Наука, 1982.
9. Mullins, B.P. Studies of the Spontaneous Ignition of Fuels Injected into a Bot Air Stream / B.P. Mullins. - NATO AGARD AG S/P2. 1952.
10. Patch, R.W. Shock-Tube Measurement of Dissociation on Rates of H2-J chem / R.W. Patch // J. Chem. Phys. - 1962. - V. 36, №7.-6 p.
Егор Евгеньевич Пигасов, аспирант, кафедра «Вычислительная механика сплошных сред>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Валерий Константинович Рябинин, кандидат физико-математических наук, доцент, кафедра «Вычислительная механика сплошных сред>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Юрий Михайлович Ковалев, доктор физико-математических наук, профессор, кафедра «Вычислительная механика сплошных сред>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Bulletin of the South Ural State University. Series «Mathematical Modelling, Programming & Computer Software>,
2013, vol. 6, no. 3, pp. 130-135.
MSC 80A25, 80A30, 97M10
Mathematical Modelling of the Adiabatic Thermal Explosion for the Hydrogen Oxidation Reaction
E.E. Pigasov, South Ural State University, Chelyabinsk, Russian Federation, [email protected],
V.K. Ryabinin, South Ural State University, Chelyabinsk, Russian Federation, [email protected],
Y.M. Kovalev, South Ural State University, Chelyabinsk, Russian Federation, [email protected]
This paper presents the results of mathematical modelling of the adiabatic reactor filled with a gas mixture of hydrogen and air. Despite the fact that the reaction of oxidation (combustion) of hydrogen is widely used for practical purposes, and there are a lot of publications, the development of kinetic schemes of hydrogen combustion is relevant at the present time. Various kinetic schemes of hydrogen combustion were tested and we have come to the conclusion that the most accurate description of the experimental data can be obtained using the complete oxidation kinetics. It is shown that the use of the reduced oxidation kinetics is relevant to describe the high-temperature asymptote, while the simplification of the mechanism to a brutto-reaction is suitable only for rough estimates.
Keywords: hydrogen, chemical kinetics, combustion and explosion, the induction period.
References
1. Bode Shrestha S.O., Karim G.A. Hydrogen as an Additive to Methane for Spark Ignition Engine Applications. Intern. J. of Hydrogen Energy, 1994, vol. 24, № 5, pp. 521-586.
2. Gelfand B.E., Popov O.E. Vodorod: Parametri Goreniya i Vzriva [Hydrogen: the Parameters of Combustion and Explosion]. Moscow, FIZMATLIT, 2008.
3. Ryabinin V.K., Kovalev Y.M. Mathematical Modeling Adiabatic Induction Period for Oxygen-Methane Mixtures in a Wide Range of Initial Temperatures and Pressures [Matematicheskoye modelirovaniye adiabaticheskogo perioda induktsii dlya metan-kiclorodnih smesey v shirokom diapazone nachalnih davleniy i temperatur]. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming & Computer software, 2013, vol. 6, no. 1, pp. 56-71.
4. Varnats Y., Mass U., Dibbl R. Goreniye. Fizicheskiye i himicheskiye aspekti, modelirovaniye, eksperimenty, obrazovanie zagryaznyayuschih veschestv [Combustion. Physical and Chemical Aspects, Modeling, Experimentation, Creation of Pollutants]. Moscow, FIZMATLIT, 2003.
5. Babushok V.I., Dakdancha A.N. Global Kinetic Parameters for High-Temperature GasPhase Reactions [Globalniye kineticheskiye parametry dlya visokotemperaturnih gazofaznih reakciy]. Combustion, Explosion, and Shock Waves, 1993, no. 4, pp. 48-80.
6. Babushok V.I., Dakdancha A.N. Testoviya primery modelirovaniya kinetiki sloznih reakciy [Test Examples of Modeling the Kinetics of Complex Reactions]. Krasnoyarsk, 1990.
7. Maksimov E.I., Merzanov A.G. Physical Simulation of Chemical Processes in the Reactor of Ideal Mixing [Fizicheskoe modelirovaniye himicheskih processov v reaktore idealnogo peremeshivaniya]. DAS USSR, 1975, vol. 6, no. 2, pp. 402-405.
8. Dimitrov V.I. Prostaya kinetika [Simple Kinetics]. Novosibirsk, 1982.
9. Mullins B.P. Studies of the Spontaneous Ignition of Fuels Injected into a Bot Air Stream. NATO AGARD AG S/P2, 1952.
10. Patch R.W. Shock Tube Measurement of Dissociation on Rates of H2-J chem. J. Chem. Phys., 1962, vol. 36, no. 7, 6 p.
Поступила в редакцию 30 мая 2013 г.
