Export citation

Export citation

Choose format for download:

Download Citation
  • Access by Xinjiang University

Relativistic magnetohydrodynamics in dynamical spacetimes: Improved electromagnetic gauge condition for adaptive mesh refinement grids

Zachariah B. Etienne*, Vasileios Paschalidis, Yuk Tung Liu, and Stuart L. Shapiro

  • Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA

  • *zetienne@illinois.edu
  • Also at Department of Astronomy and NCSA, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.

Phys. Rev. D 85, 024013 – Published 11 January, 2012

DOI: https://doi.org/10.1103/PhysRevD.85.024013

Abstract

We recently developed a new general relativistic magnetohydrodynamic code with adaptive mesh refinement that evolves the electromagnetic (EM) vector potential Ai instead of the magnetic fields directly. Evolving Ai enables one to use any interpolation scheme on refinement level boundaries and still guarantee that the magnetic field remains divergenceless. As in classical EM, a gauge choice must be made when evolving Ai, and we chose a straightforward “algebraic” gauge condition to simplify the Ai evolution equation. However, magnetized black hole–neutron star (BHNS) simulations in this gauge exhibit unphysical behavior, including the spurious appearance of strong magnetic fields on refinement level boundaries. This spurious behavior is exacerbated when matter crosses refinement boundaries during tidal disruption of the neutron star. Applying Kreiss-Oliger dissipation to the evolution of the magnetic vector potential Ai slightly weakens this spurious magnetic effect, but with undesired consequences. We demonstrate via an eigenvalue analysis and a numerical study that zero-speed modes in the algebraic gauge, coupled with the frequency filtering that occurs on refinement level boundaries, are responsible for the creation of spurious magnetic fields. We show that the EM Lorenz gauge exhibits no zero-speed modes, and as a consequence, spurious magnetic effects are quickly propagated away, allowing for long-term, stable magnetized BHNS evolutions. Our study demonstrates how the EM gauge degree of freedom can be chosen to one’s advantage, and that for magnetized BHNS simulations the Lorenz gauge constitutes a major improvement over the algebraic gauge.

Article Text

