Методы фиктивных областей и погруженных границ (Immersed Boundary Method)
Проклятие подвижных сеток в гидродинамике
Классическая вычислительная гидродинамика (CFD) требует создания согласованной (конформной) сетки. Границы сетки должны строго совпадать с твердыми стенками обтекаемого объекта. Если мы рассчитываем течение воздуха вокруг неподвижного крыла самолета, создание такой сетки — сложная, но разовая задача. Однако что делать, если геометрия задачи непрерывно меняется? Как смоделировать поток крови через бьющийся сердечный клапан, падение тысяч камней в воду, вращение турбины или полет машущего крыла насекомого?
Если использовать классический метод конечных объемов, нам придется на каждом крошечном шаге по времени заново перестраивать всю вычислительную сетку (ремешинг), чтобы она подстраивалась под новые границы объекта. Перестроение объемной 3D-сетки — это невероятно медленная, тяжелая алгоритмическая задача, которая съедает до 90% всего процессорного времени, а сильные деформации ячеек быстро разрушают точность вычислений. Для решения этой фундаментальной проблемы механики сплошных сред был разработан радикальный подход — отказ от согласования сеток. Этот подход известен как Метод фиктивных областей (Fictitious Domain Method) и его знаменитая вариация — Метод погруженных границ (Immersed Boundary Method, IBM).
Идея метода погруженных границ Чарльза Пескина
Метод погруженных границ был изобретен в 1970-х годах выдающимся математиком Чарльзом Пескиным специально для моделирования потока крови в сердце. Гениальность идеи Пескина заключается в использовании двух абсолютно разных и независимых сеток, которые наложены друг на друга.
Первая сетка — эйлерова. Она представляет собой простую, неподвижную, равномерную и строгую прямоугольную (картезианскую) сетку, натянутую на весь объем пространства. В ней решаются уравнения течения жидкости (Навье-Стокса). Она никогда не меняется, что позволяет использовать сверхбыстрые решатели на базе БПФ. Вторая сетка — лагранжева. Это эластичная, подвижная сетка, состоящая только из одномерных кривых или двумерных поверхностей, которая описывает само движущееся тело (например, створку клапана). Эта сетка «плавает» сквозь неподвижную жидкостную сетку, совершенно не обращая внимания на ее узлы.
Взаимодействие через дельта-функции Дирака
Главный математический вызов метода — как заставить неподвижную жидкость почувствовать, что сквозь нее летит твердое тело или мембрана? Здесь на помощь приходит мощный инструмент функционального анализа — дельта-функция Дирака.
Поверхностная лагранжева сетка (тело) вычисляет упругие силы, возникающие при ее деформации. Затем эти поверхностные силы с помощью специальной сглаженной дельта-функции «размазываются» (интерполируются) на ближайшие узлы неподвижной жидкостной сетки, превращаясь в объемные источники силы в уравнениях Навье-Стокса. Жидкость, получив этот силовой толчок, меняет свое поле скоростей. Затем жидкость диктует свою новую скорость точкам мембраны, заставляя их сместиться на новый шаг по времени. Этот изящный алгоритм взаимодействия (Fluid-Structure Interaction, FSI) исключает необходимость перестроения сеток. Сегодня методы погруженных границ являются мощнейшим инструментом для моделирования гемодинамики, движения клеточных мембран, создания цифровых двойников кровеносной системы человека и расчетов сложной многофазной промышленной химии.
