Рівняння Нав'є-Стокса — математика течії рідини

1. Гіпотеза суцільного середовища

Реальні рідини складаються з дискретних молекул — приблизно 2,7 × 10²⁵ молекул на кубічний метр повітря на рівні моря. Відстежувати кожну молекулу окремо обчислювально неможливо для макроскопічних течій, та й непотрібно: у масштабах, набагато більших за довжину вільного пробігу (~70 нм у повітрі за нормальних умов), статистичні середні збігаються, і рідина поводиться як суцільне середовище.

Гіпотеза суцільного середовища розглядає властивості рідини — густину ρ, швидкість u, тиск P, температуру T — як гладкі поля, визначені в кожній точці простору. Частинка рідини — це нескінченно малий об'єм, який усе ж достатньо великий, щоб містити мільйони молекул, тож термодинамічні величини добре визначені.

Ця гіпотеза порушується за дуже низьких тисків (динаміка розрідженого газу, вхід космічного апарата в атмосферу на великій висоті) або всередині наномасштабних пристроїв. У таких режимах число Кнудсена Kn = λ/L перевищує близько 0,1, де λ — довжина вільного пробігу, а L — характерний масштаб довжини, і доводиться вдаватися до кінетичної теорії або молекулярної динаміки.

2. Нестисливі рівняння Нав'є-Стокса

Для нестисливої ньютонівської рідини (стала густина, напруження пропорційне швидкості деформації) рівняння руху є формулюванням другого закону Ньютона, застосованого до частинки рідини:

Розберемо кожен доданок у рівнянні імпульсу:

• ρ ∂u/∂t — локальне прискорення (як швидкість змінюється в фіксованій точці з часом)
• ρ (u·∇u) — адвективне (конвективне) прискорення; частинка рідини, що рухається в область іншої швидкості, зазнає цього нелінійного члена. Це джерело практично всієї складності.
• −∇P — сила градієнта тиску; рідина прискорюється від високого тиску до низького
• μ∇²u — в'язка дифузія; внутрішнє тертя згладжує градієнти швидкості. μ — динамічна в'язкість (Па·с).
• f — об'ємні сили на одиницю об'єму (гравітація, електромагнітні сили тощо)

Рівняння нерозривності ∇·u = 0 стверджує, що дивергенція швидкості дорівнює нулю — рідина не створюється й не зникає в жодному елементі об'єму. Разом ці чотири скалярні рівняння (три компоненти імпульсу + нерозривність) визначають чотири невідомі u = (u, v, w) та P.

Матеріальна похідна

Комбінація ∂/∂t + u·∇ називається матеріальною (або субстанціальною) похідною D/Dt. Вона описує швидкість зміни, що відбувається разом із частинкою рідини під час її руху:

Отже, рівняння імпульсу NS — це просто: ρ Du/Dt = −∇P + μ∇²u + f — маса, помножена на прискорення, дорівнює сумарній силі, точнісінько як задумав Ньютон. Складність полягає у нелінійному адвективному члені (u·∇)u, який пов'язує компоненти швидкості й породжує хаотичну поведінку на високих швидкостях.

3. Число Рейнольдса: Re = ρUL/μ

У 1883 році Осборн Рейнольдс впорснув тонкий струмінь барвника у трубну течію й виявив різкий перехід: нижче критичної швидкості барвник тік прямою лінією, вище — барвник вибухав хаотичним перемішуванням. Він визначив ключовий безрозмірний параметр, що керує цим переходом:

де ρ — густина рідини (кг/м³), U — характерна швидкість (м/с), L — характерний масштаб довжини (м), μ — динамічна в'язкість (Па·с), а ν = μ/ρ — кінематична в'язкість (м²/с).

Число Рейнольдса вимірює відношення інерційних сил (які прагнуть підтримувати рух і створювати вихори) до в'язких сил (які чинять опір відносному руху й гасять збурення):

• Re < ~2300 — ламінарна трубна течія; домінує в'язкість, збурення гасяться
• Re ~ 2300–4000 — перехідний режим; течія чергується між ламінарною та турбулентною
• Re > ~4000 — повністю турбулентна трубна течія

