An Hermitian finite element solution of the two-dimensional saturated-unsaturated flow equation

This paper describes a Galerkin-type finite element solution of the two-dimensional saturated-unsaturated flow equation. The numerical solution uses an incomplete (reduced) set of Hermitian cubic basis functions and is formulated in terms of normal and tangential coordinates. The formulation leads to continuous pressure gradients across interelement boundaries for a number of well-defined element configurations, such as for rectangular and circular elements. Other elements generally lead to discontinuous gradients; however, the gradients remain uniquely defined at the nodes. The method avoids calculation of second-order derivatives, yet retains many of the advantages associated with Hermitian elements. A nine-point Lobatto-type integration scheme is used to evaluate all local element integrals. This alternative scheme produces about the same accuracy as the usual 9- or 16-point Gaussian quadrature schemes, but is computationally more efficient.