Elliptical Orbits: A Newtonian Simulation

Od obserwacji do prawa, następnie z prawa do siły

Johannes Kepler spędzony lata na dopasowywanie obserwacji Tycho Brahe'a Marsa, wykonanych gołym okiem, i wyprowadzanie trzech empirycznych praw ruchu planet wczesnych lat 1600., dekadę przed tym, kiedy ktokolwiek wiedział dlaczego były prawdziwe. Potrzebna była do tego zasada grawitacji uniwersalnej Newtona, F = -GMm/r², skierowana po linii łączącej dwa ciała, aby pokazać, że wszystkie trzy prawa Keplera nie są niezależnymi faktami o układzie słonecznym, lecz pojedynczym matematycznym konsekwencją jednej zasady siłowej zastosowanej do problemu dwóch ciał.

Pierwsza zasada: elipsy, a nie koła

Rozwiązanie równań Newtona dla zwartego obiegu dwóch ciał pod działaniem siły odwrotnej-kwadratowej daje cięcie skończone, w którym masa centralna znajduje się w jednym ogniskowym punkcie, a nie w centrum. Dla orbit związanych (energia całkowita ujemna) to cięcie jest zawsze elipsą; koło stanowi jedynie specjalny, dokładny przypadek z liczbą e oszczędną równą 0, wymagający precyzyjnego związku między prędkością początkową a promieniem. Zmiana prędkości początkowej nawet nieznacznie powoduje rozciągnięcie orbity w elipsę, w której planeta jest najbliżej gwiazdy na periście i najdalej na apozie – dokładnie taki kształt Kepler musiał wywnioskować z nieubelżywie nieregularnego ruchu Marsa.

Drugie prawo: równe obszary, konsekwencja momentu pędu

Drugie prawo Keplera mówi, że linia z gwiazdy na planetę wycina w równych odstępach czasu równe obszary – planeta porusza się najszybciej w perihelium i najwolniej w aphelium. Ponieważ siła grawitacji zawsze skierowana jest wzdłuż tej samej linii ("siła centralna"), ona wywiera zerowy moment bezwładności względem gwiazdy, więc zachodzi zachowanie momentu pędu L = m·r×v w trakcie orbity. Szybkość, z jaką jest wycinany obszar, wynosi dokładnie dA/dt = L/(2m), czyli stała – dlatego równe obszary w równych odstępach czasu to nic innego niż zachowanie momentu pędu, wyrażone w języku geometrycznym.

Trzecie prawo: okres kwadratowy, pół-oś obiegowa do potęgi trzeciej

T² = (4π² / G·M) · a³ T = okres orbitalny a = pół-oś obiegowa (średnia z najbliższego i najdalszego odległości) M = masa ciała centralnego (właściwa masa planety zanika, gdy jest znacznie mniejsza niż M) Trzecie prawo Keplera wypływa z równowagi siły grawitacji z wymaganiem centryfugacyjnym dla zamkniętego obiegów i integracji w jednym pełnym okresie. Jest to prawo, którego najczęściej używają astronomowie: mierzą okres planety (lub księżyca, lub gwiazdy podwójnej) oraz pół-oś obiegową, a masa centralna M wypada bezpośrednio - dosłownie tak ustalono masy Słońca, planet z księżycami i niezliczonych gwiazd karłowatych).

T² = (4π² / G·M) · a³

T  = orbital period
a  = semi-major axis (average of closest and farthest distance)
M  = mass of the central body (the planet's own mass drops out
     when it is much smaller than M)

Symulacja: dlaczego integrator ma znaczenie

Równania ruchu same są proste – przyspieszenie jest równe -GM r̂/r² - ale stabilna numeryczna ich integracja przez wiele orbit nie jest automatyczna. Niewłaściwa, jawna metoda Eulera czwartego rzędu (zwiększanie pozycji przy użyciu starej prędkości i aktualizacja prędkości) systematycznie wprowadza energię do systemu w każdym kroku; symulowany eliptyczny obieg w ten sposób widocznie spiralizuje się na zewnątrz, a efekt pogarsza się wraz z upływem czasu. Metoda Runge-Kutta czwartego rzędu (RK4) rozwiązuje to dla skal czasów istotnych tutaj poprzez próbkowanie przyspieszenia w czterech punktach w każdym kroku – na początku, dwukrotnie w połowie drogi z próbami półkrokowych i w punkcie docelowym – i łączy je w aktualizację, gdzie lokalny błąd skali z czasem krok piątym potęgą, co zapewnia znacznie lepsze zachowanie podczas dziesiątek orbit wymaganych do przekonującej symulacji.

Energia decydująca o kształcie

Jedno kolejne element rozwiązuje cały obraz: całkowita energia mechaniczna orbity, E = ½mv² - GMm/r, determinuje, jaką koniczną sekcję uzyskamy. E < 0 daje ograniczony eliptyczny orbitę (regime Keplera); E = 0 daje trajektorię paraboliczną, która ledwo ucieka w nieskończoność z zerową prędkością końcową; E > 0 daje hiperboliczny przelet, który ucieka z dużą prędkością. Dostosowywanie prędkości startowej w symulacji takiej jak ta polega głównie na ustawianiu E na wyższym lub niższym poziomie i obserwowanie, jak orbita zmienia się z ciasnej elipsy przez prawie okrąg, szeroką eliptykę i ostatecznie rozszerza się do paraboli i hiperboli – to najszybszy sposób na zobaczenie wszystkich trzech praw Keplera oraz ich ograniczeń w jednym ciągłym suwakie.

Często zadawane pytania

Wypróbuj na żywo

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