One unconditionally stable CABARET class scheme for shallow water equations

Abstract

The paper considers a new approach to the construction of implicit unconditionally stable schemes within the framework of the CABARET balance-characteristic technique in relation to the system of shallow water equations. The method is based on the idea of inverting coordinate axes in the CABARET scheme to overcome the time step limitation. The system of equations is nonlinear, since the equations include minmax limiting operations based on the maximum principle for local Riemann invariants. This limitation significantly improves the dispersion properties of the numerical scheme. The nonlinear system of equations is solved using the marching order method. The paper presents the derivation of the numerical scheme for the Courant–Friedrichs–Lewy number CFL ≤ 1 and CFL > 1. Test calculations are presented on the simplest transport equation and one-dimensional shallow water problems for the subsonic case. Conclusions are drawn about the influence of nonlinear flux correction on the solution.