A 2D parallel algorithm using MPICH for calculation of ground water flux at evaporation from water table

Abstract

This paper discusses the calculation of distribution of hydraulic head H (x, y), [m] in flux of ground water, which is described by the known steady state groundwater flow equation with nonlinear boundary condition at unknown surface boundary. A special iterative algorithm for a numerical solution of the finite-difference approximation of the above-stated problem using known Richardson's over-relaxation method has been developed. The algorithm for parallel calculation is developed using message passing interface (MPI) for two dimensional case. It has been shown by calculation that above insertion with small coefficient of filtration, a local depression of water surface has been formed. This is the area without outflow of ground water. It may substantially cause the growth of soil salinization above surface of that depression due to evaporation of ground water, as dissolved salt stays in soil layer. Dependencies of the process on evaporation rate and depth of insertion are also discussed. An approximate 58% efficiency in time was obtained using four processes in cluster using message passing interface.