A time discretization method for the saturation-based formulation of Richards’ equation is analysed numerically. Depending on the flow regime, diffusion may vanish or blow up, leading to degeneracy in the original nonlinear parabolic problem. The numerical approach is based on an Euler implicit 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.
|Name||Lecture Notes in Applied and Computational Mechanics|
|Conference||conference; Stuttgart; 2002-04-08; 2002-04-10|
|Period||8/04/02 → 10/04/02|