varrho_h>0$, i.e. $\Omega_h \subset \Omega$.

Let $u \in H^{k+1}(\Omega)$ denote the solution to (\[eq:poissoninterior\])–(\[eq:poissonbc\]). Let $u_h, \lambda_h \in V_h \times \Lambda_h$ denote the solution of (\[eq:mainproblem\]). Assume that the polynomial order of $V_h$ is $k \in \{1,2,3\}$, with enrichment on the boundary and $\Lambda_h \equiv X_h^{k-1}$. Assume that $V_h \times \Lambda_h$ satisfies (\[eq:infsup\]). Then there holds, with $\tilde \lambda = {\boldsymbol n}_h \cdot \nabla
u\vert_{\partial \Omega_h}$, $$\begin{aligned}
\label{eq:errorestenergy}
{|\mspace{-1mu}|\mspace{-1mu}|}(u - u_h, \tilde \lambda- \lambda_h) {|\mspace{-1mu}|\mspace{-1mu}|}&\lesssim 
h^{k} \|u \|_{H^{k+1}(\Omega)}
+ 
h^{-1/2} \delta_h^{2} 
\sup_{0\leq t \leq \delta_0} \| D^{k+1} u\|_{L^2(\partial \Omega_t)}
\\ \nonumber
&\qquad + 
\textcolor{black}{h^{1/2}} \delta_h^{l+1} 
\sup_{-\delta_0\leq t < 0} 
\| D_n^{l} (f + \Delta u)\|_{L^2(\partial \Omega_t)}.\end{aligned}$$

