Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

522103

.pdf
Скачиваний:
38
Добавлен:
12.05.2020
Размер:
1.09 Mб
Скачать

Y,Y*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Yн

 

 

 

 

 

 

 

 

 

 

C 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

Yк

О

X н

X 1

X к

X 2

X

X

 

Yк

 

X н

 

tк

 

н

 

 

 

X i

 

 

М

 

М

i

 

 

 

 

 

 

 

Y

X к

 

н

 

t

н

к

 

 

АВ - рабочая линия; BF - предельное положение рабочей линии, соответствующее теоретически минимальному расходу поглотителя; ОС -

равновесная линия при изменяющейся температуре жидкости ; ОС1, ОС2, ОС3 - равновесные линии для изотермической абсорбции при соответствующих температурах 1, 2, 3 .

Рисунок 2 - Схема адиабатической абсорбции

3 АБСОРБЦИЯ С РЕЦИРКУЛЯЦИЕЙ ЖИДКОСТИ

Схема абсорбции с рециркуляцией жидкости и циркуляционным отводом теплоты показана на рисунке 3 . Для приближенного метода расчета используются те же допущения, которые были сделаны в разделе 2.

Последовательность расчета материального и теплового балансов следующая:

Строят равновесную линию ОС и рассчитывают последовательно: Yк,G,M,Lмин,L,Xк как изложено в пунктах 2.1 -2.6 раздела 2.

Задаются коэффициентом рециркуляции жидкости b в пределах 0.2 - 0.4 и определяют концентрацию поглощаемого компонента на входе в

абсорбер

 

X '

по уравнению:

 

 

 

 

 

 

 

 

 

н

X н L+X к Lр

 

 

 

 

 

 

 

X

'

 

 

 

 

X

н

b X

к

,

(17)

 

 

L+L р

 

 

 

 

 

н’

 

 

 

 

 

1+b

 

 

 

11

где Lр - расход рециркулируемого поглотителя, кмоль/с; b= LLр .

Используя уравнение ( 14 ) с заменой X н на X'н’ рассчитывают

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

Уравнение рабочей линии АВ находят по формуле (13). При высоких температурах газовой смеси на входе в абсорбер

( tн к ) данный метод расчета не годится и необходимо искать

численное решение дифференциальных уравнений, описывающих процесс массо- и теплообмена между фазами [1,3].

Y,Y*

 

 

 

 

 

 

 

Yк

L', ,X'

L, í ,Xí

 

 

 

 

 

 

 

 

í

í

 

 

 

 

 

 

С

 

С

t

к

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

F

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Y н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lð, í ,Xê

Y

к

В

В

 

 

 

 

 

 

 

 

L , 'к , Xк

 

 

 

 

 

 

Yн

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

t н

 

 

 

О

 

 

 

 

 

 

 

 

 

 

X н

X н

X к

X *к

X * к

 

X

 

L', 'к

 

 

 

 

 

 

 

 

1 - абсорбер; 2 - холодильник; 3 - насос; АВ - рабочая линия без циркуляции жидкости; АВ - рабочая линия с ре-

циркуляцией жидкости; ОС - равновесная линия при изменяющейся температуре жидкости без рециркуляции последней; ОС - равновесная линия при изменяющейся температуре жидкости ' с рециркуляцией последней;

Рисунок 3 - Схема абсорбции с циркуляционным отводом теплоты (рециркуляция жидкости)

12

4 АБСОРБЦИЯ С ПРОМЕЖУТОЧНЫМ ОТВОДОМ ТЕПЛОТЫ

Схема противоточной абсорбции с промежуточным отводом теплоты для жидкой фазы показана на рис.4 . Для расчета предложенной схемы необходимо задаться числом секций n и концентрациями компонента в газовой фазе на выходе из секций Y1к Y2н,Y2к Y3н ( см. рисунок 4 ).

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

концентраций ΔYi Yн Yк по секциям, размеры их будут разными - по- n

следняя секция по высоте окажется наибольшей. Чтобы секции по высоте оказались одинаковыми, необходимо уменьшать значение ΔYi от первой

