Mechanika orbitalna: Keplar, Verlet i punkty Lagrange'a

Zasady Keplera i moment, w którym przestają być wystarczające

Dla jednej planety obiegującej pojedynczą, znacznie cięższą gwiazdę, trzy zasady Keplera opisują ruch dokładnie: orbity są elipsami z gwiazdą w jednym skupisku; linia łącząca gwiazdę z planety przeprowadza równe pola w równe czasy (więc planeta porusza się szybciej blisko perihelionu, wolniej blisko afelionu); a kwadrat okresu obiegu jest proporcjonalny do sześcianu półosi wielkiej. Problem dwóch ciał ma dokładne rozwiązanie w formie zamkniętej — można bezpośrednio obliczyć pozycję planety na dowolnym przyszłym czasie z sześciu elementów orbitalnych, bez potrzeby symulacji.

Chwilą, gdy do systemu wejdzie trzecie ciało — druga planeta, luna, statkowiec kosmiczny — system ogólnie nie ma rozwiązania w formie zamkniętej. Jest to słynny problem trzech ciał: Poincaré pokazał w latach 1890–pierwszych, że ścieżki ruchu mogą być chaotyczne i nie można je wyrazić ogólnie za pomocą skończonej kombinacji funkcji elementarnych. Praktycznym odpowiedzią jest integracja numeryczna: przesunięcie pozycji i prędkości każdego ciała o małe kroki czasowe, korzystając z siły grawitacyjnej, jaką każde ciało obecnie wykazuje na każdym innym, a następnie powtórzenie.

Dlaczego integrator musi być metoda Verleta szybkiej

Ta symulacja przyspiesza każdą sytuację użyciem metody Verleta szybkiej, a nie bardziej zaawansowanego metodu o wyższym rzędzie, ponieważ mechaniczna orbitalna dba o właściwość, która większość integratorów ignoruje: zachowanie energii w długowiecznej operacji. Metoda Verlet jest symplektyczna — zawsze zachowuje energetyczny cień bliski prawdziwej energii, dlatego orbita owalna pozostaje owalną i orbita okresowa pozostaje okresową przez miliony kroków. Metoda o wyższym rzędzie, ale nie symplektyczna jak np. RK4, jest bardziej precyzyjna na każdym kroku, ale powoli i regularnie utrata energii prowadzi do tego, że planeta integrowana przez RK4 po dostatecznie długim czasie widocznie spiraliwej w głąb lub zewnątrz, mimo że każdy indywidualny krok wygląda prawie dokuczliwie dokładny.

x += v * h + 0.5 * aOld * h * h;
const aNew = gravityFrom(allOtherBodies, x);   // sum of GM_i / r_i^2 terms
v += 0.5 * (aOld + aNew) * h;
aOld = aNew;

Punkty Lagrange'a i trojkatne asteroidy Trojańskie

W ograniczonym trójcieściennym problemie (gdzie masa jednego ciała jest niewielka w porównaniu do dwóch pozostałych, np. asterozyd orbitujący w systemie Słońce–Jowisz), istnieją pięć specjalnych pozycji, punkty Lagrange'a L1–L5, gdzie siła grawitacji obu dużych ciał i efekt centripetowy orbitalny równoważają się, co pozwala trzeciemu ciału, bezmasowemu, siedzieć w ustalonej pozycji względem innych dwóch, koordynującym się z nimi. Punkty L1–L3 leżą na linii przechodzącej przez oba duże ciała i są niestabilne punkty łukowate – przydatne dla statków kosmicznych (np. obserwatorium słonecznych w punkcie Słońce– Ziemia L1), które potrzebują często małych poprawek, nieprzydatne jednak dla ciał pozostawionych samych.

Punkty L4 i L5 są inny: leżą 60° przed i za mniejszym z dużych ciał na jego orbicie, tworząc równoległobok równoboczny z dwoma dużymi masami. Dla wystarczająco dużej proporcji masy między oboma głównymi ciałami są one dynamicznie stabilne. Asterozydy Trojańskie Jowisza są zatrzymywane dokładnie w ten sposób w punktach L4 i L5 – tysiące ich, niektóre o średnicy ponad 100 km, orbitują razem z Jowiszem Słońca bez potrzeby integracji każdego indywidualnie jako jakby były przypadkową płytką.

Zyski z grawitacji: kradniec chwilę prędkości od planety

Efekt przelotu (lub flyby) zmienia prędkość łazika w odniesieniu do Słońca bez wykorzystania żadnego paliwa, przez wymianę chwilności między łazikiem a planetą. W własnym układzie odniesienia planety spotkanie jest elastyczne — łazik opuszcza z taką samą prędkością, z jaką przybył, ale kierunek zmieniony o kąt deflekcji hiperbolicznej. Podczas przeliczenia na układ odniesienia Słońca, ta zmiana kierunku dodaje (lub odejmuje) wektor prędkości orbity planety do wektora prędkości łazika, co może spowodować znaczącą netto zmianę prędkości w odniesieniu do Słońca — mechanismum za tą prawdopodobnie każdą rzeczywistą ścieżką lotu każdego missyjnego zadania poza systemem słonecznym.

Często zadawane pytania

Wypróbuj na żywo

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