Zum Hauptinhalt springen
Autor/-in:

Külling Michaela

Programmierung einer Simulation des n-Körperproblems und Analyse mit dem Dreikörperproblem und den Planetenorbits

Begleitarbeit zur Maturitätsarbeit

Betreuer/-in:
José Abreu Castiñeira
Schule:

Kantonsschule Zürcher Unterland

Physik
Durch meine Maturitätsarbeit erlangte ich ein tieferes Verständnis im Programmieren und in der Astrophysik. Es hat mir Freude bereitet, mein Wissen zu erweitern und es anzuwenden.
Abstract

Das Ziel meiner Maturitätsarbeit war, eine Simulation für das n-Körperproblem in Python zu programmieren, welche die Bewegung der Körper richtig berechnet. Zusätzlich sollte die Simulation die Grundsätze der Energie- und Impulserhaltung berücksichtigen. Die Berechnungen für die resultierenden Kräfte auf die beteiligten Körper wurden mit dem Newtonschen Gravitationsgesetz vorgenommen. Die Bewegung der Körper wurde mit dem numerischen Integrationsalgorithmus Leapfrog berechnet.

Die Simulation wurde zum einen anhand eines von Joseph-Louis Lagrange entdeckten Spezialfalls des Dreikörperproblems, zum anderen mit den Planetenorbits in unserem Sonnensystem analysiert.

Die Positionen und Geschwindigkeiten der Körper im Spezialfall des Dreikörperproblems konnten analytisch berechnet werden. Für die Planetenbewegung im Sonnensystem wurden die Positionen und Geschwindigkeiten den Ephemeriden des Jet Propulsion Laboratory der NASA entnommen.

Bei beiden simulierten Situationen blieben die Energie und der lineare Impuls im System erhalten. Die Bewegung erfolgte, wie erwartet, regelmässig und die Endpositionen stimmten relativ gut mit den erwarteten überein. Bei den Planetenorbits im Sonnensystem stellten sich bei den inneren Planeten grössere Abweichungen ein, weil für die Simulation die Zeitberechnung des Kalenders, nach welchem sich die Ephemeriden richten, nicht berücksichtigt wurde.