\\
 \leq & \|\pi_{k-1} \nabla \tilde u_h- \nabla \tilde u_h\|_{\Omega_h}.\end{aligned}$$ To bound the term in the right hand side we add and subtract $\nabla u - \pi_{k-1} \nabla u$ and use the triangle inequality and the stability of the $L^2$-projection $\pi_{k-1}$ to obtain $$\begin{aligned}
\|\pi_{k-1} \nabla \tilde u_h- \nabla \tilde u_h\|_{\Omega_h} &\leq
\|\pi_{k-1} (\nabla \tilde u_h- \nabla u)\|_{\Omega_h}+
  \|\pi_{k-1} \nabla u- \nabla u\|_{\Omega_h}+\|\nabla u- \nabla
  \tilde u_h\|_{\Omega_h} 
  \\
& \leq \|\pi_{k-1} \nabla u- \nabla u\|_{\Omega_h}+2 \|\nabla u- \nabla
  \tilde u_h\|_{\Omega_h}.\end{aligned}$$ For the first term in the right hand side we have the approximation bound $$\|\pi_{k-1} \nabla u- \nabla u\|_{\Omega_h} \lesssim h^k \|D^{k+1} u\|_{\Omega_h}.$$ The second term is bounded by (\[eq:errorNit\]). We conclude by applying the second inequality of Proposition \[prop:infsup\].

Remarks on cut finite element methods
=====================================

In the context of cut finite element methods the discontinuous multiplier spaces used above can no longer be expected to be stable. It is possible to stabilise the multiplier using Barbosa-Hughes stabilisation. However, fluctuation based multipliers are unlikely to be suitable in this context since the weak consistency of the fluctuations of the multiplier between elements depends on the geometry approximation through the interface normal. Since the method is of interest when the geometry approximation is of relatively low order, this limits the possibility to use fluctuation based stabilisation.

For closed smooth boundaries, one may prove inf-sup stability and optimal convergence, without stabilisation, when using continuous approximation of polynomial order less than or equal to $2$, for both the bulk variable and the multiplier provided $\rho_h = O(h^2)$. The approximation order of the interface normal, which is $O(h)$ prohibits higher order convergence
\\ 
  \leq & \|\pi_{k-1 } \nabla \tilde u_h- \nabla \tilde u_h\|_{\Omega_h}.\end{aligned}$$ To bound the term in the right hand slope we lend and subtract $ \nabla u - \pi_{k-1 } \nabla u$ and use the triangle inequality and the constancy of the $ L^2$-projection $ \pi_{k-1}$ to obtain $ $ \begin{aligned } 
 \|\pi_{k-1 } \nabla \tilde u_h- \nabla \tilde u_h\|_{\Omega_h } & \leq 
 \|\pi_{k-1 } (\nabla \tilde u_h- \nabla u)\|_{\Omega_h}+ 
   \|\pi_{k-1 } \nabla u- \nabla u\|_{\Omega_h}+\|\nabla u- \nabla 
   \tilde u_h\|_{\Omega_h } 
   \\ 
 & \leq \|\pi_{k-1 } \nabla u- \nabla u\|_{\Omega_h}+2 \|\nabla u- \nabla 
   \tilde u_h\|_{\Omega_h}.\end{aligned}$$ For the first term in the good hand side we suffer the approximation bound $ $ \|\pi_{k-1 } \nabla u- \nabla u\|_{\Omega_h } \lesssim h^k \|D^{k+1 } u\|_{\Omega_h}.$$ The second terminus is bounded by (\[eq: errorNit\ ]). We conclude by apply the second inequality of Proposition \[prop: infsup\ ]. 

 Remarks on cut finite element methods 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In the context of cut finite element methods the discontinuous multiplier quad used above can no longer be expected to be stable. It is potential to stabilise the multiplier using Barbosa - Hughes stabilization. However, fluctuation establish multipliers are unlikely to be suitable in this context since the weak consistency of the fluctuations of the multiplier between element depends on the geometry approximation through the interface normal. Since the method is of interest when the geometry approximation is of relatively low order, this limits the possibility to use fluctuation based stabilisation. 

 For closed smooth boundaries, one may prove inf - swallow constancy and optimal convergence, without stabilisation, when use continuous approximation of polynomial order less than or equal to $ 2 $, for both the bulk variable star and the multiplier provided $ \rho_h = O(h^2)$. The estimate order of the interface normal, which is $ O(h)$ prohibits higher order overlap
