On the reservoir technique convergence for nonlinear hyperbolic conservation laws. I (Q1025026)

Continuing a series of papers on the numerical solution of hyperbolic systems by this relatively new method (which basically consists in updating the solution locally, for every cell, only in that moment where the accumulated local Courant number equals 1, for this needing about three times as much memory but dramatically reducing numerical diffusion), the authors consider here its convergence in the one-dimensional case, concentrating on explicit conservative three-point schemes. They prove that the method, applied to problems with strictly convex flux and without interaction of shock and rarefaction waves, converges in \(L^1\) to the exact solution, moreover, the error is estimated by the initial \(L^1\)-error plus const\,\(t_{\Delta_t}\). In numerical experiments, the 1D Burgers equation and two linear hyperbolic systems of 3 equations are considered (all with discontinuous initial data), the systems with the aim to illustrate a further theoretical result: that the discrete solution at regular time intervals is exact on the discrete level (where const\,\(t_{\Delta_t}\) is reduced to zero), for rational eigenvalues.