УДК 536.3:535.34
Д. Б. Вафин, А. М. Абдуллин
СЛОЖНЫЙ ТЕПЛООБМЕН В ТЕХНОЛОГИЧЕСКИХ ПЕЧАХ
НЕФТЕХИМИЧЕСКОЙ ПРОМЫШЛЕННОСТИ
Описывается дифференциальный метод расчета радиационно — конвективного теплообмена в топках трубчатых печей, основанный на численном решении системы двухмерных уравнений переноса излучения, энергии, движения, k —е модели турбулентности и простой модели горения газообразного топлива. Перенос излучения рассматривается в рамках S4 — приближения метода дискретных ординат. Спектр излучения продуктов сгорания описывается шестиполосной моделью. Приведены некоторые результаты расчетов. In this paper is given a differential method of calculation for Radiative and convective heat exchange in the furnaces through numerical solving a set of two dimensional equations for transfer of radiation, energy, motion, k — e turbulence models, as well as conventional gas fuel combustion model. Radiation transfer is studied using S4 — approximation of discrete ordinates method. The spectrum of radiation for combustion products is shown using a six array model. Here are also given some calculation figures.
На нефтехимических установках наиболее часто применяют печи с вертикально расположенными трубчатыми змеевиками. Имеются печи с беспламенными горелками на излучающих стенах (рис.1а), с горелками настильного (рис.1 б) и открытого сжигания газообразного топлива (рис.1 в, г). Рассматриваемые в данной работе трубчатые печи характеризуются малой шириной радиационной камеры по сравнению с ее длиной и высотой, симметричным расположением трубчатого экрана и ряда горелок. В этом случае изменение параметров потока по длине намного меньше, чем по ее ширине и высоте. Поэтому задачу теплообмена и газовой динамики продуктов сгорания можно рассматривать в двухмерной постановке. При этом трубчатый экран заменяется непрозрачной для излучения лучевоспринимающей поверхностью с эффективной степенью черноты.
Рис. 1 - Упрощенные схемы трубчатых печей
Для определения лучистых потоков к реакционным трубам необходимо решить ин-тегро-дифференциальное уравнение переноса излучения:
бУI, (М,Л) + к,(М) I, (М,Л) = ая(М) I», (0) + 4® 1I, (М,Л)у,(Л,Л&)ЛССП, (1)
где /, (М, Л) - спектральная интенсивность излучения в точке М в направлении Л; !ь, (0) - спектральная интенсивность излучения черного тела при температуре Т; к,(М)
= а, + в, - спектральный коэффициент ослабления; а, (М), в, (М) - спектральные коэффициенты поглощения и рассеяния; у,(Л,Л&~) - индикатриса рассеяния.
На граничных поверхностях необходимо учитывать излучение и отражение стенки
I,(М,Л) = £д(М,Л)I»,(0*)= П \\I,(М,Л>„(М,Л>,(М,Л ,Л)пЛ&СП (2)
Л&п>0
при Лп < 0. Здесь £я(М,Л), Г* (М,Л&) - спектральные степень черноты и отражательная способность стенки в граничной точке М; р^(М,Л,Л&) - индикатриса отражения стенки.
При использовании метода дискретных ординат уравнение (1) заменяется системой дифференциальных уравнений относительно интенсивности излучения !т вдоль ограниченного количества выделенных направлений [1].
Поле температуры определяется в результате решения уравнения энергии:
Р°р иI0 + Рср° |у = дх(^у6 Ж) + у6 °а И0 + (^ - (3)
где и, и - компоненты скорости продуктов сгорания вдоль осей х и у ; Ср - изобарная теплоемкость; Аэф = , +,т - коэффициент эффективной теплопроводности; - объемная
плотность источников тепла; -мощность плотности лучистых потоков; а = 0 - для
плоской геометрии и а = 1 - для цилиндрически симметричных задач.
Поле течения определяется в ходе решения уравнений движения:
ри % + + £Ш2ди - 2V)) + у--д-(уа„,ф(| + ди) + *. (4)
риди+р»%=-!Р + д. « §+§))+(м ди - ^))+
+ !Ь- (2-У- - 2 ^іу V) + /2 , (5)
у- х уа 3 ’ ’ 4 ’
где Аэф = Я +Ят - эффективная вязкость; f1, /2 - массовые силы. Если ось х направлена по вертикали вниз, то f1 = (1 - в(Т -- Т»)), /2 = 0, где в = ррдр _ коэффициент объемного расширения; д -ускорение свободного падения.
К уравнениям движения добавляются уравнение неразрывности и уравнение состояния газа
прт + у- ду (= °- Р = тргкт - <6>
Для замыкания системы уравнений (для определения коэффициентов турбулентного переноса) используется двухпараметрическая диссипативная к-є модель турбулентности. Уравнения для кинетической энергии турбулентных пульсаций к и скорости ее диссипации є имеют вид
д1(р иф) + (у-риф) = д^(ГфдФ) + (уаГфдф) + Бф,
доУИ о- додо ф до у- ду ф до ф
где ф = {к, а}; Гф = « + «п /оф - коэффициент переноса; Бф - источниковый член. Эмпирические константы модели взяты по рекомендациям работы [2].
Распределение источников тепловыделений в уравнении (3) зависит от характера выгорания газообразного топлива. В данной работе используется модель простой химической реакции, согласно которой горение предварительно перемешанной газовой смеси описывается уравнениями для массовой концентрации горючего тд и окислителя тш .
Уравнение для тд имеет вид:
(р ит§) + СадС(Гро т§) = (ГдтЯ) + -1 дС(Гуадт) + Б§. (8)
Такому же уравнению удовлетворяет уравнение ток Источниковый член уравнения для окислителя определяется соотношением Бш = Бд А , где А - стехиометрическое количество окислителя для сгорания 1 кг горючего. Из этих двух уравнений в предположении равенства коэффициентов переноса (А = Аш ), получается уравнение для переменной р г = т§ - тш / А с нулевым источниковым членом. Коэффициент переноса в (8) Ад = «/ ог, где ог - число Шмидта. Скорость химической реакции (Источниковый член Бд) определяется по модели «обрыва вихрей» [3].
Состав и температура получаемого целевого продукта, температура наружной стенки реакционных труб определяется условиями протекания процессов как внутри труб, так и в объеме топочной камеры. Поэтому метод теплового расчета радиантной камеры трубчатой печи должен основываться на совместном решении задач внешнего и внутреннего теплообмена. Эффективная степень черноты трубчатого экрана вычисляется по формуле:
аф = р»( 2-Р,,,К
аф е»+пР|,у(2_Р|,()г» ’ ( )
где рп,э - угловой коэффициент, зависящий от шага между трубами Ээ и внешнего диаметра трубы б. В данной работе для изучения сопряженного теплообмена рассматривается процесс паровой конверсии углеводородного сырья. Для расчета внутреннего теплообмена использована методика, описанная в работе [4]. В основу метода положена равновесная модель, учитывающая следующие реакции окисления метана:
СН4 + Н2О о- СО + ЗН2 - 206.4 кДж/моль, (10)
СО + Н2О о СО 2 + Н2 Н2 + 41.0 кДж/ моль. (11)
Для численного решения дифференциальных уравнений метода дискретных ординат (3) применяется конечно-разностный метод [5]. Решение уравнений движения (4) и (5), описывающих конвективный теплообмен и поле турбулентного течения продуктов сгорания в топочной камере прямоугольного сечения, выполняется в (и-и- р) переменных [6]. При этом возникают некоторые трудности, связанные с составлением дискретных аналогов для членов уравнений в виде первых производных. Для преодоления этих трудностей составляющие скорости определены на разностной сетке, отличной от сетки, которая используется для всех других переменных, т.е. использовали смещенную шахматную сетку. Поле течения в цилиндрических топках определяется решением уравнений движения в переменных функция тока - вихрь ( щ - о) [7]. Для получения конечно-разностного аналога уравнения энергии используется метод контрольного объема. Диффузионные члены аппроксимируются центральными разностями, для конвективных членов используется про-тивопоточная схема. В ходе итерационного процесса совместного решения уравнений энергии и переноса излучения в первых внешних итерациях поле температуры сильно «осциллирует». Для уменьшения величины «осцилляций» применяется нижняя релаксация и линеаризация источникового члена.
Алгоритм расчета сопряженного теплообмена можно выразить следующей итерационной схемой:
Для работы пакета прикладных программ необходимы следующие данные:
Выходными данными пакета прикладных программ являются:
На рис. 2 показаны поля продольной скорости и температуры для половины камеры
радиации трубчатой печи коробчатого типа для паровой конверсии метана БПК-9000. В области факела приведены также линии постоянных значений массовой концентрации горючего в процентах относительно входного сечения: 1%; 10%; 60%. На рис. 3 приведены распределения плотностей лучистых и конвективных потоков по длине реакционных труб. Для учета радиационных свойств продуктов сгорания использована «селективно-серая» модель Хоттеля и шестиполосная модель широкой полосы. Распределения лучистых потоков, рассчитанные по модели Хоттеля и шестиполосной модели на начальном участке сильно отличаются, что можно объяснить наличием окон прозрачности и сильной неизо-термичностью потока на этом участке. Необходимо отметить, что основное количество тепла трубчатому экрану передается за счет теплового излучения.
Рис. 2 - Поля продольной скорости и температуры в поперечном сечении топки: скорость на входе ивх = 1.52 м/с; Твх = 1269 К; сплошная линия - скорость; пунктирная линия - температура
Рис. 3 - Распределение плотностей тепловых потоков по длине труб, профили скоростей и температуры потока: 1 — плотность лучистых; 2 — конвективных тепловых потоков (пунктирная линия — по шестиполосной модели; сплошная линия — по модели Хоттеля)
© Д. Б. Вафин - канд. техн. наук, доц. каф. физики Нижнекамского химико-технологического института КГТУ; А. М. Абдуллин - канд. техн. наук, доцент той же кафедры.