On the Accuracy of Calculating Invariants in Centered Rarefaction Waves and in Their Influence Area

We perform a comparative analysis of the accuracy of second-order TVD (Total Variation Diminishing), third-order RBM (Rusanov–Burstein–Mirin), and fifth-order in space and third-order in time A-WENO (Alternative Weighted Essentially Non-Oscillatory) difference schemes for solving a special Cauchy problem for shallow water equations with discontinuous initial data. The exact solution of this problem contains a centered rarefaction wave, but does not contain a shock wave. It is shown that in the centered rarefaction wave and its influence area, the solutions of these three schemes converge with different orders to different invariants of the exact solution. This leads to a decrease in the accuracy of these schemes when they used to calculate the vector of base variables of the considered Cauchy problem. The P-form of the first differential approximation of the difference schemes is used for theoretical justification of these numerical results.