Сравнение двух методов распараллеливания вычислений при решении интегро-дифференциального уравнения переноса излучения

Abstract

Рассматривается задача распараллеливания вычислений при решении интегро-дифференциального уравнения переноса излучения в сильно рассеивающих средах. Распараллеливание выполняется для двухшагового итерационного алгоритма решения системы сеточных уравнений. На первом шаге применяется метод простой итерации. На втором шаге к сеточным значениям, полученным на первом шаге, прибавляется ускоряющая сходимость итераций поправка. Уравнение для поправки решается методом подпространств Крылова. Выполнено сравнение эффективности двух методов распараллеливания двухшагового итерационного алгоритма. В BJ (Block Jacobi) методе расчет простой итерации выполняется локально в каждой пространственной подобласти. В BS (Block Seidel) методе выполняется сквозной расчет по всей области. Оба метода реализованы в программе РАДУГА Т решения уравнения переноса на неструктурированных сетках. Исследование эффективности методов выполнено на модели легководного реактора.