Equations ( 1) are integrated in a periodic quadratic box of length 2Θ using adaptive mesh refinement with rectangu- The formation of current sheets in ideal incompressible magnetohydrodynamic flows in two dimensions is studied numerically using lar grids self-adjusting to the flow. In each rectang
Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
β Scribed by Dinshaw S. Balsara
- Publisher
- Elsevier Science
- Year
- 2001
- Tongue
- English
- Weight
- 567 KB
- Volume
- 174
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
β¦ Synopsis
Several physical systems, such as nonrelativistic and relativistic magnetohydrodynamics (MHD), radiation MHD, electromagnetics, and incompressible hydrodynamics, satisfy Stoke's law type equations for the divergence-free evolution of vector fields. In this paper we present a full-fledged scheme for the second-order accurate, divergence-free evolution of vector fields on an adaptive mesh refinement (AMR) hierarchy. We focus here on adaptive mesh MHD. However, the scheme has applicability to the other systems of equations mentioned above. The scheme is based on making a significant advance in the divergence-free reconstruction of vector fields. In that sense, it complements the earlier work of D. S. Balsara and D. S. Spicer (1999, J. Comput. Phys. 7, 270) where we discussed the divergence-free time-update of vector fields which satisfy Stoke's law type evolution equations. Our advance in divergence-free reconstruction of vector fields is such that it reduces to the total variation diminishing (TVD) property for one-dimensional evolution and yet goes beyond it in multiple dimensions. For that reason, it is extremely suitable for the construction of higher order Godunov schemes for MHD. Both the two-dimensional and three-dimensional reconstruction strategies are developed. A slight extension of the divergence-free reconstruction procedure yields a divergence-free prolongation strategy for prolonging magnetic fields on AMR hierarchies. Divergence-free restriction is also discussed. Because our work is based on an integral formulation, divergencefree restriction and prolongation can be carried out on AMR meshes with any integral refinement ratio, though we specialize the expressions for the most popular situation where the refinement ratio is two. Furthermore, we pay attention to the fact that in order to efficiently evolve the MHD equations on AMR hierarchies, the refined meshes must evolve in time with time steps that are a fraction of their parent mesh's time step. An electric field correction strategy is presented for use on AMR meshes. The electric field correction strategy helps preserve the divergence-free evolution of the magnetic field even when the time steps are subcycled on refined meshes. The above-mentioned innovations have been implemented in Balsara's RIEMANN framework for parallel, self-adaptive computational astrophysics, which supports both nonrelativistic and relativistic MHD. Several rigorous, three-dimensional AMR-MHD 614
π SIMILAR VOLUMES
A Cartesian adaptive mesh code for time-dependent, compressible hydrodynamics (HD) and idea1 magnetohydrodynamics (MHD) in three space dimensions has been developed. The strategy of multiple subgrid nesting is used for mesh refinement and is incorporated into an operator-split, mixed finite-differen