Метод конечных объёмов

Метод конечных объёмов (МКО, англ. Finite Volume Method, FVM) — это численный метод решения дифференциальных уравнений в частных производных, основанный на интегральной формулировке законов сохранения. В отличие от методов конечных разностей (МКР) и конечных элементов (МКЭ), МКО оперирует не значениями функции в точках, а средними значениями по контрольным объёмам, что обеспечивает выполнение законов сохранения (массы, импульса, энергии) на дискретном уровне. Метод широко применяется в вычислительной гидродинамике (CFD), аэродинамике, теплофизике, электродинамике и других областях, где требуется моделирование потоков и переноса субстанций.

Основные принципы

Интегральная форма законов сохранения

МКО базируется на интегральной форме дифференциальных уравнений, которая для произвольного контрольного объёма \(V\) с границей \(S\) записывается как:

\[ \frac{\partial}{\partial t} \int_V \rho \phi \, dV + \oint_S \rho \phi \mathbf{v} \cdot \mathbf{n} \, dS = \oint_S \Gamma \nabla \phi \cdot \mathbf{n} \, dS + \int_V S_\phi \, dV \]

где:

• \(\rho\) — плотность;
• \(\phi\) — переносимая скалярная величина (например, температура, концентрация, компонента скорости);
• \(\mathbf{v}\) — вектор скорости;
• \(\Gamma\) — коэффициент диффузии;
• \(S_\phi\) — источниковый член;
• \(\mathbf{n}\) — внешняя нормаль к границе.

Это уравнение выражает баланс: изменение во времени внутри объёма + конвективный перенос через границу = диффузионный перенос через границу + генерация/поглощение внутри объёма.

Дискретизация расчётной области

Расчётная область разбивается на непересекающиеся контрольные объёмы (ячейки). Каждая ячейка представляет собой многогранник (в 3D) или многоугольник (в 2D). Сетка может быть:

• структурированной (прямоугольные или гексаэдральные ячейки, упорядоченные по индексам);
• неструктурированной (тетраэдры, призмы, полиэдры, треугольники, произвольные многоугольники);
• гибридной (комбинация структурированных и неструктурированных зон).

Для каждого контрольного объёма вводится узловая точка (центр ячейки), в которой хранятся значения искомых переменных. Среднее значение по объёму аппроксимируется значением в центре.

Аппроксимация потоков

Ключевой этап МКО — вычисление потоков через грани ячеек. Потоки делятся на:

• конвективные (перенос за счёт движения среды) — аппроксимируются с использованием схем интерполяции (например, противопоточная схема, схема центральных разностей, TVD-схемы, схемы высокого порядка);
• диффузионные (перенос за счёт градиентов) — аппроксимируются через градиенты на гранях, часто с использованием разностных формул и коррекции на неортогональность сетки.

Для устойчивости и точности применяются различные схемы дискретизации по времени: явные (Эйлера, Рунге-Кутты) и неявные (Эйлера, Кранка-Николсона, BDF).

История развития

Ранние работы

Основы МКО были заложены в 1950-х — 1960-х годах в работах по численному решению уравнений газовой динамики. Первые публикации, явно формулирующие метод как «конечно-объёмный», относятся к началу 1970-х годов. Важную роль сыграли работы:

• С. Патанкара (S. Patankar) и Д. Сполдинга (D. Spalding) — разработка метода SIMPLE (Semi-Implicit Method for Pressure Linked Equations) для решения уравнений Навье-Стокса;
• Б. Ломакса (B. Lomax) и Т. Пуллиама (T. Pulliam) — применение МКО в аэродинамике;
• Р. Хирша (C. Hirsch) — систематическое изложение метода в учебнике «Numerical Computation of Internal and External Flows».

Развитие в 1980-1990-х годах

В этот период МКО стал доминирующим методом в коммерческих CFD-пакетах (Fluent, CFX, STAR-CD, OpenFOAM). Были разработаны:

• методы решения связанных систем уравнений (SIMPLE, PISO, SIMPLEC);
• схемы высокого порядка для конвективных потоков (MUSCL, QUICK);
• адаптивные сетки и методы измельчения/разукрупнения ячеек;
• параллельные алгоритмы для распределённых вычислений.

Современное состояние

В XXI веке МКО активно развивается в направлениях:

• гибридные методы (связь с МКЭ, методом погружённых границ);
• высокоточные схемы (WENO, DG-FVM);
• многофазные и реагирующие течения (модели VOF, Euler-Euler, Euler-Lagrange);
• аэродинамика летательных аппаратов (включая дозвуковые, трансзвуковые и сверхзвуковые режимы);
• биомеханика (кровоток, дыхание);
• метеорология и океанология (модели общей циркуляции атмосферы и океана).

Классификация методов конечных объёмов

По типу решаемых уравнений

• МКО для уравнений Навье-Стокса (несжимаемые и сжимаемые течения);
• МКО для уравнений Эйлера (идеальный газ, без вязкости);
• МКО для уравнений теплопроводности и диффузии (параболические уравнения);
• МКО для уравнений Максвелла (электродинамика);
• МКО для уравнений переноса излучения (радиационный теплообмен).

По схеме дискретизации

• Противопоточная схема (upwind) — простая и устойчивая, но с численной диффузией;
• Схема центральных разностей (central difference) — более точная, но может быть неустойчивой при больших числах Пекле;
• Схемы высокого порядка (TVD, MUSCL, WENO) — минимизируют численную диффузию и осцилляции;
• Схемы с ограничителями (limiters) — для подавления нефизичных осцилляций вблизи разрывов.

