[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,
Purpose -The purpose of this paper is to develop an efficient non-iterative model combining advanced numerical methods for solving buoyancy-driven flow problems. Design/methodology/approach -The solution strategy is based on two independent numerical procedures. The Navier-Stokes equation is solved using the non-conforming Crouzeix-Raviart (CR) finite element method with an upstream approach for the non-linear convective term. The advectiondiffusion heat equation is solved using a combination of Discontinuous Galerkin (DG) and Multi-Point Flux Approximation (MPFA) methods. To reduce the computational time due to the coupling, the authors use a non-iterative time stepping scheme where the time step length is controlled by the temporal truncation error. Findings -Advanced numerical methods have been successfully combined to solve buoyancy-driven flow problems on unstructured triangular meshes. The accuracy of the results has been verified using three test problems: first, a synthetic problem for which the authors developed a semi-analytical solution; second, natural convection of air in a square cavity with different Rayleigh numbers (103-108); and third, a transient natural convection problem of low Prandtl fluid with horizontal temperature gradient in a rectangular cavity. Originality/value -The proposed model is the first to combine advanced numerical methods (CR, DG, MPFA) for buoyancy-driven flow problems. It is also the first to use a non-iterative time stepping scheme based on local truncation error control for such coupled problems. The developed semi analytical solution based on Fourier series is also novel.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.