Extension of WAF type methods to non-homogeneous Shallow Water Equations with pollutant

This paper deals with the extension of the WAF method to discretize Shallow Water Equations with pollutants. We consider two different versions of the WAF method, by approximating the intermediate waves using the flux of HLL or the direct approach of HLLC solver. It is seen that both versions can be written under the same form with different definitions for the approximation of the velocity waves. We also propose an extension of the method to non-homogeneous systems. In the case of homogeneous systems it is seen that we can rewrite the third component of the numerical flux in terms of an intermediate wave speed approximation. We conclude that – in order to have the same relation for non-homogeneous systems – the approximation of the intermediate wave speed must be modified. The proposed extension of the WAF method preserves all stationary solutions, up to second order accuracy, and water at rest in an exact way, even with arbitrary pollutant concentration. Finally, we perform several numerical tests, by comparing it with HLLC solver, reference solutions and analytical solutions.