We consider the second-order pseudo parabolic equation of finding $u:[0,T]\mapsto H_0^1(\mathrm{\Omega })$, $\mathrm{\Omega }\subset {\mathbb{R}}^d$, such that $$\begin{array}{c}\mathcal{L}{\partial }_{t}u(t)+\mathcal{M}u(t)=F(t)\text{in }(0,T]\end{array}$$ with two second order, elliptic operators $\mathcal{L},\mathcal{M}:H_0^1(\mathrm{\Omega })\mapsto H^{-1}(\mathrm{\Omega })$ and a source function $F:[0,T]\mapsto H^{-1}(\mathrm{\Omega })$. Furthermore an initial condition $$\begin{array}{c}u(0)=u_0,u_0\in H_0^1(\mathrm{\Omega }),\end{array}$$ is given. We will assume, that the operator $\mathcal{M}$ is bounded and that $\mathcal{L}$ is both bounded and coercive.

A computable a posteriori error bound in a weighted $H^1$-norm for a full discretisation using the backward differential formula of order two (BDF-2 method) in time and ${\mathbb{P}}_2$-elements in space is derived. To do so, we leverage the $C_0$ semigroup, generated by the operator ${\mathcal{L}}^{-1}\mathcal{M}$. Furthermore we adapt elliptic reconstructions introduced by C. Makridakis and R. N. Nochetto to pseudo parabolic problems.
Given some numerical results we analyze the estimate’s order, efficiency and components. Last but not least, we show that we can apply our a posteriori error bound to other time discretisations like the backward Euler and Crank Nicolson method.

References

1. Ch. Makridakis and R. H. Nochetto. Elliptic reconstruction and a posteriori error estimates for parabolic problems. SIAM J. Numer. Anal., 41(4):1585–1594, 2003.

[link to pdf] [back to Numdiff-17]