Численный метод исследования стационарных режимов в технологических печах Текст научной статьи по специальности «Математика»

Journal of Siberian Federal University. Engineering & Technologies 1 (2014 7) 55-61

УДК [665.63+662.764]:62.50

Численный метод исследования стационарных режимов в технологических печах

Н.Д. Демиденкоа, Л.В. Кулагина6*

аСКТБ «Наука» КНЦ СО РАН Россия, 660049, Красноярск, пр. Мира, 53 бСибирский федеральный университет, Россия, 660041, Красноярск, пр. Свободный, 79

Received 18.11.2013, received in revised form 14.12.2013, accepted 12.01.2014

В статье предлагаются математические модели для стационарных и динамических процессов технологических печей как объектов с распределёнными параметрами. Эта математическая модель получена на основе законов сохранения энергии, массы и импульса и включает в себя дифференциальные уравнения в частных производных. Проведено численное исследование стационарных режимов трубчатых печей. Для этого сформулирована краевая задача для расчета концентрации горючего вещества, плотности, скорости, температуры дымовых газов и температуры нагреваемого сырья, идущего на разделение в ректификационную колонну. Представлены результаты расчета технологических параметров. Предложенный метод может быть использован при автоматизации ректификационных установок в нефтеперерабатывающей и нефтехимической промышленности.

Ключевые слова: математическое моделирование, системы с распределёнными параметрами, тепломассообмен, численный метод.

Для высокотемпературного нагрева нефти и нефтепродуктов в процессе их переработки в ректификационных установках применяют трубчатые печи. Печи работают следующим образом [1]: мазут и газ сжигают в горелках, расположенных в камере радиации. Продукты сгорания из камеры радиации поступают в камеру конвекции, затем направляются в газосборник и по дымовой трубе уходят в атмосферу. В камере конвекции расположены конвекционные трубы, воспринимающие тепло при соприкосновении дымовых газов с поверхностью нагрева путем конвекции. Нагреваемый продукт в печи последовательно проходит через конвекционные и радиантные трубы, поглощая тепло. Радиант -ная поверхность воспринимает большую часть тепла, выделяемого при сгорании топлива (рис. 1).

Горение жидкого топлива всегда протекает в паровой фазе, причем испарение обеспечивается подводом тепла из зоны горения паров. Топливо всегда сжигается в распыленном виде.

© Siberian Federal University. All rights reserved

* Corresponding author E-mail address: klvation@gmail.com

Рис. 1. Схема трубчатой печи: 1 - горелка; 2 - настильная стенка; 3 - камера радиации (топочная каме ра); 4 - камера ко нвекции; 5 -дымовая труба; 6 - змеевик конвекционных труб; 7 - змеевик радиантных труб; 8 -футеровка. Потоки: I - вход сырья; II -выход сырья; III - топливо и воздух; IV -дымовые газы

Крайним случаем горения капель является диффузионное горение, когда скорость сгорания паров очень велика по отношению к скорости диффузии паров и окислителя в зоне горения, толщина которой становится исчезающе малой. К такому режиму может приближаться горение сравнительно круглых капель. Диффузионная теория горения развита Г.А. Варшавским [2].

В математических методах анализа процесса горения имеются некоторые расхождения у разных авторов, но для стационарного сферического горения используется единый подход. Для упрощения анализ проводится при следующих предпо -ложениях:

1. Жидкая капля имеет сферическую форму.

2. Влиянием конвекции пренебрегают, пламя рассматривают как сферическую поверхность, концентрическую с каплей.

3. Пламя считают разновидностью диффузионного пламени, которое образуется в результате реакции между парами горючего и воздухом, реагирующими в стехиометрическом соотношении.

4. Рассматривают стационарное состояние при постоянном диаметре капли, хотя реально диаметр жидкой капли уменьшается по мере горения, однако это изменение происходит медленно по сравнению с изменением скорости диффузии и прочими факто-

рами.

5. Температура капли одинакова по всему объему.

6. Давление в течение в сего процеско горения считается постоянным.

7. Влияние излучения рассматривают отдельно.

Исоодя из зянинов механики сплошных сред, можно получить следующие уравнения не-

стац онарного горения:

др п

дх

др

дх х

ди

дС и д1 рд1'

дС и д1

ди

_=_ ди дГп

дС и д1 д1~

дТк=(1_ )Т—_ д7П . хд С(Тп)

дС ( У) п д1 и д1 Сут Сур

Я7Л др ' р дГ

+ КОТ1 _ Тп) + Кг(Тс2 _ Тп)'

дТ1 дТ1 1ч

~дс=™1)Г + К2(Тп_ Тс1) _ С (Тп, Тс1)'

дТс2 дТс2 2Ч , 2Ч

~д^ = _™ИГ + К(Тп_ Тс2) _ С (Тп, Тс2).

Здесь x(l,t) - концентрация горючего вещества, p(l,t), u (l,t), Tn(l,t) - плотность, скорость и тем пература дымовых газяв; T'c(l,t), Vc(l,t) - температура нисходящего и восходящего потоков нагреваемого аырья; Г, l - временная ипространственная коьрдинаты. Начальные условия

р((,0) = Ро,х(),0) = Х0,и(1, 0) = U0,Ta(i,0) = Tao,Tc1(i,0) = ТТо,Т°аО) = Тсо0. Граничные условия

p(0,t) = а1, x (0, t )= а2, u (0, t) = а3, Тп (0, t) = а4, Tc\ L, t) = а5, Tc2(0, t) = Tc\0, t). Здесь L - длина печи.

Стационарную модель, которая следует из (1)), рассмотрим боле е подробно:

др ди -и — - р — = 0, д) д)

дх х -и—-- = 0, д) г

