Maximal regularity of evolving FEMs
for parabolic equations on an evolving surface

Abstract

In this paper, we prove that spatially semi-discrete evolving finite element method for parabolic equations on a given evolving hypersurface of arbitrary dimensions preserves the maximal $L^{p}$-regularity at the discrete level. We first establish the results on a stationary surface and then extend them, via a perturbation argument, to the case where the underlying surface is evolving under a prescribed velocity field. The proof combines techniques in evolving finite element method, properties of Green’s functions on (discretised) closed surfaces, and local energy estimates for finite element methods.

Keywords:    parabolic equation, evolving surface, evolving FEM, maximal $L^{p}$-regularity, local energy estimates.

1 Introduction

Overview of maximal regularity results in flat domains. Let $\varOmega\subset{\mathbb{R}}^{d}$ be a flat, polyhedral of dimensions $d\in\{2,3\}$ or smooth domain of arbitrary dimensions, and consider the initial and boundary value problem for a linear parabolic partial differential equation (PDE)

$$\left\{\begin{array}[]{ll}\displaystyle\frac{\partial u(t,x)}{\partial t}-\sum_{i,j=1}^{d}\frac{\partial}{\partial x_{i}}\bigg{(}a_{ij}(x)\frac{\partial u(t,x)}{\partial x_{j}}\bigg{)}=f(t,x)&\mbox{for}\,\,(t,x)\in{\mathbb{R}}_{+}\times\varOmega,\\[10.0pt] u(t,x)=0&\mbox{for}\,\,(t,x)\in{\mathbb{R}}_{+}\times\partial\varOmega,\\[7.0pt] u(0,x)=u_{0}(x)&\mbox{for}\,\,x\in\varOmega,\end{array}\right.$$ (1.4)

where $a_{ij}=a_{ji}$ are real-valued bounded measurable functions satisfying the following uniform ellipticity condition for some constant $\lambda>0$:

$$\lambda^{-1}|\xi|^{2}\leq\mbox{$\sum_{i,j=1}^{d}$}a_{ij}(x)\xi_{i}\xi_{j}\leq\lambda|\xi|^{2},\quad\forall\,\xi=(\xi_{1},\dotsc,\xi_{d})\in{\mathbb{R}}^{d},\,\,\,\forall\,x\in\varOmega.$$ (1.5)

Under condition (1.5), the elliptic partial differential operator $A=\sum_{i,j=1}^{d}\frac{\partial}{\partial x_{i}}\big{(}a_{ij}(x)\frac{\partial}{\partial x_{j}}\big{)}$ generates a bounded analytic semigroup on $L^{q}(\varOmega)$, $1<q<\infty$, and the solution of (1.4) possesses the following maximal $L^{p}$-regularity in $L^{q}(\varOmega)$:

$$\|\partial_{t}u\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\varOmega))}+\|Au\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\varOmega))}\leq C\|f\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\varOmega))},\quad\mbox{if}\,\,\,u_{0}=0,\,\,1<p,q<\infty,$$ (1.6)

which is an important mathematical tool in studying the well-posedness and regularity of solutions of nonlinear parabolic PDEs; see [4, 9, 36, 27, 39, 41].

Analogously, denoting by $A_{h}$ the finite element approximation of the elliptic operator $A$ on a finite element subspace $S_{h}\subset H^{1}_{0}(\varOmega)$, i.e.,

$$(A_{h}\phi_{h},\varphi_{h})=-(a_{ij}\nabla\phi_{h},\nabla\varphi_{h}),\quad\forall\,\phi_{h},\varphi_{h}\in S_{h},$$ (1.7)

it is known that the semi-discrete finite element solutions given by

$$\left\{\begin{array}[]{ll}(\partial_{t}u_{h},v_{h})+\sum_{i,j=1}^{d}(a_{ij}\partial_{j}u_{h},\partial_{i}v_{h})=(f,v_{h}),&\forall\,v_{h}\in S_{h},\,\,\forall\,t\in(0,T),\\[5.0pt] u_{h}(0)=u_{h,0},\end{array}\right.$$ (1.10)

satisfies the following spatially discrete $h$-uniform maximal $L^{p}$-regularity [31, 35] (with a constant $C>0$ independent of the mesh size $h$):

$$\|\partial_{t}u_{h}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\varOmega))}+\|A_{h}u_{h}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\varOmega))}\leq C\|f\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\varOmega))}, \mbox{if}\,\,\,u_{h,0}=0,\,\,\,1<p,q<\infty,$$ (1.11)

which has applications in numerical analysis for semilinear parabolic equations with strong nonlinearities [18], and quasi-linear parabolic equations with nonsmooth coefficients [34]. The spatially discrete maximal $L^{p}$-regularity results were firstly proved in smooth domains with Neumann boundary condition [17, 31], and then extended to polyhedral domains [35, 32] with the Dirichlet boundary condition. The discrete maximal $L^{p}$-regularity is also closely related (in the techniques of proof) to the maximum-norm stability of finite element solutions of parabolic equations [29, 40, 42, 43, 47]:

$$\|u-u_{h}\|_{L^{\infty}(0,T;L^{\infty}(\Omega))}\leq C\ln(2+1/h)\inf_{\chi_{h}\in L^{\infty}(0,T;S_{h})}\|u-\chi_{h}\|_{L^{\infty}(0,T;L^{\infty}(\Omega))}.$$ (1.12)

The extension of maximal $L^{p}$-regularity to the time-discrete setting was established for different time discretization methods, including the backward Euler method [5], discontinuous Galerkin method [30], $\theta$-schemes [21], and A-stable multistep and Runge-Kutta methods [24]. All these methods are A-stable. The maximal $L^{p}$-regularity of A$(\alpha)$-stable backward differentiation formulae (BDF) was established in [33]. The discrete maximal $L^{p}$-regularity helps us to control the nonlinearlity as well (see [1, 26, 23]), and besides it enables us to obtain optimal-order $L^{p}$-norm error estimates without using the inverse inequality (cf. [24, 32, 23]).

Overview of maximal regularity results on surfaces. The maximal $L^{p}$-regularity of parabolic equations on an evolving surface $\Gamma\subset{\mathbb{R}}^{m+1}$, $m\in\mathbb{N}$, as well as the maximal $L^{p}$-regularity of time discretizations on an evolving surface and its application to the convergence analysis of BDF methods for nonlinear PDEs on an evolving surface, was discussed in [23].

Semi-discrete maximal regularity results on evolving surfaces. However, since the maximal $L^{p}$-regularity of spatial discretizations for parabolic equations on an evolving surface is still missing, only semi-discretization in time were considered in [23]. The aim of this article is to fill in this gap, by establishing spatially discrete maximal $L^{p}$-regularity of isoparametric finite element methods (FEMs) for parabolic equations on an evolving surface which is approximated by quasi-uniform curved triangles.

In order to prove the discrete maximal $L^{p}$-regularity for the spatially semi-discrete problems, we combine the techniques developed for evolving surface FEMs and local energy estimates. Firstly, we shall prove the discrete maximal $L^{p}$-regularity for spatially semi-discrete FEM in (2.12) on a stationary surface (Theorem 2.2). Then we use a temporal perturbation argument to extend this result to evolving surfaces (Theorem 2.3). The discrete maximal $L^{p}$-regularity results for (2.14) can be obtained analogously by a perturbation argument for the lower-order advection term.

Since we are considering a spatially semi-discrete setting, the underlying smooth surface $\Gamma(t)$ should also be replaced by the finite element surface $\Gamma_{h}(t)$. The discrepancy between $\Gamma(t)$ and $\Gamma_{h}(t)$ is the main obstacle in the proof, leading to the following technical difficulties to be addressed:

• The discrete delta functions on $\Gamma(t)$ and $\Gamma_{h}(t)$ are not simply related by the lift via the distance projection. Indeed, they are correlated via a nonlinear way which stems from the nonlinear relation of $L^{2}$ projections $P_{h}(\Gamma(t))$ and $P_{h}(\Gamma_{h}(t))$. Therefore, it is necessary for us to obtain the high-order consistency between the discrete delta functions on different surfaces (Lemma 3.3) in order to ensure the consistency of the corresponding Green’s functions, which are indispensable in the used local energy estimate (Lemma 3.5) and dyadic decomposition argument (Lemma 4.1).

• The discrepancy of $\Gamma(t)$ and $\Gamma_{h}(t)$ will introduce a bunch of additional geometric perturbation terms in the local energy estimate (Lemma 3.5) and in the dyadic decomposition argument (Lemma 4.1). We need to treat them carefully to make sure that the leading order stability is still available.

• In the temporal perturbation argument, we need to develop the norm equivalence of the discrete Laplacian (Lemma 5.1 and Remark 5.1). The super-approximation property (cf. (P3) in Section 3.2) also plays an important role in the derivation of this equivalence. Besides, as a nature of parametric finite elements, the matrix-valued coefficient $B_{h}(t,x)$ of the following change of variable

  $$\int_{\Gamma_{h}(t)}\nabla_{\Gamma_{h}(t)}u^{-\ell}\nabla_{\Gamma_{h}(t)}v^{-\ell}=\int_{\Gamma(t)}B_{h}(t,\cdot)\nabla_{\Gamma(t)}u\nabla_{\Gamma(t)}v$$

  always has jumps at the edges and thus is discontinuous. To this end, it is desirable as well to construct a globally continuous substitute $\tilde{B}(t,x)$ according to the definition of the discrete Laplacian.

• The norm equivalence of the discrete Laplacian will bring in a lower-order term $\|\nabla u\|_{L^{p}}$. To control this term by the maximal $L^{p}$-regularity, we need to use the discrete interpolation inequality (Lemma 5.2) whose proof greatly relies on the $W^{1,p}$-stability of the Ritz projection. The latter is a consequence of the Green’s function estimate on closed manifold ([13, Theorem 3.2]).

The article is organized as follows: In Section 2, we introduce the basic notations for evolving surface FEMs, the semi-discrete evolving surface FEMs for parabolic equations on an evolving surface, and the main theoretical results about discrete maximal $L^{p}$-regularity of evolving surface FEMs. In Section 3, we develop the preliminary results of geometric perturbation estimates, Green’s function estimates and local energy estimates on the stationary surface. We will prove the maximal regularity on stationary surface (Theorem 2.2) and on evolving surface via a temporal perturbation argument (Theorem 2.3) in Section4 and Section 5, respectively. In Appendix A we present the detailed proof of local energy estimates (Lemma 3.5), which is a technical lemma used in the proof of the main theorems.

2 Evolving surface FEMs and main results

We now introduce the problem to be considered in this paper and briefly recall the standard parametric evolving surface FEM, see [14, 13, 15, 22].

2.1 Parabolic equations on an evolving surface

We assume that the evolution of a closed hypersurface $\Gamma(t)\subset{\mathbb{R}}^{m+1}$ of arbitrary dimensions is given by a diffeomorphic flow map $X(t,\cdot)\colon\Gamma^{0}\to\Gamma(t)$, where $m\in\mathbb{N}$, and $\Gamma^{0}$ is a smooth $m$-dimensional initial hypersurface, with $X(0,\cdot)$ being the identity map on $\Gamma^{0}$. We assume that $X(t,y)$ is smooth with respect to $(t,y)\in[0,T]\times\Gamma^{0}$ and the inverse function $X^{-1}(t,x)$ is smooth with respect to $x\in\Gamma(t)$ uniformly for $t\in[0,T]$.

The material velocity (which is simply called velocity below) and material derivative on the surface are respectively given, for $x=X(y,t)\in\Gamma(t)$ with $y\in\Gamma^{0}$, by

$$v(t,x)= \ \partial_{t}X(t,y),$$ (2.1)

and

$$\partial_{t}^{\bullet}u(t,x)= \ \frac{{\mathrm{d}}}{{\mathrm{d}}t}u(t,X(t,y)).$$

Let $\nu$ be the unit outward normal vector to the surface $\Gamma(t)$. We denote by $\nabla_{\Gamma(t)}u$ the tangential gradient of the function $u$, and denote by $\Delta_{\Gamma(t)}u=\nabla_{\Gamma(t)}\cdot\nabla_{\Gamma(t)}u$ the Laplace–Beltrami operator acting on $u$. For more details on all these basic concepts we refer to [14, 12, 15], and the references therein. An unified abstract theory for evolving mesh FEMs for partial differential equations on evolving domains can be found in [16].

We consider the following two types of linear parabolic equations on the closed surface $\Gamma(t)$.

(i) A parabolic model problem on a given evolving surface:

