The Henry semianalytical solution for saltwater intrusion with reduced dispersion

Abstract

[1] The Henry semianalytical solution for salt water intrusion is widely used for benchmarking density dependent flow codes. The method consists of replacing the stream function and the concentration by a double set of Fourier series. These series are truncated at a given order and the remaining coefficients are calculated by solving a highly nonlinear system of algebraic equations. The solution of this system is often subject to substantial numerical difficulties. Previous works succeeded to provide semianalytical solutions only for saltwater intrusion problems with unrealistic large amount of dispersion. In this work, different truncations for the Fourier series are tested and the Levenberg-Marquardt algorithm, which has a quadratic rate of convergence, is applied to calculate their coefficients. The obtained results provide semianalytical solutions for the Henry problem in the case of reduced dispersion coefficients and for two freshwater recharge values: the initial value suggested by Henry (1964) and the reduced one suggested by Simpson and Clement (2004). The developed semianalytical solutions are compared against numerical results obtained by using the method of lines and advanced spatial discretization schemes. The obtained semianalytical solutions improve considerably the worthiness of the Henry problem and therefore, they are more suitable for testing density dependent flow codes.Citation: Zidane, A., A. Younes, P. Huggenberger, and E. Zechner (2012), The Henry semianalytical solution for saltwater intrusion with reduced dispersion, Water Resour. Res., 48, W06533,

The Role of Tectonic Structures and Density-Driven Groundwater Flow for Salt Karst Formation

“…Both, the MPFA and DG use the same type of unknowns that give the advantage of gathering it into one system matrix. The combination DG-MPFA has shown to be a robust and accurate approach for modelling density driven flow problems (Zidane et al 2012). To model the salt dissolution process, we used a technique based on the variation of the size of the mesh.…”

Simulation of rock salt dissolution and its impact on land subsidence

“…On the other hand, both the MPFA and DG use the same type of unknowns that give the advantage of gathering it into one system matrix. So far, the combination DG-MPFA has shown to be a robust and accurate approach for modelling density-driven flow problems (Konz et al, 2008;Zechner et al, 2011;Zidane et al, 2012). Non-iterative time stepping is used in this work.…”

“…The solute mass conservation equation is written in terms of mass fraction as follows (Boufadel, 2000;Cardenas et al, 2007;Younes et al, 2009;Zidane et al, 2012):…”

Groundwater withdrawal in randomly heterogeneous coastal aquifers

“…Dentz et al (2006) and Bolster et al (2007) applied perturbation techniques to solve analytically the Henry problem for a range of (small and intermediate) values of the Péclet number, which characterizes the relative strength of convective and dispersive transport mechanisms. Zidane et al (2012) solved the Henry problem for realistic (small) diffusion coefficients. Fahs et al (2014) 20 presented a semi-analytical solution for a system characterized by a square porous cavity subject to different salt concentrations at its vertical walls.…”

“…The complexity of the problem typically prevents its solution via analytical or semi-analytical methods, with few notable exceptions (Henry, 1964;Dentz et 15 al., 2006;Bolster et al 2007;Zidane et al, 2012;Fahs et al 2014). Henry (1964) developed a semi-analytical solution for a variable-density diffusion problem in a (vertical) two-dimensional homogeneous and isotropic domain.…”