Phurbas: An Adaptive, Lagrangian, Meshless, Magnetohydrodynamics Code. II. Implementation and Tests
Abstract
We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. The particles move with the fluid, so the time step is not limited by the Eulerian Courant-Friedrichs-Lewy condition. Full spatial adaptivity is required to ensure the particles fill the computational volume and gives the algorithm substantial flexibility and power. A target resolution is specified for each point in space, with particles being added and deleted as needed to meet this target. We have parallelized the code by adapting the framework provided by GADGET-2. A set of standard test problems, including 10-6 amplitude linear magnetohydrodynamics waves, magnetized shock tubes, and Kelvin-Helmholtz instabilities is presented. Finally, we demonstrate good agreement with analytic predictions of linear growth rates for magnetorotational instability in a cylindrical geometry. This paper documents the Phurbas algorithm as implemented in Phurbas version 1.1.
- Publication:
-
The Astrophysical Journal Supplement Series
- Pub Date:
- May 2012
- DOI:
- 10.1088/0067-0049/200/1/7
- arXiv:
- arXiv:1110.0836
- Bibcode:
- 2012ApJS..200....7M
- Keywords:
-
- hydrodynamics;
- magnetohydrodynamics: MHD;
- methods: numerical;
- Astrophysics - Instrumentation and Methods for Astrophysics;
- Astrophysics - Earth and Planetary Astrophysics
- E-Print:
- 14 pages, 14 figures, ApJS accepted, revised in accordance with changes to paper I (arXiv:1110.0835)