Entstehung des Mondes

[Bernhard]

Hast du eine Vermutung, warum der Integrator bei close encounters so langsam wird?

Bei Vorbeiflügen hat man stärkere Bahnkrümmungen als im Raum fern von gravitierenden Massen. Dies sollte ein guter Integrator in der Schrittweite berücksichtigen. Kleine Schrittweite bedeutet dann viel Rechenzeit und damit langsamer Programmablauf. Hast Du das mit dem SCATR-Code ausprobiert?

[SRMeister]

Die Trümmerteile des hypothetischen Urmondes sind allerdings gemäß Vorgabe so schnell, dass man für eine aussagekräftige Analyse des Vorbeifluges von Trümmerteilen an Planeten ziemlich viel Rechenzeit benötigt (>>1 Stunde).

Ich weis nicht genau ob ich dich hier richtig verstanden habe. Benötigte jeder einzelne Vorbeiflug eine Std. oder die Simulation insgesamt?

Nein ich habe das mit dem SCATR nicht ausprobiert. Ich habe momentan noch einen anderen Integrator gefunden, Swifter, basierend auf SWIFT, worauf ja auch SCATR basiert. Das Programm ist nahezu identisch strukturiert wie SCATR, bietet aber noch ein paar andere Integratoren

[Bernhard]

Ich weis nicht genau ob ich dich hier richtig verstanden habe. Benötigte jeder einzelne Vorbeiflug eine Std. oder die Simulation insgesamt?

Sämtliche Zeitangaben bezogen sich immer auf die gesamte Simulation.

[Bernhard]

Die Frage wäre jetzt, wie man daraus möglichst elegant (also mit möglichst wenigen Rechenschritten) die Bahnparameter bestimmt, also numerische Exzentrizität, große Halbachse usw.

Ich habe mittlerweile in der Literatur zum Zweikörperproblem die dazu notwendigen Gleichungen gefunden. Ich hatte den Eindruck, dass auch im SCATR-Code so eine Rechnung gemacht wird und man könnte damit dann lediglich die gravitativen Störungen der Planeten untereinander numerisch berechnen. Die grobe Bahn als Keplerbahn ist im Gegensatz dazu weitgehend analytisch berechenbar. Ich habe momentan nur wirklich wenig Zeit das auch zu programmieren.
Gruß

EDIT: Die Fehlermeldung "Danby konvergiert nicht" könnte demnach bedeuten, dass an einer bestimmten Stelle die Keplerbahn nicht ausreichend genau berechnet werden kann. Bei der Keplerbahn wird zuletzt bekanntlich die wahre Anomalie des Planeten berechnet. Das geht nur numerisch und mit einer vorgegebenen Genauigkeit. Ich vermute deswegen, dass diese Genauigkeit über die DANBY-Parameter vorgegeben wird.

[Bernhard]

Hallo zusammen,

BdW hat diesen Monat eben dieses Thema auf seiner Titelseite stehen: http://www.wissenschaft.de/wissensch...lles_heft.html

Leider gibt es scheinbar nur das gesamte Heft in digitaler Form.
Gruß