Прямое численное моделирование

Прямое численное моделирование (ПЧМ, англ. Direct Numerical Simulation, DNS) — это метод численного решения уравнений Навье — Стокса, описывающих движение вязкой жидкости или газа, при котором все пространственные и временные масштабы турбулентности, от самых крупных (энергонесущих вихрей) до самых мелких (диссипативных, определяемых вязкостью), разрешаются непосредственно на расчётной сетке без использования каких-либо моделей турбулентности. В отличие от методов осреднённых по Рейнольдсу уравнений Навье — Стокса (RANS) или методов крупных вихрей (LES), ПЧМ не требует дополнительных допущений о структуре турбулентности, что делает его наиболее точным, но и наиболее вычислительно затратным подходом к моделированию течений.

Физические основы

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

В основе ПЧМ лежит численное решение полной, нестационарной системы уравнений Навье — Стокса для несжимаемой или сжимаемой среды. В декартовых координатах для несжимаемой жидкости (плотность \(\rho = \text{const}\)) уравнения имеют вид:

\[ \frac{\partial u_i}{\partial x_i} = 0 \]

\[ \frac{\partial u_i}{\partial t} + u_j \frac{\partial u_i}{\partial x_j} = -\frac{1}{\rho} \frac{\partial p}{\partial x_i} + \nu \frac{\partial^2 u_i}{\partial x_j \partial x_j} \]

где \(u_i\) — компоненты скорости, \(p\) — давление, \(\nu\) — кинематическая вязкость, \(t\) — время. ПЧМ решает эти уравнения без осреднения или фильтрации, разрешая все вихревые структуры вплоть до колмогоровского масштаба длины \(\eta = (\nu^3 / \varepsilon)^{1/4}\), где \(\varepsilon\) — скорость диссипации турбулентной энергии.

Требования к разрешению

Для корректного разрешения всех масштабов турбулентности необходимо, чтобы шаг расчётной сетки \(\Delta x\) был меньше или порядка колмогоровского масштаба \(\eta\). В трёхмерном случае количество узлов сетки \(N\) пропорционально числу Рейнольдса \(Re\) в степени 9/4:

\[ N \sim Re^{9/4} \]

Это означает, что при увеличении числа Рейнольдса вдвое объём вычислений возрастает примерно в 4,8 раза. Для инженерных приложений с высокими числами Рейнольдса (например, обтекание самолёта, \(Re \sim 10^7\)) требуемое количество узлов может достигать \(10^{15}\)–\(10^{18}\), что на несколько порядков превышает возможности современных суперкомпьютеров.

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

Ранние работы (1970–1980-е годы)

Первые попытки прямого численного моделирования турбулентности были предприняты в 1970-х годах. В 1972 году Стивен Орзаг (США) выполнил ПЧМ однородной изотропной турбулентности на сетке \(32^3\) узлов. В 1980-х годах, с появлением более мощных вычислительных машин, стали возможны расчёты канальных течений (Джон Ким, Парвиз Моин, Роберт Мозер, 1987) и пограничных слоёв.

Современный этап (1990-е — настоящее время)

С развитием суперкомпьютеров и параллельных вычислений ПЧМ стало применяться для изучения сложных течений: турбулентных струй, пламени, течений с теплообменом и химическими реакциями. В 2010-х годах были выполнены расчёты на сетках до \(10^{10}\) узлов (например, моделирование турбулентного пограничного слоя на стенке). В России исследования в области ПЧМ ведутся в Институте проблем механики РАН, Московском физико-техническом институте и Новосибирском государственном университете.

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

Дискретизация по пространству

Для ПЧМ используются высокоточные схемы, минимизирующие численную диссипацию и дисперсию:

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

Дискретизация по времени

Используются явные схемы Рунге — Кутты (чаще всего третьего или четвёртого порядка) или полунеявные схемы (например, метод Адамса — Башфорта для конвективных членов и метод Кранка — Николсона для вязких членов). Шаг по времени \(\Delta t\) ограничивается условием Куранта — Фридрихса — Леви (CFL): \(\Delta t \leq C \cdot \Delta x / u_{\text{max}}\).

Обработка граничных условий

Важнейшая задача — корректная постановка граничных условий, не вносящая нефизических возмущений. Для входа и выхода потока применяются методы характеристик или неотражающие граничные условия. Для стенок — условие прилипания (нулевая скорость).

Применение

Фундаментальные исследования

ПЧМ является основным инструментом для изучения механизмов турбулентности:

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

Инженерные приложения

Из-за высокой вычислительной стоимости ПЧМ редко применяется для прямого проектирования, но используется для:

• Верификации моделей турбулентности — сравнение результатов DNS с данными RANS и LES для калибровки эмпирических констант.
• Аэродинамика — моделирование обтекания профилей крыла при низких числах Рейнольдса (\(Re < 10^5\)).
• Теплообмен — изучение турбулентного переноса тепла в каналах и трубах.

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

ПЧМ применяется для моделирования мелкомасштабной турбулентности в атмосферном пограничном слое и верхнем слое океана, где характерные числа Рейнольдса относительно невелики.

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

Вычислительная стоимость

Основное ограничение ПЧМ — экспоненциальный рост требуемых вычислительных ресурсов с увеличением числа Рейнольдса. Для большинства инженерных задач (Re > 10^6) ПЧМ остаётся недоступным даже на самых мощных суперкомпьютерах. Например, для моделирования турбулентного течения в трубе диаметром 10 см при скорости 10 м/с (Re ≈ 10^5) требуется сетка порядка \(10^9\) узлов, а время счёта на современном кластере может составлять недели.

Чувствительность к граничным условиям

ПЧМ крайне чувствительно к заданию начальных и граничных условий. Небольшие погрешности на входе потока могут привести к существенному искажению результатов, особенно в пристеночных областях.

Ограниченность применимости

ПЧМ эффективно только для течений с умеренными числами Рейнольдса и простой геометрией. Для сложных трёхмерных конфигураций (например, обтекание автомобиля или самолёта) используются методы LES или гибридные подходы (DES, SAS).

Сравнение с другими методами

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

Экзафлопсные вычисления

С вводом в строй суперкомпьютеров экзафлопсной производительности (10^18 операций с плавающей запятой в секунду) в 2020-х годах (например, российский суперкомпьютер «Лобачевский» в Нижнем Новгороде) стало возможным выполнение ПЧМ для течений с числом Рейнольдса до 10^5–10^6 на сетках порядка \(10^{11}\) узлов.

Машинное обучение

Развиваются методы ускорения ПЧМ с помощью нейросетей: обучение моделей для предсказания мелкомасштабной турбулентности на основе данных DNS, что позволяет снизить требования к разрешению сетки.

Гибридные подходы

Комбинирование ПЧМ с методами LES или RANS в одной расчётной области (например, DNS вблизи стенки, LES в ядре потока) позволяет снизить вычислительные затраты при сохранении высокой точности в критических зонах.

Источники

• Моин П., Ма Дж. Прямое численное моделирование турбулентности. — М.: Мир, 1992. — 320 с.
• Попов С. А. Численное моделирование турбулентных течений. — М.: Физматлит, 2010. — 256 с.
• Фрик П. Г. Турбулентность: методы и модели. — Ижевск: Институт компьютерных исследований, 2010. — 336 с.
• Orszag S. A. Analytical theories of turbulence // Journal of Fluid Mechanics. — 1970. — Vol. 41, № 2. — P. 363–386.
• Kim J., Moin P., Moser R. Turbulence statistics in fully developed channel flow at low Reynolds number // Journal of Fluid Mechanics. — 1987. — Vol. 177. — P. 133–166.