
- •Расчет скорости пара и диаметра в нижней части колонны:
- •Расчет высоты светлого слоя жидкости в верхней части колонны
- •Расчет высоты светлого слоя жидкости в нижней части колонны
- •Расчет коэффициента диффузии распределяемого компонента в паре в верхней и нижней частях колонны.
- •Расчет коэффициента диффузии распределяемого компонента в жидкости в верхней и нижней частях колонны.
лист 1
Технологический расчет
Расчет тарельчатой ректификационной колонны
Материальный баланс колонны.
Расчет материальных потоков в колонне проводится на основании уравнений материального баланса. Уравнения материального баланса колонны:
где:
F–расход исходной смеси 1,111 кг/c;
W–расход кубового остатка, кг/c;
P–расход дистиллята, кг/c;
хF –концентрация легколетучего компонента в исходной смеси;
xW–концентрация легколетучего компонента в кубовом остатке;
xP–концентрация легколетучего компонента в дистилляте;
Решая систему этих уравнений, находим расход кубового остатка и дистиллята:
Пересчет концентраций
Для расчетов по диаграмме У–Х необходимо перейти от массовых концентраций к мольным:
где хW, xF, xP–массовые концентрации низкокипящего компонента в кубовом остатке, питании, дистилляте соответственно; Мм, Мв–молекулярные массы низкокипящего и высококипящего компонентов.
Расчет минимального флегмового числа
Определяем по диаграмме
У–Х состав пара, равновесного
к составу жидкости в исходной смеси:
Рассчитываем Rmin в соответствии с формулой:
Расчет условно–оптимального флегмового числа
Нагрузки ректификационной колонны по пару и жидкости (и основные геометрические размеры) определяются рабочим флегмовым числом, найдем условно–оптимальное флегмовое число исходя из минимального объема ректификационной колонны по минимальному значению произведения N×(R+1), путем построения графика N(R+1) от R.
Для этого:
Рассчитываем флегмовые
числа по уравнению:с
различными значениями b,
где b–коэффициент
избытка флегмы.
Строим на диаграмме У–Х рабочие линии соответствующие рассчитанным флегмовым числам, вырисовываем ступени между рабочей и равновесной линиями. Считаем теоретические ступени и результаты расчетов сводим в таблицу:
Далее строим график зависимости N(R+1) от R из которого определяем условно–оптимальное флегмовое число: Rопт=1,155
лист 2
Расчет скорости пара и диаметра колонны
Диаметр колонны находим
из уравнения расхода:
где:
Q–объемный расход пара в колонне, м3/с;
G–массовый расход пара в колонне, кг/с;
D–диаметр колонны, м2;
W–скорость пара в сечении колонны, м/с;
ry–плотность пара, кг/м3.
Свойства пара в верхней и в нижней части колонны будут различны, для учета этого факта расчет свойств жидкости и пара, а также основных геометрических размеров колонны будем проводить отдельно для обеих частей колонны.
Средний мольный состав жидкости в верхней и нижней части колонны:
Средняя молекулярная масса жидкости в верхней и нижней части колонны:
Молекулярная масса исходной смеси:
МF=32×0,217+18×(1–0,217)=21,04 кг/кмоль;
Молекулярная масса дистиллята:
МP=32×0,882+18×(1–0,882)=30,35 кг/кмоль;
Средние массовые расходы жидкости в верхней и нижней частях колонны:
Средний мольный состав пара в верхней и нижней части колонны:
yср.в=(yP+yF)/2=70,4%;
yср.н=(yF+yw)/2=26,3%;
Средние мольные массы пара в верхней и нижней частях колонны:
Средние массовые потоки пара в верхней и нижней частях колонны:
В качестве контактных элементов используются ситчатые тарелки.
Рабочую скорость пара в интервале устойчивой работы ситчатых тарелок ректификационной колонны, можно определить из уравнения:
Р
асчет
скорости пара и диаметра в верхней части
колонны:
Плотность паров:
Средний массовый состав жидкости:
Плотность жидкости:
Из уравнения расхода определяем диаметр верхней части колонны:
Расчет скорости пара и диаметра в нижней части колонны:
Плотность паров:
Средний массовый состав жидкости:
Плотность жидкости:
Д
иаметр
нижней части колонны:
Таким образом рассчитанные диаметры верхней и нижней части колонн близки по величине, поэтому можно использовать колонну со стандартным внутренним диаметром: D=1000 мм.
Пересчитаем скорость в верхней и нижней частях колонны на реальный диаметр:
Техническая характеристика ситчатой тарелки типа ТС
Диаметр тарелки, мм |
1000 |
Диаметр отверстий, мм |
8 |
Шаг между отверстиями, мм |
16 |
Относительное свободное сечение тарелки, % |
13,6 |
Высота переливного порога, мм |
30 |
Периметр слива, м |
0,8 |
Рабочее сечение тарелки, м2 |
0,713 |
Масса тарелки, кг |
41,5 |
Определение числа тарелок и высоты колонны
Для определения высоты
колонны необходимо знать число
действительных тарелок. Число тарелок
рассчитывается графоаналитическим
методом – построением кинетической
линии при помощи КПД по Мэрфри, рассчитанным
через числа единиц переноса. КПД по
Мэрфри равен:
[1,
c.104, 3–46]
Для модели идеального смешения для жидкой фазы и идеального вытеснения для газовой фазы КПД по Мэрфри может быть рассчитан по уравнению:
[1, стр.239, 6–34]
где
– общее число единиц переноса. [1, cтр.239,
6–35]
Коэффициент массопередачи, отнесенный к единице рабочей площади тарелки, KYF определяется по уравнению аддитивности фазовых диффузионных сопротивлений:
,
[1, стр.239, 6–36]
где bXF и bYF – коэффициенты массоотдачи, отнесенные к единице рабочей площади тарелки, соответственно для жидкой и паровой фаз.
[1, стр.206, 5–40, 5–41]
U – плотность орошения, м3/(м2×с);
e – газосодержание барботажного слоя, м3/м3;
DX, DY – коэффициенты молекулярной диффузии распределяемого компонента соответственно в жидкости и паре, м2/с;
h0 –
высота светлого слоя жидкости на тарелке,
которая определяется из соотношения: [1,
стр.207, 5–42]
hп – высота газожидкостного барботажного слоя (пены) на тарелке, м.
Расчет высоты светлого слоя жидкости в верхней части колонны
В
ысоту
светлого слоя жидкости для ситчатых
тарелок определяют из уравнения:
Удельный расход жидкости на 1 м ширины сливной перегородки:
П
лотность
орошения U для ситчатых
тарелок:
Скорость жидкости на тарелке:
q=0,384/(814,60,8)=0,5910-3 м2/с;
wТ=1,121 м/с;
m=0,05-4,6hпер= -0,088
h0=0,0173 м;
Паросодержание барботажного слоя:
Fr – критерий Фруда;
F
r=7,40;
Расчет высоты светлого слоя жидкости в нижней части колонны
q=1,323/(917,10,8)=1,8010-3 м2/с;
wТ=1,173 м/с;
m=0,05-4,6hпер= -0,088
h0=0,0214 м;
Fr=6,55;
Расчет коэффициента диффузии распределяемого компонента в паре в верхней и нижней частях колонны.
Коэффициент диффузии в паровой фазе может быть вычислен по уравнению:
где:
Mмет , МН2О – молярные массы соответственно метанола и воды;
Р–абсолютное давление в колонне, Па;
Т–средняя температура в соответствующей части колонны, К;
vмет,vН2О – мольные объемы компонентов в жидком состоянии при нормальной температуре кипения (vмет=41,5, vН2О=28,2), см3/моль.
Для верхней части колонны:
Для нижней части колонны: