An adaptive space-time method for nonlinear poroviscoelastic flows with discontinuous porosities

This paper is concerned with a space-time adaptive numerical method for instationary porous media flows with nonlinear interaction between porosity and pressure, with focus on problems with discontinuous initial porosities. A convergent method that yields computable error bounds is constructed by a combination of Picard iteration and a least-squares formulation. The adaptive scheme permits spatially variable time steps, which in numerical tests are shown to lead to efficient approximations of solutions with localized porosity waves. The method is also observed to exhibit optimal convergence with respect to the total number of spatio-temporal degrees of freedom.