Метод спектральной коллокации Чебышева для дифференциальных уравнений

Аппроксимация непериодических краевых задач

В вычислительной гидродинамике и квантовой физике мы часто применяем спектральные методы для решения дифференциальных уравнений. Как мы обсуждали ранее, методы, основанные на рядах Фурье, идеальны для периодических задач (где правый конец области гладко сшивается с левым). Но подавляющее большинство реальных инженерных задач непериодичны. Течение жидкости в трубе ограничено твердыми стенками, на которых скорость строго равна нулю (условие прилипания). Попытка применить ряды Фурье к такой задаче приводит к катастрофическому явлению Гиббса: на границах возникают незатухающие высокочастотные осцилляции, уничтожающие точность решения.

Для непериодических задач идеальным базисом являются ортогональные полиномы Чебышева. Метод спектральной коллокации Чебышева (Chebyshev Pseudospectral Method) берет лучшее от метода конечных разностей (простота работы в узлах) и спектральных методов (экспоненциальная сходимость). Ключевой особенностью метода является расположение расчетных узлов. Мы не можем использовать равномерную сетку из-за феномена Рунге (дикой раскачки полиномов на краях). Вместо этого используются точки Гаусса-Лобатто-Чебышева, которые сгущаются к границам отрезка по закону косинуса: x_j = cos(pi * j / N). Такое сгущение идеально компенсирует поведение многочленов, обеспечивая абсолютно стабильную и фантастически точную аппроксимацию.

Матрицы спектрального дифференцирования

Гениальность метода коллокации заключается в том, как он вычисляет производные. Если в методе конечных разностей мы берем значения функции в трех соседних точках и делим их на шаг h, то в спектральном методе Чебышева мы берем значения функции во всех N узлах сетки одновременно! Это глобальный подход.

Алгоритм заранее строит так называемую матрицу спектрального дифференцирования D_N. Эта матрица полностью заполнена (плотная). Чтобы найти производную функции во всех узлах сетки, мы просто умножаем вектор значений функции на эту матрицу. Элементы этой матрицы рассчитываются аналитически (или через Быстрое преобразование Фурье) и зависят от взаимного расположения узлов Чебышева. Если нам нужна вторая производная (для уравнения Пуассона), мы просто возводим матрицу D_N в квадрат! После подстановки граничных условий задача мгновенно сводится к решению СЛАУ. Скорость сходимости потрясает воображение: если решение является гладкой (аналитической) функцией, то ошибка убывает быстрее любой степени числа узлов. Там, где методу конечных разностей потребовались бы сотни тысяч точек, методу Чебышева достаточно 30-40 узлов для достижения машинной точности (10^-15). Этот алгоритм незаменим в теории гидродинамической устойчивости (уравнение Орра-Зоммерфельда) и физике плазмы.