References (54)

  1. M. D. Duez, S. L. Shapiro, and H.-J. Yo, Phys. Rev. D 69, 104016 (2004).
  2. Z. B. Etienne, Y. T. Liu, and S. L. Shapiro, Phys. Rev. D 82, 084031 (2010).
  3. M. Shibata and T. Nakamura, Phys. Rev. D 52, 5428 (1995).
  4. T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 024007 (1998).
  5. G. Tóth, J. Comput. Phys. 161, 605 (2000).
  6. M. D. Duez, Y. T. Liu, S. L. Shapiro, M. Shibata, and B. C. Stephens, Phys. Rev. Lett. 96, 031101 (2006).
  7. M. D. Duez, Y. T. Liu, S. L. Shapiro, M. Shibata, and B. C. Stephens, Phys. Rev. D 73, 104015 (2006).
  8. B. C. Stephens, S. L. Shapiro, and Y. T. Liu, Phys. Rev. D 77, 044001 (2008).
  9. M. Shibata, Y. T. Liu, S. L. Shapiro, and B. C. Stephens, Phys. Rev. D 74, 104026 (2006).
  10. Y. T. Liu, S. L. Shapiro, and B. C. Stephens, Phys. Rev. D 76, 084017 (2007).
  11. Y. T. Liu, S. L. Shapiro, Z. B. Etienne, and K. Taniguchi, Phys. Rev. D 78, 024012 (2008).
  12. Z. B. Etienne, J. A. Faber, Y. T. Liu, S. L. Shapiro, K. Taniguchi, and T. W. Baumgarte, Phys. Rev. D 77, 084002 (2008).
  13. Z. B. Etienne, Y. T. Liu, S. L. Shapiro, and T. W. Baumgarte, Phys. Rev. D 79, 044024 (2009).
  14. V. Paschalidis, Z. Etienne, Y. T. Liu, and S. L. Shapiro, Phys. Rev. D 83, 064002 (2011).
  15. V. Paschalidis, Y. T. Liu, Z. Etienne, and S. L. Shapiro, Phys. Rev. D 84, 104032 (2011).
  16. B. D. Farris, Y. T. Liu, and S. L. Shapiro, Phys. Rev. D 81, 084008 (2010).
  17. B. D. Farris, Y. T. Liu, and S. L. Shapiro, Phys. Rev. D 84, 024024 (2011).
  18. J. U. Brackbill and D. C. Barnes, J. Comput. Phys. 35, 426 (1980).
  19. D. S. Balsara and D. S. Spicer, J. Comput. Phys. 149, 270 (1999).
  20. C. R. Evans and J. F. Hawley, Astrophys. J. 332, 659 (1988).
  21. A. Dedner, F. Kemm, D. Kröner, C. Munz, T. Schnitzer, and M. Wesenberg, J. Comput. Phys. 175, 645 (2002).
  22. M. Anderson, E. W. Hirschmann, S. L. Liebling, D. Neilsenand , Classical Quantum Gravity 23, 6503 (2006).
  23. M. Campanelli, C. O. Lousto, P. Marronetti, and Y. Zlochower, Phys. Rev. Lett. 96, 111101 (2006).
  24. J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and J. van Meter, Phys. Rev. Lett. 96, 111102 (2006).
  25. J. D. Brown, Phys. Rev. D 80, 084042 (2009).
  26. Z. B. Etienne, J. A. Faber, Y. T. Liu, S. L. Shapiro, and T. W. Baumgarte, Phys. Rev. D 76, 101503(R) (2007).
  27. D. Brown, O. Sarbach, E. Schnetter, M. Tiglio, P. Diener, I. Hawke, and D. Pollney, Phys. Rev. D 76, 081503 (2007).
  28. D. Brown, P. Diener, O. Sarbach, E. Schnetter, and M. Tiglio, Phys. Rev. D 79, 044023 (2009).
  29. D. S. Balsara, J. Comput. Phys. 174, 614 (2001).
  30. D. S. Balsara, J. Comput. Phys. 228, 5040 (2009).
  31. C. P. McNally, Mon. Not. R. Astron. Soc. 413, L76 (2011).
  32. L. Del Zanna, N. Bucciantini, and P. Londrillo, Astron. Astrophys. Rev. 400, 397 (2003).
  33. B. Giacomazzo, L. Rezzolla, and L. Baiotti, Phys. Rev. D 83, 044014 (2011).
  34. Note that the EM Lorenz gauge condition is named after Ludvig Lorenz and not after Hendrik Lorentz known for the Lorentz transformations (see also the note in p. 294 in [35]).

  35. J. D. Jackson, Classical Electrodynamics, edited by John David Jackson (Wiley-VCH, New York, 1998), 3rd ed., [Amazon][WorldCat].
  36. Z. B. Etienne, Y. T. Liu, V. Paschalidis, and S. L. Shapiro, arXiv:1112.0568.
  37. T. W. L. Baumgarte and S. L. Shapiro, Numerical Relativity (Cambridge University Press, Cambridge, England, 2010).
  38. J. R. van Meter, J. G. Baker, M. Koppitz, and D.-I. Choi, Phys. Rev. D 73, 124011 (2006).
  39. T. W. Baumgarte and S. L. Shapiro, Astrophys. J. 585, 921 (2003).
  40. B. Gustafsson, H.-O. Kreiss, and J. Oliger, Time Dependent Problems and Difference Methods (John Willey and Sons, New York, 1995).
  41. H.-O. Kreiss and J. Lorenz, Initial-Boundary Value Problems and the Navier-Stokes Equations (Academic Press, New York, 1989).
  42. Strictly speaking, for the quasilinear, partial differential equations considered here, the existence of real eigenvalues and a complete set of eigenvectors are only necessary conditions for strong hyperbolicity (see [40, 41] for complete definitions of strong hyperbolicity). However, here we are not interested in rigorously proving strong hyperbolicity. We are only interested in identifying the eigenvalues of the characteristic matrix and, in particular, whether static modes exist or not.

  43. M. Alcubierre, G. Allen, B. Brügmann, E. Seidel, and W.-M. Suen, Phys. Rev. D 62, 124011 (2000).
  44. H.-a. Shinkai and G. Yoneda, arXiv:gr-qc/0209111.
  45. V. Paschalidis, Phys. Rev. D 78, 024002 (2008).
  46. V. Paschalidis, J. Hansen, and A. Khokhlov, Phys. Rev. D 78, 064048 (2008).
  47. K. Taniguchi, T. W. Baumgarte, J. A. Faber, and S. L. Shapiro, Phys. Rev. D 77, 044003 (2008).
  48. J.-P. De Villiers, J. F. Hawley, and J. H. Krolik, Astrophys. J. 599, 1238 (2003).
  49. J. C. McKinney and C. F. Gammie, Astrophys. J. 611, 977 (2004).
  50. M. Anderson, E. W. Hirschmann, L. Lehner, S. L. Liebling, P. M. Motl, D. Neilsen, C. Palenzuela, J. E. Tohlineand , Phys. Rev. Lett. 100, 191101 (2008).
  51. S. Chawla, M. Anderson, M. Besselman, L. Lehner, S. L. Liebling, P. M. Motl, and D. Neilsen, Phys. Rev. Lett. 105, 111101 (2010).
  52. S. C. Noble, C. F. Gammie, J. C. McKinney, and L. Del Zanna, Astrophys. J. 641, 626 (2006).
  53. S. C. Noble, C. F. Gammie, J. C. McKinney, and L. Del Zanna (2006), publicly available on http://rainman.astro.illinois.edu/codelib/.
  54. K. Beckwith and J. M. Stone, Astrophys. J. Suppl. Ser. 193, 6 (2011).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation