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

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

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

17. Сазыкин А.А. Термодинамический подход к разделению изотопов, в кн. Изотопы (Свойства. Получение. Применение.) / Под ред. В.Ю. Баранова. М: ИздАТ, 2000, С. 72108.

18. Wood H.G., Borisevich V.D., Sulaberidze G.A. On a criterion efficiency for multi-isotope mixtures separation // Separ. Sci. Technol. 1999. V. 34 (3). P. 343357.

19. Палкин В.А. Обобщение решения Смородинского для потенциала разделения многокомпонентной смеси изотопов // Атом-

ная энергия. 2003. Т. 95. № 5. С. 373382.

20. Палкин В.А., Фролов Е.С. Оптимизация каскадов по критерию максимума использования разделительной способности элементов // Физико-химические процессы при селекции атомов и молекул: сб. докладов XI международной научной конференции. Звенигород, Цнииатоминформ 2006. С. 2128.

21. Song T.M., Zeng S. On the optimity of separation cascade for a binary and a multi-component case // Proc. 9th Workshop on Separation Phenomena in Liquids and Gases. September 18-21, 2006, China, Beijing. P. 132143.

21

Часть 1

ТЕОРИЯ КАСКАДОВ ДЛЯ РАЗДЕЛЕНИЯ БИНАРНЫХ

СМЕСЕЙ

Внастоящем разделе описывается теория каскадов, состоящих из дискретных элементов и предназначенных для разделения бинарных изотопных смесей, прежде всего в приложении к однофазным молекулярно-кинетическим методам разделения. Рассматриваются как малые, так и большие эффекты разделения, реализуемые в разделительном элементе, в качестве которого используются такие аппараты как газодиффузионная машина, газовая центрифуга, масс-диффузионный насос Герца и др. Излагаемая теория основана на общих принципах умножения элементарного (однократного) эффекта разделения, которые позволяют использовать её, в частности, для расчета противоточных разделительных колонн (термодиффузионных, масс-диффузионных, дистилляционных, химобменных) и каскадов, составными частями которых являются такие колонны (каскады колонн с сокращением потоков).

1.1.Разделительный элемент, разделительная ступень. Основные параметры и уравнения [1-5]

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

зывается простым разделительным элемен-

том (рис.

1.1).

В этом

L , C

 

Lэ′ =θLэ, C

случае

в

разделитель-

э

 

 

 

 

 

 

 

ный

элемент

подают

 

 

 

 

поток разделяемой сме-

 

 

 

 

си LЭ с концентрацией

(мольной долей*) ценного (целевого изотопа) c , а из него отбирают два потока: «обогащенный», с более высокой, чем во входящем потоке, кон-

Lэ′′= (1θ)Lэ, C′′

Рис. 1.1. Схема разделительного элемента

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

23

центрацией c, и «обедненный», с более низкой концентрацией c′′. В отсутствие потерь рабочего вещества, которые могут быть вызваны, например, его частичным разложением в процессе разделе-

ния, обогащенный поток будет равен LЭ′ = θLЭ , а обедненный, со-

ответственно,

L′′ = (1 θ)L . Отношение

θ =

LЭ

называют ко-

 

 

Э

Э

 

LЭ

 

 

 

 

эффициентом деления потока, а расход

LЭ производительно-

стью элемента. В разделительную ступень параллельно соединяют элементы с одинаковыми вели-

 

 

 

 

 

 

 

 

чинами θ , c,

cи c′′ (рис.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.2). Она имеет суммарную

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

производительность L и выхо-

 

 

 

 

 

 

 

 

L, C

 

 

 

 

L′ =θL , C

дящие

потоки

L′ =θL

и

 

 

 

 

 

 

L′′ = (1 θ)L . Такое представ-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ление о работе элементов в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ступени предполагает, что все

 

 

 

 

 

 

 

 

они работают

в одинаковых

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

условиях и имеют одинаковые

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

характеристики.

 

 

 

 

 

L′′ = (1θ)L, C′′

 

 

 

 

 

Между

 

приращением

 

 

 

 

 

 

 

 

δ′ = c′−c

и

уменьшением

Рис. 1.2. Схема параллельного соеди-

δ′′ = c c′′

концентрации целе-

нения разделительных элементов в

вого компонента в выходящих

ступень

 

 

 

 

 

потоках

разделительного

эле-

 

 

 

 

 

 

 

 

мента существует определенная

связь. Для стационарного состояния эта связь может быть найдена из следующих уравнений материального баланса:

L = L′+ L′′,

 

 

 

 

 

(1.1)

Lc = L(c +δ) + L′′(c δ′′).

(1.2)

Из (1.1) и (1.2) следует, что

 

 

 

 

 

 

 

δ′′ =

L

δ′ =

 

 

θ

 

δ.

(1.3)

L′′

1

