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

Борман Теория каскадов для разделения бинарных 2011

.pdf
Скачиваний:
179
Добавлен:
16.08.2013
Размер:
6.3 Mб
Скачать

ние концентрации ценного компонента на одной ступени δs можно

заменить производной dc

, а уравнение (1.234) переписать в виде

ds

 

 

 

 

 

 

dc = εt c(1c)

2P (cP c) ,

(1.235)

ds

 

 

 

 

L

 

где концентрация c = c(s)

 

является непрерывной функцией номера

теоретической ступени s, а

 

L

= L = L .

 

 

 

 

 

 

2

1

2

 

 

 

 

 

 

Уравнение (1.235) совпадает с уравнением (1.85) при L=const, т.е. с уравнением прямоугольного каскада. Таким образом, математическое описание процесса разделения в установке с «внутренним» каскадированием (т.е. в противоточной колонне) совпадает с описанием процесса разделения в прямоугольном каскаде из однофазных дискретных ступеней.

Если через h обозначить величину ВЭТС – высоту, эквивалентную теоретической ступени, через z продольную координату в еди-

ницах длины (ds = h1dz) и считать, что величина h от z не зависит, то уравнение (1.235) примет вид

h dc

=εt c(1c)

2P

(cP c) .

(1.236)

dz

 

L

 

 

Величину h для каждого метода разделения изотопов находят, решая задачу конвективной диффузии в противоточной колонне [7]. Анализ показывает, что ВЭТС определяется гидродинамическим режимом течения в колонне, физико-химическими свойствами разделяемой смеси и геометрическими параметрами колонны. Следует отметить, что при использовании 2-фазных колонн (дистилляция, химический изотопный обмен) на их торцах устанавливают специальные устройства, предназначенные для создания противоточного движения фаз путем перевода разделяемой изотопной смеси из одной фазы в другую. Эти устройства называют узлами обращения потоков (УОП) (рис. 1.19а, б).

Заметим, что изменение потоков от ступени к ступени (секционирование) в одиночной колонне осуществить невозможно. Поэтому в этом случае соединяют последовательно ряд колонн (секций),

91

в каждой из которых потоки постоянны по длине, но меняются от колонны к колонне. В случае 2-фазных методов каждая колонна (секция) имеет УОП для того, чтобы вернуть в нее ту часть потока, которая не поступает в последующую колонну (секцию).

Рис. 1.19а. РазделиРис. 1.19б. Схема прямоугольно-секционированного тельная колонна в 2- каскада из 2-фазных колонн (секций), предназначенфазных методах разных для концентрирования легкого (в случае дистилделения с устройстляции более летучего) компонента. Секцией обедневами обращения пония является нижняя часть первой колонны, концентоков трирование происходит в верхней части первой и всех

последующих колонн (секций)

Возможность сокращения потоков позволяет уменьшить объем аппаратуры, а также количество разделяемого вещества, находящегося в этом объеме. Такая установка по сути дела представляет собой прямоугольно-секционированный каскад, процесс разделения в котором описан в разделе 1.9.1. Таким образом, теории каскадов для разделения бинарных изотопных смесей оказывается весьма универсальной. То есть, с точки зрения теории, не принципиально является ли каскад установкой, состоящей из последовательно соединенных разделительных ступеней, или одиночной противоточной колонной, а также состоит ли ступень из одного или несколь-

92

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

1.9.3.Оптимизация ПК и ПСК в случае «слабого обогащения»

[3, 5, 7]

Целью оптимизации при заданных внешних параметрах каскада является определение таких его внутренних параметров, при которых КПД формы каскада имеет максимальное значение.

При оптимизации каскадов, помимо КПД формы, используют понятие КПД ступени (или локального КПД). Определить его можно как отношение приращения ценности, создаваемого ступенью, к

ее разделительной способности (мощности):

 

 

 

 

 

 

 

 

 

ηст =

P[V

(cP ,c)V

(cP ,c +

c)]

