ВОЗМОЖНОСТИ КОДАДля пространственной аппроксимации используется линейный разрывный метод конечных элементов (МКЭ) с аппроксимацией по Галеркину на тетраэдральных сеткаx.
Энергетическая зависимость рассматривается в групповом приближении, групповые уравнения решается методом простой итерации (метод итерации источника).
Анизотропия рассеяния представляется разложением в ряд Фурье по полиномам Лежандра.
Для аппроксимации угловой зависимости используется метод дискретных ординат. Реализован широкий набор угловых квадратур (треугольных и квадратур типа произведения): Карлсона, Гаусса-Лежандра, Чебышева-Лежандра.
Обеспечена возможность задания нулевых (вакуумных), зеркальных, периодических или альбедных граничных условий.
В коде реализованы методы расчета нерассеянной компоненты излучения, направленные на снижение «лучевого» эффекта (методы первого и последнего столкновения). Кроме того, метод первого столкновения позволяет корректно моделировать точечные источники, а метод последнего столкновения – рассчитывать функционалы за пределами расчетной области в вакууме (без построения дополнительной сетки).
Программа ODETTA поддерживает интерфейс OpenMP проведения параллельных вычислений для многопроцессорных систем с общей памятью. Распараллеливание кода осуществляется по угловым направлениям.
Основной объем расчетных данных хранится в формате HDF5.
Для построения неструктурированных тетраэдральных сеток и работы с CAD-моделями используется открытая интегрируемая платформа Salome.
Визуализация данных – SmartViewer или VisIt.
Программа для ЭВМ разработана на языке Fortran (с использованием современных стандартов Fortran 200x).
Функциональная схема работы кода ODETTA представлена на рисунке 1.
Рисунок 1 – Функциональная схема работы кода ODETTA
Подготовка входных данных состоит из трех блоков:
- процессинг CAD-модели, генерация неструктурированной сетки;
- подготовка групповых констант взаимодействия излучения со средой;
- подготовка конфигурационного файла, файла источника и файла на расчет линейных функционалов плотности потока излучения.
Выходные данные генерируются программой ODETTA в ходе расчета:
- дамп-файл в формате HDF5 (используется для расчета функционалов и для рестарта);
- лог-файл с информацией о ходе расчета, диагностическими сообщениями;
- out-файл с результаты расчета функционалов;
- vtk-файл для визуализации трехмерных распределений расчетных функционалов.
|