Метод конечных объемов (МКО): фундамент вычислительной гидродинамики

Законы сохранения как основа физики сплошных сред

Ранее мы рассматривали метод конечных разностей (МКР) и метод конечных элементов (МКЭ). Однако, когда дело доходит до моделирования течений газов и жидкостей (аэродинамика самолетов, течения в газопроводах, взрывы, горение в двигателях), инженеры обращаются к третьему, специализированному инструменту — методу конечных объемов (МКО, или Finite Volume Method, FVM). Именно этот алгоритм заложен в основу большинства ведущих коммерческих и открытых пакетов вычислительной гидродинамики (CFD), таких как ANSYS Fluent, CFX и OpenFOAM.

Причина такого выбора кроется в самой сути дифференциальных уравнений газовой динамики (системы уравнений Эйлера или Навье-Стокса). Эти уравнения являются математической записью фундаментальных законов природы: сохранения массы, сохранения импульса и сохранения энергии. МКР и МКЭ оперируют узлами сетки и пытаются аппроксимировать дифференциальные операторы. При грубой сетке они могут математически «потерять» часть массы или энергии. Метод конечных объемов, напротив, строится не на дифференциальных, а на интегральных формах законов сохранения, гарантируя строгий баланс физических величин на любой, даже самой грубой сетке.

Контрольные объемы и расчет потоков на гранях

Суть МКО состоит в следующем: расчетная область со сложной геометрией разбивается на множество мелких, непересекающихся ячеек (многогранников). Эти ячейки называются контрольными объемами. Искомые величины (плотность, давление, скорость) хранятся не в узлах, а в центрах этих объемов и представляют собой средние значения внутри ячейки.

Поскольку масса, импульс или тепло не могут просто исчезнуть или появиться из ниоткуда, изменение этих величин внутри контрольного объема во времени может происходить только по одной причине: из-за перетекания (потока) этих величин через границы (грани) ячейки в соседние объемы или из них. Математически интеграл по объему преобразуется в интеграл по поверхности с помощью теоремы Остроградского-Гаусса. Таким образом, главная задача МКО сводится к максимально точному расчету потоков на гранях между ячейками. Что вытекло из ячейки А, то со стопроцентной точностью втекло в ячейку Б (консервативность метода).

Улавливание ударных волн и решатели Римана

Расчет потоков на гранях — это самая сложная часть МКО, особенно при сверхзвуковых течениях, где возникают разрывы плотности и давления (ударные волны). Обычные методы интерполяции данных из центров соседних ячеек на их общую грань приводят к нефизичным осцилляциям (колебаниям) вокруг ударной волны, которые разрушают вычисления.

Для решения этой проблемы российский математик Сергей Годунов в 1959 году предложил гениальную идею. Он предложил рассматривать разрыв параметров на гране между двумя ячейками как локальную задачу о распаде произвольного разрыва (задачу Римана в газовой динамике). Решая задачу Римана на каждой грани, алгоритм получает точные, физически обоснованные потоки, которые естественным образом учитывают направление распространения волн (скорость звука). Развитие этой идеи привело к созданию современных высокоразрешающих схем с ограничителями потоков (TVD-схемы, алгоритмы Roe, HLLC). Благодаря им МКО способен рассчитывать сложнейшие сверхзвуковые течения с резкими фронтами ударных волн без численных осцилляций, обеспечивая человечество инструментами для проектирования гиперзвуковых аппаратов и турбин.