Моделирование нестационарных процессов в центробежных каскадах Текст научной статьи по специальности «Математика»

УДК 661.87.519

МОДЕЛИРОВАНИЕ НЕСТАЦИОНАРНЫХ ПРОЦЕССОВ В ЦЕНТРОБЕЖНЫХ КАСКАДАХ

В.А. Власов, В.Г. Бутов, Д.Н. Голдобин, А.А. Орлов, С.Н. Тимченко

Томский политехнический университет E-mail: richi@inet.tsu.ru

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

Введение

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

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

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

1.1. Модель противоточного симметричного каскада

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

ной изотопной смеси. Пусть в некоторую промежуточную точку каскада подается поток питания Р, а из концов каскада отбираются потоки Р (отбор) и Ж (отвал). Схематично такой каскад представлен на рис. 1.

Величины Д С, Ж, См, Р, Ср являются внешними параметрами каскада, а Ы, О+, С+, О-, С,-(/=1,п) - внутренними. Эти параметры получаются из внешних согласно их определениям и известным характеристикам ступеней. Величины Р, Р и Ж являются функциями от времени и связаны уравнением баланса:

^ = Р + Ж,

^ • СЕ = Р • Ср + Ж • Сж .

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

Система уравнений, описывающая нестационарные гидравлические процессы строится следующим образом [1]:

Газосодержание каскада, состоящего из п ступеней, разбивается на п объемов газа.

Газосодержание в каждой ступени определяется массой газа, вращающегося в роторах центрифуг [2]

Mo =pV

1 -

A =

M Q2 ra 2RT

(1)

W, С

F, C

P,Cp ►

закр

Рис. 1. Противоточный симметричный каскад: I - номер ступени; N - количество ступеней; )М, С„ - поток и концентрация отвала каскада; Р, Ср - поток и концентрация отбора каскада; F, СI - поток и концентрация внешнего питания, подаваемого в 1-ую ступень; G+, С+ -поток и концентрация отбора I-ой ступени; G-, С--поток и концентрация отвала I-ой ступени; GI, С - поток и концентрация питания 1-ой ступени

где М0 - газосодержание ротора газовой центрифуги (ГЦ) с объемом V=nraH, радиусом ra и высотой H определено для идеального газа с молярной массой M, равновесной плотностью p, давлением p, при температуре T0, вращающегося как единое целое вокруг оси z с угловой скоростью Q=Vw|ra; Vw -окружная скорость ротора. Газосодержание M¡(t) i-ой секции рассчитывается как произведение количества входящих в нее центрифуг на M0

Mt = Nt ■ M0, i = 1...n.

Перепишем это соотношение с учетом ур. (2), заменив p=Mp|RT0:

M E E N¡ VM 1 - e-A

M = Ep, E = —-----------------. (2)

i RT0 A

В общем случае входящими в i-ый объем потоками являются поток отбора предыдущей ступеней G+-1; поток отвала следующей ступеней G+1; внешний поток питания данной ступени F, являющийся заданной функцией от времени t. Кроме того, в случае закрытия соседних ступеней, в i-ый объем поступают потоки отбора G+-2 и отвала G-2 с ближа-щей работающей ступени.

Из объема ступени выходят поток отбора G+ и поток отвала G-

1.2. Основные уравнения нестационарной гидравлики

Это уравнения баланса вещества в выделенных объемах. Изменение газосодержания в выделенном объеме определяется разностью входящих и выходящих потоков.

Для i-ой ступени можно записать следующие уравнения нестационарной гидравлики:

dMi- = 8fF+s;_iG+_i +8i+ 2 +

dt

+8-+G+l +8- G+ 2 - G+- G

(3)

где SF, 8+_ъ S+_2, 5Í+1, S+2 - признаки (0 или 1) наличия потока в секции, t - время.

Из (1) и (2) следует, что изменение величины E во времени зависит от Q(t), поэтому в общем случае газосодержание газовой центрифуги зависит от частоты ее вращения:

dM ^ dp dE

----= E— + p —.

dt dt dt

Начальное условие для уравнения (3) получается из ур. (2) при Q(0)=Q0, A(0)=A и имеет вид

NMV 1 - e 0 M, (0) = -^------------------;--------pi (0),

RT0

(4)

где р(0) - начальное давление в г-ой ступени.

1.3. Уравнения для потоков