секции к последней. Предполагается, что каждая секция работает в адиабатических условиях и температуры жидкости после холодильников равны ее начальной температуре, то есть 3н 2н 1н . Расчет начинают

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

того, как будет найдена температура после последней секции (n = 3) ( см. рисунок 4 ), переходят к расчету следующих, который ведется аналогично.

13

Yк X,

Y, Y*

n=3

2

 

 

 

 

 

 

 

 

 

Y н

 

 

 

B

 

 

 

 

 

 

 

 

n=2

X,

2

Y=Y

 

 

 

E1

 

 

 

 

 

 

 

 

 

 

 

E2

 

С

 

 

Y=Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

E3

 

 

n=1

X,

 

Y к

 

 

D1

 

 

 

D2

 

 

1

 

 

D3

 

 

 

 

X,

 

 

 

 

 

 

Yн

 

О

X 3 н

Х 2 н

Х 1 н

Х к

Х

 

 

 

 

 

1 - секция ( n - номер секции); 2 - холодильники; ОС - равновесная линия при температуре 3’ ; D1E1, D2E2, D3E3, - равновесные линии при темпера-

турах жидкости в каждой секции: 1, 2, 3; АВ - рабочая линия. Рисунок 4 - Схема абсорбции с промежуточным охлаждением аб-

сорбента

5 ГИДРАВЛИЧЕСКИЙ РАСЧЕТ АБСОРБЦИОННЫХ НАСАДОЧНЫХ КОЛОНН

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

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

-определение диаметра колонны; -расчет плотности орошения и доли активной поверхности насадки;

-расчет гидравлического сопротивления слоя насадки.

5.1 ОПРЕДЕЛЕНИЕ СКОРОСТИ ГАЗА

В насадочных колоннах при противотоке газа и жидкости в зависимости от скоростей потоков наблюдаются три различных гидродинамических режима: пленочный, подвисания и эмульгирования ( захлебывания).

14

Зависимость обобщающая многочисленные экспериментальные данные для абсорбции по определению скорости газа в переходных точках (в частности, захлебывания) имеет следующий вид [3]:

w2

a

 

 

 

0.16

 

 

 

 

 

0.25

 

 

 

0.125

 

 

y

x

 

 

L

y

 

 

 

пр

 

 

 

 

 

 

 

 

 

 

 

lg

 

 

 

 

 

 

 

 

А В

 

 

 

 

 

 

 

,

(18)

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

g

x

 

 

G

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где wпр -предельная скорость газа, м/с; a- удельная поверхность

насадки, м23;

x , y - плотности жидкости и газа, кг/м3;

x - динамический коэффициент вязкости жидкости, мПа с;

g- ускорение свободного падения, м/с2;

- свободный объем насадки, м33;

L ,G - массовые расходы жидкости и газа, кг/с.

А, В - коэффициенты, значения которых для точки захлебывания приведены ниже [1]. Аналогичная формула, записанная в виде экспоненциальной зависимости [4], приведена в данном издании с ошибкой и определять по ней wпр нельзя.

Насадка

А

В

кольца Рашига (нава-

-0.073

1.75

лом)

 

 

седла, 25 мм (навалом)

-0.33

1.04

седла, 50мм (навалом)

-0.58

1.04

Затем рассчитывается рабочая скорость газа (w, м/с), по отноше-

нию:

 

w kwпр ,

(19)

где k - коэффициент, зависящий от выбранного режима работы колонны, значения k для wпр , рассчитанной в режиме подвисания приведе-

ны ниже [1]:

 

Процесс

k

абсорбция хорошо растворимых

0.75 - 0.85

газов

 

абсорбция плохо растворимых

0.45 - 0.7

газов

 

абсорбция в пенистых системах

0.25 - 0.4

Выбор рабочей скорости газа можно сделать и по представленной в [4] графической зависимости Эдулджи ( рисунок 5). Для этого

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

15

Fr w2

,

 

(20)

gd

 

 

 

Re усл

wd y

,

(21)

 

x

 

 

 

 

где d - номинальный размер насадки выбранного типа: для колец d равен наружному диаметру, м; для седел - среднему арифметическому трех основных размеров насадочного тела, м ( см. Приложение);