$$\left\{\begin{aligned} &\partial_{t}^{\bullet}u-\Delta_{\Gamma(t)}u=f&\qquad&\mbox{on}\,\,\,\Gamma(t),\,\,\,\forall\,t\in(0,T),\\ &u(0,\cdot)=u^{0}&&\mbox{on}\,\,\,\Gamma^{0},\end{aligned}\right.$$ (2.2)

where $\partial_{t}^{\bullet}$ is the material derivative associated to the flow velocity $v$. The above model problem arises in the analysis and numerical analysis of mean curvature flow, see [20] and [25], respectively, wherein the following evolution equations of normal vector $\nu$ and mean curvature $H$ play an important role:

$$\partial_{t}^{\bullet}\nu-\Delta_{\Gamma(t)}\nu=|\nabla_{\Gamma(t)}\nu|^{2}\nu\quad\mbox{and}\quad\partial_{t}^{\bullet}H-\Delta_{\Gamma(t)}H=|\nabla_{\Gamma(t)}\nu|^{2}H.$$

(ii) The heat equation on a given evolving surface (see, e.g., [14, equation (1.1)]):

$$\left\{\begin{aligned} &\partial_{t}^{\bullet}u+u\nabla_{\Gamma(t)}\cdot v-\Delta_{\Gamma(t)}u=f&\qquad&\mbox{on}\,\,\,\Gamma(t),\,\,\,\forall\,t\in(0,T),\\ &u(0,\cdot)=u^{0}&&\mbox{on}\,\,\,\Gamma^{0}.\end{aligned}\right.$$ (2.3)

The regularity of solutions to PDEs on an evolving surface is often characterized by the following Sobolev spaces on a space-time manifold $\mathcal{G}_{T}=\bigcup_{t\in(0,T)}(\{t\}\times\Gamma(t))$:

$$L_{t}^{p}(0,T;W^{k,q}(\Gamma(t))) =\bigg{\{}w\colon\mathcal{G}_{T}\rightarrow{\mathbb{R}}:w(t,\cdot)\in W^{k,q}(\Gamma(t))\mbox{ a.e. }t\in(0,T),$$

$$\hskip 121.0ptt\rightarrow\|w(t,\cdot)\|_{W^{k,q}(\Gamma(t))}\in L^{p}(0,T)\bigg{\}},$$ (2.4)

$$W_{t}^{1,p}(0,T;W^{k,q}(\Gamma(t))) =\bigg{\{}w\in L_{t}^{p}(0,T;W^{k,q}(\Gamma(t))):\partial_{t}^{\bullet}w\in L_{t}^{p}(0,T;W^{k,q}(\Gamma(t)))\bigg{\}}.$$ (2.5)

The conventional notational convention $H^{k}(\Gamma(t))=W^{2,k}(\Gamma(t))$ will also be used. For more details on these spaces we refer to [3, 2] (employing a different notation).

The weak formulation of (2.2) reads as follows. Find $u\in H^{1}(\mathcal{G}_{T})\cap L_{t}^{2}(0,T;H^{1}(\Gamma(t)))$ satisfying relation

$$\int_{\Gamma(t)}\partial_{t}^{\bullet}u\,\varphi+\int_{\Gamma(t)}\nabla_{\Gamma(t)}u\cdot\nabla_{\Gamma(t)}\varphi=\int_{\Gamma(t)}f\varphi$$ (2.6)

for all $\varphi\in H^{1}(\mathcal{G}_{T})$ and almost all $t\in(0,T)$.

The weak formulation of (2.3) reads as follows (see, e.g., [14, equation (1.2)]): Find $u\in H^{1}(\mathcal{G}_{T})\cap L_{t}^{2}(0,T;H^{1}(\Gamma(t)))$ satisfying relation

$$\frac{{\mathrm{d}}}{{\mathrm{d}}t}\int_{\Gamma(t)}u\,\varphi+\int_{\Gamma(t)}\nabla_{\Gamma(t)}u\cdot\nabla_{\Gamma(t)}\varphi=\int_{\Gamma(t)}f\varphi$$ (2.7)

for almost all $t\in(0,T)$ and all $\varphi\in H^{1}(\mathcal{G}_{T})$ with $\partial_{t}^{\bullet}\varphi=0$.

2.2 Evolving surface finite elements

Let $\Gamma_{h}(t)$ be the closed and continuous piecewise polynomial surface (of degree $k$) which approximates the smooth surface $\Gamma(t)$ evolving under the prescribed velocity $v$. Additionally, we require that the nodes of $\Gamma_{h}(t)$ stay on $\Gamma(t)$ and move with the same velocity as $\Gamma(t)$. In other words, the evolution of $\Gamma_{h}(t)$ is uniquely determined with velocity $v_{h}=I_{h}v$.

Each polynomial element $K$ of $\Gamma_{h}(t)$ is the image of an element $K^{0}\subset\Gamma_{h}(0)$ under the discrete flow map. We denote by $K_{\rm f}^{0}$ the unique flat element which has the same endpoints as $K^{0}$, and denote by $F_{K}\colon K_{\rm f}^{0}\rightarrow K$ the parametrization of $K$, i.e., $F_{K}$ is the unique polynomial of degree $k$ that maps $K_{\rm f}^{0}$ onto $K$. For a more detailed description of high-order surface triangulations we refer to [13, 22]. The finite element space of degree $k$ on the discrete surface $\Gamma_{h}(t)$ is defined as

$$S_{h}(\Gamma_{h}(t)):= \ \big{\{}\varphi_{h}\in C(\Gamma_{h}(t)):\varphi_{h}\circ F_{K}\in\mathbb{P}^{k}(K_{\rm f}^{0})\,\,\mbox{for every element}\,\,K\subset\Gamma_{h}(t)\big{\}}$$

$$= \ \textnormal{span}\big{\{}\phi_{1},\dotsc,\phi_{N}\big{\}},$$

where $\mathbb{P}^{k}(K_{\rm f}^{0})$ denotes the space of polynomials of degree $k\geq 1$ on the flat element $K_{\rm f}^{0}$, and where $\phi_{j}$ denote the standard nodal basis functions of $S_{h}(\Gamma_{h}(t))$. By construction $S_{h}(\Gamma_{h}(t))$ is an isoparametric finite element space.

We denote by $\partial_{h,t}^{\bullet}$ the discrete material derivative associate to the discrete flow $v_{h}$. The basis functions $\phi_{j}$ satisfy the following transport property (see [14, Proposition 5.4]):

$$\partial_{h,t}^{\bullet}\phi_{j}=0\qquad\text{for}\quad j=1,\dotsc,N.$$ (2.8)

Let $\delta>0$ be a fixed sufficiently small constant such that every point $x$ in the $\delta$-neighbourhood of surface $\Gamma(t)$, denoted by $D_{\delta}(\Gamma(t))=\{x\in{\mathbb{R}}^{m+1}:{\rm dist}(x,\Gamma(t))\leq\delta\}$, has a unique distance projection onto $\Gamma(t)$, denoted by $q(t,x)$, satisfying the following relation:

$$q(t,x)-x=|x-q(t,x)|\nu(t,q(t,x)),$$

where $\nu(t,\cdot)$ is the outward unit normal vector to $\Gamma(t)$. It is known that such a constant $\delta$ exists and only depends on the curvature of $\Gamma(t)$ (thus $\delta$ is independent of $t\in[0,T]$, but possibly dependent on $T$). We will use the notation $q^{-1}(t,\cdot)$ to denote the inverse of the bijective map $q(t,\cdot)|_{\Gamma_{h}(t)}\colon\Gamma_{h}(t)\rightarrow\Gamma(t)$.

We assume that each element $K^{0}\subset\Gamma_{h}(0)$ interpolates the smooth initial surface $\Gamma(0)$ and that the parametrization $F_{K^{0}}\colon K_{\rm f}^{0}\rightarrow K^{0}$ is a polynomial of degree at most $k$ with the following property:

$$\max_{K^{0}\subset\Gamma_{h}(0)}\Big{(}\|F_{K^{0}}\|_{W^{k,\infty}(K_{\rm f}^{0})}+\|\nabla_{K^{0}}F_{K^{0}}^{-1}\|_{L^{\infty}(K^{0})}\Big{)}\leq\kappa_{0},$$ (2.9)

where $\kappa_{0}$ is some constant that is independent of $h$. This property holds for standard parametric finite elements which interpolate the smooth surface $\Gamma(0)$ and guarantees the following optimal-order approximation to $\Gamma(0)$ by $\Gamma_{h}(0)$ (see [13, 15]):

$$\max_{K^{0}\subset\Gamma_{h}(0)}\|q(0)\circ F_{K^{0}}-F_{K^{0}}\|_{L^{\infty}(K_{\rm f}^{0})}\leq Ch^{k+1}.$$ (2.10)

The projection $q(0,x)$ is well defined for points $x$ in a neighbourhood of $\Gamma(0)$ and therefore well defined on $\Gamma_{h}(0)$ for sufficiently small mesh size $h$. Since the discrete flow $v_{h}=I_{h}v$ is piecewise smooth, at any time $t\in[0,T]$ the optimal-order approximation to $\Gamma(t)$ by $\Gamma_{h}(t)$ is also available

$$\max_{K\subset\Gamma_{h}(t)}\|q(t)\circ F_{K}-F_{K}\|_{L^{\infty}(K_{\rm f}^{0})}\leq Ch^{k+1}.$$ (2.11)

2.3 Lift

The lift of a finite element function $w_{h}\in S_{h}(\Gamma_{h}(t))$ onto the smooth surface $\Gamma(t)$ is defined as

$$w_{h}^{\ell}=w_{h}\circ(q(t)|_{\Gamma_{h}(t)})^{-1},$$

see [13, Section 2.4] and [25, Section 3.4]. The inverse lift of a function $f(t,\cdot)\colon\Gamma(t)\to{\mathbb{R}}$ is denoted by $f^{-\ell}(t,\cdot)\colon\Gamma_{h}(t)\to{\mathbb{R}}$, which is the function satisfying that $(f^{-\ell})^{\ell}=f$. By lifting the triangulation from $\Gamma_{h}(t)$ to $\Gamma(t)$, we obtain a triangulation and finite element space on the smooth surface $\Gamma(t)$, i.e.,

$$S_{h}(\Gamma(t))=\big{\{}w_{h}\in C(\Gamma(t)):w_{h}\circ q(t)\circ F_{K}\in\mathbb{P}^{k}(K_{\rm f}^{0})\,\,\mbox{for every element}\,\,K\subset\Gamma_{h}(t)\big{\}}.$$

For the ease of notations, if the specified domain is $\Gamma(t)$ we will also identify $v_{h}$ as its lift on $\Gamma(t)$. Namely, $v_{h}^{\ell}\in S_{h}(\Gamma(t))$ and $w_{h}\in S_{h}(\Gamma_{h}(t))$ simply mean that $v_{h}^{-\ell}\in S_{h}(\Gamma_{h}(t))$ and $w_{h}^{\ell}\in S_{h}(\Gamma(t))$. We will use the explicit superscript ^ℓ to denote the lift where we want to emphasize its underlying domain.

2.4 The semi-discrete problems

We will now formulate the evolving surface finite element semi-discretisations of the two weak formulations in (2.6) and (2.7).

The semi-discretisation of (2.6) reads: Find $u_{h}\in C^{1}_{t}([0,T];S_{h}(\Gamma_{h}(t)))$ such that

$$\left\{\begin{aligned} &(\partial_{t}u_{h}(t),\varphi_{h})_{\Gamma_{h}(t)}+(\nabla_{\Gamma_{h}(t)}u_{h}(t),\nabla_{\Gamma_{h}(t)}\varphi_{h})_{\Gamma_{h}(t)}=(f^{-\ell}(t),\varphi_{h})_{\Gamma_{h}(t)}&&\forall\,t\in(0,T],\\ &u_{h}(0)=u_{h,0}&&\mbox{on}\quad\Gamma_{h}^{0},\end{aligned}\right.$$ (2.12)

for all $\varphi_{h}\in S_{h}(\Gamma_{h}(t))$, where $u_{h}^{0}=I_{h}u(0)$, and $\partial_{t}u_{h}=\partial_{h,t}^{\bullet}u_{h}$ should be interpreted in the nodal sense. Equivalently, (2.12) can be written into the following strong form

$$\left\{\begin{aligned} &\partial_{t}u_{h}(t)-\Delta_{\Gamma_{h}(t),h}u_{h}(t)=f_{h}(t)&&\forall\,t\in(0,T],\\ &u_{h}(0)=u_{h,0}&&\mbox{on}\quad\Gamma_{h}^{0},\end{aligned}\right.$$ (2.13)

where $f_{h}(t,x)=P_{h}(\Gamma_{h}(t))f^{-\ell}(t,x)$ and $P_{h}(\Gamma_{h}(t))\colon L^{2}(\Gamma_{h}(t))\rightarrow S_{h}(\Gamma_{h}(t))$ is the $L^{2}$ projection on surface $\Gamma_{h}(t)$.

Similarly, the semi-discretisation of (2.7) reads: Find $u_{h}\in C^{1}_{t}([0,T];S_{h}(\Gamma_{h}(t)))$ such that

$$\left\{\begin{aligned} &\frac{{\mathrm{d}}}{{\mathrm{d}}t}(u_{h}(t),\varphi_{h})_{\Gamma_{h}(t)}+(\nabla_{\Gamma_{h}(t)}u_{h}(t),\nabla_{\Gamma_{h}(t)}\varphi_{h})_{\Gamma_{h}(t)}=(f^{-\ell}(t),\varphi_{h})_{\Gamma_{h}(t)}&&\forall\,t\in(0,T],\\ &u_{h}(0)=u_{h,0}&&\mbox{on}\quad\Gamma_{h}^{0},\end{aligned}\right.$$ (2.14)

for all $\varphi_{h}\in C^{1}_{t}([0,T];S_{h}(\Gamma_{h}(t)))$ such that $\partial_{h,t}^{\bullet}\varphi_{h}=0$.

2.5 Main results: Spatially discrete maximal regularity

Following [23], we recall that $u\in H^{1}(\mathcal{G}_{T})$ is a solution of (2.2) if and only if the function $U(t,y):=u(t,X(t,y))$, where $X\colon\Gamma(0)\rightarrow\Gamma(t)$ is the flow map, defines a solution $U\in H^{1}(\Gamma(0)\times(0,T))$ of the following weak formulation:

$$\int_{\Gamma(0)}a(t,y)\partial_{t}U(t,y)\psi(y)+\int_{\Gamma(0)}B(t,y)\nabla_{\Gamma(0)}U(t,y)\cdot\nabla_{\Gamma(0)}\psi(y)=\int_{\Gamma(0)}a(t,y)F(t,y)\psi(y),$$ (2.15)

for all $\psi\in H^{1}(\Gamma(0))$ and almost all $t\in[0,T]$, with $F(t,y):=f(X(t,y),t)$ and $U(0,\cdot)=u_{0}$.

Since the Riemannian metric on the evolving surface is positive definite (see [23, Appendix]), it follows that the functions $a(t,y)$ and $B(t,y)$ satisfy the following estimates:

$$C^{-1}\leq a(t,y)\leq C, \quad\forall y\in\Gamma(0),\forall t\in[0,T],$$ (2.16)

$$C^{-1}|\xi_{y}|^{2}\leq B(t,y)\xi_{y}\cdot\xi_{y}\leq C|\xi_{y}|^{2}, \quad\forall y\in\Gamma(0),\forall\xi_{y}\in T\Gamma(0)_{y},\forall t\in[0,T],$$ (2.17)

where $C$ is some positive constant (depending only on the given flow map).

The continuous maximal regularity result for the solution of (2.15) (cf. [23, Theorem 3.1, 3.2], [49, p. 211] and [41]) is summarized below.

Theorem 2.1 (Maximal regularity of continuous PDE, [23, Theorem 3.1])

If $u_{0}=0$, then the solution $U$ of the pulled-back PDE (2.15) obeys the following estimate for $p,q\in(1,\infty)$:

$$\|\partial_{t}U\|_{L^{p}(0,T;L^{q}(\Gamma^{0}))}+\|U\|_{L^{p}(0,T;W^{2,q}(\Gamma^{0}))}\leq C\|F\|_{L^{p}(0,T;L^{q}(\Gamma^{0}))}$$ (2.18)

where the constant $C>0$ only depends on $\mathcal{G}_{T}$.

The main results of this paper are the following two theorems, which concern discrete maximal regularity on discrete stationary surfaces and a discrete evolving surfaces, respectively.

Theorem 2.2 (Maximal regularity of semi-discrete surface FEMs on a stationary surface)

On a discrete stationary surface $\Gamma_{h}(s)$ (for any fixed $s\in[0,T]$), the finite element solution $u_{h}$ of weak formulation

$$\left\{\begin{aligned} &(\partial_{t}u_{h}(t),\varphi_{h})_{\Gamma_{h}(s)}+(\nabla_{\Gamma_{h}(s)}u_{h}(t),\nabla_{\Gamma_{h}(s)}\varphi_{h})_{\Gamma_{h}(s)}=(f_{h},\varphi)_{\Gamma_{h}(s)}&&\forall\,\varphi_{h}\in S_{h}(\Gamma_{h}(s)),\\ &u_{h}(0)=0&&\mbox{on}\,\,\,\Gamma_{h}(s),\end{aligned}\right.$$ (2.19)

satisfies the following estimate for any given $p,q\in(1,\infty)$, and $h\leq h_{0}$:

$$\|\partial_{t}u_{h}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\Gamma_{h}(s)))}+\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\Gamma_{h}(s)))}\leq C\|f_{h}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\Gamma_{h}(s)))}$$ (2.20)

where the constant $C>0$ is independent of $h$ and $s$, but depends on the $\mathcal{G}_{T}$ and $T$.

Theorem 2.3 (Maximal regularity of semi-discrete surface FEM on an evolving surface)

In the case $u_{h,0}=0$, the solution $u_{h}$ of (2.12) obeys the following estimate for any given $p,q\in(1,\infty)$, and $h\leq h_{0}$:

$$\|\partial_{t}u_{h}\|_{L^{p}_{t}(0,T;L^{q}(\Gamma_{h}(t)))}+\|\Delta_{\Gamma_{h}({t}),h}u_{h}\|_{L^{p}_{t}(0,T;L^{q}(\Gamma_{h}(t)))}\leq C\|f_{h}\|_{L^{p}(0,T;L^{q}(\Gamma_{h}(t)))}$$ (2.21)

where $f_{h}(t,x)=P_{h}(\Gamma_{h}(t))f^{-\ell}(t,x)$, and the constant $C>0$ is independent of $h$, but depends on $T$.

The proofs of Theorems 2.2–2.3 are presented in the following sections.

3 Preliminary results on stationary surfaces

In this section, we develop preliminary results of geometric perturbation estimates, Green’s function estimates and local energy estimates on surfaces $\Gamma(s)$ and $\Gamma_{h}(s)$ for some fixed $s\in[0,T]$. For the simplicity of notations, we omit the dependence on $s$ by using $\Gamma$ and $\Gamma_{h}$ instead to denote the surfaces. The constants in this section will not depend on $s$, but may depend on $T$.

3.1 Function spaces

We use the conventional notations of Sobolev spaces $W^{s,q}(\Gamma)$, $s\geq 0$ and $1\leq q\leq\infty$, with abbreviations $L^{q}:=W^{0,q}(\Gamma)$, $W^{s,q}:=W^{s,q}(\Gamma)$ and $H^{s}:=W^{s,2}(\Gamma)$. We denote by $H^{-s}(\Gamma)$ the dual space of $H^{s}_{0}(\Gamma)$. The latter is defined as the closure of $C^{\infty}_{0}(\Gamma)$ in $H^{s}(\Gamma)$.

For any given function $f\colon(0,T)\rightarrow W^{s,q}$ we define the following Bochner norm (for space-time functions):

$$\|f\|_{L^{p}(0,T;W^{s,q})}=\big{\|}\|f(\cdot)\|_{W^{s,q}}\big{\|}_{L^{p}(0,T)},\quad\forall\,\,1\leq p,q\leq\infty,\,\,s\in{\mathbb{R}}.$$ (3.1)

For any subdomain $D\subset\Gamma$, we define

$$\|f\|_{W^{s,q}(D)}:=\inf_{\tilde{f}|_{D}=f}\|\tilde{f}\|_{W^{s,q}(\Gamma)},\quad\forall\,\,1\leq q\leq\infty,\,\,s\in{\mathbb{R}},$$ (3.2)

where the infimum extends over all possible $\tilde{f}$ defined on $\Gamma$ such that $\tilde{f}=f$ in $D$. Similarly, for any subdomain $Q\subset{\mathcal{Q}}=(0,1)\times\Gamma$, we define

$$\|f\|_{L^{p}W^{s,q}(Q)}:=\inf_{\tilde{f}|_{Q}=f}\|\tilde{f}\|_{L^{p}(0,T;W^{s,q})},\quad\forall\,\,1\leq p,q\leq\infty,\,\,s\in{\mathbb{R}},$$ (3.3)

where the infimum extends over all possible $\tilde{f}$ defined on ${\mathcal{Q}}$ such that $\tilde{f}=f$ in $Q$. For more details on these spaces we refer to [3, 2].

In addition, the following notations of inner products will be used:

$$(\phi,\varphi):=\int_{\Gamma}\phi(x)\varphi(x){\mathrm{d}}x,\qquad[u,v]:=\int_{0}^{T}\int_{\Gamma}u(t,x)v(t,x){\mathrm{d}}x\,{\mathrm{d}}t.$$ (3.4)

For any function $w$ defined on ${\mathcal{Q}}$, We denote $w(t)=w(t,\cdot)$. The notation $1_{0<t<T}$ stands for the characteristic function of time interval $(0,T)$, i.e., $1_{0<t<T}(t)=1$ if $t\in(0,T)$ and $1_{0<t<T}(t)=0$ if $t\notin(0,T)$.

3.2 Properties of the finite element space

For any subdomain $D\subset\Gamma$, we denote by $S_{h}(D)$ the space of functions in $S_{h}(\Gamma)$ restricted to the domain $D$, and denote by $S_{h}^{0}(D)$ the subspace of $S_{h}(D)$ consisting of functions which equal zero outside $D$. For $d>0$, we denote by $B_{d}(D)=\{x\in\Gamma:{\rm dist}(x,D)\leq d\}$ a neighborhood of $D$ in $\Gamma$. On a quasi-uniform triangulation of surface $\Gamma$, there exist positive constants $K$ and $\kappa$ such that the triangulation and the corresponding finite element space $S_{h}(\Gamma)$ possess the following properties ($K$ and $\kappa$ are independent of the subset $D$ and $h$).

1. Quasi-uniformity:

   For all triangles $\tau_{l}^{h}$ in the partition, the diameter $h_{l}$ of $\tau_{l}^{h}$ and the radius $\rho_{l}$ of its inscribed ball satisfy

   $$K^{-1}h\leq\rho_{l}\leq h_{l}\leq Kh.$$

2. Inverse inequality:

   If $D$ is a union of elements in the partition, then, for $0\leq k\leq l\leq 1$ and $1\leq q\leq p\leq\infty$,

   $$\|\chi_{h}\|_{W^{l,p}(D)}\leq Kh^{-(l-k)-(m/q-m/p)}\|\chi_{h}\|_{W^{k,q}(D)},\quad\forall\,\,\chi_{h}\in S_{h}.$$

3. Local approximation and super-approximation:

   There exists an operator $I_{h}\colon H^{1}(\Gamma)\rightarrow S_{h}$ with the following properties:

   1. For $v\in H^{2}(\Gamma)$ the following estimate holds, for $\Gamma\subset{\mathbb{R}}^{m+1}$ with $m=1,2,3$, and $h\leq h_{0}$:

      $$\|v-I_{h}v\|_{L^{2}}+h\|\nabla(v-I_{h}v)\|_{L^{2}}\leq Kh^{2}\|v\|_{H^{2}}.$$

   2. If $d\geq 2h$ then the value of $I_{h}v$ in $D$ depends only on the value of $v$ in $B_{d}(D)$. If $d\geq 2h$ and supp$(v)\subset\overline{D}$, then $I_{h}v\in S_{h}^{0}(B_{d}(D))$.

   3. If $d\geq 2h$, $\chi=0$ outside $D$ and $|\partial^{\beta}\chi|\leq Cd^{-|\beta|}$ for all multi-index $\beta$, then

      $$\psi_{h}\in S_{h}(B_{d}(D))\implies I_{h}(\chi\psi_{h})\in S_{h}^{0}(B_{d}(D)),$$

      $$\|\chi\psi_{h}-I_{h}(\chi\psi_{h})\|_{L^{2}}+h\|\chi\psi_{h}-I_{h}(\chi\psi_{h})\|_{H^{1}}\leq Khd^{-1}\|\psi_{h}\|_{L^{2}(B_{d}(D))}.$$

   4. If $d\geq 2h$ and $\chi\equiv 1$ on $B_{d}(D)$, then $I_{h}(\chi\psi_{h})=\psi_{h}$ on $D$.

Properties (P1)–(P3) hold for any quasi-uniform triangulation with the standard finite element spaces consisting of globally continuous piecewise polynomials of degree $r\geq 1$ (cf. [45, Appendix]).

Remark 3.1

Properties (P1)–(P3) hold for parametric surface finite element spaces $S_{h}(\Gamma(t))$ and $S_{h}(\Gamma_{h}(t))$ for all $t\in[0,T]$, see [13, 15].

Given $\Gamma$ and $\Gamma_{h}$, we define the scalar-valued function $a_{h}(x)$ and the ${\mathbb{R}}^{(m+1)\times(m+1)}$-valued function $B_{h}(x)$ with $x\in\Gamma(s)$ to be the piecewise smooth prefactors in the following change of variables:

$$\int_{\Gamma_{h}}uv =\int_{\Gamma}a_{h}(x)u^{\ell}v^{\ell},$$ (3.5)

$$\int_{\Gamma_{h}}\nabla_{\Gamma_{h}}u\cdot\nabla_{\Gamma_{h}}v =\int_{\Gamma}B_{h}(x)\nabla_{\Gamma}u^{\ell}\cdot\nabla_{\Gamma}v^{\ell}.$$ (3.6)

Then the following estimates hold (see [13]):

Lemma 3.1

Let $\Gamma$ and $\Gamma_{h}$ be as above. Then, for sufficiently small $h\leq h_{0}$, the prefactors satisfy the geometric estimates:

$$\|a_{h}-1\|_{L^{\infty}(\Gamma)}+\|B_{h}-I\|_{L^{\infty}(\Gamma)}\leq Ch^{k+1}.$$

Besides, the basic parametric finite element theory says that the $L^{p}$ and $W^{1,p}$ norms on $\Gamma$ and $\Gamma_{h}$ are equivalent (cf. [13, p. 811]), i.e. $\|\varphi_{h}\|_{L^{p}(\Gamma)}\sim\|\varphi_{h}\|_{L^{p}(\Gamma_{h})}$ and $\|\nabla_{\Gamma}\varphi_{h}\|_{L^{p}(\Gamma)}\sim\|\nabla_{\Gamma_{h}}\varphi_{h}\|_{L^{p}(\Gamma_{h})}$ for all $\varphi_{h}\in S_{h}$ and $p\in[1,\infty]$. Here the notation $a\sim b$ means there exist a positive constant $C$ such that $C^{-1}a\leq b\leq Ca$.

3.3 Green’s functions

We consider the heat equation on a fixed closed and smooth surface $\Gamma$, i.e.,