ди _ дТа ЙТадр_^

еи-теЯ~дГ~~-де-°-

ч ди дТа х <3(Та) , 1 о дТ1

(2)

(3)

(4)

W-д^ + ХоОГ - тт1) -о (ЗаГа,Гс0) = 0, дТо

-иу-дц- -и- ад - т-0) - ¿зота,ГЖ) = 0-

Второе, пятое и шестое уравнения содержат пространственные производные то лько одной неизвестной функции.

Приведем первое, третье и чеовертое уравнения к виду, -одержащему производную толяко ят одной искомой (функции dp/dt, du/dl и dTn/dl. Для этого применим теорию о неявных функциях к системе.

Зр ди З) З)

Пи З- Й-пЗр и—иП—Й-^ = 0,

З) З+ р З)

Пи ЗГП х р(-п) ,

(LL--->- эГи и "Ж = " + + ^СРТ -- 1+ - -- -- ) = * ■

(5)

В этом случае

где

Пр _ДХ Пи _ Д2 ЗГп _ д3

Ж _ ~Д'"П) _ "Д,ПППТ _ ~д'

(6)

д _

и

Й-п

- (1 - г)-п -и

ЙТ

_ -и3 + ир-(-- иЙ(1 - г)+п = -и3 + -Йи-п - и)уД-п - и2), (7)

Д1 _

0 р о

О и Й <р (1-у)-, -и

_ Йрр, Д2_

и 0 0

Й-п

— 0Й

р

0 ^ -и

_ -Йи^>,

(8)

U

RT„

Р

U

= U2(p.

Р1-У)Т„ V

(9)

Таким образом, будем иметь:

ду Rp(

~dl = ~K = u(yRTn - и2)' ди У2 Rucp R(p

~д=У1=~ ~ u(yRT„-и2)'' yRTn - и2' дТу У3 и2( и(р

~Ж = ~К = u(yRT„ - и2) = yRТп - и2'

Тем самым система (4) приведена к следующему виду:

аг и(удтп - и2) ■

Зх х

и(- X "с^ ., ;(-п) + С„р ■К1(-с1+-с2 - 2-п))

УЙ-п -и2

-й(- X "С,г£ -Ki(^+-C2 - 2-п))

Уй-п" -и2

К—i К 1

К-2 К, „ 1 , „ ^=-^(-n--2)--Q(---:2).

р(0) = а, *(0) = а2, «(0) = а3, Э„ (0) = а4, Т1С (Ь) = а5, Э/ (0) = Т1С (0). (12)

На рис. 2-55 изображены результаты расчетов стационарных режимов при различных начальных условиях для плотности р(1), концентрации с31), скорости и(1) из температуры дымо -вых газов 7Ь(7=. При этом за начальные условия приняты сЬ0)=лСС кг/м3, х(0)=О,47, и(0)=5 м/с, Тп(0)=450 °С, Т1/¿3=270 °С. Затем начальные данные изменялись с шагом ±5 % при постоянных значениях остальных параметров.

Предлагаемый метод анализа стационарных режимов технологических печей позволяет еще на стации проектирования ректификационных установок обеспечить высокое качество разделения многокомпонентных смесей в нефтепереработке и нефтехимии.

Ii IS 1t

5 «

5 ld

■2i

Ii 2 j i t i« и ii ie

Рис. 2. Графики распределения плотности дымовых газов по длине объекта

Рис. 3. Графики изменения концентрации по длине объекта

им/с

О 2 4 б 8 10 12 14 16 18 1, м

Рис. 4. Графики изменения скорости дымовых газов по длине объекта

Тн^С

ЦП ulu

и*

HP

tft ttt

V

ist V

Мп

Hfl б I II В М it J" з

Рис. 5. Графики изменения температуры дымовых газов по длине объекта

Список литературы

[1] Скобло А.И., Молоканов Ю.К., Владимиров А.И., Щелкунов В.А. Процессы и аппараты нефтегазопереработки и нефтехими и. М.: ООО «<Недра-Бизнесцентр»», 2000. 677 с.

[2] Варшавский Г.А. II Бюро новой техники НКАП. М.: Гостехиздат, 19)456. № 6. С. 87-106.

[3] Демиденко Н.Д. II ЗЕВестник Томского государственного университета. Управление вычислительная техника и информатика. 20122. № 3 (20). С. 135-21.

[4] Демиденко Н.Д., Кулагин В.А., Шокин ЮИ. Моделирование и вычислите льные технологии распределенных систем. Новосибирск: Наука, 2012. 424 с.

[5] Демиденко Н.Д., Кулагина Л.В. Моделирование и оптимизация технических систем с распределенными параметрами: учеб. пособие для вузов. Красноярск: ИПЦ КГТУ, 2006. 210 с.

Numerical Research Method

of Stationary Modes in Technological Furnaces

Nikolay D. Demidenkoa and Ludmila V. Kulaginab

aSDTB «Nauka» KSC SB RAS 53 Mira, Krasnoyarsk, 660049, Russia bSiberian Federal University, 79 Svobodny, Krasnoyarsk, 660041, Russia

In article we propose mathematical models for stationary and dynamic processes of technological furnaces as object with distributed parameters. This mathematical model is based on the laws of conservation of energy, mass and impulse and includes differential equations in private derivatives. Computational investigation of stationary modes of tubular furnaces is conducted. Formulated boundary-value problem for calculation concentration of combustible substance, density, speed, temperature of theflue gas and temperature of the heated raw materials going on division in rectificative column. Results of calculation technological parameters are presented. The proposed method can be used at automation rectificative plant in petroleum-refining and petrochemical industry.

Keywords: mathematical modeling, distributed parameter systems, heat mass exchange, numerical method.