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

Вычислительная гидродинамика (англ. Computational Fluid Dynamics, CFD) — раздел механики сплошных сред, объединяющий математическое моделирование, численные методы и программное обеспечение для расчёта течений жидкостей и газов, а также сопутствующих процессов тепло- и массообмена, химических реакций и деформаций твёрдых тел. Основная цель CFD — получение количественного описания гидродинамических и тепловых полей в заданной области пространства путём численного решения системы уравнений Навье — Стокса (для вязких сред) или уравнений Эйлера (для невязких сред), дополненных уравнениями неразрывности, энергии и состояния.

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

Ранние этапы

Первые попытки численного решения гидродинамических задач относятся к началу XX века. В 1910 году британский метеоролог Льюис Фрай Ричардсон предложил метод численного прогноза погоды, основанный на конечно-разностной аппроксимации уравнений атмосферной динамики. Однако из-за отсутствия вычислительной техники расчёты выполнялись вручную и требовали огромных затрат времени. В 1928 году Рихард Курант, Курт Фридрихс и Ганс Леви разработали условие устойчивости (CFL-условие) для явных разностных схем, ставшее фундаментом численного моделирования.

Эра электронных вычислительных машин

С появлением первых ЭВМ в 1950-х годах началось активное развитие CFD. В 1953 году американские учёные Джон фон Нейман и Роберт Рихтмайер опубликовали работу по численному решению уравнений газовой динамики с использованием метода искусственной вязкости. В 1960-х годах были разработаны первые коммерческие коды: в 1963 году компания General Motors создала программу для расчёта обтекания автомобилей, а в 1967 году появился пакет PHOENICS, ставший предшественником современных CFD-систем.

Современный этап

С 1980-х годов CFD стала доступна для широкого круга инженеров благодаря росту производительности компьютеров и появлению графических интерфейсов. В 1990-х годах были разработаны методы крупных вихрей (LES) и моделирования отсоединённых вихрей (DES), а в 2000-х — методы прямого численного моделирования (DNS) для задач турбулентности. Сегодня CFD применяется в авиации, автомобилестроении, энергетике, медицине, метеорологии и других отраслях.

Математическая основа

Уравнения Навье — Стокса

Центральным объектом CFD является система уравнений Навье — Стокса, описывающая движение вязкой ньютоновской жидкости. В векторной форме для несжимаемой жидкости она имеет вид:

• Уравнение неразрывности: ∇·u = 0
• Уравнение импульса: ∂u/∂t + (u·∇)u = -∇p/ρ + ν∇²u + f

где u — вектор скорости, p — давление, ρ — плотность, ν — кинематическая вязкость, f — вектор внешних сил.

Уравнение энергии

Для задач с теплообменом добавляется уравнение энергии, учитывающее конвекцию, теплопроводность и источники тепла. В случае сжимаемых течений уравнение энергии связывается с уравнением состояния (например, идеального газа p = ρRT).

Модели турбулентности

Поскольку прямое численное моделирование турбулентных течений требует огромных вычислительных ресурсов, на практике используются осреднённые по Рейнольдсу уравнения Навье — Стокса (RANS) с дополнительными моделями турбулентности. Наиболее распространённые модели:

• k-ε (стандартная, RNG, реализуемая)
• k-ω (стандартная, SST)
• Модель переноса сдвиговых напряжений (SST)
• Модели напряжений Рейнольдса (RSM)

Численные методы

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

Простейший метод, основанный на аппроксимации производных разностными отношениями на регулярной сетке. Применяется для задач с простой геометрией, но ограничен в работе с криволинейными границами.

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

Наиболее популярный в коммерческих CFD-пакетах. Расчётная область разбивается на контрольные объёмы, для каждого из которых интегрируются уравнения сохранения. Метод обеспечивает консервативность потоков и хорошо работает на неструктурированных сетках.

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

Используется для задач с деформируемыми границами и многофазных течений. Основан на разбиении области на элементы и поиске приближённого решения в виде линейной комбинации базисных функций.

Спектральные методы

