Regional Quasi Three Dimensional Unsaturated Saturated Water Flow Model Based on a Vertical Horizontal Splitting Concept
Due to the high nonlinearity of the three dimensional (3 D) unsaturated saturated water flow equation, using a fully 3 D numerical model is computationally expensive for large scale applications. A new unsaturated saturated water flow model is developed in this paper based on the vertical/horizontal splitting (VHS) concept to split the 3 D unsaturated saturated Richards' equation into a two dimensional (2 D) horizontal equation and a one dimensional (1 D) vertical equation. The horizontal plane of average head gradient in the triangular prism element is derived to split the 3 D equation into the 2 D equation. The lateral flow in the horizontal plane of average head gradient represented by the 2 D equation is then calculated by the water balance method. The 1 D vertical equation is discretized by the finite difference method. The two equations are solved simultaneously by coupling them into a unified nonlinear system with a single matrix. Three synthetic cases are used to evaluate the developed model code by comparing the modeling results with those of Hydrus1D, SWMS2D and FEFLOW. We further apply the model to regional scale modeling to simulate groundwater table fluctuations for assessing the model applicability in complex conditions. The proposed modeling method is found to be accurate with respect to measurements.