High-Order Low-Dissipation Targeted ENO Schemes for Ideal Magnetohydrodynamics

AbstractThe recently proposed targeted ENO (TENO) schemes (Fu et al. J Comput Phys 305:333–359, 2016) are demonstrated to feature the controllable low numerical dissipation and sharp shock-capturing property in compressible gas dynamic simulations. However, the application of low-dissipation TENO schemes to ideal magnetohydrodynamics (MHD) is not straightforward. The complex interaction between fluid mechanics and electromagnetism induces extra numerical challenges, including simultaneously preserving the ENO-property, maintaining good numerical robustness and low dissipation as well as controlling divergence errors. In this paper, based on an unstaggered constrained transport framework to control the divergence error, we extend a set of high-order low-dissipation TENO schemes ranging from 5-point to 8-point stencils to solving the ideal MHD equations. A unique set of built-in parameters for each TENO scheme is determined. Different from the TENO schemes in Fu et al.  (2016), a modified scale-separation formula is developed. The new formula can achieve stronger scale separation, and it is simpler and more efficient than the previous version as the computation cost of high-order global smoothness measure $${\tau _K}$$τK is avoided. The performances of tailored schemes are systematically studied by several benchmark simulations. Numerical experiments demonstrate that the TENO schemes in the constrained transport framework are promising to simulate more complex MHD flows.

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

[2]  Sergey Yakovlev,et al.  Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field , 2011, J. Comput. Phys..

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

[4]  Klaus A. Hoffmann,et al.  Minimizing errors from linear and nonlinear weights of WENO scheme for broadband applications with shock waves , 2013, J. Comput. Phys..

[5]  Bertram Taetz,et al.  An unstaggered constrained transport method for the 3D ideal magnetohydrodynamic equations , 2010, J. Comput. Phys..

[6]  James A. Rossmanith,et al.  An Unstaggered, High-Resolution Constrained Transport Method for Magnetohydrodynamic Flows , 2006, SIAM J. Sci. Comput..

[7]  Hans d,et al.  Multi-dimensional upwind constrained transport on unstructured grids for 'shallow water' magnetohydrodynamics , 2001 .

[8]  Dinshaw S. Balsara.  A two-dimensional HLLC Riemann solver for conservation laws: Application to Euler and magnetohydrodynamic flows , 2012, J. Comput. Phys..

[9]  Michael Dumbser,et al.  A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems , 2016, J. Comput. Phys..

[10]  C. Munz,et al.  Hyperbolic divergence cleaning for the MHD equations , 2002 .

[11]  Feng Xiao,et al.  Boundary Variation Diminishing (BVD) reconstruction: A new approach to improve Godunov schemes , 2016, J. Comput. Phys..

[12]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[13]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[14]  Xiao Feng,et al.  A high-order positivity-preserving single-stage single-step method for the ideal magnetohydrodynamic equations , 2015, J. Comput. Phys..

[15]  Nikolaus A. Adams,et al.  Targeted ENO schemes with tailored resolution property for hyperbolic conservation laws , 2017, J. Comput. Phys..

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

[17]  Nikolaus A. Adams,et al.  A new class of adaptive high-order targeted ENO schemes for hyperbolic conservation laws , 2018, J. Comput. Phys..

[18]  Liwei Xu,et al.  Positivity-preserving DG and central DG methods for ideal MHD equations , 2013, J. Comput. Phys..

[19]  Dinshaw Balsara,et al.  Divergence-free adaptive mesh refinement for Magnetohydrodynamics , 2001 .

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

[21]  Shiyi Chen,et al.  A hybrid scheme for compressible magnetohydrodynamic turbulence , 2016, J. Comput. Phys..

[22]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[23]  Lin Fu,et al.  A low-dissipation finite-volume method based on a new TENO shock-capturing scheme , 2019, Comput. Phys. Commun..

[24]  Zhengfu Xu,et al.  Positivity-Preserving Finite Difference Weighted ENO Schemes with Constrained Transport for Ideal Magnetohydrodynamic Equations , 2015, SIAM J. Sci. Comput..

[25]  Yuxin Ren,et al.  A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws , 2003 .

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

[27]  Bruno Costa,et al.  An improved WENO-Z scheme , 2016, J. Comput. Phys..

[28]  Dinshaw Balsara,et al.  A Comparison between Divergence-Cleaning and Staggered-Mesh Formulations for Numerical Magnetohydrodynamics , 2003 .

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

[30]  Steven H. Frankel,et al.  Numerical simulations of compressible multicomponent and multiphase flow using a high-order targeted ENO (TENO) finite-volume method , 2017 .

[31]  Nikolaus A. Adams,et al.  A family of high-order targeted ENO schemes for compressible-fluid simulations , 2016, J. Comput. Phys..

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

[33]  Manuel Torrilhon,et al.  A Constrained Transport Upwind Scheme for Divergence-free Advection , 2003 .

[34]  Xiao Feng,et al.  A high-order finite difference WENO scheme for ideal magnetohydrodynamics on curvilinear meshes , 2017, 2018 IEEE International Conference on Plasma Science (ICOPS).

[35]  Chi-Wang Shu,et al.  Locally Divergence-Free Discontinuous Galerkin Methods for MHD Equations , 2005, J. Sci. Comput..

[36]  Qi Tang,et al.  Finite difference weighted essentially non-oscillatory schemes with constrained transport for ideal magnetohydrodynamics , 2013, J. Comput. Phys..

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

[38]  Soshi Kawai,et al.  Divergence-free-preserving high-order schemes for magnetohydrodynamics: An artificial magnetic resistivity method , 2013, J. Comput. Phys..

[39]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[40]  Manuel Torrilhon,et al.  Uniqueness conditions for Riemann problems of ideal magnetohydrodynamics , 2003 .

[41]  Björn Sjögreen,et al.  Adaptive filtering and limiting in compact high order methods for multiscale gas dynamics and MHD systems , 2008 .

[42]  Björn Sjögreen,et al.  Non-Linear Filtering and Limiting in High Order Methods for Ideal and Non-Ideal MHD , 2006, J. Sci. Comput..

[43]  Andree Susanto,et al.  High-Order Finite-Volume Schemes for Magnetohydrodynamics , 2014 .

[44]  J. M. Picone,et al.  Evolution of the Orszag-Tang vortex system in a compressible medium. I: Initial average subsonic flow , 1989 .