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

Борман Теория разделения изотопов 2007

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

 

c

P

c Rdi R

F

di

 

 

ciP = ciF

 

 

W

 

 

 

W

 

,

cF

cW

 

 

RWdi RP

di

 

 

c

 

c

 

 

 

 

Rdi R

 

di

(2.199)

 

 

P

 

 

 

 

 

 

 

ciW = ciF

 

 

W

 

 

 

F

 

P

 

.

c

P

c

F

 

 

 

Rdi R

 

di

 

 

 

 

 

 

 

 

W

 

F

 

 

Формулы (2.199) могут быть использованы для предварительной оценки содержания примесных изотопов в отбираемых из каскада потоках, исходя из заданных величин внешних потоков (подаваемых и отбираемых) и содержания в них основного (целевого) компонента, например, изотопа 235U при обогащении урана, а также содержания примесных изотопов

вподаваемом потоке.

2.3.4.4.Расчет R-каскада с дополнительным потоком отбора [3, 15]

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

Рассмотрим для определенности каскад, в котором имеет место несмешение относительных концентраций Rn,n+1 =

= cn cn+1 т.е. каскад, «идеальный» по компонентам n и n + 1.

Процесс разделения в таком каскаде задается следующими параметрами:

231

величинами внешних основных потоков (потоки отбора P, отвала W и поток питания F) и величинами концентраций целевого изотопа в этих потоках;

величинами внешних дополнительных потоков, подавае-

мых в каскад (или отбираемых из каскада) E j ( j =1, 2, ..., M)

здесь и далее M – число опорных точек, в которых производят дополнительный отбор или подачу продукта),

– положительный знак E j – для отбираемых потоков и от-

рицательный – для подаваемых;

– концентрации целевого изотопа cnE j в дополнительных

потоках;

– концентрации остальных компонентов ciE j (i n) в до-

полнительных потоках.

Целью расчета является определение значений концентра-

ций ciE

(i n) , числа ступеней между соседними опорными

 

j

точками, распределение потока разделяемой смеси (производительности) по длине каскада.

Для внешних потоков должны соблюдаться условия общего и покомпонентного баланса

M

 

 

 

 

E j

+ P +W F = 0 ,

(2.200)

j=1

 

 

 

M

 

 

 

 

E jciE

+ PciP +WciW FciF = 0,

(201)

j=1

j

 

 

 

 

 

 

i =1, 2, ...,

m.

 

Уравнения H-баланса

(2.141)

применительно

к случаю

k = n + 1 и при наличии M внешних дополнительных потоков записывают следующим образом

232

M

d

i + PciP (RnP,n+1)di +

 

E jciE j (RnE,nj +1)

 

 

j=1

 

 

 

+WciW (RnW,n+1)di FciF (RnF,n+1)di = 0,

 

i n, n +1;

 

di = Qi Qn .

(2.202)

Из условия несмешения относительных концентраций Rn,n+1 = cn cn+1 следует, что

Rn,n+1 = RnP,n+1 exp (Qnl),

(2.203)

где l – координата, отсчитываемая от точки основного отбора. Для участка каскада между двумя опорными точками (j-1)-й и j-й (отсчет опорных точек ведется от отбора к отва-

лу) уравнение (2.95) с учетом обозначения ϕi = exp (Qil) пе-

репишется в виде

 

 

 

1exp (Qis j )

 

 

 

L(s

 

)c

= L(0)c

 

+ 2E

c

,

(2.204)

 

 

Qi

 

 

j

iEJ

 

iE j 1

 

j1 iEJ 1

 

 

 

 

 

 

i =1, 2,

... , m,

 

 

 

 

 

где L(s j ), ciE j

– поток и концентрации компонентов в сече-

нии ввода (вывода) дополнительного потока Ej;

L(0), ciE

 

 

 

 

 

 

 

 

 

 

 

j1

поток и концентрации компонентов в сечении ввода (вывода) дополнительного потока E j1 ; s j – число ступеней каскада

между (j-1)-й и j-й опорными точками.

