$ is described in polar coordinates by $$\label{eq: 0 nbhd}
0<r<r_0\quad\text{and}\quad 0<\theta<\pi/\beta,$$ as illustrated in Figure \[fig: Omega\]. We denote the vertices of $\Omega$ by $p_0=(0,0)$, $p_1$, $p_2$, …, $p_J=p_0$, and the $j$th side by $$\Gamma_j=(p_j,p_{j+1})
    =\{\,(1-\sigma)p_j+\sigma p_{j+1}: 0<\sigma<1\,\}
    \quad\text{for $0\le j\le J-1$.}$$

Section \[sec: elliptic\] summarizes some key facts about the singular behaviour of the solution to the elliptic problem . In Section \[sec: FEM\], we describe a family of shape-regular triangulations ${\mathcal{T}}_h$ (indexed by the mesh parameter $h$) that depend on a local refinement parameter $\gamma\ge1$. The elements near the origin have sizes of order $h^\gamma$, so the ${\mathcal{T}}_h$ are quasi-uniform if $\gamma=1$ but become more highly refined with increasing $\gamma$. Our error bounds will be stated in terms of the quantity $$\label{eq: epsilon}
\epsilon(h,\gamma)=\begin{cases}
        
h^{\gamma\beta}/\sqrt{\gamma^{-1}-\beta},&1\le\gamma<1/\beta,\\
        h\sqrt{\log(1+h^{-1})},&\gamma=1/\beta,\\
        h/\sqrt{\beta-\gamma^{-1}},&\gamma>1/\beta,
\end{cases}$$ which ranges in size from $O(h^\beta)$ when $\gamma=1$ (the quasiuniform case) down to $O(h)$ when $\gamma>1/\beta$. We briefly review results for the finite element approximation of the elliptic problem, needed for our subsequent analysis: the error in $H^1(\Omega)$ is of order $\epsilon(h,\gamma)$, and the error in $L_2(\Omega)$ is of order $\epsilon(h,\gamma)^2$.

Section \[sec: time-dep\] gathers together some
$ is described in polar coordinates by $ $ \label{eq: 0 nbhd } 
 0 < r < r_0\quad\text{and}\quad 0<\theta<\pi/\beta,$$ as illustrate in Figure   \[fig: Omega\ ]. We announce the vertices of   $ \Omega$ by $ p_0=(0,0)$, $ p_1 $, $ p_2 $, …, $ p_J = p_0 $, and the $ j$th side by $ $ \Gamma_j=(p_j, p_{j+1 }) 
     = \{\,(1-\sigma)p_j+\sigma p_{j+1 }: 0<\sigma<1\,\ } 
     \quad\text{for $ 0\le j\le J-1$.}$$ 

 Section   \[sec: elliptic\ ] sum up some cardinal facts about the singular behaviour of the solution to the egg-shaped trouble  . In Section   \[sec: FEM\ ], we describe a family of shape - regular triangulations   $ { \mathcal{T}}_h$ (index by the mesh parameter   $ h$) that count on a local refinement parameter   $ \gamma\ge1$. The elements near the beginning have sizes of ordering   $ h^\gamma$, so the $ { \mathcal{T}}_h$ are quasi - uniform if $ \gamma=1 $ but become more highly refined with increasing   $ \gamma$. Our erroneousness bounds will be stated in terms of the quantity $ $ \label{eq: epsilon } 
 \epsilon(h,\gamma)=\begin{cases } 
        
 h^{\gamma\beta}/\sqrt{\gamma^{-1}-\beta},&1\le\gamma<1/\beta,\\ 
         h\sqrt{\log(1+h^{-1})},&\gamma=1/\beta,\\ 
         h/\sqrt{\beta-\gamma^{-1}},&\gamma>1/\beta, 
 \end{cases}$$ which range in size from $ O(h^\beta)$ when   $ \gamma=1 $ (the quasiuniform case) down to $ O(h)$ when   $ \gamma>1/\beta$. We briefly review results for the finite element approximation of the egg-shaped problem, needed for our subsequent analysis: the error in $ H^1(\Omega)$ is of order   $ \epsilon(h,\gamma)$, and the error in   $ L_2(\Omega)$ is of club   $ \epsilon(h,\gamma)^2$. 

 Section   \[sec: prison term - dep\ ] gathers together some