$$\left\{\begin{array}[]{ll}\partial_{t}u(t,x)-\Delta_{\Gamma}u(t,x)=f(t,x)&\mbox{in}\,\,\,(0,T)\times\Gamma,\\ u(0,x)=u_{0}(x)&\mbox{in}\,\,\,\Gamma,\end{array}\right.$$ (3.9)

and the corresponding semi-discrete finite element scheme:

$$\left\{\begin{array}[]{ll}(\partial_{t}u_{h}(t,x),v_{h})_{\Gamma_{h}}+(\nabla_{\Gamma_{h}}u_{h}(t,x),\nabla_{\Gamma_{h}}v_{h})_{\Gamma_{h}}=(f_{h}(t,x),v_{h})_{\Gamma_{h}},&\forall\,v_{h}\in S_{h}(\Gamma_{h}),\\[5.0pt] u_{h}(0,x)=u_{h,0},\end{array}\right.$$ (3.12)

where $f_{h}=P_{h}(\Gamma_{h})f^{-\ell}$. According to Lemma 3.1, we can equivalently lift the scheme onto $\Gamma$:

$$\left\{\begin{array}[]{ll}(a_{h}(x)\partial_{t}u_{h}^{\ell}(t,x),v_{h}^{\ell})_{\Gamma}+(B_{h}(x)\nabla_{\Gamma}u_{h}^{\ell}(t,x),\nabla_{\Gamma}v_{h}^{\ell})_{\Gamma}=(a_{h}(x)f_{h}^{\ell}(t,x),v_{h}^{\ell})_{\Gamma},&\forall\,v_{h}^{\ell}\in S_{h}(\Gamma),\\[5.0pt] u_{h}^{\ell}(0,x)=u_{h,0}^{\ell}.\end{array}\right.$$ (3.15)

Let $G(t,x,x_{0})$ denote the Green’s function (i.e. the heat kernel) of the parabolic equation (3.9), i.e. $G=G(\cdot,\cdot\,,x_{0})$ is the solution of

$$\left\{\begin{array}[]{ll}\partial_{t}G(\cdot,\cdot\,,x_{0})-\Delta_{\Gamma}G(\cdot,\cdot\,,x_{0})=0&\mbox{in}\,\,\,(0,T)\times\Gamma,\\ G(0,\cdot,x_{0})=\delta_{x_{0}}&\mbox{in}\,\,\,\Gamma.\end{array}\right.$$ (3.18)

The Li–Yau heat kernel estimate on closed manifold $\Gamma$ with Ricci curvature bounded from below (cf. [38, Corollary 3.1]) gives

$$G(t,x,y)\leq c_{0}t^{-\frac{m}{2}}e^{c_{1}t}e^{-\frac{d(x,y)^{2}}{c_{0}t}}\qquad t>0,$$ (3.19)

where the constant $c_{0}>0$ and $c_{1}\geq 0$ only depend on the lower bound of the Ricci curvature. Additionally, $c_{1}=0$ if the Ricci curvature of $\Gamma$ is non-negative everywhere and $c_{1}>0$ for the general case. The self-adjointness of Laplace-Beltrami operator implies the symmetry of the kernel, i.e. $G(t,x,y)=G(t,y,x)$, and the analyticity of the kernel $G(z,x,y)$ on the right complex plane $z\in\{t+is:t>0\}$ (see [11, Lemma 2]). Following [10, Theorem 3.4.8], the following pointwise bound for the complex time heat kernel $G(z,x,y)$ holds.

Lemma 3.2

The analytic complex time heat kernel on surface $\Gamma$ satisfies

$$|G(z,x,y)|\leq c_{0}(\text{\rm Re}(z))^{-\frac{m}{2}}e^{c_{1}\text{\rm Re}(z)}e^{-\text{\rm Re}\frac{d(x,y)^{2}}{c_{0}z}},$$ (3.20)

for all $\text{\rm Re}(z)>0$ and $x,y\in\Gamma$.

Then Cauchy’s integral formula says that for all real time $t>0$

$$\partial_{t}^{k}G(t,x,y)=\frac{k!}{2\pi i}\int_{|z-t|=\frac{t}{2}}\,\,\frac{G(z,x,y)}{(z-t)^{k+1}}{\mathrm{d}}z,$$ (3.21)

which, together with (3.20), yields the following Gaussian pointwise estimate for the time derivatives of Green’s function:

$$|\partial_{t}^{k}G(t,x,x_{0})|\leq\frac{C_{k}}{t^{k+m/2}}e^{C_{k}t}e^{-\frac{d(x-x_{0})^{2}}{C_{k}t}}, \forall\,x,x_{0}\in\Gamma,\,\,\,\forall\,t>0,\,\,k=0,1,2,\dotsc,$$ (3.22)

where the constant $C_{k}$ depends on $c_{0}$ and $c_{1}$ which are related to the Ricci curvature of $\Gamma$.

As constructed in [48, Lemma 2.2], for any $x_{0}\in\bar{K}^{l}\subseteq\Gamma$ (where $K^{l}$ is a lifted open triangle in the triangulation of $\Gamma$), there exists a function $\tilde{\delta}_{x_{0}}\in C_{0}^{\infty}(K^{l})$ such that

$$\chi_{h}(x_{0})=\int_{\Gamma}\chi_{h}\tilde{\delta}_{x_{0}}{\mathrm{d}}x,\quad\forall\,\chi_{h}\in S_{h}(\Gamma),$$

and

$$\|\tilde{\delta}_{x_{0}}\|_{W^{l,p}}\leq Ch^{-l-m(1-1/p)}\qquad\mbox{for}\,\,\,1\leq p\leq\infty,\,\,\,l=0,1,2,\dotsc,$$ (3.23)

$$\sup_{y\in\Gamma}\int_{\Gamma}|\tilde{\delta}_{y}(x)|{\mathrm{d}}x+\sup_{x\in\Gamma}\int_{\Gamma}|\tilde{\delta}_{y}(x)|{\mathrm{d}}y\leq C.$$ (3.24)

Let $\delta_{x_{0}}$ denote the Dirac Delta function centered at $x_{0}\in\Gamma$. In other words, $\int_{\Gamma}\delta_{x_{0}}(y)\varphi(y){\mathrm{d}}y=\varphi(x_{0})$ for arbitrary $\varphi\in C(\Gamma)$. Then the discrete Delta function associated to $\Gamma$

$$\tilde{\delta}_{h,x_{0}}:=P_{h}(\Gamma)\delta_{x_{0}}=P_{h}(\Gamma)\tilde{\delta}_{x_{0}}\qquad\forall x_{0}\in\Gamma$$

decays exponentially away from $x_{0}$ (cf. [48, Lemma 2.3] and [46, Lemma 6.1]):

$$|\tilde{\delta}_{h,x_{0}}(x)|=|P_{h}(\Gamma)\tilde{\delta}_{x_{0}}(x)|\leq Ch^{-m}e^{-\frac{d(x,x_{0})}{Kh}},\quad\forall\,x,x_{0}\in\Gamma.$$ (3.25)

Let ${H}={H}(\cdot,\cdot\,,x_{0})$, $x_{0}\in\Gamma$ be the regularized Green’s function of the parabolic equation on the closed manifold $\Gamma$, defined by

$$\left\{\begin{array}[]{ll}\partial_{t}{H}(\cdot,\cdot\,,x_{0})-\Delta_{\Gamma}{H}(\cdot,\cdot\,,x_{0})=0&\mbox{in}\,\,\,(0,T)\times\Gamma,\\ {H}(0,\cdot,x_{0})=\tilde{\delta}_{x_{0}}&\mbox{in}\,\,\,\Gamma.\end{array}\right.$$ (3.28)

The regularized Green’s function can be represented by

$${H}(t,x,x_{0})=\int_{\Gamma}G(t,y,x)\tilde{\delta}_{x_{0}}(y){\mathrm{d}}y=\int_{\Gamma}G(t,x,y)\tilde{\delta}_{x_{0}}(y){\mathrm{d}}y.$$ (3.29)

From the representation (3.29) and (3.22), one can easily derive that the regularized Green’s function ${H}$ also satisfies the Gaussian pointwise estimate in the region $\max(d(x,x_{0}),\sqrt{t})\geq 2h$:

$$|\partial_{t}^{k}{H}(t,x,x_{0})|\leq\frac{C_{k}}{t^{k+m/2}}e^{C_{k}t}e^{-\frac{d(x-x_{0})^{2}}{C_{k}t}},\qquad\forall\,x,x_{0}\in\Gamma,\,\,\forall\,t>0,$$ (3.30)

with $k=0,1,2,\dotsc$.

Similarly, for any $x_{0}\in\bar{K}\subseteq\Gamma_{h}$ (where $K$ is an open triangle in the triangulation of $\Gamma_{h}$), we denote by $\bar{\delta}_{x_{0}}\in C_{0}^{\infty}(K)$ the regularized Green’s function on $\Gamma_{h}$ (cf. [48, Lemma 2.2]), satisfying the following relations:

$$\chi_{h}(x_{0})=\int_{\Gamma_{h}}\chi_{h}\bar{\delta}_{x_{0}}{\mathrm{d}}x,\quad\forall\,\chi_{h}\in S_{h}(\Gamma_{h}),$$

and

$$\|\bar{\delta}_{x_{0}}\|_{W^{l,p}}\leq Ch^{-l-m(1-1/p)}\qquad\mbox{for}\,\,\,1\leq p\leq\infty,\,\,\,l=0,1,2,\dotsc,$$ (3.31)

$$\sup_{y\in\Gamma_{h}}\int_{\Gamma_{h}}|\bar{\delta}_{y}(x)|{\mathrm{d}}x+\sup_{x\in\Gamma_{h}}\int_{\Gamma_{h}}|\bar{\delta}_{y}(x)|{\mathrm{d}}y\leq C.$$ (3.32)

Let $\bar{\delta}_{h,x_{0}}$ be the discrete Green’s function in $S_{h}(\Gamma_{h})$, defined as

$$\bar{\delta}_{h,x_{0}}:=P_{h}(\Gamma_{h})\delta_{x_{0}}=P_{h}(\Gamma_{h})\bar{\delta}_{x_{0}}\qquad\forall x_{0}\in\Gamma_{h}.$$

By the same argument as in a flat domain (see [46, Lemma 6.1]), it is straightforward to verify that $\bar{\delta}_{h,x_{0}}$ exponentially decays away from $x_{0}$, i.e.,

$$|\bar{\delta}_{h,x_{0}}(x)|=|P_{h}(\Gamma_{h})\bar{\delta}_{x_{0}}(x)|\leq Kh^{-m}e^{-\frac{d(x,x_{0})}{Kh}},\quad\forall\,x,x_{0}\in\Gamma_{h}.$$ (3.33)

This estimate of the discrete Delta function also implies the boundedness of projection $P_{h}(\Gamma_{h})$ in $L^{q}(\Gamma_{h})$, i.e.,

$$\|P_{h}(\Gamma_{h})v\|_{L^{p}(\Gamma_{h})}\leq C\|v\|_{L^{p}(\Gamma_{h})}\quad\forall\,v\in L^{p}(\Gamma_{h}),\,\,\,1\leq p\leq\infty.$$ (3.34)

Let ${H}_{h}={H}_{h}(\cdot,\cdot,x_{0})\in S_{h}(\Gamma_{h}),x_{0}\in\Gamma_{h},$ be the discrete Green’s function associated to (3.12), i.e.,

$$\left\{\begin{array}[]{ll}(\partial_{t}{H}_{h}(t,\cdot,x_{0}),v_{h})_{\Gamma_{h}}+(\nabla_{\Gamma_{h}}{H}_{h}(t,\cdot,x_{0}),\nabla_{\Gamma_{h}}v_{h})_{\Gamma_{h}}=0,&\forall\,v_{h}\in S_{h}(\Gamma_{h}),\\[5.0pt] {H}_{h}(0,\cdot,x_{0})=\bar{\delta}_{h,x_{0}}.\end{array}\right.$$ (3.37)

Then the discrete Green’s function ${H}_{h}(t,x,x_{0})$ is symmetric with respect to $x$ and $x_{0}$. Moreover, in view of the relations in (3.5)–(3.6), the lifted discrete Green’s function ${H}_{h}^{\ell}(t,x,x_{0})={H}_{h}(t,q^{-1}(x),q^{-1}(x_{0}))$, for $x,x_{0}\in\Gamma$, satisfies the following weak formulation on $\Gamma$:

$$\left\{\begin{array}[]{ll}(a_{h}(\cdot)\partial_{t}{H}_{h}^{\ell}(t,\cdot,x_{0}),v_{h}^{\ell})_{\Gamma}+(B_{h}(\cdot)\nabla_{\Gamma}{H}_{h}^{\ell}(t,\cdot,x_{0}),\nabla_{\Gamma}v_{h}^{\ell})_{\Gamma}=0,&\forall\,v_{h}\in S_{h}(\Gamma),\\[5.0pt] {H}_{h}^{\ell}(0,\cdot,x_{0})=\bar{\delta}_{h,q^{-1}(x_{0})}^{\ell}.\end{array}\right.$$ (3.40)

By testing (3.9) and (3.12) with the backward Green’s function and discrete Green’s functions, respectively, the solutions of (3.9) and (3.12) can be represented by

$$u(t,x_{0}) =\int_{\Gamma}G(t,x,x_{0})u_{0}(x){\mathrm{d}}x+\int_{0}^{t}\int_{\Gamma}G(t-s,x,x_{0})f(s,x){\mathrm{d}}x{\mathrm{d}}s,$$ (3.41)

$$u_{h}(t,x_{0}) =\int_{\Gamma_{h}}{H}_{h}(t,x,x_{0})u_{h,0}(x){\mathrm{d}}x+\int_{0}^{t}\int_{\Gamma_{h}}{H}_{h}(t-s,x,x_{0})f_{h}(s,x){\mathrm{d}}x{\mathrm{d}}s$$

$$=\int_{\Gamma}a_{h}(x){H}_{h}^{\ell}(t,x,q(x_{0}))u_{h,0}^{\ell}(x){\mathrm{d}}x+\int_{0}^{t}\int_{\Gamma}a_{h}(x){H}_{h}^{\ell}(t-s,x,q(x_{0}))f_{h}^{\ell}(s,x){\mathrm{d}}x{\mathrm{d}}s,$$ (3.42)

where the last equality uses (3.5) and relation $f_{h}=P_{h}(\Gamma_{h})f^{-\ell}$.

Remark 3.2

The construction of the regularized Green’s function is not unique, since we can indeed choose any smooth weight function in the reference element in the proof of [48, Lemma 2.2]. By definition, $\bar{\delta}_{h,x_{0}}$ and $\tilde{\delta}_{h,x_{0}}$ are the $L^{2}$ projections of the Dirac delta functions onto the closed subspaces $S_{h}(\Gamma_{h})$ and $S_{h}(\Gamma)$ (in the sense of distribution), respectively, and therefore they do not depend on the definitions of the regularized Green’s functions $\bar{\delta}_{x_{0}}$ and $\tilde{\delta}_{x_{0}}$.

In order to see the relation between the regularized Delta functions $\bar{\delta}_{x_{0}}$ on $\Gamma_{h}$ and $\tilde{\delta}_{x_{0}}$ on $\Gamma$, respectively, we can consider a fixed point $x_{0}\in\bar{K}\subset\Gamma_{h}$ and push the regularized delta function $\bar{\delta}_{x_{0}}$ to $\Gamma$ through the distance projection $q$ (which projects any point in a neighborhood of $\Gamma$ onto $\Gamma$). This leads to the following relations:

$$\bar{\delta}_{x_{0}}=(a_{h}^{-1})^{-\ell}\tilde{\delta}_{q(x_{0})}^{-\ell}\qquad\mbox{and}\qquad\tilde{\delta}_{q(x_{0})}=a_{h}\bar{\delta}_{x_{0}}^{\ell}\qquad\mbox{for all}\,\,\,x_{0}\in\Gamma_{h}.$$ (3.43)

This relation can be used to show the following estimate for the discrete Green’s functions.

Lemma 3.3

Let $h\leq h_{0}$, then the discrete Green’s functions satisfy the estimate

$$\|\bar{\delta}_{h,x_{0}}^{\ell}-\tilde{\delta}_{h,q(x_{0})}\|_{L^{p}(\Gamma)}\leq Ch^{k+1}\min\{\|\tilde{\delta}_{q(x_{0})}\|_{L^{p}(\Gamma)},\|\bar{\delta}_{x_{0}}\|_{L^{p}(\Gamma_{h})}\},$$ (3.44)

for all $x_{0}\in\Gamma_{h}$ and $p\in[1,\infty]$.

Proof. By the construction of the regularized delta function on $\Gamma$ and $\Gamma_{h}$,

$$\|\bar{\delta}_{h,x_{0}}^{\ell}-\tilde{\delta}_{h,q(x_{0})}\|_{L^{p}(\Gamma)}=\|(P_{h}(\Gamma_{h})\bar{\delta}_{x_{0}})^{\ell}-P_{h}(\Gamma)\tilde{\delta}_{q(x_{0})}\|_{L^{p}(\Gamma)}$$

$$=\|(P_{h}(\Gamma_{h})\bar{\delta}_{x_{0}})^{\ell}-P_{h}(\Gamma)(a_{h}\bar{\delta}_{x_{0}}^{\ell})\|_{L^{p}(\Gamma)}$$

$$\leq\|(P_{h}(\Gamma_{h})\bar{\delta}_{x_{0}})^{\ell}-P_{h}(\Gamma)\bar{\delta}_{x_{0}}^{\ell}\|_{L^{p}(\Gamma)}+\|P_{h}(\Gamma)(1-a_{h})\bar{\delta}_{x_{0}}^{\ell}\|_{L^{p}(\Gamma)}$$

$$\leq\sup_{v_{h}\in S_{h}(\Gamma),\|v_{h}\|_{L^{p^{\prime}}(\Gamma)}=1}((P_{h}(\Gamma_{h})\bar{\delta}_{x_{0}})^{\ell}-P_{h}(\Gamma)\bar{\delta}_{x_{0}}^{\ell},v_{h})_{\Gamma}+C\|1-a_{h}\|_{L^{\infty}(\Gamma)}\|\bar{\delta}_{x_{0}}^{\ell}\|_{L^{p}(\Gamma)}$$

$$=\sup_{v_{h}\in S_{h}(\Gamma),\|v_{h}\|_{L^{p^{\prime}}(\Gamma)}=1}(\bar{\delta}_{x_{0}},P_{h}(\Gamma_{h})(a_{h}^{-1})^{-\ell}v_{h}-(a_{h}^{-1})^{-\ell}v_{h})_{\Gamma_{h}}+C\|1-a_{h}\|_{L^{\infty}(\Gamma)}\|\bar{\delta}_{x_{0}}^{\ell}\|_{L^{p}(\Gamma)}$$

$$\leq\sup_{v_{h}\in S_{h}(\Gamma),\|v_{h}\|_{L^{p^{\prime}}(\Gamma)}=1}(\bar{\delta}_{x_{0}},P_{h}(\Gamma_{h})(a_{h}^{-1})^{-\ell}v_{h}-v_{h})_{\Gamma_{h}}$$

$$\quad+\sup_{v_{h}\in S_{h}(\Gamma),\|v_{h}\|_{L^{p^{\prime}}(\Gamma)}=1}(\bar{\delta}_{x_{0}},v_{h}-(a_{h}^{-1})^{-\ell}v_{h})_{\Gamma_{h}}+C\|1-a_{h}\|_{L^{\infty}(\Gamma)}\|\bar{\delta}_{x_{0}}^{\ell}\|_{L^{p}(\Gamma)}$$

$$\leq\sup_{v_{h}\in S_{h}(\Gamma),\|v_{h}\|_{L^{p^{\prime}}(\Gamma)}=1}\|P_{h}(\Gamma_{h})((a_{h}^{-1})^{-\ell}-1)v_{h}\|_{L^{p^{\prime}}(\Gamma_{h})}\|\bar{\delta}_{x_{0}}\|_{L^{p}(\Gamma_{h})}$$

$$\quad+\sup_{v_{h}\in S_{h}(\Gamma),\|v_{h}\|_{L^{p^{\prime}}(\Gamma)}=1}\|((a_{h}^{-1})^{-\ell}-1)v_{h}\|_{L^{p^{\prime}}(\Gamma_{h})}\|\bar{\delta}_{x_{0}}\|_{L^{p}(\Gamma_{h})}+Ch^{k+1}\|\bar{\delta}_{x_{0}}^{\ell}\|_{L^{p}(\Gamma)}$$

$$\leq Ch^{k+1}\|\bar{\delta}_{x_{0}}\|_{L^{p}(\Gamma_{h})},$$

where we have used the $L^{\infty}$ estimate of $1-a_{h}$ (Lemma 3.1) and the $L^{p}$ stability of $P_{h}(\Gamma_{h})$ shown in (3.34). Since $\bar{\delta}_{x_{0}}^{\ell}=a_{h}^{-1}\tilde{\delta}_{q(x_{0})}$, we know that $\|\bar{\delta}_{x_{0}}\|_{L^{p}(\Gamma_{h})}$ and $\|\tilde{\delta}_{q(x_{0})}\|_{L^{p}(\Gamma)}$ are equivalent. $\square$

3.4 $W^{1,q}$-stability of Ritz projection

We consider the Ritz projection $R_{h}\colon H^{1}(\Gamma)\rightarrow S_{h}(\Gamma)$ associated to the operator $-\Delta_{\Gamma}+1$ on $\Gamma$, defined as

$$(\nabla_{\Gamma}(u-R_{h}u),\nabla_{\Gamma}\varphi_{h}^{\ell})_{\Gamma}+(u-R_{h}u,\varphi_{h}^{\ell})_{\Gamma}=0\quad\forall\varphi_{h}^{\ell}\in S_{h}(\Gamma).$$

From [13, Theorem 3.2] we know that the Ritz projection defined above is $W^{1,\infty}$ stable. Consequently the $W^{1,q}$-stability of $R_{h}$ for $q\in[1,\infty]$ follows from interpolation (cf. [6, Theorem 3.70]) and duality via Hodge decomposition (cf. [19, Eq. (10.101)]).

3.5 Dyadic decomposition of the domain ${\mathcal{Q}}=(0,1)\times\Gamma$

In the proof of Theorem 2.2, we decompose the domain ${\mathcal{Q}}=\Gamma_{T}$ with $T=1$, i.e. ${\mathcal{Q}}=(0,1)\times\Gamma$, into subdomains, and present estimates of the finite element solutions in each subdomain. The following dyadic decomposition of ${\mathcal{Q}}$ was introduced in [44]. The readers may skip this subsection if they are familiar with such dyadic decompositions.

For any integer $j$, we define $d_{j}=2^{-j}$. For a given $x_{0}\in\Gamma$, we let $J_{1}=1$, $J_{0}=0$ and $J_{*}$ be an integer satisfying $2^{-J_{*}}=C_{*}h$ with $C_{*}\geq 16$ to be determined later. If

$$h<1/(4C_{*}),$$ (3.45)

then

$$2\leq J_{*}=\log_{2}[1/(C_{*}h)]\leq\log_{2}(2+1/h).$$ (3.46)

Let

$$Q_{*}(x_{0})=\{(t,x)\in\Gamma_{T}:\max(d(x,x_{0}),t^{1/2})\leq d_{J_{*}}\},$$

$$\Gamma_{*}(x_{0})=\{x\in\Gamma:d(x,x_{0})\leq d_{J_{*}}\}\,.$$

We define

$$Q_{j}(x_{0})=\{(t,x)\in\Gamma_{T}:d_{j}\leq\max(d(x,x_{0}),t^{1/2})\leq 2d_{j}\} \mbox{for}\,\,\,j\geq 1,$$

$$\Gamma_{j}(x_{0})=\{x\in\Gamma:d_{j}\leq d(x,x_{0})\leq 2d_{j}\} \mbox{for}\,\,\,j\geq 1,$$

$$D_{j}(x_{0})=\{x\in\Gamma:d(x,x_{0})\leq 2d_{j}\} \mbox{for}\,\,\,j\geq 1,$$

and

$$Q_{0}(x_{0})={\mathcal{Q}}\big{\backslash}\big{(}\cup_{j=1}^{J_{*}}Q_{j}(x_{0})\cup Q_{*}(x_{0})\big{)},$$

$$\Gamma_{0}(x_{0})=\Gamma\big{\backslash}\big{(}\cup_{j=1}^{J_{*}}\Gamma_{j}(x_{0})\cup\Gamma_{*}(x_{0})\big{)}.$$

For $j<0$, we simply define $Q_{j}(x_{0})=\Gamma_{j}(x_{0})=\emptyset$. For all integer $j\geq 0$, we define

$$\Gamma_{j}^{\prime}(x_{0})=\Gamma_{j-1}(x_{0})\cup\Gamma_{j}(x_{0})\cup\Gamma_{j+1}(x_{0}),\quad Q_{j}^{\prime}(x_{0})=Q_{j-1}(x_{0})\cup Q_{j}(x_{0})\cup Q_{j+1}(x_{0}),$$

$$\Gamma_{j}^{\prime\prime}(x_{0})=\Gamma_{j-2}(x_{0})\cup\Gamma_{j}^{\prime}(x_{0})\cup\Gamma_{j+2}(x_{0}),\quad Q_{j}^{\prime\prime}(x_{0})=Q_{j-2}(x_{0})\cup Q_{j}^{\prime}(x_{0})\cup Q_{j+2}(x_{0}),$$

$$D_{j}^{\prime}(x_{0})=D_{j-1}(x_{0})\cup D_{j}(x_{0}),\quad D_{j}^{\prime\prime}(x_{0})=D_{j-2}(x_{0})\cup D_{j}^{\prime}(x_{0}).$$

Then we have

$$\Gamma_{T}=\bigcup^{J_{*}}_{j=0}Q_{j}(x_{0})\,\cup Q_{*}(x_{0})\quad\mbox{and}\quad\Gamma=\bigcup^{J_{*}}_{j=0}\Gamma_{j}(x_{0})\,\cup\Gamma_{*}(x_{0}).$$ (3.47)

We refer to $Q_{*}(x_{0})$ as the “innermost” set. We shall write $\sum_{*,j}$ when the innermost set is included and $\sum_{j}$ when it is not. When $x_{0}$ is fixed, if there is no ambiguity, we simply write $Q_{j}=Q_{j}(x_{0})$, $Q_{j}^{\prime}=Q_{j}^{\prime}(x_{0})$, $Q_{j}^{\prime\prime}=Q_{j}^{\prime\prime}(x_{0})$, $\Gamma_{j}=\Gamma_{j}(x_{0})$, $\Gamma_{j}^{\prime}=\Gamma_{j}^{\prime}(x_{0})$ and $\Gamma_{j}^{\prime\prime}=\Gamma_{j}^{\prime\prime}(x_{0})$.

We shall use the notations

$$\|v\|_{k,D}=\biggl{(}\int_{D}\sum_{|\alpha|\leq k}|\partial_{\Gamma}^{\alpha}v|^{2}{\mathrm{d}}x\biggl{)}^{\frac{1}{2}},\qquad{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|v\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{k,Q}=\biggl{(}\int_{Q}\sum_{|\alpha|\leq k}|\partial_{\Gamma}^{\alpha}v|^{2}{\mathrm{d}}x{\mathrm{d}}t\biggl{)}^{\frac{1}{2}},$$ (3.48)

for any subdomains $D\subset\Gamma$ and $Q\subset(0,1)\times\Gamma$. Throughout this paper, we denote by $C$ a generic positive constant that is independent of $h$, $x_{0}$ and $C_{*}$ (until $C_{*}$ is determined in Section 4). To simplify the notations, we also denote $d_{*}=d_{J_{*}}$.

3.6 Local energy estimates

Since the Gagliardo–Nirenberg interpolation inequality still holds on closed manifold (cf. [6, Theorem 3.70]), the standard energy estimates of heat equation (cf. [28, Lemma 2.1 of Chapter III]) are also valid here. Consequently we have the following local estimates of (regularized) Green’s function (cf. [32, Lemma 4.1] and [32, Eq. (4.9)] with $\alpha=1$).

Lemma 3.4

The Green’s function $G$ defined in (3.18) and the regularized Green’s function $H$ defined in (3.28) satisfy the following estimates:

$$d_{j}^{-5+m/2}\|{H}(\cdot,\cdot,x_{0})\|_{L^{\infty}(Q_{j}(x_{0}))}+d_{j}^{-5}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}(Q_{j}(x_{0}))}$$

$$+d_{j}^{-4}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}H^{2}(Q_{j}(x_{0}))}+d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}H^{2}(Q_{j}(x_{0}))}$$

$$+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}H^{2}(Q_{j}(x_{0}))}\leq Cd_{j}^{-m/2-5},$$ (3.49)

$${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}(Q_{j}(x_{0}))}+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\partial_{t}{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}(Q_{j}(x_{0}))}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}(Q_{j}(x_{0}))}$$

$$\quad+d_{j}^{3}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\partial_{tt}{H}(\cdot,\cdot,x_{0})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}(Q_{j}(x_{0}))}\leq Cd_{j}^{-m/2-1},$$ (3.50)

$$\|G(\cdot,\cdot,x_{0})\|_{L^{\infty}H^{2}(\cup_{k\leq j}Q_{k}(x_{0}))}+d_{j}^{2}\|\partial_{t}G(\cdot,\cdot,x_{0})\|_{L^{\infty}H^{2}(\cup_{k\leq j}Q_{k}(x_{0}))}\leq Cd_{j}^{-m/2-2}.$$ (3.51)

Remark 3.3

The spatial high-order counterpart of (3.51), i.e.

$$\|G(\cdot,\cdot,x_{0})\|_{L^{\infty}H^{k}(\cup_{k\leq j}Q_{k}(x_{0}))}\leq Cd_{j}^{-m/2-k},\quad k\in\mathbb{N},$$

follows from a standard induction argument by differentiating the localized equation in the spatial direction. Note that on the surface the surface derivatives do not commute with each other (cf. [8, Lemma 5.1, Item 4]). However the commutator is of lower order and therefore dose not affect the dominant stability in the energy estimates.

We also need the following local energy estimate for the error equation, which can be proved in a similar way as [32, Lemma 5.1]. The proof of Lemma 3.5 can be found in Appendix.

Lemma 3.5

Suppose that $\phi\in L^{2}(0,T;H^{1}_{0}(\Gamma))\cap H^{1}(0,T;L^{2}(\Gamma))$ and $\phi_{h}^{\ell}\in H^{1}(0,T;S_{h}(\Gamma))$ satisfy the equation

$$(\partial_{t}\phi-a_{h}\partial_{t}\phi_{h}^{\ell}),\chi_{h}^{\ell})_{\Gamma}+(\nabla_{\Gamma}\phi-B_{h}\nabla_{\Gamma}\phi_{h}^{\ell},\nabla_{\Gamma}\chi_{h}^{\ell})_{\Gamma}=0,\quad\forall\,\chi_{h}^{\ell}\in S_{h}(\Gamma),\,\,\mbox{a.e.}\,\,t>0,$$ (3.52)

with $\phi(0)=0$ in $\Gamma_{j}^{\prime\prime}$. Then

$${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(\phi-\phi_{h}^{\ell})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\phi-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}$$

$$\leq C\epsilon^{-3}\big{(}I_{j}(\phi_{h}^{\ell}(0))+X_{j}(I_{h}\phi-\phi)+Y_{j}(\phi)+d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\phi-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\big{)}$$

$$\quad+(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)\big{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(\phi-\phi_{h}^{\ell})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\phi-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}\big{)},$$ (3.53)

