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

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

История

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

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

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

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

Уравнение неразрывности: \[ \nabla \cdot \mathbf{u} = 0 \]

Уравнение движения: \[ \rho \left( \frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} \right) = -\nabla p + \mu \nabla^2 \mathbf{u} + \mathbf{f} \]

Где:

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

Для сжимаемых жидкостей и газов уравнение неразрывности принимает вид: \[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0 \] и система дополняется уравнением состояния (например, уравнением Клапейрона — Менделеева для идеального газа) и уравнением энергии.

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

• Локальная производная \(\frac{\partial \mathbf{u}}{\partial t}\) — изменение скорости в данной точке со временем.
• Конвективная производная \((\mathbf{u} \cdot \nabla) \mathbf{u}\) — нелинейный член, описывающий перенос импульса самой жидкостью. Именно он ответственен за сложные явления, такие как турбулентность.
• Градиент давления \(-\nabla p\) — сила, вызванная перепадом давления.
• Диффузионный член \(\mu \nabla^2 \mathbf{u}\) — вязкое трение, стремящееся сгладить неоднородности в поле скорости.
• Внешние силы \(\mathbf{f}\) — гравитация, электромагнитные силы и т.д.

Свойства и классификация

Нелинейность

Уравнения Навье — Стокса являются нелинейными из-за конвективного члена \((\mathbf{u} \cdot \nabla) \mathbf{u}\). Это свойство приводит к тому, что сумма двух решений не является решением, и делает систему крайне сложной для аналитического решения в общем случае.

Вязкость

Наличие вязкости различает уравнения Навье — Стокса от уравнений Эйлера. При \(\mu = 0\) (идеальная жидкость) уравнения Навье — Стокса переходят в уравнения Эйлера. Вязкость приводит к диссипации кинетической энергии, превращению её в тепло.

Турбулентность

При высоких числах Рейнольдса (отношение инерционных сил к вязким) течение становится турбулентным. Турбулентность характеризуется хаотическими, нестационарными вихревыми движениями. Уравнения Навье — Стокса в принципе описывают турбулентность, но прямое численное моделирование (DNS) турбулентных течений требует огромных вычислительных ресурсов из-за необходимости разрешать все масштабы вихрей.

Применение

Уравнения Навье — Стокса используются в огромном количестве областей науки и техники:

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

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

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

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

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

Из-за отсутствия общих аналитических решений, уравнения Навье — Стокса решаются численно. Основные подходы включают:

• Метод конечных разностей: дискретизация производных на равномерной сетке.
• Метод конечных объёмов: наиболее распространён в коммерческих CFD-пакетах (ANSYS Fluent, OpenFOAM). Основан на интегральной форме законов сохранения.
• Метод конечных элементов: гибкий для сложных геометрий.
• Метод граничных элементов: используется для потенциальных течений.
• Прямое численное моделирование (DNS): решение без каких-либо моделей турбулентности, требует разрешения всех масштабов.
• Методы моделирования турбулентности: RANS (осреднение по Рейнольдсу), LES (моделирование крупных вихрей), DES (гибридный подход).

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

• Уравнения Навье — Стокса описывают не только жидкости, но и газы при условии, что длина свободного пробега молекул мала по сравнению с характерными размерами задачи.
• Существует точное аналитическое решение для течения Пуазёйля (ламинарное течение в трубе) и течения Куэтта (течение между двумя параллельными пластинами).
• В 2014 году российский математик Александр Боровков (НИУ ВШЭ) представил доказательство отсутствия сингулярностей для трёхмерных уравнений Навье — Стокса, однако его работа была раскритикована и не признана научным сообществом.
• Уравнения Навье — Стокса являются основой для многих компьютерных игр и спецэффектов в кино, где они используются для реалистичной симуляции жидкостей и дыма.

Источники

• Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. Том 6. Гидродинамика. — М.: Наука, 1986.
• Лойцянский Л. Г. Механика жидкости и газа. — М.: Дрофа, 2003.
• Флетчер К. Вычислительные методы в динамике жидкостей. — М.: Мир, 1991.
• Temam R. Navier-Stokes Equations: Theory and Numerical Analysis. — North-Holland, 1984.
• Математический институт Клэя. Проблемы тысячелетия: Уравнения Навье — Стокса.