Relativistic magnetohydrodynamics in dynamical spacetimes: A new AMR implementation
arXiv:1007.2848 · doi:10.1103/PhysRevD.82.084031
Abstract
We have written and tested a new general relativistic magnetohydrodynamics (GRMHD) code, capable of evolving MHD fluids in dynamical spacetimes with adaptive-mesh refinement (AMR). Our code solves the Einstein-Maxwell-MHD system of coupled equations in full 3+1 dimensions, evolving the metric via the Baumgarte-Shapiro Shibata-Nakamura (BSSN) formalism and the MHD and magnetic induction equations via a conservative, high-resolution shock-capturing scheme. The induction equations are recast as an evolution equation for the magnetic vector potential, which exists on a grid that is staggered with respect to the hydrodynamic and metric variables. The divergenceless constraint div(B)=0 is enforced by the curl of the vector potential. Our MHD scheme is fully compatible with AMR, so that fluids at AMR refinement boundaries maintain div(B)=0. In simulations with uniform grid spacing, our MHD scheme is numerically equivalent to a commonly used, staggered-mesh constrained-transport scheme. We present code validation test results, both in Minkowski and curved spacetimes. They include magnetized shocks, nonlinear Alfvén waves, cylindrical explosions, cylindrical rotating disks, magnetized Bondi tests, and the collapse of a magnetized rotating star. Some of the more stringent tests involve black holes. We find good agreement between analytic and numerical solutions in these tests, and achieve convergence at the expected order.
23 pages, 27 figures, submitted to Phys. Rev. D
References in corpus (21)
- Calibration of Moving Puncture Simulations
- Accurate evolutions of inspiralling neutron-star binaries: prompt and delayed collapse to black hole
- Binary black hole merger dynamics and waveforms
- How to move a black hole without excision: gauge conditions for the numerical evolution of a moving puncture
- Simulating coalescing compact binaries by a new code SACRA
- General relativistic simulations of magnetized binary neutron star mergers
- Fully General Relativistic Simulations of Black Hole-Neutron Star Mergers
- High-spin binary black hole mergers
- Evolution of magnetized, differentially rotating neutron stars: Simulations in full general relativity
- Relativistic Simulations of Black Hole-Neutron Star Mergers: Effects of black-hole spin
- Magnetorotational collapse of massive stellar cores to neutron stars: Simulations in full general relativity
- Geometry and Regularity of Moving Punctures
- Wormholes and trumpets: the Schwarzschild spacetime for the moving-puncture generation
- Relativistic MHD with Adaptive Mesh Refinement
- Filling the holes: Evolving excised binary black hole initial data with puncture techniques
- Modeling gravitational recoil from precessing highly-spinning unequal-mass black-hole binaries
- Excision without excision: the relativistic turducken
- Magnetorotational collapse of very massive stars to black holes in full general relativity
- Evolution of near-extremal-spin black holes using the moving puncture technique
- Collapse of magnetized hypermassive neutron stars in general relativity: Disk evolution and outflows
- Constraint preserving boundary conditions for the Ideal Newtonian MHD equations