1,037 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 Second-Order Unsplit Godunov Scheme for Cell-Centered MHD: the CTU-GLM scheme
We assess the validity of a single step Godunov scheme for the solution of
the magneto-hydrodynamics equations in more than one dimension. The scheme is
second-order accurate and the temporal discretization is based on the
dimensionally unsplit Corner Transport Upwind (CTU) method of Colella. The
proposed scheme employs a cell-centered representation of the primary fluid
variables (including magnetic field) and conserves mass, momentum, magnetic
induction and energy. A variant of the scheme, which breaks momentum and energy
conservation, is also considered. Divergence errors are transported out of the
domain and damped using the mixed hyperbolic/parabolic divergence cleaning
technique by Dedner et al. (J. Comput. Phys., 175, 2002). The strength and
accuracy of the scheme are verified by a direct comparison with the eight-wave
formulation (also employing a cell-centered representation) and with the
popular constrained transport method, where magnetic field components retain a
staggered collocation inside the computational cell. Results obtained from two-
and three-dimensional test problems indicate that the newly proposed scheme is
robust, accurate and competitive with recent implementations of the constrained
transport method while being considerably easier to implement in existing hydro
codes.Comment: 31 Pages, 16 Figures Accepted for publication in Journal of
Computational Physic
An Unsplit, Cell-Centered Godunov Method for Ideal MHD
We present a second-order Godunov algorithm for multidimensional, ideal MHD.
Our algorithm is based on the unsplit formulation of Colella (J. Comput. Phys.
vol. 87, 1990), with all of the primary dependent variables centered at the
same location. To properly represent the divergence-free condition of the
magnetic fields, we apply a discrete projection to the intermediate values of
the field at cell faces, and apply a filter to the primary dependent variables
at the end of each time step. We test the method against a suite of linear and
nonlinear tests to ascertain accuracy and stability of the scheme under a
variety of conditions. The test suite includes rotated planar linear waves, MHD
shock tube problems, low-beta flux tubes, and a magnetized rotor problem. For
all of these cases, we observe that the algorithm is second-order accurate for
smooth solutions, converges to the correct weak solution for problems involving
shocks, and exhibits no evidence of instability or loss of accuracy due to the
possible presence of non-solenoidal fields.Comment: 37 Pages, 9 Figures, submitted to Journal of Computational Physic
Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
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. 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 Balsara and Spicer (1999) 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. Divergence-free restriction is also discussed. 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 sub-cycled on refined
meshes. The above-mentioned innovations have been implemented in Balsara's
RIEMANN framework for parallel, self-adaptive computational astrophysics which
supports both non-relativistic and relativistic MHD. Several rigorous, three
dimensional AMR-MHD test problems with strong discontinuities have been run
with the RIEMANN framework showing that the strategy works very well.Comment: J.C.P., figures of reduced qualit
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
WhiskyMHD: a new numerical code for general relativistic magnetohydrodynamics
The accurate modelling of astrophysical scenarios involving compact objects
and magnetic fields, such as the collapse of rotating magnetized stars to black
holes or the phenomenology of gamma-ray bursts, requires the solution of the
Einstein equations together with those of general-relativistic
magnetohydrodynamics. We present a new numerical code developed to solve the
full set of general-relativistic magnetohydrodynamics equations in a dynamical
and arbitrary spacetime with high-resolution shock-capturing techniques on
domains with adaptive mesh refinements. After a discussion of the equations
solved and of the techniques employed, we present a series of testbeds carried
out to validate the code and assess its accuracy. Such tests range from the
solution of relativistic Riemann problems in flat spacetime, over to the
stationary accretion onto a Schwarzschild black hole and up to the evolution of
oscillating magnetized stars in equilibrium and constructed as consistent
solutions of the coupled Einstein-Maxwell equations.Comment: minor changes to match the published versio
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
Positivity-Preserving Finite Difference WENO Schemes with Constrained Transport for Ideal Magnetohydrodynamic Equations
In this paper, we utilize the maximum-principle-preserving flux limiting
technique, originally designed for high order weighted essentially
non-oscillatory (WENO) methods for scalar hyperbolic conservation laws, to
develop a class of high order positivity-preserving finite difference WENO
methods for the ideal magnetohydrodynamic (MHD) equations. Our schemes, under
the constrained transport (CT) framework, can achieve high order accuracy, a
discrete divergence-free condition and positivity of the numerical solution
simultaneously. Numerical examples in 1D, 2D and 3D are provided to demonstrate
the performance of the proposed method.Comment: 21 pages, 28 figure
- âŠ