📊 Спотлайт #45 Вер 2030 • ~18 хв читання

Гідродинаміка та турбулентність — Нав'є-Стокс, число Рейнольдса та каскад енергії

Від ламінарного потоку в трубі до самоподібного хаосу повністю розвиненої турбулентності гідродинаміка пов'язує повсякденні явища з деякими з найглибших нерозв'язаних проблем математики та фізики. Цей спотлайт простежує рівняння Нав'є-Стокса, примежові шари, вихрову динаміку, аеродинаміку та каскад Колмогорова.

Рівняння Нав'є-Стокса

Уся класична гідродинаміка спирається на два закони збереження, застосовані до суцільного середовища рідини: збереження маси (рівняння неперервності) та збереження імпульсу (другий закон Ньютона для елемента рідини). Для нестисливої ньютонівської рідини це дає рівняння Нав'є-Стокса:

∇ · u = 0     (нестисливість)

ρ(∂u/∂t + (u · ∇)u) = −∇p + μ∇²u + ρg

локальне + конвективне = градієнт тиску + в'язке + об'ємна сила

Член ( u · ∇) u — це нелінійне конвективне прискорення — джерело турбулентності. Він представляє розтягування та нахил завихреності рідини, коли швидші елементи рідини випереджають повільніші. В'язкий член μ∇²u (динамічна в'язкість μ) згладжує градієнти швидкості, розсіюючи кінетичну енергію у вигляді тепла. Доведення того, що гладкі розв'язки 3D-рівнянь Нав'є-Стокса завжди існують, — одна з семи Проблем тисячоліття (Математичний інститут Клея, приз $1 мільйон).

Фізична інтерпретація

Кожен член має пряме фізичне значення. Градієнт тиску −∇p жене потік від високого до низького тиску (течія Пуазейля в трубі, метеорологічні вітри). В'язкість протидіє деформації: вода має μ ≈ 10³ × 10−³ Па·с при 20 °C; повітря має μ ≈ 1,8 × 10−&sup5; Па·с. Густина ρ множить інерцію: важчим рідинам потрібно більше сили для прискорення.

Число Рейнольдса та ламінарно-турбулентний перехід

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

Re = ρ U L / μ = U L / ν

U = характерна швидкість   L = характерна довжина
ν = μ/ρ = кінематична в'язкість (м²/с)