θ

 

 

 

 

 

 

24

 

 

 

 

Одной из основных разделительных характеристик ступени является полный коэффициент разделения, который для большинства однофазных методов разделения не зависит от состава изотопной смеси

q =

 

 

c

 

c′′

=

R

.

(1.4)

1

c

1 c′′

 

 

 

R′′

 

Здесь Rи R′′ - значения относительной концентрации ценного

(целевого) изотопа ( R = 1 c c ) в потоках обогащенной и обеднен-

ной фракции. Коэффициент разделения q характеризует эффект

разделения, достигаемый в одном элементе или ступени, и может зависеть от производительности ступени L и коэффициента деления потока θ :

q = q(L,θ),

(1.5)

где q(L,θ) – некоторая функция, которую определяют по резуль-

татам теоретических или экспериментальных исследований.

В соответствии со сказанным, к числу основных параметров ступени относятся восемь величин: L, L, L′′, c, c, c′′, q, θ , кото-

рые связаны независимыми соотношениями (1.1) –(1.5). Причём из перечисленных параметров свободными (независимыми) являются только три. Как правило, в качестве свободных параметров рассматривают величины L, c и θ . Однако в зависимости от рас-

сматриваемой задачи могут быть выбраны их различные комбинации.

Для удобства анализа эффекта разделения в элементе (ступени) могут быть введены дополнительные параметры и характеристики. Например, коэффициенты разделения по обогащённой α и обед-

нённой β фракциям, рассчитываемые как

 

 

 

 

α =

 

 

c

 

c

, β =

c

 

c′′

. (1.6)

1

c1c

1c 1c′′

 

 

 

Эти параметры характеризуют величину эффекта разделения в обогащённой и обедненной фракций по отношению к концентрации в потоке питания. Кроме того, для описания процесса разделе-

25

ния удобно пользоваться коэффициентами обогащения ε , ε, ε′′, равными

ε = q 1, ε′ =α 1, ε′′ =11 β .

(1.7)

Набор параметров ε , ε, ε′′

позволяет определить полное обо-

гащение ступени

 

 

δ = c′−c′′ = δ′+δ′′,

(1.8)

где δ′ = c′−c , δ′′ = c c′′. Если выразить cи c′′ из (1.6)

и ис-

пользовать (1.7), то в результате получим

 

δ ′ = εc(1 c)

, δ ′′ = ε′′c(1 c) .

(1.9)

1 + εc

1 ε′′c

 

Отсюда видно, что при заданных коэффициентах ε, ε′′

зави-

симости δ и δ′′ от c имеют максимумы, не совпадающие друг с другом. В случае «слабого обогащения» ( ε = q 1 <<1) наиболь-

шие значения δ и δ′′ достигаются в одной точке c = 0,5 , а фор-

мулы (1.9) существенно упрощаются:

 

δ′ = εc(1 c) , δ′′ = ε′′c(1 c) .

(1.10)

При подстановке (1.10) в (1.8) имеем

 

δ = εc(1 c) , ε′ = ε(1 θ) , ε′′ = θε .

(1.11)

Согласно (1.3) и (1.8) обогащения δ , δ , δ′′ связаны друг с

другом балансовыми соотношениями

 

δ′ = (1 θ)δ , δ′′ = θδ .

(1.12)

Следовательно, если коэффициенты разделения не зависят от параметров L и θ , можно путем уменьшения θ повысить концентрацию ценного компонента в обогащённом потоке c, произведя таким образом «перераспределение» полного обогащения δ в выходных потоках. При этом концентрация в обеднённом потоке c′′ будет приближаться к концентрации во входном потоке с . Очевидно, что возможность такого изменения обогащений δ и δ′′ связана с условиями сохранения материального баланса в разделительном элементе.

26

1.2.Разделительная способность (мощность). Работа разделения. Разделительный потенциал [1-7, 52]

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

~

Они предположили, что должна существовать функция U , с помощью которой можно охарактеризовать «ценность» изотопной смеси как в качественном, так и в количественном отношении, и которую можно представить в виде произведения экстенсивной величины – количества разделяемой смеси M – на интенсивную

величину – функцию V(c), зависящую только от концентрации

ценного изотопа и характеризующую качество смеси

 

~

(1.13)

U = MV (c) .

Функцию V(c) было предложено называть разделительным по-

~

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

Процесс разделения смеси в разделительной установке* для любого метода разделения схематически можно представить следующим образом. До начала процесса имелось некоторое количество исходной смеси F*(в ед. массы) с концентрацией ценного изотопа сF . С использованием введенных понятий изотопная ценность этой

~

 

* = F*V (cF ) . В

смеси будет определяться значением функции U

F

 

 

результате процесса разделения получают обогащённый продукт в количестве P* (в ед. массы) с концентрацией ценного изотопа сP и обеднённый продукт в количестве W* (в ед. массы) с концентраци-