Применяются для задач с высокой гладкостью решения (например, в метеорологии). Решение представляется в виде ряда по ортогональным функциям (полиномы Чебышёва, тригонометрические функции).

Программное обеспечение

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

• ANSYS Fluent (разработчик — компания ANSYS, США) — один из наиболее распространённых CFD-кодов, поддерживает широкий спектр моделей турбулентности и многофазных течений.
• STAR-CCM+ (разработчик — Siemens Digital Industries Software, Германия) — интегрированная платформа для многодисциплинарного моделирования.
• COMSOL Multiphysics (разработчик — COMSOL Group, Швеция) — пакет, объединяющий CFD с расчётами электромагнитных полей, акустики и механики деформируемого твёрдого тела.
• OpenFOAM (разработчик — OpenCFD Ltd, Великобритания) — открытая платформа на языке C++, позволяющая создавать пользовательские решатели.

Отечественные разработки

В Российской Федерации разработаны и используются несколько CFD-пакетов:

• Логос (разработчик — Российский федеральный ядерный центр — ВНИИЭФ, г. Саров) — предназначен для расчёта аэродинамики, тепломассообмена и прочности в ракетно-космической и атомной отраслях.
• FlowVision (разработчик — компания «ТЕСИС», г. Москва) — используется в авиастроении, судостроении и энергетике.
• VP2/3 (разработчик — Институт проблем механики РАН) — специализированный пакет для моделирования течений с химическими реакциями.

Применение

Авиация и космонавтика

CFD используется для расчёта обтекания крыльев, фюзеляжей, лопаток турбин и сопел ракетных двигателей. Например, при проектировании российского истребителя Су-57 применялись CFD-расчёты для оптимизации аэродинамических характеристик.

Автомобилестроение

Моделирование внешней аэродинамики (снижение лобового сопротивления), системы охлаждения двигателя, обтекания кузова и тормозных механизмов. Компании «АвтоВАЗ» и «КАМАЗ» используют CFD для доводки новых моделей.

Энергетика

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

Медицина

Моделирование кровотока в сосудах (для диагностики атеросклероза), воздушных потоков в дыхательных путях, движения лекарственных аэрозолей. В России такие расчёты проводятся в Сеченовском университете и НМИЦ им. В. А. Алмазова.

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

Численный прогноз погоды (модели WRF, COSMO-Ru), расчёт течений в Мировом океане (модель NEMO), моделирование распространения загрязнителей в атмосфере.

Ограничения и критика

Вычислительные ресурсы

Прямое численное моделирование турбулентности требует суперкомпьютеров с производительностью в десятки петафлопс. Для большинства практических задач приходится использовать осреднённые модели, которые могут давать погрешности до 20–30% в сложных течениях (отрывные зоны, закрученные потоки).

Верификация и валидация

Качество CFD-расчёта сильно зависит от выбора модели турбулентности, типа сетки, граничных условий. Для достоверности результатов требуется сравнение с экспериментальными данными или аналитическими решениями. В ряде отраслей (например, в авиации) CFD-результаты не могут быть единственным основанием для сертификации изделия.

Проблемы сходимости

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

Перспективы развития

Основные направления совершенствования CFD включают:

• Разработку гибридных моделей турбулентности (например, DES, IDDES)
• Использование методов машинного обучения для ускорения расчётов и построения суррогатных моделей
• Интеграцию CFD с расчётами прочности (FSI — Fluid-Structure Interaction) и акустики
• Переход к экзафлопсным вычислениям, позволяющим проводить DNS для промышленных задач

Источники

• Андерсон Д. А., Таннехилл Дж. К., Плетчер Р. Х. Вычислительная гидродинамика и тепломассообмен. — М.: Мир, 1990.
• Флетчер К. Вычислительные методы в динамике жидкостей. — М.: Мир, 1991.
• Ferziger J. H., Perić M. Computational Methods for Fluid Dynamics. — Springer, 2002.
• Wilcox D. C. Turbulence Modeling for CFD. — DCW Industries, 2006.
• Отчёт о научно-исследовательской работе «Развитие методов вычислительной гидродинамики для задач авиационной техники» (ЦАГИ, 2021).