Towards the ultimate understanding of MUSCL: Pitfalls in achieving third-order accuracy

Abstract We present a proof by analysis and numerical results that Van Leer's MUSCL conservative scheme with the discretization parameter κ is third-order accurate for κ = 1 / 3 . We include both the original finite-volume MUSCL family, updating cell-averaged values of the solution, and the related finite-difference version, updating point values. The presentation is needed because in the CFD literature claims have been made that not κ = 1 / 3 but κ = 1 / 2 yields third-order accuracy, or even that no value of κ can yield third-order accuracy. These false claims are the consequence of mixing up finite-difference concepts with finite-volume concepts. In a series of Pitfalls, we show how incorrect conclusions can be drawn when pointwise values of the discrete solution are interchanged with cell-averaged values. All flawed schemes presented in the Pitfalls, and some correct ones for comparison, are tested numerically and shown to behave as predicted by the analysis. We conclude with firm recommendations on how to achieve third-order accuracy at all output times, or just in a steady state.

[1]  H. Nishikawa.  The QUICK scheme is a third‐order finite‐volume scheme with point‐valued numerical solutions , 2020, International Journal for Numerical Methods in Fluids.

[2]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme III. Upstream-centered finite-difference schemes for ideal compressible flow , 1977 .

[3]  Bram van Leer,et al.  Upwind and High-Resolution Methods for Compressible Flow: From Donor Cell to Residual-Distribution Schemes , 2003 .

[4]  Panagiotis Tsoutsanis,et al.  Extended bounds limiter for high-order finite-volume schemes on unstructured meshes , 2018, J. Comput. Phys..

[5]  Bram van Leer,et al.  High-order finite-volume methods for the shallow-water equations on the sphere , 2010, J. Comput. Phys..

[6]  Chunhua Sheng,et al.  A new method towards high-order weno schemes on structured and unstructured grids , 2020 .

[7]  Stephen M. Ruffin,et al.  Higher Order Cell Centered Finite Volume Schemes for Unstructured Cartesian Grids , 2018 .

[8]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[9]  Mengping Zhang,et al.  On the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes , 2011 .

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

[11]  W. K. Anderson,et al.  Comparison of Finite Volume Flux Vector Splittings for the Euler Equations , 1985 .

[12]  Hiroaki Nishikawa.  On False Accuracy Verification of UMUSCL Scheme , 2021 .

[13]  B. P. Leonard,et al.  Beyond first‐order upwinding: The ultra‐sharp alternative for non‐oscillatory steady‐state simulation of convection , 1990 .

[14]  R. Denaro.  TOWARDS A NEW MODEL‐FREE SIMULATION OF HIGH‐REYNOLDS‐FLOWS: LOCAL AVERAGE DIRECT NUMERICAL SIMULATION , 1996 .

[15]  Huamo Wu,et al.  NON-EXISTENCE OF THIRD ORDER MUSCL SCHEMES UNIFIED CONSTRUCTION OF ENO SCHEMES AND A NEW DISCONTINUITY SHARPENNING TECHNIQUE - STIFF SOURCE TERM APPROACH , 1998 .

[16]  Liang Cheng,et al.  A review on TVD schemes and a refined flux-limiter for steady-state calculations , 2015, J. Comput. Phys..

[17]  Hong Luo,et al.  A reconstructed discontinuous Galerkin method for the compressible Navier-Stokes equations on three-dimensional hybrid grids , 2017 .

[18]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[19]  B. P. Leonard.  Order of accuracy of QUICK and related convection-diffusion schemes , 1995 .

[20]  Carl Ollivier-Gooch,et al.  Higher-order unstructured finite volume RANS solution of turbulent compressible flows , 2017 .

[21]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[22]  Minye Liu,et al.  Computational study of convective–diffusive mixing in a microchannel mixer , 2011 .

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

[24]  Roger Alexander Falconer,et al.  Modified forms of the third-order convection, second-order diffusion scheme for the advection-diffusion equation , 1994 .

[25]  Clarence O. E. Burg,et al.  Higher Order Variable Extrapolation For Unstructured Finite Volume RANS Flow Solvers , 2005 .

[26]  Hong Luo,et al.  A Reconstructed Discontinuous Galerkin Method Based on a Hierarchical Hermite WENO Reconstruction for Compressible Flows on Tetrahedral Grids , 2012 .

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

[28]  Daniela di Serafino,et al.  SParC-LES: Enabling large eddy simulations with parallel sparse matrix computation tools , 2015, Comput. Math. Appl..

[29]  Robert J. MacKinnon,et al.  Equivalent versions of the quick scheme for finite‐difference and finite‐volume numerical methods , 1992 .

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

[31]  Robert W. Walters,et al.  Implicit flux-split schemes for the Euler equations , 1985 .

[32]  J. Fromm.  A method for reducing dispersion in convective difference schemes , 1968 .

[33]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[34]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[35]  G. D. van Albada,et al.  A comparative study of computational methods in cosmic gas dynamics , 1982 .

[36]  Hong Q. Yang,et al.  Development of Vertex-Centered High-Order Schemes and Implementation in FUN3D , 2015 .