ПОВЫШЕНИЕ АЛГЕБРАИЧЕСКОГО ПОРЯДКА ТОЧНОСТИ ESn-КВАДРАТУРЫ
М. П. Пепеляев, Е. А. Ириничев Вопросы атомной науки и техники. Сер. Математическое моделирование физических процессов 2021. Вып.2. С. 11-23.
При решении задач переноса частиц в кинетическом приближении разностными методами возникает задача построения квадратурных формул на сфере (по угловым переменным). Одна из часто используемых квадратур - ESn-квадратура с равными весами. Равенство весов уменьшает погрешность квадратурной формулы. Однако ESn-квадратура обладает сравнительно невысоким алгебраическим порядком точности. Авторами разработан способ повышения точности квадратуры ESn по угловым переменным для решения уравнения переноса частиц. Равные веса квадратуры сохраняются, а направляющие косинусы полярного угла корректируются таким образом, чтобы выполнялись четные моментные соотношения. Система уравнений для вычисления направлений полета линеаризуется методом Ньютона и решается итерационно, с применением метода Гаусса на каждой итерации. Полученная угловая квадратура ESn обладает повышенным алгебраическим порядком точности по сравнению с ESn-квадратурой. Это подтверждается результатами проведенных численных исследований - расчетов интегралов от заданной функции по поверхности сферы, модельной задачи с анизотропным источником, симметричного теста К. Кобаяши (рис. - 6, табл. - 5, список лит. - 14). Ключевые слова: угловая квадратура, алгебраический порядок точности, равные веса, метод дискретных ординат, четные моментные соотношения, трехмерное уравнение переноса, декартова система координат.
Полный текст статьи
|