Strona głównaArtykułyWortoniczność cieczy

Płyty wortexowe: Jak N punktowych wortexów symulują warstwę obrotową

Dyskretny zestaw płyty wortexowej z N punktowych wortexów, każdy podporządkowany zasadom Biot-Savarta, powtarza roll-up niepokojącej instabilności Kelvin-Helmholtza bez kiedykolwiek dotykania siatki.

mysimulator teamZaktualizowano — czerwiec 2026≈ 8 min czytania▶ Otwórz symulację

Warstwa ścinająca jako łańcuch wótków punktowych

Zamiast rozwiązywać równania Eulera kontynuum na siatce, ta symulacja przedstawia płaszczyznę wótków między dwiema warstwami ścinających się cieczy jako skończoną liczbę N wótków punktowych, każdy z niezmienną okrążeniem Gamma. Każdy wótek powoduje pole prędkości dla każdego innego wótku poprzez prawo Biot-Savarta, dokładnie tak jak analogia magnetyczna używana dla przewodników przeprowadzających prąd: pole prędkości indukowane maleje jako 1/r i obraca się wokół źródłowego wótku. Integracja prędkości wszystkich wótków w przyszłość czasową stanowi całą symulację — nie ma rozwiązywania ciśnienia, ani siatki, tylko N(N-1) interakcji parowych na krok.

demo na żywo · powiązana symulacja● LIVE

Dlaczego płaszczyzna prostopadła jest niestabilna

Prawie idealnie prosta, równomiernie odstępująca linia wirów stanowi dokładne (ale niewystarczające) równowagę. Nieskończenie niewielka perturbacja pokazuje analiza klassena Kelvin-Helmholtza, że każda długość fali rosnie eksponencjalnie, z dominantnym trywialnym trybem porównywalnym z odstępem między wirami. Fizycznie: mała kreska nawiązuje wiry sąsiednie na jednej stronie bliżej siebie, co zwiększa ich wzajemną prędkość indukowaną, co w konsekwencji zwiększa kreskę — pętlę zwrotną bez siły powrotnej w warstwie nieprzewodzącej i nienapędzonej.

for each vortex i:
  v_i = sum over j != i of  Gamma_j / (2*pi*r_ij) * tangent_ij
  x_i += v_i * dt        // simple Euler, or RK4 for accuracy

Rozkład i singularność Moore'a

Podczas destabilizacji płyty, model punktowych wirów powtarza klasyczny skręcony rozkład w formie pęcherzy dymu, chmur i brzegów strumienia — płyta krzywi się w szereg dyskretnych jąder wirów połączonych cienkimi wirami. W granicy kontynuum ten rozkład jest sławnym dla rozwijania singularności krzywizny o skończonym czasie (singularność Moore'a), po której równania gładkiej płyty wirowej technicznie przestają być prawdziwe; dyskretny model N-punktowy uniknie tego zbudowany, ponieważ zawsze można zaavansować skończoną liczbę punktów, ale to zapłacenie za stabilność polega na powstaniu dymiącego szumu małoskalowego, który w końcu wymaga regularizacji lub okresowej łączenia wirów, aby pozostawać fizyczny.

Gdzie to występuje

Płyty wortexowe Kelvin-Helmholtza powstają tam, gdzie dwa płyny lub warstwy płynu przemieszczają się jedno względem drugiego z nawet małą różnicą prędkości: ostrze chmur cumulus, granica między wiatrem a wolniejszym powietrzeniem nad wodą (generując falę powierzchniową), brzegi boczne Jupitera na tle jego warstw chmur, magnetopauzę, gdzie wiatr słoneczny spotyka się z magnetycznym polem Ziemi, oraz warstwy cięciu w silnikach doświtowych i ductach mieszających. Instabylność jest również głównym mechanizmem mieszania na powierzchniach o różnej gęstości po osłabieniu stratyfikacji (warunek numeru Richarda).

Często zadawane pytania

Jak to różni się od siatkowej symulacji płynnej tej samej instabilności?

Metoda siatki dyskretyzuje przestrzeń i rozwiązuje równania Eulera lub Naviera-Stokesa w każdym komórce. Metoda punktowych wirów dyskretyzuje sam front wiru w N cząstek i oblicza prędkości analitycznie za pomocą zasady Biot-Savarta — bez siatki, bez rozwiązywania ciśnienia, dokładne zachowanie niezmienniczości, ale działa tylko dla obszarów, gdzie wirowość jest skoncentrowana na frontu.

Dlaczego w koniec wiry wyglądają szumowo zamiast czystej spirali?

Nieskończenie ciągły front wirów rozwija singularność krzywizny w skończonym czasie (singularność Moore); skończona liczba punktowych wirów przybliża to, co oznacza, że małe skalne położenia stać się chaotycznie czułe do sąsiednich wirów, co powoduje drgania longitudynowe ramion spirali, chyba że symulacja regularnie łączy lub łagodzi bliskie wirovki.

Co to jest cyrkulacja Gamma i dlaczego jest ona stała dla każdego wiru?

Cyrkulacja to liniowy całkowy prędkości wokół jądra wiru, a zgodnie z twierdzeniem Kelvina o zachowaniu cyrkulacji jest ono utrzymywane dla każdego elementu cieczy w niewiskoslowym, barotropicznym przepływie. Przypisując każdemu punktowemu wirowi stałą Gamma, dokładnie zasady zachowania cyrkulacji, co jest powodem, dla którego metoda powtarza poprawne wieloskalowe dynamy, mimo że używa tak mało stopni swobody.

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Fluid Vorticity & Kelvin–Helmholtz Instability i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.

▶ Otwórz symulację Fluid Vorticity & Kelvin–Helmholtz Instability

Co znalazłeś?

Dodaj kroki odtworzenia (opcjonalnie)