The Henry problem: New semianalytical solution for velocity‐dependent dispersion. Issue 9 (28th September 2016)

• Record Type: Journal Article
• Authors: Fahs, Marwan; Ataie‐Ashtiani, Behzad; Younes, Anis; Simmons, Craig T.; Ackerer, Philippe
• Is Part Of: Water resources research. Volume 52:Issue 9(2016:Sep.)
• Page Start: 7382
• Page End: 7407
• Publication Date: 2016-09-28
• Subjects: density‐driven flow -- velocity‐dependent dispersion -- Henry problem -- semianalytical solution -- code benchmarking; Hydrology -- Periodicals; 333.91
• DOI: 10.1002/2016WR019288
• Languages: English
• ISSNs: 0043-1397

Abstract

Abstract: A new semianalytical solution is developed for the velocity‐dependent dispersion Henry problem using the Fourier‐Galerkin method (FG). The integral arising from the velocity‐dependent dispersion term is evaluated numerically using an accurate technique based on an adaptive scheme. Numerical integration and nonlinear dependence of the dispersion on the velocity render the semianalytical solution impractical. To alleviate this issue and to obtain the solution at affordable computational cost, a robust implementation for solving the nonlinear system arising from the FG method is developed. It allows for reducing the number of attempts of the iterative procedure and the computational cost by iteration. The accuracy of the semianalytical solution is assessed in terms of the truncation orders of the Fourier series. An appropriate algorithm based on the sensitivity of the solution to the number of Fourier modes is used to obtain the required truncation levels. The resulting Fourier series are used to analytically evaluate the position of the principal isochlors and metrics characterizing the saltwater wedge. They are also used to calculate longitudinal and transverse dispersive fluxes and to provide physical insight into the dispersion mechanisms within the mixing zone. The developed semianalytical solutions are compared against numerical solutions obtained using an in house code based on variant techniques for both space and time discretization. The comparison provides better confidence on the accuracy of both numerical and semianalytical results. It shows that the new solutions are highly sensitive to the approximation techniques used in the numerical code which highlights their benefits for code benchmarking. Key Points: Analytical solutions for density‐driven flow in porous media with velocity‐dependent dispersion are missing A new semianalytical solution is developed for the Henry problem with velocity‐dependent dispersion The semianalytical solution provides physical insights into the dispersion mechanisms in the mixing zone