Ці пороги залежать від геометрії. Для потоку навколо сфери «криза опору» настає приблизно при Re ~ 3 × 10⁵, коли пограничний шар переходить у турбулентний і коефіцієнт опору різко падає — саме тому на м'ячах для гольфу є ямки (вони змушують пограничний шар перейти раніше, знижуючи ефективне Re і зменшуючи опір на швидкостях, властивих гольфу).

4. Ламінарна проти турбулентної течії

Ламінарна течія характеризується гладкими, паралельними лініями течії без бокового перемішування. Кожен шар рідини ковзає повз сусідні шари з чітко визначеними градієнтами швидкості. Профіль швидкості в круглій трубі (течія Хагена-Пуазейля) — парабола: u(r) = (ΔP/4μL)(R² − r²), де R — радіус труби, а r — радіальна відстань від центра.

Турбулентна течія має зовсім інший характер. Вона проявляє:

• Хаотичні, тривимірні коливання швидкості на багатьох масштабах довжини й часу
• Посилене перемішування — турбулентність переносить імпульс, тепло та масу значно ефективніше, ніж ламінарна дифузія
• Значну дисипацію енергії; турбулентна трубна течія вимагає набагато вищих градієнтів тиску, ніж ламінарна течія при тій самій витраті
• Чутливість до початкових умов — ефект метелика означає, що крихітні збурення зростають експоненційно

Перехід від ламінарної до турбулентної течії — не простий перемикач «увімкнено/вимкнено». Він проходить через послідовність нестійкостей: малі збурення (хвилі Толлміна-Шліхтінга) зростають, стають тривимірними, формують шпилькоподібні вихори і зрештою розпадаються на повний каскад турбулентних вихорів. Рівняння Орра-Зоммерфельда керує лінійною стійкістю паралельних зсувних течій і передбачає, які довжини хвиль збурень посиляться.

5. Каскад Колмогорова

У 1941 році Андрій Колмогоров опублікував, мабуть, найуспішнішу теорію в турбулентності: статистичний опис того, як енергія переходить через масштаби довжини. Турбулентна течія не є випадковою в неструктурований спосіб — вона має точну, універсальну статистичну архітектуру.

Картина така: каскад — великомасштабні вихори (внесені середньою течією на інтегральному масштабі L) розпадаються на дедалі менші вихори, передаючи кінетичну енергію вниз за масштабом без дисипації, доки не досягнуть мікромасштабу Колмогорова η, де в'язка дисипація нарешті перетворює кінетичну енергію на тепло:

де ε — швидкість дисипації енергії на одиницю маси (м²/с³), а ν — кінематична в'язкість. Для атмосферної турбулентності η ~ 1 мм; для океану η ~ 1 см.

У інерційному піддіапазоні — масштабах між L та η — розмірний аналіз Колмогорова передбачає відомий степеневий закон −5/3 для спектра енергії:

де k — хвильове число (1/довжина). Це означає, що вихори розміру r несуть кінетичну енергію ∝ r^(2/3). Спектр −5/3 підтверджений в аеродинамічних трубах, океанських вимірюваннях, атмосферних зондуваннях та даних сонячного вітру — вражаюча універсальність.

Відношення найбільших до найменших масштабів у турбулентності масштабується як Re^(3/4), тобто повне розв'язання турбулентної течії у прямому чисельному моделюванні (DNS) вимагає сітки з Re^(9/4) точок — обчислювально неможливо для інженерних чисел Рейнольдса. DNS реактивного двигуна при Re ~ 10⁷ вимагав би ~10¹⁶ вузлів сітки.

6. Пограничні шари та відрив потоку

За великих чисел Рейнольдса вплив в'язкості обмежується тонкою областю біля твердих поверхонь, яку називають пограничним шаром, уперше описаним Людвігом Прандтлем 1904 року. Це прозріння перетворило аеродинаміку: поза пограничним шаром течію можна вважати нев'язкою (рівняння Ейлера), що радикально спрощує аналіз.

