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

УДК 661.87.519

МОДЕЛИРОВАНИЕ НЕСТАЦИОНАРНЫХ ГИДРАВЛИЧЕСКИХ ПРОЦЕССОВ В КАСКАДАХ ЦЕНТРИФУГ ОБОГАЩЕНИЯ ИЗОТОПОВ УРАНА

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

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

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

Ключевые слова:

Разделение изотопов, газовая центрифуга, нестационарная гидравлика, математическое моделирование. Key words:

Isotope separation, gas centrifuge, nonstationary hydraulics, mathematical modeling.

Введение

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

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

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

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

В отличие от ранее известных работ [1, 2] рассматриваемая модель нестационарных процессов более точно описывает особенности построения каскада, в частности, учитывается наличие байпасных коммуникаций и наличие секционных компрессоров.

1. Описание математической модели

нестационарных процессов

1.1. Объемы, входящие в модель

В описываемой модели газосодержание /-ой ступени считается сосредоточенным в 18 объемах.

Восемь объемов роторов газовых центрифуг (ГЦ) секцийу=1,8 с газосодержанием М1уступени /. Газосодержание ^определяется давлением на носике отборника тяжелой фракции после прохождения скачка уплотнениярщу [3]:

M,.. = N V..■

\ j 1J \j

Hij '

(1)

где И} — количество ГЦ в у-ой секции, /-ой ступени, У1у, - коэффициенты зависящие от скорости вращения ротора со, температуры и геометрических размеров газовой центрифуги, Рщ - давление на носике отборника тяжелой фракции после прохождения прямого скачка.

Аналитическое выражение для газосодержания М1 у определено исходя из принятой упрощенной газодинамической модели ГЦ [4].

В газовые центрифуги секции у /-ой ступени подается поток питания g1 и выходят поток отбора g^ (легкой фракции) и поток отвала g1- (тяжелой фракции).

Траса отвала секции имеет газосодержание: М

М2 У = *Р*' (2)

где р- - давление в трассе отвала секции, У4у -объем трассы отвала. В объем входит поток отвала газовых центрифуг секции g^¡j и выходит поток отвала секции g у

Трасса отбора секции имеет газосодержание:

M + м,.. =— v, p.,

RT i

d. .rn

у у

где р+ - давление в трассе отбора секции, V/ -объем трассы отвала. В объем входит поток отвала газовых центрифуг секции g¡¡j и выходит выкачиваемый компрессором поток отбора секции g+j.

Трасса питания секции имеет газосодержание: М

М4. = ' (4)

где р - давление в трассе питания секции, ¥4/] -объем трассы питания. В объем входит поток питания ступени g2j и выходит поток питания газовых центрифуг секции gli¡.

Объем трассы отвала ступени с газосодержанием И5/ определяется давлением коллектора отвала Р- Этот объем включает в себя блочные коммуникации коллектора отвала до межступенного регулятора ВРД. Газосодержание трассы отвала ступени определяется выражением:

М _ М 5, = Р ,

(5)

I

р ., выходит суммарный поток отвала ступени

О-, определяемый расходом через межступенный регулятор ВРД.

Объем трассы отбора ступени с газосодержанием И6/ определяется давлением в коллекторе отбора ступени Р+:

ММ. = -

М КТ

К ■ р ,

61 I '

(6)

М, = -

М

ят

V ■ Р,

Ъ I '

(7)

ций

I

В объем байпасных коммуникаций линии отбора И8/ входят межблочные байпасные коммуникации ограниченные клапанами К1, К12, К5 1-го блока, и клапанами К3, К5 /-1-го блока.

Газосодержание данного участка И8/ определяется давлением в линии байпасных коммуникаций Р+■:

М„. = -

М КТ

КР-,

б, б, ~

(8)

где ¥5/ - константа зависящая от геометрических размеров трассы отвала.

В объем входит сумма потоков отвала ступеней

где ¥¡1 - константа, зависящая от геометрических размеров байпасных коммуникаций линии отбора.

В объём могут входит следующие потоки: поток отбора /-1 блока О¡м; закруточный поток Ош от клапана К12; поток приходящий с линии байпасных коммуникаций /-1-го блока О5/-1. Из объема выходит поток питания, проходящий через клапан К1.

В объем байпасных коммуникаций линии отбора И входят межблочные байпасные коммуникации ограниченный клапанами К1, К12, К5 /-го блока, и клапанами К3, К5 /-1-го блока.

Газосодержание данного участка И8 / определяется давлением в линии байпасных коммуникаций Р6+:

М„. = -

М КТ

V-Рб

б,1 б,1 '

(9)

где ¥6/ - константа зависящая от геометрических размеров трассы отбора.

В объем входит сумма потоков отвала ступеней

т

I g , выходит суммарный поток отвала ступени О+.

.=1

Объем трассы питания ступени с газосодержанием И7/ определяется давлением в коллекторе питания ступени р

где V ¡6 / - константа, зависящая от геометрических размеров байпасных коммуникаций линии отбора.

В объем могут входит следующие потоки: поток отбора /-1 блока О¡/-1; закруточный поток Ош от клапана К12; поток, приходящий с линии байпас-ных коммуникаций /-1-го блока О5/-1. Из объема выходит поток питания проходящий через клапан К1.

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

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

СМ, ■. —