$ is described in polar coorainates by $$\label{eq: 0 nbid}
0<r<r_0\quas\text{and}\duad 0<\theta<\pi/\beta,$$ as illustraved un Fiture \[fig: Omega\]. We denotd the vernices of $\Onega$ vy $p_0=(0,0)$, $p_1$, $p_2$, …, $'_N=p_0$, and bke $j$tg sidz uy $$\Gamma_j=(p_j,p_{j+1})
    =\{\,(1-\sigma)p_j+\shgma p_{j+1}: 0<\sigma<1\,\}
    \qbad\text{for $0\le j\le J-1$.}$$

Section \[sec: elliptyc\] summsrlzes some key sactx abohn uhe singular behaviour of the solhtion tm the elliptiv problem . In Section \[sec: FEL\], we describe a family of shape-retulaw triangulatiuns ${\mathcal{U}}_h$ (indexed bi the mesh parameter $h$) that depena on c local refunwmejj parameter $\gemma\ge1$. The elements near tve origon have sizes pf mrdwr $h^\gamma$, so the ${\mathcel{T}}_h$ are quasi-unifori if $\gammd=1$ yut become more highlt eefingd widh ivxrexsihg $\jamja$. Our ercor bounds sill be stared in terms of the qtqntity $$\label{es: epsijog}
\epsilon(h,\gamma)=\begin{cases}
        
h^{\gamma\btta}/\sqdt{\gamma^{-1}-\beta},&1\le\gamma<1/\beta,\\
        h\sqrt{\log(1+h^{-1})},&\gamma=1/\heta,\\
        h/\sqrt{\beta-\gamma^{-1}},&\gamma>1/\beta,
\end{cases}$$ which ranges it sizx wron $J(f^\veha)$ when $\gamma=1$ (the quasiuniform case) down to $O(h)$ shtn $\gsmma>1/\beta$. We brlefly review resulys fpt the finite euement aplroximation of the elliptyc priblem, neeqed gor our subsequent analysis: the error iu $H^1(\Imega)$ is of order $\e'silon(h,\gamma)$, and the rrror in $L_2(\Omega)$ is of orber $\epsjlon(h,\gamma)^2$.