Уравнения (2.200) – (2.204) представляют собой полную систему для расчета указанных выше параметров каскада.

Конкретизируем эту систему для частного случая каскада с двумя потоками отбора (рис. 2.8).

Уравнения (2.200) – (2.202) для рассматриваемого случая будут иметь вид

1+

P1

+W

F

= 0 ,

(2.205)

P

P

 

P

 

 

 

233

Рис. 2.8. Принципиальная схема R-каскада с дополнительным отбором

ciP

+

P1

ciP

+W ciW

F

ciF = 0,

P

 

 

 

1

P

1

P

 

 

 

i =1, 2, ... , m,

ciP

(RnP,n+1)λi +

P1

ciP (RnP,1n+1)λi +

P

 

 

 

 

 

1

 

 

+W ciW (RnW,n+1)λi

F

ciF (RnF,n+1)λi = 0,

 

P

 

 

 

 

 

P

 

 

i n, n +1.

(2.206)

(2.207)

Считая, что поток L непрерывен по каскаду, а на его концах обращается в нуль, выражение, связывающее концентрации в основном и дополнительном потоках отбора из каскада, в соответствии с (2.204) запишем в виде:

L

P

c

iP

= 2Pc

iP

1exp (Qis1)

,

(2.208)

Qi

 

1

1

 

 

 

 

 

 

 

 

 

 

где LP1 – поток в сечении вывода дополнительного потока отбора.

234

j=1
ciF (RnW,n+1)di(RnF,n+1)di (RnW,n+1)di(RnP,n+1)di +
Pc ( ) ( )d + 1 iP RW diRP i 1 + 1+
ciP = m PciP n,n 1 n,n 1
cjF (RnW,n+1)d j(RnF,n+1)d j (RnW,n+1)d j(RnP,n+1)d j +j=1
Pc ( ) ( )d 1
+ 1 jP RW d jRP j .
1 + 1+ PcjP n,n 1 n,n 1
m
c j =1, находим

Учитывая условия несмешения по выбранной паре компонентов, находим:

RnP,n+1

RnP,1n+1

= exp (Qns1), RnP,n+1

RnF,n+1

= exp (QnsP ),

 

 

 

 

(2.209)

 

 

RnW,n+1 RnF,n+1 = exp (QnsW ).

Соотношения (2.209) позволяют исключить длины отборной и отвальной секций из числа неизвестных и свести задачу о расчете каскада к определению концентраций компонентов в отводимых потоках, поскольку величины внешних потоков с помощью уравнений (2.205) и (2.206) могут быть выражены через концентрации целевого компонента в них и отношение

P1/P:

F

 

c

nP

c

nW

 

 

cnP

cnW

 

P

 

 

 

=

 

 

 

 

1

 

1

,

(2.210)

P

cnF cnW

cnF cnW

P

 

 

 

 

 

 

 

 

W

=

F

 

P1

 

1.

 

 

 

(2.211)

 

 

 

 

P

P

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решая систему уравнений (2.206)–(2.207) относительно концентраций ciP и ciW с учетом уравнений (2.210) и (2.211)

и условия

(2.212)

235

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

di

 

 

 

P

 

 

di

+

 

 

 

 

 

 

 

 

 

ciF (Rn,n+1)

 

 

 