Для течії в трубі (Хаген-Пуазейль): ламінарна при Re < 2300, перехід 2300–4000, повністю турбулентна Re > 4000. Для зовнішньої течії над плоскою пластиною: перехід близько Re ≈ 5 × 10&sup5;. Комерційний літак летить на Re ∼ 10&sup7;; бактерія плаває на Re ∼ 10−&sup4; (домінує в'язкість; інерція незначна). При низьких Re потік зворотний — знаменитий експеримент Дж. І. Тейлора з розворотом чорнила 1966 року драматично це продемонстрував.

Механізми переходу

Турбулентність не виникає спонтанно. Перехід ініціюється нестійкостями: хвилями Толлміна-Шліхтінга у примежових шарах (лінійна нестійкість, Reδ ≈ 520), вторинними нестійкостями та обхідним переходом при реальній турбулентності набігаючого потоку. Рівняння Орра-Зоммерфельда керує лінійною стійкістю паралельних зсувних течій і передбачає критичні числа Рейнольдса для течій у каналах і трубах.

🌊

Хвилі у ванні та Симулятор Бернуллі дозволяють інтерактивно регулювати швидкість потоку та спостерігати, як змінюються тиск, швидкість та режим течії залежно від числа Рейнольдса.

Рівняння Бернуллі та застосування

Для стаціонарної, нев'язкої, нестисливої течії вздовж лінії струму рівняння Бернуллі є твердженням про збереження енергії:

p + ½ρu² + ρgz = стала вздовж лінії струму

статичний тиск + динамічний тиск + гідростатичний тиск = const

Застосування повсюдні: витратоміри Вентурі, трубки Піто на літаках, карбюратори, тяга димаря, ефект Коанда, і найважливіше — аеродинамічна підйомна сила. Крило (аерофойл) має форму, яка прискорює повітря над верхньою поверхнею (коротший шлях крізь западину тиску) більше, ніж над нижньою поверхнею, створюючи різницю тиску Δp, що створює підйомну силу L = Δp × площа. Бернуллі суворо застосовується лише вздовж лінії струму; повна картина включає циркуляцію Γ (теорема Кутта-Жуковського: L = ρ U Γ на одиницю розмаху).

Пастка: помилка про однаковий час прольоту

Поширена помилкова думка стверджує, що повітря розділяється над крилом і має «знову зустрітися» на задній кромці (маючи на увазі, що верхнє повітря рухається швидше, оскільки шлях довший). Це неправда: частинки повітря не зустрічаються. Реальний механізм — умова Кутта (гладкий відхід із задньої кромки) у поєднанні з циркуляцією, забезпеченою умовою прилипання та в'язкістю.

✈️

Симулятор підйомної сили аерофойлу показує лінії струму, розподіл тиску та відношення L/D для профілів NACA за різних кутів атаки. Досліджуйте умову зриву, коли відрив потоку спричиняє колапс підйомної сили.

Примежові шари

Умова прилипання вимагає, щоб швидкість рідини на твердій стінці дорівнювала швидкості стінки. Це створює тонку зону — примежовий шар — де швидкість переходить від нуля (на стінці) до значення набігаючого потоку U∞. Розв'язок Блазіуса для ламінарного примежового шару над плоскою пластиною дає товщину примежового шару:

δ(x) ≈ 5x / √Rex    де Rex = U∞x/ν

Дотичне напруження на стінці: τw = 0,332 ρ U∞² / √Rex
Коефіцієнт тертя шкіри: Cf = 0,664 / √Rex

У турбулентному примежовому шарі профіль швидкості слідує логарифмічному закону стінки: u+ = (1/κ) ln y+ + B, зі сталою фон Кармана κ ≈ 0,41 та B ≈ 5,0. В'язкий підшар (y+ < 5) домінується молекулярною в'язкістю; логарифмічний шар (30 < y+ < 300) та зовнішня область домінують у загальному опорі. Сучасні літаки витрачають величезні інженерні зусилля на зменшення внеску турбулентного примежового шару в опір тертя шкіри, який становить приблизно 50% загального опору.

Вихрова динаміка та вихрова доріжка фон Кармана

Завихреність ω = ∇ × u вимірює локальне обертання елементів рідини. Теореми Гельмгольца (1858) показують, що вихрова трубка в нев'язкій рідині переноситься потоком, зберігає свою силу і не може закінчитися всередині рідини. Теорема циркуляції Кельвіна: DΓ/Dt = 0 у нев'язкій баротропній рідині.

Обтікаючи круговий циліндр при помірному Re (47–200), слід стає періодичним: вихори сходять почергово з верхньої та нижньої сторін, утворюючи вихрову доріжку фон Кармана. Число Струхаля St = f d / U ≈ 0,2 керує частотою сходження f для діаметра циліндра d. Ця вихрозбуджена вібрація відома тим, що спричинила обвал мосту Такома-Нарроус у 1940 році (хоча повний механізм був аеропружним флаттером, а не чистим вихровим резонансом). Морські нафтові стояки та лінії електропередач також вразливі і потребують гвинтових ребер для порушення періодичного сходження.

Струхаль: St = f d / U ≈ 0,198 (1 − 19,7/Re)   для Re = 250–2×10&sup5;
Частота сходження: f = St · U / d

Турбулентність та каскад енергії Колмогорова

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

Масштаб Колмогорова: η = (ν³/ε)¹⁄&sup4;
Часовий масштаб: τη = (ν/ε)¹⁄²
Масштаб швидкості: uη = (νε)¹⁄&sup4;

Енергетичний спектр (інерційний підінтервал): E(k) = C ε²⁄³ k−&sup5;⁄³
(закон Колмогорова −5/3, підтверджений експериментально)

Тут ε — швидкість розсіювання турбулентної кінетичної енергії на одиницю маси (Вт/кг), k — хвильове число. Степеневий закон −5/3 підтверджений для атмосферної турбулентності (масштаби 10−1000 м), океанських течій, потоку в трубах і навіть на інших планетах. Відношення L/η масштабується як Re3/4: при Re ∼ 10&sup6; діапазон охоплює чотири порядки масштабу. Саме тому пряме чисельне моделювання (DNS) турбулентності при високому Re вимагає ∼Re9/4 точок сітки і є обчислювально нездійсненним вище Re ∼ 10&sup4; навіть на сучасних суперкомп'ютерах.

Моделі турбулентності

Оскільки DNS занадто дороге, інженерія використовує моделі турбулентності. Осереднені за Рейнольдсом рівняння Нав'є-Стокса (RANS): усереднення рівнянь у часі дає тензор напружень Рейнольдса −ρ<u′u′>, що замикається вихровою в'язкістю (моделі k-ω, k-ε). Моделювання великих вихорів (LES): розв'язуються вихори, більші за масштаб фільтра, моделюються лише підсіткові напруження (модель Смагоринського). Гібридний RANS/LES (DES): дешевий поблизу стінок, точний у відірваних областях. Усі моделі вводять вільні параметри, що потребують калібрування за експериментами.

Стисливий потік та ударні хвилі

Коли швидкість потоку наближається до локальної швидкості звуку, припущення нестисливості порушується. Число Маха Ma = U/a (a = √(γRT) для ідеального газу) характеризує режим: Ma < 0,3 нестисливий; 0,3–0,8 дозвуковий стисливий; 0,8–1,2 трансзвуковий (з локальними надзвуковими зонами); 1,2–5 надзвуковий; >5 гіперзвуковий.

При Ma > 1 утворюється прямий стрибок: потік різко сповільнюється, густина і тиск стрибкоподібно змінюються, а ентропія зростає (косі стрибки в 2D/3D). Умови стрибка Ренкіна-Гюгоніо пов'язують властивості до і після стрибка:

ρ2/ρ1 = (γ+1)Ma1² / [(γ−1)Ma1² + 2]
p2/p1 = [2γMa1² − (γ−1)] / (γ+1)
Ma2² = [1 + (γ−1)/2 · Ma1²] / [γMa1² − (γ−1)/2]

Ключове застосування: надзвуковий літак має прискоритися крізь трансзвуковий звуковий бар'єр (коефіцієнт опору досягає піку при Ma ≈ 1 через хвильовий опір) і працювати на надзвуковій швидкості, де хвильовий опір керується стрілоподібними крилами та правилом площ. Особливість Прандтля-Глауерта (нескінченна підйомна сила при Ma = 1 у лінеаризованій теорії) — артефакт моделі — реальні потоки неперервні, але нелінійні.

Аналіз розмірностей та подібність

Теорема Бекінгема Π стверджує, що фізично значуще рівняння, яке пов'язує n розмірних величин з k незалежними розмірностями, може бути записане через n−k безрозмірних груп (Π-параметрів). Гідродинаміка постійно використовує це: числа Re, Ma, Fr (Фруда), We (Вебера), Eu (Ейлера) та Струхаля. Дві течії є динамічно подібними, якщо всі відповідні безрозмірні групи збігаються — основа для випробувань в аеродинамічних трубах (масштабні моделі за узгодженого Re) та випробувань у буксирувальних басейнах для кораблів (узгоджене Fr).

Ключові безрозмірні числа

Re = ρUL/μ (інерція/в'язкість) • Ma = U/a (інерція/стисливість) • Fr = U/√(gL) (інерція/гравітація) • We = ρU²L/σ (інерція/поверхневий натяг) • Pr = ν/α (імпульс/теплова дифузійність) • Nu = hL/k (конвективний/ кондуктивний теплообмін)

Надплинність та квантова турбулентність

При температурах нижче λ-точки (Tλ ≈ 2,17 К) рідкий 4He стає надплинним (конденсат Бозе-Ейнштейна). Його в'язкість тотожно зникає, і він тече без опору. Завихреність квантована: будь-яка циркуляція має бути цілим кратним κ = h/m (квант циркуляції). Ядра вихорів мають товщину ∼1 Å і несуть рівно один квант κ. Квантова турбулентність — це заплутування цих квантованих вихрових ниток; на великих масштабах вона демонструє спектр Колмогорова −5/3, натякаючи, що класичний каскад виникає з квантової вихрової динаміки.

Симуляції на цій платформі

Наша категорія гідродинаміки включає практичні інструменти по всьому спектру чисел Рейнольдса: