A finite element numerical model for unsaturated flow in layered soils

In the present study, a Galerkin finite element numerical model is developed to analyse moisture movement in layered soils. The model basically solves the Richards equation governing moisture flow through unsaturated zone using the constitutive relationships proposed by Van Genuchten. The model is applied to analyse the flow behaviour in a layered soil consisting of sand and clay during gravity drainage and infiltration. It has been observed that the movement of water in layered soils is a complex process which is influenced by the relative thickness and position of individual layers.