(Rn,n+1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pc

 

 

 

 

 

 

d

 

 

(

 

 

 

d

i

 

 

 

 

 

+

 

 

1 iP1

 

RF

 

 

 

 

i

R P1

 

 

{...}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PciP

( n,n+1)

 

 

 

 

n,n+1)

 

 

 

, (2.213)

ciW

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

F

 

 

 

d j

 

P

 

 

d j

+

 

 

 

 

 

 

c jF

(Rn,n+1)

 

 

 

(Rn,n+1)

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pc

jP1

 

 

 

 

 

d

j

 

 

 

 

 

 

d j

 

1

 

 

 

+

1

 

(

RF

 

 

 

R P1

 

 

 

 

{...}

 

 

 

 

 

 

 

 

 

 

+1

 

 

 

 

 

 

 

Pc jP

 

n,n+1)

 

 

 

 

(

 

n,n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где выражения в фигурных скобках {...} те же, что и в форму-

ле (2.212).

Из уравнения (2.208) с использованием соотношений (2.209) нетрудно получить

ciP

 

ciP

 

Q (RnP,n+1)di

(RnP,1n+1)di

1di

 

1

=

1

 

i

 

 

 

(RnP,n+1)

. (2.214)

ciP

ciP Qn

 

P

P

 

 

Rn,n+1

Rn,1n+1

 

 

Кроме того, суммируя уравнение (2.203) по всем компонентам и разделив (2.208) почленно на результат суммирования, с учетом (2.209) находим

 

 

P

P

di

 

P

P

dj

 

 

1(Rn,n+1

Rn,1n+1)

 

m

1(Rn,n+1

Rn,1n+1)

 

 

ci,

= ciP

 

 

 

c jP

 

 

 

. (2.215)

Qi

 

 

Qj

 

1

 

 

j=1

 

 

 

 

 

 

 

 

Таким образом, задавая значения концентраций целевого компонента в основном cnP и дополнительном потоках отбо-

ра cnP1 , в потоке отвала cnW (концентрации ciF считаем за-

данными), а также отношение потоков P1/P и, решая систему уравнений (2.212), (2.213), (2.215) для i = n, находим величи-

ны относительных концентраций RP

, RP1

, RW

, рас-

n,n+1

n,n+1

n,n+1

 

сматривая их как параметры задачи. После определения этих

236

параметров по соотношениям (2.212) – (2.215) рассчитывают концентрации остальных компонентов в отводимых из каскада потоках, а по соотношениям (2.209) длины отборной SP и отвальной SW секций каскада, а также длину участка каскада между двумя потоками отбора s1.

Алгоритм расчета каскада на заданные концентрации ключевого компонента cnP , cnP1 и cnW с использованием уравне-

ний (2.212) – (2.214) выглядит следующим образом.

1. Задают значения концентраций ключевого компонента в отводимых из каскада потоках cnP, cnP1 и cnW , а также отно-

шение потоков отбора P1/P.

2. Задают начальные приближения для относительных

концентраций RP

, RP1

, RW

.

n,n+1

n,n+1

n,n+1

 

3. С использованием уравнений (2.212) – (2.215) рассчи-

тывают текущие значения концентраций cnP( расч.), cnP ( расч.) и

1

cnW(расч.) .

4. Определяют величины невязок концентраций

δP

=cnP

cnP(расч.)

; δP

=cnP cnP (расч.); δW =cnW cnW(расч.) . (2.216)

 

 

 

1

1

1

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

иповторяют расчет по пунктам 3 и 4.

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

Таблица 2.4

Природный состав изотопов вольфрама

 

 

 

 

 

 

Номер компонента

1

2

3

4

5

 

 

 

 

 

 

Изотоп

180W

182W

183W

184W

186W

 

 

 

 

 

 

Концентрация, ciF , %

0,13

26,30

14,27

30,70

28,60

 

 

 

 

 

 

237

Ключевым компонентом выбран изотоп вольфрама с промежуточной атомной массой 183W. Нумерация компонентов проведена от «легкого» к «тяжелому» концу спектра масс изотопов, так что целевым является изотоп с номером n = 3

(cnF =14,3%). Решение системы уравнений выполнено мето-

дом Ньютона [16]. При решении системы (2.212)–(2.215) в качестве начальных приближений для относительных кон-

центраций

RP

, RP1

, RW

взяты их значения, получен-

 

n,n+1

n,n+1

n,n+1

 

ные из расчета каскада без дополнительного потока отбора

(P1

= 0), причем значениеRP1

было выбрано в области мак-

 

n,n+1

 

симума ключевого компонента.

При расчете каскада с дополнительным отбором на получение заданных концентраций c3P = 32%, c3P1 = 43% пара-

метром задачи являлось отношение потоков отбора P1/P, которое варьировалось в диапазоне от 0 до 1.

Результаты расчета в виде зависимостей W/P, ε0SP, ε0s1 и ε0SW ( ε0 – коэффициент обогащения, приходящийся на еди-

ницу разности массовых чисел) от величины параметра P1/P представлены в табл. 2.5.

Таблица 2.5

Параметры R-каскада с дополнительным потоком отбора для разделения изотопов вольфрама

 

 

c3P

= 32%, c3P = 43%, c3W =1, 2%

 

 

 

 

 

 

1

 

 

 

 

 

P1/P

0,042

0,060

0,080

0,100

0,120

0,160

0,200

0,204

0,208

 

 

 

 

 

 

 

 

 

 

W/P

1,449

1,487

1,531

1,575

1,619

1,707

1,795

1,805

1,814

 

 

 

 

 

 

 

 

 

 

ε0SP

5,764

6,008

6,368

6,776

7,251

8,590

12,332

13,654

22,879

 

 

 

 

 

 

 

 

 

 

ε0SW

2,331

1,732

1,551

1,439

1,358

1,246

1,168

1,162

1,155

 

 

 

 

 

 

 

 

 

 

ε0s1

5,912

5,917

5,925

5,933

5,941

5,955

5,967

5,971

5,973

 

 

 

 

 

 

 

 

 

 

238

Результаты расчета показывают возможность включением дополнительного отбора получить изотопный продукт с большим обогащением по ключевому компоненту, чем в концевом отборе каскада. Величина дополнительного отбора зависит от концентрации в нем ключевого изотопа, а также от состава исходной разделяемой смеси и концентрации ключевого компонента в основном отборе P. В рассмотренном случае величина дополнительного отбора достигала ~ 20% величины основного отбора при относительном увеличении величины концентрации целевого изотопа в нем по сравнению с c3P более чем на 30%. Анализ результатов показывает также,

что увеличение относительного отбора приводит к увеличению числа ступеней в каскаде и суммарного потока. Важно отметить, что отношение P2 P1 ограничено сверху. Сущест-

вование предельной величины P2 P1 объясняется тем, задан-

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

2.3.4.5.Аппроксимация каскадов непрерывного профиля прямоугольно-секционированным каскадом

Проектировочный расчет прямоугольно-секционирован- ного каскада по заданному отбору и ограничениям диапазона концентраций целевого изотопа в потоках отбора и отвала заключается в определении следующих его параметров: число, длина и потоки секций, отношение потоков отвала и отбора W/P. При этом, как правило, значения параметров должны соответствовать значениям, оптимальным по тому или иному критерию. Использование методов непосредственного расчета целевых каскадов для концентрирования заданного компонента представляется неэффективным, поскольку концентрации компонентов в отборе многокомпонентного каскада существенно зависят от распределения потоков, вследствие че-

239

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

До настоящего времени вопрос о модельном каскаде, обеспечивающем в случае «слабого» разделения условие

L = min , остаётся открытым. На практике в качестве мо-

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

Принципиальная возможность замены Q-каскада каскадом из секций постоянной ширины впервые была продемонстрирована в работе [12]. В [17] предлагается один из возможных критериев аппроксимации и разрабатывается методика расчета ПСК, распределение концентраций в котором наиболее близко к распределению концентраций в модельном Q-каска- де. Задача о замене участка Q-каскада, работающего в диапазоне концентраций от ci,нач до ci,кон ( (i =1,K, m , где m – число компонентов смеси), секцией постоянной ширины сформулирована следующим образом [17]: найти такие значения потока L и числа ступеней S в секции, при которых отклонения полученных в результате расчета концентраций в конце прямоугольной секции ci от заданных значений ci,кон

минимальны. При этом предполагается, что концентрации в начале секции и в отборе каскада ciP совпадают с соответст-

вующими концентрациями Q-каскада. В качестве критерия аппроксимации удобно использовать сумму относительных отклонений концентраций компонентов в конце секции

m

c

i,кон

c

 

 

φ =

 

i

 

.

(2.216)

 

c

 

 

i=1

 

i,кон

 

 

 

 

 

 

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

240