Adaptive time discretization and linearization based on a posteriori estimates for the Richards equation

Abstract : We derive some a posteriori error estimates for the Richards equation, based on the dual norm of the residual. This equation is nonlinear in space and in time, thus its resolution requires fixed-point iterations within each time step. We propose a strategy to decrease the computational cost relying on a splitting of the error terms in three parts: linearization, time discretization, and space discretization. In practice, we stop the fixed-point iterations after the linearization error becomes negligible, and choose the time step in order to balance the time and space errors.