_ gui _ , . = 1m, (10)

л

где ¥7/ - константа, зависящая от геометрических размеров трассы питания.

В объем могут входить следующие потоки: поток отбора предыдущей ступени О¡/-1; поток отвала последующей ступени Овнешний поток питания Т0; внешний поток воздуха Т„; поток закрутки Озак{.

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

Ж

СМ

gUj к\ g_ ,

С

СМ..

= gi _ к3 g+,

с

■ = к 2 g 2 _

2У

СМ 5, а

= I к,■ р_■. _ К2.в-,

/ > 21. 2,1 I 5

(11)

(12)

(13)

(14)

.=1

.=1

сМ6.

_6.

Ж СМ 7.

_7.

л

= У к 2,

„ 2 - К3 в+

и 21. 3,1 I

= К40+ К 0+ +

4, г+1 1, г+1

(15)

+<51Т0 + 82Т -Ук2 . .,

0 и / , . о 2. '

и=1

СМ„

с

■ = К2 в + к6 в + 2, 6,

+к„.в -к41в.. -к6.-о°-,

12,1 I 4,г-1 г-1 6,г4 г4 5

Ж

= К3 в' + к5.о -

3,1 5, I

-кш-1°-1 - к1„+1°'+1 - кы +1°5+1'

Работу секционного компрессора можно описать следующим выражением:

Ап2 + Б^ 2+,. + СЯ2+,и2 = р+ (е-1), (22)

где А, В, С - константы, получаемые экспериментально и учитывающие геометрические размеры компрессора, его гидравлическое сопротивление, п

(16) - частота вращения компрессора, е - показатель степени сжатия компрессора. Решение этого квадратного уравнения дает величину потока отбора g у-секции с учетом работы секционного компрессора.

Для определения текущего значения расхода га-

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

в- = А • (Нф)Р-, (23)

(18)

где к1$, к2р к30- состояние секционных клапанов на трассах отвала, питания, отбора (0 - клапан закрыт, 1 - клапан открыт); Ки-К6;, Кш состояние блочных клапанов на трассах отвала, питания, отбора (0 - клапан закрыт, 1 - клапан открыт).

Ур. (10) составляется для каждой секции ступени /. Ур. (11-18) составляются для каждой ступени.

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

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

а) турбулентного

в =

Р2 - Р2

С

(19)

турб

где Стурб - постоянный коэффициент; б) ламинарного

Р - Р2

в =

С

(20)

= N..

Ь1.. + Ь2.. + Ь3..

N.. N.. ,

+Ь4 р + Ь5 р 2 + Ь6 р

. . . . . .

(21)

пС

где А =-и , ик, - критическая скорость, а -

ЯТ ^ р

диаметр клапана, Н - положение клапана регулятора, зависящее от конструкции регулятора, ф - коэффициент расхода, зависящий от положения клапана регулятора и описываемый с помощью дифференциального уравнения:

С П (АР. ~+АР) +ц

— ге! '-, (24)

Ж Т

где /л = -

АН Н

- относительное изменение положе-

ния заслонки регулятора; Трег - постоянная времени регулятора; Прег - коэффициент усиления регулятора.

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

цЦ) = я (АР- +АР)(1 -е"

(25)

где Сшм - постоянный коэффициент.

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

С (~ \ (~ \2 ^

+

где Ь19...Ь6^ - коэффициенты гидравлической характеристики для определенного типа газовых центрифуг в рассчитываемой ступени. Эти коэффициенты зависят от частоты вращения ротора и определяются экспериментально.

2. Схема решения уравнений нестационарной гидравлики

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

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

и=1

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

Таким образом, задача моделирования нестационарных процессов в промышленных центрифужных каскадах сводится к решению системы дифференциальных уравнений первого порядка с заданными начальными условиями. Задача Коши за редким исключением не имеет аналитических методов решения. В численном методе вместо отыскания непрерывных функций времени рассчитываются значения искомых функций в дискретные моменты времени, а система дифференциальных уравнений заменяется, тем или иным способом, разностными уравнениями. В применяемом алгоритме использована неявная схема Эйлера с пересчетом [5], которая обладает третьим порядком точности на шаге и вторым на интервале.

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

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

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

Результаты расчетов нестационарных режимов работы ТС, для которых имеются экспериментальные данные (изменение величин внешних потоков, закрытие блоков [1]), позволяют судить об адекватности описанной модели поведению реального объекта. Расчетное время распространения

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

1. Филимонов С.В., Скорынин Г.М., Орлов А.А., Голдобин Д.Н. Моделирование нестационарных гидравлических процессов в промышленных центрифужных каскадах обогащения урана // Известия Томского политехнического университета - 2006. -Т. 309. - № 3. - С. 85-88.

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

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

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

Заключение

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

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

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

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

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

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

Работа выполнена при финансовой поддержке гранта 09-08-9915-рофи.

3. Дейч М.Е. Гидрогазодинамика. - М.: Энергия, 1984. - 384 с.

4. Идельчик И.Е. Справочник по гидравлическим сопротивлениям / Под ред. М. О. Штейнберга. - 3-е изд., перераб. и доп.-

М.: Машиностроение, 1992. - 672 с.

5. Пирумов У.Г. Численные методы. - М.: Дрофа, 2004. - 224 с.

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