High-order central ENO finite-volume scheme for ideal MHD

A high-order accurate finite-volume scheme for the compressible ideal magnetohydrodynamics (MHD) equations is proposed. The high-order MHD scheme is based on a central essentially non-oscillatory (CENO) method combined with the generalized Lagrange multiplier divergence cleaning method for MHD. The CENO method uses k-exact multidimensional reconstruction together with a monotonicity procedure that switches from a high-order reconstruction to a limited low-order reconstruction in regions of discontinuous or under-resolved solution content. Both reconstructions are performed on central stencils, and the switching procedure is based on a smoothness indicator. The proposed high-order accurate MHD scheme can be used on general polygonal grids. A highly sophisticated parallel implementation of the scheme is described that is fourth-order accurate on two-dimensional dynamically-adaptive body-fitted structured grids. The hierarchical multi-block body-fitted grid permits grid lines to conform to curved boundaries. High-order accuracy is maintained at curved domain boundaries by employing high-order spline representations and constraints at the Gauss quadrature points for flux integration. Detailed numerical results demonstrate high-order convergence for smooth flows and robustness against oscillations for problems with shocks. A new MHD extension of the well-known Shu-Osher test problem is proposed to test the ability of the high-order MHD scheme to resolve small-scale flow features in the presence of shocks. The dynamic mesh adaptation capabilities of the approach are demonstrated using adaptive time-dependent simulations of the Orszag-Tang vortex problem with high-order accuracy and unprecedented effective resolution.

[1]  P. Roe,et al.  Divergence- and curl-preserving prolongation and restriction formulas , 2002 .

[2]  Chang-Hsien Tai,et al.  Design of optimally smoothing multistage schemes for the Euler equations , 1992 .

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

[4]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[5]  Eitan Tadmor,et al.  Non-oscillatory central schemes for one- and two-dimensional MHD equations. II: high-order semi-discrete schemes. , 2006 .

[6]  T. Mexia,et al.  Author ' s personal copy , 2009 .

[7]  E. Tadmor,et al.  Non-oscillatory central differencing for hyperbolic conservation laws , 1990 .

[8]  Andrea Mignone,et al.  High-order conservative finite difference GLM-MHD schemes for cell-centered MHD , 2010, J. Comput. Phys..

[9]  Kenneth G. Powell,et al.  AN APPROXIMATE RIEMANN SOLVER FOR MAGNETOHYDRODYNAMICS (That Works in More than One Dimension) , 1994 .

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

[11]  V. Venkatakrishnan.  On the accuracy of limiters and convergence to steady state solutions , 1993 .

[12]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

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

[14]  K. Germaschewski,et al.  Three-dimensional MHD high-resolution computations with CWENO employing adaptive mesh refinement , 2004 .

[15]  Scott Northrup,et al.  Parallel Implicit Adaptive Mesh Refinement Scheme for Body-Fitted Multi-Block Mesh , 2005 .

[16]  Clinton P. T. Groth,et al.  High-Order Central ENO Finite-Volume Scheme with Adaptive Mesh Refinement , 2007 .

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

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

[19]  H. De Sterck,et al.  Stationary two-dimensional magnetohydrodynamic flows with shocks: characteristic analysis and grid convergence study , 2001 .

[20]  Timothy J. Barth,et al.  Recent developments in high order K-exact reconstruction on unstructured meshes , 1993 .

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

[22]  C. Ollivier-Gooch,et al.  A high-order-accurate unstructured mesh finite-volume scheme for the advection-diffusion equation , 2002 .

[23]  Carl Ollivier-Gooch,et al.  Accuracy preserving limiter for the high-order accurate solution of the Euler equations , 2009, J. Comput. Phys..

[24]  Andrea Mignone,et al.  A second-order unsplit Godunov scheme for cell-centered MHD: The CTU-GLM scheme , 2009, J. Comput. Phys..

[25]  Clinton P. T. Groth,et al.  Parallel Adaptive Mesh Refinement Scheme for Turbulent Non-Premixed Combusting Flow Prediction , 2006 .

[26]  George Em Karniadakis,et al.  A Discontinuous Galerkin Method for the Viscous MHD Equations , 1999 .

[27]  Clinton P. T. Groth,et al.  A Mesh Adjustment Scheme for Embedded Boundaries , 2006 .

[28]  Lucian Ivan,et al.  Development of High-order CENO Finite-volume Schemes with Block-based Adaptive Mesh Refinement (AMR) , 2011 .

[29]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

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

[31]  Scott Northrup,et al.  Parallel solution-adaptive method for two dimensional non-premixed combusting flows , 2011 .

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

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

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

[35]  Andrea Lani,et al.  A finite volume implicit time integration method for solving the equations of ideal magnetohydrodynamics for the hyperbolic divergence cleaning approach , 2011, J. Comput. Phys..

[36]  W. Gropp,et al.  Using MPI-2nd Edition , 1999 .

[37]  Clinton P. T. Groth,et al.  International Journal of Computational Fluid Dynamics a Parallel Adaptive Mesh Refinement Algorithm for Predicting Turbulent Non-premixed Combusting Flows a Parallel Adaptive Mesh Refinement Algorithm for Predicting Turbulent Non-premixed Combusting Flows , 2022 .

[38]  A. Harten,et al.  Multi-Dimensional ENO Schemes for General Geometries , 1991 .

[39]  Clinton P. T. Groth,et al.  High-Order Solution-Adaptive Central Essentially Non-Oscillatory (CENO) Method for Viscous Flows , 2011 .

[40]  C. Angelopoulos.  High resolution schemes for hyperbolic conservation laws , 1992 .

[41]  Dimitri J. Mavriplis,et al.  Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes , 2003 .

[42]  Jean-Marc Hérard,et al.  A LOCAL TIME-STEPPING DISCONTINUOUS GALERKIN ALGORITHM FOR THE MHD SYSTEM , 2009 .

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

[44]  Harish Kumar,et al.  On the role of Riemann solvers in Discontinuous Galerkin methods for magnetohydrodynamics , 2010, J. Comput. Phys..

[45]  Clinton P. T. Groth,et al.  Parallel High-Order Anisotropic Block-Based Adaptive Mesh Refinement Finite-Volume Scheme , 2011 .

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

[47]  Dinshaw S. Balsara.  Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics , 2009, J. Comput. Phys..

[48]  J. Sachdev,et al.  A parallel solution-adaptive scheme for multi-phase core flows in solid propellant rocket motors , 2005 .

[49]  Claus-Dieter Munz,et al.  Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model , 2000 .

[50]  Michael Dumbser,et al.  Efficient, high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics , 2008, Journal of Computational Physics.

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

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

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

[54]  Eitan Tadmor,et al.  Nonoscillatory Central Schemes for One- and Two-Dimensional Magnetohydrodynamics Equations. II: High-Order SemiDiscrete Schemes , 2006, SIAM J. Sci. Comput..

[55]  Robert B. Ross,et al.  Using MPI-2: Advanced Features of the Message Passing Interface , 2003, CLUSTER.

[56]  Clinton P. T. Groth,et al.  A parallel solution - adaptive method for three-dimensional turbulent non-premixed combusting flows , 2010, J. Comput. Phys..