Fourier series solution for an anisotropic and layered configuration of the dispersive Henry Problem

¹ Laboratoire d’Hydrologie et Géochimie de Strasbourg, University of Strasbourg/EOST/ENGEES CNRS, Strasbourg, France
² Department of Civil Engineering, Sharif University of Technology, Tehran, Iran
³ National Centre for Groundwater Research & Training and College of Science & Engineering, Flinders University, Adelaide, Australia

Contact Information: Marwan FAHS, Laboratoire d’Hydrologie et Géochimie de Strasbourg, University of Strasbourg/EOST/ENGEES, CNRS, 1 rue Blessig 67084 Strasbourg, France, Phone : +33-3-68 85 04 48, Email : fahs@unistra.fr

Abstract

Henry Problem (HP) still plays an important role in benchmarking numerical models of seawater intrusion (SWI) as well as being applied to practical and managerial purposes. The popularity of this problem is due to having a closed-form semi-analytical (SA) solution. The early SA solutions obtained for HP were limited to extensive assumptions that restrict its application in practical works. Several further studies expended the generality of the solution by assuming lower diffusion coefficients or including velocity-dependent dispersion in the results. However, all these studies are limited to homogeneous and isotropic domains. The present work made an attempt to improve the reality of the SA solution obtained for dispersive HP by considering anisotropic and stratified heterogeneous coastal aquifers. The solution is obtained by defining Fourier series for both stream function and salt concentration, applying a Galerkin treatment using the Fourier modes as trial functions and solving the flow and the salt transport equations simultaneously in the spectral space. In order to include stratified heterogeneity, a special depth-hydraulic conductivity model is applied that can be solved analytically without significant mathematical complexity. Several examples are proposed and studied. The results show excellent agreement between the SA and numerical solutions obtained with an in-house advanced finite element code.