Турбулентність простими словами — від Рейнольдса до Колмогорова

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

1. Що таке турбулентність?

Турбулентність — це не випадковий шум. Це детермінований, хаотичний стан руху рідини, що характеризується вихорами на багатьох масштабах, які нелінійно взаємодіють. Її визначають три властивості:

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

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

2. Число Рейнольдса і перехід

Осборн Рейнольдс (1883) уприскував барвник у потік у трубі й виявив різкий перехід. Нижче критичного числа Рейнольдса барвник тік гладкими паралельними смужками (ламінарно). Вище нього барвник миттєво перемішувався по всій трубі (турбулентно):

Re = ρ \cdot U \cdot L / μ = U \cdot L / ν ρ = густина, U = швидкість, L = діаметр труби μ = динамічна в'язкість, ν = кінематична в'язкість = μ/ρ Потік у трубі: Re < 2300: ламінарний, стабільний 2300-4000: перехідний, переривчаста турбулентність («пориви») Re > 4000: повністю турбулентний

Фізична інтерпретація: Re — це відношення сил інерції до сил в'язкості. За низького Re в'язкість гасить будь-яке збурення, перш ніж воно встигне зрости. За високого Re інерція долає в'язке гасіння, і малі збурення підсилюються до турбулентних флуктуацій через нестійкість Кельвіна-Гельмгольца та інші види нестійкості.

Перехід не миттєвий — він відбувається через стадії: хвилі Толлміна-Шліхтінга → тривимірні нестійкості → турбулентні плями → повністю турбулентний потік. За сприятливого градієнта тиску (потік, що прискорюється над крилом літака поблизу передньої кромки) перехід затримується; несприятливі градієнти його прискорюють.

3. Енергетичний каскад Колмогорова

Андрій Колмогоров (1941) сформулював сучасне розуміння структури турбулентності. Енергія входить у турбулентний потік на інтегральному масштабі L (найбільші вихори, порядку розміру границі/геометрії). Ці великі вихори нестійкі, розпадаються на менші вихори, що розпадаються на ще менші, безперервно передаючи енергію вниз за масштабом — це й є каскад.

Вірш Річардсона 1922 (у переказі): «У великих вихорів є малі вихори, що живляться їхньою швидкістю; а в малих вихорів — ще менші, і так далі аж до в'язкості». Енергетичний спектр в інерційному піддіапазоні (Колмогоров 1941): E(k) \propto ε^(2/3) \cdot k^(-5/3) E(k) = енергія на хвильовому числі k (обернений масштаб довжини) ε = швидкість дисипації турбулентної кінетичної енергії (m²/s³) k = хвильове число = 2π/λ Цей закон степеня -5/3 експериментально підтверджено на 8 порядках величини за хвильовим числом.

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

4. Мікромасштаби Колмогорова

З двох параметрів, що керують малими масштабами — кінематичної в'язкості ν та швидкості дисипації ε — Колмогоров вивів найменші турбулентні масштаби довжини, часу та швидкості:

Масштаб довжини Колмогорова: η = (ν³/ε)^(1/4) Масштаб часу Колмогорова: τ_η = (ν/ε)^(1/2) Масштаб швидкості Колмогорова: u_η = (νε)^(1/4) Зв'язок з інтегральним масштабом L: η/L \approx Re^(-3/4) \to Re = 10⁴: η/L \approx 1/1000 \to 1000 масштабів довжини між великими та малими вихорами \to Re = 10⁸: η/L \approx 1/10^6 \to потрібні DNS-сітки з 10^18 комірок — неможливо!

Саме тому пряме чисельне моделювання (DNS) — розв'язання всіх масштабів — обмежене Re ≲ 10⁴. Інженерні потоки (крила літаків за Re ≈ 10⁷, корпуси кораблів за Re ≈ 10⁹) потребують моделей турбулентності, щоб уникнути явного розв'язання астрономічної кількості потрібних комірок.

5. Доріжка вихорів Кармана

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

Число Струхаля: St = f \cdot D / U \approx 0.198 \cdot (1 - 19.7/Re) f = частота сходження вихорів (Hz) D = діаметр циліндра U = швидкість потоку Перетворено для частоти: f = St \cdot U / D Приклад: циліндр діаметром 0.01 m у вітрі 5 m/s: f \approx 0.2 \times 5 / 0.01 = 100 Hz (чутний вихровий тон!)

Доріжки вихорів Кармана спричиняють резонанс конструкцій. Міст Такома-Нарровз (1940) обвалився, коли сходження вихорів збіглося з власною частотою конструкції — це явище тепер називають флатером. Сучасні мости використовують аеродинамічно сформовані перерізи та обтічники, щоб розладнати зв'язок частот вихор-тіло.

