959 research outputs found
Relativistic MHD with Adaptive Mesh Refinement
This paper presents a new computer code to solve the general relativistic
magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh
refinement (AMR). The fluid equations are solved using a finite difference
Convex ENO method (CENO) in 3+1 dimensions, and the AMR is Berger-Oliger.
Hyperbolic divergence cleaning is used to control the
constraint. We present results from three flat space tests, and examine the
accretion of a fluid onto a Schwarzschild black hole, reproducing the Michel
solution. The AMR simulations substantially improve performance while
reproducing the resolution equivalent unigrid simulation results. Finally, we
discuss strong scaling results for parallel unigrid and AMR runs.Comment: 24 pages, 14 figures, 3 table
A multidimensional grid-adaptive relativistic magnetofluid code
A robust second order, shock-capturing numerical scheme for multi-dimensional
special relativistic magnetohydrodynamics on computational domains with
adaptive mesh refinement is presented. The base solver is a total variation
diminishing Lax-Friedrichs scheme in a finite volume setting and is combined
with a diffusive approach for controlling magnetic monopole errors. The
consistency between the primitive and conservative variables is ensured at all
limited reconstructions and the spatial part of the four velocity is used as a
primitive variable. Demonstrative relativistic examples are shown to validate
the implementation. We recover known exact solutions to relativistic MHD
Riemann problems, and simulate the shock-dominated long term evolution of
Lorentz factor 7 vortical flows distorting magnetic island chains.Comment: accepted for publication in Computer Physics Communication
High-order conservative finite difference GLM-MHD schemes for cell-centered MHD
We present and compare third- as well as fifth-order accurate finite
difference schemes for the numerical solution of the compressible ideal MHD
equations in multiple spatial dimensions. The selected methods lean on four
different reconstruction techniques based on recently improved versions of the
weighted essentially non-oscillatory (WENO) schemes, monotonicity preserving
(MP) schemes as well as slope-limited polynomial reconstruction. The proposed
numerical methods are highly accurate in smooth regions of the flow, avoid loss
of accuracy in proximity of smooth extrema and provide sharp non-oscillatory
transitions at discontinuities. We suggest a numerical formulation based on a
cell-centered approach where all of the primary flow variables are discretized
at the zone center. The divergence-free condition is enforced by augmenting the
MHD equations with a generalized Lagrange multiplier yielding a mixed
hyperbolic/parabolic correction, as in Dedner et al. (J. Comput. Phys. 175
(2002) 645-673). The resulting family of schemes is robust, cost-effective and
straightforward to implement. Compared to previous existing approaches, it
completely avoids the CPU intensive workload associated with an elliptic
divergence cleaning step and the additional complexities required by staggered
mesh algorithms. Extensive numerical testing demonstrate the robustness and
reliability of the proposed framework for computations involving both smooth
and discontinuous features.Comment: 32 pages, 14 figure, submitted to Journal of Computational Physics
(Aug 7 2009
- …