ODETTA

CODE FEATURES

Для пространственной аппроксимации используется линейный разрывный метод конечных элементов (МКЭ) с аппроксимацией по Галеркину на тетраэдральных сеткаx.

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

Анизотропия рассеяния представляется разложением в ряд Фурье по полиномам Лежандра.

Для аппроксимации угловой зависимости используется метод дискретных ординат. Реализован широкий набор угловых квадратур (треугольных и квадратур типа произведения): Карлсона, Гаусса-Лежандра, Чебышева-Лежандра.

Обеспечена возможность задания нулевых (вакуумных), зеркальных, периодических или альбедных граничных условий.

В коде реализованы методы расчета нерассеянной компоненты излучения, направленные на снижение «лучевого» эффекта (методы первого и последнего столкновения). Кроме того, метод первого столкновения позволяет корректно моделировать точечные источники, а метод последнего столкновения – рассчитывать функционалы за пределами расчетной области в вакууме (без построения дополнительной сетки).

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

Основной объем расчетных данных хранится в формате HDF5.

Для построения неструктурированных тетраэдральных сеток и работы с CAD-моделями используется открытая интегрируемая платформа Salome.

Визуализация данных – SmartViewer или VisIt.

Программа для ЭВМ разработана на языке Fortran (с использованием современных стандартов Fortran 200x).

Функциональная схема работы кода ODETTA представлена на рисунке 1.

Рисунок 1 – Функциональная схема работы кода ODETTA

Подготовка входных данных состоит из трех блоков:

• процессинг CAD-модели, генерация неструктурированной сетки;
• подготовка групповых констант взаимодействия излучения со средой;
• подготовка конфигурационного файла, файла источника и файла на расчет линейных функционалов плотности потока излучения.

Выходные данные генерируются программой ODETTA в ходе расчета:

• дамп-файл в формате HDF5 (используется для расчета функционалов и для рестарта);
• лог-файл с информацией о ходе расчета, диагностическими сообщениями;
• out-файл с результаты расчета функционалов;
• vtk-файл для визуализации трехмерных распределений расчетных функционалов.