Sfction \[sec: gime-dep\] gathers gogvthes some
$ is described in polar coordinates by nbhd} 0<\theta<\pi/\beta,$$ as in Figure \[fig: of by $p_0=(0,0)$, $p_1$, …, $p_J=p_0$, and $j$th side by $$\Gamma_j=(p_j,p_{j+1}) =\{\,(1-\sigma)p_j+\sigma p_{j+1}: \quad\text{for $0\le j\le J-1$.}$$ Section \[sec: elliptic\] summarizes some key facts about the behaviour of the solution to the elliptic problem . In Section \[sec: FEM\], describe family shape-regular ${\mathcal{T}}_h$ (indexed by the mesh parameter $h$) that depend on a local refinement parameter $\gamma\ge1$. The near the origin have sizes of order $h^\gamma$, the ${\mathcal{T}}_h$ are quasi-uniform $\gamma=1$ but become more highly with $\gamma$. Our bounds be in terms of quantity $$\label{eq: epsilon} \epsilon(h,\gamma)=\begin{cases} h^{\gamma\beta}/\sqrt{\gamma^{-1}-\beta},&1\le\gamma<1/\beta,\\ h\sqrt{\log(1+h^{-1})},&\gamma=1/\beta,\\ h/\sqrt{\beta-\gamma^{-1}},&\gamma>1/\beta, \end{cases}$$ which ranges in size from $O(h^\beta)$ when $\gamma=1$ (the case) down when $\gamma>1/\beta$. briefly results the finite element the elliptic problem, needed for our error in $H^1(\Omega)$ is of order $\epsilon(h,\gamma)$, and error in is of order $\epsilon(h,\gamma)^2$. Section \[sec: gathers together some
$ is described in polar coordinAtes by $$\labeL{eq: 0 nbHd}
0<r<R_0\quAd\Text{And}\qUad 0<\theta<\pi/\beta,$$ AS illUstrated in Figure \[fig: OmeGa\]. We dEnOTe thE VeRticeS of $\OmegA$ By $P_0=(0,0)$, $P_1$, $P_2$, …, $p_J=P_0$, aNd The $J$tH SiDe by $$\GAmmA_j=(p_j,p_{j+1})
    =\{\,(1-\sIgma)p_j+\sigmA p_{j+1}: 0<\SiGma<1\,\}
    \quad\text{fOR $0\lE j\le J-1$.}$$

SectiOn \[sEc: elliptic\] suMmaRizes sOmE keY Facts AboUt the SingulAR behavIour of the SoLUtion tO The elliPTIc ProbLem . In Section \[sec: FEm\], We DEscribe a family Of shapE-rEGuLAR trIanGulations ${\mAtHcal{T}}_H$ (Indexed BY tHE MEsh PArameter $h$) that Depend on a loCAl rEfinemEnT paRAmeter $\Gamma\Ge1$. tHe eLements near The oRigin have Sizes oF Order $h^\gAMma$, so thE ${\mathcAl{T}}_H$ arE quaSI-uNiForM iF $\GamMA=1$ bUt bEComE more higHlY rEfineD witH INCReasIng $\GammA$. Our eRror bounds wilL be StatED in Terms Of the QuanTiTy $$\labEl{eq: epSilon}
\EpSilon(h,\gamma)=\begiN{casEs}
        
h^{\gamma\bEta}/\SqRt{\gAmMa^{-1}-\betA},&1\Le\gammA<1/\beTa,\\
        h\Sqrt{\log(1+H^{-1})},&\gamma=1/\bETa,\\
        h/\SqRT{\BEtA-\gamma^{-1}},&\gamma>1/\beta,
\end{CaSES}$$ wHich rangEs in siZE fRoM $o(h^\beta)$ whEn $\GamMa=1$ (thE QUasiuNifoRM cAse) down tO $O(h)$ wheN $\GaMmA>1/\beta$. We BrIefly rEvIew ResUlts fOR the Finite Element aPproxIMation of the ellIPtic problem, neEDeD FOr OUr suBseQuent analysIs: thE ErroR in $H^1(\oMeGa)$ iS Of ordEr $\epsIlON(h,\GAmma)$, and the error in $L_2(\OMeGa)$ is of Order $\Epsilon(h,\gamma)^2$.

section \[sec: TIME-dep\] gathErs tOGeTHer some
$ is described in polar co ordinatesby $$ \la bel {e q: 0  nbh d}
0<r<r_0\qua d \tex t{and}\quad 0<\theta<\ pi/\b et a ,$$a sillus tratedi nF i gur e\[ fig :O me ga\].  We  denote  the verti ces  o f $\Omega$ b y  $ p_0=(0,0)$ , $ p_1$, $p_2$,  …,  $p_J= p_ 0$, and t he$j$th  sideb y $$\G amma_j=(p _j , p_{j+1 } )
    = \ { \, (1-\ sigma)p_j+\sigmap _{ j +1}: 0<\sigma< 1\,\}  \ q u ad\ tex t{for $0\l ej\leJ -1$.}$$ 
S e c t ion \[sec: ellipt ic\] summar i zes  someke y f a cts ab out t he sin gular behav iour  of the s olutio n  to the ellipti c prob lem  .In S e ct io n \ [s e c:F EM \], wedescribe  a  f amily  ofs h a p e-re gul ar t riang ulations ${\m ath cal{ T }}_ h$ (i ndexe d by  t he me sh par amete r$h$) that depen d on  a localref in eme nt  para m eter $ \ga mma \ge1$.The ele m ent sn e a rthe origin have si ze s of  order $ h^\gam m a$ ,s o the ${ \m ath cal{ T } }_h$areq ua si-unifo rm if$ \g am ma=1$ b ut  becom emor e h ighly refi ned wi th incre asing $\gamma$. Oure rror bounds w i ll b es tate d i n terms oftheq uant ity$ $\ lab e l{eq:  epsi lo n }\ epsilon(h,\gamma)=\ be gin{ca ses}        
h^{\ gamma\beta } / \ sqrt{\ga mma^ { -1 } -\beta},&1\le\ gamma <1/\beta,\ \ 
       h\sq rt{\log( 1+h^{-1}) } , &\gamma= 1/\ bet a,\ \
    h/\sqrt{\b e t a-\g am ma^{-1} },& \gamma> 1/\ bet a,\en d{ cases}$$which ra ng es  i nsiz e fro m  $O(h^\b et a)$  w hen  $\ga m ma=1$(thequas iu ni f orm  case)d ow n to $ O( h) $ wh en$\ gamma >1/\ b eta $. We b riefly re vie w  res ul ts  for th e finite elem en t approxim at ion  of th e elliptic  problem, needed for ou r  subseq uen t ana lysi s: the er ror  in $H ^1( \ Omega) $ is o f ord er  $\ e p silon ( h ,\ gam ma )$, and th e err or in  $ L_2( \Omega) $ is of order $\ep s ilo n(h,\gamma)^2 $.
Sec t i on  \[ s ec :  ti me - dep \ ]  gathers togeth er some
$ is_described in_polar coordinates by $$\label{eq:_0 nbhd}
0<r<r_0\quad\text{and}\quad_0<\theta<\pi/\beta,$$_as illustrated_in_Figure \[fig: Omega\]. We_denote the vertices_of $\Omega$ by $p_0=(0,0)$, $p_1$,_$p_2$, …, $p_J=p_0$,_and_the $j$th side by $$\Gamma_j=(p_j,p_{j+1})
    =\{\,(1-\sigma)p_j+\sigma p_{j+1}: 0<\sigma<1\,\}
   _\quad\text{for_$0\le j\le_J-1$.}$$

Section \[sec:_elliptic\]_summarizes some key facts about_the singular behaviour of the_solution to_the elliptic problem . In Section \[sec: FEM\], we describe_a_family of shape-regular_triangulations ${\mathcal{T}}_h$ (indexed by the mesh parameter $h$) that depend on_a local refinement parameter $\gamma\ge1$. The elements_near the origin_have_sizes_of order $h^\gamma$, so the_${\mathcal{T}}_h$ are quasi-uniform if $\gamma=1$ but_become more highly refined with increasing $\gamma$._Our error bounds will be stated in_terms of the quantity $$\label{eq: epsilon}
\epsilon(h,\gamma)=\begin{cases}
_    _ _
h^{\gamma\beta}/\sqrt{\gamma^{-1}-\beta},&1\le\gamma<1/\beta,\\
    _  _h\sqrt{\log(1+h^{-1})},&\gamma=1/\beta,\\
 _  _   h/\sqrt{\beta-\gamma^{-1}},&\gamma>1/\beta,
\end{cases}$$ which ranges_in size from_$O(h^\beta)$ when $\gamma=1$ (the quasiuniform case) down_to_$O(h)$ when $\gamma>1/\beta$. We_briefly_review_results for_the finite element_approximation_of the_elliptic_problem, needed for our subsequent analysis:_the_error in $H^1(\Omega)$ is of order $\epsilon(h,\gamma)$, and_the error in $L_2(\Omega)$ is_of_order $\epsilon(h,\gamma)^2$.

Section \[sec: time-dep\] gathers together_some