Error Estimates for a Time Discretization Method for the Richards’ Equation

We present a numerical analysis of an implicit time discretization method applied to Richards' equation. Written in its saturation-based form, this nonlinear parabolic equation models water flow into unsaturated porous media. Depending on the soil parameters, the diffusion coefficient may vanish or explode, leading to degeneracy in the original parabolic equation. The numerical approach is based on an implicit Euler time discretization scheme and includes a regularization step, combined with the Kirchhoff transform. Convergence is shown by obtaining error estimates in terms of the time step and of the regularization parameter