PL 
EnglishDomyślny
Українська🌊 Konwekcja Bénarda
Nieściśliwe równania Naviera-Stokesa z wyporem w przybliżeniu Boussinesqa. Obserwuj powstawanie komórek konwekcyjnych, gdy ogrzany płyn unosi się, a schłodzony opada w cienkiej poziomej warstwie.
Co pokazuje
Konwekcja Rayleigha-Bénarda zachodzi, gdy płyn jest ogrzewany od dołu. Poniżej krytycznej liczby Rayleigha Rac ≈ 1708 ciepło przepływa wyłącznie przez przewodzenie. Powyżej Rac niestabilność wyporu napędza tworzenie się przeciwbieżnie wirujących wałków konwekcyjnych (komórek Bénarda). Ta symulacja rozwiązuje dwuwymiarowe nieściśliwe równania Naviera-Stokesa w przybliżeniu Boussinesqa, wykorzystując sformułowanie funkcji prądu-wirowości (ψ-ω).
Jak korzystać
- Zwiększ Buoyancy, aby wywołać silniejszą konwekcję — powstanie więcej i szybszych wałków
- Zwiększ Viscosity, aby wytłumić przepływ — komórki staną się ospałe i mogą zniknąć
- Przełącz strzałki prędkości, by zobaczyć kierunek cyrkulacji
- Kliknij Perturb, aby dodać losowy szum temperatury i wywołać tworzenie nowych komórek
Czy wiesz, że...
Henri Bénard po raz pierwszy zaobserwował te komórki konwekcyjne w 1900 roku, używając wosku spermacetowego. Lord Rayleigh wyprowadził kryterium stabilności w 1916 roku. Dziś konwekcja Rayleigha-Bénarda leży u podstaw granulacji słonecznej, oceanicznej cyrkulacji termohalinowej i konwekcji płaszcza Ziemi — wszystkie napędzane tą samą niestabilnością.
O konwekcji Rayleigha-Bénarda
Ta symulacja pokazuje warstwę płynu ogrzewaną od dołu i chłodzoną od góry — klasyczny układ Rayleigha-Bénarda. Rozwiązuje dwuwymiarowe nieściśliwe równania Naviera-Stokesa w przybliżeniu Boussinesqa, wykorzystując sformułowanie funkcji prądu-wirowości (ψ-ω) na siatce 120×72. Temperatura i wirowość są adwekowane metodą semi-Lagrange'a, dyfundowane jawnie, a równanie Poissona ∇²ψ = -ω jest relaksowane metodą kolejnych nadrelaksacji. Wypór wchodzi jako człon źródłowy proporcjonalny do poziomego gradientu temperatury.
Suwak Buoyancy skaluje siłę napędową (efektywnie liczbę Rayleigha), Viscosity ν tłumi pęd, a Diffusivity κ wygładza temperaturę. Powyżej krytycznej liczby Rayleigha Ra ≈ 1708 nieruchoma, przewodząca warstwa staje się niestabilna i spontanicznie powstają przeciwbieżnie wirujące wałki konwekcyjne. Ta sama niestabilność rządzi granulacją słoneczną, cyrkulacją atmosferyczną i oceaniczną oraz powolną konwekcją płaszcza Ziemi, co czyni ją fundamentem dynamiki płynów.
Najczęściej zadawane pytania
Czym jest konwekcja Rayleigha-Bénarda?
To ruch napędzany wyporem, który powstaje, gdy pozioma warstwa płynu jest ogrzewana od dołu i chłodzona od góry. Ciepły, mniej gęsty płyn w pobliżu gorącej dolnej płyty unosi się, podczas gdy chłodny płyn opada, organizując się w regularne wałki konwekcyjne lub heksagonalne komórki Bénarda, gdy ogrzewanie jest wystarczająco silne.
Co oznacza liczba Rayleigha?
Liczba Rayleigha Ra = gβΔTH³/(νκ) porównuje siłę napędową wyporu z tłumiącymi efektami lepkości i dyfuzji cieplnej. Gdy Ra przekracza wartość krytyczną wynoszącą około 1708 dla sztywnych płyt, samo przewodzenie nie jest już w stanie odprowadzić ciepła i rozpoczyna się konwekcja.
Co kontrolują trzy suwaki?
Buoyancy skaluje siłę napędową termiczną, pełniąc rolę liczby Rayleigha; jej zwiększenie daje bardziej intensywne wałki. Viscosity tłumi pole prędkości, spowalniając lub tłumiąc komórki, natomiast Diffusivity kontroluje, jak szybko wyrównują się różnice temperatur.
Jakiej metody numerycznej używa symulacja?
Wykorzystuje formułę funkcji prądu i wirowości dla dwuwymiarowych równań Naviera-Stokesa z przybliżeniem Boussinesqa. Temperatura i wirowość są adwekowane schematem semi-Lagrange'a, dyfundowane jawnie, a funkcja prądu jest odtwarzana w każdym kroku przez rozwiązanie równania Poissona metodą kolejnych nadrelaksacji.
Dlaczego wałki pojawiają się dopiero powyżej pewnego progu?
Poniżej krytycznej liczby Rayleigha lepkość i dyfuzja cieplna wygaszają każde małe zaburzenie, zanim zdąży ono urosnąć, więc ciepło przemieszcza się wyłącznie przez przewodzenie. Powyżej progu siła wyporu wygrywa, drobne zaburzenia się wzmacniają, a warstwa reorganizuje się w stały ruch konwekcyjny. To podręcznikowy przykład niestabilności hydrodynamicznej i bifurkacji widłowej.
Czym jest formuła funkcji prądu i wirowości?
Zamiast śledzić bezpośrednio ciśnienie i dwie składowe prędkości, przepływ opisywany jest przez wirowość ω (lokalny obrót) i funkcję prądu ψ, której gradienty dają prędkość. To automatycznie wymusza nieściśliwość i usuwa ciśnienie z równań, dlatego jest to wygodne podejście dla takich dwuwymiarowych symulacji jak ta.
Co pokazują strzałki prędkości i linie prądu?
Strzałki prędkości oznaczają lokalny kierunek i względną prędkość płynu, ujawniając cyrkulację wewnątrz każdego wałka. Opcja linii prądu rysuje kontury funkcji prądu ψ; płyn przepływa wzdłuż tych linii, więc zamknięte pętle wyznaczają komórki konwekcyjne bez ich przecinania.
Czy ta symulacja jest fizycznie dokładna?
Odwzorowuje poprawną jakościową fizykę: przejście od przewodzenia do konwekcji, przeciwbieżnie wirujące wałki oraz zależność od wyporu, lepkości i dyfuzyjności. Jednak siatka jest zgrubna, wartości suwaków są skalowane, a nie wyrażone w prawdziwych jednostkach SI, a symulacja jest dwuwymiarowa, więc to narzędzie dydaktyczne, a nie solver klasy badawczej.
Co robi przycisk Perturb?
Perturb wprowadza małe losowe fluktuacje temperatury do wnętrza warstwy. W pobliżu lub powyżej progu niestabilności te zarodki rosną i mogą wywołać nowe komórki konwekcyjne lub przeorganizować istniejący wzór, pozwalając obserwować, jak układ wybiera nową strukturę przepływu z szumu.
Gdzie w przyrodzie występuje konwekcja Rayleigha-Bénarda?
Ta sama niestabilność napędza granulację na powierzchni Słońca, pasma chmur i prądy wznoszące w burzach w atmosferze, cyrkulację termohalinową w oceanach oraz powolną konwekcję płaszcza Ziemi, która przesuwa płyty tektoniczne. Henri Bénard po raz pierwszy zaobserwował te komórki w 1900 roku, a Lord Rayleigh wyjaśnił kryterium w 1916 roku.
Jakich warunków brzegowych używa model?
Dolna płyta jest utrzymywana w cieple przy T = 1, a górna w chłodzie przy T = 0, co ustala spadek temperatury w warstwie. Płyty są sztywnymi ścianami bez poślizgu, narzucanymi przez standardową relację wirowości na ścianie, natomiast kierunek poziomy jest traktowany jako okresowy, dzięki czemu wałki mogą zawijać się wokół obszaru.