,

 

(1.237)

 

 

 

 

 

 

Lε 2

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c)(12c)

 

 

 

 

 

 

c

P

 

 

 

 

 

c

 

(c

P

,(1.238)

V (cP

,c) = P (2cP

1) ln

 

 

 

 

(2cP

1) ln

 

 

 

+

 

 

 

 

 

1

cP

1

c

 

 

c(1c)

 

 

 

 

 

 

 

 

 

 

 

а c - приращение концентрации на ступени.

Разлагая разность в квадратных скобках (1.237) в ряд вблизи c , ограничиваясь первым членом разложения и учитывая, что прира-

щение концентрации на одной ступени c =

dc

s = dc

(по-

скольку s =1), найдем:

ds

ds

 

 

 

 

 

 

ηст

= −

4P

 

 

V (cp ,c)

 

dc

 

 

 

(1.239)

Поскольку

Lε 2

c

 

 

ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εc(1c)

 

 

L*

 

 

dc

= εc(1

c)

2P(cP c)

=

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ds

 

 

 

L

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

93

 

 

 

 

 

 

 

 

 

 

 

 

*

4P(cP c)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где L =

εc(1c)

, а производная разделительного потенциала по

концентрации с равна V (cP ,c)

= −

cP c

 

, то КПД ступени

c2 (1c)2

принимает следующий вид:

c

 

 

 

 

 

 

 

L*

 

 

 

L*

 

 

 

 

 

 

 

 

 

 

 

 

 

ηст

=

 

 

2

 

 

 

.

 

(1.240)

 

 

 

 

L

 

 

 

 

 

L

 

 

 

 

 

 

Нетрудно видеть, что

ηст = 0

при

 

L*

 

= 0 , т.е. при очень

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

больших потоках ( L ∞) и при L* = 2 , т.е. при значении потока,

L

при котором dcds =0.

Решение задачи об определении оптимальных условий разделения начнем с простейшего случая – прямоугольного каскада (ПК).

Математический анализ показывает, что чем «короче» ПК, тем выше его КПД формы. Выясним, величину максимального КПД, которую может иметь ПК. Для «короткого» ПК можно принять c(1c) const = cF (1cF ) , и, следовательно, для обогатительной

части идеального каскада

*

 

4P(cP c)

 

L

=

 

(1.241)

εcF (1cF )

для обеднительной части идеального каскада

*

 

4W (c cW )

 

L

=

 

.

(1.242)

εcF (1cF )

Соответственно, вычисляя суммарные потоки отборной и отвальной частей «короткого» идеального каскада, имеем:

94

 

 

 

 

 

 

 

 

 

 

SP

*

 

 

CP

 

L* (c)dc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

IP =

L ds

=

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dc

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

CF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

(1.243)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8P

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

P (cP

c)dc =

 

 

(L* (cF ))2

 

 

 

 

 

 

 

ε

2

2

(1cF )

2

 

 

 

4P

 

 

 

 

 

 

 

 

 

 

 

 

cF

 

 

 

CF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

SW

*

 

 

CF

 

 

L* (c)dc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

IW =

L ds

=

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dc

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

CW

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

(1.244)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8W

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

F (c cW )dc =

 

 

 

 

(L* (cF ))2

 

 

 

 

 

 

 

ε

 

2

2

(1cF )

2

 

 

 

4W

 

 

 

 

где

 

 

 

 

cF

 

 

 

CW

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dc

*

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

4P(cP

cF )

 

4W (cF

cW )

 

 

 

 

=

 

 

 

εcF (1 cF ) , L (cF ) =

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

.

 

2

εcF (1cF )

εcF

(1

cF )

ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Суммарный поток «короткого» идеального каскада составит

 

 

 

 

 

 

 

I = L*

= IP + IW =

 

1

 

(

 

1

 

+

 

1

)(L* (cF ))2

 

(1.245)

 

 

 

 

 

 

4

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

Решения (1.218) и (1.219), перепишем в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

SP =

 

 

 

L

 

ln

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

,

 

 

 

 

 

 

 

(1.246)

 

 

 

 

 

 

 

 

 

 

 

2P

 

 

 

 

 

 

 

 

 

*

(cF )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

=

 

 

 

 

ln

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

.

 

 

 

(1.247)

 

 

 

 

 

 

 

 

 

 

 

 

 

2W

 

 

 

 

 

 

 

L* (cF )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

Тогда величина КПД формы «короткого» ПК с учетом (1.245)-

(1.247) запишется в виде

 

 

 

 

 

η =

L*

=

ς 2

,

(1.248)

L SP + L SW

2ln

2

 

 

 

 

 

 

 

2 ς

 

 

 

95

 

 

 

 

 

где

ς = L* (cF ) .

L

а)