x - динамический коэффициент вязкости жидкости, Па с;y - плотность газа, кг/м3.

Далее необходимо рассчитать комплексы Y и X :

0.1

 

 

2

y

0.85

 

 

 

 

 

Y CFrRe усл

 

воды

 

 

 

,

 

 

x

 

0

 

 

 

 

 

 

Χ Vx , Vy

(22)

(23)

где С - константа, зависящая от типа насадки, для колец Рашига (навалом) С=1, для седел Берля С=0.471; - Инталокс С=0.37 ;

воды, 0 - плотности воды и воздуха при давлении 760 мм рт. ст. и 20 С, кг/м3;

Vx,Vy - объемные расходы жидкости и газа, м3/с.

Для определения рабочей скорости газа данным способом задаются коэффициентом k = 0.5 0.85, определяют скорость газа, находят на графике Эдулджи положение рабочей точки, и если она находится выше линии захлебывания или значение гидравлического сопротивления не соответствует значениям гидравлического сопротивления для данного типа абсорберов (см. соответствующий раздел), то изменяют скорость газа и повторяют расчет.

Для определения скорости газа в насадочных колоннах, работающих в режиме начала подвисания, приводится также зависимость [3]:

 

 

 

 

 

 

 

 

 

 

0.43

 

 

Reг 0.045Ar

0.57

G

,

(24)

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

где Reг

wdэ y

 

- критерий Рейнольдса для газа, условный;

 

y

 

 

 

 

 

 

 

 

 

 

 

dэ

4

- эквивалентный диаметр насадки, м;

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

d 3

 

y

 

x

 

Ar

э

 

 

- критерий Архимеда.

 

2

 

 

 

 

 

 

 

 

 

 

y

 

 

 

Однако значения скорости, определенные этим способом, значительно ( до 2 раз) отличаются от полученных по формулам (18 - 19) и по графику Эдулджи и не подтверждаются опытными данными.

Y

 

 

 

101

 

 

 

8

 

 

 

6

1200

 

 

 

4

800

 

 

Захлебывание

 

 

 

2

600

 

 

 

 

Подвисание

 

300

0

 

 

 

2

 

100

 

 

 

8

100

Па/м

 

6

 

 

 

4

 

 

 

2

 

 

 

10-1

4

6

8 10-3

2

4

6

8 10-2

2

4

6

8 10-1

X

Рисунок 5 - Графическая зависимость Эдулджи для определения гидродинамического режима работы насадочной колонны и сопротивле-

ния насадки

17

5.2 ОПРЕДЕЛЕНИЕ ДИАМЕТРА КОЛОННЫ Диаметр колонны D (м) может быть рассчитан как:

D

Vy

,

(25)

0.785w

где V y - объемный расход газа, м3/с;

w - рабочая скорость газа, м/с.

По расчетному диаметру колонны выбирается ближайший стандартный аппарат, по диаметру которого пересчитывается рабочая скорость.

5.3 РАСЧЕТ ПЛОТНОСТИ ОРОШЕНИЯ И ДОЛИ АКТИВНОЙ ПОВЕРХНОСТИ НАСАДКИ

Плотность орошения (U, м3/(м2 с)) определяется по зависимости:

 

 

 

 

 

 

 

U

L

 

.

(26)

 

 

 

 

 

 

 

 

 

0.785D

2r

 

 

 

 

 

x

 

Желательно, чтобы полученное значение находилось в пределах 0.002-0.005 м3/(м2 с). Большие значения характерны для нерегулярных насадок. При малой плотности орошения аппарат работает неэффективно из-за плохой смоченности насадки. В этом случае следует использовать более крупную насадку или отказаться от применения насадочной колонны.

Долю активной поверхности насадки ( ya ), участвующую в массопереносе, находят по соотношению:

a

U

 

 

 

,

 

( 27 )

p+qU a

 

где p,q - постоянные, зависящие от типа и размера насадки.

Значения p,q для нерегулярных насадок следующие:

 

Тип насадки

p 105

q

кольца Рашига, 25 мм

1.019

0.0086

кольца Рашига, 35 мм

0.843

0.0103

кольца Рашига, 50 мм

0.667

0.012

18

