Numerical analysis of a finite volume scheme for a seawater intrusion model with cross-diffusion in an unconfined aquifer

We consider a degenerate parabolic system modeling the flow of fresh and saltwater in a porous medium in the context of seawater intrusion. We propose and analyze a finite volume scheme based on two-point flux approximation with upwind mobilities. The scheme preserves at the discrete level the main features of the continuous problem, namely the nonnegativity of the solutions, the decay of the energy and the control of the entropy and its dissipation. Based on these nonlinear stability results, we show that the scheme converges toward a weak solution to the problem. Numerical results are provided to illustrate the behavior of the model and of the scheme.