За низького Re (40–190) слід ламінарний, але осцилюючий (2D-доріжка Кармана). За вищого Re (>190) вихори розпадаються на 3D турбулентні структури у сліді.

6. Інженерні моделі турбулентності

Моделі турбулентності RANS замінюють поле пульсуючої швидкості турбулентною в'язкістю ν_t, що доповнює молекулярну в'язкість:

Гіпотеза Буссінеска: турбулентне напруження τ_ij = ρν_t (\partialuᵢ/\partialxⱼ + \partialuⱼ/\partialxᵢ) Завдання: як обчислити ν_t?

Модель k-ε (Лаундер і Сполдінг, 1974)

Два додаткові рівняння переносу: k = турбулентна кінетична енергія = ½⟨u'²⟩ (m²/s²) ε = швидкість дисипації (m²/s³) ν_t = C_μ \cdot k² / ε (C_μ = 0.09, емпірична) Рівняння переносу: Dk/Dt = P_k - ε + \nabla\cdot[(ν + ν_t/σ_k)\nablak] Dε/Dt = C₁ε(ε/k)P_k - C₂ε(ε²/k) + \nabla\cdot[(ν + ν_t/σ_ε)\nablaε] P_k = породження = 2ν_t S_ij : S_ij (S = тензор швидкості деформації) Стандартні константи: C₁ε=1.44, C₂ε=1.92, σ_k=1.0, σ_ε=1.3

Сильні сторони k-ε: надійна, добре перевірена для приєднаних потоків і потоків у трубах. Модель за замовчуванням у багатьох промислових пакетах CFD (Fluent, OpenFOAM). Слабкість: неточна для відривних потоків (рециркуляція за тілами), несприятливих градієнтів тиску та потоків зі значною кривиною ліній течії.

Модель k-ω SST (Ментер, 1994)

Використовує k-ω (краща поведінка біля стінки) у пограничному шарі та переходить до k-ε у вільному потоці. Модифікація перенесення дотичних напружень (SST) також обмежує ν_t у несприятливих градієнтах тиску, поліпшуючи прогнози відривних потоків. Найпоширеніша робоча модель в аерокосмічному CFD сьогодні.

Спаларт-Аллмарас (1992)

Одне рівняння переносу для ν̃ (модифікованої турбулентної в'язкості). Дуже швидка (одне додаткове рівняння проти двох для k-ω). Розроблена спеціально для аеродинамічних потоків; менш точна для внутрішніх потоків і рециркуляції.

7. DNS, LES і RANS — практичний погляд

Метод	Розв'язує	Комірок сітки (3D, Re=10⁶)	Час	Хто це використовує
DNS	Усі масштаби	~10¹⁴ (неможливо)	Роки процесорного часу	Фундаментальні дослідження турбулентності, Re < 10⁴
LES	Великі вихори; малі моделюються	~10⁸–10¹⁰	Дні–тижні	Аероакустика, горіння, погода, складні потоки
RANS (k-ε)	Лише осереднений потік	~10⁵–10⁷	Години	Промисловість (авто, літаки, ОВК, турбіни)
LES з моделюванням стінки	Великі вихори; стінка моделюється	~10⁷–10⁹	Дні	LES за високих Re в аеродинаміці

8. Нерозв'язана проблема

Попри понад століття досліджень, турбулентність лишається однією з Задач тисячоліття Інституту Клея (приз 1 мільйон доларів): довести або спростувати, що гладкі розв'язки 3D рівнянь Нав'є-Стокса завжди існують (без сингулярності за скінченний час).

Практично ми можемо точно моделювати турбулентність за допомогою DNS за низьких Re та модельних підходів за вищих Re. Але ми не можемо передбачити турбулентність із перших принципів без запуску моделювання. Рівняння Нав'є-Стокса нестійкі за Ляпуновим за високих Re — крихітні початкові збурення зростають експоненційно (хаос). Саме тому погода передбачувана лише на ~10 днів, і саме тому слід кожного літака унікальний.

Машинне навчання для турбулентності: нещодавні роботи (2020-ті) використовують нейронні мережі, щоб навчати замикання турбулентності безпосередньо з даних DNS, замінюючи вручну налаштовані константи k-ε навчуваними моделями. Ці «моделі турбулентності на основі даних» краще узагальнюються на нерозрахункові умови, ніж класичний RANS, але потребують навчальних даних DNS — які самі вимагають умов низького Re. Фундаментальна проблема «курки і яйця» лишається.