Псевдоспектральные методы решения нелинейных эволюционных уравнений
Нелинейные волны: от солитонов до оптоволокна
В современной физике существуют удивительные явления, поведение которых невозможно предсказать с помощью линейных теорий. Самым известным примером являются солитоны — уединенные волны, которые распространяются в нелинейной среде, не меняя своей формы, и проходят друг сквозь друга без разрушения (как частицы). Солитоны описываются нелинейным уравнением Кортевега-де Фриза (КдФ) на мелкой воде, а распространение лазерных импульсов в трансатлантических оптоволоконных кабелях (основа современного интернета) подчиняется нелинейному уравнению Шредингера (НУШ).
Эти эволюционные уравнения содержат две конкурирующие физические силы: линейную дисперсию (которая заставляет волну расплываться в пространстве) и нелинейность (которая заставляет волну сжиматься, делая фронт более крутым). Точный баланс этих двух сил и порождает стабильный солитон. Попытка решить такие уравнения классическим методом конечных разностей часто заканчивается провалом: искусственная численная дисперсия разностных схем (ошибки аппроксимации производных) вмешивается в тонкий физический баланс уравнения, и численный солитон либо разваливается, либо порождает нефизичную рябь.
Псевдоспектральный подход: лучшее из двух миров
Для достижения необходимой экспоненциальной точности в расчете пространственных производных (чтобы свести численную дисперсию к машинному нулю) математики применяют спектральные методы Фурье. Однако чистый спектральный метод Галеркина сталкивается с непреодолимой алгоритмической стеной при обработке нелинейных членов (например, когда функция умножается сама на себя). Вычисление произведения двух функций в спектральном пространстве требует операции свертки, которая имеет вычислительную сложность O(N^2), что делает метод катастрофически медленным.
Гениальным выходом из ситуации стали псевдоспектральные методы (Pseudospectral methods, или методы Фурье-коллокации). Их философия проста: «Выполняй каждую математическую операцию в том пространстве, где она алгоритмически дешевле всего». В этом методе вычисление пространственных производных происходит в частотном пространстве (где это тривиальное и точное умножение на волновое число). Но как только дело доходит до нелинейного члена, алгоритм вызывает Быстрое преобразование Фурье (БПФ) и мгновенно переносит данные в физическое пространство!
Расщепление по физическим процессам (Split-Step)
В физическом (координатном) пространстве нелинейная операция — это простейшее поточечное умножение чисел, которое выполняется за O(N) операций. Затем алгоритм снова применяет обратное БПФ, возвращаясь в спектральный мир для продолжения интегрирования. Этот феноменально быстрый перескок между двумя пространствами благодаря БПФ снижает общую сложность решения до O(N log N).
Особую популярность приобрел Псевдоспектральный метод с расщеплением по физическим процессам (Split-Step Fourier Method). Уравнение искусственно разбивается на две независимые части: чисто линейную (отвечающую за дисперсию) и чисто нелинейную. Алгоритм на каждом временном шаге делает микро-шаг интегрирования только с учетом линейной части в спектральном пространстве, затем перепрыгивает в физическое пространство и делает микро-шаг с учетом только нелинейности, после чего складывает результаты (схема Ли-Троттера или симметричная схема Стрэнга). Этот подход стал абсолютным мировым стандартом в нелинейной волоконной оптике, фотонике и физике плазмы, позволяя рассчитывать динамику сотен взаимодействующих солитонов с точностью до 15 знаков.