🦋 Lorenz Attractor — 3D Chaos

Symulacja wykorzystuje klasyczną czterorzędową metodę Rungego–Kutty (RK4) ze stałym krokiem czasowym dt = 0,005 i 6 podkrokami na renderowaną klatkę (efektywny postęp 0,03 na klatkę przy 60 fps). RK4 zapewnia dokładność czwartego rzędu — lokalny błąd obcięcia O(dt⁵) — co jest istotne dla sztywnego układu chaotycznego, jakim jest atraktor Lorenza, gdzie prostsza metoda Eulera szybko gromadzi błędy w przestrzeni fazowej. Trajektoria dzięki temu ściśle śledzi prawdziwy matematyczny atraktor przez tysiące kroków, zanim chaotyczna rozbieżność od jakiejkolwiek „prawdziwej” orbity staje się nieunikniona.

Dziwny atraktor to atraktor w przestrzeni fazowej o strukturze fraktalnej — zajmuje zerową objętość, a mimo to trajektorie pozostają w nim uwięzione na zawsze. Atraktor Lorenza ma wymiar Hausdorffa (fraktalny) wynoszący w przybliżeniu 2,06, co oznacza, że jest nieco bardziej złożony niż powierzchnia, ale znacznie mniej wypełnia przestrzeń niż bryła 3D. Sąsiednie trajektorie na atraktorze rozbiegają się ze średnim tempem mierzonym maksymalnym wykładnikiem Lapunowa λ₁ ≈ 0,906 s⁻¹, potwierdzającym chaos deterministyczny: przewidywalność zmniejsza się o połowę mniej więcej co 0,76 jednostki czasu (≈ln 2 / λ₁).

Parametr ρ kontroluje siłę napędową (analogiczną do liczby Rayleigha w konwekcji). Dla ρ < 1 wszystkie trajektorie zbiegają do początku układu (brak konwekcji). Dla 1 < ρ < 24,74 układ ma dwa stabilne punkty stałe (ustalone wałki konwekcyjne). Przy ρ ≈ 24,74 zachodzi podkrytyczna bifurkacja Hopfa i pojawia się chaos. Przy ρ = 28 istnieje słynny dziwny atraktor. W pobliżu ρ = 99,96 ponownie pojawia się orbita okresowa — czasem nazywana oknem okresowym wewnątrz reżimu chaotycznego. Menu Presetów pozwala przeskakiwać między tymi jakościowo różnymi zachowaniami bez ręcznego dostrajania.

σ to liczba Prandtla, stosunek dyfuzyjności pędu (lepkości kinematycznej) do dyfuzyjności cieplnej. Dla powietrza σ ≈ 0,72; Lorenz użył σ = 10, co odpowiada wodzie i daje solidny chaos. Zmniejszenie σ poniżej 1 ma tendencję do tłumienia chaosu. β to czynnik geometryczny związany ze stosunkiem wymiarów komórki konwekcyjnej: dla klasycznego przybliżenia okrągłego wałka β = 8/3 ≈ 2,667. Zmiana β przesuwa położenie skrzydeł punktów stałych i nieznacznie zmienia wymiar fraktalny atraktora, ale rzadko eliminuje chaos przy standardowym ρ=28.

W każdym kroku całkowania chwilowa prędkość jest obliczana jako wielkość wektora prędkości w przestrzeni fazowej: |v| = √(ẋ² + ẏ² + ż²). Jest to mapowane liniowo z zakresu 0–30 na wartość barwy 240°–0° (niebieski dla wolnego, czerwony dla szybkiego). Ślad zazwyczaj pokazuje więc niebieski (wolno) blisko wewnętrznych spirali każdego płata — gdzie trajektoria owija się ciasno — i czerwony (szybko) blisko obszaru siodłowego przejścia między dwoma skrzydłami, czyniąc strukturę dynamiczną natychmiast widoczną.

Bufor kołowy (pierścieniowy) o długości równej wartości suwaka śladu przechowuje tylko N najnowszych punktów całkowania. Gdy bufor jest pełny, każdy nowy punkt nadpisuje najstarszy, więc zużycie pamięci jest stałe niezależnie od tego, jak długo trwa symulacja. Geometria Three.js jest odbudowywana co klatkę poprzez odczyt bufora w porządku chronologicznym i przesłanie wynikowych uporządkowanych tablic zmiennoprzecinkowych do GPU za pomocą BufferAttribute. Takie podejście zapobiega nieograniczonemu wzrostowi pamięci, jednocześnie utrzymując dokładnie żądaną długość śladu.

Tak. Każda orbita żyje w trójwymiarowej przestrzeni fazowej (x, y, z). Sam atraktor jest nieskończenie cienkim zbiorem o wymiarze fraktalnym ≈ 2,06 — jest zasadniczo złożoną dwuwymiarową powierzchnią, sfałdowaną nieskończenie wiele razy jak ciasto francuskie. W dwóch wymiarach widziałbyś tylko projekcję: klasyczny widok z boku pokazuje dwie połączone spirale, widok z góry pokazuje skrzydła motyla. Ta symulacja pozwala swobodnie orbitować, aby docenić prawdziwą geometrię 3D, w tym spiralne helisy każdego płata i wąski mostek je łączący.

Efekt motyla jest często mylnie rozumiany jako oznaczający, że motyl dosłownie powoduje tornado. To, co faktycznie stwierdza, to że w układzie chaotycznym mała zmiana warunków początkowych prowadzi do wykładniczo rozbieżnych wyników w czasie — czyniąc długoterminową prognozę fundamentalnie ograniczoną, a nie że motyl jest rzeczywistą fizyczną przyczyną. Drugim powszechnym błędnym przekonaniem jest to, że chaos oznacza losowość. Układ Lorenza jest całkowicie deterministyczny: przy identycznych warunkach startowych trajektoria jest za każdym razem identyczna. Pozorna losowość wynika z wykładniczej czułości na drobną niedokładność nieodłączną każdemu rzeczywistemu pomiarowi.