A Slope Constrained 4th Order Multi-Moment Finite Volume Method with WENO Limiter

This paper presents a new and better suited formulation to implement the limiting projection to high-order schemes that make use of high-order local reconstructions for hyperbolic conservation laws. The scheme, so-called MCV-WENO4 (multimoment Constrained finite Volume with WENO limiter of 4th order) method, is an extension of the MCV method of Ii a Xiao (2009) by adding the 1st order derivative (gradient or slope) at the cell center as an additional constraint for the cell-wise local reconstruction. The gradient is computed from a limiting projection using the wrii7,No (weighted essentially non-oscillatory) reconstruction that is built from the nodal values at 5 solution points within 3 neighboring cells. Different from other existing methods where only the cell-average value is used in the WINO reconstruction, the present method takes account of the solution structure within each mesh cell, and thus minimizes the stencil for reconstruction. The resulting scheme has 4th-order accuracy and is of significant advantage in algorithmic simplicity and computational efficiency. Numerical results of one and two dimensional benchmark tests for scalar and Euler conservation laws are shown to verify the accuracy and oscillation-less property of the scheme.

[1]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[2]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[3]  Wai-Sun Don,et al.  An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws , 2008, J. Comput. Phys..

[4]  Feng Xiao,et al.  High order multi-moment constrained finite volume method. Part I: Basic formulation , 2009, J. Comput. Phys..

[5]  Feng Xiao,et al.  CIP/multi-moment finite volume method for Euler equations: A semi-Lagrangian characteristic formulation , 2007, J. Comput. Phys..

[6]  Feng Xiao,et al.  Global shallow water models based on multi-moment constrained finite volume method and three quasi-uniform spherical grids , 2014, J. Comput. Phys..

[7]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[8]  Feng Xiao,et al.  A note on the general multi-moment constrained flux reconstruction formulation for high order schemes , 2012 .

[9]  Feng Xiao,et al.  Constructing a multi-dimensional oscillation preventing scheme for the advection equation by a rational function , 1996 .

[10]  T. Yabe,et al.  Completely conservative and oscillationless semi-Lagrangian schemes for advection transportation , 2001 .

[11]  Feng Xiao,et al.  A Multimoment Constrained Finite-Volume Model for Nonhydrostatic Atmospheric Dynamics , 2013 .

[12]  Zhi J. Wang,et al.  High-Order Multidomain Spectral Difference Method for the Navier-Stokes Equations , 2006 .

[13]  Chi-Wang Shu,et al.  Runge-Kutta Discontinuous Galerkin Method Using WENO Limiters , 2005, SIAM J. Sci. Comput..

[14]  Feng Xiao,et al.  Shallow water model on cubed-sphere by multi-moment finite volume method , 2008, J. Comput. Phys..

[15]  Steven J. Ruuth,et al.  A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods , 2002, SIAM J. Numer. Anal..

[16]  A. Patera.  A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[17]  Zhi J. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids. Basic Formulation , 2002 .

[18]  Feng Xiao,et al.  A multi-moment finite volume method for incompressible Navier-Stokes equations on unstructured grids: Volume-average/point-value formulation , 2014, J. Comput. Phys..

[19]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[20]  Z. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids , 2002 .

[21]  Feng Xiao,et al.  Constructing oscillation preventing scheme for advection equation by rational function , 1996 .

[22]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[23]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .

[24]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

[25]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[26]  Feng Xiao,et al.  Two and three dimensional multi-moment finite volume solver for incompressible Navier–Stokes equations on unstructured grids with arbitrary quadrilateral and hexahedral elements , 2014 .

[27]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one-dimensional case , 2004 .

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

[29]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[30]  R. LeVeque.  Approximate Riemann Solvers , 1992 .

[31]  Eleuterio F. Toro,et al.  ADER: Arbitrary High Order Godunov Approach , 2002, J. Sci. Comput..

[32]  John H. Kolias,et al.  A CONSERVATIVE STAGGERED-GRID CHEBYSHEV MULTIDOMAIN METHOD FOR COMPRESSIBLE FLOWS , 1995 .

[33]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

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

[35]  Zhi Jian Wang,et al.  A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids , 2009, J. Comput. Phys..

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

[37]  Eleuterio F. Toro,et al.  Derivative Riemann solvers for systems of conservation laws and ADER methods , 2006, J. Comput. Phys..

[38]  Chi-Wang Shu,et al.  A simple weighted essentially nonoscillatory limiter for Runge-Kutta discontinuous Galerkin methods , 2013, J. Comput. Phys..

[39]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

[40]  T. Yabe,et al.  Conservative and oscillation-less atmospheric transport schemes based on rational functions , 2002 .

[41]  Gecheng Zha,et al.  Improvement of weighted essentially non-oscillatory schemes near discontinuities , 2009 .

[42]  Feng Xiao,et al.  Fifth Order Multi-moment WENO Schemes for Hyperbolic Conservation Laws , 2015, J. Sci. Comput..

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

[44]  Mengping Zhang,et al.  A simple weighted essentially non-oscillatory limiter for the correction procedure via reconstruction (CPR) framework , 2015 .

[45]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

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

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