Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows , I : Basic Theory

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls t o limit the amount of and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvilinear grids. The four advantages of the present approach over existing MHD schemes reported in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion MHD flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike existing schemes for the conservative MHD equations which suffer from ill-conditioned eigendecompositions, the present scheme makes use of a well-conditioned eigen-decomposition obtained from a minor modification of the eigenvectors of the non-conservative MHD equations t o solve the conservative form of the MHD equations. Third, this approach of using the non-conservative eigensystem when solving the conservative equations also works well in the context of standard shock-capturing schemes for the MHD equations. Fourth, a new approach to minimize the numerical error of the divergence-free magnetic condition for high order schemes is introduced. Numerical experiments with typical MHD model problems revealed the applicability of the newly developed schemes for the MHD equations.

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

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

[3]  H. C. Yee Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows , 2003 .

[4]  Neil D. Sandham,et al.  Low-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters , 1999 .

[5]  H. C. Yee,et al.  Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations , 2001 .

[6]  K. Yee Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

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

[8]  H. C. Yee,et al.  Entropy Splitting and Numerical Dissipation , 2000 .

[9]  Björn Sjögreen,et al.  Designing Adaptive Low Dissipative High Order Schemes , 2003 .

[10]  Björn Sjögreen,et al.  Low dissipative high‐order numerical simulations of supersonic reactive flows , 2003 .

[12]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[13]  Gérard Gallice,et al.  Roe Matrices for Ideal MHD and Systematic Construction of Roe Matrices for Systems of Conservation Laws , 1997 .

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

[15]  Paul R. Woodward,et al.  A Simple Finite Difference Scheme for Multidimensional Magnetohydrodynamical Equations , 1998 .

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