Аппроксимации, сохраняющие локальный баланс массы и нейтронов деления в расчётах радиационной защиты

• Автор: Руссков, Александр Алексеевич
• Шифр специальности: 05.13.18
• Научная степень: Кандидатская
• Год защиты: 2010
• Место защиты: Москва
• Количество страниц: 168 с. : ил.

ОГЛАВЛЕНИЕ
Введение
Глава 1. Постановка задачи и обзор литературы
Глава 2. Использование метода трейсинга для конвертации комбинаторного представления геометрии в растровое (сеточное)
2.1. Принцип метода трейсинга
2.2. Реализация метода трейсинга для двумерной х,у геометрии
2.3. Реализация метода трейсинга для двумерной г,2 геометрии
2.4. Реализация метода трейсинга для двумерной г, 3 геометрии
2.5. Реализация метода трейсинга для трёхмерной х,у,г геометрии
2.6. Реализация метода трейсинга для трёхмерной г,3,г геометрии
2.7. Особенности вычисления материального состава ячеек
2.8. Корректировка луча при невозможности его проведения
2.9. Корректировка плотности материалов
2.10. Проверка корректности задания геометрии и средства визуализация геометрии
Глава 3. Подготовка комбинаторного источника на основе потвэльных и
покассетных данных о выгорании топлива
Введение
3.1. Постановка задачи
3.2. Учёт изменения множественности в процессе выгорания топлива
3.3. Алгоритм формирования комбинаторного источника
3.4. Расчёт плотности нейтронов деления, усреднённой по нескольким кампаниям
Г лава 4. Использование метода трейсинга для конвертации комбинаторного источника на сетку задачи
4.1. Введение и постановка задачи
4.2. Метод решения

Глава 5. Численные результаты использования подготовленных моделей геометрии и источника в расчётах активной зоны и радиационной защиты
РУ с ВВЭР
5.1 Расчёт нейтронных полей в активной зоне реактора с использованием VF метода. Эффект суперсходимости
5.2 Расчёт нейтронных полей в радиационной защите реактора с использованием VF метода
Глава 6. Реализация метода CADIS с использованием программ КАТРИН (детерминистический метод) и MCU (метод Монте-Карло)
6.1. Теоретические основы метода CADIS
6.2. Реализация алгоритма CADIS
6.3. Подготовка сопряжённого решения программы КАТРИН
6.4 Реализация и применение методики CADIS
Заключение
Благодарности
Литература
Приложение 1. Пакет начальных данных для 2D варианта программы
ConDat
Приложение 2. Пакет начальных данных для 3D варианта программы ConDat
Приложение 3. Описание бинарного формата пихтар для 20 геометрий

Приложение 4. Описание бинарного формата пихтар для 30 геометрий

Приложение 5. Формат текстового файла Ш геометрического описания

Приложение 6. Пример пускового пакета для 20 версии программы
СопОа!
Приложение 7. Пример пускового пакета для 30 версии программы ConDat

Приложение 8. Пакет начальных данных для программы ВигпОа!
Приложение 9. Пример пускового пакета для программы ВигпЕ)а1 для РУ
ВВЭР-1
Приложение 10. Пример пускового пакета для программы ВигпОа! для РУ
ВВЭР-
Приложение 11. Пакет начальных данных для конвертера СопЗоигсе... 154 Приложение 12. Пример пакета начальных данных для конвертера СопБоигсе
Приложение 13. Аппроксимация источника по аксиальной переменной г

Расчётной областью является множество точек с координатами, принадлежащими множествам точек разбиений. В случае декартовой системы расчётная область является прямоугольником или прямоугольным параллелепипедом, в случае полярной - кругом, сектором круга, кольцом, сектором кольца, в случае цилиндрической — аналогично полярной с отрезком по третьей координате. В случае г,2-системы расчётная область, по сути, является множеством точек из прямоугольника, образованного в плоскости Ш со всеми возможными значениями угла 3.
Методом трейсинга будем определять материальный состав для каждой ячейки, то есть доли простых материалов, из которых она состоит. В двумерном случае декартовой и полярной систем эти доли являются поверхностными долями, в трёхмерном случае — объёмными. В случае г,г-системы доли являются объёмными, но результат получается аналогично двумерному случаю с незначительным отличием, состоящим в введении множителей.
В методе трейсинга параллельно координатной оси Ох для декартовых х,у и х,у,г геометрий, и оси г для криволинейных г,г, г,3 и г,3,2 геометрий с необходимой плотностью проводятся лучи и вычисляются длины отрезков, относящиеся к различным материалам задачи. Исходя из прямолинейных участков траектории, лежащих в пределах области с одним материалом, определяется вклад каждого отрезка. Исходя из этой информации рассчитывается вклад каждого отрезка в пересекаемые ячейки разностной сетки, что иллюстрируется на Рис. 2.1.1 а) и б) для двумерных декартовой и полярной геометрий, соответственно. Отметим, что все вычисления выполняются с помощью арифметики двойной точности. Достаточность выбранной плотности траекторий алгоритма трассирования контролируется путем расчёта объемов каждого из материалов расчётной области.