Для ламінарного пограничного шару над плоскою пластиною (розв'язок Блазіуса) товщина зростає як:

Шар починається тонким на передній кромці й потовщується вниз за потоком. Всередині нього швидкість зростає від нуля на стінці (умова прилипання) до значення вільного потоку U на відстані δ. Напруження зсуву на стінці τ_w = μ(∂u/∂y)|_{y=0} створює опір тертя шкіри.

Несприятливі градієнти тиску та відрив

Коли тиск зростає за напрямком течії (несприятливий градієнт тиску, ∂P/∂x > 0), рідина біля стінки — вже сповільнена в'язкістю — може зупинитися, а потім піти назад. Пограничний шар відривається від поверхні, створюючи велику зону рециркуляції та різке зростання опору тиску (опору форми).

Відрив відповідає за звалювання крил літака на великих кутах атаки, слід за тупими тілами та шум, що виникає при обтіканні порожнин. Керування відривом пограничного шару — за допомогою шорсткості поверхні, відсмоктування, видування або вихорогенераторів — одна з центральних задач аеродинамічного проєктування.

7. Проблема тисячоліття

У 2000 році Математичний інститут Клея перерахував сім проблем тисячоліття, кожна з яких несла нагороду в 1 000 000 доларів США. Проблема існування та гладкості Нав'є-Стокса — одна з них — і станом на 2026 рік вона лишається невирішеною.

Проблема, сформульована неформально: за гладких початкових даних (добре поведеного початкового поля швидкості у 3D) чи існують гладкі розв'язки нестисливих рівнянь NS для всіх майбутніх часів? Чи можуть розв'язки розвинути сингулярності — точки, де швидкість або тиск стають нескінченними за скінченний час?

У двох вимірах глобальна регулярність доведена: 2D-розв'язки NS залишаються гладкими назавжди. У трьох вимірах задача відкрита. Нелінійний член (u·∇)u в принципі може концентрувати енергію в дедалі тонших структурах швидше, ніж в'язкість здатна її розсіяти — потенційно спричиняючи «вибух за скінченний час».

Що відомо:

• Локальне існування: гладкі розв'язки існують протягом короткого часу після будь-якої гладкої початкової умови (Лерей, 1934)
• Слабкі розв'язки Лерея-Гопфа існують глобально, але можуть бути неєдиними або негладкими
• Якщо вибух стається, завихреність ‖ω‖_{L∞} має розходитися (критерій Біла-Като-Мажди)
• Умовна регулярність: розв'язки гладкі, якщо швидкість задовольняє певні умови інтегровності (критерії Ладиженської-Проді-Серріна)

8. Числові методи: CFD на практиці

Оскільки аналітичні розв'язки існують лише для найпростіших геометрій, практична гідродинаміка спирається на обчислювальну гідродинаміку (CFD). Основні числові підходи:

• Метод скінченних об'ємів (FVM) — область поділяється на контрольні об'єми; закони збереження застосовуються до кожної комірки. Більшість промислових CFD-кодів (OpenFOAM, ANSYS Fluent) використовують FVM для прямого забезпечення збереження.
• Метод скінченних різниць (FDM) — похідні апроксимуються різницями між значеннями у вузлах сітки. Простий у реалізації; використовується в прогнозуванні погоди й моделюванні океану.
• Спектральні методи — швидкість розкладається в ряд Фур'є або поліноми Чебишова. Дуже точні для гладких течій у простих геометріях; метод вибору для DNS турбулентності.
• Метод ґраткового Больцмана (LBM) — замість прямого розв'язання NS, LBM моделює спрощене кінетичне рівняння на ґратці. NS виникає в макроскопічній межі. Чудово підходить для складних геометрій і багатофазних течій.

Моделювання турбулентності додає ще один рівень складності. DNS розв'язує всі масштаби, але здійсненне лише за низьких Re. Метод великих вихорів (LES) розв'язує великі вихори й моделює малі. Осереднені за Рейнольдсом рівняння Нав'є-Стокса (RANS) моделюють усю турбулентність із рівняннями замикання (k-ε, k-ω, SST) і є робочим конем промислової CFD, попри свої фундаментальні наближення.

Сучасні підходи машинного навчання — фізично-інформовані нейронні мережі (PINN) та нейронні оператори — починають прискорювати CFD, навчаючись відображенням розв'язків безпосередньо з навчальних даних, хоча повна заміна першопринципних розв'язувачів залишається віддаленою метою.