where

$$I_{j}(\phi_{h}^{\ell}(0))=\|\phi_{h}^{\ell}(0)\|_{1,\Gamma_{j}^{\prime}}+d_{j}^{-1}\|\phi_{h}^{\ell}(0)\|_{\Gamma_{j}^{\prime}}\,,$$

$$X_{j}(I_{h}\phi-\phi)=d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(I_{h}\phi-\phi)\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(I_{h}\phi-\phi)\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}$$

$$\qquad\qquad\qquad\,\,+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}\phi-\phi\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}\phi-\phi\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\,,$$

$$Y_{j}(\phi)=h^{k+1}\bigg{(}d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\partial_{t}\phi\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}\phi\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\phi\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\phi\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\bigg{)}.$$

The positive constant $C$ is independent of $h$, $j$ and $C_{*}$; the norms ${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\cdot\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{k,Q_{j}^{\prime}}$ and ${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\cdot\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{k,\Gamma_{j}^{\prime}}$ are defined in (3.48).

Remark 3.4

The term $Y_{j}(\phi)$ is generated by approximating $\Gamma$ with $\Gamma_{h}$.

4 Discrete maximal regularity on a stationary surface
(Proof of Theorem 2.2)

We shall prove the following key lemma using the local energy estimates in Section 3.6.

Lemma 4.1

The functions ${H}_{h}(t,x,x_{0})$, ${H}(t,x,x_{0})$ and $F(t,x,x_{0}):={H}_{h}^{\ell}(t,x,x_{0})-{H}(t,x,x_{0})$ satisfy

$$\sup_{x_{0}\in\Gamma}\|\partial_{t}F(\cdot,\cdot,x_{0})\|_{L^{1}((0,+\infty)\times\Gamma)}+\sup_{x_{0}\in\Gamma}\|t\partial_{tt}F(\cdot,\cdot,x_{0})\|_{L^{1}((0,+\infty)\times\Gamma)}\leq C,$$ (4.1)

$$\sup_{x_{0}\in\Gamma_{h}}\|\partial_{t}{H}_{h}(t,\cdot,x_{0})\|_{L^{1}(\Gamma_{h})}\leq Ce^{-\lambda_{0}t},\qquad\forall\,t\geq 1,$$ (4.2)

where the constants $C>0$ and $\lambda_{0}\geq 0$ are independent of $h$.

To simplify the notation, in this subsection we continue to omit the dependence of surface $\Gamma(s)$ on $s\in[0,T]$. Additionally, we relax the dependence of the Green’s functions on $x_{0}\in\Gamma$ by denoting

$${H}_{h}^{\ell}(t)={H}_{h}^{\ell}(t,\cdot,x_{0}),\qquad{H}(t)={H}(t,\cdot,x_{0})\quad\mbox{and}\quad F(t)={H}_{h}^{\ell}(t)-{H}(t).$$

Since the coefficient $a_{h}$ and $B_{h}$ are bounded everywhere, the same proofs in [32, Section 4.3–4.4] using Lemma 4.1 can be carried over almost verbatim here. For example, by differentiating (3.42) in time, we obtain

$$\partial_{t}u_{h}(t,q^{-1}(x_{0})) =\int_{0}^{t}\int_{\Gamma}a_{h}(x)\partial_{t}F(t-s,x,x_{0})f_{h}^{\ell}(s,x){\mathrm{d}}x{\mathrm{d}}s$$

$$\quad+\int_{0}^{t}\int_{\Gamma}a_{h}(x)\partial_{t}{H}(t-s,x,x_{0})f_{h}^{\ell}(s,x){\mathrm{d}}x{\mathrm{d}}s+a_{h}(x_{0})f_{h}^{\ell}(t,x_{0})$$

$$=:\mathcal{M}_{h}(a_{h}f_{h}^{\ell})+\mathcal{K}_{h}(a_{h}f_{h}^{\ell})+a_{h}f_{h}^{\ell},\quad x_{0}\in\Gamma,$$

where $\mathcal{M}_{h}$ and $\mathcal{K}_{h}$ are the linear operators whose kernel are $\partial_{t}F$ and $\partial_{t}{H}$ respectively. Therefore, by using Lemma 4.1 and the arguments in [32, Section 4.3–4.4], we can obtain the following result for $p,q\in[1,\infty)$:

$$\|\partial_{t}u_{h}^{\ell}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\Gamma))}\leq C\|a_{h}f_{h}^{\ell}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\Gamma))}\leq C\|f_{h}^{\ell}\|_{L^{p}({\mathbb{R}}_{+};L^{q}(\Gamma))}.$$

This proves the maximal regularity result in Theorem 2.2 (by choosing $\Gamma=\Gamma(s)$ and $\Gamma_{h}=\Gamma_{h}(s)$ for any fixed $s\in[0,T]$).

Remark 4.1

Analogous to [32, Section 4.2], we also have the analyticity and maximum-norm stability (or discrete maximum principle, for the continuous maximum principle on closed manifold $\Gamma$; see [37, Theorem 15.2]) results for the discrete heat semigroup on $\Gamma_{h}$: The semi-discrete finite element solution $u_{h}=E_{h}(t)u_{h,0}$ of equation

$$\left\{\begin{aligned} &(\partial_{t}u_{h}(t),\varphi_{h})_{\Gamma_{h}}+(\nabla_{\Gamma_{h}}u_{h}(t),\nabla_{\Gamma_{h}}\varphi_{h})_{\Gamma_{h}}=0&&\forall\,\varphi_{h}\in S_{h},\\ &u_{h}(0)=u_{h,0}&&\mbox{on}\,\,\,\Gamma,\end{aligned}\right.$$

satisfies the following estimates for $q\in[1,\infty)$:

$$\sup_{t>0}(\|E_{h}(t)u_{h,0}\|_{L^{q}(\Gamma_{h})}+t\|\partial_{t}E_{h}(t)u_{h,0}\|_{L^{q}(\Gamma_{h})})\leq C\|u_{h,0}\|_{L^{q}(\Gamma_{h})},$$

$$\sup_{t>0}\|E_{h}(t)u_{h,0}\|_{L^{\infty}(\Gamma_{h})}\leq C\|u_{h,0}\|_{L^{\infty}(\Gamma_{h})}.$$

In order to prove Lemma 4.1, we apply the local energy estimate in Lemma 3.5 to estimate $\|\partial_{t}F\|_{L^{1}((0,+\infty)\times\Gamma)}+\|t\partial_{tt}F\|_{L^{1}((0,+\infty)\times\Gamma)}$. The estimation consists of two parts: The first part concerns estimates for $t\in(0,1)$, and the second part concerns estimates for $t\geq 1$, which is a consequence of the smoothing property of parabolic equations.

Part I.    First, we present estimates in the domain ${\mathcal{Q}}=(0,1)\times\Gamma$ with the restriction $h<1/(4C_{*})$; see (3.45). In this case, the basic energy estimate gives

$$\|\partial_{t}{H}\|_{L^{2}({\mathcal{Q}})}+\|\partial_{t}{H}_{h}^{\ell}\|_{L^{2}({\mathcal{Q}})}\leq C(\|{H}(0)\|_{H^{1}}+\|{H}_{h}^{\ell}(0)\|_{H^{1}})\leq Ch^{-1-m/2},$$ (4.4)

$$\|{H}\|_{L^{\infty}L^{2}({\mathcal{Q}})}+\|{H}_{h}^{\ell}\|_{L^{\infty}L^{2}({\mathcal{Q}})}\leq C(\|{H}(0)\|_{L^{2}}+\|{H}_{h}^{\ell}(0)\|_{L^{2}})\leq Ch^{-m/2},$$ (4.5)

$$\|\nabla_{\Gamma}{H}\|_{L^{2}({\mathcal{Q}})}+\|\nabla_{\Gamma}{H}_{h}^{\ell}\|_{L^{2}({\mathcal{Q}})}\leq C(\|{H}(0)\|_{L^{2}}+\|{H}_{h}^{\ell}(0)\|_{L^{2}})\leq Ch^{-m/2},$$ (4.6)

$$\|\partial_{tt}{H}\|_{L^{2}({\mathcal{Q}})}+\|\partial_{tt}{H}_{h}^{\ell}\|_{L^{2}({\mathcal{Q}})}\leq C(\|\Delta_{\Gamma}{H}(0)\|_{H^{1}}+\|\Delta_{\Gamma,h}{H}_{h}^{\ell}(0)\|_{H^{1}})\leq Ch^{-3-m/2},$$ (4.7)

$$\|\nabla_{\Gamma}\partial_{t}{H}\|_{L^{2}({\mathcal{Q}})}+\|\nabla_{\Gamma}\partial_{t}{H}_{h}^{\ell}\|_{L^{2}({\mathcal{Q}})}\leq C(\|\Delta_{\Gamma}{H}(0)\|_{L^{2}}+\|\Delta_{\Gamma,h}{H}_{h}^{\ell}(0)\|_{L^{2}})\leq Ch^{-2-m/2},$$ (4.8)

where we have used (3.23) and (3.25) to estimate ${H}(0)$ and ${H}_{h}(0)$, respectively. Hence, we have

$${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}+\||{H}_{h}^{\ell}\||_{Q_{*}}\leq Cd_{*}\|{H}\|_{L^{\infty}L^{2}(Q_{*})}+Cd_{*}\|{H}_{h}^{\ell}\|_{L^{\infty}L^{2}(Q_{*})}\leq Cd_{*}h^{-m/2}\leq CC_{*}h^{1-m/2}.$$ (4.9)

Since the volume of $Q_{j}$ is $Cd_{j}^{2+m}$, we can decompose $\|\partial_{t}F\|_{L^{1}(\mathcal{Q})}+\|t\partial_{tt}F\|_{L^{1}(\mathcal{Q})}$ in the following way:

$$\|\partial_{t}F\|_{L^{1}(\mathcal{Q})}+\|t\partial_{tt}F\|_{L^{1}(\mathcal{Q})}$$

$$\leq\|\partial_{t}F\|_{L^{1}(Q_{*})}+\|t\partial_{tt}F\|_{L^{1}(Q_{*})}+\sum_{j}\big{(}\|\partial_{t}F\|_{L^{1}(Q_{j})}+\|t\partial_{tt}F\|_{L^{1}(Q_{j})}\big{)}$$

$$\leq Cd_{*}^{1+m/2}\big{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}+d_{*}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}\big{)}+\sum_{j}Cd_{j}^{1+m/2}\big{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}\big{)}$$

$$\leq CC_{*}^{3+m/2}+{\mathscr{K}},$$ (4.10)

where we have used (4.4), (4.6) and (4.7) to estimate

$$Cd_{*}^{1+m/2}\big{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}+d_{*}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}\big{)},$$

and introduced the following notation:

$${\mathscr{K}}: =\sum_{j}d_{j}^{1+m/2}(d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}).$$ (4.11)

It remains to estimate ${\mathscr{K}}$. To this end, in view of (3.28) and (3.40), we set “$\phi_{h}=H_{h}$, $\phi=H$, $\phi_{h}(0)=P_{h}(\Gamma_{h})\bar{\delta}_{q^{-1}(x_{0})}$ and $\phi(0)=\tilde{\delta}_{x_{0}}$” and “$\phi_{h}=\partial_{t}H_{h}$, $\phi=\partial_{t}H$, $\phi_{h}^{\ell}(0)=\Delta_{\Gamma,h}(P_{h}(\Gamma_{h})\bar{\delta}_{q^{-1}(x_{0})})^{\ell}$ and $\phi(0)=\Delta_{\Gamma}\tilde{\delta}_{x_{0}}$” in Lemma 3.5, respectively. Then we obtain

$$d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}} \leq C\epsilon^{-3}(\widehat{I_{j}}+\widehat{X_{j}}+\widehat{Y_{j}}+d^{-2}_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q^{\prime}_{j}})$$ (4.12)

$$\quad+(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)\big{(}d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\big{)}\qquad\quad\,\,$$

and

$$d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}} \leq C\epsilon^{-3}(\overline{I_{j}}+\overline{X_{j}}+\overline{Y_{j}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q^{\prime}_{j}})$$ (4.13)

$$\quad+(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)\big{(}d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\big{)}\,,$$

respectively. By using interpolation error estimate, (3.33) (exponential decay of $P_{h}(\Gamma_{h})\bar{\delta}_{q^{-1}(x_{0})}$) and Lemma 3.4 (local estimates of regularized Green’s function), we have

$$\widehat{I_{j}} =\|P_{h}(\Gamma_{h})\bar{\delta}_{q^{-1}(x_{0})}\|_{1,\Gamma_{j}^{\prime}}+d_{j}^{-1}\|P_{h}(\Gamma_{h})\bar{\delta}_{q^{-1}(x_{0})}\|_{\Gamma_{j}^{\prime}}\leq Ch^{2}d_{j}^{-3-m/2},$$ (4.14)

$$\widehat{X_{j}} =d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(I_{h}{H}-{H})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(I_{h}{H}-{H})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}$$

$$\quad\,+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}{H}-{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}{H}-{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}$$

$$\leq Cd_{j}h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}H^{2}(Q_{j}^{\prime\prime})}+Cd_{j}^{-1}h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}H^{2}(Q_{j}^{\prime\prime})}$$

$$\leq Chd_{j}^{-2-m/2},$$ (4.15)

$$\widehat{Y_{j}} =h^{k+1}\bigg{(}d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\bigg{)}$$

$$\leq Ch^{k+1}d_{j}^{-1-m/2},$$ (4.16)

and

$$\overline{I_{j}} =d_{j}^{2}\|\Delta_{\Gamma,h}(P_{h}(\Gamma_{h})\bar{\delta}_{q^{-1}(x_{0})})^{\ell}\|_{1,\Gamma_{j}^{\prime}}+d_{j}\|\Delta_{\Gamma,h}(P_{h}(\Gamma_{h})\bar{\delta}_{q^{-1}(x_{0})})^{\ell}\|_{\Gamma_{j}^{\prime}}\leq Ch^{2}d_{j}^{-3-m/2},$$ (4.17)

$$\overline{X_{j}} =d_{j}^{3}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}\partial_{tt}{H}-\partial_{tt}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}\partial_{tt}{H}-\partial_{tt}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}$$

$$\quad\,+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}\partial_{t}{H}-\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|I_{h}\partial_{t}{H}-\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}$$

$$\leq C(d_{j}^{3}h+d_{j}^{2}h^{2}){\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}H^{2}(Q_{j}^{\prime\prime})}+C(d_{j}h+h^{2}){\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}H^{2}(Q_{j}^{\prime\prime})}$$

$$\leq Chd_{j}^{-2-m/2},$$ (4.18)

$$\overline{Y_{j}} =h^{k+1}\bigg{(}d_{j}^{3}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\partial_{tt}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\bigg{)}$$

$$\leq Ch^{k+1}d_{j}^{-1-m/2}.$$ (4.19)

By substituting (4.12)-(4) into the expression of ${\mathscr{K}}$ in (4.11), we have

$${\mathscr{K}} =\sum_{j}d_{j}^{1+m/2}(d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}})$$

$$\leq C\sum_{j}d_{j}^{1+m/2}\epsilon^{-3}\big{(}h^{2}d_{j}^{-3-m/2}+hd_{j}^{-2-m/2}+d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q^{\prime}_{j}}\big{)}$$

$$\quad+\sum_{j}(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)d_{j}^{1+m/2}(d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}})$$

$$\leq C+C\epsilon^{-3}\sum_{j}d_{j}^{-1+m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}$$

$$\quad+\sum_{j}(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)d_{j}^{1+m/2}(d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}).$$ (4.20)

Since ${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}\leq C({\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j-1}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j+1}})$, we can convert the $Q_{j}^{\prime}$-norm in the inequality above to the $Q_{j}$-norm:

$${\mathscr{K}} \leq C+C\epsilon^{-3}d_{*}^{-1+m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}+C\epsilon^{-3}\sum_{j}d_{j}^{-1+m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}$$

$$\quad+\sum_{j}(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)d_{*}^{1+m/2}(d_{*}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{*}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}+d_{*}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{*}}+d_{*}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}})$$

$$\quad+\sum_{j}(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)d_{j}^{1+m/2}(d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+d_{j}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}+d_{j}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}})$$

$$\leq C_{\epsilon}+C_{\epsilon}C_{*}^{3+m/2}+C_{\epsilon}\sum_{j}d_{j}^{-1+m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+C(C_{\epsilon}C_{*}^{-1/2}+CC_{*}^{-1}+\epsilon){\mathscr{K}},$$ (4.21)

where we have used $d_{j}\geq C_{*}h$ and (4.4)-(4.8) to estimate

$${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{*}},\qquad{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}},\qquad{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{*}},\qquad\text{and}\qquad{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{tt}F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}},$$

and used the expression of ${\mathscr{K}}$ in (4.11) to bound the terms involving $Q_{j}$. Since $\epsilon_{*}=\epsilon+\epsilon^{-1}/C_{*}$, we can make $\epsilon_{*}$ sufficiently small by first choosing $\epsilon$ small enough and then choosing $C_{*}$ large enough ($\epsilon$ can be fixed now and $C_{*}$ will be determined later). Then the last term on the right-hand side of (4) can be absorbed by the left-hand side. Therefore, we obtain

$${\mathcal{K}} \leq C+CC_{*}^{3+m/2}+\sum_{j}d_{j}^{-1+m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}.$$ (4.22)

It remains to estimate ${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}$. We apply the parabolic duality argument: Let $w$ be the solution of the following backward parabolic equation on the domain $\mathcal{Q}$:

