High-order regularised symplectic integrator for collisional planetary systems
| Autor: | Jacques Laskar, Gwenaël Boué, Mickael Gastineau, Antoine Petit |
|---|---|
| Přispěvatelé: | Hôpital Avicenne [AP-HP], Assistance publique - Hôpitaux de Paris (AP-HP) (AP-HP), Institut de Mécanique Céleste et de Calcul des Ephémérides (IMCCE), Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire de Paris, Université Paris sciences et lettres (PSL)-Université Paris sciences et lettres (PSL)-Université de Lille-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS) |
| Jazyk: | angličtina |
| Rok vydání: | 2019 |
| Předmět: |
Solar System
planets and satellites: dynamical evolution and stability FOS: Physical sciences Astrophysics 01 natural sciences Stability (probability) methods: numerical 0103 physical sciences Statistical physics Symplectic integrator 010303 astronomy & astrophysics Earth and Planetary Astrophysics (astro-ph.EP) Physics [PHYS]Physics [physics] 010308 nuclear & particles physics Astronomy and Astrophysics Radius celestial mechanics Planetary system Celestial mechanics Space and Planetary Science Integrator Astrophysics::Earth and Planetary Astrophysics [PHYS.ASTR]Physics [physics]/Astrophysics [astro-ph] Astrophysics - Earth and Planetary Astrophysics Symplectic geometry |
| Zdroj: | Astronomy and Astrophysics-A&A Astronomy and Astrophysics-A&A, 2019, 628, pp.A32. ⟨10.1051/0004-6361/201935786⟩ Astronomy and Astrophysics-A&A, EDP Sciences, 2019, 628, pp.A32. ⟨10.1051/0004-6361/201935786⟩ |
| ISSN: | 0004-6361 |
| DOI: | 10.1051/0004-6361/201935786⟩ |
| Popis: | We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators, while also reducing the effective step size whenever two planets encounter. We use a high order MVS scheme such that it is possible to integrate with large time steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e with a mutual separation of a fraction of the physical radius) while using the same time-step as in weakly perturbed problem such as the Solar System. We demonstrate the long term behaviour on systems of six super-Earths experiencing strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost we obtain much better energy conservation. Comment: Accepted in Astronomy and Astrophysics, updated to match the published version |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|