Содержание к диссертации
ВВЕДЕНИЕ 3
На защиту выносятся результаты. 14
ГЛАВА I. Задачи математического моделирования ядерных реакторов и их защиты. 15
Математическое моделирование ядерных реакторов. 15
Постановка задач расчета ядерных реакторов и их защиты
в многогрупповом приближении. 18
1.3. Приближение индикатрисы рассеяния полиномами Лежандра
( Рт - приближение). 20
1.4. Угловые квадратуры Sn и ESn метода дискретных ординат. 22
ГЛАВА II. Итерационные методы решения. 26
2.1. Ускорение сходимости внешних итераций. 26
2.1.1. 8 - процесс ускорения сходимости итераций в задачах расчета
значений Keff и источника деления. 26
Исследование скорости сходимости 8 - процесса в задаче поиска наибольшего собственного значения неотрицательной неразложимой квадратной матрицы А . 32
Результаты тестовых расчетов критической сборки GODIVA
в одномерной сферической геометрии. 38
2.1.4. Результаты расчета значений Kcjf и источника деления исходного
состояния критической сборки BZD/1 в экспериментах «ZEBRA»
в Х-Y-Z геометрии. 41
2.2.Ускорение сходимости внутренних итераций. 47
Дискретизация уравнения переноса и метод решения. 47
Метод пространственного ребаланса. 49
Фурье анализ устойчивости метода пространственного ребаланса совместно с «алмазной» схемой на примере одногрупповой плоской
задачи для изотропного рассеяния в бесконечной среде. 54
2.2.4. Результаты тестовых расчетов железоводной борированной защиты
в одномерной геометрии. 61
ГЛАВА III. Реализация метода дискретных ординат и эффективных методов ускорения
сходимости в трехмерных модулях KIN3D, KINRTZ и KIN3D6 пакета
«РЕАКТОР» для моделирования ядерных реакторов и их защиты. 63
Описание модуля KIN3D. 64
Описание модуля KINRTZ. 65
Описание модуля KIN3D6. 65
Входная информация к модулям. 65
Выходная информация модулей. 70
Информация для продолжения работы модулей. 71
Общая блок-схема управляющей программы KIN модулей KIN3D,
KINRTZ и KIN3D6 и основной подпрограммы DSNPN. 73
ГЛАВА IV. Проверка эффективности предлагаемых алгоритмов
на математических моделях реальной ядерно-энергетической установки. 76
Расчет значений Keff и источника деления реактора СВБР 75/100. 76
Расчет радиационных полей в защите реактора СВБР 75/100. 79
ЗАКЛЮЧЕНИЕ 87
ЛИТЕРАТУРА 88
Приложение I. Метод верхней релаксации для нахождения поправок метода ребаланса. 95
Приложение И. DD0 схема в HEX-Z геометрии. 99
Приложение III. 0-WDD схема в X-Y-Z и R-
Введение к работе
Развитие безопасной ядерной энергетики является одной из актуальных задач современной технологии. Использование ядерной энергетики обеспечивает энергетическую независимость страны и дает явные экономические преимущества. Одним из новых направлений развития ядерной энергетики является создание ядерных реакторов малой и средней мощности на быстрых нейтронах с тяжелым жидкометаллическим теплоносителем, которые необходимы для развития регионов крайнего севера и дальнего востока. Продолжается дальнейшее развитие быстрых реакторов различного типа и назначения. При проектировании таких реакторов необходимо обеспечить выполнение всех норм ядерной и радиационной безопасности, правильно рассчитать дозы облучения и оценить надежность конструкционных материалов. Проектирование невозможно без эффективного решения задач математического моделирования ядерных реакторов, максимально приближенных к реальности.
Одной из таких задач является проблема численного решения уравнения переноса нейтронов и гамма-квантов с детальным описанием геометрии ядерного реактора и с подробной зависимостью от энергетической переменной. Возникающие при этом системы конечно-разностных уравнений обладают высокой размерностью. Итерационные методы решения таких систем, как правило, очень медленно сходятся. Проблема разработки и использования эффективных методов решения является весьма актуальной задачей, которой посвящена диссертация.
В течение более 20-ти лет в Институте Прикладной Математики им. М.В. Келдыша РАН под руководством А.В. Воронкова проводились работы по созданию численных методов и программ решения стационарных и нестационарных систем многогрупповых уравнений переноса частиц в приближении метода дискретных ординат [1] с учетом структуры расчетной области в различных геометриях [2] - [37]. Использование таких программ позволило в 1994 году впервые в нашей стране провести трехмерные расчеты уравнения переноса [8] для получения натриевого пустотного коэффициента реактивности критической сборки BZD/1 [9] - [10] в экспериментах «ZEBRA» [38] в SA приближении метода дискретных ординат, в Р{ приближении индикатрисы рассеяния и в 26 групповом разбиении по энергии. Численное решение получено методом итерации источника [1] без использования алгоритмов ускорения сходимости итераций в рамках пакета прикладных программ «РЕАКТОР» [39]. Пакет «РЕАКТОР» разработан по модульному принципу и использовался ранее для математического моделирования ядерных реакторов на быстрых нейтронах в диффузионном приближении. Поэтому все программы решения
многогрупповых систем уравнений переноса оформлены в виде модулей. В состав пакета «РЕАКТОР» входят модули различного типа и назначения. Обмен информацией между модулями осуществляется через общий Архив. Системное обеспечение пакета служит для создания Архива данных, ввода информации и сопровождения работы каждого модуля. Имеются геометрические модули, модуль ввода ядерно-физических данных, модули расчета констант, модуль задания параметров рассчитываемой задачи, модуль обработки результатов расчета и др. Большинство модулей используются до расчета решения систем многогрупповых уравнений переноса частиц, подготавливая и записывая в Архив исходные данные.
Для решения уравнения переноса нейтронов и гамма-квантов в одномерной геометрии в пакете «РЕАКТОР» использовалась положительная AWDD схема метода дискретных ординат совместно с согласованной PXSA [11] схемой ускорения сходимости внутренних итераций. Пакет программ KINRZ дополнен двумя методами коррекции «алмазной» схемы [1] для получения положительного решения: схемой с «нулевой» коррекцией потока (DD0 схема) и методом коррекции по шаговой схеме (DDST схема) [1], программами для расчета индикатрисы рассеяния в Рт приближении и DSA — методом
ускорения внутренних итераций [16]. Разработан алгоритм [17] и создана программа решения многогруппового стационарного уравнения переноса нейтронов и гамма-квантов на сетках, согласованных со структурой расчетной области в двумерной R-Z геометрии. Разработана нестационарная методика решения уравнения переноса в R-Z геометрии на сетках, состоящих из произвольных выпуклых четырехугольников [18], [20] и получены результаты тестовых расчетов. Для проведения расчетов активной зоны ядерного реактора на быстрых нейтронах с гексагональной решеткой ТВЭЛов разработан метод и создана программа решения уравнения переноса в HEX-Z геометрии [21] с использованием DD0 схемы на сетках, максимально согласованных с конфигурацией расчетной области. Первоначально написанные в однопроцессорном виде трехмерные программы решения уравнения переноса перенесены на параллельные ЭВМ кластерного типа [22] - [26] в пакете «REACTOR-P».
Настоящая работа посвящена разработке и исследованию новых эффективных методов ускорения сходимости итераций при решении однородных и неоднородных систем линейных уравнений, возникающих в задачах расчета ядерных реакторов и их защиты. При решении однородных задач расчета Kejf - эффективного коэффициента
размножения и источника деления ядерного реактора в многогрупповом приближении различают внешние и внутренние итерации. Внешними итерациями называются итерации
по источнику деления, а внутренними - итерации по внутригрупповому рассеянию с фиксированным источником, состоящим из источника за счет переходов из других групп и источника деления. При решении неоднородных задач защиты ядерных реакторов внешние итерации отсутствуют, если не учитывается рассеяние в верхние энергетические группы.
Наиболее полный обзор современных методов ускорения сходимости итераций для решения уравнений переноса частиц методом дискретных ординат дан в работе [40]. Плохую сходимость внешних итераций обычно связывают с близостью двух наибольших
по модулю собственных значений матрицы Т системы уравнений TQ = АО, когда доминантное отношение модуля второго собственного значения к наибольшему первому близко к единице. Одним из известных методов ускорения сходимости внешних итераций является метод полиномов Чебышева [41], [42]. Этот метод организует внешние итерации
вида Q(p+]) =Cn(T)Q{p) с помощью полиномов Чебышева Сп(Т) степени п, построенных
на основе информации о двух наибольших по модулю собственных значениях матрицы Т. Для его использования надо предварительно оценить величину доминантного отношения. Получение оценки сравнимо по сложности с решением исходной задачи. Кроме того, метод полиномов Чебышева слишком чувствителен к вычислительной погрешности [1], [43].
В методе сдвига [40] для ускорения сходимости внешних итераций предлагается подобрать положительную константу к > KeJf для того, чтобы доминантное отношение у
новой задачи стало меньше, чем у исходной. Новая задача получается путем добавления к обеим частям системы уравнений отрицательного слагаемого, записанного на разных итерациях в виде источника деления, деленного на к. В этом случае приходится переопределять источник деления на каждой внутренней итерации, что приводит к усложнению процесса сходимости внутренних итераций в каждой группе. Кроме того, в задачах с доминантным отношением близким к единице трудно заранее угадать наиболее подходящую константу к.
В методе обратных итераций Виланда [41] предлагается выбрать константу «/о> близкую к наибольшему собственному значению и затем организовать итерационный процесс, как в методе сдвига, чтобы получить соответствующий собственный вектор. Полученную задачу нельзя свести к последовательному решению одногрупповых уравнений [41].
В методе ребаланса на грубой сетке [43] предлагается ускорять внешние и внутренние итерации одновременно. Схема решения - двухступенчатая: расчеты чередуются на основной и на грубой сетке.
В итерационном методе подпространств Крылова [44] ищутся несколько собственных значений, близких по модулю к наибольшему собственному значению, и им соответствующие собственные векторы. Этот метод тем эффективнее, чем больше собственных значений и соответствующих им собственных векторов будут найдены, что, как правило, усложняет исходную задачу, т.к. вместо одного собственного значения ищутся несколько, и поэтому увеличивается трудоемкость решения реальных задач.
В методе квазидиффузии [45], [46] в основном решается многогрупповая система квазидиффузионных уравнений, а решение системы уравнений переноса используется для поправки диффузионного дробно-линейного функционала. Для расчета Кс1Г используется
система кназидиффузионных уравнений, полученная суммированием исходной системы по энергетическим группам. Расчет Keff происходит быстрее за счет того, что более
трудоемкое решение системы уравнений переноса осуществляется редко. Метод квазидиффузии хорошо зарекомендовал себя в одномерной геометрии и в двумерной R-Z геометрии. Неизвестно, насколько этот метод эффективен в трехмерной геометрии.
На этом список известных методов ускорения сходимости внешних итераций практически заканчивается. Среди других методов ускорения выделяются методы, основанные на идее Л.А. Люстерника [47], предложенной для решения систем неоднородных линейных уравнений
х = Вх + с (в.1)
методом простой итерации:
х(р) =вх{р-х) +д.
Установлено [47], что при большом числе итераций р векторы х{р) -х(р+1) и х(р+,) - х(р+2) практически пропорциональны, а отношение компонент этих векторов приблизительно равно наибольшему собственному значению \ матрицы В. Если известны
последовательные приближения х(/7),Зс(р+1),х(р+2), то можно оценить Я, по какой-либо i-той компоненте:
х(Р+1) _ х(Р+2)
а затем найти х, используя приближенное равенство:
x*xlp)+-l—(xlp+l)-xip)). (в.З)
1-я,
Метод Люстерника (в.2) - (в.З) был модифицирован и использован для ускорения сходимости итераций при решении неоднородной одногрупповой задачи (в.1) переноса нейтронов [48]. Величина параметра Я(р) < 1 вычислялась по формуле:
^p)-(p(p-X)\dV
Л(р)=^г
jy"4) V'~2)
где ф{р) - приближение скалярного потока на итерации р , G - пространственная область, в которой рассматривается процесс переноса нейтронов. Линейная экстраполяция выполнялась через каждые четыре итерации по формуле:
<Р = <Р(Р)+~^(<Р{Р)-<Р(Р~Х)) (в.4)
с последующей нормировкой полученного результата для того, чтобы сохранить полный баланс частиц. Сделан вывод о том, что метод Люстерника может быть более эффективным, если его использовать в ходе итерационного процесса [48].
За рубежом метод Люстерника известен, как метод асимптотической экстраполяции источника ASE (Asymptotic Source Extrapolation) [40]. Этот метод применяется после выполнения нескольких простых итераций по источнику рассеяния. Оценка собственного значения Я с наибольшей амплитудой гармоники ошибки делается по формуле отношения двух норм:
:(/>+) _ ;:(/>)
Л*
Г, (в.5)
Wx^-x^l
а линейная экстраполяция по формуле:
x = xip+l,2) +-^(xlp+l,1) -xip)). 1-Я
Сделан вывод о том, что метод ASE не препятствует неприемлемо медленной сходимости
внутренних итераций [40], если имеется много гармоник.
При использовании метода Люстерника для ускорения сходимости итераций
возникает проблема разработки критериев его эффективного применения. В [49] описан
близкий к методу Люстерника алгоритм ускорения сходимости итераций при решении
неоднородной системы линейных уравнений (в.1), который называется «д2 - процесс».
Этот метод рассмотрен для случая, когда матрица В имеет простую структуру. В 82 -
процессе компоненты приближения 5с(р) уточняются в ходе итерационного процесса после выполнения условия пропорциональности двух векторов [49]:
juip) =
х^-х1р-2)\\\\х1р)-х1р-1)
1, (в.6)
где в числителе стоит скалярное произведение векторов, а в знаменателе их евклидовы нормы. Если условие (в.6) выполнено, то вычисляется собственное значение Л\р) матрицы В задачи (в.1) по формуле:
(x{p)-xll-l),xlp)-xlp-») (jc^-jc^-2',^-^"-0)
^ ^—f/3-11 -Гп-21 -(п) -C/j-lK - \В-')
Это собственное значение используется для уточнения приближения х к решению х :
у(Р) =х-м +^L—(x(p) -х(р~])). (в.8)
1 - Л\р)
В [49] получены следующие оценки /^^=^+0(^/^), х = х^ +0(\At\P) и х = у(р) + 0(\Л2\Р), где х - точное решение (в.1). Сделан вывод о том, что у(р) является
лучшим начальным условием для последующих итераций по сравнению с х(р). Производя время от времени уточнения (в.7) и (в.8), иногда удается существенно уменьшить общее число итераций. Для того чтобы делать такие уточнения, необходимо чтобы в разложении в ряд по собственным векторам п -мерной матрицы В невязки начального приближения:
х(р)-х = ^Л?а,ё, i=i
одно из слагаемых преобладало над другими [49]. Условие (в.6) называется условием
практической применимости линейной экстраполяции (в.7), (в.8). Кроме того, в [49]
предложена возможная схема метода простой итерации с применением S2 - процесса
ускорения сходимости. Отмечено наличие ряда моментов, затрудняющих реализацию
метода и требующих серьезной математической подготовки и проведения большой серии
численных экспериментов [49].
b(/>+i) -х{р+2)
В [50] предложено использовать метод Люстерника (в.З), где Л1 = -у-
xlp)-xip+l)\
для ускорения сходимости итераций при решении задачи Ах = х, чередуя выполнение нескольких обычных итераций х(р) = Ах(р'Х) с расчетом уточненного приближения, которое каждый раз принимается за новое начальное приближение.
В данной работе предложен разработанный автором новый метод ускорения сходимости итераций, который называется 5 - процесс, для решения задачи нахождения наибольшего по модулю собственного значения матрицы системы линейных уравнений. В задачах расчета Кс// этот метод служит для ускорения сходимости внешних итераций в
случае их замедления и заключается в использовании линейной экстраполяции полученного приближения источника деления и соответствующих ему скалярных потоков вблизи положительных корней характеристического многочлена итерируемой матрицы, а также пересчета моментов угловых потоков. Этот метод является дальнейшим развитием метода Люстерника [47] и S2 - процесса Бахвалова [49] для решения однородных задач.
Проблема ускорения сходимости внутренних итераций возникает при решении уравнений переноса в одной энергетической группе методом итераций по столкновениям. Этот метод сходится быстро в задачах с оптически тонкой средой или с высокими скоростями реакции захвата. Однако он сходится медленно в задачах с оптически толстыми средами и в тех случаях, когда отношения c-~Z0S /Ег для разных веществ
близки к единице, где 0 s. - нулевой момент сечения рассеяния и Еу. - полное сечение.
Одним из способов ускорения сходимости внутренних итераций является метод
ребаланса [51], [52], дополнительно требующий на каждой итерации относительно
небольшого количества арифметических действий. В тех случаях, когда этот метод
устойчив, он значительно ускоряет расчеты. Это нелинейный метод, основанный на
дополнительном расчете на каждой итерации по столкновениям решения системы
диффузионно-подобных разностных уравнений для мультипликативных
пространственных поправок к нулевому и первому угловым моментам решения в Р] приближении индикатрисы рассеяния. Этот метод привлекателен еще и тем, что если решается задача в Рт приближении индикатрисы рассеяния, то для сохранения баланса
частиц при ускорении все угловые моменты решения умножаются на соответствующие поправки в каждой пространственной точке. Такой возможности лишены DSA метод [53] и PXSA метод [54] ускорения сходимости внутренних итераций, т.к. в них используются аддитивные пространственные поправки. В DSA методе поправки ищутся к нулевому моменту решения, а в PXSA методе - к нулевому и первому угловым моментам решения.
Различают два способа ускорения методом ребаланса: на мелкой сетке (fine-mesh rebalance) и на грубой сетке (coarse-mesh rebalance). В методе ребаланса на мелкой сетке [55] уравнение для поправок записывается на сетке, заданной для решения уравнения переноса. В методе ребаланса на грубой сетке для расчета поправок используется более
крупная сетка, объединяющая две или несколько ячеек мелкой сетки. Автором проведено теоретическое исследование устойчивости метода ребаланса в обоих случаях [56], [57] с помощью Фурье-анализа линеаризованного метода пространственного ребаланса совместно с «алмазной» схемой Sn метода дискретных ординат в одномерной плоской
геометрии на равномерной сетке для случая изотропного рассеяния частиц на специальном классе задач: с постоянными сечениями и постоянным изотропным источником в бесконечной области. Результаты анализа показали, что для любого фиксированного 0 < с < 1 область устойчивости метода зависит от величины параметра Y.Th, где h - шаг сетки. Для достаточно малых значений с метод ребаланса устойчив и эффективно ускоряет в обоих случаях. При приближении значений с к 1 оба метода могут быть неустойчивыми. Скорость сходимости метода ребаланса становится медленной и метод даже расходится на сетках, образованных очень большим числом мелких ячеек [57]. Метод ребаланса на мелкой сетке устойчив в очень узкой области изменения параметра Е7./г ж 1 для значений с « 1 [56], [40].
В [58] предложен способ повышения устойчивости метода ребаланса на мелкой сетке в одномерной геометрии путем введения фиктивных граничных токов на сторонах ячейки. Подтверждением устойчивости этого метода служат результаты численных расчетов [58]. Этот усовершенствованный метод называется методом пространственного ребаланса или методом частичных токов. Этот алгоритм обобщен на случай двумерной геометрии и предложена другая формула для фиктивных токов [59]. На основе результатов расчетов сделан вывод о том, что новая формула предпочтительнее для значений « с », близких к единице. Дальнейшее усовершенствование методики ускорения состоит из двух модификаций [60]: получены граничные условия для нахождения мультипликативных поправок, согласованные с граничными условиями уравнений переноса и способствующие ускорению и разработан алгоритм увеличения фиктивных токов с целью достижения эффективности ускорения в наиболее трудных задачах. Устойчивость окончательного варианта метода пространственного ребаланса исследована автором [61] с помощью Фурье-анализа линеаризованного метода в одномерной плоской геометрии.
Усовершенствованы имеющиеся программы в X-Y-Z геометрии (модуль KIN3D) и в HEX-Z геометрии (модуль KIN3D6) [26], [28] - [31] и разработаны новые программы (модуль KINRTZ) в трехмерной R-(p-Z геометрии [27] для проведения трехмерных расчетов эффективного коэффициента размножения Kcff ядерных реакторов на быстрых
нейтронах и источника деления в однородных задачах, а также потоков нейтронов и
гамма-квантов в задачах их защиты. В модулях KIN3D и KINRTZ реализована взвешенная 0 - схема [62] (0-WDD схема) метода дискретных ординат для расчета задач защиты. Во все трехмерные модули добавлены программы расчета индикатрисы рассеяния в Рт
приближении. Для ускорения сходимости внутренних итераций реализован метод пространственного ребаланса, а для ускорения сходимости внешних итераций разработан и реализован новый метод 8 - процесс [36]. Результаты трехмерных расчетов реакторной установки СВБР 75/100 [63] (Свинцово-Висмутовый Быстрый Реактор) эквивалентной энергетической мощностью 75/100 МВт и его радиационной защиты, полученные по этим программам, доложены на конференциях [32] - [35], [37]. Данная работа построена следующим образом.
В первой главе описаны две основные задачи математического моделирования ядерных реакторов и их защиты. Раздел 1 посвящен проблеме математического моделирования ядерных реакторов на основе использования линеаризованного уравнения Больцмана для описания переноса нейтронов и гамма-квантов. В разделе 2 приведена постановка этих задач в многогрупповом приближении. В разделе 3 описано приближение индикатрисы рассеяния в виде ряда по полиномам Лежандра Рт до степени т, т.е. в Рт -
приближении. В разделе 4 описаны угловые квадратуры метода дискретных ординат, используемые в программах.
Во второй главе рассмотрены итерационные методы решения основных задач. Первый раздел посвящен новому методу 8- процессу ускорения сходимости внешних итераций. В параграфе 1 описан степенной метод решения задачи расчета значений Kcff и
источника деления, приведено свойство включения, позволяющее оценить на каждой внешней итерации верхнюю и нижнюю границы величины Keff, перечислены критерии
окончания итерационного процесса. Обоснована возможность рассмотрения эквивалентной задачи, имеющей наибольшее собственное значение, равное единице, путем разложения начального приближения в ряд по собственным и присоединенным векторам итерируемой матрицы. Затем описан предлагаемый алгоритм ускорения сходимости внешних итераций и приведены его формулы. В параграфе 2 проведено исследование скорости сходимости 8 - процесса на примере решения вспомогательной задачи Ах = х поиска наибольшего собственного значения неотрицательной неразложимой квадратной матрицы А. Показано, что формулу линейной экстраполяции 8 - процесса можно интерпретировать, как шаг метода Ньютона [49], сделанный вблизи положительного корня характеристического уравнения итерируемой матрицы А. Представлена блок-схема 8 - процесса. В параграфе 3 приведены результаты тестовых
расчетов Kejr критической сборки GODIVA [64] в одномерной сферической геометрии в SgP, приближении с различными значениями параметра д и без ускорения. Получена оценка эффективности ускорения. Показаны графики поведения значений параметров К(р) и Я(р) на внешних итерациях без ускорения и с использованием д - процесса ускорения для 8 = 0.9. В параграфе 4 рассмотрены результаты расчетов Kcjf и источника
деления исходного состояния критической сборки BZD/1 в экспериментах «ZEBRA» [38] в X-Y-Z геометрии в 30-ти групповом приближении по энергии для S4Pt и SSP3
приближений методом дискретных ординат. Сделан вывод о возможности выбора эффективного значения параметра S на основе анализа результатов предварительных расчетов в более низком S4P{ приближении.
Второй раздел посвящен методу пространственного ребаланса. В параграфе 1 описан алгоритм решения уравнения переноса в одной группе на одной внутренней итерации с использованием «взвешенной» схемы [1] в X-Y-Z геометрии. В параграфе 2 приведены формулы метода пространственного ребаланса и описана эффективная методика его использования. В параграфе 3 проводится исследование устойчивости метода пространственного ребаланса с «алмазной» схемой [61] с помощью Фурье-анализа линеаризованного метода в одномерной плоской геометрии на равномерной сетке. Рассмотрен случай изотропного рассеяния частиц на специальном классе задач: с постоянными сечениями и постоянным изотропным источником в бесконечной области. Устойчивость этого метода ранее не исследовалась. Результаты расчетов показали, что Фурье-анализ корректно предсказывает свойства устойчивости и сходимости нелинейного метода пространственного ребаланса [61]. В параграфе 4 приведены результаты тестовых расчетов железоводной борированной защиты [11] в одномерной геометрии 0-WDD схемой совместно с методом пространственного ребаланса. Для сравнения эффективности ускорения приведены результаты расчетов этой задачи AWDD схемой совместно с P^SA методом, выполненные ранее в пакете «РЕАКТОР».
Третья глава посвящена реализации метода дискретных ординат и эффективных методик ускорения сходимости внутренних и внешних итераций [35] - [37], [61] в трехмерных программных модулях пакета «РЕАКТОР», предназначенных для расчета значений Kcff и источника деления и задач защиты. Приведено описание каждого модуля.
В первом разделе описан модуль KIN3D расчета потока частиц в X-Y-Z геометрии в выпуклых областях, горизонтальное сечение которых задается на равномерной квадратной сетке. Во втором разделе описан модуль KINRTZ, предназначенный для
расчета потока частиц в R-(p-Z геометрии в цилиндрических областях с неравномерной сеткой. В третьем разделе описан модуль KIN3D6 расчета уравнения переноса в HEX-Z геометрии в выпуклых областях, горизонтальное сечение которых задается на сетке, состоящей из правильных шестиугольников. В четвертом разделе приведена общая блок-схема модулей и их основной подпрограммы DSNPN.
В четвертої! главе приведены результаты решения задач математического моделирования проектируемого в настоящее время энергетического ядерного реактора СВБР 75/100 [63] и его защиты в трехмерной геометрии. Продемонстрирована эффективность предложенных алгоритмов и созданных программ. В первом разделе приведены результаты расчетов Keff и источника деления реактора СВБР 75/100 в HEX-Z
геометрии в 30-ти групповом iSgP, приближении для различных значений параметра д от
0.2 до 0.9 по DD0 схеме. Эта же задача рассчитана в более высоком S&P^ приближении
индикатрисы рассеяния с параметром S = 0.7, при котором результаты расчетов в SgP^
приближении были в 2 раза эффективнее по времени. По результатам расчетов разработаны рекомендации о выборе значения параметра д. Во втором разделе рассмотрены результаты расчета радиационных полей в защите реакторной установки СВБР 75/100 (30 нейтронных групп и 19 гамма групп) в X-Y-Z и R-q>-Z геометрии с использованием 0-WDD схемы, в S& приближении метода дискретных ординат и в Р3 -
приближении индикатрисы рассеяния. Результаты расчетов сопоставлены с аналогичными результатами, полученными по программе TORT [65] при использовании одной и той же системы констант. Сделан вывод о том, что на эффективность ускорения внутренних итераций методом пространственного ребаланса влияют в основном два фактора: выбор номера итерации, с которой начинается ускорение итераций, и выбор 0 параметра схемы 0-WDD. Приведены рекомендации по оптимальному выбору этих параметров.
В Приложении I описан метод верхней релаксации, используемый для нахождения поправок метода пространственного ребаланса.
В Приложении II приведены основные формулы DD0 схемы в HEX-Z геометрии.
В Приложении III описана 0-WDD схема и приведены ее основные формулы в X-Y-Z и R-cp-Z геометрии.
На защиту выносятся следующие результаты:
Предложен новый метод (S- процесс) ускорения сходимости внешних итераций при решении задач расчета эффективного коэффициента размножения и источника деления, возникающих при математическом моделировании ядерных реакторов методом дискретных ординат в трехмерной геометрии.
Исследована устойчивость линеаризованного приближения метода пространственного ребаланса на одногрупповой задаче с изотропным рассеянием в бесконечной однородной среде с помощью Фурье-анализа. Доказана его устойчивость.
Реализованы эффективные методы ускорения внешних итераций (S- процесс) и внутренних итераций (метод пространственного ребаланса) в программных модулях пакета «РЕАКТОР» для решения уравнения переноса нейтронов и гамма-квантов, используемого при математическом моделировании ядерных реакторов и их защиты в трехмерной геометрии. Реализованы взвешенные схемы метода дискретных ординат с коррекцией.
Проведены численные исследования на тестовых и реальных задачах эффективности реализованных методов ускорения сходимости внутренних и внешних итераций. Разработаны рекомендации по выбору параметра 8 для эффективного ускорения внешних итераций при расчетах в SnPm приближении. Разработаны рекомендации по
оптимальному выбору параметра в взвешенной 0-WDD схемы.
5. Выполнены трехмерные расчеты проектируемого в настоящее время энергетического
ядерного реактора СВБР 75/100, которые были использованы в проектных работах.
