PHANTOM's fast multipole gravity solver violates total linear momentum conservation by computing leaf-cell-to-supercell forces without reverse pairs; in neutron star simulations this shifts the center of mass by up to about one stellar radius.
Octo-Tiger: A New, 3D Hydrodynamic Code for Stellar Mergers that uses HPX Parallelisation
1 Pith paper cite this work. Polarity classification is still indexing.
abstract
OCTO-TIGER is an astrophysics code to simulate the evolution of self-gravitating and rotat-ing systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. OCTO-TIGER is currently optimised to simulate the merger of well-resolved stars that can be approximated by barotropic structures, such as white dwarfs or main sequence stars. The gravity solver conserves angular momentum to machine precision, thanks to a correction algorithm. This code uses HPX parallelization, allowing the overlap of work and communication and leading to excellent scaling properties, allowing for the computation of large problems in reasonable wall-clock times. In this paper, we investigate the code performance and precision by running benchmarking tests. These include simple problems, such as the Sod shock tube, as well as sophisticated, full, white-dwarf binary simulations. Results are compared to analytic solutions, when known, and to other grid based codes such as FLASH. We also compute the interaction between two white dwarfs from the early mass transfer through to the merger and compare with past simulations of similar systems. We measure OCTO-TIGERs scaling properties up to a core count of 80,000, showing excellent performance for large problems. Finally, we outline the current and planned areas of development aimed at tackling a number of physical phenomena connected to observations of transients.
fields
astro-ph.SR 1years
2025 1verdicts
CONDITIONAL 1representative citing papers
citing papers explorer
-
Non-conservation of linear momentum in widely used hierarchical methods in gravitational gas dynamics
PHANTOM's fast multipole gravity solver violates total linear momentum conservation by computing leaf-cell-to-supercell forces without reverse pairs; in neutron star simulations this shifts the center of mass by up to about one stellar radius.