Improved multislope MUSCL reconstruction on unstructured grids for shallow water equations. (12th March 2018)

• Record Type: Journal Article
• Authors: Zhao, J.; Özgen, I.; Liang, D.; Hinkelmann, R.
• Is Part Of: International journal for numerical methods in fluids. Volume 87:Number 8(2018)
• Page Start: 401
• Page End: 436
• Publication Date: 2018-03-12
• Subjects: finite‐volume method -- MUSCL scheme -- shallow water -- slope limiters -- total variation diminishing -- unstructured mesh; Fluid dynamics -- Mathematics -- Periodicals; 532
• DOI: 10.1002/fld.4499
• Languages: English
• ISSNs: 0271-2091

Abstract

Summary: In shallow water flow and transport modeling, the monotonic upstream‐centered scheme for conservation laws (MUSCL) is widely used to extend the original Godunov scheme to second‐order accuracy. The most important step in MUSCL‐type schemes is MUSCL reconstruction, which calculate‐extrapolates the values of independent variables from the cell center to the edge. The monotonicity of the scheme is preserved with the help of slope limiters that prevent the occurrence of new extrema during reconstruction. On structured grids, the calculation of the slope is straightforward and usually based on a 2‐point stencil that uses the cell centers of the neighbor cell and the so‐called far‐neighbor cell of the edge under consideration. On unstructured grids, the correct choice for the upwind slope becomes nontrivial. In this work, 2 novel total variation diminishing schemes are developed based on different techniques for calculating the upwind slope and the downwind slope. An additional treatment that stabilizes the scheme is discussed. The proposed techniques are compared to 2 existing MUSCL reconstruction techniques, and a detailed discussion of the results is given. It is shown that the proposed MUSCL reconstruction schemes obtain more accurate results with less numerical diffusion and higher efficiency. Abstract : This work presents two novel MUSCL schemes compared with two existing schemes, novel schemes can provide accuracy, stable and computational efficient result; 2D TVD conditions are discussed which give a new condition for preserving the TVD condition with searching for the far‐neighboring cell; an additional limitation for the MUSCL scheme is proposed to avoid the overestimated velocity at wetting and drying interface.