Dec 3 – 4, 2020
Europe/Paris timezone

Convergence of an MPFA-O finite volume scheme for a seawater intrusion model with sharp-diffuse interfaces

Dec 4, 2020, 10:00 AM
Zoom (Virtuel)



Salle 1 : Salle 2 : Salle 3 :


Khadija Talali (LMAP, Universite de Pau et des Pays de l’Adour, CNRS, E2S UPPA, 64000 Pau, France---L2M3S-ENSAM, Université Moulay Ismaïl, Meknès 50000, Morocco)


We consider a sharp-diffuse interfaces seawater intrusion model [3] in
coastal aquifers. This process leads to a coupled system of two nonlinear parabolic partial differential equations simulating two immiscible fluids considering the dynamics of transition zones.

To the discretization in time, we apply a cell-centred Multi-Point Flux Approximation (MPFA-O) finite volume scheme [1] for the nonlinear system on an unstructured mesh and for the time discretization, we use an implicit Euler scheme, which allows the use of large time steps and then the reduction of CPU time. It is shown that this scheme ensures the non-negativity of the discrete solution of the freshwater and seawater thicknesses, taking into account the anisotropy and heterogeneity of the coastal aquifer. Based on an apriori estimate and using a fixed point theorem, we have established the existence of the discrete solution issued from the MPFA-O scheme. After that, we proved the strong convergence of the numerical solution to the weak solution of the continuous problem due to some recent compactness arguments.

We have developed and implemented a new module in the open-source platform DuMu$^X$[4]. The efficiency and the robustness of our method are studied numerically by comparing our results with an analytical solution presented in [2]. Moreover, the comparison between our $2D$ numerical simulations and those of [3] have shown that this approach yields performant results. Our new module based on recent numerical tools is accurate, efficient and able to solve a $2 D$ seawater intrusion model tacking into account the dynamics of the transition zones. Some numerical results will be presented. Future works will mainly focus on the validation of this approach for more realistic test cases.

 [1] I. Aavatsmark, "An introduction to multipoint flux approximations for quadrilateral grids, Computational Geosciences, vol. 6, no. 3-4, pp. 405--432, 2002.

 [2] Abudawia, A.; Mourad, A.; Rodrigues, J.H.; Rosier, C. A finite element method for a seawater intrusion problem in unconfined aquifers. Appl. Numer. Math. 2018, 127, 349--369

  [3] Choquet, C.; Diédhiou, M.; Rosier, C. Derivation of a sharp-diffuse interfaces model for seawater intrusion in a free aquifer, Numerical simulations. SIAM J. Appl. Math. 2016, 76, 138--158.

  [4] DuMu$^X$: Dune for multi- Phases, Component, Scale, Physics, ... flow and transport in porous media,, last accessed December 14, 2019.

Primary authors

Prof. Brahim Amaziane (Universite de Pau et des Pays de l’Adour, E2S UPPA, CNRS, LMAP, Pau, France) Prof. Mustapha El Ossmani (Université Moulay EMMACS-ENSAM Ismail, Meknès, Morocco) Khadija Talali (LMAP, Universite de Pau et des Pays de l’Adour, CNRS, E2S UPPA, 64000 Pau, France---L2M3S-ENSAM, Université Moulay Ismaïl, Meknès 50000, Morocco)

Presentation materials