$$\left\{\begin{array}[]{ll}-\partial_{t}w-\Delta_{\Gamma}w=v,\\ w(1)=0,\end{array}\right.$$

with ${\rm supp}(v)\subseteq Q_{j}$ and $\|v\|_{L^{2}(Q_{j})}=1$. After testing the above equation by $F$, we get

$$[F,v]=(F(0),w(0))+[F_{t},w]+[\nabla_{\Gamma}F,\nabla_{\Gamma}w],$$ (4.23)

where

$$(F(0),w(0))=(\bar{\delta}_{h,q^{-1}(x_{0})}^{\ell}-\tilde{\delta}_{x_{0}},w(0))$$

$$=(P_{h}\tilde{\delta}_{x_{0}}-\tilde{\delta}_{x_{0}},w(0)-I_{h}w(0))+(\bar{\delta}_{h,q^{-1}(x_{0})}^{\ell}-\tilde{\delta}_{h,x_{0}},w(0))$$

$$=(P_{h}\tilde{\delta}_{x_{0}},w(0)-I_{h}w(0))_{\Gamma_{j}^{\prime\prime}}+(P_{h}\tilde{\delta}_{x_{0}}-\tilde{\delta}_{x_{0}},w(0)-I_{h}w(0))_{(\Gamma_{j}^{\prime\prime})^{c}}+(\bar{\delta}_{h,q^{-1}(x_{0})}^{\ell}-\tilde{\delta}_{h,x_{0}},w(0))$$

$$=:\mathcal{I}_{1}+\mathcal{I}_{2}+\mathcal{I}_{3}.$$ (4.24)

Both $\mathcal{I}_{1}$ and $\mathcal{I}_{2}$ can be estimated in the same way as [32, Eq. (5.21), (5.23) and (5.24)]. The only difference is the replacement of flat domain $\Omega$ in [32] by surface $\Gamma$ here. These estimates of $\mathcal{I}_{1}$ and $\mathcal{I}_{2}$ can be written as follows:

$$|\mathcal{I}_{1}|\leq Ch^{2}d_{j}^{-1-m/2}\quad\mbox{and}\quad\,\,|\mathcal{I}_{2}| \leq Ch^{2}\|w(0)\|_{W^{2,\infty}((\Gamma_{j}^{\prime})^{c})},$$

$$\leq Ch^{2}\sup_{y\in\Gamma}\|G(\cdot,\cdot,y)\|_{L^{\infty}W^{2,\infty}(\cup_{k\leq j+\log_{2}C}Q_{k})}\|v\|_{L^{1}(Q_{j})},$$

$$\leq Ch^{2}d_{j}^{-2-m}\|v\|_{L^{1}(Q_{j})},$$

$$\leq Ch^{2}d_{j}^{-1-m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|v\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}(Q_{j})},$$

where in the second-to-last inequality we have used the following pointwise estimate of Green’s function (cf. (3.22)):

$$\sup_{y\in\Gamma}\|G(\cdot,\cdot,y)\|_{L^{\infty}W^{2,\infty}(\cup_{k\leq j+\log_{2}C}Q_{k})}$$

$$\leq C\sup_{y\in\Gamma}\|G(\cdot,\cdot,y)\|_{L^{\infty}L^{\infty}(\cup_{k\leq j+\log_{2}C}Q_{k})}^{1-\theta}\|G(\cdot,\cdot,y)\|_{L^{\infty}H^{M}(\cup_{k\leq j+\log_{2}C}Q_{k})}^{\theta}$$

$$\leq Cd_{j}^{-m(1-\theta)}d_{j}^{-(M+m/2)\theta}=Cd_{j}^{-2-m},$$

where in the first and second inequality we have used the Gagliardo–Nirenberg interpolation inequality with $\theta=\frac{4}{2M-m}$ for any integer $M>\frac{m}{2}$ (cf. [6, Theorem 3.70]) and the local estimates of Green’s function (Lemma 3.4 and Remark 3.3) respectively.

Similarly, by using Lemma 3.3 and (3.22), the argument in [32] leads to the following estimate:

$$|\mathcal{I}_{3}| \leq\|\bar{\delta}_{h,q^{-1}(x_{0})}^{\ell}-\tilde{\delta}_{h,x_{0}}\|_{L^{1}(\Gamma)}\|w(0)\|_{L^{\infty}(\Gamma)}$$

$$\leq Ch^{k+1}\sup_{y\in\Gamma}\|G(\cdot,\cdot,y)\|_{L^{\infty}L^{\infty}(\cup_{k\leq j+\log_{2}C}Q_{k})}\|v\|_{L^{1}(Q_{j})}$$

$$\leq Ch^{k+1}d_{j}^{-m}\|v\|_{L^{1}(Q_{j})}$$

$$\leq Ch^{k+1}d_{j}^{1-m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|v\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{L^{2}(Q_{j})}.$$

We decompose the second and third terms on the right-hand side of (4.23) as follows:

$$[F_{t},w]+[\nabla_{\Gamma}F,\nabla_{\Gamma}w]$$

$$=[\partial_{t}({H}_{h}^{\ell}-{H}),w]+[\nabla_{\Gamma}({H}_{h}^{\ell}-{H}),\nabla_{\Gamma}w]$$

$$=[a_{h}\partial_{t}{H}_{h}^{\ell}-\partial_{t}{H},w]+[B_{h}\nabla_{\Gamma}{H}_{h}^{\ell}-\nabla_{\Gamma}{H},\nabla_{\Gamma}w]$$

$$\quad+[(1-a_{h})\partial_{t}{H}_{h}^{\ell},w]+[(1-B_{h})\nabla_{\Gamma}{H}_{h}^{\ell},\nabla_{\Gamma}w]$$

$$=[a_{h}\partial_{t}{H}_{h}^{\ell}-\partial_{t}{H},(w-I_{h}w)]+[B_{h}\nabla_{\Gamma}{H}_{h}^{\ell}-\nabla_{\Gamma}{H},\nabla_{\Gamma}(w-I_{h}w)]$$

$$\quad+[(1-a_{h})\partial_{t}{H},w]+[(1-B_{h})\nabla_{\Gamma}{H},\nabla_{\Gamma}w]$$

$$\quad+[(1-a_{h})\partial_{t}({H}_{h}^{\ell}-{H}),w]+[(1-B_{h})\nabla_{\Gamma}({H}_{h}^{\ell}-{H}),\nabla_{\Gamma}w]$$

$$=[a_{h}\partial_{t}({H}_{h}^{\ell}-{H}),(w-I_{h}w)]+[B_{h}\nabla_{\Gamma}({H}_{h}^{\ell}-{H}),\nabla_{\Gamma}(w-I_{h}w)]$$

$$\quad+[(a_{h}-1)\partial_{t}{H},(w-I_{h}w)]+[(B_{h}-1)\nabla_{\Gamma}{H},\nabla_{\Gamma}(w-I_{h}w)]$$

$$\quad+[(1-a_{h})\partial_{t}{H},w]+[(1-B_{h})\nabla_{\Gamma}{H},\nabla_{\Gamma}w]$$

$$\quad+[(1-a_{h})\partial_{t}({H}_{h}^{\ell}-{H}),w]+[(1-B_{h})\nabla_{\Gamma}({H}_{h}^{\ell}-{H}),\nabla_{\Gamma}w]$$

$$=:\sum_{i=4}^{11}\mathcal{I}_{i}.$$ (4.25)

By [32, Eq. (5.26)], we know

$$|\mathcal{I}_{4}|+|\mathcal{I}_{5}|\leq C\sum_{*,i}(h^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\|w\|_{L^{2}H^{2}(Q_{i}^{\prime})}.$$ (4.26)

Using the interpolation error estimate,

$$|\mathcal{I}_{6}|+|\mathcal{I}_{7}| \leq Ch^{k+1}\sum_{*,i}(h^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\|w\|_{L^{2}H^{2}(Q_{i}^{\prime})},$$

and from Lemma 3.1, Sobolev embedding and Hölder’s inequality, it holds

$$|\mathcal{I}_{8}|+|\mathcal{I}_{9}| \leq Ch^{k+1}\sum_{*,i}(d_{i}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+d_{i}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\|w\|_{L^{2}H^{2}(Q_{i}^{\prime})},$$

$$|\mathcal{I}_{10}|+|\mathcal{I}_{11}| \leq Ch^{k+1}\sum_{*,i}(d_{i}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+d_{i}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\|w\|_{L^{2}H^{2}(Q_{i}^{\prime})}.$$

According to the estimates above and [32, Eqs. (5.31)–(5.32)] (with $\alpha=1$ therein), we get

$${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}} \leq Ch^{2}d_{j}^{-1-m/2}+Ch^{k+1}\sum_{*,i}(d_{i}^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+d_{i}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|{H}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}$$

$$\quad+C\sum_{*,i}(h^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}$$

$$\leq Ch^{2}d_{j}^{-1-m/2}+Ch^{k+1}\sum_{*,i}d_{i}^{1-m/2}\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}$$

$$\quad+C\sum_{*,i}(h^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}.$$ (4.27)

Note that (cf. [32, Eq. (5.33)] with $\alpha=1$ therein)

$$\sum_{j}d_{j}^{-1+m/2}\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}\leq Cd_{i}^{-1-m/2}.$$ (4.28)

By substituting (4)–(4.28) into (4.22) we obtain

$${\mathscr{K}} \leq C+CC_{*}^{3+m/2}+\sum_{j}Cd_{j}^{-1+m/2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}$$

$$\leq C+CC_{*}^{3+m/2}+C\sum_{j}\left(\frac{h}{d_{j}}\right)^{2}\qquad\mbox{(here we have used \eqref{FL2Qj})}$$

$$\quad+C\sum_{j}d_{j}^{-1+m/2}h^{k+1}\sum_{*,i}d_{i}^{1-m/2}\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}$$

$$\quad+C\sum_{j}d_{j}^{-1+m/2}\sum_{*,i}(h^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}.$$

Consequently,

$${\mathscr{K}} \leq C+CC_{*}^{3+m/2}+CC_{*}^{-2}\qquad\mbox{(here we exchange the order of summation)}$$

$$\quad+C\sum_{*,i}h^{k+1}d_{i}^{1-m/2}\sum_{j}d_{j}^{-1+m/2}\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}$$

$$\quad+C\sum_{*,i}(h^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})\sum_{j}d_{j}^{-1+m/2}\Big{(}\frac{\min(d_{i},d_{j})}{\max(d_{i},d_{j})}\Big{)}$$

$$\leq C+CC_{*}^{3+m/2}+CC_{*}^{-2}+C\sum_{*,i}h^{k+1}$$

$$\quad+C\sum_{*,i}(h^{2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+h{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}})d_{i}^{\frac{m}{2}-1}\qquad\mbox{(here we use \eqref{lalnu})}$$

$$=C+CC_{*}^{3+m/2}+CC_{*}^{-2}+C(\log C_{*}+\log\frac{1}{h})h^{k+1}$$

$$\quad+C\sum_{*,i}d_{i}^{1+m/2}\Big{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}\Big{(}\frac{h}{d_{i}}\Big{)}^{2}+d_{i}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}}\Big{(}\frac{h}{d_{i}}\Big{)}\Big{)}$$

$$\leq C+CC_{*}^{3+m/2}+CC_{*}^{-2}+C(\log C_{*}+\log\frac{1}{h})h^{k+1}+Cd_{*}^{1+m/2}\left({\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{*}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{*}}\right)$$

$$\quad+C\sum_{i}d_{i}^{1+m/2}\left({\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F_{t}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{i}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|F\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{i}}\right)\bigg{(}\frac{h}{d_{i}}\bigg{)}$$

$$\leq C+CC_{*}^{3+m/2}+CC_{*}^{-2}+C\Big{(}\log C_{*}+\log\frac{1}{h}\Big{)}h^{k+1}+CC_{*}^{-1}{\mathscr{K}}.$$

By choosing $C_{*}$ to be large enough ($C_{*}$ is determined now), the term $CC_{*}^{-1}{\mathscr{K}}$ will be absorbed by the left-hand side of the inequality above. In this case, the inequality above implies

$${\mathscr{K}}\leq C.$$ (4.29)

Substituting the last inequality into (4) yields

$$\|\partial_{t}F\|_{L^{1}((0,1)\times\Gamma)}+\|t\partial_{tt}F\|_{L^{1}((0,1)\times\Gamma)} \leq C.$$ (4.30)

Part II.    For $t\in[1,+\infty)$, we first define the finite element space with vanishing average: $\mathring{S}_{h}(\Gamma_{h})=\{\varphi_{h}\in S_{h}(\Gamma_{h}):\int_{\Gamma_{h}}\varphi_{h}=0\}$ and $\mathring{S}_{h}(\Gamma)=\{\varphi_{h}^{\ell}\in S_{h}(\Gamma):\int_{\Gamma}\varphi_{h}^{\ell}=0\}$ on which we can define the inverse of the discrete Laplacian operators $\Delta_{\Gamma_{h},h}^{-1}$ and $\Delta_{\Gamma,h}^{-1}$ respectively.

Lemma 4.2

For any $f_{h}\in\mathring{S}_{h}(\Gamma_{h})$, we define $f_{h}^{\rm avg,\ell}=f_{h}^{\ell}-|\Gamma|^{-1}\int_{\Gamma}f_{h}^{\ell}\in\mathring{S}_{h}(\Gamma)$. Then we have the following properties

$$\|(\Delta_{\Gamma_{h},h}^{-1}f_{h})^{\ell}\|_{L^{2}(\Gamma)} \sim\|\Delta_{\Gamma,h}^{-1}f_{h}^{\rm avg,\ell}\|_{L^{2}(\Gamma)},$$

$$\|\Delta_{\Gamma,h}^{-1}f_{h}^{\rm avg,\ell}\|_{L^{2}(\Gamma)} \leq C\|f_{h}^{\rm avg,\ell}\|_{L^{1}(\Gamma)}\leq 2C\|f_{h}^{\ell}\|_{L^{1}(\Gamma)}.$$

Proof. The first result is a consequence of consistency of $\Gamma$ and $\Gamma_{h}$. We define $u_{h}=\Delta_{\Gamma_{h},h}^{-1}f_{h}\in\mathring{S}_{h}(\Gamma_{h})$ and $\tilde{u}_{h}^{\ell}=\Delta_{\Gamma,h}^{-1}f_{h}^{\rm avg,\ell}\in\mathring{S}_{h}(\Gamma)$ with $\int_{\Gamma_{h}}u_{h}=\int_{\Gamma}\tilde{u}_{h}^{\ell}=0$ and define

$$u_{h}^{\rm avg,\ell}=u_{h}^{\ell}-|\Gamma|^{-1}\int_{\Gamma}u_{h}^{\ell}\in\mathring{S}_{h}(\Gamma).$$

Since $\int_{\Gamma}u_{h}^{\ell}=\int_{\Gamma_{h}}(a_{h}^{-1})^{-\ell}u_{h}=\int_{\Gamma_{h}}((a_{h}^{-1})^{-\ell}-1)u_{h}$ and $|(a_{h}^{-1})^{-\ell}-1|\leq Ch^{k+1}$, it follows that

$$|u_{h}^{\rm avg,\ell}-u_{h}^{\ell}|\leq Ch^{k+1}\|u_{h}\|_{L^{1}(\Gamma_{h})},$$ (4.31)

and, similarly

$$|f_{h}^{\rm avg,\ell}-f_{h}^{\ell}|\leq Ch^{k+1}\|f_{h}\|_{L^{1}(\Gamma_{h})}.$$ (4.32)

By definition,

$$(\nabla_{\Gamma_{h}}u_{h},\nabla_{\Gamma_{h}}v_{h})_{\Gamma_{h}} =(f_{h},v_{h})_{\Gamma_{h}}\quad\forall v_{h}\in S_{h},$$

$$(\nabla_{\Gamma}\tilde{u}_{h}^{\ell},\nabla_{\Gamma}v_{h}^{\ell})_{\Gamma} =(f_{h}^{\rm avg,\ell},v_{h}^{\ell})_{\Gamma}\quad\forall v_{h}\in S_{h}.$$

Applying change of variables and subtraction, we obtain

$$(\nabla_{\Gamma}(u_{h}^{\ell}-\tilde{u}_{h}^{\ell}),\nabla_{\Gamma}v_{h}^{\ell})_{\Gamma} =-((B_{h}-1)\nabla_{\Gamma}\tilde{u}_{h}^{\ell},\nabla_{\Gamma}v_{h}^{\ell})_{\Gamma}+((a_{h}-1)f_{h}^{\ell},v_{h}^{\ell})_{\Gamma}+((f_{h}^{\ell}-f_{h}^{\rm avg,\ell}),v_{h}^{\ell})_{\Gamma}.$$

Then we test with $v_{h}^{\ell}=u_{h}^{\rm avg,\ell}-\tilde{u}_{h}^{\ell}\in\mathring{S}_{h}(\Gamma)$ and use (4.31)–(4.32), Lemma 3.1, Poincaré inequality and the inverse inequality on $\Gamma$ to derive

$$\|\nabla_{\Gamma}(u_{h}^{\rm avg,\ell}-\tilde{u}_{h}^{\ell})\|_{L^{2}(\Gamma)}=\|\nabla_{\Gamma}(u_{h}^{\ell}-\tilde{u}_{h}^{\ell})\|_{L^{2}(\Gamma)}$$

$$\leq Ch^{k+1}\|\nabla_{\Gamma}\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)}+Ch^{k+1}\|f_{h}^{\ell}\|_{L^{2}(\Gamma)}$$

$$\leq Ch^{k+1}\|\nabla_{\Gamma}\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)}+Ch^{k+1}\min\{\|f_{h}^{\ell}\|_{L^{2}(\Gamma)},\|f_{h}^{\rm avg,\ell}\|_{L^{2}(\Gamma)}\}$$

$$\leq Ch^{k}\|\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)}+Ch^{k-1}\min\{\|u_{h}^{\ell}\|_{L^{2}(\Gamma)},\|\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)}\},$$

and from the triangle inequality and Poincaré inequality again

$$\|u_{h}^{\ell}-\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)} \leq\|u_{h}^{\rm avg,\ell}-u_{h}^{\ell}\|_{L^{2}(\Gamma)}+C\|\nabla_{\Gamma}(u_{h}^{\rm avg,\ell}-\tilde{u}_{h}^{\ell})\|_{L^{2}(\Gamma)}$$

