Konwekcja Bénarda: Porządek z Rozgrzanego Warstwy Płynnej

Podgrzej cienką warstwę płynu z dołu

Weź cienką warstwę płynu, podgrzej ją równomiernie od dołu i schłodź od góry, a dla małych różnic temperatur nie obserwuje się nic widocznego – ciepło po prostu przewodzi w górę przez nieruchomy płyn. Doprowadź różnicę temperatur ponad krytyczny próg, a płyn spontanicznie uorganizuje się w regularny wzór obrotowych ciał powietrznych lub heksagonalne komórki: gorący, mniej gęsty płyn unosi się w niektórych obszarach, a zimny, bardziej gęsty płyn opada w innych, tworząc zaskakująco uporządkowany wzorzec, który wyłonił się z całkowicie jednolitego stanu początkowego. Jest to konwekcja Bénarda Rayleigha, jeden z klasycznych przykładów spontanicznej formacji wzorców w fizyce.

Liczba Rayleigh’a: Jedna skala, która decyduje o wszystkim

Czy zachodzi konwekcja w ogóle jest regulowana przez jedną bezwymiarową wielkość, liczbę Rayleigh’a, która porównuje siły wypornościowe prowadzące do przelewania z dwoma efektami, które je opóźniają — lepkość, która tłumi ruch płynu, i dyfuzja termiczna, która usuwa różnice temperatur przed mogą one napędzać przepływ:

Ra = ( g * alpha * DeltaT * L^3 ) / ( nu * kappa ) g = przyspieszenie ziemskie alpha = współczynnik rozszerzalności cieplnej DeltaT = różnica temperatur w kierunku pionowym L = głębokość warstwy nu = lepkość kinematyczna kappa = dyfuzja termiczna Ra czyste przewodzenie, brak ruchu Ra > Ra_c -> konwekcja w postaci rulonów pojawia się Krytyczna wartość Ra_c ≈ 1708 dla warstwy z dwoma nieruchomymi granicami została pierwotnie wyprowadzona analitycznie przez fizyka Lorda Rayleigha w 1916 roku, bazując na obserwacjach eksperymentalnych Henri Bénarda, który w 1900 roku podgrzał cienką warstwę wosku spermaceti i sfotografował powstałe sześciokątne komórki. Poniżej Ra_c przewodzenie zawsze wygrywa — każdy początkujący ruch konwekcyjny jest tłumiony szybciej niż siły wypornościowe mogą go utrzymać. Przekrocz próg i płyn nagle reorganizuje się w rulony, ponieważ ten wzór przepływu konwekcyjnego znacznie bardziej wydajnie przenosi ciepło do góry niż przewodzenie samo w sobie.

Ra = ( g * alpha * DeltaT * L^3 ) / ( nu * kappa )
  g      = gravitational acceleration
  alpha  = thermal expansion coefficient
  DeltaT = temperature difference top-to-bottom
  L      = layer depth
  nu     = kinematic viscosity
  kappa  = thermal diffusivity

Ra < Ra_c (~1708 for rigid boundaries) -> pure conduction, no motion
Ra > Ra_c                              -> convection rolls appear

Dlaczego się formują wałki i heksagonalne komórki

Tuż nad Ra_c, preferowanym wzorem w idealnie uwarstwionym płynie z ustalonymi temperaturami na brzegach jest zestaw równoległych wałków konwekcyjnych: płyn unosi się po jednej stronie wałka, przemieszcza się po wierzchu, opada po drugiej stronie i wraca wzdłuż spodu, a sąsiednie wałki przeciwstawnie obracają się, tak że wspólny brzeg ma płyn poruszający się w tym samym kierunku. Jeśli natomiast górną powierzchnię pozostawia się wolną (otwartą na powietrze) a napięcie powierzchniowe jest istotne, ta sama niestabilność częściej organizuje się w komórki heksagonalne – oryginalne komórki Bénarda napędzane napięciem powierzchniowym (technicznie efekt Marangoni) różniły się od czysto napędzanego konwekcją ciągu wypornościowego, a ta rozróżnienie zostało sprecyzowane dopiero dziesiątki lat po jego oryginalnych eksperymentach.

Przybliżenie Boussinesqa i równania Navier-Stokesa

Symulowanie tego właściwie oznacza rozwiązywanie niezmienniczych równań Navier-Stokesa połączonych z równaniem transportu ciepła, ale przy użyciu jednej sztuczki upraszczającej: przybliżenie Boussinesqa traktuje gęstość płynu jako stałą wszędzie, z wyjątkiem członu siły wypornościowej, gdzie niewielkie zmiany gęstości spowodowane rozszerzalnością cieplną są tym, co faktycznie napędza przepływ. Pozwala to utrzymać prosty warunek zachowania masy w równaniu ruchu (przepływ niezmienniczy, zerowa dywergencja) przy jednoczesnym umożliwieniu pracy siłom wypornościowym naprowadzonym przez temperaturę.

Boussinesq momentum equation (2D, simplified):
du/dt = -grad(p)/rho0 + nu * laplacian(u) + alpha * g * (T - T0) * y_hat
dT/dt = kappa * laplacian(T) - u . grad(T)
div(u) = 0                                    // incompressibility retained exactly

Beyond onset: chaos at high Rayleigh number

Push Ra far beyond Ra_c and the orderly rolls themselves become unstable — first developing wavy, oscillating distortions, then time-periodic and eventually chaotic motion, a route to turbulence that Edward Lorenz’s famous three-variable truncation of these very equations (the Lorenz system) was originally derived to study. At very high Rayleigh number the flow becomes fully turbulent convection, the regime relevant to the Sun’s convective zone, Earth’s mantle, and the atmosphere’s convective boundary layer — the same instability mechanism operating across vastly different scales and fluids.

Co pokazuje symulacja

Ta symulacja rozwiązuje 2D model Boussinesqa Navier-Stokes płynu podgrzanego, pozwalając na podniesienie efektywnego liczby Rayleigh'a ponad próg krytyczny i obserwację nagłego przejścia od przewodzenia do uporządkowanej konwekcji w postaci rulonów, a następnie obserwację destabilizacji tych rulonów w bardziej złożone, zmienne w czasie przepływy wraz ze wzrostem dodatkowej różnicy temperatur.

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Bénard Convection i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.