Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method
We present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field (∇⋅B=0) on adaptively refined, conformally moving meshes. The method relies on evolving the magnetic vector potential and then...
Main Authors: | , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Elsevier
2019-03-01
|
Series: | Journal of Computational Physics: X |
Online Access: | http://www.sciencedirect.com/science/article/pii/S2590055219300368 |
id |
doaj-d92887276f4f46529037a28d8738c894 |
---|---|
record_format |
Article |
spelling |
doaj-d92887276f4f46529037a28d8738c8942020-11-25T00:21:02ZengElsevierJournal of Computational Physics: X2590-05522019-03-012Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential methodP. Chris Fragile0Daniel Nemergut1Payden L. Shaw2Peter Anninos3Department of Physics and Astronomy, College of Charleston, Charleston, SC 29424, USA; Corresponding author.Department of Physics and Astronomy, College of Charleston, Charleston, SC 29424, USADepartment of Physics and Astronomy, College of Charleston, Charleston, SC 29424, USALawrence Livermore National Laboratory, P.O. Box 808, Livermore, CA 94550, USAWe present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field (∇⋅B=0) on adaptively refined, conformally moving meshes. The method relies on evolving the magnetic vector potential and then using it to reconstruct the magnetic fields. The advantage of this approach is that the vector potential is not subject to a constraint equation in the same way the magnetic field is, and so can be refined and moved in a straightforward way. We test this new method against a wide array of problems from simple Alfvén waves on a uniform grid to general relativistic MHD simulations of black hole accretion on a nested, spherical-polar grid. We find that the code produces accurate results and in all cases maintains a divergence-free magnetic field to machine precision. Keywords: Magnetohydrodynamics, Adaptive mesh refinement, Divergence constraint, Astrophysics, Magnetic fieldshttp://www.sciencedirect.com/science/article/pii/S2590055219300368 |
collection |
DOAJ |
language |
English |
format |
Article |
sources |
DOAJ |
author |
P. Chris Fragile Daniel Nemergut Payden L. Shaw Peter Anninos |
spellingShingle |
P. Chris Fragile Daniel Nemergut Payden L. Shaw Peter Anninos Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method Journal of Computational Physics: X |
author_facet |
P. Chris Fragile Daniel Nemergut Payden L. Shaw Peter Anninos |
author_sort |
P. Chris Fragile |
title |
Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method |
title_short |
Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method |
title_full |
Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method |
title_fullStr |
Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method |
title_full_unstemmed |
Divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method |
title_sort |
divergence-free magnetohydrodynamics on conformally moving, adaptive meshes using a vector potential method |
publisher |
Elsevier |
series |
Journal of Computational Physics: X |
issn |
2590-0552 |
publishDate |
2019-03-01 |
description |
We present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field (∇⋅B=0) on adaptively refined, conformally moving meshes. The method relies on evolving the magnetic vector potential and then using it to reconstruct the magnetic fields. The advantage of this approach is that the vector potential is not subject to a constraint equation in the same way the magnetic field is, and so can be refined and moved in a straightforward way. We test this new method against a wide array of problems from simple Alfvén waves on a uniform grid to general relativistic MHD simulations of black hole accretion on a nested, spherical-polar grid. We find that the code produces accurate results and in all cases maintains a divergence-free magnetic field to machine precision. Keywords: Magnetohydrodynamics, Adaptive mesh refinement, Divergence constraint, Astrophysics, Magnetic fields |
url |
http://www.sciencedirect.com/science/article/pii/S2590055219300368 |
work_keys_str_mv |
AT pchrisfragile divergencefreemagnetohydrodynamicsonconformallymovingadaptivemeshesusingavectorpotentialmethod AT danielnemergut divergencefreemagnetohydrodynamicsonconformallymovingadaptivemeshesusingavectorpotentialmethod AT paydenlshaw divergencefreemagnetohydrodynamicsonconformallymovingadaptivemeshesusingavectorpotentialmethod AT peteranninos divergencefreemagnetohydrodynamicsonconformallymovingadaptivemeshesusingavectorpotentialmethod |
_version_ |
1725364279145660416 |