Fala 33: Tsunami, Konstelacja Satelitarna i Wulkan

Liczby platformy

Symulacje Fali 33

🌊 Tsunami

Równania płytkiej wody

Tsunami zasadniczo różnią się od fal oceanicznych generowanych przez wiatr. Ich długość fali — zazwyczaj 100–500 km — znacznie przekracza głębokość oceanu (średnio 3,8 km), więc kolumna wody oscyluje w zasadzie w fazie od powierzchni do dna morskiego. To jest reżim płytkiej wody, opisywany przez uśrednione po głębokości równania płytkiej wody (SWE):

∂η/∂t = −∇·(H·u)

∂u/∂t = −g∇η

Tu η oznacza wzniesienie swobodnej powierzchni, u to uśredniona po głębokości prędkość pozioma, H to niezakłócona głębokość wody, a g to przyspieszenie grawitacyjne. Połączenie tych dwóch równań daje liniowe równanie falowe ze zmienną przestrzennie prędkością fali c = √(gH): fale poruszają się szybciej na głębokiej wodzie i gwałtownie zwalniają w miarę zbliżania się do płytkich wybrzeży, powodując spiętrzanie i wzmocnienie.

Implementacja: schemat leapfrog na siatce przesuniętej

Symulacja wykorzystuje siatkę 220×220 komórek ze schematem krokowania czasowego leapfrog (Störmer-Verlet). Wysokość powierzchni η jest przechowywana w środkach komórek; składowe prędkości u i v są przechowywane odpowiednio na wschodnich i północnych ściankach komórek (przesunięta siatka Arakawy typu C). Równania aktualizacji to:

• Krok prędkości: u += −g·dt · (η[i+1] − η[i]) / dx
• Krok wysokości: η_next = damp · (2η − η_prev) + dt² · ∇·(c²∇η)

Zmienna prędkość fali c = √(gH) jest obliczana na każdej ściance komórki przez uśrednienie głębokości dwóch sąsiadujących komórek. Zapewnia to, że operator Laplace'a poprawnie obsługuje nieciągłości batymetrii. Absorbujące warunki brzegowe tłumią wychodzące fale na wszystkich czterech krawędziach, aby zapobiec zanieczyszczaniu symulacji odbiciami od granicy siatki.

Presety batymetrii

Cztery profile głębokości ilustrują różne środowiska, w których tsunami zachowują się odmiennie:

• ⛰️ Kontynent — szelf kontynentalny z łagodnym nachyleniem przechodzącym od głębokiego oceanu (~4000 m) do płytkich wód przybrzeżnych (~50 m). Kliknij w strefie głębokiej, aby wywołać trzęsienie ziemi w stylu strefy subdukcji.
• 🏝️ Wyspa — góra podwodna wznosząca się z głębokiego oceanu, modelująca łańcuch wysp wulkanicznych. Fale załamują się wokół wyspy i mogą ulegać wzmocnieniu po przeciwnej stronie.
• 🌊 Basen — zamknięty basen z płytkimi krawędziami, wytwarzający wielokrotne odbicia i rezonansowe wzory fal stojących podobne do oscylacji w porcie (sejsze).
• 🏔️ Rów oceaniczny — rów oceaniczny otoczony górami podwodnymi. Rów wychwytuje energię fali i prowadzi ją wzdłuż swojej długości — zjawisko znane jako uwięzienie fal krawędziowych.

Wyzwalanie trzęsienia ziemi

Kliknięcie lub przeciąganie po siatce dodaje gaussowskie wypiętrzenie dna morskiego o konfigurowalnej amplitudzie (1–15 m) i promieniu (4–30 komórek). Początkowe przemieszczenie jest dodawane bezpośrednio do pola η, symulując impulsywną deformację powierzchni morza spowodowaną przez trzęsienie ziemi typu uskoku odwróconego. Powstały pakiet falowy rozchodzi się następnie na zewnątrz z prędkością c = √(gH), widocznie zwalniając z prędkości głębokiego oceanu (~200 m/s przy głębokości 4000 m) do prędkości płytkiego wybrzeża (~22 m/s przy głębokości 50 m).

Co warto wypróbować

• W presecie Kontynent kliknij w głębokim oceanie i obserwuj, jak czoło fali przyspiesza nad głęboką równiną abisalną, a następnie gwałtownie zwalnia i stromieje po wejściu na szelf kontynentalny.
• W presecie Basen wywołaj wiele źródeł trzęsienia ziemi i obserwuj wzory interferencji — konstruktywne wzmocnienie i destruktywne wygaszanie.
• Obniż suwak absorpcji blisko 1,000, aby zezwolić na odbicia brzegowe i zobaczyć energię fali uwięzioną w domenie, budującą chaotyczną interferencję.
• Ustaw dużą wysokość wypiętrzenia (15 m), aby symulować zdarzenie typu megathrust, i obserwuj szczyt odczytu maksymalnego wzniesienia powierzchni w miarę spiętrzania się fali.

🛰️ Konstelacja Satelitarna

Konstelacje Walker Delta

Konstelacja Walker Delta jest opisywana przez trzy parametry: T/P/F — T to łączna liczba satelitów, P to liczba płaszczyzn orbitalnych, a F to współczynnik fazowania określający kątowe przesunięcie między płaszczyznami. W obrębie każdej płaszczyzny satelity są rozmieszczone równomiernie; płaszczyzny są równomiernie rozłożone pod względem rektascensji węzła wstępującego (RAAN) wokół równika przy tej samej inklinacji. Ta symetria maksymalizuje jednorodność pokrycia dla danej łącznej liczby satelitów.

Cztery presety obejmują główne architektury operacyjnych konstelacji:

• 🔵 Klasa Starlink — 6 płaszczyzn × 10 satelitów przy inklinacji 53°, wysokość 550 km. Pokrycie średnich szerokości geograficznych, krótki okres orbitalny (~95 min), niskie opóźnienie.
• 🟠 GPS — 6 płaszczyzn × 4 satelity przy inklinacji 55°, wysokość 20 200 km (MEO). Klasyczna konstelacja nawigacyjna; okres orbitalny ~12 h, globalne pokrycie poniżej 55° szerokości geograficznej.
• 🟡 Galileo — 3 płaszczyzny × 8 satelitów przy inklinacji 56°, wysokość 23 222 km. Europejska konstelacja GNSS na średniej orbicie okołoziemskiej.
• 🟢 Polarna — 8 płaszczyzn × 6 satelitów przy inklinacji 90°, wysokość 800 km. Orbity polarne w stylu heliosynchronicznym zapewniające pełne globalne pokrycie, w tym biegunów.

Renderowanie 3D w Three.js

Wizualizator wykorzystuje Three.js r160 ze wzorem pozycji orbity kepplerowskiej. Każdy satelita jest umieszczany w:

x = a(cosΩcosν − sinΩsinνcosi)
y = a·sinν·sini
z = a(sinΩcosν + cosΩsinνcosi)

gdzie a to półoś wielka (promień orbity), Ω to RAAN, ν to anomalia prawdziwa, a i to inklinacja. Początkowa anomalia prawdziwa każdego satelity jest przesunięta tak, aby satelity w obrębie płaszczyzny były rozmieszczone równomiernie, a współczynnik fazowania Walker F określa przesunięcie anomalii między płaszczyznami. Dla konstelacji Walker Delta RAAN dla p-tej płaszczyzny wynosi p·2π/P, a przesunięcie anomalii wynosi p·F·2π/(T).

Okręgi zasięgu są rysowane bezpośrednio na sferze Ziemi jako pierścienie EllipseCurve, nachylone tak, aby były zgodne z torem naziemnym satelity. Połowa kąta zasięgu ρ jest obliczana z:

ρ = arccos(R_e/a · cosε_min) − ε_min

gdzie ε_min to minimalny kąt elewacji (domyślnie 5°), a R_e = 6371 km. Procent pokrycia jest szacowany jako ułamek sfery Ziemi pokryty przez sumę śladów satelitów.

Co warto wypróbować

• Przełącz ze Starlink na GPS i obserwuj, jak okręgi zasięgu stają się znacznie większe przy wyższej wysokości, osiągając globalne pokrycie przy znacznie mniejszej liczbie satelitów.
• Ustaw inklinację na 90° (Polarna) i obserwuj, jak płaszczyzny orbitalne blokują się przez bieguny; przełącz nakładkę zasięgu, aby potwierdzić pełne globalne pokrycie.
• Zwiększ prędkość czasu do 1000×, aby obserwować, jak konstelacja wykonuje pełne orbity w ciągu kilku sekund; zauważ, że fazowanie Walker zapewnia, że żadne dwa satelity nigdy nie znajdują się w tej samej części swojej orbity.
• Przeciągnij, aby obracać, przewiń, aby przybliżyć — 3D Ziemia pozwala sprawdzić zasięg na dowolnej szerokości geograficznej, w tym czapy polarne, gdzie konstelacje o nachyleniu mają luki.

🌋 Erupcja Wulkanu

Style erupcji wulkanicznych

Erupcje wulkaniczne obejmują ogromny zakres intensywności i stylów, determinowany głównie przez lepkość magmy i zawartość rozpuszczonych substancji lotnych (gazów). Nisko lepka magma bazaltowa umożliwia łatwą ucieczkę gazu, wytwarzając łagodne erupcje efuzywne i słabo eksplozywne. Wysoko lepka magma ryolitowa uwięzia gaz, aż ciśnienie wywoła katastrofalną fragmentację. Wskaźnik eksplozywności wulkanicznej (VEI) kwantyfikuje to na skali logarytmicznej; symulacja modeluje cztery klasyczne style:

• 🌋 Strombolijski — łagodne, rytmiczne wybuchy wyrzucające rozżarzone bomby lawowe i rozpryski przy niskiej prędkości. Nazwany od wulkanu Stromboli we Włoszech, który wybucha niemal nieprzerwanie od 2000 lat.
• 🔥 Hawajski — trwałe fontanny lawy zasilające potoki lawy. Niska lepkość bazaltu umożliwia spokojną efuzję przerywaną dramatycznymi fontannami ognia sięgającymi do 500 m wysokości.
• 💥 Wulkański — odrębne, armatnie eksplozje wyrzucające gęste bloki i popiół. Korek z zastygłej lawy lub lepkiej magmy jest okresowo wystrzeliwany, gdy ciśnienie narasta poniżej.
• ☁️ Pliniański — megawybuchy podtrzymywane przez intensywną fragmentację magmy. Kolumna erupcyjna wznosi się dziesiątki kilometrów w stratosferę, wtryskując aerozole siarkowe, które mogą ochładzać globalny klimat. Nazwany od Pliniusza Starszego, który zginął, obserwując erupcję Wezuwiusza w 79 r. n.e.

System cząsteczkowy

Symulacja utrzymuje pulę do 3000 cząsteczek, z których każda należy do jednego z czterech typów fizycznych:

• Bomba lawowa (pomarańczowo-czerwona) — duża masa, tor balistyczny, przekształca się w przepływ piroklastyczny przy uderzeniu o ziemię. Podlega grawitacji i znoszeniu przez wiatr.
• Popiół (szary) — unosi się na dużej wysokości, powoli opada, szeroko się rozprzestrzenia. Dryfuje z symulowanym wiatrem i stopniowo opada w postaci deszczu.
• Chmura gazu (żółto-zielona) — szybko się wznosi, blaknie podczas rozpraszania. Modeluje pióropusze SO&sub2; i pary wodnej.
• Przepływ piroklastyczny (ciemnopomarańczowy) — gęsta gorąca mieszanina gazu i cząstek, która przylega do ziemi i przyspiesza w dół zbocza. Najbardziej śmiertelne zagrożenie wulkaniczne; modelowane jako niskokątowy, szybki strumień cząstek przylegający do ziemi.

Prędkości pionowego wystrzelenia są pobierane z rozkładu zależnego od presetu i mnożone przez suwak Siły Erupcyjnej ustawiany przez użytkownika. Słabe poziome pole wiatru (0,008 px/klatkę) znosi wszystkie cząsteczki w jedną stronę, wytwarzając charakterystyczny asymetryczny pióropusz popiołu widoczny na zdjęciach z podwiatru.

Narastanie ciśnienia magmy