Let $\tilde u_h \in V_h$ denote the solution to (\[eq:Nitform\]). We recall from [@BrDuTh72; @BuHaLa18] that the following error bound holds $$\begin{aligned}
\nonumber
{|\mspace{-1mu}|\mspace{-1mu}|}(u - \tilde u_h,0) {|\mspace{-1mu}|\mspace{-1mu}|}&+\|h^{\frac12} {\boldsymbol n}_h \cdot \nabla (u -
  \tilde u_h)\|_{\partial \Omega_h}+ \|h^{-\frac12
varrho_h>0 $, i.e. $ \Omega_h \subset \Omega$. 

 Let $ u \in H^{k+1}(\Omega)$ denote the solution to (\[eq: poissoninterior\])–(\[eq: poissonbc\ ]). Let $ u_h, \lambda_h \in V_h \times \Lambda_h$ announce the solution of (\[eq: mainproblem\ ]). wear that the polynomial order of $ V_h$ is $ k \in \{1,2,3\}$, with enrichment on the boundary and $ \Lambda_h \equiv X_h^{k-1}$. Assume that $ V_h \times \Lambda_h$ satisfies (\[eq: infsup\ ]). Then there hold, with $ \tilde \lambda = { \boldsymbol n}_h \cdot \nabla 
 u\vert_{\partial \Omega_h}$, $ $ \begin{aligned } 
 \label{eq: errorestenergy } 
 { |\mspace{-1mu}|\mspace{-1mu}|}(u - u_h, \tilde \lambda- \lambda_h) { |\mspace{-1mu}|\mspace{-1mu}|}&\lesssim 
 h^{k } \|u \|_{H^{k+1}(\Omega) } 
 + 
 h^{-1/2 } \delta_h^{2 } 
 \sup_{0\leq t \leq \delta_0 } \| D^{k+1 } u\|_{L^2(\partial \Omega_t) } 
 \\ \nonumber 
 & \qquad + 
 \textcolor{black}{h^{1/2 } } \delta_h^{l+1 } 
 \sup_{-\delta_0\leq deoxythymidine monophosphate < 0 } 
 \| D_n^{l } (f + \Delta u)\|_{L^2(\partial \Omega_t)}.\end{aligned}$$ 

 Let $ \tilde u_h \in V_h$ denote the solution to (\[eq: Nitform\ ]). We remember from [ @BrDuTh72; @BuHaLa18 ] that the following error oblige holds $ $ \begin{aligned } 
 \nonumber 
 { |\mspace{-1mu}|\mspace{-1mu}|}(u - \tilde u_h,0) { |\mspace{-1mu}|\mspace{-1mu}|}&+\|h^{\frac12 } { \boldsymbol n}_h \cdot \nabla (u - 
   \tilde u_h)\|_{\partial \Omega_h}+ \|h^{-\frac12
vargho_h>0$, i.e. $\Omega_h \subset \Omtga$.

Let $u \in H^{k+1}(\Omeyq)$ denove the aolution to (\[eq:poissoninterior\])–(\[eq:poisslnvc\]). Leu $u_h, \lambda_h \in V_h \gimes \Lamhda_h$ denite uhe solution of (\[eq:mainpronjem\]). Zdsumz vhat the polynokial order of $V_h$ is $k \in \{1,2,3\}$, dich enrichment on the boundary and $\Laibda_h \eaulv X_h^{k-1}$. Assume jhat $N_h \timss \Lambda_h$ satisfies (\[eq:infsup\]). Then there iolds, with $\tildr \lambda = {\boldsymbol n}_h \cdlt \nwbla
u\vert_{\partial \Olega_h}$, $$\begin{qliggwd}
\label{eq:errurestenergj}
{|\kspace{-1mu}|\mskace{-1mu}|}(u - u_h, \tilde \lambda- \lambda_h) {|\mspaee{-1mu}|\mspace{-1mo}|}&\lssdvim 
h^{k} \|u \|_{H^{k+1}(\Imega)}
+ 
h^{-1/2} \delta_h^{2} 
\sui_{0\keq t \neq \delya_0} \| D^{k+1} u\|_{L^2(\partisl \Mmeta_t)}
\\ \nonumber
&\qquad + 
\teetcolor{black}{h^{1/2}} \delta_h^{j+1} 
\sup_{-\deltd_0\lzq t < 0} 
\| D_n^{l} (f + \Delta y)\|_{L^2(\partidl \Okega_g)}.\wnd{xlifnxd}$$

Lst $\tilfe n_h \in V_h$ dehote the sooution to (\[eq:Nitform\]). Wq recall from [@GrDuTh72; @BtHaLa18] that the following error bound honds $$\begin{aligned}
\nonumber
{|\mspace{-1mu}|\mspace{-1mu}|}(u - \tildg u_h,0) {|\mspacq{-1mu}|\mspace{-1mu}|}&+\|h^{\frac12} {\boldsymbol n}_h \cdot \nabla (u -
  \tilge u_h)\|_{\'aftico \Omeeq_h}+ \|h^{-\frac12
varrho_h>0$, i.e. $\Omega_h \subset \Omega$. Let $u denote solution to Let $u_h, \lambda_h the of (\[eq:mainproblem\]). Assume the polynomial order $V_h$ is $k \in \{1,2,3\}$, with on the boundary and $\Lambda_h \equiv X_h^{k-1}$. Assume that $V_h \times \Lambda_h$ satisfies Then there holds, with $\tilde \lambda = {\boldsymbol n}_h \cdot \nabla u\vert_{\partial \Omega_h}$, \label{eq:errorestenergy} - \tilde \lambda_h) {|\mspace{-1mu}|\mspace{-1mu}|}&\lesssim h^{k} \|u \|_{H^{k+1}(\Omega)} + h^{-1/2} \delta_h^{2} \sup_{0\leq t \leq \delta_0} \| D^{k+1} u\|_{L^2(\partial \Omega_t)} \nonumber &\qquad + \textcolor{black}{h^{1/2}} \delta_h^{l+1} \sup_{-\delta_0\leq t < \| D_n^{l} (f + u)\|_{L^2(\partial \Omega_t)}.\end{aligned}$$ Let $\tilde u_h V_h$ the solution (\[eq:Nitform\]). recall [@BrDuTh72; @BuHaLa18] that following error bound holds $$\begin{aligned} \nonumber {|\mspace{-1mu}|\mspace{-1mu}|}(u - \tilde u_h,0) {|\mspace{-1mu}|\mspace{-1mu}|}&+\|h^{\frac12} {\boldsymbol n}_h \cdot \nabla (u - u_h)\|_{\partial \Omega_h}+
varrho_h>0$, i.e. $\Omega_h \subset \OmegA$.

Let $u \in H^{k+1}(\OMega)$ dEnoTe tHe SoluTion To (\[eq:poissonintERior\])–(\[Eq:poissonbc\]). Let $u_h, \lambda_H \in V_h \TiMEs \LaMBdA_h$ denOte the sOLuTIOn oF (\[eQ:mAinPrOBlEm\]). AssUme That the Polynomial OrdEr Of $V_h$ is $k \in \{1,2,3\}$, witH EnRichment on The Boundary and $\LAmbDa_h \equIv x_h^{k-1}$. aSsume ThaT $V_h \tiMes \LamBDa_h$ satIsfies (\[eq:iNfSUp\]). Then THere holDS, WiTh $\tiLde \lambda = {\boldsymbOL n}_H \Cdot \nabla
u\vert_{\PartiaL \OMEgA_H}$, $$\BegIn{aLigned}
\labeL{eQ:erroREstenerGY}
{|\mSPACe{-1mU}|\Mspace{-1mu}|}(u - u_h, \tiLde \lambda- \laMBda_H) {|\mspacE{-1mU}|\msPAce{-1mu}|}&\lEsssiM 
h^{K} \|U \|_{H^{k+1}(\omega)}
+ 
h^{-1/2} \delta_H^{2} 
\sup_{0\Leq t \leq \deLta_0} \| D^{k+1} u\|_{l^2(\Partial \oMega_t)}
\\ \noNumber
&\QquAd + 
\tExtcOLoR{bLacK}{h^{1/2}} \DEltA_H^{l+1} 
\Sup_{-\DEltA_0\leq t < 0} 
\| D_n^{l} (F + \DElTa u)\|_{L^2(\pArtiAL \oMEga_t)}.\End{AligNed}$$

LeT $\tilde u_h \in V_h$ dEnoTe thE SolUtion To (\[eq:NItfoRm\]). we recAll froM [@BrDuth72; @buHaLa18] that the foLlowIng error bOunD hOldS $$\bEgin{aLIgned}
\nOnuMbeR
{|\mspace{-1Mu}|\mspacE{-1Mu}|}(u - \TiLDE U_h,0) {|\Mspace{-1mu}|\mspace{-1mu}|}&+\|h^{\fRaC12} {\BOlDsymbol n}_H \cdot \nABlA (u -
  \TIlde u_h)\|_{\paRtIal \omegA_H}+ \|H^{-\frac12
varrho_h>0$, i.e. $\Omega_ h \subset\Omeg a$. 

L et  $u\inH^{k+1}(\Omega ) $ de note the solution to ( \[eq: po i sson i nt erior \])–(\[ e q: p o iss on bc \]) .L et  $u_h , \ lambda_ h \in V_h\ti me s \Lambda_h$ de note the s olu tion of (\[e q:m ainpro bl em\ ] ). As sum e tha t thep olynom ial order  o f  $V_h$ is $k \ i n  \ {1,2 ,3\}$, with enric h me n t on the bound ary an d$ \L a m bda _h\equiv X_h ^{ k-1}$ .  Assume th a t $V_ h  \times \Lamb da_h$ satis f ies  (\[eq :i nfs u p\]).Thenth e reholds, with  $\t ilde \lam bda ={ \boldsy m bol n}_ h \cdo t \ nab la
u \ ve rt _{\ pa r tia l  \ Ome g a_h }$, $$\b eg in {alig ned} \ l a bel{ eq: erro reste nergy}
{|\msp ace {-1m u }|\ mspac e{-1m u}|} (u  - u_ h, \ti lde \ la mbda- \lambda_h ) {| \mspace{- 1mu }| \ms pa ce{-1 m u}|}&\ les ssi m 
h^{k } \|u \ | _{H ^{ k + 1 }( \Omega)}
+ 
h^{-1/ 2} \ de lta_h^{2 } 
\su p _{ 0\ l eq t \le q\de lta_ 0 }  \| D ^{k+ 1 }u\|_{L^2 (\part i al  \ Omega_t )} 
\\ \n on umb er&\qqu a d +
\text color{bl ack}{ h ^{1/2}} \delta _ h^{l+1} 
\sup _ {- \ d el t a_0\ leq  t < 0} 
\|  D_n ^ {l}(f + \D elt a  u)\| _{L^2 (\ p ar t ial \Omega_t)}.\end {a ligned }$$
Let $\tilde u _h \in V_h $ d enote th e so l ut i on to (\[eq:Ni tform \]). We re c all from  [@Br DuTh72;@BuHaLa18 ] that the  fo llo win g e r r or  bound holds$ $ \beg in {aligne d}\nonumb er{|\ msp ace {- 1mu}|\msp ace{-1mu }| }( u-\ti lde u _ h,0) {|\ ms pac e{ -1m u}|\m s pace{- 1mu}| }&+\ |h ^{ \ fra c12} {\ b ol d s ymbo ln} _h \ cdo t\nabl a (u -
 \tilde  u_h)\|_{ \pa r tial  \ Om ega_h}+  \|h^{-\frac1 2
varrho_h>0$, i.e._$\Omega_h \subset_\Omega$.

Let $u \in H^{k+1}(\Omega)$_denote the_solution_to (\[eq:poissoninterior\])–(\[eq:poissonbc\])._Let_$u_h, \lambda_h \in_V_h \times \Lambda_h$_denote the solution of_(\[eq:mainproblem\]). Assume that_the_polynomial order of $V_h$ is $k \in \{1,2,3\}$, with enrichment on the boundary and_$\Lambda_h_\equiv X_h^{k-1}$._Assume_that_$V_h \times \Lambda_h$ satisfies (\[eq:infsup\])._Then there holds, with $\tilde_\lambda =_{\boldsymbol n}_h \cdot \nabla
u\vert_{\partial \Omega_h}$, $$\begin{aligned}
\label{eq:errorestenergy}
{|\mspace{-1mu}|\mspace{-1mu}|}(u - u_h,_\tilde_\lambda- \lambda_h) {|\mspace{-1mu}|\mspace{-1mu}|}&\lesssim_
h^{k} \|u \|_{H^{k+1}(\Omega)}
+ 
h^{-1/2} \delta_h^{2} 
\sup_{0\leq t \leq \delta_0}_\| D^{k+1} u\|_{L^2(\partial \Omega_t)}
\\ \nonumber
&\qquad +_
\textcolor{black}{h^{1/2}} \delta_h^{l+1} 
\sup_{-\delta_0\leq_t_<_0} 
\| D_n^{l} (f_+ \Delta u)\|_{L^2(\partial \Omega_t)}.\end{aligned}$$

Let $\tilde u_h_\in V_h$ denote the solution to_(\[eq:Nitform\]). We recall from [@BrDuTh72; @BuHaLa18] that_the following error bound holds $$\begin{aligned}
\nonumber
{|\mspace{-1mu}|\mspace{-1mu}|}(u_- \tilde u_h,0) {|\mspace{-1mu}|\mspace{-1mu}|}&+\|h^{\frac12} {\boldsymbol_n}_h \cdot_\nabla (u -
  \tilde_u_h)\|_{\partial \Omega_h}+ \|h^{-\frac12