Теория переноса излучений Ф8-01Н Московский инженерно-физический

уравнений в матричной форме.
Организация итерационного процесса.
Проблема сходимости численных

схем.
Улучшенные итерационные методы.
Внутренние и внешние итерации.

в матричной форме в виде:

=

Диагональные компоненты матрицы

положительны, в то время как недиагональные члены - отрицательны или равны нулю. Сумма недиагональных элементов в любом данном ряду меньше, чем диагональный элемент. Таким образом, матрица является неприводимой диагонально преобладающей. Следовательно, для матрицы существует обратная матрица , и решение уравнения можно записать в виде:

=

суммы трех матриц:
=
где – диагональная матрица

(отличные от нуля элементы находятся только на основной диагонали), – верхняя треугольная матрица (отличные от нуля элементы находятся только выше основной диагонали) и – нижняя треугольная матрица (отличные от нуля элементы находятся ниже основной диагонали).
= ( + ) +
Итерационный процесс можно определить следующим образом:

= ( + ) +

пор, пока разность между потоками и

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

в правой части уравнения будут использоваться только значения потока

из последней итерации, т. е. . Может оказаться, что после того, как рассчитана новая компонента , более предпочтительно использовать именно ее, а не для определения последующих компонент :

( – ) = +

Так как матрица ( – ) треугольная, включая основную диагональ, то можно легко найти обратную ей или решить уравнение относительно .

и внешние итерации, основана на идее вычисления компонент

на базе только вычисленных компонент и на внутренних итерациях. На внешних итерациях производится пересчет источника с учетом всех вычисленных . Часто на внутренних итерациях решается уравнение с фиксированным источником деления, а полученное решение в итерациях по рассеянию используется для пересчета источника деления.