Целью гидравлического расчета является определение общего гидравлического сопротивления змеевика печи или давление сырья на входе в змеевик, которое, в свою очередь, необходимо для выбора сырьевого насоса.
Давление сырья на входе в печь складывается из следующих составляющих:
P 0 = Pk + Δ Pи + Δ Рн + Δ Рк – Рст., (59)
где Pk, Δ Pи, Δ Рн, Δ Рк, Рст. – соответственно давление сырья на выходе из змеевика печи, значение которого приводится в исходных данных; потери напора: на участке испарения, на участке нагрева радиантных труб, в конвекционных трубах; статический напор.
Расчет необходимо начинать с определения потерь напора на участке испарения:
Δ Pи = Рн - Рк., (60)
где Рн – давление в начале участка испарения, которое, в свою очередь, рассчитывается методом последовательного приближения, используя уравнение Бакланова:
, (61)
где А и В – расчетные коэффициенты:
; (62)
(63)
где λ, L 1, , 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)∙ S 1 + 0,5∙ S 1 + 2∙0,25 (77)
hT = (37 – 1)∙0,300 + 0,5∙0,300 + 2∙0,25 = 11,45 м.
.
Тогда:
.
Давление сырья на входе в печь:
.
График зависимости давления в начале участка испарения от температуры представлен на рисунке 4.
Схема, поясняющая расчет змеевика трубчатой печи, представлена на рисунке 5.
Рисунок 4 – график зависимости давления в начале участка испарения от температуры.
Рисунок 5 – схема, поясняющая расчет змеевика трубчатой печи.