An iterative scaling function procedure for solving scalar non-linear hyperbolic balance laws (2007.09290v1)

Abstract: The scaling of the exact solution of a hyperbolic balance law generates a family of scaled problems in which the source term does not depend on the current solution. These problems are used to construct a sequence of solutions whose limiting function solves the original hyperbolic problem. Thus this gives rise to an iterative procedure. Its convergence is demonstrated both theoretically and analytically. The analytical demonstration is in terms of a local in time convergence and existence theorem in the $L^2$ framework for the class of problems in which the source term $s(q)$ is bounded, with $s(0) = 0$, is locally Lipschitz and belongs to $C^2(\mathbb{R}) \cap H^1 (\mathbb{R}) $. A convex flux function, which is usual for existence and uniqueness for conservation laws, is also needed. For the numerical demonstration, a set of model equations is solved, where a conservative finite volume method using a low-dissipation flux is implemented in the iteration stages. The error against reference solutions is computed and compared with the accuracy of a conventional first order approach in order to assess the gaining in accuracy of the present procedure. Regarding the accuracy only a first order scheme is explored because the development of a useful procedure is of interest in this work, high-order accurate methods should increase the computational cost of the global procedure. Numerical tests show that the present approach is a feasible method of solution.