Między erupcjami ciśnienie magmy narasta w tempie ustawionym przez suwak Tempo Ciśnienia. Gdy ciśnienie osiąga 100%, wulkan automatycznie wybucha — wyrzucając wybuch cząsteczek, których liczba i rozkład typów zależą od presetu stylu erupcji. Ten automatyczny mechanizm modeluje cykl regeneracji i uwalniania w zamkniętym kanale wulkanicznym. Podczas aktywnej erupcji ciśnienie szybko spada; pozostaje na zerze, dopóki cząsteczki są jeszcze w powietrzu, a następnie zaczyna się odbudowywać, gdy komin się oczyszcza.

Co warto wypróbować

• Wybierz Pliniański i ustaw Tempo Ciśnienia na maksimum: wulkan wybucha często z masywnymi kolumnami popiołu. Obniż Tempo Ciśnienia i obserwuj rzadsze, ale bardziej intensywne pojedyncze erupcje.
• Przełącz na Strombolijski i ustaw wysoką Lepkość: erupcje stają się bardzo krótkie, wyrzucając za każdym razem jedynie garść bomb lawowych — dokładnie jak prawdziwy Stromboli.
• Naciśnij przycisk Erupcja ręcznie w dowolnym trybie, aby wywołać erupcję na żądanie niezależnie od aktualnego poziomu ciśnienia.
• Obserwuj przepływy piroklastyczne w presetach Wulkański i Pliniański — gęste strumienie cząsteczek przylegają do terenu i przyspieszają w dół zbocza, zawsze skoncentrowane po bardziej stromej stronie stożka.

Uwagi techniczne

• Absorpcja brzegowa tsunami — warstwa absorbująca wykorzystuje funkcję gąbczastą: komórki w odległości 10 komórek od każdej krawędzi siatki stosują współczynnik tłumienia damp < 1 na krok. Współczynnik jest obliczany jako 1 − 0,01 × (1 − d/10)² gdzie d to odległość od krawędzi, dając kwadratowe zwężenie od pełnej absorpcji na granicy do braku absorpcji przy d = 10.
• Fazowanie Walker Delta — RAAN p-tej płaszczyzny (indeksowanej od 0) wynosi p × 2π / P. Początkowa anomalia prawdziwa s-tego satelity w płaszczyźnie p wynosi (s / S) × 2π + (p × F / T) × 2π, gdzie F to parametr fazowania Walker (0 w presecie klasy Starlink, 1 dla GPS, 0 dla Galileo, 1 dla Polarnej). Zapewnia to równomierne rozmieszczenie satelitów zarówno w obrębie płaszczyzny, jak i między płaszczyznami.
• Pula cząsteczek wulkanu — martwe cząsteczki (ptype = 0) są przetwarzane na początek aktywnej listy cząsteczek, aby uniknąć alokacji pamięci podczas wybuchów erupcji. Przy 3000 aktywnych cząsteczek pula jest pełna; nowe cząsteczki erupcji nadpisują najstarsze aktywne, zapobiegając nieograniczonemu zużyciu pamięci.
• Wszystkie symulacje Fali 33 mają pełne strony EN + UK i są zarejestrowane w simulations.json.

Tagi

Tsunami Równania płytkiej wody Propagacja fal Batymetria Trzęsienie ziemi Konstelacja Satelitarna Walker Delta Orbity kepplerowskie Three.js GPS Starlink Wulkan Erupcja Przepływ piroklastyczny Symulacja cząsteczkowa Fala 33

Zapowiedź Fali 34

Trzy symulacje planowane na Falę 34:

• Cykl skalny — animowany diagram geologicznego cyklu skalnego: wypiętrzanie, erozja, sedymentacja, pogrzebanie, metamorfizm i topienie magmowe wraz z interaktywnymi kontrolkami skali czasu.
• Kominy hydrotermalne — symulacja głębokomorskiego czarnego dymiarza z ekosystemem chemosyntetycznym: wzrost kominów mineralnych, konwekcja przegrzanego płynu i fauna kominowa reagująca na gradienty temperatury.
• Płyn MPM — symulacja płynu metodą punktów materiałowych: hybrydowy solwer cząsteczka-siatka obsługujący duże deformacje, swobodne powierzchnie i interakcję płyn-ciało stałe w jednej zunifikowanej strukturze.

Wszystkie symulacje Fali 34 zostaną wydane ze stronami EN + UK w dniu premiery.