б)

Рис. 1.20. Определение оптимальных параметров каскада при следующих ис-

ходных данных: cp=3,5%, cF=0,719%, cw=0,2%, F/P=6,353, W/P=5,353 [4]

96

Величина КПД имеет максимум в зависимости от безразмерного параметра ς , а именно η = ηmax = 0,814 (т.е. 81,4%) при

ςопт =1,43. Величина 81,4% представляет верхний предел КПД, ко-

торый может быть достигнут в ПК.

В общем случае оптимизацию ПК целесообразно проводить следующим образом. По формулам (1.217) с учетом (1.77) и (1.78)

рассчитываем зависимость ε(SP + SW ) от ε2PL .

Далее, с учетом этих результатов и построив зависимость КПД формы от той же переменной, находим оптимальные значения ε2Lp

и ε(SP + SW ) (рис. 1.20, а и б).

На рис. 1.21 для заданных значений cp , cF и cw представлены

распределения потоков питания ступеней в идеальном каскаде и аппроксимирующем его оптимизированном прямоугольном каскаде.

Рис. 1.21. Профили идеального и оптимизированного прямоугольного каскада (потоки выражены в относительных единицах). Исходные условия: cp=3,5%, cF=0,719%, cw=0,2%, F/P=6,353, W/P=5,353 [4]

97

Здесь же приведено распределение Lмин (s), соответствующее

минимальному потоку через каждую ступень. Очевидно, что величина потока прямоугольного каскада (в том числе и оптимизиро-

ванного) должна превышать величину потока Lмин(cF ) , равную

L

(c

F

) =

2P(cp cF )

=

2W (c

F

c

w

)

.

(1.249)

 

 

 

 

 

 

 

мин

 

 

εcF

(1cF )

 

εcF (1cF )

 

 

 

 

 

 

 

 

 

Рис. 1.22. Распределение КПД ступени (локального КПД) по длине прямоугольного каскада. Исходные данные: cp=3,5%, cF=0,719%, cw=0,2%,

F/P=6,353, W/P=5,353, ε2PL = 5, 45 [4]

