Гидравлический расчет змеевика трубчатой печи

 

Целью гидравлического расчета является определение общего гидравлического сопротивления змеевика печи или давление сырья на входе в змеевик, которое, в свою очередь, необходимо для выбора сырьевого насоса.

Давление сырья на входе в печь складывается из следующих составляющих:

P0 = Pk + ΔPи + ΔРн + ΔРкРст., (59)

где Pk, ΔPи, ΔРн, ΔРк, Рст. – соответственно давление сырья на выходе из змеевика печи, значение которого приводится в исходных данных; потери напора: на участке испарения, на участке нагрева радиантных труб, в конвекционных трубах; статический напор.

Расчет необходимо начинать с определения потерь напора на участке испарения:

ΔPи = Рн - Рк., (60)

где Рн – давление в начале участка испарения, которое, в свою очередь, рассчитывается методом последовательного приближения, используя уравнение Бакланова:

, (61)

где А и В – расчетные коэффициенты:

; (62)

(63)

где λ, L1, , dвн, e, ρп – соответственно коэффициент гидравлического сопротивления (для атмосферных печей λ = 0,02 – 0,024), секундный расход сырья по одному потоку, плотность сырья при средней температуре на участке испарения (определяется по уравнению (30)), внутренний диаметр труб, доля отгона сырья на выходе из змеевика, средняя плотность паров при давлении 9,1 Па (при нагреве нефти 1/ρп = 3500);

lи – длина участка испарения.

, (64)

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

(65)

tср.и – средняя температура на участке испарения.

; (66)

lрад. – эквивалентная длина радиантных труб

, (67)

np – число радиантных труб, приходящихся на один поток.

(68)

n – число потоков, Np – общее число радиантных труб.

(69)

lэ – эквивалентная длина печного двойника (ретурбента);

d – наружный диаметр трубы;

lp – рабочая длина одной трубы.

lэ = 50∙d.

lэ = 50∙0,179 = 8,95 м

.

Принимаем Np = 74.

.

 

Тогда:

.

.

Суть метода расчета (метод итераций) Рн по уравнению Бакланова заключается в следующем: задаются значения Рн; по заданной зависимости

Рн = f(tн)

определяют соответствующую этому tн; рассчитывают коэффициенты А и В, длину участка испарения и получают расчетное значение давления в начале участка испарения; если это значение не совпадает с определенной точностью с заданным значением Рн, то расчет возобновляется; при достижении заданной точности фиксируется значение Рн и определяются потери напора на участке испарения.

Таблица 4 – Результаты расчет Рн методом итераций.

Рн, задан. Па tи, оС tср.и, оС , кг/м3 А qtн, кДж/кг lи, м В Рн, расч. Па
900000,0000 268,7643 314,3822 711,5954 101,7209 605,8630 792,3507 1286017,0937 917592,8164
917592,8164 270,2101 315,1051 711,1328 101,7871 609,8191 783,9177 1299851,3949 912832,2301
912832,2301 269,8222 314,9111 711,2569 101,7693 608,7567 786,1823 1296107,2628 914113,0549
914113,0549 269,9268 314,9634 711,2234 101,7741 609,0431 785,5717 1297114,6149 913767,9016
913767,9016 269,8986 314,9493 711,2324 101,7728 608,9660 785,7362 1296843,1556 913860,8729
913860,8729 269,9062 314,9531 711,2300 101,7731 608,9868 785,6918 1296916,2765 913835,8271
913835,8271 269,9042 314,9521 711,2307 101,7731 608,9812 785,7038 1296896,5782 913842,5741
913842,5741 269,9047 314,9524 711,2305 101,7731 608,9827 785,7006 1296901,8846 913840,7565
913840,7565 269,9046 314,9523 711,2305 101,7731 608,9823 785,7014 1296900,4551 913841,2461
913841,2461 269,9046 314,9523 711,2305 101,7731 608,9824 785,7012 1296900,8402 913841,1142

 

Приведем пример расчета Рн для последней итерации.

Принимаем Рн = 913841,2461 Па. По зависимости Рн от tн Рн = f(tн) определяем температуру в начале участка испарения:

tи = 269,9046 оС.

Средняя температура на участке испарения:

.

Плотность сырья при средней температуре на участке испарения:

.

Коэффициент А:

Теплосодержание сырья при температуре начала испарения:

.

Длина участка испарения:

.

 

Коэффициент В:

.

Расчетное значение Рн:

При полученном значении Рн определяем потери напора на участке испарения:

.

Потери напора на участке нагрева радиантных труб:

(70)

где λ2 – коэффициент гидравлического сопротивления.

lн – эквивалентная длина участка нагрева радиантных труб по одному потоку:

(71)

.

ρж – плотность продукта при средней температуре на участке нагрева радиантных труб:

(72)

.

.

U – массовая скорость продукта в радиантных трубах (в случае одинакового размера труб) на один поток:

(73)

.

Следовательно:

.

Потери напора в конвекционных трубах для одного потока:

(74)

где lк – эквивалентная длина конвекционных труб.

(75)

.

ρж – плотность продукта при средней температуре в конвекционных трубах:

.

Тогда:

.

Определение статического напора в змеевике печи:

, (76)

где hT, hK, g – соответственно высота камеры радиации, высота камеры конвекции, плотность продукта при средней температуре:

hT = (np – 1)∙S1 + 0,5∙S1 + 2∙0,25 (77)

hT = (37 – 1)∙0,300 + 0,5∙0,300 + 2∙0,25 = 11,45 м.

.

Тогда:

.

Давление сырья на входе в печь:

.

График зависимости давления в начале участка испарения от температуры представлен на рисунке 4.

Схема, поясняющая расчет змеевика трубчатой печи, представлена на рисунке 5.


 

Рисунок 4 – график зависимости давления в начале участка испарения от температуры.


 

 

 

 

Рисунок 5 – схема, поясняющая расчет змеевика трубчатой печи.