Well-balanced hybrid compact-WENO scheme for shallow water equations

Abstract

We investigate the performance of the high order well-balanced hybrid compact-weighted essentially non-oscillatory (WENO) finite difference scheme (Hybrid) for simulations of shallow water equations with source terms due to a non-flat bottom topography. The Hybrid scheme employs the nonlinear fifth order characteristic-wise WENO-Z finite difference scheme to capture high gradients and discontinuities in an essentially non-oscillatory manner, and the linear spectral-like sixth order compact finite difference scheme to resolve the fine scale structures in the smooth regions of the solution efficiently and accurately. The high order multi-resolution analysis is employed to identify the smoothness of the solution at each grid point. In this study, classical one- and two-dimensional simulations, including a long time two-dimensional dam-breaking problem with a non-flat bottom topography, are conducted to demonstrate the performance of the hybrid scheme in terms of the exact conservation property (C-property), good resolution and essentially non-oscillatory shock capturing of the smooth and discontinuous solutions respectively, and up to 2–3 times speedup factor over the well-balanced WENO-Z scheme.