PL 
EnglishDomyślny
Українська🌊 Niestabilności płynów
Poznaj niestabilności Rayleigha-Taylora i Kelvina-Helmholtza powstające na styku płynów o różnej gęstości.
O tej symulacji
Główny panel renderuje granicę między dwiema warstwami płynu: lżejszą, szybciej płynącą warstwą na górze oraz cięższą, wolniejszą warstwą na dole, rozdzielonymi cienkim pasem przejściowym. Pole barwnika obliczane na GPU jest adwekowane przez analityczny profil prędkości niestabilności Kelvina-Helmholtza — gładki profil ścinania w kształcie tangensa hiperbolicznego plus narastające zaburzenie falowe — dzięki czemu można obserwować, jak granica zwija się w charakterystyczne wiry typu „koci oczy" w miarę rozwoju niestabilności.
Pod animacją symulator śledzi klasyczną relację dyspersji KH γ(k): destabilizujący człon pochodzący od ścinania prędkości rywalizuje ze stabilizującymi członami grawitacji (poprzez różnicę gęstości między warstwami) i napięcia powierzchniowego. Dolny wykres przemiata zakres liczb falowych, pokazując, która długość fali rośnie najszybciej przy aktualnych ustawieniach — ten szczyt to wartość „Most unstable k" (najbardziej niestabilne k) podana w panelu statystyk.
Ścinanie prędkości ΔU, stosunek gęstości ρ₂/ρ₁, liczbę falową k i napięcie powierzchniowe σ można zmieniać w czasie rzeczywistym; przycisk Pauza zatrzymuje animację do analizy, a Reset odtwarza od nowa granicę międzyfazową i zeruje zegar. Gdy człony stabilizujące przeważają nad ścinaniem, tempo wzrostu wynosi zero, a granica pozostaje płaska zamiast się zwijać.
Najczęściej zadawane pytania
Co decyduje o tym, czy granica między płynami jest stabilna, czy niestabilna?
Stabilność granicy między płynami zależy od rywalizacji między siłami destabilizującymi (grawitacją działającą na różnice gęstości lub ścinaniem wywołującym różnice ciśnienia) a siłami stabilizującymi (napięciem powierzchniowym i lepkością). W przypadku niestabilności Rayleigha-Taylora każda inwersja gęstości (cięższy płyn nad lżejszym) jest bezwarunkowo niestabilna przy braku napięcia powierzchniowego. Napięcie powierzchniowe stabilizuje krótkie długości fali (wysokie liczby falowe), więc rosną tylko zaburzenia dłuższe niż krytyczna długość fali (długość kapilarna).
Jak szybko rosną niestabilności Rayleigha-Taylora?
Tempo wzrostu niestabilności Rayleigha-Taylora wynosi σ = √(Agk), gdzie A = (ρ₂−ρ₁)/(ρ₂+ρ₁) to liczba Atwooda (kontrast gęstości), g to przyspieszenie grawitacyjne, a k = 2π/λ to liczba falowa zaburzenia. Większy kontrast gęstości i krótsze długości fali rosną szybciej w fazie liniowej (wczesnej). W pełni nieliniowej fazie charakterystyczna prędkość wznoszących się i opadających kolumn płynu skaluje się jak √(AgL), gdzie L to wysokość warstwy mieszania.
Gdzie w przyrodzie obserwujemy niestabilności Kelvina-Helmholtza?
Niestabilności Kelvina-Helmholtza pojawiają się wszędzie tam, gdzie stykają się płyny poruszające się z różnymi prędkościami. Charakterystyczne fale-chmury (chmury warkoczowe) powstają tam, gdzie szybki górny przepływ powietrza napotyka wolniejsze, wilgotniejsze powietrze poniżej. Pasma chmur Jowisza pokazują wiry KH wzdłuż granic prądów strumieniowych. W oceanie niestabilności KH na granicach termokliny mieszają ciepłą wodę powierzchniową z zimną wodą głębinową, napędzając wypływ składników odżywczych ku górze. Na magnetopauzie Ziemi ścinanie wiatru słonecznego wywołuje fale KH, które transportują plazmę słoneczną do magnetosfery.
Jakiego wzoru na tempo wzrostu używa ten symulator niestabilności Kelvina-Helmholtza?
Panel sterowania pokazuje klasyczną relację dyspersji Kelvina-Helmholtza dla dwóch warstw: γ = k·ΔU/2 − √(k²g·Δρ/ρ̄ − k²ΔU²/4·ρ₁ρ₂/ρ̄²), gdzie ΔU to ścinanie prędkości między warstwami, Δρ = ρ₂−ρ₁ i ρ̄ = ρ₁+ρ₂ to odpowiednio różnica i suma gęstości, g to przyspieszenie grawitacyjne, a k to liczba falowa zaburzenia. Pod pierwiastkiem destabilizujący człon ścinania (∝ ΔU²) rywalizuje ze stabilizującymi członami grawitacji i napięcia powierzchniowego; gdy człony stabilizujące przeważają, wyrażenie pod pierwiastkiem nie jest dodatnie i granica jest stabilna (tempo wzrostu podawane jako zero). Gdy wygrywa ścinanie, mod rośnie wykładniczo, a symulator animuje odpowiadające mu zwijanie się wirów.
Co kontrolują suwaki ścinania prędkości, stosunku gęstości, liczby falowej i napięcia powierzchniowego?
Ścinanie prędkości ΔU ustala, o ile szybciej porusza się górna (lżejsza) warstwa względem dolnej (cięższej) — większe ścinanie wtłacza więcej energii kinetycznej w granicę i napędza szybsze zwijanie się wirów. Stosunek gęstości ρ₂/ρ₁ ustala, o ile cięższa jest dolna warstwa od górnej; większe stosunki zwiększają stabilizujący człon grawitacyjno-wyporowy, czyniąc granicę bardziej odporną na ścinanie. Liczba falowa k wybiera, którą długość fali zaburzenia (2π/k) śledzi widok główny i odczyt tempa wzrostu; wykres dyspersji poniżej pokazuje tempo wzrostu w całym zakresie liczb falowych, dzięki czemu widać, która rośnie najszybciej. Napięcie powierzchniowe σ dodaje dodatkowy człon stabilizujący, który tłumi zwłaszcza krótkie długości fali (wysokie k), podobnie jak rzeczywiste napięcie powierzchniowe przeciwdziała ostrej krzywiźnie granicy.
Dlaczego wskaźnik „Stabilne?" czasem pokazuje „Tak", a zwijanie się wirów przestaje rosnąć?
Symulacja zgłasza stan „stabilny", gdy destabilizujący człon ścinania nie jest wystarczająco duży, by przezwyciężyć łączny stabilizujący wpływ grawitacji (wyporu wynikającego z różnicy gęstości) i napięcia powierzchniowego dla wybranej liczby falowej. Dzieje się tak, gdy zmniejszysz ścinanie prędkości, zwiększysz stosunek gęstości, zwiększysz napięcie powierzchniowe lub przesuniesz suwak liczby falowej poza pasmo najszybszego wzrostu widoczne jako szczyt krzywej dyspersji. W takim reżimie tempo wzrostu γ jest podawane jako zero, a granica pozostaje w zasadzie płaska, zamiast zwijać się w wiry.
Jaką rolę odgrywa napięcie powierzchniowe w stabilizacji warstwy ścinania?
Napięcie powierzchniowe przeciwdziała każdemu wzrostowi krzywizny granicy, więc szczególnie hamuje zaburzenia o krótkiej długości fali (wysokiej liczbie falowej), które przy danej amplitudzie mają najostrzejszą krzywiznę. We wzorze na tempo wzrostu wchodzi ono jako dodatkowy człon stabilizujący rosnący z k³, więc zwiększanie suwaka napięcia powierzchniowego tłumi głównie drobne, gęsto rozmieszczone zmarszczki, pozostawiając względnie nienaruszonymi mody o dłuższej fali. Odzwierciedla to zachowanie rzeczywistych granic między płynami, gdzie napięcie powierzchniowe wyznacza krótkofalowe odcięcie, poniżej którego niestabilności napędzane ścinaniem lub wyporem nie mogą rosnąć, niezależnie od siły destabilizującej.