$$\leq Ch^{k+1}\|u_{h}^{\ell}\|_{L^{1}(\Gamma)}+Ch^{k}\|\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)}+Ch^{k-1}\min\{\|u_{h}^{\ell}\|_{L^{2}(\Gamma)},\|\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)}\}.$$

Therefore $\|\tilde{u}_{h}^{\ell}\|_{L^{2}(\Gamma)}\sim\|u_{h}^{\ell}\|_{L^{2}(\Gamma)}$ , and this proves the first result.

The second result follows from [7, Lemma 4.3] where we use the elliptic regularity theory and the $W^{1,q}$-stability of Ritz projection on the smooth surface $\Gamma$ (see Section 3.4). $\square$

Denote by $e^{t\Delta_{\Gamma_{h},h}}$ the analytical semigroup generated by the linear operator $\Delta_{\Gamma_{h},h}\colon\mathring{S}_{h}(\Gamma_{h})\rightarrow\mathring{S}_{h}(\Gamma_{h})$. Then, by the previous lemma with $f_{h}$ replaced by $\bar{\delta}_{h,x_{0}}^{\rm avg}$ therein for any fixed $x_{0}\in\Gamma_{h}$, if we define $\bar{\delta}_{h,x_{0}}^{\rm avg}=\bar{\delta}_{h,x_{0}}-|\Gamma_{h}|^{-1}\int_{\Gamma_{h}}\bar{\delta}_{h,x_{0}}\in\mathring{S}_{h}(\Gamma_{h})$ then

$$\|\partial_{t}{H}_{h}(t,\cdot,x_{0})\|_{L^{2}(\Gamma_{h})} =\|\partial_{t}e^{t\Delta_{\Gamma_{h},h}}\bar{\delta}_{h,x_{0}}^{\rm avg}\|_{L^{2}(\Gamma_{h})}$$

$$=t^{-2}\|(t\Delta_{\Gamma_{h},h})^{2}e^{t\Delta_{\Gamma_{h},h}}\Delta_{\Gamma_{h},h}^{-1}\bar{\delta}_{h,x_{0}}^{\rm avg}\|_{L^{2}(\Gamma_{h})}$$

$$=t^{-2}\|\Delta_{\Gamma_{h},h}^{-1}\bar{\delta}_{h,x_{0}}^{\rm avg}\|_{L^{2}(\Gamma_{h})}$$

$$\leq Ct^{-2}\|\bar{\delta}_{h,x_{0}}^{\rm avg}\|_{L^{1}(\Gamma_{h})}$$

$$\leq Ct^{-2}.$$

From the norm equivalence on $\Gamma$ and $\Gamma_{h}$, it is straightforward to show the constant of Poincaré inequality (with vanishing average) on $\Gamma_{h}$ is bounded from above by the constant of Poincaré inequality (with vanishing average) on $\Gamma$. As a consequence, the smallest eigenvalue of $-\Delta_{\Gamma_{h}}\colon\mathring{S}_{h}(\Gamma_{h})\rightarrow\mathring{S}_{h}(\Gamma_{h})$ has a lower bound, denoted by $\lambda_{0}$, which is uniform w.r.t $h$ and only depends on $\Gamma$. Hence (the average vanishes due to the differential $\partial_{t}$),

$$\|\partial_{t}{H}_{h}(t,\cdot,x_{0})\|_{L^{2}(\Gamma_{h})}\leq Ce^{-\lambda_{0}(t-1)}\|\partial_{t}{H}_{h}(1,\cdot,x_{0})\|_{L^{2}(\Gamma_{h})}\leq Ce^{-\lambda_{0}(t-1)}.$$

Similarly, we also have

$$\|\partial_{tt}{H}_{h}(t,\cdot,x_{0})\|_{L^{2}(\Gamma_{h})}+\|\partial_{t}{H}(t,\cdot,x_{0})\|_{L^{2}(\Gamma_{h})}+\|\partial_{tt}{H}(t,\cdot,x_{0})\|_{L^{2}(\Gamma_{h})}\leq Ce^{-\lambda_{0}(t-1)}\quad t\geq 1.$$

The proof of Lemma 4.1 is complete.

5 Perturbation arguments for an evolving surface
(Proof of Theorem 2.3)

Following usual notational conventions, we will always identify a finite element function and the corresponding vector collecting its nodal values. Such representation is unique if we have specified the underlying domain. For example, given any $t,s\in[0,T]$, the two integrands of

$$\int_{\Gamma_{h}(t)}v_{h}\quad\mbox{and}\quad\int_{\Gamma_{h}(s)}v_{h}$$

have the same vector of nodal values, denoted by ${\bf v}$, but are defined on different domains $\Gamma_{h}(t)$ and $\Gamma_{h}(s)$. When the underlying domain is specified, ${\bf v}$ is automatically substantialized to a finite element function $v_{h}$ on that domain. Since all of the quantitative computations in this paper involve either integrals or norms, our notations for finite element functions will always have a unique and clear meaning. For another example, $\|v_{h}\|_{\Gamma_{h}(t)}$ and $\|v_{h}\|_{\Gamma_{h}(s)}$ denote the norms of a finite element function (a nodal vector) on the two different surfaces $\Gamma_{h}(t)$ and $\Gamma_{h}(s)$, respectively. Correspondingly, the interpolation operator $I_{h}$ should be interpreted as the determination of the nodal vector which uniquely corresponds to a finite element function after specifying the underlying surface.

Given any $s,t\in[0,T]$, we first pull back the scheme onto the piecewise polynomial approximate surface $\Gamma_{h}(s)$ and then carry out temporal perturbation argument on $\Gamma_{h}(s)$.

For the fixed time $s\in[0,T]$, we define the scalar function $a(t,x)$ and the ${\mathbb{R}}^{(m+1)\times(m+1)}$-valued function $B(t,x)$ with $x\in\Gamma(s)$ to be the smooth prefactor of the following change of variables via the smooth flow parametrized on $\Gamma(s)$, i.e. $F_{s}(t)\colon\Gamma(s)\rightarrow\Gamma(t)$,

$$\int_{\Gamma(t)}uv =\int_{\Gamma(s)}a(t,x)(u\circ F_{s})\,(v\circ F_{s}),$$

$$\int_{\Gamma(t)}\nabla_{\Gamma(t)}u\cdot\nabla_{\Gamma(t)}v =\int_{\Gamma(s)}B(t,x)\nabla_{\Gamma(s)}(u\circ F_{s})\cdot\nabla_{\Gamma(s)}(v\circ F_{s}).$$

Similarly, we define the piecewise smooth scalar function $\bar{a}(t,x)$ and matrix-valued function $\bar{B}(t,x)$ with $x\in\Gamma_{h}(s)$ to be the prefactor of the following change of variables via the smooth flow parametrized on $\Gamma_{h}(s)$, i.e. $F_{s}(t)\circ q(s)\colon\Gamma_{h}(s)\rightarrow\Gamma(t)$,

$$\int_{\Gamma(t)}uv =\int_{\Gamma_{h}(s)}\bar{a}(t,x)(u\circ F_{s}\circ q(s))\,(v\circ F_{s}\circ q(s)),$$

$$\int_{\Gamma(t)}\nabla_{\Gamma(t)}u\cdot\nabla_{\Gamma(t)}v =\int_{\Gamma_{h}(s)}\bar{B}(t,x)\nabla_{\Gamma_{h}(s)}(u\circ F_{s}\circ q(s))\cdot\nabla_{\Gamma_{h}(s)}(v\circ F_{s}\circ q(s)),$$

and define the piecewise smooth $\bar{a}_{h}(t,x)$ and $\bar{B}_{h}(t,x)$ with $x\in\Gamma_{h}(s)$ to be the prefactor of the following change of variables via the discrete flow parametrized on $\Gamma_{h}(s)$, i.e. $F_{h,s}(t):=I_{h}(F_{s}(t)\circ q(s))\colon\Gamma_{h}(s)\rightarrow\Gamma_{h}(t)$,

$$\int_{\Gamma_{h}(t)}uv =\int_{\Gamma_{h}(s)}\bar{a}_{h}(t,x)u\circ F_{h,s}\,v\circ F_{h,s},$$

$$\int_{\Gamma_{h}(t)}\nabla_{\Gamma_{h}(t)}u\cdot\nabla_{\Gamma_{h}(t)}v =\int_{\Gamma_{h}(s)}\bar{B}_{h}(t,x)\nabla_{\Gamma_{h}(s)}(u\circ F_{h,s})\cdot\nabla_{\Gamma_{h}(s)}(v\circ F_{h,s}).$$

Since $F_{h,s}(t)=I_{h}(F_{s}(t)\circ q(s))$, by the interpolation error estimates, it follows that for any $s,t\in[0,T]$

$$\|\bar{a}_{h}(t,\cdot)-\bar{a}(t,\cdot)\|_{L^{\infty}(\Gamma_{h}(s))} \leq Ch^{k+1},$$

$$\|\bar{B}_{h}(t,\cdot)-\bar{B}(t,\cdot)\|_{L^{\infty}(\Gamma_{h}(s))} \leq Ch^{k+1},$$

and moreover from the geometric perturbation error estimates (cf. [13])

$$\|\bar{a}^{\ell}(t,\cdot)-a(t,\cdot)\|_{L^{\infty}(\Gamma(s))} \leq Ch^{k+1},$$

$$\|\bar{B}^{l}(t,\cdot)-B(t,\cdot)\|_{L^{\infty}(\Gamma(s))} \leq Ch^{k}.$$

Hence $\tilde{B}(t,x):=(B(t)\circ q(s))(x)\approx\bar{B}_{h}(t,x),x\in\Gamma_{h}(s),s,t\in[0,T]$ is a good approximation which is globally continuous and piecewise smooth. This global continuity will allow us to take advantage of the definition of discrete Laplacian. With the definitions above we have, for any $v_{h}\in S_{h}$,

$$\quad(\partial_{t}u_{h},v_{h})_{\Gamma_{h}(s)}+(\nabla_{\Gamma_{h}(s)}u_{h},\nabla_{\Gamma_{h}(s)}v_{h})_{\Gamma_{h}(s)}$$

$$=(\bar{a}_{h}(s,\cdot)\partial_{t}u_{h}(t,\cdot),v_{h})_{\Gamma_{h}(s)}+(\bar{B}_{h}(s,\cdot)\nabla_{\Gamma_{h}(s)}u_{h}(t,\cdot),\nabla_{\Gamma_{h}(s)}v_{h})_{\Gamma_{h}(s)}$$

$$=(\bar{a}_{h}(t,\cdot)f_{h}(t,\cdot),v_{h})_{\Gamma_{h}(s)}+((\bar{a}_{h}(s,\cdot)-\bar{a}_{h}(t,\cdot))\partial_{t}u_{h}(t,\cdot),v_{h})_{\Gamma_{h}(s)}$$

$$\quad+((\bar{B}_{h}(s,\cdot)-\bar{B}_{h}(t,\cdot))\nabla_{\Gamma_{h}(s)}u_{h}(t,\cdot),\nabla_{\Gamma_{h}(s)}v_{h})_{\Gamma_{h}(s)}$$

$$=:(\bar{a}_{h}(t,\cdot)f_{h}(t,\cdot),v_{h})_{\Gamma_{h}(s)}+((\bar{a}_{h}(s,\cdot)-\bar{a}_{h}(t,\cdot))\partial_{t}u_{h}(t,\cdot),v_{h})_{\Gamma_{h}(s)}$$

$$\quad+(\Delta_{\Gamma_{h}(s),h}^{B_{t}}u_{h}(t,\cdot),v_{h})_{\Gamma_{h}(s)},$$ (5.1)

where $\Delta_{\Gamma_{h}(s),h}^{B_{t}}u_{h}\in S_{h}(\Gamma_{h}(s))$ is defined as the Riesz representative of the following linear functional on $S_{h}(\Gamma_{h}(s))$:

$$((\bar{B}_{h}(s)-\bar{B}_{h}(t))\nabla_{\Gamma_{h}(s)}u_{h}(t),\nabla_{\Gamma_{h}(s)}\cdot)_{\Gamma_{h}(s)}\colon S_{h}(\Gamma_{h}(s))\rightarrow S_{h}(\Gamma_{h}(s)).$$

The estimate of $\Delta_{\Gamma_{h}(s),h}^{B_{t}}u_{h}$ is given in the lemma below.

Lemma 5.1

For any $q\in[1,\infty]$, we have the estimate

$$\|\Delta_{\Gamma_{h}(s),h}^{B_{t}}u_{h}(t)\|_{L^{q}(\Gamma_{h}(s))}\leq C\|\nabla_{\Gamma_{h}(s)}u_{h}(t)\|_{L^{q}(\Gamma_{h}(s))}+C|s-t|\|\Delta_{\Gamma_{h}(s),h}u_{h}(t)\|_{L^{q}(\Gamma_{h}(s))}.$$

Proof. We apply change of variables, Leibniz rule, super-approximation estimate (cf. (P3) in Section 3.2) and the inverse inequality to get

$$\quad(\Delta_{\Gamma_{h}(s),h}^{B_{t}}u_{h},v_{h})_{\Gamma_{h}(s)}$$

$$=(\nabla_{\Gamma_{h}(s)}u_{h},(\tilde{B}(s)-\tilde{B}(t))\nabla_{\Gamma_{h}(s)}v_{h})_{\Gamma_{h}(s)}$$

$$\quad+\bigg{(}\nabla_{\Gamma_{h}(s)}u_{h},\bigg{(}(\bar{B}_{h}(s)-\bar{B}_{h}(t))-(\tilde{B}(s)-\tilde{B}(t))\bigg{)}\nabla_{\Gamma_{h}(s)}v_{h}\bigg{)}_{\Gamma_{h}(s)}$$

$$=-(\nabla_{\Gamma_{h}(s)}u_{h},\nabla_{\Gamma_{h}(s)}(\tilde{B}(s)-\tilde{B}(t))v_{h})_{\Gamma_{h}(s)}+(\nabla_{\Gamma_{h}(s)}u_{h},\nabla_{\Gamma_{h}(s)}I_{h}((\tilde{B}(s)-\tilde{B}(t))v_{h}))_{\Gamma_{h}(s)}$$

$$\quad+(\nabla_{\Gamma_{h}(s)}u_{h},\nabla_{\Gamma_{h}(s)}(1-I_{h})((\tilde{B}(s)-\tilde{B}(t))v_{h}))_{\Gamma_{h}(s)}$$

$$\quad+\bigg{(}\nabla_{\Gamma_{h}(s)}u_{h},\bigg{(}(\bar{B}_{h}(s)-\bar{B}_{h}(t))-(\tilde{B}(s)-\tilde{B}(t))\bigg{)}\nabla_{\Gamma_{h}(s)}v_{h}\bigg{)}_{\Gamma_{h}(s)}$$

$$\leq C\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{p}(\Gamma_{h}(s))}\|v_{h}\|_{L^{q}(\Gamma_{h}(s))}+C|s-t|\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{p}(\Gamma_{h}(s))}\|v_{h}\|_{L^{q}(\Gamma_{h}(s))}$$

$$\quad+C(h+h^{k})\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{p}(\Gamma_{h}(s))}\|\nabla_{\Gamma_{h}(s)}v_{h}\|_{L^{q}(\Gamma_{h}(s))}$$

$$\leq C\big{(}\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{p}(\Gamma_{h}(s))}+|s-t|\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{p}(\Gamma_{h}(s))}\big{)}\|v_{h}\|_{L^{q}(\Gamma_{h}(s))}\quad\forall v_{h}\in S_{h},$$

where in the second to last inequality we have used the consistency estimate $\|\bar{B}_{h}(t)-\tilde{B}(t)\|_{L^{\infty}(\Gamma_{h}(s))}\leq Ch^{k}$ for any $t\in[0,T]$, and all of the constants are independent of $s,t\in[0,T]$ (possibly depend on $T$). Thus we conclude the lemma by duality. $\square$

Remark 5.1

The above proof additionally leads to, for any $s\in[0,T]$,

$$\|\Delta_{\Gamma(s),h}u_{h}\|_{L^{q}(\Gamma(s))}\leq C\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{q}(\Gamma_{h}(s))}+C\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{q}(\Gamma_{h}(s))},$$

and, for any $s^{\prime}\in[0,T]$,

$$\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{q}(\Gamma_{h}(s))}\leq C\|\nabla_{\Gamma_{h}(s^{\prime}),h}u_{h}\|_{L^{q}(\Gamma_{h}(s^{\prime}))}+C\|\Delta_{\Gamma_{h}(s^{\prime}),h}u_{h}\|_{L^{q}(\Gamma_{h}(s^{\prime}))}.$$

Lemma 5.2

The following discrete interpolation inequality holds on the polynomial surface $\Gamma_{h}(s)$, for all $q\in(1,\infty)$,

$$\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{q}(\Gamma_{h}(s))} \leq C\epsilon^{-1}\|u_{h}\|_{L^{q}(\Gamma_{h}(s))}+C\epsilon\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{q}(\Gamma_{h}(s))}.$$

Proof. We define function $u$ to be the solution of the following elliptic equation on $\Gamma(s)$

$$-(\Delta_{\Gamma(s)}-1)u=-(\Delta_{\Gamma(s),h}-1)u_{h}^{\ell},$$

or equivalently $R_{h}u=u_{h}^{\ell}$ where $R_{h}$ is the Ritz projection associated to the finite element space $S_{h}(\Gamma(s))$ the elliptic operator $-\Delta_{\Gamma(s)}+1$ (see Section 3.4). Then the elliptic regularity theory says

$$\|u\|_{W^{2,q}(\Gamma(s))}\leq C\|(\Delta_{\Gamma(s),h}-1)u_{h}^{\ell}\|_{L^{q}(\Gamma(s))}.$$ (5.2)

Moreover from the inverse inequality

