|
c |
P |
−c R−di − R |
F |
−di |
|
|
ciP = ciF |
|
|
W |
|
|
|
W |
|
, |
cF |
−cW |
|
|
RW−di − RP |
−di |
|
|
c |
|
−c |
|
|
|
|
R−di − R |
|
−di |
(2.199) |
|
|
P |
|
|
|
|
|
|
|
ciW = ciF |
|
|
W |
|
|
|
F |
|
P |
|
. |
c |
P |
−c |
F |
|
|
|
R−di − R |
|
−di |
|
|
|
|
|
|
|
|
W |
|
F |
|
|
Формулы (2.199) могут быть использованы для предварительной оценки содержания примесных изотопов в отбираемых из каскада потоках, исходя из заданных величин внешних потоков (подаваемых и отбираемых) и содержания в них основного (целевого) компонента, например, изотопа 235U при обогащении урана, а также содержания примесных изотопов
вподаваемом потоке.
2.3.4.4.Расчет R-каскада с дополнительным потоком отбора [3, 15]
Задачи разделения многокомпонентных смесей в R-кас- кадах, как правило, связаны с определенными практическими соображениями. Так, практическая польза от результатов расчета каскада для разделения многокомпонентных изотопных смесей с несколькими отборами связана с тем обстоятельством, что распределения концентраций компонентов с промежуточными массами имеют максимальные значения внутри каскада и, следовательно, существует возможность получить в дополнительном отборе концентрацию промежуточного компонента, большую, чем в потоке основного отбора на концах каскада.
Рассмотрим для определенности каскад, в котором имеет место несмешение относительных концентраций Rn,n+1 =
= cn
cn+1 т.е. каскад, «идеальный» по компонентам n и n + 1.
Процесс разделения в таком каскаде задается следующими параметрами:
–величинами внешних основных потоков (потоки отбора 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 внешних дополнительных потоков записывают следующим образом
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) пе-
репишется в виде |
|
|
|
1−exp (−Qis j ) |
|
|
|
L(s |
|
)c |
= L(0)c |
|
+ 2E |
c |
, |
(2.204) |
|
|
Qi |
|
|
j |
iEJ |
|
iE j −1 |
|
j−1 iEJ −1 |
|
|
|
|
|
|
i =1, 2, |
... , m, |
|
|
|
|
|
где L(s j ), ciE j |
– поток и концентрации компонентов в сече- |
нии ввода (вывода) дополнительного потока Ej; |
L(0), ciE |
– |
|
|
|
|
|
|
|
|
|
|
|
j−1 |
поток и концентрации компонентов в сечении ввода (вывода) дополнительного потока E j−1 ; 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 |
|
|
|
Рис. 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.
Считая, что поток L непрерывен по каскаду, а на его концах обращается в нуль, выражение, связывающее концентрации в основном и дополнительном потоках отбора из каскада, в соответствии с (2.204) запишем в виде:
|
L |
P |
c |
iP |
= 2Pc |
iP |
1−exp (−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 −di− RP 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 j− RP 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)
|
|
|
|
|
|
|
|
|
|
|
|
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 |
1−di |
|
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 |
|
сматривая их как параметры задачи. После определения этих
параметров по соотношениям (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 |
|
|
|
|
|
|
Ключевым компонентом выбран изотоп вольфрама с промежуточной атомной массой 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 |
|
|
|
|
|
|
|
|
|
|
Результаты расчета показывают возможность включением дополнительного отбора получить изотопный продукт с большим обогащением по ключевому компоненту, чем в концевом отборе каскада. Величина дополнительного отбора зависит от концентрации в нем ключевого изотопа, а также от состава исходной разделяемой смеси и концентрации ключевого компонента в основном отборе P. В рассмотренном случае величина дополнительного отбора достигала ~ 20% величины основного отбора при относительном увеличении величины концентрации целевого изотопа в нем по сравнению с c3P более чем на 30%. Анализ результатов показывает также,
что увеличение относительного отбора приводит к увеличению числа ступеней в каскаде и суммарного потока. Важно отметить, что отношение P2
P1 ограничено сверху. Сущест-
вование предельной величины P2
P1 объясняется тем, задан-
ная концентрация ключевого компонента не может быть достигнута ни в одном сечении каскада.
2.3.4.5.Аппроксимация каскадов непрерывного профиля прямоугольно-секционированным каскадом
Проектировочный расчет прямоугольно-секционирован- ного каскада по заданному отбору и ограничениям диапазона концентраций целевого изотопа в потоках отбора и отвала заключается в определении следующих его параметров: число, длина и потоки секций, отношение потоков отвала и отбора W/P. При этом, как правило, значения параметров должны соответствовать значениям, оптимальным по тому или иному критерию. Использование методов непосредственного расчета целевых каскадов для концентрирования заданного компонента представляется неэффективным, поскольку концентрации компонентов в отборе многокомпонентного каскада существенно зависят от распределения потоков, вследствие че-
го профиль целевого каскада предугадать невозможно. Поэтому проектировочные расчеты целесообразно проводить на основе модельного каскада непрерывного профиля.
До настоящего времени вопрос о модельном каскаде, обеспечивающем в случае «слабого» разделения условие
∑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,кон |
|
|
|
|
|
|
Расчет отдельных секций наиболее целесообразно проводить методом непосредственного интегрирования системы