CC BY
17Бюллетень науки и практики /Bulletin of Science and Practice Т. 6. №6. 2020
https://www.bulletennauki.com https://doi.org/10.33619/2414-2948/55
УДК 536.3.535.34 https://doi.org/10.33619/2414-2948/55/05
ВЛИЯНИЕ СПЕКТРАЛЬНОЙ МОДЕЛИ ИЗЛУЧЕНИЯ НА РАСЧЕТНЫЕ ХАРАКТЕРИСТИКИ СЛОЖНОГО ТЕПЛООБМЕНА В ПЛАМЕННЫХ ПЕЧАХ НЕФТЕХИМИЧЕСКОЙ ПРОМЫШЛЕННОСТИ
©Абдуллин А. М., ORCID: 0000-0003-1223-0535, канд. техн. наук, Казанский национальный исследовательский технологический университет, г. Нижнекамск, Россия, amabdullin@mail.ru
INFLUENCE OF THE SPECTRAL MODEL OF RADIATION ON THE CALCULATED CHARACTERISTICS OF COMPLEX HEAT EXCHANGE IN FLAME FURNACES
OF THE PETROCHEMICAL INDUSTRY
©Abdullin A., ORCID: 0000-0003-1223-0535, Ph.D., Kazan National Research Technological University, Nizhnekamsk, Russia, amabdullin@mail.ru
Аннотация. Анализируется влияние спектральной модели излучения на тепловые потоки и температуру продуктов сгорания в радиантных камерах трубчатых печей нефтехимической промышленности. Рассматриваются модель широкой полосы и серая модель Хоттеля. Показано, что спектральная модель излучения топочной среды слабо влияет на расчетные характеристики суммарной теплопередачи.
Abstract. The influence of the spectral model of radiation on heat fluxes and the temperature of combustion products in the radiant chambers of tube furnaces of the petrochemical industry is analyzed. A wide-band model and a Hottel gray model are considered. It is shown that the spectral model of the combustion medium radiation weakly affects the calculated characteristics of the total heat transfer.
Ключевые слова: интенсивность излучения, лучистый поток, степень черноты, спектр, конвективный поток.
Keywords: radiation intensity, radiant flux, degree of blackness, spectrum, convective flux.
В тепловом балансе радиантных камер трубчатых печей определяющая роль принадлежит лучистому теплообмену, доля которого может быть на порядок больше по сравнению с другими механизмами теплопередачи. Точность расчета лучистого теплообмена в значительной мере определяется корректностью модели спектра излучения компонентов, заполняющих камеру сгорания. К ним относятся, в первую очередь, пары воды ^O и двуокись углерода CO2, а также частицы сажи, образующиеся в результате неполного сгорания топлива. Излучение паров воды и двуокиси углерода происходит в определенных интервалах спектра, излучение частиц сажи и стенок камеры характеризуется сплошным спектром. Одним из наиболее удобных методов учета селективности излучения является модель широкой полосы, параметры которой приведены, например, в работе [1].
В данной работе излагаются результаты сравнительного анализа применимости спектральных моделей широкой полосы и «серой» модели Хоттеля для теплового расчета радиантных камер трубчатых печей. За исходную принимается ступенчатая модель,
Бюллетень науки и практики /Bulletin of Science and Practice Т. 6. №6. 2020
https://www.bulletennauki.com https://doi.org/10.33619/2414-2948/55
учитывающая полосы 15; 4,3; 2,7 мкм излучения двуокиси углерода и 10; 6,3; 2,7; 1,5 мкм излучения водяного пара.
Лучистый теплообмен описывается системой двухмерных дифференциальных уравнений Р1-приближения метода сферических гармоник:
Чп Чп (1)
^xihb(T,-dd-i ;
h
Сх, = -Ъх^Ух.- (2)
Условия однозначности на ограничивающих поверхностях с учетом диффузного отражения и излучения имеет вид:
й(х 1 ( С*2 \ (3)
КГ = 2ЁХ-(ГТТ)(4^I - (1 -г)<рх.)
Здесь С[х1 — вектор плотности потока излучения в ьтой спектральной полосе; , ах1— коэффициенты диффузии излучения и поглощения в ьтой спектральной полосе; ¡хь — функция Планка; е,г — степень черноты и отражательная способность поверхности.
Распределение температуры в расчетной области определяется путем решения уравнения энергии:
(4)
д Т д Т д (4)
PCpU^ + PCPV^,=TX
дТ
а+лт)-
д
ду
д Т
(Л + ЛТ)~
+ Q - div^qx. i=l
В качестве краевых условий на поверхности нагрева и входном сечении задается температура, на выходном сечении используется свойство «односторонних» координат [2]. Температура футеровки определяется из уравнения теплового баланса, составленного с учетом тепловых потерь через нее:
Ь ^ (5)
-ф (Tw - То) = + ^ Axi.
i=l
Конвективный поток тепла к стенке вычисляется по методу пристеночных функций [3]. В уравнениях (4) и (5) р,ср,и,р — плотность, изобарная теплоемкость и составляющие вектора скорости продуктов сгорания; Л, Лт — коэффициенты молекулярной и турбулентной диффузии тепла; Q — объемная плотность тепловыделения в результате горения топлива; — количество излучающих полос в спектральной модели; Лф, й — коэффициент теплопроводности и толщина футеровки, Т/, Т — температура внутренней и внешней поверхности футеровки.
Конечно-разностная аппроксимация уравнений (1)-(3) производится методом конечных элементов, а уравнения (4) методом интегрирования по контрольному объему [2].
Полученная в результате система алгебраических уравнений относительно узловых значений искомых величин решается методом трехдиагонального матричного алгоритма (ТДМА). Проведенные расчеты показывают, что используемая в данной работе явная схема совместного интегрирования системы уравнений (1)-(5) при тепловых нагрузках, характерных для высокотемпературных камер сгорания, является неустойчивой. Для
Бюллетень науки и практики / Bulletin of Science and Practice https://www.bulletennauki.com
Т. 6. №6. 2020 https://doi.org/10.33619/2414-2948/55
подавления этой неустойчивости к дивергенции лучистых потоков тепла применяется нижняя релаксация с параметром о = 0,2.
Анализ проводился на примере прямоугольной радиантной камеры трубчатой печи конверсии углеводородов БПК-15. Топливо (природный газ) подается в камеру через горелки, расположенные на своде в одном ряду симметрично относительно поверхности нагрева, продукты сгорания отводятся через тоннель, расположенный на поду камеры. Наличие осевой симметрии позволяет рассматривать только одну половину радиантной камеры. Полагалось, что поверхности футеровки и нагрева диффузно излучающие и отражающие со степенями черноты 0,6 и 0,87 соответственно. Для получения замкнутой системы уравнений в расчетах входное и выходное сечения заменяются полностью отражающими поверхностями. Поля скоростей и коэффициентов турбулентного обмена определяются из решения уравнений газовой динамики и (к-е) — модели турбулентности [4].
В модели спектра интервалы, где происходит перекрывание полос излучения водяного пара и двуокиси углерода, выделялись как отдельные полосы с суммарным коэффициентом поглощения:
аЕ = а^ + асРс, (6)
где а^, ас — коэффициенты поглощения водяного пара и двуокиси углерода в данном интервале спектра, соответственно, Рмп Рс — парциальные давления каждой из компонент. Учитывалось перекрывание спектральных полос 2,7 мкм; 10 и 15 мкм.
Спектральный коэффициент поглощения частиц сажи вычислялся по формуле [5]:
3яФ(Я)
ux,s
21
■fv
(7)
где функция Ф(Я), описывающая дисперсию оптических констант сажи, по данным работы [6] в интервале спектра до 10 мкм может быть описана зависимостью:
Ф(Л) =
1,75
(8)
1 + 0,3Л
Так называемая объемная фракция сажи определялась по эмпирической формуле [7]:
С
fv = 0,068(2 - aT)-f-.
пр
(9)
Здесь ат — коэффициент избытка воздуха; Ср/Нр — относительный вес углерода в топливном газе.
Для вычисления интегральной степени черноты топочной среды использовалась модель Хоттеля набором четырех серых газов [6]:
= ^ ап(Т){1 - ехр\-(апР1 + а5)/эф]}
(10)
п=0
где 1эф — эффективная толщина топочной среды, (атм.) — суммарное давление водяного пара и двуокиси углерода, интегральный коэффициент поглощения частиц сажи определяется по зависимости [6]:
(11)
С
as = 0,03(0,0167 - 0,5)(2 -aT)-f.
Нр
£
Бюллетень науки и практики /Bulletin of Science and Practice Т. 6. №6. 2020
https://www.bulletennauki.com https://doi.org/10.33619/2414-2948/55
На Рисунках 1-2 приведены расчетное распределение плотности интегральных по спектру лучистых потоков тепла к поверхности нагрева по длине труб и профиль температуры на оси потока по высоте радиантной камеры, определенные с использованием спектральных моделей широкой полосы и Хоттеля. Наибольшее отличие в распределениях лучистых потоков, рассчитанных по модели Хоттеля и полосовой модели без учета влияния сажи, имеет место на начальном участке течения, что можно объяснить существенной неизотермичностью потока, а также наличием «окон» прозрачности в полосовой модели. Ниже по течению, где поперечный градиент температуры становится незначительным, расчет по модели Хоттеля хорошо согласуется с результатами расчета по полосовой модели. Учет излучения частиц сажи приводит к увеличению лучистых потоков тепла и соответственно уменьшению максимальной температуры газов в топке и ее поперечного градиента. В этом случае распределение лучистых потоков, определенное по модели Хоттеля, удовлетворительно согласуется с результатами расчетов, сделанными по модели широкой полосы с учетом трех «окон» прозрачности, в которых лучистая теплопередача происходит в основном за счет излучения частиц сажи.
О 2 4 6 8 L..1
1 — расчет без учета, 2 — с учетом излучения сажи, сплошная линия — модель широкой полосы, пунктиром — серая модель Хоттеля.
Рисунок 1. Распределение плотности Рисунок 2. Профиль температуры
лучистых потоков тепла к поверхности нагрева по продуктов сгорания по высоте радиальной длине труб. камеры.
Наибольшее влияние спектральная модель излучения оказывает на профиль температуры дымовых газов в объеме радиантной камеры. В то же время это влияние на температуру газов, покидающих топку, незначительно. Наличие частиц сажи в составе продуктов сгорания снижает температуру газов на выходе из топки на 16%.
Учет излучения сажи приводит к уменьшению расчетных значений конвективных потоков тепла (Рисунок 3). Сплошной спектр излучения интенсифицирует теплообмен между стенкой и средой, вследствие чего градиент температуры вблизи стенки уменьшается. Расчеты также показывают, что в месте присоединения рециркуляционного и основного течений коэффициент конвективной теплоотдачи имеет минимум. Ниже по течению вблизи выходного сечения увеличивается, что вызвано сужением основного потока и его ускорением.
Бюллетень науки и практики /Bulletin of Science and Practice Т. 6. №6. 2020
https://www.bulletennauki.com https://doi.org/10.33619/2414-2948/55
ЧкКВТ/к2
1.5
1.0
0,5
/ / ^ 2 \_ v\ / \ * — 4X 1
0 2 4 6 а и м
1 — расчет без учета, 2 — с учетом излучения сажи, сплошная линия — модель широкой полосы, пунктиром — серая модель Хоттеля.
Рисунок 3. Распределение плотности конвективных потоков тепла к поверхности нагрева по длине труб.
Анализ суммарной теплопередачи показывает, что интегральное тепловосприятие поверхностью нагрева, рассчитанное по модели Хоттеля, на 13% выше по сравнению с рассчитанным по модели широкой полосы. Этому соответствует снижение температуры дымовых газов на выходе из топки, что качественно согласуется и с данными других авторов [7]. Вместе с тем следует заметить, что спектральная модель излучения топочной среды слабо влияет на расчетные характеристики суммарной теплопередачи в радиантной камере рассмотренной трубчатой печи.
Список литературы:
1. Абдуллин А. М. Анализ спектра излучения продуктов сгорания газообразного топлива // Вестник Казанского технологического университета. 2013. Т. 16. №12. С. 67-70.
2. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Энергоатомиздат, 1984. 152 с.
3. Белов И. А., Кудрявцев Н. А. Теплоотдача и сопротивление пакетов труб. Л.: Энергоатомиздат, 1987. 223 с.
4. Абдуллин А. М. Численный метод определения температуры излучающей стенки в трубчатых печах // Известия высших учебных заведений. Проблемы энергетики. 2011. №1112. С. 30-39.
5. Hubbard G. L., Tien C. L. Infrared mean absorption coefficients of luminous flames and smoke. 1978. https://doi.org/10.1115/1.3450788
6. Руководящий технический материал 26-02-40-77. Нормативная методика теплового расчета трубчатых печей. Введена 01.01.78. 360 с.
7. Блох А. Г. Теплообмен в топках паровых котлов. Л.: Энергоатомиздат, 1984. 240 с.
References:
1. Abdullin, A. M. (2013). Analiz spektra izlucheniya produktov sgoraniya gazoobraznogo topliva. VestnikKazanskogo tekhnologicheskogo universiteta, 16(12), 67-70. (in Russian).
2. Patankar, S. (1984). Chislennye metody resheniya zadach teploobmena i dinamiki zhidkosti. Moscow, Energoatomizdat, 152. (in Russian).
Бюллетень науки и практики / Bulletin of Science and Practice Т. 6. №6. 2020
https://www.bulletennauki.com https://doi.org/10.33619/2414-2948/55
3. Belov, I. A., & Kudryavtsev, N. A. (1987). Teplootdacha i soprotivlenie paketov trub. Leningrad, Energoatomizdat, 223. (in Russian).
4. Abdullin, A.M. (2011). Numerical method for determining the temperature of radiating wall in tube furnaces. Izvestiya vysshikh uchebnykh zavedenii. Problemy energetiki, (11-12), 30-39. (in Russian).
5. Hubbard, G. L., & Tien, C. L. (1978). Infrared mean absorption coefficients of luminous flames and smoke. (in Russian). https://doi.org/10.n15/L3450788
6. Rukovodyashchii tekhnicheskii material 26-02-40-77. Normativnaya metodika teplovogo rascheta trubchatykh pechei. Vvedena 01.01.78. (in Russian).
7. Blokh, A. G. (1984). Teploobmen v topkakh parovykh kotlov. Leningrad, Energoatomizdat, 240. (in Russian).
Работа поступила Принята к публикации
в редакцию 08.05.2020 г. 11.05.2020 г.
Ссылка для цитирования:
Абдуллин А. М. Влияние спектральной модели излучения на расчетные характеристики сложного теплообмена в пламенных печах нефтехимической промышленности // Бюллетень науки и практики. 2020. Т. 6. №6. С. 42-47. https://doi.org/10.33619/2414-2948/55/05
Cite as (APA):
Abdullin, A. (2020). Influence of the Spectral Model of Radiation on the Calculated
Characteristics of Complex Heat Exchange in Flame Furnaces of the Petrochemical Industry. Bulletin of Science and Practice, 6(6), 42-47. (in Russian). https://doi.org/10.33619/2414-2948/55/05
