Numerical solution of the vertical infiltration problem in bounded profiles

There is presented a numerical solution of the one-dimensional infiltration problem in bounded profiles. The soil is assumed to have constant water diffusivity and linear dependence between the hydraulic conductivity and the water content. Then, the vertical infiltration problem is modeled as an initial boundary value problem for a diffusion equation. We combine the finite difference scheme for the time variable with the fundamental sequence method for the spatial variable. The derived numerical scheme is applied to both flooding and rainfall scenarios. The convergence of the numerical approximated solution to the analytical one justifies the applicability of the method.