We present a novel hp-adaptive space-time discontinuous Galerkin (hp-STDG) method for the numerical solution of the nonstationary Richards equation equipped with Dirichlet, Neumann and seepage face boundary conditions. The hp-STDG method presented in this paper is a generalization of a hp-STDG method which was developed for time dependent non-linear convective-diffusive problems.
We describe the method and the single ring experiment, and then we present a numerical experiment which clearly demonstrates the superiority of the hp-STDG method over a discontinuous Galerkin method based on a static fine mesh.