Ме́тод подві́йної сфе́ри Фур'є́ (англ. Double Fourier Sphere method, DFS) — у математиці це метод, що дає змогу перетворити функцію, задану на сферичній поверхні, на функцію, в прямокутній області із збереженням періодичності як за довготою, так і за широтою.

Область застосування Представлення та чисельне маніпулювання функціями на різних геометричних многовидах є важливим для таких галузей

Прогнозування погоди та глобальне кліматичне моделювання; Геофізика, геодезія та науки про Землю; Астрофізика; Молекулярна біологія (зокрема алгоритм білкового докінгу); Вивчення динаміки активних рідин і біологічних системах.

Математичні підходи стикались з труднощами при апроксимації на многовидах до появи методу DFS:

Атласні методи вимагали складного гладкого "зшивання" локальних карт (поверхню розрізали на шматки, які потім було надзвичайно важко гладко з'єднати назад без помилок); Методи об'ємної аппроксимації збільшували розмірність задачі (рахували не тільки поверхню, а й увесь простір навколо неї, що сильно перевантажувало комп'ютер); Спеціалізовані базиси (наприклад, сферичні гармоніки) призводили до чисельнних сингулярностей на полюсах (створювали математичні аномалії й помилки на полюсах, де лінії координат сходяться в точку) і поступались за швидкістю ШПФ (Швидкого перетворення Фур'є). Метод DFS усуває ці недоліки, зводячи задачу до класичного аналізу Фур'є на торі.

Вступ

Спочатку функція на сфері радіуса записується як , використовуючи сферичні координати, тобто через зв'язок із декартовими координатами:

Для таких обчислювальних задач стандартним підходом є використання координат довготи та широти , де та позначають азимутальний (довгота) та полярний (коширота) кути відповідно. Тобто:

За такого представлення функція є -періодичною за змінною , але не є періодичною за , через що періодичність у напрямку широти втрачається. Щоб відновити її, функцію «подвоюють» і визначають пов’язану функцію на розширеній прямокутній області , яка компактно задається як:

де , якщо , та , якщо .

Для детального аналізу в кожному з чотирьох квадрантів розширеної області отриману функцію можна деталізувати через допоміжні функції і , визначені на початковому домені :

Завдяки такому штучному розширенню області функція стає повністю -періодичною за обома координатами ( та ), а також є сталою вздовж ліній і , що відповідають полюсам сфери.

Це дозволяє розкласти функцію у стандартний подвійний (двовимірний) ряд Фур'є:

Історія Метод був запропонований Філіпом Мерілісом (Philip E. Merilees) у 1973 році для розв’язання рівнянь мілкої води на сфері. Значний розвиток отримав у 1974 році завдяки Стівену Орзагу (Steven A. Orszag) у класичній статті «Fourier Series on Spheres».

У перші десятиліття метод використовувався відносно мало через домінування спектральних методів на основі сферичних гармонік, для яких існували добре розвинені обчислювальні бібліотеки. Ситуація змінилася у 1990-х роках завдяки роботам Бенгта Форнберга (Bengt Fornberg), який детально проаналізував умови гладкості та періодичності, необхідні для коректної поведінки методу на полюсах.

У 2010–2020-х роках метод DFS пережив новий підйом завдяки зростанню обчислювальних потужностей і потребі в моделях надвисокої роздільної здатності. Особливо активно він застосовується у чисельній відносності, моделюванні чорних дір, сучасному кліматичному моделюванні та спектральних методах. Головною перевагою залишається можливість використання швидкого перетворення Фур'є (FFT), що робить метод надзвичайно ефективним на сучасних багатоядерних і GPU-архітектурах.

Посилання