scholarly journals Maximum-principle-satisfying and positivity-preserving high-order schemes for conservation laws: survey and new developments

Author(s):  
Xiangxiong Zhang ◽  
Chi-Wang Shu

In an earlier study (Zhang & Shu 2010 b J. Comput. Phys. 229 , 3091–3120 ( doi:10.1016/j.jcp.2009.12.030 )), genuinely high-order accurate finite volume and discontinuous Galerkin schemes satisfying a strict maximum principle for scalar conservation laws were developed. The main advantages of such schemes are their provable high-order accuracy and their easiness for generalization to multi-dimensions for arbitrarily high-order schemes on structured and unstructured meshes. The same idea can be used to construct high-order schemes preserving the positivity of certain physical quantities, such as density and pressure for compressible Euler equations, water height for shallow water equations and density for Vlasov–Boltzmann transport equations. These schemes have been applied in computational fluid dynamics, computational astronomy and astrophysics, plasma simulation, population models and traffic flow models. In this paper, we first review the main ideas of these maximum-principle-satisfying and positivity-preserving high-order schemes, then present a simpler implementation which will result in a significant reduction of computational cost especially for weighted essentially non-oscillatory finite-volume schemes.

2015 ◽  
Vol 18 (4) ◽  
pp. 901-930 ◽  
Author(s):  
Ziyao Sun ◽  
Honghui Teng ◽  
Feng Xiao

AbstractThis 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 (multi-moment Constrained finite Volume with WENO limiter of 4th order) method, is an extension of the MCV method of Ii & 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 WENO (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 WENO 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.


Mathematics ◽  
2021 ◽  
Vol 9 (16) ◽  
pp. 1885
Author(s):  
Shijian Lin ◽  
Qi Luo ◽  
Hongze Leng ◽  
Junqiang Song

We propose a new multi-moment numerical solver for hyperbolic conservation laws by using the alternating polynomial reconstruction approach. Unlike existing multi-moment schemes, our approach updates model variables by implementing two polynomial reconstructions alternately. First, Hermite interpolation reconstructs the solution within the cell by matching the point-based variables containing both physical values and their spatial derivatives. Then the reconstructed solution is updated by the Euler method. Second, we solve a constrained least-squares problem to correct the updated solution to preserve the conservation laws. Our method enjoys the advantages of a compact numerical stencil and high-order accuracy. Fourier analysis also indicates that our method allows a larger CFL number compared with many other high-order schemes. By adding a proper amount of artificial viscosity, shock waves and other discontinuities can also be computed accurately and sharply without solving an approximated Riemann problem.


Sign in / Sign up

Export Citation Format

Share Document