$$\|(\Delta_{\Gamma(s),h}-1)u_{h}^{\ell}\|_{L^{q}(\Gamma(s))} =\sup_{v_{h}\in S_{h},\|v_{h}^{\ell}\|_{L^{q^{\prime}}(\Gamma(s))}=1}(\nabla_{\Gamma(s)}u_{h}^{\ell},\nabla_{\Gamma(s)}v_{h}^{\ell})_{\Gamma(s)}+(u_{h}^{\ell},v_{h}^{\ell})_{\Gamma(s)}$$

$$\leq Ch^{-2}\|u_{h}^{\ell}\|_{L^{q}(\Gamma(s))}.$$ (5.3)

Define $A=-\Delta_{\Gamma(s)}+1$ and we get

$$(u,v)_{\Gamma(s)}=(Au,A^{-1}v)_{\Gamma(s)}=((-\Delta_{\Gamma(s),h}+1)u_{h}^{\ell},A^{-1}v)_{\Gamma(s)}$$

$$=((-\Delta_{\Gamma(s),h}+1)u_{h}^{\ell},R_{h}A^{-1}v)_{\Gamma(s)}+((-\Delta_{\Gamma(s),h}+1)u_{h}^{\ell},(1-R_{h})A^{-1}v)_{\Gamma(s)}$$

$$=((\nabla_{\Gamma(s),h}+1)u_{h}^{\ell},(\nabla_{\Gamma(s),h}+1)A^{-1}v)_{\Gamma(s)}+((-\Delta_{\Gamma(s),h}+1)u_{h}^{\ell},(1-R_{h})A^{-1}v)_{\Gamma(s)}$$

$$=(u_{h}^{\ell},v)_{\Gamma}+((-\Delta_{\Gamma(s),h}+1)u_{h}^{\ell},(1-R_{h})A^{-1}v)_{\Gamma(s)},$$

and therefore

$$\|u\|_{L^{q}(\Gamma(s))}=\sup_{\|v\|_{L^{q^{\prime}}(\Gamma(s))}=1}(u,v)_{\Gamma(s)}$$

$$\leq\|u_{h}^{\ell}\|_{L^{q}(\Gamma(s))}+\|(-\Delta_{\Gamma(s),h}+1)u_{h}^{\ell}\|_{L^{q}(\Gamma(s))}\sup_{\|v\|_{L^{q^{\prime}}(\Gamma(s))}=1}\|(1-R_{h})A^{-1}v\|_{L^{q^{\prime}}(\Gamma(s))}$$

$$\leq\|u_{h}^{\ell}\|_{L^{q}(\Gamma(s))}+Ch^{2}\|(-\Delta_{\Gamma(s),h}+1)u_{h}^{\ell}\|_{L^{q}(\Gamma(s))}\sup_{\|v\|_{L^{q^{\prime}}(\Gamma(s))}=1}\|A^{-1}v\|_{W^{2,q^{\prime}}(\Gamma(s))}$$

$$\leq C\|u_{h}^{\ell}\|_{L^{q}(\Gamma(s))},$$ (5.4)

where in the last line we have used (5) and the elliptic regularity theory $\|A^{-1}v\|_{W^{2,q^{\prime}}(\Gamma(s))}\leq C\|v\|_{L^{q^{\prime}}(\Gamma(s))}=C$.

Hence the norm equivalence, Remark 5.1, (5.2), (5) and the $W^{1,q}$-stability of $R_{h}$ (see Section 3.4) imply

$$\|\nabla_{\Gamma(s)}u_{h}^{\ell}\|_{L^{q}(\Gamma(s))} \leq\|\nabla_{\Gamma(s)}(R_{h}u-I_{h}u)\|_{L^{q}(\Gamma(s))}+\|\nabla_{\Gamma(s)}(I_{h}-1)u\|_{L^{q}(\Gamma(s))}+\|\nabla_{\Gamma(s)}u\|_{L^{q}(\Gamma(s))}$$

$$\leq Ch\|u\|_{W^{2,q}(\Gamma(s))}+\|u\|_{L^{q}(\Gamma(s))}^{1/2}\|u\|_{W^{2,q}(\Gamma(s))}^{1/2}$$

$$\leq C\epsilon^{-1}\|u_{h}\|_{L^{q}(\Gamma_{h}(s))}+C\epsilon\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{q}(\Gamma_{h}(s))}+C\epsilon\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{q}(\Gamma_{h}(s))}.$$

We complete the proof by absorbing $\epsilon\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{q}(\Gamma_{h}(s))}$ into the left-hand side. $\square$

Applying the discrete maximal regularity on the stationary surface $\Gamma_{h}(s)$ with $p,q\in(1,\infty)$ (Theorem 2.19) to (5) and using the norm equivalence between $\Gamma_{h}(s)$ and $\Gamma_{h}(s)$, we get

$$\quad\|\partial_{t}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}+\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}$$

$$\leq C\|f_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}+C\|(\bar{a}_{h}(s,\cdot)-\bar{a}_{h}(t,\cdot))\partial_{t}u_{h}(t,\cdot)\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}$$

$$\quad+C\|\Delta_{\Gamma_{h}(s),h}^{B_{t}}u_{h}(t,\cdot)\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}$$

$$\leq C\|f_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}+C|s-t|\|\partial_{t}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}$$

$$\quad+C\|\nabla_{\Gamma_{h}(s)}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))}+C|s-t|\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(s)))},$$

where in the second inequality we have used Lemma 5.1. Then we apply the norm equivalence, Remark 5.1 and Lemma 5.2 to the inequality above to change the domain from $\Gamma_{h}(s)$ to $\Gamma_{h}(0)$:

$$\quad\|\partial_{t}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}+\|\Delta_{\Gamma_{h}(s),h}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}$$

$$\leq C\|f_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}+C\|u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}+C|s-t|\|\partial_{t}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}$$

$$\quad+C|s-t|\|\Delta_{\Gamma_{h}(0),h}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}.$$ (5.5)

We define (cf. the proof of Theorem 3.1 and (4.22) in [23])

$$L_{p,q}(t):=\|\partial_{t}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}^{p}+\|\Delta_{\Gamma_{h}(0),h}u_{h}\|_{L^{p}(0,t;L^{q}(\Gamma_{h}(0)))}^{p}$$

with

$$\partial_{t}L_{p,q}(t)=\|\partial_{t}u_{h}(t,\cdot)\|_{L^{q}(\Gamma_{h}(0)))}^{p}+\|\Delta_{\Gamma_{h}(0),h}u_{h}(t,\cdot)\|_{L^{q}(\Gamma_{h}(0)))}^{p}.$$

Then (5) implies that

$$L_{p,q}(s) =\|\partial_{t}u_{h}\|_{L^{p}(0,s;L^{q}(\Gamma_{h}(0)))}^{p}+\|\Delta_{\Gamma_{h}(0),h}u_{h}\|_{L^{q}(0,s;L^{2}(\Gamma_{h}(0)))}^{p}$$

$$\leq C\|f_{h}\|_{L^{p}(0,s;L^{q}(\Gamma_{h}(0)))}^{p}+C\|u_{h}\|_{L^{p}(0,s;L^{q}(\Gamma_{h}(0)))}^{p}$$

$$\quad+C\int_{0}^{s}(s-t)^{p}\big{(}\|\partial_{t}u_{h}(t,\cdot)\|_{L^{q}(\Gamma_{h}(0))}^{p}+\|\Delta_{\Gamma_{h}(0),h}u_{h}(t,\cdot)\|_{L^{q}(\Gamma_{h}(0))}^{p}\big{)}{\mathrm{d}}t$$

$$\leq C\|f_{h}\|_{L^{p}(0,s;L^{q}(\Gamma_{h}(0)))}^{p}+CT^{p}\int_{0}^{s}L_{p,q}(t){\mathrm{d}}t+CpT^{p-1}\int_{0}^{s}L_{p,q}(t){\mathrm{d}}t,$$

where we have used the following estimates in the last line (cf. [23, p. 10])

$$\|u_{h}\|_{L^{p}(0,s;L^{q}(\Gamma_{h}(0)))}^{p}\leq T^{p}\int_{0}^{s}L_{p,q}(t){\mathrm{d}}t,$$

$$\int_{0}^{s}(s-t)^{p}\big{(}\|\partial_{t}u_{h}(t,\cdot)\|_{L^{q}(\Gamma_{h}(0))}^{p}+\|\Delta_{\Gamma_{h}(0),h}u_{h}(t,\cdot)\|_{L^{q}(\Gamma_{h}(0))}^{p}\big{)}{\mathrm{d}}t\leq pT^{p-1}\int_{0}^{s}L_{p,q}(t){\mathrm{d}}t.$$

Thus, by Grönwall’s inequality, it holds that

$$L_{p,q}(s)\leq C\|f_{h}\|_{L^{p}(0,s;L^{q}(\Gamma_{h}(0)))}^{p}\quad\forall s\in[0,T].$$

Therefore, the proof of Theorem 2.3 is complete.

Appendix A Proof of Lemma 3.5

The following lemma is proved in [32, Lemma A.1], which can be easily extended to our current scenario.

Lemma A.1

Suppose that, for any $t\in(0,T)$, $\phi_{h}^{\ell}(t)\in S_{h}(\Gamma)$ satisfies

$$(\partial_{t}\phi_{h}^{\ell},\chi_{h}^{\ell})_{\Gamma}+(\nabla_{\Gamma}\phi_{h}^{\ell},\nabla_{\Gamma}\chi_{h}^{\ell})_{\Gamma}=0\quad\forall\chi_{h}^{\ell}\in S_{h}^{0}(\Gamma_{j}^{\prime\prime}),t\in(0,d_{j}^{2}),$$

$$(\partial_{t}\phi_{h}^{\ell},\chi_{h}^{\ell})_{\Gamma}+(\nabla_{\Gamma}\phi_{h}^{\ell},\nabla_{\Gamma}\chi_{h}^{\ell})_{\Gamma}=0\quad\forall\chi_{h}^{\ell}\in S_{h}^{0}(D_{j}^{\prime\prime}),t\in(d_{j}^{2}/4,2d_{j}^{2}).$$

Then we have

$${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}$$

$$\leq(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)\big{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\nabla_{\Gamma}\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}\big{)}$$

$$\quad+C\epsilon^{-1}\bigg{(}d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}\|\phi_{h}^{\ell}(0)\|_{\Gamma_{j}^{\prime}}+\|\nabla_{\Gamma}\phi_{h}^{\ell}(0)\|_{1,\Gamma_{j}^{\prime}}\bigg{)},$$

where the constant $C>0$ is independent of $h$, $j$ and $C_{*}$.

It is easy the construct the space-time cut-off function $\tilde{\omega}$ which is constant one on $Q_{j}^{\prime}$ and zero outside $Q_{j}^{\prime\prime}$. We then define $\tilde{\phi}=\tilde{\omega}\phi$ and $\tilde{\phi}_{h}^{\ell}\in S_{h}(\Gamma)$ to be the finite element solution of

$$(\partial_{t}\tilde{\phi}-a_{h}\partial_{t}\tilde{\phi}_{h}^{\ell},\chi_{h}^{\ell})_{\Gamma}+(\nabla_{\Gamma}\tilde{\phi}-B_{h}\nabla_{\Gamma}\tilde{\phi}_{h}^{\ell},\chi_{h}^{\ell})_{\Gamma}=0\quad\forall\chi_{h}^{\ell}\in S_{h}(\Gamma),$$ (A.1)

with the initial condition $\tilde{\phi}_{h}^{\ell}(0)=\tilde{\phi}(0)=0$. Obviously, $\tilde{\phi}_{h}^{\ell}$ is zero outside $Q_{j}^{\prime\prime}$ as well.

Then,

$$\quad(a_{h}\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell}),\chi_{h}^{\ell})_{\Gamma}+(B_{h}\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell}),\chi_{h}^{\ell})_{\Gamma}$$

$$=(a_{h}\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}),\chi_{h}^{\ell})_{\Gamma}+(B_{h}\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}),\chi_{h}^{\ell})_{\Gamma}$$

$$\quad+((a_{h}-1)\partial_{t}\tilde{\phi},\chi_{h}^{\ell})_{\Gamma}+((B_{h}-1)\nabla_{\Gamma}\tilde{\phi},\chi_{h}^{\ell})_{\Gamma}\quad\forall\chi_{h}^{\ell}\in S_{h}(\Gamma).$$

Testing the above equation with $\chi_{h}^{\ell}=I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell}\in S_{h}(\Gamma)$ and using Lemma 3.1, temporal integration by parts and Hölder’s inequality (cf. [32, p. 37]), we get

$$\quad\|I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell}\|_{L^{\infty}L^{2}(\mathcal{Q})}^{2}+\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}^{2}$$

$$\leq\|I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell}\|_{L^{2}L^{2}(\mathcal{Q})}\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}+\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Ch^{k+1}\|\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}+Ch^{k+1}\|\nabla_{\Gamma}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})},$$ (A.2)

and, similarly, testing the above equation with $\chi_{h}^{\ell}=\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\in S_{h}(\Gamma)$ leads to

$$\quad\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}^{2}+\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{\infty}L^{2}(\mathcal{Q})}^{2}$$

$$\leq\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}+\|\nabla_{\Gamma}\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Ch^{k+1}\|\partial_{t}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}+Ch^{k+1}\|\nabla_{\Gamma}\partial_{t}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}.$$ (A.3)

Applying Young’s inequality to (A) and (A), we obtain

$$\quad\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}+\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{\infty}L^{2}(\mathcal{Q})}$$

$$\quad+d_{j}^{-1}\|I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell}\|_{L^{\infty}L^{2}(\mathcal{Q})}+d_{j}^{-1}\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\leq C\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}\|\nabla_{\Gamma}\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Cd_{j}^{-2}\|I_{h}\tilde{\phi}-\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}^{-1}\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Ch^{k+1}\|\partial_{t}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}h^{k+1}\|\nabla_{\Gamma}\partial_{t}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Cd_{j}^{-2}h^{k+1}\|\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}^{-1}h^{k+1}\|\nabla_{\Gamma}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}.$$

Using the triangle inequality and Hölder’s inequality in the temporal direction,

$$\quad\|\partial_{t}(\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}+d_{j}^{-1}\|\nabla_{\Gamma}(\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+d_{j}^{-2}\|\tilde{\phi}-\tilde{\phi}_{h}^{\ell}\|_{L^{2}L^{2}(\mathcal{Q})}+d_{j}^{-1}\|\nabla_{\Gamma}(\tilde{\phi}-\tilde{\phi}_{h}^{\ell})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\leq C\|\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}\|\nabla_{\Gamma}\partial_{t}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Cd_{j}^{-2}\|I_{h}\tilde{\phi}-\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}^{-1}\|\nabla_{\Gamma}(I_{h}\tilde{\phi}-\tilde{\phi})\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Ch^{k+1}\|\partial_{t}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}h^{k+1}\|\nabla_{\Gamma}\partial_{t}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}$$

$$\quad+Cd_{j}^{-2}h^{k+1}\|\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}+Cd_{j}^{-1}h^{k+1}\|\nabla_{\Gamma}\tilde{\phi}\|_{L^{2}L^{2}(\mathcal{Q})}.$$ (A.4)

By subtracting (3.52) from (A.1) and using suitably arranged cut-off function, we know that for $\chi_{h}^{\ell}\in S_{h}^{0}(\Gamma_{j}^{\prime\prime}),t\in(0,d_{j}^{2})$ and $\chi_{h}^{\ell}\in S_{h}^{0}(D_{j}^{\prime\prime}),t\in(d_{j}^{2}/4,2d_{j}^{2})$

$$(a_{h}\partial_{t}(\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell}),\chi_{h}^{\ell})_{\Gamma}+(B_{h}\nabla_{\Gamma}(\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell}),\nabla_{\Gamma}\chi_{h}^{\ell})_{\Gamma}=0.$$

Then we apply Lemma A.1 to $\tilde{\phi}_{h}-\phi_{h}$ and get

$${\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}}$$

$$\leq(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)\bigg{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}\bigg{)}$$

$$\quad+C\epsilon^{-1}\bigg{(}d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}\|\phi_{h}^{\ell}(0)\|_{Q_{j}^{\prime}}+\|\phi_{h}^{\ell}(0)\|_{1,Q_{j}^{\prime}}\bigg{)}$$

$$\leq(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)\bigg{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(\tilde{\phi}-\phi_{h}^{\ell})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\tilde{\phi}-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}\bigg{)}$$

$$\quad+C\epsilon^{-1}\bigg{(}d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\tilde{\phi}-\phi_{h}^{\ell}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}\|\phi_{h}^{\ell}(0)\|_{Q_{j}^{\prime}}+\|\phi_{h}^{\ell}(0)\|_{1,Q_{j}^{\prime}}\bigg{)}$$

$$\quad+(Ch^{1/2}d_{j}^{-1/2}+C\epsilon^{-1}hd_{j}^{-1}+\epsilon)\bigg{(}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\partial_{t}(\tilde{\phi}_{h}^{\ell}-\tilde{\phi})\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}}+d_{j}^{-1}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\tilde{\phi}_{h}^{\ell}-\tilde{\phi}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{1,Q_{j}^{\prime}}\bigg{)}$$

$$\quad+C\epsilon^{-1}d_{j}^{-2}{\left|\kern-1.07639pt\left|\kern-1.07639pt\left|\tilde{\phi}_{h}^{\ell}-\tilde{\phi}\right|\kern-1.07639pt\right|\kern-1.07639pt\right|}_{Q_{j}^{\prime}},$$ (A.5)

where in the last inequality, we use the splitting $\tilde{\phi}-\phi_{h}^{\ell}=(\tilde{\phi}-\tilde{\phi}_{h}^{\ell})+(\tilde{\phi}_{h}^{\ell}-\phi_{h}^{\ell})$.

Finally, we complete the proof of Lemma 3.5 by applying the triangle inequality to (A) and (A) and using the fact that both $\tilde{\phi}$ and $\tilde{\phi}_{h}^{\ell}$ are zero outside $Q_{j}^{\prime\prime}$.