ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ И ПРОГРАММИРОВАНИЕ
Архив статей журнала
Разработан алгоритм высокоточного численного решения эллиптического уравнения второго порядка при наличии в области нескольких интерфейсов, в том числе пересекающихся и невыпуклых. Для аппроксимации задачи в окрестности интерфейсов используются нерегулярные ячейки (н-ячейки), отсекаемые ими от регулярных ячеек прямоугольной сетки, и законтурные части этих ячеек. Для построения приближенного решения предложено: 1) выписывать дополнительные условия согласования в н-ячейках на интерфейсах, увеличивая количество согласуемых ячеек вблизи интерфейсов; 2) уменьшать общую часть интерфейса, заключенную в соседних ячейках и используемую для записи условий. Для решения краевой задачи Дирихле реализован hp-вариант метода коллокации и наименьших квадратов (hp-МКНК) в сочетании с современными алгоритмами ускорения итерационного процесса: предобуславливание; распараллеливание с помощью OpenMP; ускорение, основанное на подпространствах Крылова; многосеточный алгоритм. При решении различных тестовых задач исследованы сходимость hp-МКНК и обусловленность возникающих переопределенных систем линейных алгебраических уравнений (СЛАУ). Проведено сравнение результатов, полученных МКНК, с результатами других авторов, использовавших метод MIB (англ. matched interface and boundary).