По способу решения системы уравнений

• Явные схемы — просты в реализации, но требуют малого шага по времени (условие Куранта-Фридрихса-Леви);
• Неявные схемы — позволяют использовать большие шаги, но требуют решения систем линейных алгебраических уравнений (СЛАУ) на каждом шаге;
• Методы расщепления (SIMPLE, PISO, SIMPLEC) — итерационные процедуры для связанных систем «скорость-давление».

Применение

Вычислительная гидродинамика (CFD)

МКО является основным методом в большинстве коммерческих и открытых CFD-пакетов. Применяется для:

• моделирования обтекания автомобилей, самолётов, кораблей;
• расчёта течений в трубопроводах, теплообменниках, реакторах;
• анализа аэродинамики зданий и сооружений (ветровые нагрузки);
• моделирования процессов горения, смешения, химической кинетики.

Теплофизика и энергетика

• Расчёт температурных полей в твёрдых телах и жидкостях;
• Моделирование тепловых режимов электронных устройств (теплоотводы, радиаторы);
• Анализ систем охлаждения ядерных реакторов;
• Проектирование тепловых насосов и солнечных коллекторов.

Аэродинамика и авиакосмическая техника

• Расчёт подъёмной силы и лобового сопротивления крыла;
• Моделирование обтекания фюзеляжа, мотогондол, пилонов;
• Анализ трансзвуковых и сверхзвуковых течений (ударные волны, скачки уплотнения);
• Моделирование работы двигателей (внутренние течения, камеры сгорания).

Метеорология и океанология

• Модели общей циркуляции атмосферы (например, WRF, COSMO);
• Модели океанических течений (ROMS, MITgcm);
• Прогноз погоды и климатические модели.

Биомеханика и медицина

• Моделирование кровотока в сосудах (аневризмы, стенозы);
• Анализ дыхания в лёгких (вентиляция, перенос аэрозолей);
• Моделирование движения спинномозговой жидкости.

Преимущества и недостатки

Преимущества

• Сохранение законов — интегральная форма гарантирует выполнение законов сохранения на дискретном уровне (свойство консервативности);
• Гибкость — работает на неструктурированных и адаптивных сетках, позволяет моделировать сложные геометрии;
• Устойчивость — при правильном выборе схем дискретизации метод устойчив даже при больших градиентах;
• Широкая применимость — от дозвуковых до гиперзвуковых течений, от однофазных до многофазных сред.

Недостатки

• Численная диффузия — особенно при использовании простых противопоточных схем на грубых сетках;
• Сложность реализации — требуется аккуратная аппроксимация потоков, коррекция на неортогональность, обработка граничных условий;
• Вычислительная стоимость — для неявных схем требуется решение больших СЛАУ, что может быть ресурсоёмко;
• Ограничения для задач с разрывами — требует специальных схем (TVD, WENO) для подавления осцилляций.

Связь с другими методами

Метод конечных разностей (МКР)

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

Метод конечных элементов (МКЭ)

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

Метод погружённых границ (IBM)

IBM позволяет моделировать движение тел в потоке без перестроения сетки, накладывая граничные условия через источниковые члены. МКО часто используется в комбинации с IBM для задач с подвижными границами.

Программные реализации

Коммерческие пакеты

• ANSYS Fluent — один из наиболее распространённых CFD-пакетов, использующий МКО;
• ANSYS CFX — также основан на МКО, специализируется на турбомашинах;
• STAR-CCM+ — пакет от Siemens, использует МКО на полиэдральных сетках;
• COMSOL Multiphysics — использует МКЭ, но включает модули на основе МКО для некоторых задач.

Открытые пакеты

• OpenFOAM — открытая библиотека для CFD, полностью основанная на МКО, поддерживает неструктурированные сетки, параллельные вычисления, широкий набор физических моделей;
• SU2 — открытый пакет для аэродинамики и оптимизации, использует МКО;
• Code_Saturne — разработан EDF (Франция), применяется в энергетике;
• Gerris (ныне Basilisk) — открытый пакет с адаптивными сетками, основан на МКО.

Интересные факты

• МКО является основой для большинства CFD-расчётов в авиационной промышленности, включая проектирование самолётов Boeing 787 и Airbus A350.
• В OpenFOAM реализовано более 200 решателей для различных физических задач (от несжимаемых течений до магнитогидродинамики).
• МКО используется в моделях общей циркуляции атмосферы, таких как WRF (Weather Research and Forecasting), которые применяются для прогноза погоды в России и других странах.
• В России МКО активно применяется в Институте проблем механики РАН, МФТИ, МГУ, ЦАГИ, а также в отраслевых институтах (например, ВИАМ, НИЦ «Курчатовский институт»).

Источники

• Patankar S. V. Numerical Heat Transfer and Fluid Flow. — Hemisphere Publishing Corporation, 1980.
• Ferziger J. H., Perić M. Computational Methods for Fluid Dynamics. — 3rd ed. — Springer, 2002.
• Hirsch C. Numerical Computation of Internal and External Flows. — 2nd ed. — Butterworth-Heinemann, 2007.
• Versteeg H. K., Malalasekera W. An Introduction to Computational Fluid Dynamics: The Finite Volume Method. — 2nd ed. — Pearson, 2007.
• OpenFOAM User Guide. — The OpenFOAM Foundation, 2023.
• Weller H. G., Tabor G., Jasak H., Fureby C. A tensorial approach to computational continuum mechanics using object-oriented techniques // Computers in Physics. — 1998. — Vol. 12, No. 6. — P. 620–631.