На рис. 1.22 представлено распределение локального КПД по длине оптимизированного прямоугольного каскада. На двух ступенях этого каскада величина локального КПД достигает 100% (на входе в каждую из этих ступеней концентрации c таковы, что

L = L* (c) . Монотонное уменьшение ηст у концов каскада, очевидно, связано с тем, что на концевых ступенях каскада ( c cp , c cw ) и, соответственно, L* 0 (см. формулу 1.240).

98

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

Таблица 1.4

Параметры ПК, оптимизированных по величине КПД формы (разделяемая смесь 235UF6/238UF6) [4]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cP ,%

 

3,5

 

5,0

 

20

 

20

 

93,5

 

93,5

 

93,5

 

 

 

 

 

 

 

 

 

 

 

 

 

cF ,%

 

0,7

 

0,7

 

0,7

 

3,5

 

0,7

 

3,5

 

5,0

 

 

 

 

 

 

 

 

 

 

 

 

 

cW ,%

 

0,2

 

0,2

 

0,2

 

0,7

 

0,2

 

0,7

 

0,7

 

 

 

 

 

 

 

 

 

 

 

 

 

U / P

 

5,356

 

8,765

 

45,397

 

8,261

 

235,285

 

26,583

 

45,286

 

 

 

 

 

 

 

 

 

 

 

 

 

F / P

 

6,353

 

9,241

 

38,121

 

6,934

 

179,630

 

33,367

 

21,675

 

 

 

 

 

 

 

 

 

 

 

 

 

εL / 2P

 

5,443

 

8,241

 

35,094

 

6,869

 

155,644

 

33,468

 

24,091

 

 

 

 

 

 

 

 

 

 

 

 

 

εSP

 

2,438

 

2,948

 

5,092

 

2,860

 

9,432

 

7,716

 

7,386

 

 

 

 

 

 

 

 

 

 

 

 

 

εSW

 

2,672

 

2,628

 

2,8562

 

3,268

 

3,334

 

3,894

 

4,688

 

 

 

 

 

 

 

 

 

 

 

 

 

ε(SP + SW )

 

5,010

 

5,576

 

7,764

 

6,128

 

12,766

 

11,610

 

11,966

 

 

 

 

 

 

 

 

 

 

 

 

 

η,%

 

78,6

 

76,3

 

66,7

 

78,5

 

47,4

 

58,2

 

62,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача оптимизации ПСК для заданных величин P , cP , cF , cW , ε , числа прямоугольных секций в обогатительной и обедни-

тельной частях каскада сводится к определению оптимальных параметров, соответствующих максимальному значению КПД формы, а именно потоков во всех секциях, концентраций на границах («стыках») секций и соответственно чисел ступеней во всех секциях.

Оптимальные значения «стыковых» концентраций находят следующим образом. Для ПСК, состоящего из k прямоугольных участков, сумму потоков можно записать виде:

99

 

 

k

k

ci

 

 

 

 

 

 

K

dc

 

 

 

 

Li Si

= Li

 

,

(1.250)

 

 

 

dc ds

 

 

i =1

i =1

c

i

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

где ci

, ci

– концентрация ценного изотопа в начале и конце i

H

k

 

 

 

 

 

 

 

секции, соответственно. Для нахождения оптимума (минимума) соотношения (1.250) продифференцируем его по cKi , учитывая, что cKi входит в уравнение дважды: как верхний предел интеграла для

секции i и как нижний – для секции i +1 (так как cHi +1 = cki ). По-

скольку производная интеграла по нижнему пределу отрицательна, условие оптимума (минимума) будет иметь вид:

 

 

 

 

 

 

 

k

 

L

 

 

L

 

 

 

 

 

 

 

 

 

 

 

[Li Si ] =

 

i

 

i +1

= 0

(1.251)

 

 

 

 

 

ci

dc i

dc i +1

 

 

 

 

 

i =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

ds

 

ds

 

 

dc i

dc

i +1

 

 

 

 

 

 

 

 

 

где

 

 

,

 

 

 

– значения производных в конце i -й секции и

 

 

ds K

ds

H

 

 

 

 

 

 

 

 

 

в начале i +1-й секции соответственно.

Учитывая соотношения (1.85) и (1.144), условие (1.251) можно преобразовать к виду

 

ηCTi

,K =ηCTi+1,H

(1.252)

где ηCTi

,K ,ηCTi+1,H – локальные КПД в конце секции i

и вначале сек-

ции i +1 соответственно. Полученное условие (1.252) означает, что при оптимальном соединении секций в ПСК локальные КПД в точке перехода от предыдущей секции к последующей должны быть равны.

Подставляя в (1.252) значенияηCTi

,K ,ηCTi+1,H ,

после несложных

преобразований получим:

 

2Li

 

 

 

*

 

 

 

 

L

=

 

 

,

(1.253)

σi +

1

i

 

 

 

где

100