Abstract We propose and analyse a general framework for space-time finite element methods that is based on a first-order system least-squares reformulation of the thick parabolic obstacle problem. Discretizations based on simplicial and prismatic meshes are studied and we show a priori error estimates for both. Convergence rates are derived for sufficiently smooth solutions. Reliable a posteriori bounds are provided and used to steer an adaptive algorithm. Numerical experiments including a one-phase Stefan problem and an American option pricing problem are presented.
Carvajal et al. (Tue,) studied this question.