Strona główna Dynamika Płynów Wirowość płynu i niestabilność Kelvina–Helmholtza

🌊 Wirowość płynu i niestabilność Kelvina–Helmholtza

Dynamika dwuwymiarowej warstwy wirowej i niestabilność Kelvina–Helmholtza. N wirów punktowych oddziałuje wg prawa Biota–Savarta, zwijając powierzchnię ścinania w charakterystyczne struktury spiralne.

Dynamika Płynów3DŚredni60 FPS
fluid-vorticity ↗ Otwórz osobno
Interfejs samej symulacji jest w języku angielskim.

O niestabilności Kelvina-Helmholtza

Ta symulacja modeluje dwuwymiarową nielepką warstwę wirową na granicy między dwoma strumieniami cieczy poruszającymi się z różnymi prędkościami. Ciągła warstwa jest dyskretyzowana na N punktowych wirów, z których każdy niesie cyrkulację Γ i porusza się pod wpływem pola prędkości indukowanego przez wszystkie pozostałe poprzez regularyzowane prawo Biota-Savarta: vx = −(Γ/2π)·dy/r², vy = +(Γ/2π)·dx/r², przy czym skończony rdzeń ε zapobiega osobliwościom. Okresowość w kierunku x jest wymuszana za pomocą kolumn obrazów.

Suwaki ustawiają liczbę wirów N (8–64), cyrkulację Γ na wir, promień rdzenia ε, jawny krok czasowy Eulera dt oraz amplitudę początkowego zaburzenia sinusoidalnego. Większe zaburzenia zapoczątkowują szybsze zwijanie. Przyciski Pojedyncza, Podwójna i Losowa warstwa zmieniają konfigurację początkową. Ta niestabilność ścinająca napędza turbulencję w chmurach kłębiastych, falach oceanicznych, atmosferze Jowisza i wietrze słonecznym, co czyni ją fundamentalną dla mechaniki płynów.

Najczęściej zadawane pytania

Czym jest niestabilność Kelvina-Helmholtza?

To niestabilność powstająca, gdy dwie warstwy cieczy przesuwają się względem siebie z różnymi prędkościami. Każde małe zaburzenie na wspólnej granicy rośnie wykładniczo, zwijając granicę w rząd spiralnych wirów. Jest to jedna z najczęstszych dróg do turbulencji w przyrodzie.

Jak ta symulacja oblicza ruch?

Warstwa wirowa jest reprezentowana przez N dyskretnych wirów punktowych. W każdym kroku każdy wir jest przemieszczany przez łączną prędkość indukowaną przez wszystkie pozostałe poprzez nielepkie prawo Biota-Savarta, całkowane w przód prostym jawnym schematem Eulera z krokiem czasowym dt.

Dlaczego istnieje promień rdzenia ε?

Czysty wir punktowy indukuje prędkość, która rozbiega się jak 1/r, gdy dwa wiry się zbliżają, powodując niestabilność numeryczną. Dodanie ε² do mianownika regularyzuje jądro, nadając każdemu wirowi gładki, skończony rdzeń. Utrzymuje to stabilność symulacji, ledwie wpływając na zwijanie w dużej skali.

Co robią suwaki Liczba wirów N i Cyrkulacja Γ?

N (8–64) ustala, jak dokładnie ciągła warstwa wirowa jest próbkowana przez dyskretne wiry; więcej wirów daje ostrzejsze spiralne rdzenie. Γ to cyrkulacja niesiona przez każdy wir, skalująca siłę indukowanego pola prędkości, a tym samym szybkość rozwoju niestabilności.

Jaka jest szybkość wzrostu niestabilności?

Dla warstwy wirowej zaburzenie o liczbie falowej k rośnie wykładniczo z szybkością σ = |Γ|k/4π w granicy nielepkiej. Krótsze długości fal rosną najszybciej, dlatego zdyskretyzowana warstwa szybko tworzy małoskalowe zwoje, które później łączą się w większe struktury poprzez parowanie wirów.

Co pokazuje tryb Podwójna warstwa?

Podwójna warstwa dodaje drugi, przeciwnie wirujący rząd wirów około 0,08 jednostki nad pierwszym. Przeciwstawne cyrkulacje tworzą układ przypominający ścieżkę wirową Kármána, powstającą za opływanym ciałem, ilustrując interakcję sparowanych warstw ścinających.

Czy symulacja jest fizycznie dokładna?

Jest to dokładny dwuwymiarowy nielepki model wirów punktowych, więc wiernie oddaje wczesne zwijanie i parowanie przewidywane przez klasyczną teorię. Pomija jednak lepkość, efekty trójwymiarowe i ściśliwość, więc bardzo późne chaotyczne etapy są jakościowe, a nie ilościowo precyzyjne.

Dlaczego przepływ staje się chaotyczny z czasem?

Układy wirów punktowych z więcej niż trzema wirami są na ogół niecałkowalne i wykazują wrażliwą zależność od warunków początkowych. Gdy wiry zwijają się i parują, drobne różnice narastają, tworząc chaotyczny, zdominowany przez mieszanie stan charakterystyczny dla rozwiniętej turbulencji ścinającej.

Dlaczego wiry są kolorowane inaczej?

Niebieskie wiry niosą dodatnią cyrkulację (Γ > 0), a czerwone — ujemną (Γ < 0), wskazując przeciwne kierunki rotacji. W trybie Pojedyncza warstwa wszystkie wiry mają jeden znak, podczas gdy tryby Podwójna i Losowa mieszają oba znaki, co widać na pozostawianych przez nie śladach.

Gdzie w rzeczywistym świecie występuje niestabilność Kelvina-Helmholtza?

Powoduje ona przypominające łamiące się fale chmury Kelvina-Helmholtza na niebie, fale wewnętrzne w oceanie, warstwy mieszania w gwiazdach, niestabilności plazmowe przy magnetopauzie Ziemi oraz pasmowe przepływy ścinające na Jowiszu. Ta sama matematyka leży u podstaw odrywania wirów za przeszkodami i mieszania w strugach.

Podobne symulacje