The numerical entropy production (NEP) for shallow water equations (SWE) is discussed and implemented as a smoothness indicator. We consider SWE in three different dimensions, namely, one-dimensional, one-and-a-half-dimensional, and two-dimensional SWE. An existing numerical entropy scheme is reviewed and an alternative scheme is provided. We prove the properties of these two numerical entropy schemes relating to the entropy steady state and consistency with the entropy equality on smooth regions. Simulation results show that both schemes produce NEP with the same behaviour for detecting discontinuities of solutions and perform similarly as smoothness indicators. An implementation of the NEP for an adaptive numerical method is also demonstrated.

Standard compact scheme and upwinding compact scheme have high order accuracy and high resolution, but cannot capture the shock which is a discontinuity. This work developed a modified upwinding compact scheme which uses an effective shock detector to block compact scheme to cross the shock and a control function to mix the flux with WENO scheme near the shock. The new scheme makes the original compact scheme able to capture the shock sharply and, more importantly, keep high order accuracy and high resolution in the smooth area which is particularly important for shock boundary layer and shock acoustic interactions. Numerical results show the scheme is successful for 2-D Euler and 2-D Navier-Stokes solvers. The examples include 2-D incident shock, 2-D incident shock and boundary layer interaction. The scheme is robust, which does not involve case related parameters.