paper

Improved rates for a space-time FOSLS of parabolic PDEs

arXiv:2208.10824 · doi:10.1007/s00211-023-01387-3

Abstract

We consider the first-order system space-time formulation of the heat equation introduced in [Bochev, Gunzburger, Springer, New York (2009)], and analyzed in [Führer, Karkulik, Comput. Math. Appl. 92 (2021)] and [Gantner, Stevenson, ESAIM Math. Model. Numer. Anal.} 55 (2021)], with solution components . The corresponding operator is boundedly invertible between a Hilbert space and a Cartesian product of -type spaces, which facilitates easy first-order system least-squares (FOSLS) discretizations. Besides -norms of and , the (graph) norm of contains the -norm of . When applying standard finite elements w.r.t. simplicial partitions of the space-time cylinder, estimates of the approximation error w.r.t. the latter norm require higher-order smoothness of . In experiments for both uniform and adaptively refined partitions, this manifested itself in disappointingly low convergence rates for non-smooth solutions . In this paper, we construct finite element spaces w.r.t. prismatic partitions. They come with a quasi-interpolant that satisfies a near commuting diagram in the sense that, apart from some harmless term, the aforementioned error depends exclusively on the smoothness of , i.e., of the forcing term . Numerical results show significantly improved convergence rates.

References in corpus (5)

Cited by in corpus (6)