Тема: Сонячна система. Власна пісочниця. (Прочитано 624 раз)

zhyvchyk · « : 30 Січня 2016, 15:47:35 »

Просто для власного розвитку написав на С++ симулятор сонячної системи по мануалам системи HORIZONS.
Цікаво, хто небудь іще таким займався?
Ось зараз думаю, над цікавими задачами, де б його можна було застосувати. Може у кого будуть ідеї?
Варіанти "проінтегрувати комету" чи "подивитися як літають троянці і греки" давно вже проходили. Апофіс і всілякі інші тіла, які зближаються з планетами Сонячної системи, також прораховані до дірок.

Orion

Тоже как курсовой проект сделал на Delphi ПО для решения задачи N тел, взаимодействующих гравитационно. Составлял модели солнечной системы, тел падающих под углом к горизонту, двойная звезда, орбиты комет, перехват малого тела в системе двух масивных тел, симметричные системы. Использовал метод Рунге-Кутта. Редактор графический, для каждого объекта можно установить массу, размер, координаты и скорость в двумерном пространстве.

Orion

Еще фрагмент из презентации, сделанной на конкурс.

zhyvchyk

Orion, ну у мене з програмуванням GUI завжди було туго, тож все своє візуалізую в Origin або IDL.
Мій інтегратор консольний. Вивід даних безпосередньо в файли.
Інтегратор побудований по схемі предиктор-коректор, предиктором виступають явні схеми адамса-башфорта, коректором - неявні схеми адамса-мултона. Порядок інтегратора сталий і дорівнює 11, крок змінний.
Окрім "звичайного ньютона", при розрахуванні взаємодії між тілами, враховуються релятивістські поправки а також несферичності Сонця, Землі та Місяця.
Загалом кістяк інтегратора складають 11 масивних тіл: Сонце, планети, Місяць та Плутон і 343 астероїди головного поясу.
А в їх гравітаційному полі уже можна крутити які завгодно тестові безмасові частинки. Якщо ця тестова частинка є кометою, то також враховується вплав негравітаційних додатків.
Поділюся і своїми результатами.
Ось знайшов графіки датовані на комп*ютері початком 2014 року. На них по осі ОХ - час інтергування, по осі OY - різниця між модулем радіус-вектора i-того тіла взятого з системи horizons та модулем радіус-вектора того ж тіла порахованого моїм софтом. Корочше кажучи delta_r = |r_jpl| - |r_my|. Різниця виражена в кілометрах.