CC BY
33
УДК 681.5
С.Л. Морева, О.И. Золотов, А.Л. Ляшенко
РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ ТЕМПЕРАТУРНЫХ ПОЛЕЙ АКТИВНОЙ ЗОНЫ РЕАКТОРА РБМК-1000
Высокий уровень надежности современных атомных станций - доказанный наукой и временем факт. Тем не менее, обеспечение безопасности является приоритетом ядерной промышленности, а ее совершенствованию уделяется самое пристальное внимание. Работа с высокими технологиями и огромная ответственность персонала атомных станций (в первую очередь перед населением) требуют, чтобы учитывались даже маловероятные риски возникновения аварии. Для этого проводится постоянный мониторинг состояния реакторной установки. Диагностирование физического состояния энергоблока позволяет выявить возможные нарушения функций безопасности, вызванные отказами оборудования и систем, оценить приоритетность восстановления этих функций и определить оптимальную последовательность действий оператора, направленных на обеспечение безопасности и перевод энергоблока в нормальный режим работы.
Описание объекта и постановка задачи
В качестве примера рассмотрим реактор АЭС РБМК-1000. Данный реактор является канальным реактором с графитовым замедлителем и водяным кипящим теплоносителем. В качестве топлива используется диоксид урана и02 с добавкой оксида эрбия Ег203. Активная зона представляет собой цилиндрическую кладку диаметром 11,8 и высотой 7 м, состоящую из вертикальных графитовых колонн [1].
Графитовая кладка цилиндрической формы состоит из вертикальных графитовых блоков сечением 250*250 мм. Блоки по вертикальной оси имеют сквозные отверстия диаметром 114 мм, предназначенные для размещения технологических каналов и каналов для контроля и управления. Часть технологического канала, размещенная в активной зоне, представляет собой циркониевую трубу диаметром 88 мм с толщиной стенки 4 мм. Для обеспечения теплового контакта с блоками кладки на трубу надеты графитовые кольца. Каждое чередующееся кольцо имеет не-
посредственный контакт по боковой поверхности либо с трубой, либо с внутренней поверхностью графитового блока, а также одно с другим по торцам. Это расположение определено из условия недопустимости заклинивания канала в кладке вследствие ее радиационно-термической усадки, а также термической и радиационной ползучести циркониевой трубы канала. В технологические каналы устанавливаются тепловыделяющие сборки (ТВС). Каждая ТВС содержит 18 твэлов (топливных стержней), дистанционированных специальными стальными решетками, и опирается на стальную концевую решетку. Твэлы диаметром 13,6 мм расположены по двум окружностям. На первом радиусе расположено 6 твэлов, на втором - 12.
Геометрическое расположение графитовой кладки, образующей вместе с топливными каналами активную зону, представлено на рис. 1.
Технологические каналы предназначены для размещения ТВС, дополнительных поглотителей и организации потока теплоносителя. Теплоноситель (вода) поступает в топливные каналы снизу, при температуре 270° С. Участок канала, на котором вода нагревается до температуры насыщения, имеет высоту около 2,5 м от низа активной зоны. На остальной ее части имеет место процесс развитого кипения. Расход теплоносителя по каждому топливному каналу может регулироваться независимо индивидуальным запорно-регулирующим клапаном.
Приблизительно 95 % энергии, выделяющейся в результате реакции деления ядер урана, передается теплоносителю, омывающему ТВС. Около 5 % мощности реактора выделяется в графите за счет замедления нейтронов и поглощения гамма-излучения. Эта доля тепла также передается теплоносителю через стенку технологического канала. Для снижения термического сопротивления, предотвращения окисления графита и улучшения теплопередачи от графита к технологическим и специальным каналам полость кладки заполнена циркулирующей газовой смесью гелия и азота, ко-
4
Рис. 1. Активная зона реактора
торая используется одновременно и для контроля целостности каналов путем измерения влажности и температуры газа [1].
Мониторинг параметров реакторной установки связан с необходимостью анализа непрерывных физических процессов, протекающих в активной зоне. С этой целью разработан ряд программных продуктов, рассмотренных в [2, 7]. Важная роль при исследовании ядерш>1х установок отводится расчету тепловыделяющей системы, математической моделью которой является уравнение теплопроводности, выведенное на основе первого закона термодинамики [3].
Составим математическую модель температурных полей для графитовой кладки.
Разработка математической модели
При моделировании рассматривались тепловые процессы, протекающие в графитовых кладках, набранных из графитовых блоков. На рис. 2 показан фрагмент кладки в количестве девяти блоков. В технологическом канале располагается тепловыделяющая сборка, омываемая теплоносителем. Высота графитовой кладки 8330 мм. Поверхности ТВС будем рассматривать как функции входа, а контрольные точки во внутренних углах графитовых кладок - как функции выхода. На рисунке представлены рассмаориваемые ячейки, ориентированные в горизонтальной плоскости осей пространственных координат X-Z, ось Y в данном случае является осью ординат.
Ниже представлен фрагмент математической модели, составленный для ячейки 5 (см. рис. 2).
Запишем дифференциальное уравнение теплопроводности применительно к графитовой кладке:
дХ.1 (х, н, г, t) dt
■ = вг
С аХм( х, н, г, t)
сх2
+ С Х5у(х, н, г, t) + С Х5у (х, н, г,t)
Сн
Сг2
(1)
BrWr
где аГ - коэффициент температуропроводности графита; х, у, х - пространственные координаты; !¥Г - мощность, выделяемая при торможении нейтронов; ХГ - коэффициент теплопроводности графита; УГ - объем графитовой кладки; Х2< х < ХЪ, У1< у <У2, 22 << <г3.
При анализе тепловых процессов в технологических каналах, уравнение теплопроводности удобнее записывать в сферических координатах.
Запишем дифференциальное уравнение теплопроводности применительно к циркониевой трубе:
СХ52 (^ Ф, R 0) dt
■ = а.
са н, ф, r, t)
Cr2
+ у сх .а (н, ф, R, t) + d 2а (н, Ф, R, t)
R
CR
Сн2
(2)
где ац - коэффициент температуропроводности циркония; .1 < у < У2, Л2 < Я<Я3 .
Запишем дифференциальное уравнение те-
Рис. 2. Схема расположения рассматриваемых ячеек технологических каналов
плопроводности применительно к теплоносителю в канале [5]:
дЫ у, ф, Я, х)
дХ
■ = а
д Т5з( у, ф, я, х)
дЯ2
1 дг5з(у,ф,я,х) , д2г,з(у,ф,я,х)
Я
дЯ
-»(Я, X) •
ду2
дг5з( у, ф, я, х) дф
- (3)
где В - скорость движения теплоносителя в канале; аТ - коэффициент температуропроводности теплоносителя; 71 < у < 72, Я1 < Я < Я2.
Функции входа и заданы в виде выражения:
и = Ттвс (у, ф, R, 0, (4)
где ТТВС - температура на поверхносаи ТВС; 71 < у < 72,Ф1 < ф < Ф2,.Я = Я1.
Яисленш>1е значеня температуры на поверхности ТВ С будем искать из значений энергораспределения, полученных с помощью информационно-вычислительной системы (ИВС) «Скала-микро».
Аналогично составлены дифференциальные уравнения теплопроводности для остальных ячеек.
В разработанной математической модели присутствуют уравнения, записанные в прямоугольной и цилиндрической системе координат. Подробно вопрос стыковки данных уравнений рассмотрен в [4]. Для их стыковки выделены зоны цилиндрической формы с радиусом R4 в графитовых блоках, они совпадают с границами сквозных отверстий в графитовых блоках (см. рис. 2). Для этих зон составлено дополнительно дифференциальное уравнение теплопроводности.
Запишем дифференциальное уравнение теплопроводности применительно к выбранной ячейке графитовой кладки:
дгп( у, ф, я, х) дХ
= аг
д %(у, ф, я, х)
дЯ2
+ дгп( у, Ф, Я, х) + д2Гп( у, Ф, Я, х) +
Я
дЯ
Я2
(5)
д2гп( у, ф, я, х)
ду2
агЩ,
71 < у < 72, Ф1 < ф < Ф2,Я3 < Я < Я4.
Данное дифференциальное уравнение отражает общий характер процесса. Для получения полной математической модели объекта, соответствующей конкретной задаче, необходимо за-
дание условий однозначности.
Граничные условия для графитовой кладки зададим соотношениями:
тп( у, ф, Я3, г) = т52( у, ф, Я3, г);
. эти(у,ф,яз,г) зт52(у,ф,яз,о
Аг - = АЦ-;
г дЯ „ дЯ
т51 (X, у, 2 2, г) = т21 (х, у, 2 2, г); . дТ51( х, у, 2 2, г) . дт21( х, у, 2 2, г)
ЛГ - = ЛГ - ;
г д2 г д2
т51( X 3, у, х, г) = Т61( X 3, у, х, г); . дт51(X3, у, х, г) А дТб1(X3, у, X, г)
Ар - = Ар - ;
г дХ г дX
т51 (х, у, 2 3, г) = тм( х, у, 2 3, г); . дт51( х, у, 2 3, г) А дт81( х, у, 2 3, г)
АГ-= Аг-;
г д2 г д2
тя( X 2, у, х, г) = т41( X 2, у, х, г); . дтм( X 2, у, х, г) . дТд{ X 2, у, х, г)
Аг-= Аг-
дX
дX
(6)
(7)
(8) (9)
(10) (11) (12)
(13)
(14)
(15)
где т21, тб1, т81, т41 - температура в графите ячеек 2, 6, 8 и 4 соответственно; ХГ - коэффициент теплопроводности графита; Хц - коэффициент теплопроводности циркония.
Граничные условия для циркониевой трубы зададим соотношениями:
т52( у, ф, Я3, г) = т51( у, ф, Я3, г); . дт52(у, ф, Я3, г) А дт51(у, ф, Я3, г)
А„ -= АГ-.
ц дЯ г дЯ
т52( у, ф, Я2, г) = тя( у, ф, Я 2, г); . дт52(у, ф, Я2, г) А дтя(у, ф, Я2, г)
Ац-^-= А Т-
(16)
(17)
(18) (19)
дЯ 1 дЯ
Граничные условия для теплоносителя зададим соотношениями:
тя(у, ф, Я 2, г) = т52(у, ф, Я 2, г); (20)
А,
дт^у, ф,Я2,0 , дт52(у,ф,Я2,г )
■=А „■
(21)
дЯ „ дЯ
Данная математическая модель получилась сложной. Решить полученную систему дифференциальных уравнений в частных производных аналитически (выделить передаточную функцию) не представляется возможным. Данные дифференциальные уравнения решались методом «сеток» [6].
Для численного анализа рассматриваемого объекта управления были составлены дискретная
модель уравнений и вычислительный алгоритм. Упрощенная блок-схема вычислительного алгоритма решения дифференциальных уравнений представлена на рис. 3.
В процессе составления дискретных моделей были решены задачи «стыковки» граничных условий, обеспечения устойчивости вычислительной схемы и выбраны шаги дискретизации по пространственным переменным.
Для проверки адекватности математической
Начало программы
Ввод данных
1/нгсльг по
, коорочнспналл
Обнуление массивов данных
ф = 1-\ Нет 1...7 >-—-
ичкпы ' бремени
1111КЛЫ но
Да
коороннсина.м
Уравнена я теплопроводности
Стыковка дифференциальных уравнений
1/НКЛЫ но
координации*
Граничные условия
Вывод данных
{ Конец программы у
Рис. 3. Блок-схема вычислительного алгоритма
Моделирование тепловым полей в активной зо>
О программе Справка Сохражггь Графики Печать
Ввод исходных данных
Температура теплоносителя на входе в канал
270
Лага 23 06 12 Текущее время 10 06 20 № контрольной точки
f^l
Скорость теплоносителя
38,4 мЗ/ч
Температура поверхности ТВС
950
Табличные константы AfP-GJ
B(P,G)
0,575
Расчет температуры в контрольной точке
141 Щ»
*4 ®
® О О
Расчет критической мощности ТК
Ввод новых данных
Время счета 5100
Отмена
Рис. 4. Окно ввода данных
модели и компьютерного моделирования специально разработано программное обеспечение. В качестве системы программирования выбран язык С++ и среда для разработки Borland C++ Builder. Выбор языка программирования С++ обусловлен тем, что он применяется на практике во многих отраслях промышленности, в т. ч. и на АЭС.
На рис. 4 представлен общий интерфейс программы.
С помощью разработанного ПО рассчитаны значения температур в контрольных точках (рис. 5), соответствующих местам расположения термопар, установленных в графитовой кладке на реакторе.
На рис. 5 графики А соответствуют данным, полученным с ИВС «Скала-микро», графики В соответствуют данным, полученным расчетным путем.
Результаты моделирования проверены с помощью информационно-вычислительной системы «Скала-микро», обеспечивающей контроль над работой реактора.
Значения, полученные путем моделирования, совпали с показаниями датчиков, установленных в графитовой кладке на реакторе. Полученные результаты свидетельствуют об адекватности разработанной математической модели.
Рис. 5. Графики переходных процессов:
а - график переходного процесса в контрольной точке 2; б - график переходного процесса в теплоносителе на выходе из канала
Результаты вычислительных экспериментов позволяют сделать вывод о том, что предложенное ПО точно моделирует тепловые процессы, протекающие в выделенном участке реактора, и может быть расширено с целью моделирования и мониторинга тепловых процессов во всем реакторе.
СПИСОК Л
1. Абрамов, М.А. Канальный ядерный энергетический реактор РБМК [Текст] / М.А. Абрамов, В.И. Авдеев, Е.О. Адамов; Под общ. ред. Ю.М. Черкашова. -М.: ГУП НИКИЭТ, 2006. -632с.
2. Киселев, А.Е. Моделирование аварий со снижением подачи теплоносителя для канала реактора РБМК-1000 [Текст] / А.Е. Киселев, А.В. Лемеха, В.Ф. Стрижов // Вопросы атомной науки и техники. Сер. Физика ядерных реакторов. -1998. -Вып. 1. -С. 48-65. -М.: Российский научный центр «Курчатовский институт».
3. Лыков, А.В. Теория теплопроводности [Текст] / А.В. Лыков. -М.: Высш. школа, 1967. -599 с.
4. Морева, С.Л. Математическое моделирование системы с распределенными параметрами на примере активной зоны реактора РМБК-1000 [Текст] / С.Л. Мо-
Разработанная методика моделирования позволит проводить диагностирование физического состояния энергоблока, выявлять возможные нарушения в его функционировании и прогнозировать возможные варианты работы при различных, в т. ч. и нештатных ситуациях.
ГЕРАТУРЫ
рева, А.Л. Ляшенко // Матер. IV Междунар. науч. конф. Системный синтез и прикладная синергетика. -Таганрог: Изд-во ТТИ ЮФУ, 2011. -C. 218-225.
5. Морева, С.Л. Разработка программного комплекса для моделирования тепловых процессов в активной зоне реакторных установок [Текст] / С.Л. Морева, А.Л. Ляшенко // Народное хозяйство. Вопросы инновационного развития. -М.: МИИ Наука, 2012. -С. 243-247.
6. Самарский, А.А. Численные методы математической физики [Текст] / А.А. Самарский, А.В. Гулин. -М.: Научный мир , 2000. -316 с.
7. Allison, C.M. SCDAP/RELAP5/MOD2 Code Manual [Text] / C.M. Allison, E.C. Johnson [et al.] // EG&G Idaho, INC. Sept. 1989.
УДК 004.8
В.В. Дармахеев
ЭФФЕКТИВНОСТЬ ПРИМЕНЕНИЯ РЯДА ГЕНЕТИЧЕСКИХ ОПЕРАТОРОВ СЕЛЕКЦИИ
Эволюционные вычисления для генерации моделей автоматов и преобразователей используются давно и успешно [6]. Однако при разработке генетических алгоритмов возникают проблемы. К одной из них относится трудность достижения сходимости алгоритма. Эта проблема тесно связана с выбором генетических операторов, включая операторы селекции. Настоящая статья посвящена вопросу выбора операторов селекции.
Генетический алгоритм разработан для решения прикладной задачи, заключающейся в генерации моделей преобразователей продукционных правил, представленных на естественном языке, в формулы логики предикатов первого порядка.
Данный генетический алгоритм (генератор) основан на следующих положениях [5]:
• автомат (особь) представляется в виде хромосомы;
• хромосома состоит из динамического набора молекул ДНК и имеет графовое представление, при этом одна молекула ДНК представляет собой описание одного состояния автомата;
• молекула ДНК состоит из фиксированного количества генов (минимальная неделимая составная часть молекулы ДНК);
• первая популяция Р(0) генерируется случайным образом, при этом значение каждого гена выбирается из соответствующей области допустимых значений.
Во время проведения вычислительных экспериментов была обнаружена проблема выбора особей, к которым впоследствии будут при-
