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

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

¶История

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

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

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

Уравнение неразрывности (сохранение массы):

∇·u = 0 где u — вектор скорости жидкости.

Уравнение движения (сохранение импульса):

∂u/∂t + (u·∇)u = - (1/ρ) ∇p + ν ∇²u + F где:

t — время;
ρ — плотность жидкости;
p — давление;
ν — кинематическая вязкость (коэффициент пропорциональности между напряжением и скоростью деформации);
F — вектор внешних сил (например, сила тяжести).

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

¶Варианты записи

Для сжимаемой жидкости: уравнение неразрывности принимает вид ∂ρ/∂t + ∇·(ρu) = 0, а в уравнение движения добавляются члены, связанные с объёмной вязкостью.
В проекциях на оси координат: в декартовой системе координат (x, y, z) каждое из уравнений распадается на три скалярных уравнения для компонент скорости (u, v, w).
В безразмерной форме: вводится число Рейнольдса (Re = UL/ν, где U — характерная скорость, L — характерный размер), которое определяет соотношение инерционных и вязких сил. При малых Re (ламинарное течение) вязкие силы доминируют, при больших Re — инерционные, что приводит к турбулентности.

¶Физический смысл и свойства

Уравнения Навье-Стокса описывают баланс сил в каждой точке жидкости. Нелинейный член (u·∇)u отвечает за перенос импульса движущейся жидкостью и является причиной сложного, хаотического поведения течений при высоких числах Рейнольдса. Вязкий член ν∇²u описывает диссипацию энергии за счёт внутреннего трения, превращающую механическую энергию в тепло.

Ключевые свойства:

Нелинейность: делает систему чрезвычайно сложной для аналитического решения, за исключением некоторых частных случаев (например, течение Пуазёйля, течение Куэтта).
Эллиптико-параболический тип: сочетает черты параболических (диффузия) и эллиптических (несжимаемость) уравнений.
Чувствительность к начальным условиям: при больших числах Рейнольдса малые возмущения могут приводить к кардинально разным сценариям развития течения (турбулентность).

¶Применение

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

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

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

Одной из важнейших нерешённых проблем современной математики является вопрос о существовании и гладкости (отсутствии разрывов) решений трёхмерных уравнений Навье-Стокса для всех моментов времени при произвольных гладких начальных данных. Эта задача входит в список «Проблем тысячелетия» Математического института Клэя (США), за её решение объявлена премия в 1 миллион долларов США. Формально проблема формулируется так: доказать, что для любого гладкого начального поля скорости и любой внешней силы существует единственное гладкое (бесконечно дифференцируемое) решение уравнений Навье-Стокса, определённое для всех t > 0, либо привести контрпример, демонстрирующий возникновение особенности (разрыва) за конечное время. Для двумерного случая существование и гладкость решений доказаны. Для трёхмерного случая доказаны лишь локальные по времени результаты (существование на малом интервале времени) или при малых начальных данных.

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

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

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

Для решения нелинейной системы уравнений на каждом шаге по времени используются итерационные методы (например, метод SIMPLE, PISO). Особую сложность представляет моделирование турбулентности, для чего применяются либо прямое численное моделирование (DNS — требует огромных вычислительных ресурсов), либо модели турбулентности (RANS, LES, DES), которые вводят дополнительные уравнения для осреднённых характеристик.

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

Уравнения Навье-Стокса являются частным случаем общих уравнений движения сплошной среды — уравнений Коши.
В 2000 году задача о существовании и гладкости решений уравнений Навье-Стокса была включена в список семи «Проблем тысячелетия».
Аналитическое решение уравнений Навье-Стокса для течения в трубе круглого сечения (течение Пуазёйля) было получено ещё в 1840-х годах.
Для моделирования турбулентных течений с помощью DNS требуются суперкомпьютеры, способные выполнять вычисления с числом узлов сетки порядка Re^(9/4).

¶Источники

Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. Том 6. Гидродинамика. — М.: Наука, 1986.
Лойцянский Л. Г. Механика жидкости и газа. — М.: Дрофа, 2003.
Флетчер К. Вычислительные методы в динамике жидкости. — М.: Мир, 1991.
Проблемы тысячелетия: Уравнения Навье-Стокса. — Математический институт Клэя, 2000.