High-Order Finite-Volume Method with Block-Based AMR for Magnetohydrodynamics Flows

A high-order central essentially non-oscillatory (CENO) finite volume scheme combined with a block-based adaptive mesh refinement (AMR) algorithm is proposed for the solution of the ideal magnetohydrodynamics equations. The high-order CENO finite-volume scheme is implemented with fourth-order spatial accuracy within a flexible multi-block, body-fitted, hexahedral grid framework. An important feature of the high-order adaptive approach is that it allows for anisotropic refinement, which can lead to large computational savings when anisotropic flow features such as isolated propagating fronts and/or waves, shocks, shear surfaces, and current sheets are present in the flow. This approach is designed to handle complex multi-block grid configurations, including cubed-sphere grids, where some grid blocks may have degenerate edges or corners characterized by missing neighboring blocks. A procedure for building valid high-order reconstruction stencils, even at these degenerate block edges and corners, is proposed, taking into account anisotropic resolution changes in a systematic and general way. Furthermore, a non-uniform or heterogeneous block structure is used where the ghost cells of a block containing the solution content of neighboring blocks are stored directly at the resolution of the neighbors. A generalized Lagrange multiplier divergence correction technique is applied to achieve numerically divergence-free magnetic fields while preserving high-order accuracy on the anisotropic AMR grids. Parallel implementation and local grid adaptivity are achieved by using a hierarchical block-based domain partitioning strategy in which the connectivity and refinement history of grid blocks are tracked using a flexible binary tree data structure. Physics-based refinement criteria as well as the CENO smoothness indicator are both used for directing the mesh refinement. Numerical results, including solution-driven anisotropic refinement of cubed-sphere grids, are presented to demonstrate the accuracy and efficiency of the approach.

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

[2]  Michael Williamschen,et al.  Parallel Anisotropic Block-based Adaptive Mesh Refinement Algorithm For Three-dimensional Flows , 2013 .

[3]  Lars Hernquist,et al.  A constrained transport scheme for MHD on unstructured static and moving meshes , 2014, 1402.5963.

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

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

[6]  Luiz Tobaldini Neto,et al.  A High-Order Finite-Volume Scheme for Large-Eddy Simulation of Turbulent Premixed Flames , 2014 .

[7]  James A. Rossmanith,et al.  Finite difference weighted essentially non-oscillatory schemes with constrained transport for 2D ideal Magnetohydrodynamics , 2013, ICOPS 2013.

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

[9]  Randall J. LeVeque,et al.  Adjoint Methods for Guiding Adaptive Mesh Refinement in Tsunami Modeling , 2016, Pure and Applied Geophysics.

[10]  Xiao-Chuan Cai,et al.  Parallel fully implicit two‐grid methods for distributed control of unsteady incompressible flows , 2013 .

[11]  Chi-Wang Shu,et al.  Preface to the Republication of “Uniformly High Order Essentially Non-oscillatory Schemes, III,” by Harten, Engquist, Osher, and Chakravarthy , 1997 .

[12]  Hans De Sterck,et al.  High-order central ENO finite-volume scheme for hyperbolic conservation laws on three-dimensional cubed-sphere grids , 2015, J. Comput. Phys..

[13]  Hans De Sterck,et al.  A Parallel High-Order CENO Finite-Volume Scheme with AMR for Three-Dimensional Ideal MHD Flows , 2017 .

[14]  Phillip Colella,et al.  High-order, finite-volume methods in mapped coordinates , 2010, J. Comput. Phys..

[15]  Carsten Burstedde,et al.  p4est: Scalable Algorithms for Parallel Adaptive Mesh Refinement on Forests of Octrees , 2011, SIAM J. Sci. Comput..

[16]  Yuxi Chen,et al.  A fifth-order finite difference scheme for hyperbolic equations on block-adaptive curvilinear grids , 2016, J. Comput. Phys..

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

[18]  H. De Sterck,et al.  A High-Order Finite-Volume Method with Anisotropic AMR for Ideal MHD Flows , 2017 .

[19]  Bernd Einfeld.  On Godunov-type methods for gas dynamics , 1988 .

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

[21]  Rony Keppens,et al.  Parallel, grid-adaptive approaches for relativistic hydro and magnetohydrodynamics , 2012, J. Comput. Phys..

[22]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[23]  Clinton P. T. Groth,et al.  Anisotropic Non-Uniform Block-Based Adaptive Mesh Refinement for Three-Dimensional Inviscid and Viscous Flows , 2015 .

[24]  Quentin F. Stout,et al.  High performance computer methods applied to predictive space weather simulations , 2000 .

[25]  Quentin F. Stout,et al.  An adaptive MHD method for global space weather simulations , 2000 .

