🌠 Symulacja grawitacyjna N ciał — orbity, zderzenia i asysty grawitacyjne
Interaktywna symulacja grawitacyjna N ciał z całkowaniem numerycznym RK4. Dodawaj własne ciała, wybieraj gotowe konfiguracje (Układ Słoneczny, gwiazdy podwójne, zderzenie galaktyk) i obserwuj rezonanse orbitalne oraz manewry asysty grawitacyjnej.
Podobne symulacje
O symulacji
Ta symulacja rozwiązuje klasyczny problem N ciał: każda para punktowych mas przyciąga się wzajemnie zgodnie z prawem Newtona F = G·m₁·m₂/r², a przyspieszenia parami są sumowane po wszystkich ciałach w każdej klatce. Równania ruchu są całkowane metodą czwartego rzędu Rungego–Kutty (RK4), a nie prostym krokiem Eulera, co daje znacznie lepszą dokładność przy tym samym kroku czasowym. Do wyrazu odległości dodawana jest długość zmiękczająca ε, dzięki czemu siła nigdy nie eksploduje, gdy dwa ciała zbliżą się bardzo blisko siebie. Energia kinetyczna i potencjalna są śledzone na bieżąco jako przybliżony sprawdzian jakości całkowania – dla trzech lub więcej ciał nie istnieje ogólne rozwiązanie w postaci zamkniętej, więc drobne zmiany położeń lub prędkości początkowych mogą skierować układ na zupełnie inne, chaotyczne tory.
🔬 Co pokazuje
Cztery gotowe konfiguracje: pięcioobiektowy układ Słoneczny, parę podwójną okrążaną przez lekkie trzecie ciało, dokładną trzyobiektową orbitę Ósemkową (Figure-8) oraz losową chmurę ciał. Ślady ujawniają kształty orbit, rezonanse, bliskie spotkania i odchylenia typu manewru procy grawitacyjnej, gdy integrator RK4 aktualizuje położenie i prędkość każdego ciała w każdej klatce.
🎮 Jak korzystać
Kliknij na płótnie, aby dodać nowe ciało o masie ustawionej suwakiem Masa nowego ciała. Dostosuj G, aby wzmocnić lub osłabić grawitację, dt, aby zamienić szybkość całkowania na dokładność, Długość śladu, aby zobaczyć więcej lub mniej historii orbity, oraz Zmiękczenie ε, aby kontrolować, jak łagodnie wygładzane są bliskie spotkania. Wstrzymaj, aby zatrzymać scenę, lub Wyczyść, aby zacząć od pustego układu.
💡 Czy wiesz, że…
Trzyobiektowa orbita ósemkowa – w której trzy równe masy nieustannie gonią się nawzajem po pojedynczym torze w kształcie lemniskaty – została odkryta numerycznie dopiero w 1993 roku, a jej istnienie rygorystycznie udowodnili Chenciner i Montgomery w 2000 roku. Wymaga dokładnego stosunku mas i konfiguracji prędkości; wersja użyta tutaj jest jedynie przybliżeniem, więc obserwuj, jak szybko rozpada się pod wpływem całkowania RK4 i błędów zaokrągleń.
Najczęściej zadawane pytania
Dlaczego problem N ciał nie ma ogólnego rozwiązania?
Dla dwóch ciał prawo grawitacji Newtona daje dokładne orbity w postaci zamkniętej (elipsy Keplera). Gdy dodane zostaje trzecie ciało, układ sprzężonych równań różniczkowych z reguły nie ma rozwiązania algebraicznego – Poincaré udowodnił w latach 90. XIX wieku, że problem trzech ciał jest ogólnie niecałkowalny. Ten symulator omija to, całkując równania ruchu numerycznie krok po kroku, zamiast rozwiązywać je analitycznie.
Dlaczego stosowana jest metoda RK4 zamiast prostego kroku Eulera?
Metoda Eulera aktualizuje położenie i prędkość, wykorzystując wyłącznie przyspieszenie z początku każdego kroku, co szybko kumuluje błąd i „wycieka" energię. Metoda Rungego–Kutty czwartego rzędu (RK4) próbkuje natomiast przyspieszenie w czterech punktach w obrębie kroku (etapy k1–k4) i łączy je w ważoną średnią, co daje błąd malejący znacznie szybciej wraz ze zmniejszaniem kroku czasowego dt – dużo lepszą długoterminową wierność orbit przy tym samym koszcie obliczeniowym.
Co właściwie robi parametr zmiękczający ε?
Siła grawitacji Newtona 1/r² dąży do nieskończoności, gdy dwa ciała zbliżają się do r = 0, co wysłałoby wartości liczbowe symulacji w nieskończoność podczas bliskiego spotkania. Zmiękczenie zastępuje r² wyrazem r²+ε² w formule przyspieszenia, ograniczając maksymalną siłę na małych odległościach. Większe ε sprawia, że niemal-zderzenia są łagodniejsze, ale mniej realistyczne fizycznie; mniejsze ε jest dokładniejsze, ale grozi numerycznymi eksplozjami przy dużej prędkości lub dużym dt.
Czym jest orbita ósemkowa i dlaczego ostatecznie się rozpada?
Orbita ósemkowa to okresowe rozwiązanie, w którym trzy równe masy gonią się nawzajem po pojedynczej krzywej w kształcie ósemki, odkryte numerycznie przez Cristophera Moore'a w 1993 roku i rygorystycznie udowodnione przez Chenciner i Montgomery w 2000 roku. Istnieje wyłącznie dla precyzyjnej kombinacji położeń początkowych, prędkości i równych mas. Ponieważ ten symulator wykorzystuje dyskretny krok czasowy i arytmetykę o skończonej precyzji, drobne odchylenia od tej dokładnej konfiguracji rosną wykładniczo – co jest cechą charakterystyczną chaosu – dlatego orbita widocznie rozpada się po wystarczającej liczbie obiegów.
Jak obliczana jest siła grawitacji między każdą parą ciał?
W każdej klatce symulator przechodzi przez każdą unikalną parę (i, j) i oblicza siłę Newtona F = G·mⁱ·m⁅/r² wzdłuż linii je łączącej, sumując wynikowe przyspieszenie aⁱ = Σ⁅ G·m⁅·(r⁅−rⁱ)/(|r⁅−rⁱ|²+ε²)^(3/2) na każdym ciele. Jest to bezpośrednia, dokładna suma parami (obliczenie rzędu O(n²)), a nie przybliżenie w rodzaju Barnes–Hut, więc jest dokładna, ale praktyczna jedynie dla skromnej liczby ciał używanej w tej symulacji.