Well-Balanced High-Order Discontinuous Galerkin Methods for Systems of Balance Laws

This work introduces a general strategy to develop well-balanced high-order Discontinuous Galerkin (DG) numerical schemes for systems of balance laws. The essence of our approach is a local projection step that guarantees the exactly well-balanced character of the resulting numerical method for smooth stationary solutions. The strategy can be adapted to some well-known different time marching DG discretisations. Particularly, in this article, Runge–Kutta DG and ADER DG methods are studied. Additionally, a limiting procedure based on a modified WENO approach is described to deal with the spurious oscillations generated in the presence of non-smooth solutions, keeping the well-balanced properties of the scheme intact. The resulting numerical method is then exactly well-balanced and high-order in space and time for smooth solutions. Finally, some numerical results are depicted using different systems of balance laws to show the performance of the introduced numerical strategy.

Existence and uniqueness of generalized solutions to hyperbolic systems with linear fluxes and stiff sources

Motivated by the modeling of boiling two-phase flows, we study systems of balance laws with a source term defined as a discontinuous function of the unknown. Due to this discontinuous source term, the classical theory of partial differential equations (PDEs) is not sufficient here. Restricting to a simpler system with linear fluxes, a notion of generalized solution is developed. An important point in the construction of a solution is that the curve along which the source jumps, which we call the boiling curve, must never be tangent to the characteristics. This leads to exhibit sufficient conditions which ensure the existence and uniqueness of a solution in two different situations: first when the initial data is smooth and such that the boiling curve is either overcharacteristic or subcharacteristic; then with discontinuous initial data in the case of Riemann problems. A numerical illustration is given in this last case.

Collocation Methods for High-Order Well-Balanced Methods for Systems of Balance Laws

In some previous works, two of the authors introduced a technique to design high-order numerical methods for one-dimensional balance laws that preserve all their stationary solutions. The basis of these methods is a well-balanced reconstruction operator. Moreover, they introduced a procedure to modify any standard reconstruction operator, like MUSCL, ENO, CWENO, etc., in order to be well-balanced. This strategy involves a non-linear problem at every cell at every time step that consists in finding the stationary solution whose average is the given cell value. In a recent paper, a fully well-balanced method is presented where the non-linear problems to be solved in the reconstruction procedure are interpreted as control problems. The goal of this paper is to introduce a new technique to solve these local non-linear problems based on the application of the collocation RK methods. Special care is put to analyze the effects of computing the averages and the source terms using quadrature formulas. A general technique which allows us to deal with resonant problems is also introduced. To check the efficiency of the methods and their well-balance property, they have been applied to a number of tests, ranging from easy academic systems of balance laws consisting of Burgers equation with some non-linear source terms to the shallow water equations—without and with Manning friction—or Euler equations of gas dynamics with gravity effects.

Convergence rate from systems of balance laws to isotropic parabolic systems, a periodic case

It is proved that partially dissipative hyperbolic systems converge globally-in-time to parabolic systems in a slow time scaling, when initial data are smooth and sufficiently close to constant equilibrium states. Based on this result, we establish the global-in-time error estimates between the smooth solutions to the partially dissipative hyperbolic systems and those to the isotropic parabolic limiting systems in a three dimensional torus, rather than in the one dimensional whole space (Appl. Anal. 100(5) (2021) 1079–1095). This avoids the condition raised for the strong connection between the flux and the source term and make the result obtained more generalized. In the proof, we provide a similar stream function technique which is valid for the three dimensional periodic case. Similar method is provided for the one-dimensional periodic case. As applications of the results, we give several examples arising from physical models at the end of the paper.

Fast Convergence Methods for Hyperbolic Systems of Balance Laws with Riemann Conditions

In this paper, we develop an accurate technique via the use of the Adomian decomposition method (ADM) to solve analytically a 2 × 2 systems of partial differential equation that represent balance laws of hyperbolic-elliptic type. We prove that the sequence of iteration obtained by ADM converges strongly to the exact solution by establishing a construction of fixed points. For comparison purposes, we also use the Sinc function methodology to establish a new procedure to solve numerically the same system. It is shown that approximation by Sinc function converges to the exact solution exponentially, also handles changes in type. A numerical example is presented to demonstrate the theoretical results. It is noted that the two methods show the symmetry in the approximate solution. The results obtained by both methods reveal that they are reliable and convenient for solving balance laws where the initial conditions are of the Riemann type.