Поток питания отдельной ГЦ проходит через критическую шайбу. Выражение для потока питания ГЦ г-ой ступени можно записать в виде:

С = N • к • р,

где N - количество ГЦ в г-ой ступени; К - коэффициент расхода; Р1 - давление питания на входе в ГЦ.

Исходя из ламинарного характера течения газа в отборнике тяжелой фракции ГЦ, для потока отвала ГЦ г-ой секции О- имеем [3]:

с-д-

Ni

= PW - (P- )2

где С,- - коэффициент гидравлического сопротивления.

Поток отбора ГЦ г-ой ступени О+ определяем с помощью формализованной универсальной гидравлической характеристики ГЦ, которая представляется в виде:

(

G+ = N

bl + b2G + b3.

i i Ni i

А

V Ni У

Gi

+b4 P-+ b5. (Pr )2 + b6.-i-P-

i i i V i ' i jN i

где Ь1-Ь61 - коэффициенты гидравлической характеристики, зависящие в общем случае от типа ГЦ, температуры ротора, концентрации легких примесей, частоты вращения О и т. п.

1.4. Основные уравнения нестационарного переноса

Перенос легкого изотопа от отвала в сторону отбора в обогатительной и регенеративной частях каскада соответственно записывается в следующем виде [2, 4, 5]:

+ G (1 ч 9c G

т, =tc + Gsc (1 - c )-----------------,

1 on 2

r¡ = -r’c + Gsc (1 -c) -—G, 9п 2

(5)

(б)

где е - коэффициент полного обогащения.

Здесь в силу плавного изменения обогащения в каскаде номер ступени можно считать переменной величиной, тогда концентрация е=е(п,1), О=О(п,1), т=т(п,/) и т'=т'(п,/).

После подстановки т и т' из (5) и (6) в уравнение переходного процесса в центробежном каскаде д(Мс) дт, д(Мс) дт/

д( дп’ д( дп

для обогатительной части получим

D(Mc) _ G д 2c д

Dt

= — т-у - —[Gsc(1 - c) + tc], 2 an 9п

(7)

для регенеративной части -

D(Mc) G д2c д rG (1 ) ¡ п

—------ =----2-~ rGsc (1 - c)-т c]. (S)

дt

2 9п 9п

Для решения уравнений (7) и (8) необходимо задать начальное условие и граничные условия на отборном и отвальном концах каскада. Принима-

ется, что распределение концентрации при (=0 соответствует некоторому стационарному состоянию каскада

с(п,0) = с0(п). (9)

Граничное условие на отборном конце каскада (п=Ы) при безотборном режиме работы определяется как

т( N, t ) = 0; т, =

тс + Gec(1 - c) -

dc_G_

dn 2

= 0.

n=N

Отсюда получим dc dn

= 2ec(1 - c).

(10)

n=N

На отвальном конце каскада (п=1) при безотвальном режиме работы граничное условие определяется как

т(1, t) = 0; t¡ =

' + G (i ) dcG '

-т c + Gs c(1 - c)-------------

dn 2

= 0.

n=1

Отсюда получим dc dn

= 2s c (1 - c).

(11)

n=1

В случае работы в безотборно-безотвальном режиме с резервуарами, имеющими газосодержание MT и МЛ на тяжелом и легком концах каскада соответственно, граничные условия имеют вид:

при

n = 1;

при n = N ;

* - 2sc = 2 Mpi*.

dn G dt

^ - 2sc = -2 MP- *.

dn G dt

(12)

На стыке ступеней разной производительности должны выполняться следующие условия непрерывности

т(+пк, 0 =т(-пк, 0, т, (+пк, 0 =т1(-пк, 0,

с(+пк, г) = с(-пк, 0- (13)

Знаками «+» и «-» отмечены значения параметров справа и слева от к-го стыка ступеней. После подстановки (5) и (6) в (12) получим С(-пк, г) дс~ С(+пк, г) дс+

2 дп 2 дп

+е[С(+пк, г) - С(-пк, г)] • с*(1 - с*) = 0. (14)

Здесь знаки «+» и «-» у производных имеют тот же смысл, а с=с(+пк,()=с(-пк,().

Если на «стыке» ступеней (п=пк) подается поток питания, то вместо уравнений (13) выполняются следующие уравнения

т(+пк, г)-т(-пк, г) = ^,

Т (+пк, г) ~т, (-пк, г) = с^,

с/ = с(+пк, г) = с(-пк, г).

После подстановки (10) и (11) в (14) получим

G(-nk, t) dc~ G(+nk, t) dc+

2 dn 2 dn

+s[G(+nk,t) - G(-nk,t)] ■ cf (1- cf) = cfF. (15)

2.1. Численный расчет нестационарных процессов

в разделительном каскаде

Расчет системы уравнений (3) для определения газосодержания и (7), (8) нестационарного переноса проводится следующим образом:

Рассчитывается начальное газосодержание в ступенях согласно (4) и начальная концентрация вдоль каскада по условию (9). Затем рассчитывается минимальный шаг по времени Дтдля уравнений (3) и (7), (8). Для уравнений (7), (8) вводится также шаг Дп по переменной п.

Для численного решения уравнения (3) используется модифицированный метод Эйлера [6] второго порядка точности по времени. При решении уравнений (7) и (8) используется двухслойная схема расчета, обеспечивающая второй порядок точности по переменным t и п [6].

Также как и в работах [7] для обеспечения расчета с общим вторым порядком точности при аппроксимации условий на отборном (10) и отвальном (11) концах каскада, на стыке ступеней разной производительности без подвода (13) и с подводом питания (15) использовались односторонние разностные схемы.

dc~

— = (-2 ■ d -з + 9 ■ cl -2 -18 ■ cl -! +11 ■ cl ) / 6Дп,

dn k k k k

dc+ ■ ■ ■ ■

— = (2 ■ c1 -3 - 9 ■ cn _2 +18 ■ < _! -11-cn )/6Дп,

СП k k k k

которые получаются из суперпозиции разложения функции с в ряды Тейлора в точках с индексами к-1, к-2, к—3, к+1, к+2, к+3 соответственно.

2.2. Результаты расчетов

Был проведен расчет нестационарного процесса разделения в симметричном противоточном каскаде, состоящем из 6 разделительных ступеней. В качестве начального условия задавалась начальная концентрация 0=0,00711. Считали, что все ступени каскада в начальный момент времени заполнены газовой смесью на заданную величину Мн. Задан внешний поток питания F=F(s3,t), как функция от времени и точки подачи потока питания в каскад. Результаты расчетов представлены на рис. 2-4.

На рис. 2 представлено распределение по вре-

мени т = ■

GS

Mn

(Gn, Mn- поток и газосодержание

последней ступени на стационарном режиме работы каскада) относительной концентрации на от-

С

борном конце каскада , где Св - требуемая кон-

Сп

центрация целевого продукта. На рис. 3 представлено распределение относительной концентра-

Ср

ции —— по каскаду при выходе на стационарный

Сп

режим работы. Относительное распределение газоМ,

содержания

Mn

по ступеням каскада при выходе

на стационарный режим работы показано на рис. 4.

s /\

Ступень

Рис. 2. Изменение относительной концентрации во времени на отборном конце разделительного каскада

Рис. 3. Распределение относительной концентрации по каскаду при выходе на стационарный режим

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

Заключение

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

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Филимонов С.В., Скорынин ГМ., Голдобин Д.Н. Моделирование нестационарных гидравлических процессов // Известия вузов. Физика. - 2004. - Т. 47. - № 11. - С. 173-177.

2. Обогащение урана / Под ред. С. Виллани. Пер. с англ. под ред. И.К. Кикоина. - М.: Энергоатомиздат, 1983. - 320 с.

3. Идельчик И.Е. Справочник по гидравлическим сопротивлениям. - М.: Госэнергоиздат, 1960. - 464 с.

4. Cohen K. The Theory of Isotope Separation as Applied to the Large Scale Production of U235. - N.Y.: McGraw-Hill, 1951.

5. Рабинович ГД., Гуревич Р.Я., Бобров ГИ. Термодиффузионное разделение жидких смесей. - Минск: Наука и техника, 1971. - 207 с.

6. Самарский А.А. Теория разностных схем. - М.: Наука, 1977. -656 с.

7. Левин Е.В., Лагунцов Н.И., Сулаберидзе Г.А. Об одном методе численного исследования нестационарных процессов в многосекционных установках для разделения многокомпонентных изотопных смесей // Инженерный физический журнал. -1982. - Т. 43. - № 3. - С. 456-462.

Поступила 20.12.2006 г.