945 research outputs found
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
A Constrained Transport Method for the Solution of the Resistive Relativistic MHD Equations
We describe a novel Godunov-type numerical method for solving the equations
of resistive relativistic magnetohydrodynamics. In the proposed approach, the
spatial components of both magnetic and electric fields are located at zone
interfaces and are evolved using the constrained transport formalism. Direct
application of Stokes' theorem to Faraday's and Ampere's laws ensures that the
resulting discretization is divergence-free for the magnetic field and
charge-conserving for the electric field. Hydrodynamic variables retain,
instead, the usual zone-centred representation commonly adopted in
finite-volume schemes. Temporal discretization is based on Runge-Kutta
implicit-explicit (IMEX) schemes in order to resolve the temporal scale
disparity introduced by the stiff source term in Ampere's law. The implicit
step is accomplished by means of an improved and more efficient Newton-Broyden
multidimensional root-finding algorithm. The explicit step relies on a
multidimensional Riemann solver to compute the line-averaged electric and
magnetic fields at zone edges and it employs a one-dimensional Riemann solver
at zone interfaces to update zone-centred hydrodynamic quantities. For the
latter, we introduce a five-wave solver based on the frozen limit of the
relaxation system whereby the solution to the Riemann problem can be decomposed
into an outer Maxwell solver and an inner hydrodynamic solver. A number of
numerical benchmarks demonstrate that our method is superior in stability and
robustness to the more popular charge-conserving divergence cleaning approach
where both primary electric and magnetic fields are zone-centered. In addition,
the employment of a less diffusive Riemann solver noticeably improves the
accuracy of the computations.Comment: 25 pages, 14 figure
A Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations
A discontinuous Galerkin method for the ideal 5 moment two-fluid plasma
system is presented. The method uses a second or third order discontinuous
Galerkin spatial discretization and a third order TVD Runge-Kutta time stepping
scheme. The method is benchmarked against an analytic solution of a dispersive
electron acoustic square pulse as well as the two-fluid electromagnetic shock
and existing numerical solutions to the GEM challenge magnetic reconnection
problem. The algorithm can be generalized to arbitrary geometries and three
dimensions. An approach to maintaining small gauge errors based on error
propagation is suggested.Comment: 40 pages, 18 figures
- …