[26]  Poedts,et al.  Intermediate shocks in three-dimensional magnetohydrodynamic bow-shock flows with multiple interacting shock fronts , 2000, Physical review letters.

[27]  S. Osher,et al.  Regular ArticleUniformly High Order Accurate Essentially Non-oscillatory Schemes, III , 1997 .

[28]  Lucie Fréret,et al.  Block-Based Anisotropic AMR withA PosterioriAdjoint-Based Error Estimation for Three-Dimensional Inviscid and Viscous Flows , 2017 .

[29]  Clinton P. T. Groth,et al.  A computational framework for predicting laminar reactive flows with soot formation , 2010 .

[30]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[31]  F. Olsson,et al.  Stability of interpolation on overlapping grids , 1996 .

[32]  Phillip Colella,et al.  High-order finite-volume methods for hyperbolic conservation laws on mapped multiblock grids , 2015, J. Comput. Phys..

[33]  M. Berger.  ON CONSERVATION AT GRID INTERFACES. , 1987 .

[34]  David L. Darmofal,et al.  Anisotropic Adaptation for Functional Outputs of Viscous Flow Simulations , 2003 .

[35]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

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

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

[38]  Bram van Leer,et al.  Design of Optimally Smoothing Multi-Stage Schemes for the Euler Equations , 1989 .

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

[40]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[41]  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 .

[42]  Peter MacNeice,et al.  Paramesh: A Parallel Adaptive Mesh Refinement Community Toolkit , 2013 .

[43]  Bertram Taetz,et al.  A High-Order Unstaggered Constrained-Transport Method for the Three-Dimensional Ideal Magnetohydrodynamic Equations Based on the Method of Lines , 2013, SIAM J. Sci. Comput..

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

[45]  Phillip Colella,et al.  A HIGH-ORDER FINITE-VOLUME METHOD FOR CONSERVATION LAWS ON LOCALLY REFINED GRIDS , 2011 .

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

[47]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[48]  C. Groth,et al.  Application of Gaussian Moment Closure to Microscale Flows with Moving Embedded Boundaries , 2014 .

[49]  C. P. T. Groth,et al.  A Parallel Adaptive 3D MHD Scheme for Modeling Coronal and Solar Wind Plasma Flows , 1999 .

[50]  Chaopeng Shen,et al.  Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations , 2011, J. Comput. Phys..

[51]  D. Venditti,et al.  Anisotropic grid adaptation for functional outputs: application to two-dimensional viscous flows , 2003 .

[52]  D. D. Zeeuw,et al.  Global three‐dimensional MHD simulation of a space weather event: CME formation, interplanetary propagation, and interaction with the magnetosphere , 2000 .

[53]  Mark F. Adams,et al.  Chombo Software Package for AMR Applications Design Document , 2014 .

[54]  Quentin F. Stout,et al.  Adaptive numerical algorithms in space weather modeling , 2012, J. Comput. Phys..

[55]  Clinton P. T. Groth,et al.  A High-Order Central ENO Finite-Volume Scheme for Three-Dimensional Turbulent Reactive Flows on Unstructured Mesh , 2013 .

[56]  Björn Sjögreen,et al.  Conservative and Non-Conservative Interpolation between Overlapping Grids for Finite Volume Solutions of Hyperbolic Problems , 1994 .

[57]  Hans De Sterck,et al.  High-order central ENO finite-volume scheme for ideal MHD , 2013 .

[58]  Marsha Berger,et al.  Three-Dimensional Adaptive Mesh Refinement for Hyperbolic Conservation Laws , 1994, SIAM J. Sci. Comput..

[59]  Hans De Sterck,et al.  Multi-dimensional finite-volume scheme for hyperbolic conservation laws on three-dimensional solution-adaptive cubed-sphere grids , 2013, J. Comput. Phys..

[60]  Louis A. Povinelli,et al.  Large-scale computation of incompressible viscous flow by least-squares finite element method , 1994 .

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

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

[63]  V. Springel,et al.  A moving mesh unstaggered constrained transport scheme for magnetohydrodynamics , 2016, 1606.02310.

[64]  Michael Dumbser,et al.  Efficient conservative ADER schemes based on WENO reconstruction and space-time predictor in primitive variables , 2015, Computational astrophysics and cosmology.

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

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

[67]  Clinton P. T. Groth,et al.  High-order solution-adaptive central essentially non-oscillatory (CENO) method for viscous flows , 2011, J. Comput. Phys..

[68]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[69]  Clinton P. T. Groth,et al.  A High-Order Central ENO Finite-Volume Scheme for Three-Dimensional Low-Speed Viscous Flows on Unstructured Mesh , 2015 .