This paper studies a family of finite volume schemes for the hyperbolic scalar conservation law $u_t+{\nabla }_g·f(x,u)=0$ on a closed Riemannian manifold M. For an initial value in BV(M) we will show that these schemes converge with a $h^{\frac{1}{4}}$ convergence rate towards the entropy solution. When M is 1-dimensional the schemes are TVD and we will show that this improves the convergence rate to $h^{\frac{1}{2}}.$

We study a two-grid scheme fully discrete in time and space for solving the Navier-Stokes system. In the first step, the fully non-linear problem is discretized in space on a coarse grid with mesh-size H and time step k. In the second step, the problem is discretized in space on a fine grid with mesh-size h and the same time step, and linearized around the velocity uH computed in the first step. The two-grid strategy is motivated by the fact that under suitable assumptions, the contribution of uH...

We propose a 1D adaptive numerical scheme for hyperbolic conservation laws based on the numerical density of entropy production (the amount of violation of the theoretical entropy inequality). This density is used as an a posteriori error which provides information if the mesh should be refined in the regions where discontinuities occur or coarsened in the regions where the solution remains smooth. As due to the Courant-Friedrich-Levy stability condition the time step is restricted and leads to...

In this paper, we first construct a model for free surface flows that takes into account the air entrainment by a system of four partial differential equations. We derive it by taking averaged values of gas and fluid velocities on the cross surface flow in the Euler equations (incompressible for the fluid and compressible for the gas). The obtained system is conditionally hyperbolic. Then, we propose a mathematical kinetic interpretation of this system to finally construct a two-layer kinetic scheme...

We consider the system of partial differential equations governing the one-dimensional flow of two superposed immiscible layers of shallow water. The difficulty in this system comes from the coupling terms involving some derivatives of the unknowns that make the system nonconservative, and eventually nonhyperbolic. Due to these terms, a numerical scheme obtained by performing an arbitrary scheme to each layer, and using time-splitting or other similar techniques leads to instabilities in...

This study concerns some asymptotic models used to compute the flow outside and inside fractures in a bidimensional porous medium. The flow is governed by the Darcy law both in the fractures and in the porous matrix with large discontinuities in the permeability tensor. These fractures are supposed to have a small thickness with respect to the macroscopic length scale, so that we can asymptotically reduce them to immersed polygonal fault interfaces and the model finally consists in a coupling between...