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

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

История

Уравнения были выведены в первой половине XIX века. Французский инженер и физик Клод-Луи Навье в 1822 году опубликовал работу, в которой впервые представил уравнения движения вязкой жидкости, основываясь на молекулярной теории. В 1845 году британский математик и физик сэр Джордж Габриель Стокс дал строгий вывод уравнений, используя континуальную механику и гипотезу о линейной связи между тензором напряжений и тензором скоростей деформации (закон Ньютона для вязкости). Работы Стокса завершили формирование классической формы уравнений.

Математическая формулировка

Уравнения Навье-Стокса представляют собой запись второго закона Ньютона для элементарного объёма жидкости. В векторной форме для несжимаемой жидкости они имеют вид:

\[ \rho \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} \right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \mathbf{f} \]

где:

• \(\rho\) — плотность жидкости,
• \(\mathbf{v}\) — вектор скорости,
• \(t\) — время,
• \(p\) — давление,
• \(\mu\) — динамическая вязкость,
• \(\mathbf{f}\) — вектор плотности объёмных сил (например, силы тяжести),
• \(\nabla\) — оператор набла,
• \(\nabla^2\) — оператор Лапласа.

Дополняется уравнение уравнением неразрывности (законом сохранения массы): \[ \nabla \cdot \mathbf{v} = 0 \]

Для сжимаемых жидкостей (газов) уравнение неразрывности содержит производную плотности по времени, а к уравнениям добавляется уравнение состояния и уравнение энергии.

Физический смысл членов уравнения

• Левая часть (\(\rho \frac{D\mathbf{v}}{Dt}\)) — субстанциальная (полная) производная скорости по времени, описывающая ускорение частицы жидкости. Включает локальное ускорение (изменение скорости во времени в данной точке) и конвективное ускорение (изменение скорости при перемещении частицы в пространстве).
• Правая часть:
• \(-\nabla p\) — градиент давления (сила, обусловленная разностью давлений),
• \(\mu \nabla^2 \mathbf{v}\) — вязкое трение (диффузия импульса),
• \(\mathbf{f}\) — внешние объёмные силы.

Свойства и сложность

Уравнения Навье-Стокса являются нелинейными из-за конвективного члена \((\mathbf{v} \cdot \nabla) \mathbf{v}\). Нелинейность приводит к тому, что:

• Решения могут быть неустойчивыми и чувствительными к начальным условиям.
• Возможны турбулентные режимы течения, характеризующиеся хаотическим изменением скорости и давления.
• Аналитические решения существуют лишь для ограниченного числа простых случаев (например, течение Пуазёйля в трубе, течение Куэтта между параллельными пластинами).

В большинстве практических задач уравнения решаются численно методами вычислительной гидродинамики (CFD).

Применение

Уравнения Навье-Стокса используются в:

• Авиа- и ракетостроении: расчёт обтекания крыльев, фюзеляжей, сопел двигателей.
• Автомобилестроении: моделирование аэродинамики автомобиля, охлаждения двигателя.
• Метеорологии и океанологии: прогноз погоды, моделирование течений в океане, атмосферной циркуляции.
• Биомеханике: расчёт кровотока в сосудах, движения воздуха в лёгких.
• Энергетике: проектирование турбин, насосов, теплообменников.
• Нефтегазовой отрасли: моделирование течения нефти и газа в пластах и трубопроводах.

Проблема существования и гладкости решений

Вопрос о том, всегда ли для трёхмерных уравнений Навье-Стокса существует гладкое (бесконечно дифференцируемое) решение при любых гладких начальных данных, является одной из семи «задач тысячелетия» Математического института Клэя. За доказательство или опровержение этого утверждения назначен приз в 1 миллион долларов США. На данный момент (2024 год) проблема остаётся нерешённой. Для двумерного случая существование и гладкость решений доказаны. Для трёхмерного случая известны результаты о существовании слабых решений (в смысле Лере — Хопфа), но их гладкость не установлена.

Численные методы решения

Для практического применения разработаны методы дискретизации уравнений:

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

Для учёта турбулентности применяются:

• Прямое численное моделирование (DNS) — решение без моделей турбулентности, требует огромных вычислительных ресурсов.
• Моделирование крупных вихрей (LES) — моделирование крупномасштабных структур, мелкие масштабы моделируются.
• Усреднённые по Рейнольдсу уравнения Навье-Стокса (RANS) — осреднение по времени, введение моделей турбулентности (например, k-ε, k-ω SST).

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

• Уравнения Навье-Стокса являются основой для симуляции жидкостей в компьютерной графике и спецэффектах в кино.
• При нулевой вязкости (\(\mu = 0\)) уравнения переходят в уравнения Эйлера для идеальной жидкости.
• В 2014 году математик Теренс Тао показал, что для предотвращения образования сингулярностей в трёхмерных уравнениях достаточно, чтобы скорость роста решения была ограничена определённым образом, но полного доказательства не получил.
• В 2022 году группа исследователей из Университета Брауна опубликовала работу, в которой утверждалось, что обнаружена новая особенность в решении уравнений, связанная с «квази-сингулярностями», но окончательного вывода о существовании гладких решений сделано не было.

Источники

• Ландау Л. Д., Лифшиц Е. М. Гидродинамика. — М.: Наука, 1986.
• Лойцянский Л. Г. Механика жидкости и газа. — М.: Дрофа, 2003.
• Batchelor G. K. An Introduction to Fluid Dynamics. — Cambridge University Press, 2000.
• Fefferman C. L. Existence and smoothness of the Navier–Stokes equation // The Millennium Prize Problems. — Clay Mathematics Institute, 2006.
• Pope S. B. Turbulent Flows. — Cambridge University Press, 2000.