Piecewise parabolic method on a local stencil for magnetized supersonic turbulence simulation

Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new, higher order-accurate, low dissipation numerical method which requires no additional dissipation or local ''fixes'' for stable execution. We describe PPML, a local stencil variant of the popular PPM algorithm for solving the equations of compressible ideal magnetohydrodynamics. The principal difference between PPML and PPM is that cell interface states are evolved rather that reconstructed at every timestep, resulting in a compact stencil. Interface states are evolved using Riemann invariants containing all transverse derivative information. The conservation laws are updated in an unsplit fashion, making the scheme fully multidimensional. Divergence-free evolution of the magnetic field is maintained using the higher order-accurate constrained transport technique of Gardiner and Stone. The accuracy and stability of the scheme is documented against a bank of standard test problems drawn from the literature. The method is applied to numerical simulation of supersonic MHD turbulence, which is important for many problems in astrophysics, including star formation in dark molecular clouds. PPML accurately reproduces in three-dimensions a transition to turbulence in highly compressible isothermal gas in a molecular cloud model. The low dissipation and wide spectral bandwidth of this method make it an ideal candidate for direct turbulence simulations.

[1]  Paul R. Woodward Numerics for iles: The ppm compressible gas dynamics scheme , 2007 .

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

[3]  Richard Sanders,et al.  High resolution staggered mesh approach for nonlinear hyperbolic systems of conser-vation laws , 1992 .

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

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

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

[7]  Alan Weiser,et al.  A high order staggered grid method for hyperbolic systems of conservation laws in one space dimension , 1989 .

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

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

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

[11]  Michael L. Norman,et al.  Simulations of Supersonic Turbulence in Molecular Clouds: Evidence for a New Universality , 2009 .

[12]  A. E. Wright,et al.  THEORY OF STAR FORMATION. , 1970 .

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

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

[15]  A. Schekochihin,et al.  Turbulence and Magnetic Fields in Astrophysical Plasmas , 2005, astro-ph/0507686.

[16]  J. Stone,et al.  An unsplit Godunov method for ideal MHD via constrained transport , 2005, astro-ph/0501557.

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

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

[19]  Björn Sjögreen,et al.  Development of low dissipative high order filter schemes for multiscale Navier-Stokes/MHD systems , 2006, J. Comput. Phys..

[20]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[21]  J. Quirk A Contribution to the Great Riemann Solver Debate , 1994 .

[22]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[23]  K. Kusano,et al.  A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics , 2005 .

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

[25]  Michael L. Norman,et al.  Scaling laws and intermittency in highly compressible turbulence , 2007, 0706.0739.

[26]  L. Margolin,et al.  Implicit Large Eddy Simulation: Computing Turbulent Fluid Dynamics , 2011 .

[27]  Richard Sanders,et al.  A third-order accurate variation nonexpansive difference scheme for single nonlinear conservation laws , 1988 .

[28]  William J. Rider,et al.  Accurate monotonicity- and extrema-preserving methods through adaptive nonlinear hybridizations , 2007, J. Comput. Phys..

[29]  Phillip Colella,et al.  A limiter for PPM that preserves accuracy at smooth extrema , 2008, J. Comput. Phys..

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

[31]  M. V. Popov,et al.  Piecewise parabolic method on a local stencil for ideal magnetohydrodynamics , 2008 .

[32]  Dongwook Lee,et al.  An unsplit staggered mesh scheme for multidimensional magnetohydrodynamics , 2009, J. Comput. Phys..

[33]  M. Carbonaro,et al.  von Karman Institute for Fluid Dynamics , 2004 .

[34]  M. Norman,et al.  The Statistics of Supersonic Isothermal Turbulence , 2007, 0704.3851.

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

[36]  Dongsu Ryu,et al.  Numerical magnetohydrodynamics in astrophysics: Algorithm and tests for multidimensional flow , 1995 .

[37]  Paul R. Woodward,et al.  An approximate Riemann solver for ideal magnetohydrodynamics , 1994 .

[38]  Vincent Guinot,et al.  High resolution Godunov‐type schemes with small stencils , 2004 .

[39]  Paul R. Woodward,et al.  Extension of the Piecewise Parabolic Method to Multidimensional Ideal Magnetohydrodynamics , 1994 .