Метод функций уровня (Level Set Method) для отслеживания движущихся границ
Моделирование капель, пузырей и топологических парадоксов
В вычислительной гидродинамике, физике плазмы и материаловедении одной из самых завораживающих и вычислительно тяжелых задач является отслеживание многофазных течений с четкими подвижными границами. Как математически описать слияние двух мыльных пузырей, распад струи воды на тысячи капель, рост дендритов снежинки или распространение фронта лесного пожара? Трудность заключается не просто в движении поверхности, а в изменении ее топологии: одна область распадается на множество частей, а разные области сливаются воедино.
Исторически для этого применялись маркерные методы (MAC): на границу жидкости сажали миллионы виртуальных точек (частиц) и двигали их вместе с потоком. Однако при попытке двух пузырей слиться, точки пересекались, образуя математические «петли» и хаос. Алгоритмам приходилось запускать тяжелейшие эвристические проверки для удаления пересекающихся линий, что обрушивало скорость вычислений. Деформируемые сетки (Arbitrary Lagrangian-Eulerian) также терпели крах из-за полного вырождения ячеек в моменты слияния. Физике требовался принципиально новый математический инструмент, который обрабатывал бы изменение топологии автоматически и элегантно.
Повышение размерности: граница как нулевой срез
В 1988 году математики Стэнли Ошер и Джеймс Сетьян совершили революционный прорыв, разработав Метод функций уровня (Level Set Method). Их гениальная идея основывалась на геометрическом принципе повышения размерности. Вместо того чтобы пытаться явно отслеживать движущуюся двумерную поверхность в трехмерном пространстве (или 1D кривую на 2D плоскости), давайте введем вспомогательную, непрерывную функцию, размерность которой на единицу больше!
Эта глобальная функция (обозначается Phi) распределена по всему пространству. Внутри капли воды она принимает положительные значения (расстояние до границы), снаружи капли (в воздухе) — отрицательные. Сама же поверхность раздела сред, за которой мы охотимся, определяется исключительно как нулевая поверхность уровня (изоповерхность, где Phi = 0). Эта функция со знаком расстояния (Signed Distance Function) образует красивый непрерывный градиент во всем объеме.
Слияние без проблем и уравнение Гамильтона-Якоби
Движение границы раздела теперь сводится к решению дифференциального уравнения в частных производных конвективного типа (уравнения Гамильтона-Якоби) для функции Phi на неподвижной эйлеровой сетке. Скорость потока переносит значения функции Phi. Эволюция поверхности уровня происходит абсолютно естественно. Если две растущие капли соприкасаются, их нулевые поверхности просто сливаются в общую изоповерхность на уровне математики, без единой строчки специального кода для отслеживания столкновений!
Конечно, алгоритм таит в себе свои численные вызовы. При решении конвективного уравнения необходимы специальные односторонние (Upwind) разностные схемы или методы WENO высокого порядка, чтобы предотвратить нефизичные осцилляции. Кроме того, в процессе эволюции функция Phi теряет свойство точного геометрического расстояния. Поэтому периодически приходится запускать сложный алгоритм переинициализации (Reinitialization), который «расправляет» градиенты, сохраняя положение нулевого контура на месте. Наконец, стандартный метод Level Set страдает от небольшой потери физической массы при вычислениях. Для борьбы с этим его часто комбинируют с методом VOF (Volume of Fluid), получая гибридный алгоритм CLSVOF, который сегодня доминирует в промышленных CFD-кодах (Fluent, OpenFOAM) для симуляции кипения, волн в океане и расплавов металлов.
