Switching integrators reversibly in the astrophysical N-body problem
We present a simple algorithm to switch between N-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader. Upgrading an ordinary non-reversible switching integrator to...
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Article (Journal) |
| Language: | English |
| Published: |
July 2023
|
| In: |
Monthly notices of the Royal Astronomical Society
Year: 2023, Volume: 522, Issue: 3, Pages: 4639-4648 |
| ISSN: | 1365-2966 |
| DOI: | 10.1093/mnras/stad657 |
| Online Access: | Verlag, lizenzpflichtig, Volltext: https://doi.org/10.1093/mnras/stad657
|
| Author Notes: | David M. Hernandez and Walter Dehnen |
| Summary: | We present a simple algorithm to switch between N-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader. Upgrading an ordinary non-reversible switching integrator to a reversible one is straightforward and introduces no appreciable computational burden in our tests. Our method checks whether the integrator during the time-step violates a time-symmetric selection condition and redoes the step if necessary. In our experiments, a few per cent of steps would have violated the condition without our corrections. By eliminating them, the algorithm avoids long-term error accumulation, of several orders of magnitude in some cases. |
|---|---|
| Item Description: | Veröffentlicht: 2. März 2023 Gesehen am 21.07.2023 |
| Physical Description: | Online Resource |
| ISSN: | 1365-2966 |
| DOI: | 10.1093/mnras/stad657 |