No abstract

The accurate description of the flow of water in porous media is of the greatest importance due to its numerous applications in several areas (groundwater, soil mechanics, etc.). The nonlinear Richards equation is often used as the governing equation that describes this phenomenon and a large number of research studies aimed to solve it numerically. However, due to the nonlinearity of the constitutive expressions for permeability, it remains a challenging modeling problem. In this paper, the stationary form of Richards’ equation used in saturated soils is solved by two numerical methods: generalized finite differences, an emerging method that has been successfully applied to the transient case, and a finite element method, for benchmarking. The nonlinearity of the solution in both cases is handled using a Newtonian iteration. The comparative results show that a generalized finite difference iteration yields satisfactory results in a standard test problem with a singularity at the boundary.