Wave 53: Przejście laminarno-turbulentne, geometria islamska i warstwa przyścienna Blasiusa

Wave 53 w skrócie

〜 Przejście od laminarnego do turbulentnego

W 1883 roku Osborne Reynolds wstrzyknął cienką nitkę barwnika do rury i zmieniał prędkość przepływu. Poniżej prędkości krytycznej barwnik tworzył idealnie prostą nić — przepływ laminarny. Powyżej pewnego progu nić nagle rozpadała się i mieszała w całym przekroju rury — przepływ turbulentny. Bezwymiarowym parametrem rządzącym tym przejściem jest liczba Reynoldsa:

Laminarny profil prędkości

Dla stałego, w pełni rozwiniętego przepływu laminarnego w rurze okrągłej (przepływ Hagena-Poiseuille'a) profil prędkości jest dokładnie paraboliczny:

u(r) = u_max (1 − r²/R²)

gdzie r to odległość radialna od osi rury, a R to promień rury. Maksymalna prędkość na osi jest dwukrotnością średniej przekrojowej: u_max = 2 Ū. Wszystkie elementy płynu poruszają się po prostych, równoległych torach — nie zachodzi żadne mieszanie.

Turbulentny profil prędkości

W w pełni turbulentnym przepływie w rurze (Re ≫ 10 000) uśredniony w czasie profil prędkości jest znacznie płaszy. Empiryczne prawo potęgi 1/n daje dobre przybliżenie:

u(r)/u_c ≈ (1 − r/R)^1/n

gdzie n ≈ 7 dla Re ~ 10⁵. Na ten uśredniony profil nakładają się turbulentne fluktuacje prędkości u′(t) i v′(t). Ich pierwiastek średniokwadratowy w stosunku do średniej to intensywność turbulencji I = u_rms/Ū, zwykle 5–15% w w pełni rozwiniętej turbulencji.

Symulacja używa znaczników cząsteczkowych: kolorowe strumienie wchodzą z lewej strony. W trybie laminarnym ich tory pozostają proste, a pasma kolorów nigdy się nie mieszają. Gdy Re przekracza próg przejścia, fluktuacje składowej v rozpraszają cząstki w poprzek kanału, odtwarzając mieszanie barwnika obserwowane w oryginalnym eksperymencie Reynoldsa.

✦ Islamskie wzory geometryczne

Islamska sztuka geometryczna — spotykana w meczetach, pałacach i manuskryptach od Andaluzji po Azję Środkową — zbudowana jest z niewielkiego zestawu przenikających się gwiazd wielokątnych, rozet i kafelków łączących. Wzory wykazują idealną symetrię translacyjną i (lokalnie) obrotową, nie pozostawiając luk ani nakładek.

Gwiazdy wielokątne {n/k}

Regularna gwiazda wielokątna {n/k} jest rysowana przez umieszczenie n równo rozmieszczonych punktów na okręgu i połączenie co k-tego punktu. Notację wprowadził Schläfli. Kluczowe przykłady:

• {6/2} — Gwiazda Dawida (heksagram); gcd(6,2) = 2 daje dwa oddzielne trójkąty równoboczne.
• {8/3} — gwiazda 8-ramienna (oktagram); gcd(8,3) = 1 daje pojedynczą zamkniętą gwiazdę 8-punktową — wszechobecną w osmańskich kafelkach i wzorach girih.
• {10/4} — 10-ramienna gwiazda girih; kanoniczny element systemu kafelków girih znaleziony w sanktuarium Darb-i Imam w Isfahanie (1453 r.). Pozostałe kształty luk (kokarda, pięciokąt, wydłużony sześciokąt, romb) wypełniają płaszczyznę bez nakładek.
• {12/5} — 12-ramienna gwiazda osmańska; częsta w tureckich kafelkach İznik i marokańskich mozaikach zellige.

Logika kafelkowania

Generator umieszcza środki gwiazd na siatce kwadratowej (odstęp d) lub heksagonalnej (dx = d, dy = d√3⁄2 z naprzemiennym przesunięciem rzędów). Dla każdego środka rysowane są ścieżki wielokąta {n/k}, opcjonalnie wypełniane. Mniejszy regularny n-kąt w środku (promień ≈ 0,45 R) tworzy wewnętrzny obszar wypełnienia o zmniejszonej krycie, naśladując kontrast między korpusem gwiazdy a otaczającymi kafelkami łączącymi.

Dostępnych jest sześć historycznych palet kolorystycznych: Islamic Gold (bursztyn i ochra), Alhambra (leśna zieleń), Ottoman Blue (kobalt i błękit nieba), Moroccan Terracotta (spalona sjena), Mughal White (biel na czerni) i Neon Night (fiolet na niemal czarnym tle). Gotowy wzór można wyeksportować jako bezstratny PNG.

〰 Warstwa przyścienna Blasiusa

Gdy lepki płyn opływa powierzchnię stałą, warunek braku poślizgu wymusza zerową prędkość płynu przy ścianie. Powstaje cienka warstwa — warstwa przyścienna — w której prędkość rośnie od zera do wartości swobodnego strumienia U∞. Poza tą warstwą przepływ jest zasadniczo nielepki.

Równanie Blasiusa

W latach 1907–1908 Heinrich Blasius (student Prandtla) pokazał, że dla stałego przepływu laminarnego nad półnieskończoną płaską płytą równania Naviera-Stokesa redukują się, przy przybliżeniach warstwy przyściennej, do pojedynczego równania różniczkowego zwyczajnego trzeciego rzędu poprzez zmienną podobieństwa η = y √(U∞ / νx):

f‴ + ½ f f″ = 0   z warunkami   f(0) = f′(0) = 0, f′(∞) = 1

Tutaj f′(η) = u/U∞ to bezwymiarowy profil prędkości. Rozwiązanie jest uniwersalne — wszystkie płaskie laminarne warstwy przyścienne pokrywają się na tej jednej krzywej, gdy wykreślone są względem η.

Metoda strzałów i RK4

Warunek brzegowy f′(∞) = 1 nie jest dany przy η = 0, więc nie można go bezpośrednio całkować. Zamiast tego stosuje się metodę strzałów: zgaduje się wartość początkową α = f″(0) ≈ 0,4696 (dobrze znaną stałą naprężenia ścinającego Blasiusa na ścianie), całkuje w przód metodą RK4 z krokiem Δη = 0,02 i weryfikuje, że f′ → 1 przy η → 12. Powstała tabela rozwiązań jest wykorzystywana w czasie działania metodą interpolacji liniowej do obliczenia u/U∞ w dowolnym punkcie (x, y) na canvasie.

Kluczowe wyniki

• Grubość warstwy przyściennej: δ₉₉(x) ≈ 5,0 √(νx/U∞)
• Lokalne tarcie powierzchniowe: C_f(x) = 0,664 / √Re_x
• Średnie tarcie powierzchniowe na długości płyty L: C̄_f = 1,328 / √Re_L
• Naprężenie ścinające na ścianie: τ_w(x) = 0,332 ρ U∞² / √Re_x

Canvas pokazuje płaską płytę u dołu z obwiednią warstwy przyściennej δ₉₉(x) (pomarańczowe wypełnienie), profilami prędkości u/U∞ w czterech stacjach x (białe krzywe) oraz strzałkami swobodnego strumienia powyżej warstwy. Przełączanie między wodą (ν = 1 × 10⁻⁶ m²/s) a powietrzem (ν = 1,5 × 10⁻⁵ m²/s) natychmiast aktualizuje grubość warstwy przyściennej i wartości tarcia powierzchniowego.

Pytania lub sugestie? Napisz przez stronę kontaktową.