CC BY
86
УДК 621.18.01
РАЗРАБОТКА ТЕОРЕТИЧЕСКИХ ОСНОВ СНИЖЕНИЯ ОБРАЗОВАНИЯ ОКСИДОВ АЗОТА В ТОПКАХ ГАЗОТРУБНЫХ КОТЛОВ
UDC 621.18.01
DEVELOPMENT OF THE THEORETICAL BASIS OF REDUCED FORMATION OF NITROGEN OXIDES IN THE FURNACE OF BOILERS
Михайлов Андрей Гаррьевич к.т.н., доцент
Батраков Петр Андреевич аспирант
Теребилов Сергей Викторович аспирант
Омский государственный технический университет, Омск, Россия
Рассмотрены вопросы образования оксидов азота при горении органического топлива в газотрубном котле, изложенные в отечественной и зарубежной литературе за последние годы. Предлагаются химические реакции и математические модели, сопровождающие появление термического, быстрого и топливного N0,,. Приведена оценка влияния формы топки газотрубного котла на состав продуктов сгорания газообразного топлива
Ключевые слова: ГОРЕНИЕ, ОКСИДЫ АЗОТА, ТЕРМИЧЕСКИЙ, БЫСТРЫЙ, ТОПЛИВНЫЙ
Mikhailov Andrey Garrievich Cand.Tech.Sci., associate professor
Batrakov Petr Andreevich postgraduate student
Terebilov Sergey Victorovich postgraduate student
Omsk State Technical University, Omsk, Russia
The problems of formation of nitrogen oxides during combustion of organic fuels in the boiler set out in the national and international literature in recent years. Chemical reactions and mathematical models that accompany the emergence of thermal, prompt and fuel NOx are discussed in this article. Impact assessment forms furnace boiler composition of the products of combustion of gaseous fuels cited by the authors
Keywords: COMBUSTION, OXIDES OF NITROGEN, THERMAL, FAST, FUEL
В России все более широкое применение находит автономное (децентрализованное) теплоснабжение. Оно особенно развито в тех регионах России, где нет недостатка в сетевом (природном) газе и имеется достаточно развитая система газоснабжения. Эти системы обеспечивают более высокий уровень теплового комфорта и создают дополнительные возможности для энергосбережения.
Одной из важнейших частей автономного теплоснабжения является газотрубные котлы малой мощности, которые нужны для получения пара и горячей воды. Интенсифицировать процессы переноса энергии и улучшения характеристик котла возможны с помощью изменения геометрических характеристик топки котла.
Моделирование процессов, происходящих в топке газотрубного котла, осуществляется с помощью расчетного комплекса АЫБУБ [1,2].
Для описания турбулентных течений реагирующих газов используется модель турбулентности с двумя уравнениями. Эта модель турбулентности получила название
к-б с двумя уравнениями (к - турбулентная кинетическая энергия, е - величина рассеивания кинетической энергии).
Рассмотрим основные уравнения, которые описывают реагирующую газовую смесь при следующих основных допущениях [3]: газовая смесь, заполняющая топочный объем - серое тело; теплота от факела к стенке в основном переносится излучением, конвекцией и молекулярной теплопроводностью; внутри пограничного слоя давление не изменяется вдоль нормали к поверхности тела и равно соответственному давлению на внешней границе пограничного слоя; внутри температурного пограничного слоя члены, характеризующие изменение энергии вследствие конвекции и изменения времени одного порядка с членами, характеризующими изменение энергии вследствие молекулярной теплопроводности; суммарный перенос теплоты на границе раздела газовой смеси - стенка осуществляется за счет адвекции и излучения; реагирующий газ СН4 - 100 %, окислитель - воздух.
1. Неразрывности для всей смеси:
= = о (1)
где р - плотность газовой смеси, II - вектор скорость, / - время.
2. Неразрывности для каждого компонента:
(2)
Sj - скорость образования /-компонента, - концентрация вещества
/-компонента, Р/ - плотность каждого /-компонента, - коэффициент диффузии, // -коэффициент диффузии для /-компонента, /лг - турбулентная составляющая
V V
Зсг = — Бс*. = — г г
динамическои вязкости, вязкость.
3. Моментов:
дри
- число Шмидта, V — кинематическая
дг
+ ч-(р11 (8> У)- V. = -*Р’ -ь V * (їріеГГУи}іт + в
(3)
где В - сумма всех сил, действующих на объем газа, цед- - эффективная турбулентная вязкость, Р - давление.
В формулах (1) - (3) используются следующие обозначения:
и =
Щ
и,
дих дЩ ди
' Л + д + д и -
дх ду дг
иу
и7
дих дЩ ди2
* а + л + я
дх о у о г
V • (р[/ ® 17) =
дх
д
V * (рі7 ® 17) =
дх
—(р их и у ) + щ (р и у и у ) + ^ (р и2 иу)
£ ^(риуиг), §^(р%и2) §-к(риА)+щ(риуих)+^(ри2их)
д^(Рихиу)+ду,(риуиу)+д^(Ригиу)
А(р ихиг) і %_[рпу и*)*
4. Энергии и диссипации:
+ ч >{рик) = V • +^?1^ + Рй-/э£ -5я
+ v. (pUk) = v. 4- £0vfc] + Рк - рг - SR
aM„.fpUQ &)«] +1;^_ - С£2р£)
°т^(риф^. [(м+gi)v£j♦^
где SR - источниковый член, Csl,Cs2, с а- , с е - справочные константы [1, 2\, Рк -параметр турбулентности, характеризует соотношение между силами вязкости и силами выталкивающими [1,2]:
Рк = №v. (vV + ?[/?)- | v. l/(3rtv. [/ + pfc) + Рк„
Pk = pitvU.(vU + vVT)-|v • I/(3ptv .U+рк) + Р№ (6)
5. Определение энтальпии:
Общая удельная энтальпия htot определяется следующим выражением:
htot=h + k, (7)
где И - удельная энтальпия неподвижной газовой смеси.
6. Определения вязкости:
к-8 модель основывается на концепции турбулентной вязкости, поэтому:
Meff=M + Mt, (8)
где // - динамическая вязкость. В данной модели предполагается, что турбулентная вязкость связана с турбулентной кинетической энергией и диссипацией через выражение:
(9)
где С ft— справочная константа [1,2].
Переменные к и е являются результатом решения дифференциальных транспортных уравнений для турбулентных кинетической энергии и диссипации:
7. Уравнение состояния:
Уравнение состояния, предложенное Редпихом-Квонгом, выглядит следующим образом [1,2]:
ЯТ а СП И.Т а(П
Р =--------_--------Р =--------_--------
у-Ь + с у(у + Ь) у-Ь + с гКи + ЬЭ, (10)
где V - удельный объем. Величины а, Ь, с - константы, зависящие от конкретного
вещества.
8. Начальные условия:
Принимаются значения всех параметров, входящих в систему уравнений, при времени ( = 0 и при начальной температуре Т = 300 К.
9. Граничные условия:
Для математического моделирования конкретных течений многокомпонентного реагирующего газа необходимо поставить соответствующие граничные условия [1,2]:
- условия прилипания на непроницаемой поверхности 1¥
Чг"0®^0. (11)
- задаются условия сложного теплообмена на стенке
Ч = Чс + Чк, (12)
где д - плотность теплового потока через пограничный слой от реагирующей газовой смеси к стенке, дс - плотность теплового потока, переносимого объемом газа, со стороны реагирующей газовой смеси к стенке, д^ - плотность радиационного теплового потока со стороны реагирующей газовой смеси к стенке (определяется методом Монте-Карло).
- непроницаемость;
- градиент давления в направлении нормали к поверхности.
На проницаемой границе, которой является вход в топке, задаются:
- массовая скорость;
- концентрации газов;
- интенсивность турбулентности, которая является приближенным значением для внутреннего течения в трубе;
На проницаемой границе, которой является выход в топке, задается:
- давление Р, которое является характеристикой системы для удаления
продуктов сгорания.
Уравнения (1 - 12), входящие в к-є модель, при соответствующих граничных и начальных условиях в трехмерной постановке численно решаются с использованием программы ANSYS CFX [1,2].
Расчетная область представляет собой горизонтальный цилиндр длиной 1,0 м, с площадью поперечного сечения 0,66 м2 (рисунок 1) при массовом расходе топливо-воздушной смеси на входе 0,061 кг/с.
Total Temperature
(temp)
- 1.6Ме+003
- 1. 156Є+003
- 7.8326+002
[К]
0-500 (m)
Рисунок 1. Распределение температур в топке
Основными результатами расчетов явилось определение значений температур, скоростей и состава продуктов сгорания. В том числе и концентрации соединений азота с кислородом [4].
При сжигании органических топлив в топках котлов азот, содержащийся в воздухе, взаимодействуя с кислородом, образует оксиды:
шх=ш + ш2+м2о,
Основная доля образовавшихся в продуктах сгорания паровых и водогрейных котлов ЫОх (95...99 %) приходится на монооксид (оксид) азота N0. Диоксид N02 и гемиоксид N20 азота образуются в значительно меньших количествах.
Образование монооксида (оксида) азота при сжигании органическго топлива СН4 происходит за счет окисления азота воздуха N2. В настоящее время известны три
механизма, по которым происходит образование оксидов азота: термический, быстрый и топливный. При образовании термических и быстрых N0 - источником азота является воздух, а в случае образования топливных N0 азотсодержащие составляющие топлива.
Механизм образования термических оксидов азота при соответствующих граничных и начальных условиях в трехмерной постановке решаются с использованием программы СБХ [2]. Он был предложен Я.Б. Зельдовичем [5,
6, 7, 8] и включает следующие реакции:
О + N2 -► N0 Ш,
N + 02 -► N0 + О,
Позднее он был дополнен реакцией атомарного азота с гидроксилом и получил название расширенного механизма Я.Б. Зельдовича:
ОН + N —>N0 +Н,
Реакции образования термических N0 характеризуются высокой энергией активации, поэтому образование оксидов азота происходит в области высоких температур, превышающих 1800 К. Концентрация термических N0 интенсивно возрастает от начала зоны горения и достигает наибольших значений непосредственно за зоной максимальных температур. Далее по длине факела концентрация оксидов азота практически не изменяется. Выражения для констант скоростей к каждой из трех реакций представлены в работе [2] и выглядят следующим образом:
к7 = (6,4 ■ 1о9)ехР^- ка = (6,4 ■ 1о9)ехР^-
къ=3,0-1013.
Образование термического N0 определяет следующие основные факторы: температура в зоне горения, коэффициент избытка воздуха и время пребывания продуктов сгорания в зоне высоких температур. Количество этого компонента
SnO,thermal Определяется Выражением:
*^N0 ,therma
где кГ},егта1=2к\, Ж^о - молярная масса термического N0; [О], [N2] - мольные концентрации кислорода и азота.
Исследования по сжиганию углеводородных топлив, проведенные Фенимором [2, 6, 7, 8], показали, что во фронте пламени за весьма короткий промежуток времени происходит образование оксидов азота по механизму, отличному от предложенного Я.Б. Зельдовичем. Обнаруженный оксид азота был назван быстрым из-за достаточно большой скорости его образования в корневой части факела. При этом вблизи зоны горения наблюдались значительные количества цианида водорода НОЧ, что объясняется реагированием молекулярного азота с углеводородными радикалами:
Реакции образования быстрых оксидов N0 протекают достаточно энергично при температурах 1200... 1600 К, когда образование термических оксидов азота практически не происходит.
Количество этого компонента 5к0_ рГотрт определяется выражением согласно [2]:
где Ж^о обозначает молярную массу N0, АрготрГ - число Аррениуса. Приведенные формулы для вычисления скоростей реакций применимы только к течениям с небольшими значениями Яе (критерий Рейнольдса). В турбулентных системах колебания параметров могут оказать доминирующее влияние на скорость формирования N0. В этом случае для расчетов используются статистические методы
В качестве объектов для расчетных исследований авторами рассмотрены топки
СН + N2 -► HCN +N, HCN + 02^N0 + ...,
/
= prompt 10 2] j
[2].
газотрубного котла с профилями окружность и эллипс. Расчетные модели имеют следующие геометрические характеристики (Рисунок 2).
Установлено, что с изменением формы профиля с окружности на эллипс наблюдается рост конвективных явлений. Это сопровождается снижением относительной концентрации оксидов азота ЫОэл / Ы0окр (N0^ - массовая концентрация оксидов азота на выходе из топки с формой поперечного сечения эллипс, Шокр - массовая концентрация оксидов азота на выходе из топки с формой поперечного сечения окружность) на выходе из топки газотрубного котла (Рисунок 3), уменьшением максимальной и средней температуры Т (Рисунок 4) газовой смеси внутри топки и процессами рециркуляции в соответствии с уравнением:
где v - стехиометрический коэффициент; Fuel - разновидность топлива.
Все выше перечисленные явления приводят к уменьшению относительной концентрации NOOT/ N0OKp на выходе из топки до 22 % в диапазоне а/Ь от 1,2 до 1,4.
а)
б)
Рисунок 2. Геометрические характеристики расчетных моделей: а - окружность: а - диаметр окружности; б - эллипс вертикальный: а - большая полуось, Ъ - малая полуось
1
NO + vFFuel --N2 +vCOzC02-bVllj0H1o
(И)
Рисунок 3. Зависимость относительной концентрации ЫОэл / Ы0ОКр от отношения большей полуоси к малой полуоси эллипса
Рисунок 4. Зависимость температуры газовой смеси осредненной по объему от отношения большей полуоси к малой полуоси эллипса
Рисунок 5. Зависимость температуры газовой смеси осредненной по объему от
коэффициента оребрения Аналогичные эффекты - уменьшение средней температуры в газовом объеме (Рисунок 5), явление рециркуляции достигаются с использованием цилиндрической топки с внутренним поперечным оребрение (Рисунок 6).
Рисунок 6. Расчетная схема цилиндрической топки с внутренним поперечным оребрением
На рисунке 7 приведена зависимость относительной концентрации ЫОор / Ы0окр от коэффициента оребрения 8, где ЫОор - массовая концентрация оксида азота на выходе из цилиндрической топки с оребренной внутренней поверхностью топки, 8 -отношение площади поверхности оребренной топки к неоребренной. Откуда следует, что данная конструкция топки позволяет снизить концентрацию оксидов азота до 50 % в диапазоне значений коэффициента оребрения от 1,05 до 1,5.
Рисунок 7. Зависимость относительной концентрации Ж)ор/Ж)окр
от коэффициента оребрения
Выводы
Показано, что процессы горения в топках газотрубных котлов всегда сопровождаются движением газов - воздуха, газообразного топлива, продуктов сгорания и являются совокупностью взаимообусловленных аэродинамических, тепловых и химических процессов. Так же при сжигании органического топлива происходит образование монооксида азота (оксида) за счет окисления азота воздуха.
Изменение формы поперечного сечения топочной камеры с окружности на
эллипс, а также оребрение внутренней поверхности цилиндрической топки приводит к:
1. Интенсификации конвективных явлений.
2. Увеличению рециркуляции продуктов сгорания.
3. Уменьшению концентраций оксидов азота на выходе из топки.
Литература
1. Михайлов А.Г., Батраков П.А., Теребилов С.В. Численное моделирование процессов тепломассопереноса при горении газообразного топлива в топочном объеме // Естественные и технические науки. - 2011.- №5(55).- С. 354-358.
2. ANSYS CFX-Solver Theory Guide. ANSYS CFX Release 11.0 / ANSYS, Inc. // Southpointe 275 Technology Drive. - Canonsburg : PA 15317, 2006. - 312 p.
3. Алексеев Б. В. Физическая газодинамика реагирующих сред. - М.: Высшая школа, 1985. -
464 с.
4. Михайлов А.Г. Вопросы образования оксидов азота при сжигании газообразных и жидких топлив // Омский научный вестник. - 2009. - № 3 (83). - С. 103 - 106.
5. Пашков Л. Т. Основы теории горения. - М. : МЭИ, 2002. - 136 с.
6. Росляков П.В., Закиров И.А. Нестехиометрическое сжигание природного газа и мазута на тепловых электростанциях. - М. : МЭИ, 2001. - 144 с.
7. Fenimore С.Р. Formation of nitric oxide in premixed hydrocarbon flame // 13-th Sympos. (Interm.) On Combustion. 1971. P. 373.
8. Hayhurts A.N., Vince I.M. Nitric oxide formation from N2 in flames: The importance of «prompt» NO // In Progress in Energy and Combustion Science. Oxford, 1980. Vol. 6. P. 35-51.
References
1. Mikhailov A.G., Batrakov P. A., Terebilov S.V. Numerical modeling of heat and mass transfer in combustion of fuel in the combustion volume // Natural and Technical Sciences. - 2011. - № 5 (55). - P. 354-358.
2. ANSYS CFX-Solver Theory Guide. ANSYS CFX Release 11.0 / ANSYS, Inc. // Southpointe 275 Technology Drive. - Canonsburg: PA 15317, 2006. - 312 p.
3. Alexeev B.V. Physical Gas Dynamics reacting media. - М.: Higher School, 1985. - 464 p.
4. Mikhailov A.G. Questions of nitrogen oxide combustion of gaseous and liquid fuels // Omsk Scientific Bulletin. - 2009. - № 3 (83). - P. 103 - 106.
5. Pashkov L.T. The fundamental theory of combustion. - М.: MEI, 2002. - 136 p.
6. Roslyakov P.V., Zakirov I.A. Non-stoichiometric combustion of natural gas and fuel oil for thermal power plants. - М.: MEI, 2001. - 144 p.
7. Fenimore C.P. Formation of nitric oxide in premixed hydrocarbon flame // 13-th Sympos. (Interm.) On Combustion. 1971. P. 373.
8. Hayhurts A.N., Vince I.M. Nitric oxide formation from N2 in flames: The importance of «prompt» NO // In Progress in Energy and Combustion Science. Oxford, 1980. Vol. 6. P. 35-51.