5.4 РАСЧЕТ ГИДРАВЛИЧЕСКОГО СОПРОТИВЛЕНИЯ СЛОЯ НАСАДКИ

Общее гидравлическое сопротивление колонны находят после определения высоты колонны, но на данном этапе проектирования необходимо убедиться в том, что выбранные ( рассчитанные) рабочая скорость газа, диаметр колонны и плотность орошения обеспечивают необходимый режим работы аппарата. Для этого определяют гидравлическое сопротивление одного метра слоя насадки по графику Эдулджи ( см. выше ). Рекомендуются следующие значения сопротивлений 1м насадки : для абсорберов, работающих под атмосферным или повышенном давлении 200-400 Па/м, для вакуумных абсорберов - 80-400 Па/м Если полученное значение значительно больше рекомендованного для данного типа абсорбера, то увеличивают диаметр аппарата, уменьшая тем самым скорость газа и плотность орошения, или применяют насадку с большим свободным объемом. В работе [1] описан метод определения соотношения сопротивлений орошаемой и сухой насадок из графической зависимости указанных сопротивлений от комплекса К, включающего числа Рейнольдса по жидкости и газу и число Фруда. Однако в диапазон данных, для которых построена зависимость, практически не входят режимы работы реальных насадочных колонн, поэтому использовать этот график, как правило, невозможно.

6 РАСЧЕТ ВЫСОТЫ КОЛОННЫ

6.1 ОПРЕДЕЛЕНИЕ КИНЕТИЧЕСКИХ ПАРАМЕТРОВ И ВЫСОТЫ НАСАДКИ

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

6.1.1 РАСЧЕТ ПО УРАВНЕНИЮ МАССОПРЕДАЧИ

M K y F

Y ср K y SHa

а

ΔYср ,

(28)

 

 

 

 

где М - количество поглощаемого компонента в единицу времени, кмоль/с;

K y - коэффициент массопередачи, отнесенный к средней движу-

щей силе в газовой (паровой) фазе, кмоль/(м2 с ед. средней движущей силы);

F - поверхность контакта фаз, м2;

19

ΔY- средняя движущая сила процесса по газовой (паровой)

кмоль/кмоль;

S - площадь поперечного сечения колонны, м2; Н - высота слоя насадки, м.

Или используют относительные массовые концентрации:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M=K y F YK y SHa

а

Y

,

 

 

 

 

 

 

 

 

 

 

 

фазе,

(29)

где M - количество поглощаемого компонента в единицу времени,

кг/с;

K y - коэффициент массопередачи, отнесенный к средней движу-

щей силе в газовой (паровой) фазе, кг/(м2 с ед. средней движущей силы); Y- средняя движущая сила процесса по газовой (паровой) фазе,

кг/кг.

Средняя движущая сила определяется в зависимости от кривизны равновесной линии.

Если равновесная линия на отрезке, ограниченном начальной и конечной концентрацией поглощаемого компонента в жидкости, прямая, то:

 

 

ΔY -

 

Yм

 

 

 

 

 

 

 

 

 

 

ΔYср=

 

 

 

 

б

 

 

 

 

,

 

 

 

 

 

(30)

 

 

 

ln

ΔYб

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ΔYм

 

 

 

 

 

 

 

 

 

 

 

 

где ΔY

 

 

Y

 

Y*

 

при

X=X

к

, кмоль/кмоль;

б

 

 

 

н

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ΔY Y

 

Y*

при X=X

н

, кмоль/кмоль.

м

 

 

к

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

Если линия равновесия криволинейная в указанном выше интерва-

ле, то:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ΔYср

 

 

Yн

Yк

 

 

 

.

 

 

 

 

 

(31)

Yн

dY

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y Y*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Такие же уравнения можно записать и для случая, когда равновесная линия задана в относительных массовых концентрациях. Уравнения, аналогичные уравнениям ( 28 - 30 ), будут и в том случае, если вести расчет по средней движущей силе и коэффициенту массопередачи, рассчитанным по жидкой фазе.

Коэффициенты массопередачи рассчитываются по уравнениям:

1

 

1

 

m yx

,

(32)

K y

y

x

 

 

 

 

20