Dynamika płynów i turbulencja — Navier-Stokes, liczba Reynoldsa i kaskada energii

Równania Naviera-Stokesa

Cała klasyczna dynamika płynów opiera się na dwóch prawach zachowania stosowanych do płynu ciągłego: zachowaniu masy (równanie ciągłości) i zachowaniu pędu (druga zasada dynamiki Newtona dla elementu płynu). Dla nieściśliwego płynu newtonowskiego dają one równania Naviera-Stokesa:

Człon ( u · ∇) u to nieliniowe przyspieszenie konwekcyjne — źródło turbulencji. Reprezentuje on rozciąganie i przechylanie wirowości płynu, gdy szybsze elementy płynu wyprzedzają wolniejsze. Człon lepkościowy μ∇²u (lepkość dynamiczna μ) wygładza gradienty prędkości, rozpraszając energię kinetyczną jako ciepło. Udowodnienie, że gładkie rozwiązania trójwymiarowych równań Naviera-Stokesa zawsze istnieją, jest jednym z siedmiu problemów milenijnych (Clay Mathematics Institute, nagroda 1 mln USD).

Interpretacja fizyczna

Każdy człon ma bezpośrednie znaczenie fizyczne. Gradient ciśnienia −∇p napędza przepływ od wysokiego do niskiego ciśnienia (przepływ Poiseuille'a w rurze, wiatry meteorologiczne). Lepkość opiera się deformacji: woda ma μ ≈ 10³ × 10−³ Pa·s przy 20°C; powietrze ma μ ≈ 1,8 × 10−&sup5; Pa·s. Gęstość ρ mnoży bezwładność: cięższe płyny wymagają więcej siły do przyspieszenia.

Liczba Reynoldsa i przejście laminarno-turbulentne

To, czy przepływ jest gładki (laminarny), czy chaotyczny (turbulentny), zależy od rywalizacji między siłami bezwładności a lepkości, kwantyfikowanej przez bezwymiarową liczbę Reynoldsa:

Dla przepływu w rurze (Hagen-Poiseuille): laminarny dla Re < 2300, przejściowy 2300–4000, w pełni turbulentny Re > 4000. Dla przepływu zewnętrznego nad płaską płytą: przejście przy Re ≈ 5 × 10&sup5;. Samolot pasażerski lata z prędkością odpowiadającą Re ∼ 10&sup7;; bakteria pływa przy Re ∼ 10−&sup4; (zdominowana przez lepkość; bezwładność jest pomijalna). Przy niskim Re przepływ jest odwracalny — słynny eksperyment G.I. Taylora z odwracaniem atramentu z 1966 roku dramatycznie to zademonstrował.

Mechanizmy przejścia

Turbulencja nie pojawia się spontanicznie. Przejście jest wywoływane przez niestabilności: fale Tollmiena-Schlichtinga w warstwach przyściennych (niestabilność liniowa, Re_δ ≈ 520), niestabilności wtórne oraz przejście typu bypass pod wpływem rzeczywistej turbulencji swobodnego strumienia. Równanie Orra-Sommerfelda rządzi liniową stabilnością równoległych przepływów ścinających i przewiduje krytyczne liczby Reynoldsa dla przepływów w kanałach i rurach.

Równanie Bernoulliego i jego zastosowania

Dla ustalonego, nielepkiego, nieściśliwego przepływu wzdłuż linii prądu równanie Bernoulliego jest stwierdzeniem zachowania energii:

Zastosowania są wszechobecne: przepływomierze Venturiego, rurki Pitota w samolotach, gaźniki, ciąg kominowy, efekt Coandy, a przede wszystkim siła nośna aerodynamiczna. Skrzydło (profil aerodynamiczny) jest ukształtowane tak, by przyspieszać powietrze nad górną powierzchnią (krótsza droga przez dolinę ciśnienia) bardziej niż pod dolną, tworząc różnicę ciśnień Δp, która wytwarza siłę nośną L = Δp × powierzchnia. Bernoulli stosuje się ściśle tylko wzdłuż linii prądu; pełna historia obejmuje cyrkulację Γ (twierdzenie Kutty-Żukowskiego: L = ρ U Γ na jednostkę rozpiętości).

Pułapka: fałszywy mit równego czasu przelotu

Powszechne błędne przekonanie głosi, że powietrze rozdziela się nad skrzydłem i musi „ponownie się spotkać" na krawędzi spływu (co sugeruje, że górne powietrze porusza się szybciej, ponieważ droga jest dłuższa). To fałsz: elementy powietrza nie spotykają się ponownie. Rzeczywistym mechanizmem jest warunek Kutty (gładkie opuszczenie krawędzi spływu) w połączeniu z cyrkulacją wymuszoną przez warunek braku poślizgu i lepkość.

Warstwy przyścienne

Warunek braku poślizgu wymaga, aby prędkość płynu przy ściance stałej była równa prędkości ścianki. Tworzy to cienką strefę — warstwę przyścienną — w której prędkość przechodzi od zera (przy ściance) do wartości swobodnego strumienia U∞. Rozwiązanie Blasiusa dla laminarnej warstwy przyściennej nad płaską płytą daje grubość warstwy przyściennej:

Wewnątrz turbulentnej warstwy przyściennej profil prędkości podlega prawu logarytmicznemu ściany: u+ = (1/κ) ln y+ + B, ze stałą von Kármána κ ≈ 0,41 i B ≈ 5,0. Podwarstwa lepka (y+ < 5) jest zdominowana przez lepkość molekularną; warstwa logarytmiczna (30 < y+ < 300) i obszar zewnętrzny dominują nad całkowitym oporem. Nowoczesne samoloty poświęcają ogromny wysiłek inżynierski na redukcję wkładu turbulentnej warstwy przyściennej w opór tarcia powierzchniowego, który odpowiada za mniej więcej 50% całkowitego oporu.

Dynamika wirów i ulica wirowa von Kármána

Wirowość ω = ∇ × u mierzy lokalną rotację elementów płynu. Twierdzenia Helmholtza (1858) pokazują, że rurka wirowa w nielepkim płynie jest konwekowana wraz z przepływem, zachowuje swoją siłę i nie może się kończyć w płynie. Twierdzenie o cyrkulacji Kelvina: DΓ/Dt = 0 w nielepkim płynie barotropowym.

Za walcem kołowym przy umiarkowanym Re (47–200) ślad staje się okresowy: wiry są zrzucane naprzemiennie z górnej i dolnej strony, tworząc ulicę wirową von Kármána. Liczba Strouhala St = f d / U ≈ 0,2 rządzi częstotliwością zrzucania f dla średnicy walca d. To wywoływane wirami drganie słynnie spowodowało zawalenie się mostu Tacoma Narrows w 1940 roku (choć pełny mechanizm to trzepotanie aeroelastyczne, nie czysty rezonans wirowy). Morskie ryzery ropociągów i linie energetyczne są również podatne i wymagają śrubowych żeber zakłócających okresowe zrzucanie.

Turbulencja i kaskada energii Kołmogorowa

W pełni turbulentny przepływ zawiera wiry obejmujące ogromny zakres skal. Andrey Kołmogorow (1941) zaproponował, że turbulencja jest samopodobna w podzakresie inercyjnym: energia dostarczana w dużych skalach (L, skala całkowa) kaskaduje w dół poprzez rozszczepianie wirów bez dyssypacji, aż do mikroskali Kołmogorowa η, gdzie lepkość w końcu zamienia energię kinetyczną w ciepło.

Tutaj ε to tempo dyssypacji turbulentnej energii kinetycznej na jednostkę masy (W/kg), k to liczba falowa. Prawo potęgowe −5/3 zostało zweryfikowane w turbulencji atmosferycznej (skale 10–1000 m), prądach oceanicznych, przepływie w rurach, a nawet na innych planetach. Stosunek L/η skaluje się jak Re^3/4: przy Re ∼ 10&sup6; zakres obejmuje cztery dekady skali. Dlatego bezpośrednia symulacja numeryczna (DNS) turbulencji przy wysokim Re wymaga ∼Re^9/4 punktów siatki i jest obliczeniowo niewykonalna powyżej Re ∼ 10&sup4; nawet na dzisiejszych superkomputerach.

Modele turbulencji

Ponieważ DNS jest zbyt kosztowna, inżynieria stosuje modele turbulencji. Uśrednione po Reynoldsie równania Naviera-Stokesa (RANS): uśredniają w czasie równania, dając tensor naprężeń Reynoldsa −ρ<u′u′>, zamykany przez lepkość wirową (modele k-ω, k-ε). Symulacja dużych wirów (LES): rozwiązuje wiry większe niż skala filtra, modeluje tylko naprężenia podsiatkowe (model Smagorinsky'ego). Hybrydowe RANS/LES (DES): tanie blisko ścian, dokładne w obszarach oderwania. Wszystkie modele wprowadzają wolne parametry wymagające kalibracji względem eksperymentów.

Przepływ ściśliwy i fale uderzeniowe

Gdy prędkość przepływu zbliża się do lokalnej prędkości dźwięku, założenie nieściśliwości przestaje obowiązywać. Liczba Macha Ma = U/a (a = √(γRT) dla gazu doskonałego) charakteryzuje reżim: Ma < 0,3 nieściśliwy; 0,3–0,8 poddźwiękowy ściśliwy; 0,8–1,2 transoniczny (z lokalnymi kieszeniami naddźwiękowymi); 1,2–5 naddźwiękowy; >5 hipersoniczny.

Przy Ma > 1 tworzy się fala uderzeniowa normalna: przepływ gwałtownie zwalnia, gęstość i ciśnienie skacze nieciągle, a entropia rośnie (fale ukośne w 2D/3D). Warunki skoku Rankine'a-Hugoniota wiążą właściwości przed i za falą:

Kluczowe zastosowanie: samolot naddźwiękowy musi przyspieszyć przez transoniczną barierę dźwięku (współczynnik oporu osiąga szczyt przy Ma ≈ 1 z powodu oporu falowego) i operować naddźwiękowo tam, gdzie opór falowy jest opanowany przez skrzydła skośne i regułę powierzchni. Osobliwość Prandtla-Glauerta (nieskończona siła nośna przy Ma = 1 w teorii linearyzowanej) jest artefaktem modelu — rzeczywiste przepływy są ciągłe, ale nieliniowe.

Analiza wymiarowa i podobieństwo

Twierdzenie Buckinghama Π stwierdza, że fizycznie sensowne równanie wiążące n wielkości wymiarowych z k niezależnymi wymiarami można zapisać za pomocą n−k bezwymiarowych grup (parametrów Π). Dynamika płynów wykorzystuje to nieustannie: liczby Re, Ma, Fr (Froude'a), We (Webera), Eu (Eulera) i Strouhala. Dwa przepływy są dynamicznie podobne, jeśli wszystkie istotne grupy bezwymiarowe się zgadzają — to podstawa testów w tunelu aerodynamicznym (modele w skali przy dopasowanym Re) i testów w basenie holowniczym statków (dopasowane Fr).

Nadciecze i turbulencja kwantowa

W temperaturach poniżej punktu λ (Tλ ≈ 2,17 K) ciekły ⁴He staje się nadcieczą (kondensatem Bosego-Einsteina). Jego lepkość znika całkowicie i płynie bez oporu. Wirowość jest skwantowana: każda cyrkulacja musi być całkowitą wielokrotnością κ = h/m (kwantu cyrkulacji). Rdzenie wirowe mają grubość ∼1 Å i niosą dokładnie jeden kwant κ. Turbulencja kwantowa to plątanina tych skwantowanych włókien wirowych; w dużych skalach wykazuje widmo −5/3 Kołmogorowa, co sugeruje, że klasyczna kaskada wyłania się z kwantowej dynamiki wirów.

Symulacje na tej platformie

Nasza kategoria dynamiki płynów obejmuje praktyczne narzędzia w całym spektrum liczby Reynoldsa:

• Symulator Bernoulliego — przepływomierz Venturiego, rurka Pitota, ciśnienie profilu aerodynamicznego
• Profil aerodynamiczny i siła nośna — profile NACA, kąt przeciągnięcia, biegun L/D
• Fale w wannie — dyspersja i odbicia fal płytkowodnych
• Przepływ krwi — profil Poiseuille'a, zwężenie, tętniak
• Formowanie chmur — konwekcja wilgotna, niestabilność gradientu temperatury
• Cyklon i Coriolis — równowaga geostroficzna, komórki Hadleya
• Konwekcja Bénarda — komórki Rayleigh-Bénarda, liczba Nusselta
• Dynamika baniek — równanie Rayleigha-Plesseta, kawitacja