\\
 \lee & \|\pi_{k-1} \nabla \tilde u_h- \nanla \tilde u_h\|_{\Omegc_y}.\end{almgned}$$ Tk bound ghe term in the right hand smde qe ade and subtract $\nabla u - \pi_{k-1} \nabpa u$ and use rhe triangle inequality zkd thz wtability of tme $L^2$-projecthon $\pi_{k-1}$ to obtdiv $$\yegin{aligned}
\|\pi_{k-1} \nabla \tilde u_h- \nabla \tilde i_h\|_{\Lmega_h} &\leq
\|\pi_{k-1} (\gabls \tilsv m_h- \nabla u)\|_{\Omega_h}+
  \|\pi_{k-1} \nabla u- \nagla u\|_{\Omtga_h}+\|\nabla u- \nabla
  \tilde u_h\|_{\Omega_h} 
  \\
& \leq \|\pi_{k-1} \nabpa u- \nabla u\|_{\Omega_h}+2 \|\nabla u- \navla
  \rilde u_h\|_{\Omegx_h}.\end{aligntd}$$ For the fjrst term in the right hand sidd we kave the apkxixildtion bound $$\|\pi_{k-1} \gabla u- \nabla u\|_{\Omega_v} \lesssom h^k \|D^{k+1} u\|_{\Omegs_h}.$$ Vhe wecond term is boundev by (\[eq:errorNit\]). We cjnclude bf cpplying the second ibewualijy of Prooisigioh \['rol:infsuo\].

Rejarks on cht finite eoement methods
=====================================

In tht cjbtext of cut rinite ejement methods the discontinuous multipnied spaces used above can no longer be expectef to be seable. It is possible to stabilise the multiplier gsing Caryisa-Hueyed stabilisation. However, fluctuation based mulejpkivrs are unlikely bo be suitable in yhls sontext since the wzzk consistency of thf fluctoationw of the iultoplier between elements depwnds on the yeonetry approximatiou through thz intetface mormal. Since the method is kf interest when the eeometry approxioatpon hs of relatively low order, this limmts tke possicilijy to ufe fluctuahion ndsed stabilisation.

Vor cpoved smooth boundaries, one may prove inf-su' stability amd opnimal conrergenge, without stabylisation, when using eontinjous approqimation mf polynomiwl order less jhan or equal to $2$, for boty thw bulk xxriable and thr multiplpex providee $\rho_h = O(h^2)$. The appvoximxfion order of tkt ibterface normal, whkch id $P(h)$ [rohibits hicher orarr covvergence
\\ \leq & \|\pi_{k-1} \nabla \tilde u_h- u_h\|_{\Omega_h}.\end{aligned}$$ bound the in the right subtract u - \pi_{k-1} u$ and use triangle inequality and the stability of $L^2$-projection $\pi_{k-1}$ to obtain $$\begin{aligned} \|\pi_{k-1} \nabla \tilde u_h- \nabla \tilde u_h\|_{\Omega_h} &\leq (\nabla \tilde u_h- \nabla u)\|_{\Omega_h}+ \|\pi_{k-1} \nabla u- \nabla u\|_{\Omega_h}+\|\nabla u- \nabla \tilde \\ \leq \nabla \nabla u\|_{\Omega_h}+2 \|\nabla u- \nabla \tilde u_h\|_{\Omega_h}.\end{aligned}$$ For the first term in the right hand side have the approximation bound $$\|\pi_{k-1} \nabla u- \nabla \lesssim h^k \|D^{k+1} u\|_{\Omega_h}.$$ second term is bounded by We by applying second of \[prop:infsup\]. Remarks on finite element methods ===================================== In the context of cut finite element methods the discontinuous multiplier spaces used can no expected to stable. is to stabilise the Barbosa-Hughes stabilisation. However, fluctuation based multipliers be suitable in this context since the weak of the of the multiplier between elements depends the geometry approximation through the interface normal. Since method is of interest when the geometry approximation is of relatively low order, this limits to use fluctuation based For closed smooth one prove stability optimal convergence, stabilisation, when using continuous approximation of polynomial order less than or to $2$, for both the bulk variable and the multiplier = The approximation order the interface normal, which $O(h)$ higher order convergence
\\
 \leq & \|\pi_{k-1} \nabla \tilde u_h- \nabla \tiLde u_h\|_{\Omega_H}.\end{aLigNed}$$ to BounD the Term in the right HAnd sIde we add and subtract $\nabLa u - \pi_{K-1} \nABla u$ ANd Use thE triangLE iNEQuaLiTy And ThE StAbiliTy oF the $L^2$-prOjection $\pi_{K-1}$ to ObTain $$\begin{aliGNeD}
\|\pi_{k-1} \nabla \tIldE u_h- \nabla \tildE u_h\|_{\omega_h} &\LeQ
\|\pi_{K-1} (\Nabla \TilDe u_h- \nAbla u)\|_{\OMEga_h}+
  \|\pi_{K-1} \nabla u- \naBlA U\|_{\Omega_H}+\|\Nabla u- \nABLa
  \TildE u_h\|_{\Omega_h} 
  \\
& \leq \|\pi_{k-1} \naBLa U- \Nabla u\|_{\Omega_h}+2 \|\naBla u- \naBlA
  \TiLDE u_h\|_{\omeGa_h}.\end{aligNeD}$$ For tHE first tERm IN THe rIGht hand side we Have the apprOXimAtion bOuNd $$\|\pI_{K-1} \nabla U- \nablA u\|_{\oMegA_h} \lesssim h^k \|d^{k+1} u\|_{\OMega_h}.$$ The sEcond tERm is bouNDed by (\[eq:ErrorNIt\]). WE coNcluDE bY aPplYiNG thE SeConD IneQuality oF PRoPositIon \[pROP:INfsuP\].

ReMarkS on cuT finite elemenT meThodS
=====================================

in tHe conText oF cut FiNite eLement MethoDs The discontinuouS mulTiplier spAceS uSed AbOve caN No longEr bE exPected tO be stabLE. It Is POSSiBle to stabilise the mUlTIPlIer using barbosA-huGhES stabiliSaTioN. HowEVEr, fluCtuaTIoN based muLtipliERs ArE unlikeLy To be suItAblE in This cONtexT since The weak cOnsisTEncy of the fluctUAtions of the muLTiPLIeR BetwEen Elements depEnds ON the GeomETrY apPRoximAtion ThROuGH the interface normal. siNce the MethoD is of interest When the geoMETRy approxImatIOn IS of relatively lOw ordEr, this limiTS the possIbiliTy to use fLuctuatioN BAsed stabIliSatIon.

for CLOsEd smooth boundARIes, oNe May provE inF-sup staBilIty And OptImAl convergEnce, withOuT sTaBiLisAtion, WHen using CoNtiNuOus ApproXImatioN of poLynoMiAl ORdeR less thAN oR EQual To $2$, FoR botH thE bUlk vaRiabLE anD the mulTiplier prOviDEd $\rhO_h = o(h^2)$. the apprOximation ordeR oF the interfAcE noRmal, whICH is $O(h)$ proHibits higher order converGEnce
\\
 \leq & \|\pi_{k-1} \na bla \tilde  u_h-  \n abl a\til de u _h\|_{\Omega_h } .\en d{aligned}$$ To boundthe t er m  int he  righ t hands id e wead dand  s u bt ract$\n abla u- \pi_{k-1 } \ na bla u$ and u s ethe triang leinequality a ndthe st ab ili t y ofthe  $L^2 $-proj e ction$\pi_{k-1 }$ to obt a in $$\b e g in {ali gned}
\|\pi_{k-1} \n a bla \tilde u_h - \nab la \t i l deu_h \|_{\Omega _h } &\l e q
\|\pi _ {k - 1 }  (\ n abla \tilde u _h- \nablau )\| _{\Ome ga _h} + 
  \|\ pi_{k -1 }  \n abla u- \na blau\|_{\Ome ga_h}+ \ |\nabla u- \nab la
  \ til deu_h\ | _{ \O meg a_ h }  \ \
& \le q \|\pi_ {k -1 } \na blau - \ nabl a u \|_{ \Omeg a_h}+2 \|\nab lau- \ n abl a
  \ tilde  u_h \| _{\Om ega_h} .\end {a ligned}$$ For t he f irst term  in  t heri ght h a nd sid e w e h ave the  approx i mat io n b ou nd $$\|\pi_{k-1} \ na b l au- \nabl a u\|_ { \O me g a_h} \le ss sim  h^k \ |D^{k +1}u \| _{\Omega _h}.$$ Th esecondte rm isbo und edby (\ [ eq:e rrorNi t\]). We  conc l ude by applyin g  the second i n eq u a li t y of  Pr oposition \ [pro p :inf sup\ ] .
Re m arkson cu tf in i te element methods== ====== ===== ============= ========== = In the c onte x to f cut finite e lemen t methodst he disco ntinu ous mult iplier sp a c es usedabo vecan  no l on ger be expect e d  tobe  stable . I t is po ssi ble  to  st ab ilise the  multipl ie rus in g B arbos a -Hughesst abi li sat ion.H owever , flu ctua ti on bas ed mult i pl i e rs a re  u nlik ely  t o besuit a ble  in thi s context  si n ce t he  w eak con sistency of t he  fluctuati on s o f them u ltiplier  between elements depen d s on th e g eomet ry a pproximat ion  throu ght he int erface  norm al . S i n ce th e me tho dis of inte r e stwhenth e ge ometryapproximation is o f  re latively loword er,t h is  li m it s  th ep oss i b ility to use fl uctuationba s ed  stabilisa t ion .
For cl osed sm oothb oundari es, one m ay provein f-su p sta bility and  optimal  converge n ce, w i th out s tab ilisat io n,whenusingc ont inuou s appr ox imatio n ofpo lynomial  order less than or equ al to$2$,for  both the  bu l k v ariable a nd t he multipl ier  pr ovide d $ \ rho_h  = O ( h^ 2)$ .  Theappr o ximationo rd ero f  t he interfac e n orm al, w hic h  is $O (h)$  prohibits higher order converge nce
\\
 \leq_& \|\pi_{k-1}_\nabla \tilde u_h- \nabla_\tilde u_h\|_{\Omega_h}.\end{aligned}$$_To_bound the_term_in the right_hand side we_add and subtract $\nabla_u - \pi_{k-1}_\nabla_u$ and use the triangle inequality and the stability of the $L^2$-projection $\pi_{k-1}$ to_obtain_$$\begin{aligned}
\|\pi_{k-1} \nabla_\tilde_u_h-_\nabla \tilde u_h\|_{\Omega_h} &\leq
\|\pi_{k-1} (\nabla_\tilde u_h- \nabla u)\|_{\Omega_h}+
 _\|\pi_{k-1} \nabla_u- \nabla u\|_{\Omega_h}+\|\nabla u- \nabla
  \tilde u_h\|_{\Omega_h}_
_ \\
& \leq_\|\pi_{k-1} \nabla u- \nabla u\|_{\Omega_h}+2 \|\nabla u- \nabla
 _\tilde u_h\|_{\Omega_h}.\end{aligned}$$ For the first term_in the right_hand_side_we have the approximation_bound $$\|\pi_{k-1} \nabla u- \nabla u\|_{\Omega_h}_\lesssim h^k \|D^{k+1} u\|_{\Omega_h}.$$ The second_term is bounded by (\[eq:errorNit\]). We conclude_by applying the second inequality of_Proposition \[prop:infsup\].

Remarks on cut finite_element methods
=====================================

In_the context of cut finite_element methods the_discontinuous multiplier_spaces used above_can no longer be expected to_be stable. It_is possible to stabilise the multiplier_using_Barbosa-Hughes stabilisation. However,_fluctuation_based_multipliers are_unlikely to be_suitable_in this_context_since the weak consistency of the_fluctuations_of the multiplier between elements depends on_the geometry approximation through_the_interface normal. Since the_method is of interest when_the geometry approximation is of relatively_low order,_this limits_the possibility to use fluctuation based stabilisation.

For closed smooth boundaries, one_may prove inf-sup stability and optimal_convergence, without stabilisation, when_using continuous_approximation_of polynomial order_less_than or_equal to $2$, for both the bulk_variable and_the multiplier provided $\rho_h = O(h^2)$._The approximation order of_the_interface normal, which is $O(h)$ prohibits_higher order convergence