Применение дробно-рациональных интерполяций для решения краевых задач с особенностями

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