CC BY
199
УДК 621.181:519.876
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ СЖИГАНИЯ РЕЗЕРВНОГО ТОПЛИВА В ТОПКЕ КОТЛА БКЗ-210-140
Гиль Андрей Владимирович,
канд. техн. наук, доцент каф. парогенераторостроения и парогенераторных установок Энергетического института ФГАОУ ВО «Национальный исследовательский Томский политехнический университет», Россия, 634050, г. Томск, пр. Ленина, д. 30. E-mail: andgil@tpu.ru
Заворин Александр Сергеевич,
д-р техн. наук, профессор каф. парогенераторостроения и парогенераторных установок Энергетического института ФГАОУ ВО «Национальный исследовательский Томский политехнический университет», Россия, 634050, г. Томск, пр. Ленина, д. 30. E-mail: zavorin@tpu.ru
Лебедь Дмитрий Викторович,
директор филиала «Приморская генерация» ОАО «Дальневосточная генерирующая компания», Россия, 690034, г. Владивосток, ул. Фадеева, д. 47А. E-mail: lebed-dv@dvgk.rao-esv.ru
Старченко Александр Васильевич,
д-р физ.-мат. наук, профессор каф. вычислительной математики и компьютерного моделирования механико-математического факультета ФГБОУ ВПО «Национальный исследовательский Томский государственный университет», Россия, 634050, г. Томск, пр. Ленина, д. 36.
E-mail: starch@math.tsu.ru
Актуальность работы обусловлена необходимостью оценки эффективности и надежности работы котельных агрегатов на резервном топливе после реконструкции связанной с изменением базового топлива и применяемых подходов численного анализа. Цель работы: анализ физико-химических процессов в топочном объеме котла БКЗ-210-140 при работе на резервном топливе и номинальной нагрузке при использовании Эйлерова и совмещенного Эйлерова и Лагранжева подходов к моделированию. Построение математических моделей объекта исследования с учетом рационального количества расчетных областей и применения современных численных алгоритмов и подходов. Сравнение численных результатов, полученных с использованием различных программных комплексов и различного представления горения жидкого топлива в камерной топке энергетического котла. Оценка температурных полей в объеме топочной камеры, концентраций химических элементов продуктов сгорания, формирования аэродинамической структуры потоков и эффективности выгорания резервного топлива.
Методы исследования: численное моделирование на основе Эйлерова и Лагранжева подходов, метод конечного объема, пристеночных функций, двухпараметрическая k-s модель турбулентности, модель протекания химической реакции в турбулентном потоке Eddy-Break-Up - EBU, Р1-приближение метода сферических гармоник, теоретический расчет на основе нормативного метода теплового расчета котлов.
Результаты: численное представление физико-химических процессов в топке котла БКЗ-210-140 при организации сжигания резервного топлива; сравнение результатов, полученных по различным численным моделям; оценка надежности работы испарительных поверхностей и эффективности эксплуатации котельного агрегата на резервном виде топлива; рекомендации по совершенствованию топочного режима с целью повышения эффективности эксплуатации котла БКЗ-210-140 при сжигании мазута.
Ключевые слова:
Топочная камера, горелочное устройство, параметр крутки, газодинамика, температура, мазут.
Введение
На тепловых электрических станциях, работающих как на твердом топливе, так и на природном газе, в качестве резервного топлива обычно используется мазут. Поскольку при эксплуатации котельных агрегатов запас резервного топлива должен обеспечивать бесперебойную работу в течение 3,5 суток, в ходе выполнения пуско-наладочных работ основное внимание уделяется работоспособности на основном топливе, а эффективность работы котельного агрегата на резервном топливе из-за краткосрочности рассматриваемого периода оценивается в меньшей степени.
Помимо того, что жидкие органические топлива являются самыми дорогими среди энергетических топлив, при нарушении режимов их сжигания кроме повышения расхода топлива могут наблюдаться и другие негативные последствия в виде снижения экономичности, надежности работы котла, повышения выбросов вредных веществ в атмосферу [1, 2].
Известно, что оптимальная для каждого вида топлива организация топочной аэродинамики обеспечивает устойчивость горения, позволяет поддерживать требуемые скорости и температуры
в топочном объеме, снизить уровень выбросов токсичных веществ в окружающую среду. Неравномерное и неустойчивое движение газов в топке, напротив, ведет к противоположным результатам и может наряду с этим вызвать превышение допустимых локальных тепловых нагрузок, что отразится на эффективности работы котельного агрегата [3, 4].
Реконструкция котельных агрегатов БКЗ-210-140, предназначенных по проекту для факельного сжигания чихезского бурого угля, проведена с целью перевода на сжигание природного газа в качестве базового топлива. Резервным топливом по-прежнему остался мазут, поэтому для оценки эффективности работы новой схемы сжигания на резервном топливе необходимо проанализировать физико-химические процессы в объеме топочной камеры.
Поскольку натурные многофакторные исследования процесса горения в топках котельных агрегатов, базирующиеся на экспериментальном изучении трехмерной аэродинамики горящего факела, в том числе с вихревыми течениями, являются очень дорогостоящими и трудозатратными, в настоящее время для исследования и анализа газодинамики и теплофизических процессов в котлах эффективным инструментом является численное моделирование с использованием как разработанных за рубежом (ANSYS CFD (Fluent & CFX), STAR-CD, FLOW3D, Open-FOAM и др.), так и отечественных (Flow Vision, VP2/3, SigmaFlow, FIRE 3D и т. п.), но в разной мере апробированных для этих целей пакетов прикладных программ. Обилие физических моделей в перечисленных пакетах позволяет с хорошей сходимостью предсказывать ламинарные и турбулентные течения, различные режимы теплопереноса, химические реакции, многофазные потоки и другие явления на основе гибкого построения сеток и их адаптации к полученному решению [5].
Процесс горения мазута представляется последовательным протеканием следующих фаз: мелкодисперсное распыление, испарение и термическое расщепление, концентрирование газовоздушной смеси, воспламенение и ее горение [6, 7]. Хотя при этом в факеле образуются и твердые соединения сажи и кокса, в большинстве работ исследуются процессы испарения, воспламенения и горения газовой фазы и не учитывается наличие твердой фазы, горение которой происходит аналогично гетерогенному горению частиц твердого топлива [3, 4]. Поэтому из анализа литературных источников сложно отдать предпочтение тому или иному численному подходу к решению задачи по горению жидкого топлива в факеле.
В связи с изложенным в данной работе были использованы ресурсы пакета прикладных программ FIRE 3D и универсального программного комплекса ANSYS. При этом расчет в программном комплексе ANSYS проведен с учетом только газовой фазы, а в пакете прикладных программ FIRE 3D -с учетом твердой фазы.
Объект моделирования
Объектом исследования является котельный агрегат БКЗ-210-140 - вертикально-водотрубный, од-нобарабанный, с естественной циркуляцией, П-об-разной компоновки, с твердым шлакоудалением. Номинальная паропроизводительность 210 т/ч, давление перегретого пара 14 МПа, расчетная температура перегретого пара 560 °С. Из-за состояния паропроводов высокого давления и необходимости продления паркового ресурса металла парового тракта котел переведен на пониженную температуру перегретого пара 545 °С. Топочная камера объемом 1093 м3до реконструкции была оборудована четырьмя блоками тангенциально установленных прямоточных пылеугольных горелок и четырьмя мельницами-вентиляторами.
В результате выполненной реконструкции топочное устройство котла оснащено шестью низкоэмиссионными вихревыми комбинированными газомазутными горелками ГМВИ (III) - 30 производства ОАО «ТКЗ - Красный котельщик» тепловой мощностью 30 МВт каждая с центральной и периферийной раздачей природного газа и двумя потоками воздуха. Горелки расположены на боковых стенах топки по схеме треугольника с вершиной вниз (рис. 1). Воздух в газомазутные горелки подается из воздушного тракта после воздухоподогревателя и распределяется на центральный и периферийный каналы посредством регулирующей арматуры.
В качестве расчетного топлива принят малосернистый мазут марки М100 [8].
Для подавления образования оксидов азота предусмотрена система двухступенчатого сжигания, суть которой состоит в том, что основная часть топлива сжигается при недостатке воздуха, а дожигание происходит за счет подачи дополнительного воздуха (около 20 %) в сопла третичного дутья в области более низких температур факела. Сопла сечением 608x208 мм установлены по два на боковых стенах топки выше верхнего яруса горелок на отметке 11,5 м.
Для сжигания мазута в центральную часть газомазутных горелок установлены паромеханические форсунки типа «Вулкан» (рис. 2). Перед форсунками обеспечивается давление топлива не менее 2,0 кгс/см2, давление пара не менее 3,0 кгс/см2, температура распыливающей среды (перегретый пар) 200...250 °С, удельный расход пара не более 0,02 кг/кг, вязкость топлива не более 6 °ВУ согласно ГОСТ 6258-85 [9].
Паромеханические форсунки довольно экономичны в работе в широком диапазоне регулирования, позволяющем достигать глубины регулирования до 10 % номинальной производительности без ухудшения распыления. Они имеют такую конструкцию, что при нагрузке выше 80 % работают как механические, при более низкой - как паровые форсунки. В паромеханическом режиме задействованы два канала - для подачи мазута и пара в топочное пространство. При работе котла на низ-
Рис. 1. Эскиз топочной камеры котла БКЗ-210-140: а) боковая стена; б) трехмерная модель
Fig. 1. Rough drawing of combustion chamber of boiler BKZ-210-140: a) side wall; b) three-dimensional model
кой нагрузке используется пар, который, пройдя паровой канал, систему отверстий в корпусе форсунки и паровой завихритель-кавитатор, встречается с мазутом, распыляя его.
106,5
89,5
////////yS///S////////
in 1Л У lliiiB^ гч X (N
to , л Жч^Щ^Ш 3 LD 5
W4
Г Г Ч\\ХХЧ\\у\\\>К\\ЧЧПьУ 5/ 4/ 1 / \ 3 \б \ 2 \J_
Рис. 2. Форсунка паромеханическая «Вулкан»: 1) корпус;
2) кавитатор мазутный; 3) сопло; 4) гайка поджимная;
5) излучатель паровой; 6) прокладка; 7) основание
Fig. 2. Steam-assisted burner «Vulkan»: 1) body; 2) fuel oil ca-
vitator; 3) nozzle; 4) press screw nut; 5) steam radiator;
6) shim piate; 7) base
Расчетная сетка и параметры моделирования
При работе в программном комплексе ANSYS построенная геометрическая модель топочной камеры посредством Ansys Workbench v. 12.1.4 преобразована в расчетную структурированную тетрагональную сетку (рис. 3, а) с размерами стороны ячеек в диапазоне от 1,40710-4 до 1,4072 10-1м, что позволяет производить численное моделирование с достаточно высокой производительностью и приемлемой точностью. Полученная сетка содержит 133777 контрольных объемов и 26336 расчетных узлов [10].
После импортирования в расчетный модуль универсального пакета прикладных программ Ansys Fluent v. 12.1.4 было произведено преобразование тетрагональной структуры сетки в полиэдральную (рис. 3, б). В результате оптимизации количество контрольных объемов снизилось до 44291, что благоприятно отразилось на производительности расчета и на количестве «плохих» ячеек. Количество расчетных узлов полиэдральной сетки увеличилось до 183649, что, в свою очередь, повысило сходимость и точность численного моделирования.
При работе с программным комплексом FIRE 3D с использованием препроцессора была построена декартовая сеточная модель топочной камеры котла БКЗ-210-140 с количеством ячеек 169200.
Граничные и начальные условия математической модели задавались в соответствии с результатами предварительно проведенного теплового расчета котла, в качестве исходных данных для которого принимались: топливо - пары мазута (С^Що) при расчете в программном комплексе Ansys Fluent v. 12.1.4 и малосернистый мазут (с элементарным составом (%) на рабочую массу: влажность 0,2, зольность 0,03, содержание серы 0,85, углерода 86,58, водорода 12,04, кислорода 0,3) при расчете с использованием пакета прикладных программ FIRE 3D; массовый расход топлива на номинальной нагрузке 4,057 кг/с; окислителем является атмосферный воздух, теоретически необходимый объем которого составляет 10,9059 м3/кг; коэффициент избытка воздуха на выходе из топки 1,2 (вследствие повышенных присосов воздуха); доля первичного воздуха по отношению к общему количеству 0,175, вторичного воздуха, подаваемого в вихревые горелки 0,525, вторичного воздуха через тангенциальные
Рис. 3. Структурированная расчетная сетка для математического моделирования: а) тетрагональная; б) полиэдральная Fig. 3. Structured computational mesh for mathematical simulation: a) tetragonal; b) polyhedral
горелки 0,10, третичного дутья через сбросные сопла 0,20; температура горячего воздуха 189 °С. Для повышения производительности расчетов воздух условно представлен в виде смеси кислорода с азотом в объемных долях 21 и 79 % соответственно.
Математическая модель
Методом конечных элементов моделировались следующие физико-химические процессы: горение паров мазута в кислороде воздуха, турбулентность в реагирующей среде, конвективный и лучистый теплообмен и т. д. В расчете принято упрощение процесса окисления горючего - уравнение брутто реакции горения паров мазута протекает необратимо в одну стадию: С19Н30+26,502=19С02+15Н20. Моделирование газовой фазы проводилось в приближении Эйлера. В соответствии с этим математическая модель включает в себя уравнения (запись для краткости приведена только для одной координаты) [11]: • неразрывности
д(рЦ) 3%:
состояния
p_-
RT
CH
MO mN Mc H MCO
O2 N2 C19H30 CO2
M
= 0;
движения
Здесь х1 - декартовы координаты; и1 - компоненты вектора скорости; р,р, Т, с, /, g - плотность, давление, температура, удельная теплоемкость, динамическая вязкость смеси, ускорение свободного падения; 8 -символ Кронекера; т^ - компоненты тензора рей-нольдсовых напряжений; й - число Шмидта; Я - коэффициент теплопроводности; С - концентрация компонентов; М] - молекулярный вес компонентов смеси.
Использовалась к-е модель турбулентности, которая хорошо зарекомендовала себя в задачах с сильной кривизной линий тока и справедлива для полностью развитого турбулентного течения, т. е. для высоких чисел Рейнольдса, когда прямое влияние вязкости на структуру турбулентности пренебрежимо мало [12, 13].
Для описания турбулентных характеристик газа используется двухпараметрическая к-е модель турбулентности Лаундера и Сполдинга
dpUUj dp
pu.
dxj
dx,
dk
dxj
■ = Tj
du, dx.
dx.
(dU, du,^ 2 dUk _
"UU"2
dx, dx, J 3 dxk j _ 1,2,3;
+ pgi,
ds
u,-
1 dx,
энергии
dUjpCpT _ dqL +dTvU
dx,
dxi dxj . dT
q< ;
dx,
+ S,
-- ps + Л 5k dxj
s du. s2 -1,92 — + k dxj k
1,3 j ds
dxj ,
k2 0,09p—. s
dx,
В данных уравнениях к, е - кинетическая энергия турбулентности и скорость ее диссипации; / -турбулентная вязкость.
В общем балансе энергии перенос тепла излучением в камерах сгорания играет важную роль. При численном исследовании радиационного теплопе-реноса в излучающих, поглощающих и рассеивающих средах, таких как топочные камеры, используется сложное интегро-дифференциальное уравнение для интенсивности термического излучения. Для приближенного решения этого уравнения в данной работе применяется метод сферических гармоник (его Р1-приближение). Использование такого подхода для моделирования термического излучения в топочной камере обусловлено его хорошей совместимостью с конечноразностны-ми методами.
Для математического моделирования протекания химической реакции в турбулентном потоке одной из часто используемых моделей является модель Сполдинга (Eddy-Break-Up - EBU), отраженная в работах [13-15]. Согласно данной модели скорость химической реакции пропорциональна интенсивности турбулентного перемешивания. Однако влияние температуры потока на скорость протекающих химических реакций EBU-модель не учитывает, что не позволяет моделировать течения с переходными режимами горения, определяющими условия существования пламени.
Модифицированная модель, основанная на совместном использовании модели Сполдинга и кинетической модели, расширяет возможности моделирования и позволяет учесть различные режимы течения.
Кинетическая скорость реакции определяется следующей формулой [16]
Wш = Ap2~cf~coz [-E /(Rf)],
где Ak - постоянный коэффициент, м3/(кг-с); E -энергия активации химической реакции, Дж/моль; R - универсальная газовая постоянная; cf - массовая концентрация горючего; cax - массовая концентрация окислителя; Т - температура газа.
Скорость турбулентного горения, согласно EBU-модели, вычисляется по формуле
w = aebuр(е/k)min{Cf,сох / v},
где е - скорость диссипации кинетической энергии турбулентности; k - удельная кинетическая энергия турбулентности; AEBU - коэффициент модели; v - стехиометрический коэффициент.
При использовании модели EBU совместно с методом крупных вихрей отношение вычисляется [17, 18] по зависимости:
W = aebubebuР I S I min{Cf , c0x / v}
е / k = у[2Щ,
где Бц - тензор скоростей деформации.
Тогда формула скорости турбулентного горения, соответственно, принимает вид
™, = Аьи|т1п{с/,сох / V}.
Для учета зависимости скорости реакции от температуры потока формула модифицируется следующим образом
ВЕви = 6,386 х 103ехр .
Вид связи ВЕВи и температуры выбран по аналогии с законом Аррениуса. Коэффициент ВЕВи изменяется в диапазоне от 0 до 1. Своё максимальное значение (ВЕВи=1) он принимает при Гши=2100 К (то есть при наибольшей температуре горения). Следовательно, при температуре, близкой к наибольшей, значения, полученные по двум последним формулам, становятся фактически одинаковыми. При уменьшении температуры коэффициент ВЕВи также начинает уменьшаться, обуславливая тем самым зависимость скорости реакции от температуры [19].
После вычисления кинетической и турбулентной скоростей горения скорость реакции окончательно определяется выражением, аналогичным формуле Щёлкина [20]
w = \ /W L + Vv2.
Использование этой формулы дает возможность моделировать переход как от турбулентного к ламинарному режиму горения при уменьшении интенсивности турбулентного перемешивания, так и, наоборот, от ламинарного режима к турбулентному, соответственно, при увеличении интенсивности.
Таким образом, предложенная модель является универсальной для расчета течений реагирующего газа в широком диапазоне чисел Рейнольдса [21].
В программном продукте FIRE 3D теплообмен и горение в газовой фазе представляются на основе Эйлерова способа описания, т. е. используются стационарные пространственные уравнения баланса массы, импульса, концентраций газовых компонентов и энергии для газовой смеси. Для описания движения и тепломассообмена одиночных частиц топлива и золы вдоль их траекторий с учетом обратного влияния дисперсной фазы на несущую среду применяется Лагранжев подход. Турбулентные характеристики газа рассчитываются с использованием двухпараметрической k-е модели турбулентности, также учитывающей влияние движущихся частиц. Радиационный теплообмен в двухфазном потоке представляется в рамках Р1-приближения метода сферических гармоник [5].
Соответственно основные уравнения запишутся в следующей форме:
• уравнение изменения массы газовой смеси вследствие испарения влаги, газификации и горения углеродных частиц
д(ри,)
дх,
■ = J + J + J, ,
evap vap char '
где р - плотность газовой смеси; и1 (¿=1,2,3) - компоненты скорости газа; х1 (¿=1,2,3) - декартовы координаты; Jevap, Jvap, Jchar - массовые скорости испарения влаги, топлива и догорания частиц;
уравнение количества движения для газовой фазы
dpUUj dp -j = -± + Pgj + F +
dx,
j
dx,
JdU. dU, ^
(ß+ß ^+dxrJ-
- 2 3 j
л dUk oxh
; j = 1,2,3,
где р - давление; лТ - турбулентная вязкость; л -молекулярная вязкость; к - энергия турбулентности; 81/ - символы Кронеккера (8Й=1; 8/=0 (Ц/)); gj -проекции ускорения свободного падения на декартовы оси; - компоненты силы межфазного взаимодействия между частицами и несущей средой; • уравнение движения частицы-маркера ¿-й группы (¿=1...Д),
dUi (U, - U'.)
—j = ^—J-+gj; j = 1,2,3.
dt t
Здесь
T =
4ppdp
3cD p| U' - U Г
24
-(1 + 0,15 Re 0'687); Re,. <1000;
Re
0,44;
Re,. > 1000;
Re =
X (Uj - Uj )2 dp,
j=1
ß
где и/ - компоненты скорости частиц; и/ - компоненты скорости газа; g¡ - проекции ускорения свободного падения на декартовы оси; рр - плотность частиц; "р1 - размер частиц 1-й группы; сВ1 - коэффициент сопротивления частиц; • уравнение энергии для газа
dpU.pT дх:
дх:
ß + ß Pr Pr
+ QvolJvol
+ QCO JCO
dq'r
дх,
dt) +
дх
- + 0,
где с - теплоемкость смеси газов; Рг, Рг( - молекулярное и турбулентное (Рг(=0,9) числа Прандтля; Т - температура газа; 2„„;, QCO - тепловые эффекты реакций горения углеводородов и угарного газа; д™" - компоненты вектора радиационного теплового потока; 0 - интенсивность теплообмена между несущей средой и частицами. • уравнение состояния
Р
Р = -
RT
Cn
C„
Co
Cc.
C
C
O
MO2 MN, Kol MCO MCO 2
MH O
где Ш} - молекулярные веса компонентов газовой смеси; Я0 - универсальная газовая постоянная, р -давление.
Результаты численного моделирования
На рис. 4 представлены результаты визуализации численных расчетов температурного уровня в горизонтальном сечении по оси верхнего яруса го-релочных устройств.
В обоих вариантах расчета наблюдается аналогичное расположение температурных полей с незначи-
a/a
-500 -450
-350 -300
0 1 2 3 4 5 6 7
б/b
Рис. 4. Температурные поля в горизонтальном сечении топочной камеры по оси верхнего яруса горелочных устройств при использовании вычислительных пакетов: а) ANSYS Fluent; б) FIRE 3D
Fig. 4. Temperature fields in horizontal section of a combustion chamber along the axis of upper layer of a burner device while using computational packages: a) ANSYS Fluent; b) FIRE 3D
D: =
I
1823.00 1737.35 1651.69 1566.04 1480.39 1394.73 1309.08 1223.42 1137.77 1052.11 966.46 880.80 795.15 709.50 623.84 538.19 452.53 366.88 281.22 195.57 109.91
i>1800 1750 У, м 1700
1650 24 1600
1550 22 ¡яц-1500 -1450 20 Ц-1400 Щ-1350 1-1300 1250 1200 1150 1100 1050 1000
а/а
бЬЬ
Рис. 5. Температурные поля в вертикальном сечении топочной камеры по продольной оси топочной камеры при использовании вычислительных пакетов: а) ANSYS Fluent; б) FIRE 3D
Fig. 5. Temperature fields in vertical section of a combustion chamber along the longitudinal axis of upper layer of a burner device while using computational packages: a) ANSYS Fluent; b) FIRE 3D
тельными отклонениями по температурным значениям, свидетельствующее об отсутствии единого, симметричного относительно центральной оси топочной камеры, ядра горения. Это может негативно отразиться на гидродинамических режимах работы экранных элементов в отдельных контурах циркуляции.
Визуальные картины распределения температурных полей в вертикальном сечении по продольной оси топочной камеры имеют отличия (рис. 5). При расчете в программе ANSYS Fluent максимальный уровень температур наблюдается в области расположения верхнего яруса горелочных устройств, а результирующая картина температурного уровня при расчете в программном комплексе FIRE 3D имеет максимум в области между ярусами горелок. При этом также можно отметить, что именно в данном сечении уровень температур до аэродинамического выступа при расчете в программном комплексе ANSYS Fluent выше (~100-200 °С), чем при расчете программой FIRE 3D.
Дальнейшие графические результаты представлены относительно пересечения продольной и поперечной оси топочной камеры, т. е. по центральной оси топки.
Анализируя изменение температуры по центральной оси топочной камеры, можно отметить в целом аналогичность профилей кривых (рис. 6, а). Но при этом значения температур на уровне середины холодной воронки и в области активного выгорания топливно-воздушной смеси существенно отличаются (~200 °С). Данное различие в большей степени связанно с локальными значениями в конкретных сечениях, поскольку применены различные численные модели расчета процессов горения, которые влияют на итерационный расчет аэродинамической структуры потока в топке, что вносит несоот-
ветствия при рассмотрении в конкретных точках. Температура на выходе из топки составляет 1093 °С при расчете с гомогенными реакциями горения и 992 °С при расчете гетерогенного процесса горения.
На рис. 6, б изображены два графика зависимости изменения концентрации кислорода по высоте топки, полученные из математического моделирования в Ansys Fluent и FIRE 3D. При сравнении кривых можно отметить аналогичность профилей, но при этом на уровне расположения горелочных устройств и третичного дутья при расчете гетерогенного горения мазута с использованием программного продукта FIRE 3D наблюдаются пики повышенной концентрации, из этого следует, что воспламенение топливовоздушной смеси при расчете в программе Ansys происходит быстрее и интенсивнее, что и соответствует гомогенному горению. В целом наблюдается более высокая концентрация кислорода при использовании программного комплекса FIRE 3D примерно на 2 %.
Скоростной режим восходящего газо-воздуш-ного потока имеет близкие значения в обоих исследованных вариантах за исключением зоны на уровне расположения нижнего яруса горелочных устройств (рис. 6, в). При расчете с использование программы FIRE 3D закрученные потоки, истекающие из горелочных устройств, сталкиваются, образуя зону с пониженным давлением, что приводит к падению движения топочной среды. Максимальные значения скоростей наблюдаются на уровне аэродинамического выступа и соответственно сужения топочной камеры, а также вследствие дополнительного ввода объема воздуха через сопла третичного дутья. В данном сечении скорость достигает 9 м/с при расчете в программе Ansys Fluent и 10,5 м/с при использовании FIRE 3D.
Температура, °С Концентрация кислорода, % Скорость, м/с
a/a б/b в/с
1) ANSYS Fluent; 2) FIRE 3D
Рис. 6. Результаты численного расчета топочной среды по высоте топочной камеры: а) изменение температуры; б) изменение концентрации кислорода; в) изменение скорости газо-воздушного потока
Fig. 6. Results of numerical computation of furnace atmosphere by a combustion chamber height: a) temperature change; b) oxygen concentration change; c) change in gas-air flow speed
Заключение
Какая из рассмотренных моделей более адекватна, заключить достаточно сложно, поскольку картины визуализации и профили кривых на представленных графиках по всем рассмотренным параметрам достаточно близки. Отличия в абсолютных значениях можно отнести к тому, что они снимались по центральной оси топки, и незначительное смещение температурных, концентрационных полей рассматриваемых параметров приводит к отличию результатов в одной и той же точке при расчете по разным моделям. При сравнении значений температур на выходе из топки при расчете в принятых расчетных комплексах и значения температуры, рассчитанного по нормативному методу теплового расчета, разница не превышает 100 °С, что является допустимым при расчетных исследованиях процессов в котельных агрегатах.
Анализ совокупности данных, полученных по двум принятым к расчету моделям, показывает, что работа котельного агрегата на резервном топливе представляется надежной, поскольку темпе-
ратурный уровень и его распределение в пределах топочного объема соответствуют допустимому тепловому напряжению и достаточно равномерен в сечениях по высоте топочной камеры. Соответственно, как процесс циркуляции пароводяной смеси, так и экраны циркуляционных контуров, судя по этим косвенным данным, в процессе эксплуатации имеют надежные параметры работы. Наличие высокоскоростных локальных зон не должно негативно отразиться на надежности работы поверхностей испарения за счет воздействия твердых частиц, поскольку наличие их в газовом потоке при сжигании мазута минимально.
Поскольку сжигание мазута по большей части является гомогенным горением в турбулентном потоке, то согласно рис. 6, б процесс выгорания топлива будет проходить полностью. Имеющий место после реконструкции котла повышенный коэффициент избытка воздуха, несомненно, при работе на мазуте влияет на снижение температуры на выходе из топки.
Работа выполнена в рамках ГЗ НИР № 2069 (2.13222014).
СПИСОК ЛИТЕРАТУРЫ
1. Ghorbani A., Bazooyar B. Optimization of the combustion of SOME (soybean oil methyl ester), B5, B10, B20 and petrodiesel in a semi industrial boiler // Energy. - 2012 August. - V. 44. -Iss. 1. - P. 217-227.
2. Daho T., Vaitilingom G., Sanogo O. Optimization of the combustion of blends of domestic fuel oil and cottonseed oil in a non-modified domestic boiler // Fuel. - 2009 July. - V. 88. - Iss. 7. -P. 1261-1268.
3. Белоусов В.Н., Смородин С.Н., Смирнов О.С. Топливо и теория горения. Ч. II. Теория горения. - СПб.: СПбГТУРП, 2011. -142 с.
4. Лебедев Б.В., Карякин С.К. Технология сжигания органических топлив. - Томск: Изд-во Томского политехнического университета, 2012. - 148 с.
5. Гиль А.В., Старченко А.В. Математическое моделирование физико-химических процессов сжигания углей в камерных топках котельных агрегатов на основе пакета прикладных программ FIRE 3D // Теплофизика и аэромеханика. - 2012. -Т. 19. - №5. - C. 655-671.
6. Аскарова А.С., Рыспаева М.Ж., Волошина И.Э. Численное моделирование горения и самовоспламенения двухфазных химических реагирующих течений с впрысками // Известия Томского политехнического университета. - 2009. - Т. 315. -№ 4. - С. 5-9.
7. González-Espinosa A., Lozano A., García J.A., Barreras F., Lincheta E. Design modification of the air diffuser in the burners of a fuel oil power plant. P. II: Interaction with the liquid spray // Fuel. - 2014 May. - V. 124. - P. 141-150.
8. Тепловой расчет котлов (нормативный метод). - СПб.: Изд-во НПО ЦКТИ, 1998. - 256 с.
9. ГОСТ 6258-85 Нефтепродукты. Метод определения условной вязкости. - М.: Изд-во ФГУП «Стандартинформ», 2006. - 7 с.
10. Гиль А.В., Старченко А.В., Заворин А.С. Применение численного моделирования топочных процессов для практики перевода котлов на непроектное топливо: монография. - Томск: STT, 2011. - 184 с.
11. Бубенчиков А.М., Старченко А.В. Численный анализ аэродинамики и горения турбулентной пылеугольной горелочной струи //Физика горения и взрыва. - 1997. - Т. 33. - № 1. -С. 51-59.
12. Бубенчиков А.М., Старченко А.В. Численные модели динамики и горения аэродисперсных смесей в каналах. - Томск: Изд-во Том. ун-та, 1998. - 236 с.
13. Spalding D.B. Mathematical models of turbulent flames: a review // Combustion Science and Technology. - 1976. -V. 13. - № 1-6. - P. 3-22.
14. Yaxin Su, Cuiwu Chen, Along Su. Simulation of High Temperature Air Combustion with modified Eddy-Break-Up combustion model // Energy Procedia. - 2012. - V. 14. - P. 127-132.
15. Abolfazl Irannejad, Farhad Jaberi. Large eddy simulation of turbulent spray breakup and evaporation // International Journal of Multiphase Flow. - 2014 May. - V. 61. - P. 108-128.
16. Westbrook C.K., Dryer F.L. Simplified reaction mechanisms for the oxidation of hydrocarbon fuels in flames // Combust. Sci. Technol. - 1981. - V. 27. - № 1-2. - P. 31-34.
17. Zhou L.X., Hu L.Y., Wang F. Large eddy simulation of turbulent combustion using different combustion models // Fuel. - 2008. -V. 87. - № 13-14. - P. 3123-3131.
18. Murrone A., Scherrer D. Large eddy simulation of turbulent pre-mixed flame stabilized by a backward facing step // 1st INCA Workshop. - Villaroche, France, 2005. - P. 1-9.
19. Шумихин А.А., Карпов А.И. Численное моделирование турбулентного диффузионного пламени на основе метода крупных вихрей // Вычислительная механика сплошных сред. -2012.- Т. 5. - №2. - С. 199-207.
20. Седов Л.И. Механика в СССР за 50 лет. М.: Наука, 1970. -Т. 2. - 880 с.
21. Бояршинов М.Г., Балабанов Д.С. Вычислительное моделирование движения сжимаемой среды, генерируемой точечным источником // Вычислительная механика сплошных сред. -2010. - Т. 3. - №3. - С. 18-31.
Поступила 14.07.2014 г.
UDC 621.181:519.876
NUMERICAL INVESTIGATION OF BURNING RESERVE FUEL IN COMBUSTION CHAMBER OF BOILER BKZ-420-140
Andrey V. Gil,
Cand. Sc., Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia. E-mail: andgil@tpu.ru
Aleksandr S. Zavorin,
Dr. Sc., Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia. E-mail: zavorin@tpu.ru
Dmitry V. Lebed,
Primorskaya Generatsiya, Far East Generating Company, 47A, Fadeev street, Vladivostok, 690034, Russia. E-mail: lebed-dv@dvgk.rao-esv.ru
Aleksandr V. Starchenko,
Dr. Sc., Tomsk State University, 36, Lenin Avenue, Tomsk, 634050, Russia. E-mail: starch@math.tsu.ru
Relevance of the work is caused by the need to assess the effectiveness and reliability of the boiler units on reserve fuel after reconstruction associated with a change in the base fuel and approaches of numerical analysis.
The main aim of the research is to analyze physical and chemical processes in the furnace volume of boiler BKZ-210-140 operating on reserve fuel and rated load when using the Euler and Euler combined and Lagrangian modeling approaches; to construct the mathematical models of the object of study, taking into account the number of rational settlement areas and application of modern numerical algorithms and approaches; to compare numerical results obtained by using different software packages and different representations of the combustion of liquid fuel in the power boiler furnace; to evaluate temperature fields in the volume of the combustion chamber, the concentration of chemical elements of the combustion products, formation of a wind flow patterns and efficiency of backup fuel burn. Methods: numerical simulation based on the Euler and Lagrangian approaches, the finite volume method, wall functions, the two-parameter k-e turbulence model, the model of the chemical reaction in a turbulent flow Eddy-Break-Up - EBU, P1 approximation of spherical harmonics method, the theoretical calculation based on the standard method of thermal calculation of boilers. Results: numerical representation of physical and chemical processes in the furnace of the boiler BKZ-210-140 at reserve fuel combustion; comparison of the results obtained by different numerical models; assessment of the reliability of evaporator surface functioning and efficient of operation of the reserve fuel boiler unit; recommendations for improving the combustion mode in order to increase the operational efficiency of the boiler BKZ-210-140 when combusting fuel oil.
Key words:
Combustion chamber, burner device, twist parameter, gas dynamics, temperature, fuel oil. The paper was carried out within the R&D State Task no. 2069 (2.1322.2014).
REFERENCES
1. Ghorbani A., Bazooyar B. Optimization of the combustion of SOME (soybean oil methyl ester), B5, B10, B20 and petrodiesel in a semi industrial boiler. Energy, 2012 August, vol. 44, Iss. 1, pp. 217-227.
2. Daho T., Vaitilingom G., Sanogo O. Optimization of the combustion of blends of domestic fuel oil and cottonseed oil in a non-modified domestic boiler. Fuel, 2009 July, vol. 88, Iss. 7, pp. 1261-1268.
3. Belousov V.N., Smorodin S.N., Smirnov O.S. Toplivo i teoriya goreniya [Fuel and combustion theory]. V. 2. Teoriya goreniya [Combustion theory]. St. Petersburg, SPbGTURP Publ., 2011. 142 p.
4. Lebedev B.V., Karyakin S.K. Tekhnologiya szhiganiya orga-nicheskikh topliv [Technology of burning fossil fuels]. Tomsk, Tomsk Polytechnic University Publ. House, 2012. 148 p.
5. Gil A.V., Starchenko A.V. Mathematical modeling of physical and chemical processes of coal combustion in chamber furnaces of boiler aggregates based on the package of applied programs FIRE 3D. Thermophysics and Aeromechanics, 2012, vol. 19, no. 3, pp. 503-519.
6. Askarova A.S., Ryspaeva M.Zh., Voloshina I.E. Chislennoe mode-lirovanie goreniya i samovosplameneniya dvukhfaznykh khi-micheskikh reagiruyushchikh techeniy s vpryskami [Numerical simulation of combustion and self-ignition of two-phase reacting chemical flows]. Bulletin of the Tomsk Polytechnic University, 2009. vol. 315, no. 4, pp. 5-9.
7. González-Espinosa A., Lozano A., García J.A., Barreras F., Lincheta E. Design modification of the air diffuser in the burners of a fuel oil power plant. P. II: Interaction with the liquid spray. Fuel, 2014 May, vol. 124, pp. 141-150.
8. Teplovoy raschet kotlov [Thermal computation of boilers] (Normative method). St. Peterburg, 3'd ed., NPO TsKTI publisher, 1998. 256 p.
9. GOST 6258-85. Nefteproduktyi. Metod opredeleniya uslovnoy vy-azkosti [State Standard 6258-85. Petroleum products. The method of determining the relative viscosity]. Moscow, Standar-tinform Publ., 2007. 10 p.
10. Gil A.V., Starchenko A.V., Zavorin A.S. Primenenie chislennogo modelirovaniya topochnykh protsessov dlya praktiki perevoda kot-lov na neproektnoe toplivo [Application of numerical simulation of furnace processes for the practice of transferring fuel boilers non-project]. Tomsk, STT, 2011. 184 p.
11. Bubenchikov A.M., Starchenko A.V. Chislennyy analiz aerodina-miki i goreniya turbulentnoy pyleugolnoy gorelochnoy strui [Numerical analysis of the aerodynamics and combustion of pulverized coal burner turbulent flow]. Fizika goreniya i vzryva, 1997, vol. 33, no. 1, pp. 51-59.
12. Bubenchikov A.M., Starchenko A.V. Chislennye modeli dinamiki i goreniya aerodispersnykh smesey v kanalakh [Numerical models of the dynamics and combustion of aerodisperse mixtures in channels]. Tomsk, Tomsk University Press, 1998. 236 p.
13. Spalding D.B. Mathematical models of turbulent flames: a review. Combustion Science and Technology, 1976, vol. 13, no. 1-6, pp. 3-22.
14. Yaxin Su, Cuiwu Chen, Along Su. Simulation of High Temperature Air Combustion with modified Eddy-Break-Up combustion model. Energy Procedia, 2012, vol. 14, pp. 127-132.
15. Irannejad A., Jaberi F. Large eddy simulation of turbulent spray breakup and evaporation. International Journal of Multiphase Flow, 2014 May, vol. 61, pp. 108-128.
16. Westbrook C.K., Dryer F.L. Simplified reaction mechanisms for the oxidation of hydrocarbon fuels in flames. Combust. Sci. Technol, 1981, vol. 27, no. 1-2, pp. 31-34.
17. Zhou L.X., Hu L.Y., Wang F. Large eddy simulation of turbulent combustion using different combustion models. Fuel, 2008, vol. 87, Iss. 13-14, pp. 3123-3131.
18. Murrone A., Scherrer D. Large eddy simulation of turbulent pre-mixed flame stabilized by a backward facing step. 1" INCA Workshop. Villaroche, France, 2005. pp. 1-9.
19. Shumikhin A.A., Karpov A.I. Chislennoe modelirovanie turbu-lentnogo diffuzionnogo plameni na osnove metoda krupnykh vikhrey [Numerical modeling of turbulent diffusion flame on the basis of large-eddy]. Computational continuum mechanics, 2012, vol. 5, no. 2, pp. 199-207.
20. Sedov L.I. Mekhanika v SSSR za 50 let [Mechanics in the USSR for 50 years]. Moscow, Nauka Publ., 1970, vol. 2, 880 p.
21. Boyarshinov M.G., Balabanov D.S. Vychislitelnoe modelirovanie dvizheniya szhimaemoy sredy, generiruemoy tochechnym is-tochnikom [Computational modeling of motion of a compressible medium, generated by a point source]. Computational continuum mechanics, 2010, vol. 3, no. 3, pp. 18-31.
Received: 14 July 2014.