~

 

* = P*V (cP )

ей сW. Изотопная ценность этих продуктов будет U

P

 

 

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

27

~

* =W *V (cW ) соответственно. В результате проведения

и U

W

 

 

 

 

~

 

~

процесса разделения функция U F *

изменится на величину

U :

~

~

~

~

 

U

=U P* +UW *

U F * =

(1.14)

 

 

 

.

= P*V (cP ) +W *V (cW ) F *V (cF )

 

Следовательно, «ценность» смеси будет выражаться соотношением

F*V (cF ) +

~

= P*V (cP ) +W *V (cW ) .

(1.15)

U

~

Из соотношения (1.15) следует, что величина U характеризует меру усилий, которую необходимо затратить на получение из первоначальной бинарной смеси изотопов два новых продукта – обогащенный и обедненный одним из изотопов этой смеси. При-

~

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

Отметим также, что величина работы разделения не дает ответ на вопрос о том, за какое время эта работа может быть выполнена на той или иной разделительной установке. Для получения ответа на этот вопрос необходимо знать разделительную мощность (способность) установки U , то есть работу разделения, выполняе-

мую установкой в единицу времени. Для перехода от

~

U

U к

достаточно в соотношении (1.14) вместо F * , P* ,W *

подставить

величины входящего (F) и выходящих (P и W) в установку потоков соответственно.

U = PV (cP ) +WV (cW ) FV (cF ) .

(1.16)

Для вычисления работы разделения и разделительной способности необходимо знать явный вид функции V(с). Пайерлс и Дирак

28

решили этот вопрос, рассматривая разделительную способность ступени (элемента), которую в соответствии с вышесказанным можно выразить следующим образом

δU = θLV (c) +(1−θ)LV (c′′) LV (c) .

(1.17)

В случае слабого разделения (q~1) функции V (с)

и V (с′′)

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

V (с) V (с) +

dV

δ′+

1 d 2V

 

(δ)2

,

 

(1.18)

 

 

 

 

 

 

2 2

 

 

 

 

 

 

V (с′′) V (с)

dV

δ′′+

1 d 2V

(δ′′)2 .

 

(1.19)

 

 

 

 

 

 

2 2

 

 

 

 

 

 

Подставив их в (1.17), находим

 

 

 

 

 

 

 

 

 

′′]dV

 

δU = [θL + (1 θ)L L]V (с) + [θLδ ′ − (1 θ)Lδ

+

+ 1 [θL(δ )2 + (1 θ)L(δ

′′)2 ]

 

 

 

 

 

 

 

(1.20)

 

d 2V

 

 

 

 

.

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

Из условий баланса (1.1), (1.2) и (1.12) коэффициенты при V (с)

и dV dс равны нулю. Поэтому с учётом (1.11) и (1.12) имеем

1

 

 

 

2

d 2V

 

с

2

(1 с)

2

.

(1.21)

δU = 2 θ(1 θ)Lε

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пайерлс и Дирак ввели условие, согласно которому разделительная способность (мощность) ступени или элемента определяется только его разделительными характеристиками L, θ, ε и не

должна зависеть от состава питающей его смеси. Это условие обосновывается следующим соображением. Если разделительная способность отдельной ячейки установки – разделительного элемента – не зависит от концентрации и все элементы работают в идентичных условиях, т.е. с одинаковыми L, θ, ε , то суммарная

разделительная способность элементов, составляющих эту установку, будет равна произведению Z δU эл , где δU эл – раздели-

тельная способность одного элемента, Z – число элементов в установке. Если процесс разделения организован так, что потери работы разделения (разделительной способности) отсутствуют, то ве-

29

личина Z δU эл будет равна разделительной способности всей ус-

тановки U , вычисляемой по формуле (1.16), т.е.

 

U = Z δU эл .

(1.22)

Откуда число разделительных элементов в каскаде будет определяться по формуле

Z =

U

.

(1.23)

 

 

δU эл

 

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

Условие независимости разделительной способности ступени (элемента) от концентрации приводит к следующему уравнению

δU =

1

θ(1θ)Lε 2

 

d 2V

c2 (1c)2

= const .

(1.24)

2

 

 

 

 

 

 

dc2

 

 

Если в уравнении (1.24) выбрать постоянную в виде

 

 

 

const =

1

θ(1 θ)Lε 2 ,

 

(1.25)

 

 

 

 

 

 

2

 

 

 

 

то для определения потенциала V(c) получается дифференциальное уравнение

 

d 2V

=

1

 

 

,

(1.26)

 

 

 

c2 (1c2 )

 

dc2

 

 

 

общее решение которого имеет вид

 

 

c

 

 

V (c) = (2c

1)ln

 

 

+ ac +b .

(1.27)

1

c

 

 

 

 

 

 

 

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

30