Piecewise parabolic method on a local stencil for ideal magnetohydrodynamics

A numerical scheme based on the piecewise parabolic method on a local stencil (PPML) is proposed for solving the ideal magnetohydrodynamic (MHD) equations. The method makes use of the conservation of Riemann invariants along the characteristics of the MHD equations. As a result, a local stencil can be used to construct a numerical solution. This approach improves the dissipative properties of the numerical scheme and is convenient in the case of adaptive grids. The basic stages in the design of the scheme are illustrated in the two-dimensional case. The conservation of the solenoidal property of the magnetic field is discussed. The scheme is tested using several typical MHD problems.

[1]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[2]  P. Roe,et al.  A Solution-Adaptive Upwind Scheme for Ideal Magnetohydrodynamics , 1999 .

[3]  H. Huynh,et al.  Accurate Monotonicity-Preserving Schemes with Runge-Kutta Time Stepping , 1997 .

[4]  G. Tóth.  The ∇·B=0 Constraint in Shock-Capturing Magnetohydrodynamics Codes , 2000 .

[5]  Guang-Shan Jiang,et al.  A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics , 1999 .

[6]  J. Brackbill,et al.  The Effect of Nonzero ∇ · B on the numerical solution of the magnetohydrodynamic equations☆ , 1980 .

[7]  S. Orszag,et al.  Small-scale structure of two-dimensional magnetohydrodynamic turbulence , 1979, Journal of Fluid Mechanics.

[8]  Dinshaw Balsara,et al.  Second-Order-accurate Schemes for Magnetohydrodynamics with Divergence-free Reconstruction , 2003, astro-ph/0308249.

[9]  T. J. Barth,et al.  On unstructured grids and solvers , 1990 .

[10]  M. V. Popov,et al.  Piecewise parabolic method on local stencil for gasdynamic simulations , 2007 .

[11]  M. Brio,et al.  An upwind differencing scheme for the equations of ideal magnetohydrodynamics , 1988 .

[12]  P. Roe.  CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .

[13]  D. Balsara,et al.  A Staggered Mesh Algorithm Using High Order Godunov Fluxes to Ensure Solenoidal Magnetic Fields in Magnetohydrodynamic Simulations , 1999 .

[14]  Jianqiang Han,et al.  An adaptive moving mesh method for two-dimensional ideal magnetohydrodynamics , 2007, J. Comput. Phys..

[15]  Kenneth G. Powell,et al.  Axisymmetric modeling of cometary mass loading on an adaptively refined grid: MHD results , 1994 .

[16]  J. Hawley,et al.  Simulation of magnetohydrodynamic flows: A Constrained transport method , 1988 .

[17]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .