ANALYSIS AND APPROXIMATION OF A SCALAR CONSERVATION LAW WITH A FLUX FUNCTION WITH DISCONTINUOUS COEFFICIENTS

Abstract

We study here a model of conservative nonlinear conservation law with a flux function with discontinuous coefficients, namely the equation ut + (k(x)u(1 - u))x = 0. It is a particular entropy condition on the line of discontinuity of the coefficient k which ensures the uniqueness of the entropy solution. This condition is discussed and justified. On the other hand, we perform a numerical analysis of the problem. Two finite volume schemes, the Godunov scheme and the VFRoe-ncv scheme, are proposed to simulate the conservation law. They are compared with two finite volume methods classically used in an industrial context. Several tests confirm the good behavior of both new schemes, especially through the discontinuity of permeability k (whereas a loss of accuracy may be detected when industrial methods are performed). Moreover, a modified MUSCL method which accounts for stationary states is introduced.