A Semi-Analytical Solution for the Reactive Henry Saltwater Intrusion Problem

Because of the existence of a semi-analytical solution, the Henry saltwater intrusion problem has been widely used for benchmarking non-reactive density-driven flow models. In this work, we extend the semi-analytical solution of Henry to reactive transport in variable-density fluid flow. Accurate semi-analytical solutions are provided for three test cases dealing with saltwater transport including dissolution and degradation reactions. About 6,195 terms are required in the Fourier series to obtain a stable solution for these test cases instead of the 78 initially used by Henry (Sea Water in Coastal Aquifers 1613-C:70–84, 1964) for the non-reactive problem. The resolution of the highly non-linear system is made possible due to the modified Powell hybrid algorithm with an analytical evaluation of the Jacobian. Numerical simulations are performed using different numerical methods and grid sizes to evaluate the benefits of these new test cases for benchmarking reactive density-driven flow models.