A unified theory for continuous in time evolving finite element space approximations to partial differential equations in evolving domains
Charles M. Elliott, Thomas Ranner

TL;DR
This paper develops a comprehensive theoretical framework for finite element methods applied to partial differential equations on evolving domains, including surfaces and bulk systems, with proven optimal bounds and numerical validation.
Contribution
It introduces a unified abstract variational theory for continuous in time finite element discretizations on evolving domains, including curved spaces and coupled surface-bulk systems.
Findings
Proved optimal a priori error bounds for the discretizations.
Developed evolving finite element spaces on flat and curved geometries.
Numerical experiments confirm theoretical convergence rates.
Abstract
We develop a unified theory for continuous in time finite element discretisations of partial differential equations posed in evolving domains including the consideration of equations posed on evolving surfaces and bulk domains as well coupled surface bulk systems. We use an abstract variational setting with time dependent function spaces and abstract time dependent finite element spaces. Optimal a priori bounds are shown under usual assumptions on perturbations of bilinear forms and approximation properties of the abstract finite element spaces. The abstract theory is applied to evolving finite elements in both flat and curved spaces. Evolving bulk and surface isoparametric finite element spaces defined on evolving triangulations are defined and developed. These spaces are used to define approximations to parabolic equations in general domains for which the abstract theory is shown to…
| error | (eoc) | ||
|---|---|---|---|
| — | |||
| error | (eoc) | ||
|---|---|---|---|
| — | |||
| error | (eoc) | ||
|---|---|---|---|
| — | |||
| error | (eoc) | ||
|---|---|---|---|
| — | |||
| error | (eoc) | error | (eoc) | ||
|---|---|---|---|---|---|
| — | — | ||||
| error | (eoc) | error | (eoc) | ||
|---|---|---|---|---|---|
| — | — | ||||
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
A unified theory for continuous in time evolving finite element space
approximations to partial differential equations in evolving domains
C. M. Elliott1
and
T. Ranner2
Abstract.
We develop a unified theory for continuous in time finite element discretisations of partial differential equations posed in evolving domains including the consideration of equations posed on evolving surfaces and bulk domains as well as coupled surface bulk systems. We use an abstract variational setting with time dependent function spaces and abstract time dependent finite element spaces. Optimal a priori bounds are shown under usual assumptions on perturbations of bilinear forms and approximation properties of the abstract finite element spaces. The abstract theory is applied to evolving finite elements in both flat and curved spaces. Evolving bulk and surface isoparametric finite element spaces defined on evolving triangulations are defined and developed. These spaces are used to define approximations to parabolic equations in general domains for which the abstract theory is shown to apply. Numerical experiments are described which confirm the rates of convergence.
Key words and phrases:
evolving finite element spaces, abstract error analysis, parabolic equations on moving domains, advection-diffusion on evolving surfaces, bulk-surface parabolic equations, evolving surface finite element methods, evolving bulk finite element methods
1Mathematics Institute, Zeeman Building, University of Warwick, Coventry. CV4 7AL. UK. Email. [email protected]
The work of CME was supported by a Royal Society Wolfson Research Merit Award.
2School of Computing, University of Leeds, Leeds. LS2 9JT. UK. Email. [email protected]
The work of TR was supported by Engineering and Physical Sciences Research Council (EPSRC EP/J004057/1) and a Leverhulme Trust Early Career Fellowship.
Dedicated to the memory of John W. Barrett
Contents
-
2.2 Abstract formulation of the partial differential equation
-
3.1.3 Abstract discrete variational problem and well posedness
-
6.2 Triangulated hypersurface and surface finite element spaces
-
8 Application I: Parabolic equation on an evolving bulk domain
1. Introduction
In this paper, we develop a unified theory for finite element discretisations of partial differential equations posed in evolving domains including the consideration of equations posed on evolving surfaces and bulk domains as well as coupled surface bulk systems. The discretisation is based on evolving finite element spaces defined on evolving triangulations using isoparametric elements. Optimal order a priori error bounds are proven. This unification is achieved by using an abstract variational setting with time dependent abstract function spaces and time dependent abstract finite element spaces. Given a time dependent Hilbert triple
[TABLE]
the abstract strong formulation is: Find such that
[TABLE]
where is an appropriate time dependent Hilbert space with dual , is an appropriate abstract material derivative arising from the evolution of the spaces, is an abstract elliptic operator satisfying suitable coercivity properties and is a lower order term arising from evolution of the space. Similar to the case of time independent function spaces this equation may be written in variational form as
Problem 1.1**.**
Given , find with such that for almost every ,
[TABLE]
subject to the initial condition .
Here and are generalisations of the Bochner spaces and , where is a time-independent Hilbert triple. The bilinear form is associated with the elliptic operator and the bilinear forms and are associated with the -inner product and its time derivative (i.e. the operator ). The above problem is also well posed for where we seek . Here we seek error bounds for sufficiently smooth solutions.
We formulate and analyse an abstract finite element discretisation based on a Galerkin ansatz with perturbations of the bilinear forms. Under assumptions on the approximation of geometry and the approximation of function spaces by abstract finite element spaces optimal order error bounds are proved. Then we construct realisations of this abstract setting in the context of partial differential equations on evolving domains. This is achieved by means of the construction of evolving finite element spaces on evolving triangulations of time dependent surfaces and bulk domains from first principals. We give a concrete realisation of these spaces based on evolving Lagrange isoparametric finite elements.
This approach is applied to three model problems: a linear parabolic problem in an evolving, bounded bulk domain in , a linear parabolic problem on an evolving compact -dimensional surface embedded in and a linear parabolic problem coupling problems in an evolving, bounded bulk domain in to a problem on its boundary. In each case, we assume that the evolution of the problem domain is prescribed. The abstract approach is applicable to other situations including nonlinear equations, coupled equations and problems with dynamic boundary conditions.
1.1. Background
Partial differential equations posed on complex evolving domains arise in numerous settings such as surfactant transport on fluid interfaces, receptor ligand dynamics on cell surfaces and phase separation on dissolving alloy surfaces (Deckelnick et al., 2015; Elliott et al., 2017; Barrett et al., 2015; Alphonse et al., 2018; Torres-Sánchez et al., 2019; Zimmermann et al., 2019). Numerical approaches to solve these problems include surface finite elements, implicit surface formulations, diffuse interface approximations, trace finite elements, unfitted finite elements, finite volume schemes and mesh free methods. See the works of Dziuk (1988); Dziuk and Elliott (2007); Deckelnick et al. (2009); Dziuk and Elliott (2010); Deckelnick et al. (2014); Deckelnick and Styles (2018); Olshanskii and Reusken (2017); Burman et al. (2016); Lehrenfeld et al. (2018); Lehrenfeld and Olshanskii (2019); Giesselmann and Müller (2014); Deckelnick et al. (2019); Suchde and Kuhnert (2019) and the review of Dziuk and Elliott (2013a).
The problem of solving parabolic problems in evolving bulk domains has been studied for many years. In particular, we mention the ALE (arbitrary Lagrangian-Eulerian) approach first proposed by Hirt et al. (1974) in the context of finite difference methods and by Donea et al. (1982); Hughes et al. (1981) for finite element methods. Analysis of a similar problem considering both spatial and temporal discretisation is given by Badia and Codina (2006); Boffi and Gastaldi (2004); Formaggia and Nobile (1999, 2004); Gastaldi (2001); Nobile (2001); Gawlik and Lew (2015). The analysis by Bonito et al. (2013b, a) provides optimal order convergence for a discrete Galerkin in time approach.
The study of finite element methods for partial differential equations posed on surfaces started with the seminal work of Dziuk (1988). Evolving surface finite elements were introduced and analysed for an advection diffusion equation posed on evolving surfaces by Dziuk and Elliott (2007, 2012, 2013a). In these works optimal error bounds were proved for piecewise linear finite elements on triangulated surface evolved using the normal and advective velocity. In this work we consider a more general parabolic equation on surfaces and discretisations which cover the case of higher-order schemes (Heine, 2005; Demlow, 2009; Kovács, 2018) and arbitrary Lagrangian-Eulerian methods (Elliott and Styles, 2012; Elliott and Venkataraman, 2015). The discretisations presented in this paper can be combined with different time stepping schemes (Dziuk et al., 2011; Dziuk and Elliott, 2012; Lubich et al., 2013; Kovács and Guerra, 2018) to provide a fully discrete scheme. We also mention the analysis of Kovács et al. (2017) who study a diffusion equation on the surface which drives the a priori unknown evolution of the surface.
Our abstract approach will be applied to equations posed on evolving bulk domains and coupled bulk surface systems. See the work of Elliott and Ranner (2013); Gross et al. (2015); Burman et al. (2016) for approaches to stationary surface problems. Kovacs and Lubich (2017) extended the results for piecewise linear elements to a coupled bulk-surface system. The functional analytic setting will be the product of spaces over the bulk domain, , and the surface, .
1.2. Some partial differential equations on evolving domains
Let . For let denote an -dimensional bounded, connected, open subset (domain) in , for . We denote by the boundary assumed to be a sufficiently smooth compact -dimensional orientable hypersurface with unit normal pointing outward from . We write , .
We assume that there exists a sufficiently smooth mapping (called the flow map) satisfying
- (1)
is a diffeomorphism of onto for each ; 2. (2)
.
We will write for the inverse of . We define the material velocity by
[TABLE]
which may be written as
[TABLE]
We use the terminology that is an evolving bulk domain and is an evolving surface domain. In order to define the evolution of the domains as sets we need only to specify the normal velocity of . The tangential components of allows for an arbitrary parametrisation of the domain.
Illustrative examples of the boundary value problems we wish to consider are:
- (1)
Find a time-dependent scalar field such that on an evolving Cartesian bulk domain
[TABLE]
where is a smooth diffusion tensor , a smooth vector field and a smooth scalar field. 2. (2)
Seek a time-dependent scalar surface field such that
[TABLE]
where is a smooth diffusion tensor which maps the tangent space of into itself and is a smooth scalar field. Here is a tangential vector field. We use the notation for the normal time derivative (Cermelli et al., 2005; Dziuk and Elliott, 2013b), which denotes the time derivative of a function along a trajectory on moving in the normal direction. This equation is supposed to hold pointwise on on trajectories evolving from with velocity . 3. (3)
Find a time-dependent pair with a scalar volumetric field and a scalar surface field such that
[TABLE]
where and are positive constants. Equation 1.6b couples the equations posed on the domain and its boundary .
Remark 1.2*.*
- (a)
We do not explain the models in the above equations other than to comment that examples may be derived as conservation laws for scalar quantities subject to diffusive and advective fluxes with linear reactions. Note that in these equations the only part of the material velocity that appears is the normal velocity of . 2. (b)
Our formulation of these initial boundary value problems allows, by means of a tangential velocity field, for a reparametrisation of the evolving domains. In our definition of material velocity the velocity component tangential to is used to define a parametrisation of the domains and . This is an Arbitrary Lagrangian Eulerian (ALE) approach. A consequence is that the differential operators and the function spaces in the abstract setting depend on this ALE velocity. In the discretisations the meshes we use are transported using this velocity. 3. (c)
Working with evolving triangulations on evolving domains has advantages when the domain is unknown and is to be found. For example the hypersurface may be a free boundary or a surface evolving via a geometric evolution equation coupled to the parabolic equation on the bulk or surface domain, Pozzi and Stinner (2017); Kovács et al. (2017); Barrett et al. (2017); Pozzi and Stinner (2018). Evolving triangulated surfaces are often computed as approximations of geometric evolution flows, Deckelnick et al. (2005); Barrett et al. (2019); Kovács et al. (2019). 4. (d)
The analysis applies to the use of evolving finite element spaces in the case of domains which are time independent. This may be useful when devising schemes which evolve the mesh in time in order to adjust to the solution.
1.3. Contributions
The primary contributions of this work which generalises the evolving surface finite element continuous in time analysis of Dziuk and Elliott (2007, 2013a) are:-
- •
The analysis is carried out in a generalised abstract setting extending earlier work by allowing for arbitrary parametrisations and higher order approximations together with domains with boundaries.
- •
We construct evolving finite element spaces and derive approximation results for quasi-uniform evolving surface and bulk triangulations using isoparametric elements.
- •
We show that the resulting finite element methods satisfy the necessary geometric and interpolation estimates for the abstract theory to apply.
- •
We give a complete presentation of the theory for the numerical approximation of parabolic equations on evolving domains using evolving finite element spaces and give three examples of concrete realisations to initial boundary value problems involving time dependent domains in flat and curved space. In particular the approach to evolving bulk domains with boundary requires a more complex treatment than that of surfaces without boundary and the abstract theory is structured to allow this.
- •
The analysis provides error bounds for higher order isoparametric approximations of the case of evolving curved hypersurfaces, evolving bulk domains and evolving coupled surface- bulk systems.
1.4. Extensions
- (a)
All the theory presented in this paper will be applicable with the addition of sufficiently smooth right hand side functionals under natural appropriate assumptions. 2. (b)
Although we will apply our abstract formulation to settings where we approximate a known continuous domain, our intention is that this framework may also be applied to situations where the evolution of the domain is a priori unknown. One way to interpret this work would be to determine the minimum requirements of an evolving finite element method in order to have a sensible, well-posed method. 3. (c)
In this paper we do not consider the example of a parabolic PDE posed on an evolving -dimensional curved sub-manifold in with a moving boundary on which there is an appropriate boundary condition. On the other hand, it is possible to take the perspective in our first example on an evolving bulk domain that is an evolving flat sub-manifold with moving boundary and work with the definitions from Sec. 6 (see also Ex. 6.7(c)), though we choose not to follow this approach.
1.5. Outline
The article is set out in three parts. The first part comprises the abstract theory. In Sec. 2, we introduce the abstract functional analytic setting in which we pose the continuous partial differential equations. An abstract analysis of evolving finite element methods is provided in Sec. 3. The second part comprises the construction of and approximation theory for evolving finite element spaces. We introduce our basic theory for evolving bulk and evolving surface finite elements in Sec. 4 and Sec. 6. Sec. 5 and 7 give technical details on relating the finite elements used in our computational methods to the structures in the underlying continuous problem. The third part comprises applications to three PDE settings. Sec. 8 to Sec. 10 apply these ideas to tackle three model problems. In Sec. 11 results of numerical experiments are given confirming the proven error bounds. Finally, there is an appendix (App. A) with a technical result concerning a Faá di Bruno formula for parametric surfaces.
Part I Abstract theory
2. Abstract formulation
2.1. Evolving function spaces
We introduce an abstract functional analytic setting formulated by Alphonse et al. (2015a) generalising the surface parabolic PDE setting of Vierling (2014). One of the key novelties of these works is to provide the basic theory for evolving Bochner-like spaces for evolving Hilbert spaces such as in order make a definition similar to “”. Using this formulation, we can pose partial differential equations on evolving domains in a fully rigorous setting.
Remark 2.1*.*
The work of Alphonse et al. (2015a) uses a Lagrangian formulation where the evolving domain is parametrised over the initial domain. This matches well with the arbitrary Lagrangian-Eulerian finite element methods we will consider. The setting may be applied to evolving parametrisations of fixed domains. A different functional analytic setting may be more appropriate for different discretisation approaches such as the trace finite element method (Olshanskii et al., 2014; Olshanskii and Reusken, 2017) or the implicit surface approach (Dziuk and Elliott, 2010).
Definition 2.2* (Compatibility).*
For , let be a separable Hilbert space and denote by . Let be a family of invertible, linear homeomorphisms, with inverse , such that there exists such that for every
[TABLE]
and such that the map is continuous for all . Under these circumstances, we call the pair compatible. We call the map the push-forward operator and the pull-back operator.
Remark 2.3*.*
If is a closed subspace in for each and maps , then form a compatible pair.
Definition 2.4* (Evolving Hilbert triple).*
For each , let and be real, separable Hilbert spaces with and such that inclusion is continuous and dense. We will write and for the norms on and , for the inner product on and for the pairing of with its dual. We assume there exists a family of linear homeomorphisms such that and are compatible. We will write for . It follows that continuously and densely. Under these assumptions, we say that is an evolving Hilbert triple.
For a compatible pair, we can define an equivalent structure to Bochner spaces in an evolving context. For a compatible pair, we define to be
[TABLE]
with norm
[TABLE]
One can show that the space is a separable Hilbert space (Alphonse et al., 2015a, Cor. 2.12), is isomorphic to and
[TABLE]
For , we also define the space of smoothly evolving in time functions by
[TABLE]
where is the space of -valued infinitely differentiable functions compactly supported in the interval .
For , we can define a strong material derivative which we denote by by
[TABLE]
This is a temporal derivative which takes into account that fact that is changing as well as the function .
Remark 2.5*.*
If and are both compatible pairs, with , then the spaces and induced using each push forward map are the same vector space with distinct but equivalent norms. On the other hand, in general, the spaces and , when they are non-trivial, are different spaces for . In such cases where ambiguity will occur we will add the push-forward map to the subscript of so no confusion will occur.
We can define a weak material derivative for which an integration by parts in time formula holds. This takes into account the evolution of the space and is often called a transport formula when applied to the derivative of the time dependent inner product. It generalises the notion of the Reynold’s transport formula, Cermelli et al. (2005). In order to provide this definition, we require a further assumption on .
Assumption 2.6*.*
We shall assume the following for all :
[TABLE]
with the constant independent of .
We define by
[TABLE]
Then we have a bilinear form by
[TABLE]
It can be shown that the map is measurable for and we have the following bound independently of :
[TABLE]
We say a function has a weak material derivative if
[TABLE]
for all .
Remark 2.7*.*
The weak and strong material derivatives coincide for .
Lemma 2.8** (Abstract transport formula).**
For all with weak material derivatives we have
[TABLE]
Proof.
See Alphonse et al. (2015a, Thm. 2.40). ∎
The natural solution space is the space given by
[TABLE]
which we equip with the inner product
[TABLE]
Assumption 2.9*.*
Let be the space given by
[TABLE]
equipped with the inner product
[TABLE]
We assume there is an evolving space equivalence between and that is
[TABLE]
and there exists such that
[TABLE]
2.2. Abstract formulation of the partial differential equation
Let . We assume we are in the setting that we have an evolving Hilbert triple and Ass. 2.9 and 2.6 hold so that we have a weak material derivative, which we denote, for appropriate , and a transport formula for the -inner product.
We assume that we have three time dependent bilinear forms , and
[TABLE]
We consider problems of the following form:
Problem 2.10**.**
Given , find such that for almost every ,
[TABLE]
subject to the initial condition .
Remark 2.11*.*
The abstract formulation of Alphonse et al. (2015a) allows for a weaker formulation with solutions in 2.8 and initial condition in . However in order to obtain optimal order error bounds we require more smoothness for the solution. We leave the consideration of error analysis for solutions which are less regular in time to a future work.
In order to make sense of this formulation we restrict to the following assumptions on the bilinear forms holding for each .
**Assumptions on : **
First, we assume that is symmetric:
[TABLE]
We assume that there exists such that for all , we have
[TABLE]
**Assumptions on : **
We assume the bilinear form satisfies
[TABLE]
and such that there exists such that for all
[TABLE]
**Assumptions on : **
- **(a): **
We assume that the map
[TABLE]
is measurable and can be decomposed
[TABLE]
where and are both measurable bilinear forms:
[TABLE]
We assume that is symmetric and allow . 2. **(b): **
We assume there exists such that
[TABLE]
We note that together As4 and An4 imply
[TABLE] 3. **(c): **
We assume the existence of bilinear forms and such that
[TABLE]
and that there exists such that
[TABLE] 4. **(d): **
We define to be
[TABLE]
and note that Bs2 and Bn2 together imply that
[TABLE]
Remark 2.12*.*
We allow for the case that is non-symmetric. This is similar to the choice of Alphonse et al. (2015a) but in contrast to much of the finite element literature (e.g. Elliott and Venkataraman (2015) use different bilinear forms for diffusion and advection terms in a parabolic operator).
Assumption 2.13*.*
We assume there exists a basis of and a sequence with for each such that there exists (which do not depend on or ) with
[TABLE]
Remark 2.14*.*
In the case of being compactly embedded in , Ass. 2.13 follows by Hilbert-Schmidt theory
Theorem 2.15**.**
Let Assumptions M1, M2, G1, G2, A1, A2, A3, As4, Bs1, Bn1, Bs2, Bn2 and Ass. 2.13 hold. Then problem 2.9 has a unique solution which satisfies the stability bound
[TABLE]
and if then
[TABLE]
Proof.
Employing Ass. 2.13 existence may be proved by a Galerkin argument similar to Alphonse et al. (2015a, Thm. 3.6 and 3.13). The a priori estimates can be shown for Galerkin approximations and then used in standard compactness arguments. We do not show the details here but indicate the arguments in deriving estimates in the continuous setting.
Testing 2.9 with :
[TABLE]
G1 gives after integrating in time
[TABLE]
[TABLE]
We infer that
[TABLE]
and applying a Grönwall inequality we see the first stability bound.
The bound 2.11 uses the decomposition A2. We test 2.9 with
[TABLE]
[TABLE]
so we have
[TABLE]
Integrating forwards in time gives
[TABLE]
Applying A3, An4, and a Young’s inequality we observe:
[TABLE]
[TABLE]
From the previous two inequalities with M2, G2, Bs2, Bn2 and An4, we infer
[TABLE]
Rearranging and applying a Grönwall inequality gives the desired result. ∎
Remark 2.16*.*
Note that also satisfies the variational form of 2.9
[TABLE]
3. Abstract discretisation analysis
In this section we present an abstract discretisation and a numerical analysis. We assume that for each there are Hilbert spaces and for all , for some value of fixed throughout this section. These spaces are used to help with the stability and error analysis. We assume that all constants are independent of unless indicated. Throughout this section the integer will denote a further discretisation parameter denoting the order of approximation. Our method will be based in a finite dimensional subspace .
- •
We assume that the evolving Hilbert space is continuously embedded in uniformly for : For the two norms and there exists a constant such that for all and all , we have
[TABLE]
- •
We assume we have a push forward map such that and () are compatible pairs (Def. 2.2) uniformly in : That is there exists such that, for all ,
[TABLE]
- •
This allows us to define the spaces and as in 2.1 and 2.2. For , we denote by the (strong) material derivative (c.f. 2.4) with respect to the push-forward map defined by
[TABLE]
3.1. Abstract discrete method
Let and . Let be an evolving, finite-dimensional space which is a subspace of at each and satisfies (where ). Since is a closed subspace of it is a Hilbert space and forms a compatible pair (Rem. 2.3). In particular, we have well defined spaces 2.1 and 2.2 and the material derivative is well defined for (2.4 and 3.3).
3.1.1. Basis functions
Let the dimension of be for all . We write for a basis of and push-forward to construct a time dependent basis of by
[TABLE]
The following important transport properties of the basis functions hold.
Lemma 3.1**.**
The material derivative of a basis function is zero,
[TABLE]
Furthermore, any function , which can be written as , satisfies
[TABLE]
Proof.
By definition, it follows that
[TABLE]
so that for a decomposition
[TABLE]
we compute that
[TABLE]
3.1.2. Discrete bilinear forms
Let and be two time dependent bilinear forms:
[TABLE]
We make assumptions similar to those in Sec. 2.2.
We assume that is symmetric:
[TABLE]
We assume that there exists such that for all , we have
[TABLE]
We assume we have a transport formula for the bilinear form : there exists a bilinear form such that
[TABLE]
We assume that there exists such that for all
[TABLE]
We assume that the map
[TABLE]
is measurable, and there exists constants such that for all , we have
[TABLE]
We assume a transport formula for the bilinear form, that there exists a bilinear form such that
[TABLE]
and that there exists such that for all we have
[TABLE]
3.1.3. Abstract discrete variational problem and well posedness
Motivated by the variational form 2.12, we consider semi-discrete problems of the following form:
Problem 3.2**.**
Given , find such that and
[TABLE]
We note that due to Assumption Gh1, the discrete scheme 3.7 can be re-written as
[TABLE]
The solution may be written as a decomposition into the time dependent basis functions of ,
[TABLE]
where . Using this notation Prob. 3.2 is equivalent to finding a solution of the (finite dimensional) system of ordinary differential equations:
[TABLE]
where
[TABLE]
Here, we have used the fact that for .
We first wish to show that there exists a solution to our discrete scheme satisfying a stability bound similar to 2.10 for the continuous case. Due to our abstract formulation the calculations follow in a similar way to Thm. 2.15.
Theorem 3.3** (Existence and stability of finite element method).**
Let Assumptions Mh1, Mh2, Gh1, Gh2, Ah1, Ah2, Ah3, Bh1 and Bh2 hold with constants independent of . Then 3.7 has a unique solution with and there exists a constant independent of such that
[TABLE]
Proof.
We consider the problem in the matrix form 3.10. Since Gh1 and is invertible Mh2, this is equivalent to
[TABLE]
This is a linear system of ordinary equations with coefficients (easily verified). Standard theory implies there exists a unique solution , which can be translated as .
To show the energy bound, we start by testing 3.7 with :
[TABLE]
The transport equality Gh1 implies that
[TABLE]
The desired stability bound follows using the same calculations as used in deriving 2.10.
∎
3.2. Abstract lifted finite element spaces
The discrete space is not assumed to be contained in the continuous space . This is an example of a “variational crime” (Strang and Fix, 2008). However it is convenient to prove error bounds in the spaces and . To do this we use lifted discrete spaces using an embedding map . The error analysis will relate the solution of 2.9 with a so-called lift of the discrete solution.
We require two further time dependent Hilbert spaces and which satisfy for each with the inclusions uniformly continuous. Furthermore, we assume that and are both compatible pairs Def. 2.2 so that we may define the spaces 2.1 and 2.2. These spaces are abstractions of spaces of smoother functions.
The link between all the function spaces used in the section is shown in Fig. 3.1.
3.2.1. Lifting operator
Let . We assume there is a continuous, bijective, linear function with inverse denoted by such that is also a bijection onto . We will denote and , i.e. for and for .
We assume that the lifting map is bounded and that there exists such that for all
[TABLE]
Remark 3.4*.*
As a convention, variables denoted by correspond to lifted objects and to inverse lifted objects. The construction of each depends on the parameter .
3.2.2. Lifted push-forward and pull back maps
The lift of the discrete space push-forward map induces a new push-forward map on . Let by given by
[TABLE]
with inverse given by
[TABLE]
so that
[TABLE]
Note that is a time dependent operator. Our assumptions imply that both pairs and are compatible uniformly in (Def. 2.2). For example, applying that is compatible and L1 gives
[TABLE]
We use the notation and for the spaces of smoothly evolving functions, each with respect to the pull back (c.f. 2.2). We recall that the definition of and only does not depend on the choice of push-forward map up to norm equivalence (Rem. 2.5). We assume the following inclusions hold:
[TABLE]
3.2.3. Lifted material derivative
We denote by the material derivative for the push-forward map 2.4:
[TABLE]
This is a different material derivative to the material derivative defined with respect to the push-forward map 3.3. However, as observed in Dziuk and Elliott (2013a), our construction implies that first taking material derivatives and then lifting is the same as first lifting and then taking material derivatives:
Lemma 3.5** (Commutation of material derivative and lifting).**
The following hold:
- •
For all
[TABLE]
- •
* if, and only if, , and if, and only if, .*
Proof.
Indeed, applying the definitions 3.13 and 3.14, we see
[TABLE]
since the lift at time and time derivative commute and and are inverses. By similar reasoning, we have that
[TABLE]
In particular, applying L1 and L2, we have the second statement. ∎
3.2.4. Abstract lifted transport formulae
We assume that we have a transport formula for functions in and for the and bilinear forms. We assume that there exists bilinear forms and such that
[TABLE]
We assume that for these bilinear forms there exists such that for all and all ,
[TABLE]
3.2.5. Lifted finite element space
Let . The lifting process allows us to introduce a new space by
[TABLE]
This is a subspace of for each .
The lifted discrete space is a closed subspace of and so we may infer that is a compatible pair (Rem. 2.3). Hence, the spaces 2.1 and 2.2 are well defined and we will write for the material derivative of . From our assumptions, in general it does not hold that .
3.2.6. Approximation property of
For each , we assume that there exists a well defined interpolation operator such that there exists a constant such that for all
[TABLE]
3.2.7. Assumptions on the geometric approximation
Finally, we assume we have the following relations between continuous and discrete bilinear forms. We assume that there exists constants such that for all the following holds for all with lifts we have
[TABLE]
For with inverse lifts , we have
[TABLE]
For and , with inverse lifts and , we have
[TABLE]
Finally, we assume
[TABLE]
3.3. Ritz projection
It is convenient to introduce a Ritz projection which is a standard approach in the finite element analysis of evolution equations (Thomée, 2006) also applied to problems on evolving surfaces (e.g. Dziuk and Elliott, 2013a). The Ritz projection is defined with respect to modified positive definite bilinear forms and .
3.3.1. A new bilinear form
We know from Assumptions M2, Mh2, A3, An4 and Ah2, there exists such that there exists such that for all and all
[TABLE]
We now take fixed in the sequel. Thus, we infer that the bilinear forms:
[TABLE]
are uniformly coercive for all .
3.3.2. The projection
Definition 3.6*.*
The Ritz projection is an operator . For , is given as the unique solution of
[TABLE]
We denote by .
We will further assume that there exists a constant such that for all , for and all we have
[TABLE]
Remark 3.7*.*
In an application it may be possible to prove that a simpler version of B3 is sufficient. For example in the case of of surfaces without boundary for all , for and all we have
[TABLE]
holds.
3.3.3. A dual problem
We introduce the dual problem: Given , find such that
[TABLE]
Assumptions on along with the previous assumptions from Sec. 2.2 imply that 3.21 has a unique solution and we assume the regularity condition that there exists such that
[TABLE]
where the constant is independent of and time .
3.3.4. Ritz error analysis
Lemma 3.8**.**
For each , there exists a unique solution of 3.20. There exists a constant such that for all and all we have
[TABLE]
Furthermore, there exists a constant such that for all and ,
[TABLE]
Proof.
Since is uniformly coercive 3.19 and bounded (Ah3 and Mh2) and is bounded (A4 and M2), standard Lax-Milgram theory gives that there exists a unique solution that satisfies the stability bound 3.22.
To show the error bound, we consider the functional given by
[TABLE]
First, note that for , we can use the definition of 3.20 to see that
[TABLE]
Then the perturbation estimates P1 and P4 and the stability bound 3.22 imply that
[TABLE]
Next, we consider for . Then, again using 3.20 we have
[TABLE]
Using the boundedness of (A4 and M2) and the interpolation bounds I1, we have
[TABLE]
We split so that together with the perturbation estimates P1, P4 and P4’ and the interpolation result we have
[TABLE]
Then combining the above estimates with the stability bound 3.22, we see that
[TABLE]
To show the -norm error bound, we have
[TABLE]
Applying the boundedness (A4 and M2) and coercivity 3.18 of , the interpolation bound I2 and the first bound on 3.24 gives
[TABLE]
Using the interpolation bound I2 and rearranging using a Young’s inequality gives
[TABLE]
For the -norm bound, we consider the dual problem 3.21 with . Then there exists a unique such that
[TABLE]
Furthermore, and satisfies R2
[TABLE]
Then we have from M2 that
[TABLE]
Then the second bound on 3.25 together with the -norm bound 3.26 and the dual regularity estimate 3.27 imply that
[TABLE]
Rearranging this inequality provides the -norm bound. ∎
3.3.5. Time derivative of Ritz projection
Since in general the material derivative and Ritz projection do not commute, we must provide a further estimate for this material derivative of the error . First, we derive an equation for .
Lemma 3.9**.**
Let then, for each , we have the equation:
[TABLE]
Proof.
First consider and take the time derivative of 3.20. We apply the discrete transport formulae Gh1 and Bh1 on the left hand side and the lifted transport formulae Bℓ1 and Gℓ1 on the right hand side and rearrange to see:
[TABLE]
Using the fact that for , we have and , we can apply 3.20 once more to see that
[TABLE]
We can expand this result to arbitrary by considering the function given by which satisfies and . ∎
Lemma 3.10**.**
For , we have that and there exists a constants such that for all ,
[TABLE]
Furthermore, if , then
[TABLE]
Proof.
For the stability bound, we see that satisfies the discrete elliptic problem 3.28. This tells us that and, combined with the boundedness of (A4 and M2), (Bℓ2 and Gℓ2), and (Bh2 and Gh2) we see that
[TABLE]
Then applying the perturbation estimate P9 and the stability of estimate 3.22, we see 3.29.
To show the error bound, we proceed in a similar fashion to Lem. 3.8, we introduce the functional given by
[TABLE]
First, for , we can use 3.28 to see
[TABLE]
Then, using the perturbation estimates on (P1 and P4), with the fact that the two discrete material derivatives and lifting commute 3.16, and (P2 and P5), the boundedness of (Gℓ2 and Bℓ2), the error bound 3.23 and the stability estimates 3.22 and 3.29 gives
[TABLE]
Secondly, for , we have using 3.28
[TABLE]
We split the four terms using the smooth functions and so that
[TABLE]
Using the boundedness of (A4 and M2) and the interpolation estimate I1, we have
[TABLE]
Using the perturbation errors for P1, P4 and P4’ (with the fact that discrete material derivatives and lifting commute 3.16, the inclusions shown in L3 and • ‣ 3.5), as well as the estimate with material derivatives P7, together with the interpolation bound I1 and the error in material derivatives P9, we have
[TABLE]
Using the simple and improved perturbation estimate for P2, P5 and P5’, the interpolation result I1, and the Ritz -norm error bound 3.23, we have
[TABLE]
Using the boundedness of (Bℓ2 and Gℓ2), the perturbation estimates P6 and P3 on , the boundedness of (B3 and G2) and the Ritz and -norm error bounds 3.23 we have
[TABLE]
Combining the previous four bounds with the stability estimates for 3.22 and 3.29 gives
[TABLE]
To show the -norm error bound, we start with
[TABLE]
Noting that the final term on the right hand side is , the bounds on (A4 and M2), the perturbation estimate P9, the interpolation estimate I2 and the first bound on 3.31 gives
[TABLE]
Again, using the interpolation bound I2 and the coercivity of 3.18 and rearranging using a Young’s inequality gives
[TABLE]
To show the -norm bound, we consider the dual problem 3.21 with . Then, there exists such that
[TABLE]
Furthermore, and satisfies the bound
[TABLE]
Then we have
[TABLE]
The second bound on 3.32, the -norm error bound 3.33 and the dual regularity result 3.34 give
[TABLE]
Rearranging this inequality gives the desired -norm bound. ∎
3.4. Abstract error bound
To show the error bound we make the following assumption on the smoothness of the continuous problem. We assume that and that there exists a constant such that satisfies that regularity estimate
[TABLE]
Theorem 3.11**.**
Let all the assumptions listed in Sec. 3.2 hold as well as B3 and R1. Denote by the solution of 2.9 and by the solution of 3.7 with lift . Then, there exists constant such that for we have the error estimate
[TABLE]
To show the error bound, we start by rescaling both solutions. Let and , which satisfy
[TABLE]
Our assumptions imply
[TABLE]
We define to be the lift of . We will decompose the error as:
[TABLE]
We already have bounds on from Lem. 3.8 and 3.10, thanks to assumption R1, so it remains to show a bound for . We will denote by , and by our assumptions, we know .
Lemma 3.12**.**
Let and denote by . Then satisfies
[TABLE]
where
[TABLE]
Proof.
We transform 3.36 into variational form using G1, then together with the definition of the Ritz projection 3.20, we see that
[TABLE]
We use the transport formulae Gh1 for and Gℓ1 for to see
[TABLE]
Subtracting this equation from 3.37 and rearranging gives 3.40. ∎
Lemma 3.13**.**
For , the consistency terms and satisfy
[TABLE]
Proof.
For , we use M2 and 2.6 together with the error bounds from 3.23 and 3.30 to see
[TABLE]
For , we use the perturbation estimates P1, P2 and P8 together with the stability bounds on the Ritz projection 3.22 and 3.29 to see
[TABLE]
Lemma 3.14**.**
The following bound holds for
[TABLE]
Proof.
We test 3.40 with to see
[TABLE]
The transport formula for Gh1 tells us that
[TABLE]
hence, applying the bound on and 3.41 we infer that
[TABLE]
Applying the boundedness and coercivity estimates form and Mh2, 3.19 and Gh2 with a Young’s inequality and integrating in time gives
[TABLE]
Finally, we use a Grönwall inequality to see the desired result. ∎
Finally, we can show the result of Thm. 3.11.
Proof of Thm. 3.11.
We apply the splitting 3.39, the bounds on from Lem. 3.8, the bounds on from Lem. 3.14 and the estimate on from 3.38 to see
[TABLE]
The final line follows since the norm is bounded by the norm. ∎
Part II Evolving finite element spaces
4. Evolving bulk finite element spaces
In this section, we will define families of evolving bulk finite element spaces on families of evolving triangulated bulk domains consisting of unions of elements. We have in mind approximating an open bounded domain . We will use the word bulk in this section to emphasise the difference to the surface case considered in Sec. 6 but more common terminology would simply remove this word.
Our work extends from standard bulk finite element theory (Ciarlet, 1978) and the work of Ciarlet and Raviart (1972) and Bernardi (1989) for Cartesian bulk domains with curved boundaries to the evolving case. Throughout this section we will denote global discrete quantities with a subscript , which is related to element size. We assume implicitly that these structures exist for each in this range. See also Rem. 4.11. For ease of exposition we begin with definitions without the time parameter, .
4.1. Reference finite element
Definition 4.1* (Reference finite element).*
The triple is a reference finite element if:
- (a)
the element domain is the closure of an open domain with Lipschitz piecewise smooth boundary, 2. (b)
the set of shape functions is a finite dimensional space of functions over , 3. (c)
the nodal variables or degrees of freedom are a basis of the dual space to ,
and determines , that is if for with for all , we have .
As part of this definition, we are implicitly assuming that the nodal variables live in the dual to a larger function space than . We will see that this usually requires further smoothness or continuity of finite element functions. We give an example of a simplical finite element, but this definition includes other examples such as isoparametric finite elements and brick finite elements.
Recall that a (non-degenerate) -simplex in , is the convex hull of distinct points , called the vertices of the -simplex, which are not contained in a common -dimension hyperplane. More precisely, we have
[TABLE]
For each , we call barycentric coordinates. For any integer with , an -facet of an -simplex is any -simplex whose vertices are also vertices of . We call an -facet a boundary facet. We will also use the term boundary facet for any boundary polytopes (union of simplicies) of a polytope . For each , we shall denote by the space of all polynomials of degree in the variables in . For any set , we let
[TABLE]
Example 4.2* (Example of reference finite element).*
The standard piecewise linear finite element is obtained by choosing to be a non-degenerate -simplex in and . We can also define higher order spaces , for , by including extra evaluation points in (see, for example, Ciarlet (1978, Section 2.2)). The key property of the extra evaluation points is that they determine the particular function in . It is also true that the restriction of to any facet determines the restriction of functions in on that facet.
Given the reference element and a function for which the nodal variables can be computed (e.g. in the case of Lagrange element we require to be continuous), we define the nodal interpolation of , written , as the unique function in which has the same nodal values as . Let be the basis of dual to then we can characterise as
[TABLE]
Lemma 4.3** (Bramble-Hilbert Lemma, Ciarlet 1978, Thm. 3.1.5).**
Let the following inclusions hold for , and ,
[TABLE]
Under the above assumptions on the reference finite element we have that there exists a constant such that for all functions ,
[TABLE]
4.2. Bulk finite element
We start by defining a single bulk finite element in . Our definition of bulk finite element combines Def. 2.1 and 2.3 from (Bernardi, 1989).
Definition 4.4* (Bulk element reference map and bulk finite element).*
Let be a reference finite element (Def. 4.1) with .
- (a)
Let satisfy
-
(1)
-
(a)
; 2. (b)
; 3. (c)
is a bijection onto its image; 2. (2)
can be decomposed into an affine part and smooth part
[TABLE]
such that is an invertible matrix, , and
[TABLE]
where denotes the two-norm of the matrix.
In this situation we call a bulk element reference map. 2. (b)
Let be a bulk element reference map and be the triple given by
[TABLE]
Under the above assumptions, we call a bulk finite element, the associated reference finite element.
With the bulk reference element map we can compute integrals and derivatives over the reference element using the transformation identity:
[TABLE]
We denote by the outward pointing normal to .
Definition 4.5* (-bulk finite element, Def. 2.4 Bernardi 1989).*
Let and be the bulk element reference map for a bulk finite element .
- (a)
We say that is a -bulk finite element reference map if
- (i)
the bulk element reference map ; 2. (ii)
for , there exists constants such that
[TABLE] 2. (b)
We say that is a -bulk finite element if is a -bulk finite element reference map and
- (i)
the space contains the functions for all ; 2. (ii)
the space is contained in .
Remark 4.6*.*
The properties of allow us to define the Sobolev spaces for , . Since and is finite dimensional, we clearly see that is a closed subspace of for , .
Example 4.7* (Bulk finite elements).*
We are thinking of two particular examples. The first is a standard Lagrange finite element and the second is an isoparametric finite element. Examples of each of these cases are shown in Fig. 4.1.
- (a)
Let be a reference finite element. Consider the affine map given by . If is non-singular, then this defines a bulk finite element . In the standard way, the element domain is defined by the location of its vertices. For a simplex reference element domain , we are thinking of line segments in , triangles in and tetrahedra in . 2. (b)
Let be a reference finite element. Let be a bulk finite element which is the image of under a map which satisfies . We call an isoparametric bulk finite element. We note that the functions in will not necessarily consist of polynomials over even if consists of polynomials over , however this leads to a practical scheme where integrals are computed over reference elements. This example is the basis for the method in Sec. 8.
The definition of bulk finite element (Def. 4.4), in particular 4.2, is constructed to allow the following result:
Lemma 4.8** (Lem. 2.1, Bernardi 1989).**
Let be a bulk finite element reference map then is a -diffeomorphism and satisfies
[TABLE]
and also for all
[TABLE]
To help us understand the geometry of the new element domains , we introduce a new element domain defined by the affine part of the parametrisation: .
Lemma 4.9**.**
Let be an element domain 4.3a parametrised by a bulk element reference map over . Denote by
[TABLE]
We will also write and for the diameter of and diameter of the maximum inscribed ball in . Then we have that
[TABLE]
Proof.
See (Ciarlet, 1978, Thm. 3.1.3). ∎
Remark 4.10*.*
We note that the volume of an element can be estimated by and by
[TABLE]
Here the positive constants depend on the volume of the unit ball in and the constant .
Remark 4.11*.*
In the sequel, we will assume implicitly that results hold for sufficiently small () for some particular value of . In general this is always possible by subdividing a particular element using a refinement procedure and applying the result to the subdivided, smaller elements.
The choice of mapping allows us to relate functions defined on to functions on .
Lemma 4.12** (Lem. 2.3, Bernardi 1989).**
Let be a -bulk finite element reference map (Def. 4.5). Let and , then implies belongs to . We have for any that
[TABLE]
for a constant which depends on . We also have for any that and
[TABLE]
where the constant here depends on and the product .
Given a bulk finite element (Def. 4.4), let be the basis dual to . This is the set of basis functions of the finite element. If is a function for which all , is well defined, then we define the local interpolant by
[TABLE]
We can think of as the unique shape function that has the same nodal values as so that, in particular, for .
Theorem 4.13** (Local interpolation estimate).**
Let be a -bulk finite element (Def. 4.5) with reference element by which satisfies the assumptions of Lem. 4.3 for some , . Then there exists a constant such that for all functions
[TABLE]
Proof.
We re-scale 4.1 using Lem. 4.12 and the estimates from 4.9:
[TABLE]
The last line holds if (note that so this statement is true for small enough). ∎
4.3. Triangulated bulk domain and spaces
Next, we bring together a finite family of bulk finite elements in order to create a bulk finite element space.
Definition 4.14*.*
- (a)
A triangulated (bulk) domain is a set equipped with an admissible subdivision consisting of bulk finite element domains such that , for with . 2. (b)
The maximum subdivision diameter is defined by:
[TABLE] 3. (c)
Let be a discrete bulk domain equipped with an admissible subdivision such that each set is an element domain for a bulk finite element parametrised over the same polygonal reference finite element . We say that is a facet if is the image of a boundary facet of . 4. (d)
We say that is a conforming subdivision of if any facet of an element domain is either a facet of another element domain , in which case we say and are adjacent, or a portion of the boundary . 5. (e)
For a conforming subdivision, we denote by the set of facets between adjacent elements and by the set of boundary facets. For a internal facet between adjacent elements , we make a choice of the two elements and denote by the outward normal to and by the outward normal to . The particular choice will not affect our calculations.
Remark 4.15*.*
- (a)
We recall that our elements domains (Def. 4.1(a)) are closed so that is also a closed set. 2. (b)
When putting together bulk finite elements in order to form a discrete domain , we will be generally thinking of the case that only elements with more than one vertex on the boundary are curved (). See Sec. 8 for more details. We allow for the more general case here.
Definition 4.16* (Broken Sobolev spaces and norms).*
Let be a subdivision of consisting of -bulk finite elements. Then for , , we define the broken Sobolev space by
[TABLE]
with norm
[TABLE]
Remark 4.17*.*
This space is often used in the context of discontinuous Galerkin finite element methods. See, for example, Arnold et al. (2002).
Lemma 4.18**.**
The space is complete.
Proof.
Consider a Cauchy sequence . This implies
- •
is Cauchy in so there exists such that in .
- •
is Cauchy for all so there exists such that in .
It is clear from the triangle inequality that :
[TABLE]
since the right hand-side converges to [math] as . Hence, we have shown that converges to a function in the norm. ∎
Let be a triangulated domain with conforming subdivision . For , we denote the trace of a function by and recall that there exists a constant such that
[TABLE]
We define the space by
[TABLE]
We equip this space with the broken norm .
Lemma 4.19**.**
The space is a closed subspace of so is complete.
Proof.
Take a sequence which converges to . Then for any pair of adjacent elements we have
[TABLE]
Clearly the right hand side converges to [math] as so we have the traces of from adjacent elements coincide and . ∎
We will use the notation for which is a Hilbert space when equipped with the obvious broken inner product.
Lemma 4.20**.**
Let be a conforming subdivision of then .
Proof.
First, let . Then we have that and it is left to show that has a weak derivative in . We have a candidate given element-wise by . It is clear that and for and , we have
[TABLE]
We note that we can write
[TABLE]
where is the set of facets between adjacent elements in , since the traces of and to coincide for any adjacent pair . We note that the sum over edges is zero since the normals from adjacent elements are equal and opposite in a conforming triangulation and the integral over boundary facets is zero since is zero here. Thus, we see that is the weak derivative of .
Second, let then it is clear that , and by the trace theorem has common trace between adjacent elements. ∎
Remark 4.21*.*
The proof is based on ideas from Ciarlet (1978, Thm. 2.1.1) which showed that appropriate finite element spaces defined in that work are contained in (in our notation).
4.3.1. Bulk finite element space
We restrict to Lagrangian finite elements over a polygonal reference finite element. More precisely, we assume that the degrees of freedom for each element are given by
[TABLE]
where is a finite set of nodes in . We call the set of Lagrange nodes of .
The set of degrees of freedom of adjacent bulk finite elements will be related as follows. Let and be two bulk finite elements such that and are adjacent with and . Then, we have
[TABLE]
We denote the global set of Lagrange nodes by
[TABLE]
For each , let be the local neighbourhood of elements for which .
Definition 4.22* (Bulk finite element space).*
- (a)
Let be a discrete bulk domain equipped with a conforming subdivision with each domain equipped with a bulk finite element (Def. 4.4) which satisfy 4.19. A bulk finite element space is a (generally proper) subset of the product space given by
[TABLE] 2. (b)
The bulk finite element space is determined by the global degrees of freedom
[TABLE]
In this definition, an element is not, in general a “function” defined over , since we do not necessarily have a good definition of over element boundaries: The “function” may be double-valued.
If it happens, however, that for each element , the restrictions and coincide along the common face of any adjacent elements and , then the function can be identified with a function defined over the set . In this case, we call the elements bulk finite element functions. Examples of bulk finite element functions are shown in Fig. 4.2.
We enumerate the nodes so that and take to be the basis of dual to . Since, we have a finite basis of we note that we can identify any with a vector so that
[TABLE]
Lemma 4.23**.**
Let be a bulk finite element space consisting of bulk finite elements over a conforming subdivision of . Assume further that for each , the corresponding reference finite element is a Lagrange element of order (Ex. 4.2). Then we can identify elements of as functions in . Furthermore is a closed subspace in .
Proof.
Consider two adjacent elements and . The functions and when restricted to the appropriate edges in are polynomials of degree which agree at the Lagrange points on this edge from 6.20 and the definition of . The Lagrange points in the reference element determine polynomials of degree so we have that on . Since is a conforming subdivision we can define a global function such that for each which is globally continuous. Indeed restricted to each element is continuous, as the composition of a polynomial (element of ) and a smooth surface finite element reference map , and is single valued on the facets where any two elements meet.
In fact the restriction to each element is a -function hence so it is clear that . The space is closed since it is finite dimensional. ∎
Remark 4.24*.*
The proof is based on ideas from Ciarlet (1978, Thm. 2.2.3) which show that appropriate finite element spaces defined in that work are contained in (in our notation).
The approximation property of the finite element space will be defined through an interpolation operator:
Definition 4.25* (Interpolation).*
If is a function on for which all , , is well defined (in case of Lagrangian finite elements, suffices), then we can define a global interpolant by
[TABLE]
Note that our construction implies that
[TABLE]
and for all .
In order to prove estimates on the global interpolant, we will first define three further properties of our subdivision .
Definition 4.26* (Regular and quasi-uniform subdivisions).*
For , let be a triangulated bulk domain (Def. 4.14(a)) equipped with a conforming subdivision (Def. 4.14(d)).
- (a)
[Def. 3.1, Bernardi 1989] The family is said to be non-degenerate or regular if there exists such that for all and all ,
[TABLE]
and there exists a constant such that
[TABLE] 2. (b)
[Def. 3.2, Bernardi 1989] The family is said to be -regular if it is regular, if for all and all , , and if, there exists a constant such that
[TABLE] 3. (c)
A regular family is said to be quasi-uniform if there exists such that
[TABLE]
Remark 4.27*.*
We note that:
- •
for a regular subdivision there exists a constant depending on the global quantities and
[TABLE]
- •
for a quasi-uniform subdivision there exists a constant depending on the global quantities and
[TABLE]
Theorem 4.28** (Global interpolation estimates, c.f. Cor. 4.1, Bernardi 1989).**
For , let be a triangulated bulk domain (Def. 4.14(a)) equipped with a -regular (Def. 4.26(b)), quasi-uniform (Def. 4.26(c)), conforming (Def. 4.14(d)) subdivision . Let each be equipped with a -bulk finite element (Def. 6.4) parametrised over a reference finite element which satisfies the assumptions of Lem. 4.3 for some , . Then there exists a constant such that for all functions ,
[TABLE]
Proof.
The proof follows by piecing together Thm. 4.13 using the fact that is quasi-uniform. ∎
Remark 4.29*.*
The approximation property shown in (Bernardi, 1989, Cor. 4.1) is a result for an projection for a more general class of finite element spaces.
4.4. Evolving bulk finite elements
Let denote time. We consider families of bulk finite elements, spaces and triangulated domains parametrised by .
Definition 4.30* (Evolving bulk finite element).*
- (a)
Let be a time dependent family of bulk finite elements (Def. 4.4) parametrised over a common reference element . If the constant is uniformly bounded away from 1,
[TABLE]
we say that is an evolving bulk finite element. 2. (b)
Let where . We say that if is such that
[TABLE]
then is the flow defining the evolution of the element domain and that is the evolving bulk element reference map. 3. (c)
The element velocity of is defined by
[TABLE] 4. (d)
If each is a -bulk finite element for each and the constants are uniformly bounded:
[TABLE]
then we say that is an evolving -bulk finite element and is an evolving -bulk element reference map. 5. (e)
We say that an evolving bulk finite element is temporally quasi-uniform, if there exists such that
[TABLE] 6. (f)
The family of element push forward maps denoted by for , indexed by , is defined to be the linear bijections defined by:
[TABLE]
Lemma 4.31**.**
Let , for be an evolving -bulk finite element reference map (Def. 4.30(d)) for a temporally quasi-uniform element domain (Def. 4.30(e)) and the family of element push forward maps (Def. 4.30(f)). Then there exists constants , which depend only on the reference element and the constants and , such that for all , if and only if and
[TABLE]
Proof.
From Lem. 4.12 and 4.9, we have
[TABLE]
and
[TABLE]
It can be easily seen that for a quasi-uniform evolving surface finite element that these constants only depend on allowed quantities. ∎
This result implies that is a compatible pair (Def. 2.2). Furthermore , equipped with the norm is also a compatible pair since is a closed subspace of (Rem. 2.3 and 4.6).
4.5. Evolving bulk triangulations and spaces
We now derive definitions of an evolving bulk finite element space which is part of a compatible pair (in the sense of Sec. 2, Def. 2.2). For each , we are given a family of discrete bulk domains and each equipped with a bulk finite element space . Furthermore, we are interested in under what assumptions does the compatibility hold independently of the element diameter .
Definition 4.32* (Evolving bulk domain).*
For , let be a family of triangulated bulk domains (Def. 4.14(a)) each equipped with a conforming subdivision (Def. 4.14(d)) such that each element domain is equipped with an element flow map (Def. 4.30(b)).
- (a)
We call an evolving conforming subdivision if for each element and each facet of either is a facet of another element , in which case is a common facet between and for all or is a portion of the boundary , in which case is a portion of the boundary for all . 2. (b)
An evolving triangulated (bulk) domain is defined to be a family of triangulated bulk domains equipped with an evolving conforming subdivision. In this case, we define the mesh parameter to be
[TABLE] 3. (c)
We define a global discrete flow element-wise by
[TABLE]
Our assumptions imply that is piecewise smooth and . 4. (d)
We define a global discrete velocity given by
[TABLE] 5. (e)
The family of linear bijections induced by the flow and called the global push forward map is denoted by for and defined by
[TABLE]
In order to bring together a collection of elements, we again restrict to Lagrangian finite elements over a polygonal reference finite element. We note that our construction implies that for each element each node is the trajectory of a point under the flow - that is .
We make the following extra requirement on how adjacent elements will be related. For each , we denote by the global set of Lagrange nodes of 4.20 and for any , is the set of elements such that is a node of . We make the restriction that the global flow is single-valued at each Lagrange point: for all we have
[TABLE]
Remark 4.33*.*
We note that our construction does not imply that the global flow map is indeed a function: Its restrictions to a facet from adjacent elements may not coincide. The assumption 4.25 imposes that the global flow map should coincide at Lagrange points along element boundaries.
Definition 4.34* (Evolving bulk finite element space).*
- (a)
Let be an evolving triangulated bulk domain (Def. 4.32(b)) equipped with an evolving conforming subdivision (Def. 4.32(a)). For , let be a bulk finite element space (Def. 4.22(a)) over . If each is equipped with an evolving bulk finite element (Def. 4.30(a)) which together satisfy 4.25 then we say is an evolving bulk finite element space. 2. (b)
For each , we will write for the set of global nodal variables (Def. 4.22(b)). We will use the convention that
[TABLE]
where is the trajectory of a Lagrange point under the global flow . We will denote by the global basis of finite element functions such that for and all . This implies that .
Definition 4.35* (Uniformly regular and uniformly quasi-uniform evolving subdivisions).*
For , let be a family of evolving conforming subdivisions (Def. 4.32(a)).
- (a)
We say that the family is uniformly regular if there exists such that for all and all times , we have
[TABLE]
and there exists such that
[TABLE] 2. (b)
We say that the family is uniformly -regular if it is uniformly regular, if for each time , the family is -regular and if there exists a constant such that
[TABLE] 3. (c)
We say that the family is uniformly quasi-uniform if there exists such that for all and all times , we have
[TABLE]
Note that a uniformly quasi-uniform subdivision consists of element domains for temporally quasi-uniform evolving bulk finite elements.
Lemma 4.36**.**
For , let be a uniformly -regular (Def. 4.35(b)), uniformly quasi-uniform (Def. 4.35(c)), evolving, conforming subdivision (Def. 4.32(a)) and let be the global push-forward map (Def. 4.32(e)). Let , . Then, if and only if for all . Furthermore, there exists independent of and such that for all
[TABLE]
Proof.
We simply sum the element-wise result from Lem. 4.31. The constants are independent of and due to the uniform quasi-uniformity of . ∎
Remark 4.37*.*
In particular, the pair is compatible with respect to the broken Sobolev norm (Def. 4.16). Furthermore, the pairs , equipped with the norm, and are also both compatible (Rem. 2.3, 4.23 and 4.19). Note that this result implies that the spaces and are well defined when equipped with the appropriate norms (c.f. 2.1 and 2.2).
5. Lifted bulk finite element spaces
This section sets out a procedure for relating functions on the discrete bulk domain to the smooth domain via the construction of from the space defined in the previous section. We start by defining a lifted bulk finite element in using a mapping . This process is called lifting. In the following we will provide the appropriate assumptions on to allow us to relate structures on to their lifted counterparts on .
5.1. Lifted bulk finite element
Consider a single bulk finite element (Def. 4.4), with element reference map , where the element domain approximates a portion of a domain with smooth boundary in . Let be a -map which is a diffeomorphism onto its image. We define . For a function , we call the lift of which is given by
[TABLE]
We assume that we can decompose into
[TABLE]
where is an invertible matrix, and . We will assume that does not affect the affine part of the parametrisation: for each vertex .
Definition 5.1* (Lifted bulk finite element).*
We call the triple defined by
[TABLE]
the lift of and the lifting map. If forms a bulk finite element over then we say that is the lifted bulk finite element associated with . In this case, we call the lifted bulk finite element reference map.
The next two results show under what assumptions on is a bulk finite element (Def. 4.4) or a -bulk finite element (Def. 4.5).
Lemma 5.2**.**
If satisfies that
[TABLE]
then is a bulk finite element. Furthermore, for , we have that there exists such that
[TABLE]
where the constants depend on , and the ratio .
Proof.
To show that is a bulk finite element (Def. 4.4), the conditions on the element reference map are clear and we are left to check the curvedness condition 4.2. Using the expansion of , we see
[TABLE]
So that
[TABLE]
Applying the curvedness condition for we see
[TABLE]
The curvedness condition is shown by applying 5.1.
To show 5.2a the result is clear for . For , we will apply Lem. 6.12. Then we see
[TABLE]
where depends on and the bound 5.1. For the bound 5.2b, we note that
[TABLE]
so that we infer
[TABLE]
Then applying Lem. 4.12 once more we see
[TABLE]
The final result is given by applying 4.9a and 4.9b. ∎
Lemma 5.3**.**
If is a -bulk finite element and satisfies and , then is also a -bulk finite element. Furthermore, for and , we have that there exists such that
[TABLE]
where the constants depend on , and the ratio .
Proof.
To see that is a -bulk finite element (Def. 4.5), the first three conditions are clear from the smoothness assumption on . We must show 4.4 holds: for , there exists a constant such that
[TABLE]
Computing directly using (Bernardi, 1989, Eq. (2.9)), we have
[TABLE]
where are constants and
[TABLE]
But applying the fact that is a -bulk finite element, the definition of , Lem. 4.12 and the smoothness assumption on , we see that
[TABLE]
Finally, applying 5.3, we have that
[TABLE]
This shows that is a -surface finite element.
To show the norm equivalence we again apply Lem. 4.12, recognising the geometric progression, to see
[TABLE]
where
[TABLE]
The final result is given by applying 4.9a, 4.9b and 5.4c. ∎
Remark 5.4*.*
Considering the case that fixes the vertices of then we can write as
[TABLE]
That is that and . Then the assumptions of Lem. 5.2 can be replaced by
[TABLE]
and the assumptions of Lem. 5.3 can be replaced by the assumption that is bounded.
We can also use to define an inverse lift. Given an element domain , lifted element domain and lifting map . We know that is invertible onto its image, namely . So for , we denote its inverse lift by defined by
[TABLE]
Lemma 5.5**.**
If the assumptions of Lem. 5.2 and 5.3 hold then, for and , we have that there exists such that for all , we have
[TABLE]
Proof.
The same proof can be applied to the case of the inverse lift as well as the lift. ∎
We next relate the geometry of the base and lifted element domains.
Lemma 5.6**.**
Using the decomposition of , we have
[TABLE]
Proof.
For 5.4a, we show the second inequality. The first follows by the same reasoning applied to the inverse of . Since is compact, there exists such that
[TABLE]
But since is invertible, there exists such that
[TABLE]
Then, we can compute that
[TABLE]
Similarly, for 5.4b, we only show the first inequality. For each , there exists a ball in such that
[TABLE]
Denote by and the centre and radius of respectively. Consider the affine map given by
[TABLE]
It is clear that maps balls to balls and is mapped to a ball centred at with radius . We claim is contained in . Indeed, take , then there exists such that . Denote by
[TABLE]
then
[TABLE]
so and
[TABLE]
so that .
Therefore, we have found a ball with radius , thus we can infer that
[TABLE]
Since the proof holds for all , we see the desired result.
Given 5.4a and 5.4b, the final result 5.4c follows directly from Rem. 4.10. ∎
5.2. Lifted bulk triangulations and spaces
Let be a smooth bulk domain and for and let be a triangulated bulk domain (Def. 4.14(a)) equipped with a conforming subdivision (Def. 4.14(d)) and a bulk finite element space (Def. 4.22(a)). Let each be associated with a lifted finite element with lifted map .
Definition 5.7*.*
- (a)
We denote by the set of all lifted element domains
[TABLE]
If the global map is single valued and forms a conforming subdivision of , we say that is an exact subdivision of the domain . 2. (b)
We define a global lifting map by . We define the inverse lift in a similar element-wise fashion by . 3. (c)
We denote by , for , the global lift given by
[TABLE]
and by , for ,, the global inverse lift given by
[TABLE]
We will also use the notations and . 4. (d)
Let be a bulk finite element space. If for each bulk finite element there is an associated lifted bulk finite element , then, we define a lifted bulk finite element space (c.f. 3.17) by
[TABLE]
Proposition 5.8**.**
Assume additionally that the family of triangulations is regular (Def. 4.26(a)), satisfies the assumptions of Lem. 5.2 and there exists such that
[TABLE]
for all for all . Then is regular and is a bulk finite element space and there exists constants , which are independent of such that
[TABLE]
Furthermore, if satisfies the assumptions of Lem. 5.3 for all and all , then for , there exists independent of such that
[TABLE]
Proof.
The regularity of the family of triangulations follows from 5.4a and 5.4b from assumption 5.5.
The results follow by combining the previous results for each element in . We achieve bounds independently of since the regularity of the subdivisions implies that is bounded independently of . ∎
We next define interpolation estimate which interpolates smooth functions over the continuous surface into the lifted surface finite element space. We denote by an interpolation operator defined by
[TABLE]
Theorem 5.9** (Global lifted interpolation theorem).**
For , let be a -regular (Def. 4.26(b)), quasi-uniform (Def. 4.26(c)) family of subdivisions of triangulated bulk domains equipped with a bulk finite element space (Def. 4.22(a)) consisting of -bulk finite elements (Def. 4.5) over a reference element which satisfies Lem. 4.3 for some , . Let each element be equipped with a lifting map such that and 5.5 and 5.1 hold each uniformly for . Then is a -regular, quasi-uniform family of subdivisions of and is a bulk finite element space consisting of -bulk finite elements. Let be a continuous function, then is well defined. Furthermore, if the assumptions of Thm. 4.13 hold for the reference element , there exists a constant such that for all functions ,
[TABLE]
Proof.
The result follows by combining previous lemmas in the appropriate way. We see the lifted triangulation is quasi-uniform by applying the results in Lem. 5.6 and the assumption 5.5. The interpolation result then follows by applying Thm. 4.28. ∎
Corollary 5.10**.**
Let , then the interpolant of into , denoted by , given by
[TABLE]
Furthermore, there exists a constant such that for all functions , and all , we have
[TABLE]
Proof.
We apply the inverse lift result to the estimates in the theorem. ∎
5.3. Evolving lifted bulk triangulations
For , let be a smoothly evolving bulk domain with flow map defined for the closure : i.e. , and . For , let be an evolving triangulated bulk domain (Def. 4.32(b)) with global discrete flow (Def. 4.32(c)) and equipped with an evolving conforming subdivision (Def. 4.32(a)) and an evolving bulk finite element space (Def. 4.34(a)). We assume that we are given a global lifting map (Def. 5.7(b)) which gives an exact subdivision of (Def. 5.7(a)).
Definition 5.11* (Lifted discrete flow map, material velocity and pushed forward map).*
- (a)
The lifted flow map of the smooth bulk domain is defined by
[TABLE]
We note that . 2. (b)
The lifted discrete material velocity on is defined by
[TABLE] 3. (c)
The family of lifted push forward maps denoted by , for , indexed by , are the linear bijections defined by
[TABLE]
Remark 5.12*.*
Note that in general is different to , but each describes a different parametrisation of the same evolving domain. Also , the material velocity of , and define the same domain evolving from , so have the same normal components on the boundary .
Proposition 5.13**.**
For , let be an evolving bulk finite element space over a -regular (Def. 4.26(b)), uniformly quasi-uniform (Def. 4.35(c)), evolving conforming subdivision consisting of -evolving bulk finite elements (Def. 4.30(d)) and the global push-forward map (Def. 4.32(e)). For each and each , let each element be equipped with a lifting map such that and 5.5 and 5.1 hold each uniformly for and . Then is also an evolving bulk finite element space over a -regular, uniformly quasi-uniform, evolving conforming subdivision consisting of -evolving bulk finite elements. Furthermore, for , , there exists independent of and such that
[TABLE]
In particular, the pair is compatible (Def. 2.2). Furthermore, , equipped with the norm , and are also compatible pairs (Rem. 2.3, 4.23 and 4.19).
Proof.
The properties of follow for the same reasoning as Prop. 5.8 since the constants are bounded uniformly in time. The bounds 5.12 and the compatibility of follow since the assumptions imply that are compatible and we can simply lift this result with uniform bounds. ∎
6. Evolving surface finite element spaces
In this section, we will give precise definitions concerning the evolving surface finite element spaces we use. The ideas follow in a similar manner to Sec. 4. The upshot will be a notion of a discrete surface consisting of a union of elements and a finite element space defined on this discrete domain. Our extensions from standard bulk finite element theory, presented in Sec. 4, to surface finite elements build on the work of Nedelec (1976), Dziuk (1988) and Heine (2005) for surfaces. The subscript parametrises the constructions and will be related to the size of elements used in our computational domain. Implicitly it is assumed that these structures exist for all for some fixed value of . See also Rem. 6.11.
6.1. Surface finite elements
We next define a surface finite element which takes its inspiration from the notion of curved finite elements studied by Ciarlet and Raviart (1972) and Bernardi (1989). The key idea here is that a surface finite element is an -dimensional parametrised surface with boundary embedded in . In the following for a matrix , we use to denote the pseudo-inverse. For any matrix of full column rank the pseudo-inverse is given by
[TABLE]
In this case is the left inverse of : .
Definition 6.1* (Surface finite element and surface element reference map).*
Let be a reference finite element with .
- (a)
Let satisfy
-
(1)
-
(a)
; 2. (b)
; 3. (c)
is a bijection onto its image; 2. (2)
can be decomposed into an affine part, and smooth part
[TABLE]
such that has full column rank, and
[TABLE]
where denotes the two-norm in this context.
In this case, we call a surface element reference map. 2. (b)
Let be a surface element reference map and be the triple given by
[TABLE]
Under the above assumptions, we call a surface finite element, and the associated reference finite element.
Remark 6.2*.*
- (a)
We note that our assumptions imply that both and are full column rank. 2. (b)
The first three assumptions in the definition of surface finite element imply that is a parametrised surface and the fourth 6.1 that is not too curved. The final assumption allows the case that is a flat simplical domain and is curved.
Remark 6.3*.*
- (a)
We denote by the unit normal vector field to . It is the unique (up to sign) unit vector orthogonal to the partial derivatives of for and is given by
[TABLE]
Here denotes the wedge product. The sign of the normal vector field is chosen by fixing a permutation of the barycentric coordinates of the reference element. By swapping any two elements, we reverse the sign of . For a simplex reference element, the orientation can be fixed by ordering the labels of vertices so that where are the vertices of the reference element domain. 2. (b)
We also define the outward pointing unit conormal on the boundary of the element domain . This is the unique (up to sign) vector which is orthogonal to the boundary and the normal . 3. (c)
Our assumptions imply that (or ) is a left inverse but not a right inverse of (respectively, ). We can compute that
[TABLE]
where denotes projection onto the tangent plane to . One way to interpret this result is to note that is the identity operator when restricted to the tangent plane so that in some sense is a right inverse of when we restrict to the tangent plane of .
Definition 6.4* (-surface finite element and -surface element reference map).*
Let and a surface element reference map for a surface finite element .
- (a)
We say that is a -surface finite element reference map if
- (i)
the surface element reference map (i.e. is a -hypersurface); 2. (ii)
for , there exists constants such that
[TABLE] 2. (b)
We say that is a -surface finite element if is a -surface finite element reference map and
- (i)
the space contains the functions for all ; 2. (ii)
the space is contained in .
Remark 6.5*.*
- (a)
The -surface finite element is a generalisation of a curved finite element of order given by Bernardi (1989). See also Def. 4.5. 2. (b)
For a -surface finite element , we have that the Sobolev space is well defined for and (Hebey, 2000; Dziuk and Elliott, 2013b). 3. (c)
Since and is finite dimensional, is a closed subspace of for and . 4. (d)
We note that is a -surface finite element if is a surface finite element (Def. 6.1), the map and contains all affine functions on . We see that the constant (see Lem. 6.8). 5. (e)
Our applications (Sec. 9 and 10) will use in the computational method but we allow the more general case here.
Using transformation formulae, we have for
[TABLE]
where and , and
[TABLE]
and finally, are the components of the inverse . Note that the surface gradient on can also be written as:
[TABLE]
where is the gradient of an arbitrary extension of away from .
Remark 6.6*.*
If the surface finite element map is a function then we can define the extended Weingarten map by
[TABLE]
Example 6.7* (Surface finite elements).*
We are thinking of three particular examples. The first is due to Dziuk (1988) and the second due to Heine (2005). Examples of each of these first two cases are shown in Fig. 6.1.
- (a)
Let be a reference Lagrangian finite element. Consider the affine map given by . If is non-degenerate, then this defines a surface finite element . The element domain is determined by its vertices and consists of affine functions over . This is the surface finite element introduced by Dziuk (1988) which we will call an affine finite element. We think of a simplex for so that the domains are either line segments embedded in , triangles embedded in , and tetrahedra embedded in . 2. (b)
Let be a reference finite element. Let be a surface finite element which the image of under a map which satisfies . We call an isoparametric (surface) finite element. This construction is a generalisation of an affine finite element and was introduced by Heine (2005). We note that the functions in will not necessarily consist of polynomials over even if consists of polynomials over , however this leads to a practical scheme where integrals are computed over reference elements. This example is the basis for the method in Sec. 9. 3. (c)
Let be a reference finite element. Then can be thought of a surface finite element by defining the parametrisation by
[TABLE]
Note that but . In general, we could consider flat surface finite elements to be surface finite elements to have parametrisation such that .
Lemma 6.8**.**
Let be a surface finite element reference map then is a -diffeomorphism and satisfies
[TABLE]
and also for all
[TABLE]
Proof.
The proof of 6.5, 6.6 and 6.7 follows immediately from 6.1 by writing as
[TABLE]
For 6.7, we use the fact that the determinant is an -linear continuous form. ∎
To help us understand the geometry of the new elements, given an element domain , we introduce a new affine element domain defined by the affine part of the parametrisation: .
Lemma 6.9**.**
Let be an element domain 6.2a parametrised by a surface element reference map over . Denote by
[TABLE]
We will also write and for the diameter of and diameter of the maximum inscribed ball in . Then we have that
[TABLE]
Proof.
To show 6.9a, we start by noticing that
[TABLE]
From the definition of we know that for all , , there exists such that . Then noting that , we have that
[TABLE]
Since the choice of was arbitrary, we have shown 6.9a.
For 6.9b, we proceed in a similar fashion with
[TABLE]
Let . We note that has a non-trivial kernel so we decompose where and is in the tangent plane to . Then, we see that there exists such that , and noting that (we see this using the definition of and that is a left inverse of ), we see that
[TABLE]
Again, since the choice of was arbitrary, we have shown 6.9b.
To see 6.9c we apply each of the previous two bounds with the result of Lem. 6.8 to see
[TABLE]
Remark 6.10*.*
We note that the volume of an element can be estimated by and by
[TABLE]
Here the positive constants depend on the volume of the unit ball in and the constant .
Remark 6.11*.*
In the sequel, we will assume implicitly that results hold for sufficiently small () for some particular value of . In general this is always possible by subdividing a particular element using a refinement procedure and applying the result to the subdivided, smaller elements.
This scaling property allows us to characterise Sobolev spaces over a surface finite element and calculate norms over (see Rem. 6.5(b)).
Lemma 6.12**.**
Let be a -reference finite element map Def. 6.4(a). Let and , then implies belongs to . We have for any that
[TABLE]
for a constant which depends on . We also have for any that and
[TABLE]
where the constant here depends on and the product .
Proof.
In this proof we use the results and notation for the Faá di Bruno result presented in App. A.
For 6.10 using A.1, we see that
[TABLE]
Then
[TABLE]
We see that from assumption 6.3
[TABLE]
Then, we have that
[TABLE]
The inequality 6.10 then follows from integration and a Minkowski inequality.
[TABLE]
We bound each of and in turn. We see that for , applying 6.6
[TABLE]
[TABLE]
Combining the above estimates and integrating over the domain, we see
[TABLE]
Define . Then, a simple induction argument shows that
[TABLE]
where satisfies
[TABLE]
Given a surface finite element (Def. 6.1), let be the basis dual to . This is the set of basis functions of the finite element. If is a function for which all , is well defined, then we define the local interpolant by
[TABLE]
We can think of as the unique shape function that has the same nodal values as so that, in particular, for .
Theorem 6.13** (Local interpolation estimate).**
Let be a -surface finite element (Def. 6.4) with reference element which satisfies the assumptions of Lem. 4.3 for some , . Then there exists a constant such that for all functions ,
[TABLE]
Proof.
We re-scale 4.1 using Lem. 6.12 and the estimates from 6.9:
[TABLE]
The last line holds if (note that so this statement is true for small enough). ∎
6.2. Triangulated hypersurface and surface finite element spaces
We will next bring together several surface finite elements in order to define as a collection of finite element domains.
Definition 6.14*.*
- (a)
A triangulated hypersurface is a set equipped with an admissible subdivision consisting of surface finite element domains 6.2a such that , for with . 2. (b)
The maximum subdivision diameter is defined by:
[TABLE] 3. (c)
Let be a discrete hypersurface equipped with an admissible subdivision such that each set is an element domain for a surface finite element parametrised over the same polygonal reference finite element . We say that is a facet if is the image of a boundary facet of . 4. (d)
We say that is a conforming subdivision of if any facet of an element domain is either a facet of another element domain , in which case we say and are adjacent, or a portion of the boundary (if such a boundary exists). 5. (e)
For a conforming subdivision of , we denote by the set of facets between adjacent elements and any boundary facets. For a common facet between elements , we make a fixed choice that the conormal on will be denoted and the conormal on will be denoted . The choice of which element is on which side of the facet is not important in our considerations.
Remark 6.15*.*
- •
In the above definition we do not impose any global assumptions on the connectivity or smoothness of . Thus there may not be an underling smooth surface.
- •
We orient a discrete hypersurface which is equipped with a conforming subdivision by choosing a particular sign to the element-wise definition of normal. We restrict that the induced orientation of the intersection of adjacent element domains are opposite. For example, for a simplex reference element, the vertices in facets between two elements should be ordered oppositely in each element.
For any triangulated hypersurface , we may define spaces of Lebesgue integrable functions with the usual norms for .
Definition 6.16* (Broken Sobolev spaces and norms).*
Let be a subdivision of consisting of -surface finite elements. Then for , , we define the broken Sobolev space by
[TABLE]
with norm
[TABLE]
Lemma 6.17**.**
The space is complete.
Proof.
Consider a Cauchy sequence . This implies
- •
is Cauchy in so there exists such that in .
- •
is Cauchy for all so there exists such that in .
It is clear from the triangle inequality that :
[TABLE]
since the right hand-side converges to [math] as . Hence, we have shown that converges to a function in the norm. ∎
We also want more connectivity between elements. This is achieved in the following space.
Let be a triangulated hypersurface with conforming subdivision . For , we denote the trace of a function by and recall that there exists a constant such that
[TABLE]
We define the space by
[TABLE]
We equip this space with the broken norm .
Lemma 6.18**.**
The space is a closed subspace of so is complete.
Proof.
Take a sequence which converges to . Then for any pair of adjacent elements we have
[TABLE]
Clearly the right hand side converges to [math] as so we have the traces of from adjacent elements coincide and . ∎
We will use the notation for which is a Hilbert space when equipped with the obvious broken inner product.
6.2.1. Surface finite element space
We now restrict to Lagrangian finite elements over a polygonal reference finite element. Here we assume that the degrees of freedom for each element are given by
[TABLE]
where is a finite set of nodes in . We call the set of Lagrange nodes of . This restriction avoids difficulties in defining the edge of elements and how to effectively bring elements together to form a global finite element space. Extensions to other element types such as Hermite elements are left to future work.
Finally, the set of degrees of freedom of adjacent surface finite elements will be related as follows. Let and be two surface finite elements such that and are adjacent with and . Then, we have
[TABLE]
We denote the global set of Lagrange nodes by
[TABLE]
For each , let be the local neighbourhood of elements for which .
Definition 6.19* (Surface finite element space).*
- (a)
Let be a discrete hypersurface equipped with a conforming subdivision with each domain equipped with a surface finite element (Def. 6.1) which satisfy 6.20. A surface finite element space is a (generally proper) subset of the product space given by
[TABLE] 2. (b)
The surface finite element space is determined by the global degrees of freedom
[TABLE]
In this definition, an element is not, in general a “function” defined over , since we do not necessarily have a good definition of over element boundaries: The “function” may be double-valued.
If it happens, however, that for each element , the restrictions and coincide along the common face of any adjacent elements and , then the function can be identified with a function defined over the set . In this case, we call the elements surface finite element functions. Examples of surface finite element functions are shown in Fig. 6.2.
Lemma 6.20**.**
Let be a surface finite element space consisting of surface element elements over a conforming subdivision of . Assume further that for each , the corresponding reference finite element is a Lagrange element of order (Ex. 4.2). Then we can identify elements of as functions in . Furthermore is a closed subspace in .
Proof.
Consider two adjacent elements and a element . The functions and when restricted to the appropriate edges in are polynomials of degree which agree at the Lagrange points on this edge from 6.20 and the definition of . The Lagrange points in the reference element determine polynomials of degree so we have that on . Since is a conforming subdivision we can define a global function such that for each which is globally continuous. Indeed restricted to each element is continuous, as the composition of a polynomial (element of for some ) and a smooth surface finite element reference map , and is single valued on the facets where any two elements meet.
In fact the restriction to each element is a -function hence so it is clear that . The space is closed since it is finite dimensional. ∎
We enumerate the nodes so that and take to be the basis of dual to . Since, we have a finite basis of we note that we can identify any with a vector so that
[TABLE]
Definition 6.21* (Interpolation).*
If is a function on for which all , , is well defined (in case of Lagrangian finite elements, suffices), then we can define a global interpolant by
[TABLE]
Note that our construction implies that
[TABLE]
and for all .
In order to prove estimates on the global interpolant, we will first define three further properties of our subdivision .
Definition 6.22* (Regular and quasi-uniform subdivisions).*
For , let be a triangulated hypersurface (Def. 6.14(a)) equipped with a conforming subdivision (Def. 6.14(d)).
- (a)
The family is said to be non-degenerate or regular if there exists such that for all and all ,
[TABLE]
and there exists a constant such that
[TABLE] 2. (b)
The family is said to be -regular if it is regular, if for all and all , , and if, there exists a constant such that
[TABLE]
This implies that is a -reference surface finite element map for each . 3. (c)
A regular family is said to be quasi-uniform if there exists such that
[TABLE]
Remark 6.23*.*
We note that:
- •
for a regular subdivision there exists a constant depending on the global quantities and
[TABLE]
- •
for a quasi-uniform subdivision there exists a constant depending on the global quantities and
[TABLE]
Theorem 6.24** (Global interpolation estimates).**
For , let be a triangulated hypersurface (Def. 6.14(a)) equipped with a -regular (Def. 6.22(b)), quasi-uniform (Def. 6.22(c)), conforming (Def. 6.14(d)) subdivision . Let each be equipped with a -surface finite element (Def. 6.4) parametrised over a reference finite element which satisfies the assumptions of Thm. 6.13 for some , . Then there exists a constant such that for all functions ,
[TABLE]
Proof.
The proof follows by piecing together Thm. 6.13 using the fact that is quasi-uniform. ∎
6.3. Evolving surface finite elements
Let . We consider families of surface finite elements, spaces and triangulated hypersurfaces parametrised by .
Definition 6.25* (Evolving surface finite element).*
- (a)
Let be a time dependent family of surface finite elements (Def. 6.1) parametrised over a common reference element . If the constant is uniformly bounded away from 1,
[TABLE]
we say that is an evolving surface finite element. 2. (b)
Let where . We say that if such that
[TABLE]
then is the flow defining the evolution of the element domain and that is the evolving surface element reference map. 3. (c)
The element velocity of is defined by
[TABLE] 4. (d)
If each is a -surface finite element for each and the constants are uniformly bounded:
[TABLE]
then we say that is an evolving -surface finite element and is the evolving -surface element reference map. 5. (e)
We say that an evolving finite element domain is temporally quasi-uniform, if there exists such that
[TABLE] 6. (f)
The family of element push forward maps , indexed by , is defined to be the linear invertible map given for by where
[TABLE]
Lemma 6.26**.**
Let , for , be an evolving -surface finite element reference map (Def. 6.25(d)) for a temporally quasi-uniform element domain (Def. 6.25(e)) and the family of element push forward maps (Def. 6.25(f)). Then there exists constant which depend only on the reference element domain , and the constant and , such that for and all , if, and only if, and
[TABLE]
Proof.
From Lem. 6.9 and 6.12, we have
[TABLE]
and
[TABLE]
It can be easily seen that for a quasi-uniform evolving surface finite element that these constants only depend on allowed quantities. ∎
This result implies that is a compatible pair (Def. 2.2) and in particular is a compatible pair when equipped with the -norm (as is a closed subspace of , Rem. 2.3 and 6.5(c)).
6.4. Evolving surface finite element spaces
We now formulate an evolving surface finite element space forming part of a compatible pair (in the sense of Sec. 2, Def. 2.2). For each , we are given a family of discrete hypersurfaces and each equipped with a surface finite element space . Furthermore, we are interested in under what assumptions does the compatibility hold independently of the mesh size .
Definition 6.27* (Evolving triangulated hypersurface).*
For , let be a family of triangulated hypersurfaces (Def. 6.14(a)) each equipped with a conforming subdivision (Def. 6.14(d)) such that each element domain is equipped with an element flow map (Def. 6.25(b)).
- (a)
We call an evolving conforming subdivision if for each element and each facet of either is a facet of another element , in which case is a common facet between and for all or is a portion of the boundary , in which case is a portion of the boundary for all . 2. (b)
An evolving triangulated hypersurface is defined to be a family of triangulated hypersurfaces equipped with an evolving conforming subdivision. In this case, we define the mesh parameter to be
[TABLE] 3. (c)
We define a global discrete flow element-wise by
[TABLE]
Our assumptions imply that is piecewise smooth and . 4. (d)
We define a global discrete velocity given by
[TABLE] 5. (e)
The family of linear homeomorphisms induced by the flow and called the global push forward map is denoted by for and defined by
[TABLE]
Again, in order to bring together a collection of evolving surface finite elements we restrict to Lagrangian surface finite elements. For each , we denote by the global set of Lagrange nodes of 6.21 and for any , is the set of elements such that is a node of . We make the further restriction that the global flow is single-valued at each Lagrange point: for all we have
[TABLE]
Definition 6.28* (Evolving surface finite element space).*
- (a)
Let be an evolving triangulated hypersurface (Def. 6.27(b)) equipped with an evolving conforming subdivision (Def. 6.27(a)). For , let be a surface finite element space (Def. 6.19(a)) over . If each is equipped with an evolving surface finite element (Def. 6.25(a)) then we say is an evolving surface finite element space. 2. (b)
For each , we will write for the set of global nodal variables (Def. 6.19(b)). We will use the convention that
[TABLE]
where is the trajectory of a Lagrange point under the global flow . We will denote by the global basis of finite element functions such that for and all . This implies that . See also 3.4.
Definition 6.29* (Uniformly regular and uniformly quasi-uniform evolving subdivisions).*
For , let be a family of evolving conforming subdivisions (Def. 6.27(a)).
- (a)
We say that the family is uniformly regular if there exists such that for all and all times , we have
[TABLE]
and there exists such that
[TABLE] 2. (b)
We say that the family is uniformly -regular if it is uniformly regular, if for each time , the family is -regular and if there exists a constant such that
[TABLE]
This implies that is an evolving -reference finite element map for each . 3. (c)
We say that the family is uniformly quasi-uniform if there exists such that for all and all times , we have
[TABLE]
Note that a uniformly quasi-uniform subdivision consists of element domains for temporally quasi-uniform evolving surface finite element domains.
Lemma 6.30**.**
For , let be a uniformly -regular (Def. 6.29(b)), uniformly quasi-uniform (Def. 6.29(c)), evolving, conforming subdivision (Def. 6.27(a)) and let be the global push-forward map (Def. 6.27(e)). Let , . Then if and only if for all . Furthermore, there exists independent of and such that for all
[TABLE]
Remark 6.31*.*
This implies that is a compatible pair. In particular the pairs , equipped with the broken Sobolev norm (Def. 6.16), and are compatible since each is a closed subspace (Rem. 2.3, 6.20 and 6.18).
Proof.
We simply sum the element-wise result from Lem. 6.26. The constants are independent of and due to the uniform quasi-uniformity of . ∎
Note that this result implies that the spaces and are well defined when equipped with the appropriate norms (c.f. 2.1 and 2.2).
7. Lifted surface finite element spaces
So far we have only defined surface finite elements without relation to approximation of a smooth hypersurface and function spaces on through the definition of a lifted finite element space . In general, due to the curvature of the surface, the computational domain will be an approximation to . We will identify surface finite elements on with a corresponding curved surface finite elements on using a mapping . We call this process lifting. This provides a convenient way to compare the smooth functions (solutions to PDEs) on with the lift of discrete functions on . We will also require an inverse lift that maps functions on the smooth domain to the computational domain.
In the following, we will answer the question what assumptions on must be made so that is a -surface finite element and how we can put several such lifted surface finite elements together in order to recover and finite element spaces on . We will also explore some interpolation results when these assumptions hold. Finally, we explore how we can lift the discrete push forward map.
7.1. Lifted surface finite element
We consider the situation of a surface finite element (Def. 6.1), with element reference map , where the element domain approximates a portion of .
Let be a -map which is a diffeomorphism onto its image. For a function , we call defined by the lift of . We will assume that we can decompose into
[TABLE]
where is an invertible matrix, and . We will assume that does affect the affine part of the parametrisations:
[TABLE]
Definition 7.1* (Lifted surface finite element).*
We call the triple defined by
[TABLE]
the lift of and the lifting map. If forms a surface finite element over then we say that is the lifted surface finite element associated with . In this case we call the lifted surface finite element reference map.
The next two results show under what assumptions on is a surface finite element (Def. 6.1) or a -surface finite element (Def. 6.4).
Lemma 7.2**.**
If satisfies that
[TABLE]
then is a surface finite element. Furthermore, for , we have that there exists such that
[TABLE]
where the constants depend on , and the ratio .
Proof.
To show that is a surface finite element (Def. 6.1), the conditions on the element reference map are clear and we are left to check the curvedness condition 6.1. Using the expansion of , we see
[TABLE]
So that
[TABLE]
Applying the curvedness condition for and the fact that is a projection (Rem. 6.3(c) applied to the flat element ) we see
[TABLE]
The curvedness condition is shown by applying 7.1.
To show 7.2 the result is clear for . For , we will apply Lem. 6.12. Then we see
[TABLE]
where depends on and the bound 7.1. For the higher order bound, we note that
[TABLE]
so that we infer
[TABLE]
Then applying Lem. 6.12 once more we see
[TABLE]
The final result is given by applying 6.9a, 6.9b and 7.4c. ∎
Lemma 7.3**.**
If is a -surface finite element and satisfies and , then is also a -surface finite element. Furthermore, for and , we have that there exists such that
[TABLE]
where the constants depend on , and the ratio .
Proof.
To see that is a -surface finite element (Def. 6.4), the first three conditions are clear from the smoothness assumption on . We must show 6.3 holds: for , there exists a constant such that
[TABLE]
Computing directly using the Faá di Bruno formula A.1, we have
[TABLE]
But applying the fact that is a -surface finite element and the smoothness assumption on , we see that
[TABLE]
Finally, applying 7.3, we have that
[TABLE]
This shows that is a -surface finite element.
To show the norm equivalence we again apply Lem. 6.12, recognising the geometric progression, to see
[TABLE]
where
[TABLE]
The final result is given by applying 6.9a, 6.9b and 7.4c. ∎
Remark 7.4*.*
Considering the case that fixes the vertices of then we can write as
[TABLE]
That is that and . Then the assumptions of Lem. 7.2 can be replaced by
[TABLE]
and the assumptions of Lem. 7.3 can be replaced by the assumption that is uniformly bounded.
We next relate the geometry of the base and lifted element domains.
Lemma 7.5**.**
Using the decomposition of , we have
[TABLE]
Proof.
For 7.4a, we show the second inequality. The first follows by the same reasoning applied to the inverse of . Since is compact, there exists such that
[TABLE]
But since is invertible, there exists such that
[TABLE]
Then, we can compute that
[TABLE]
Similarly, for 7.4b, we only show the first inequality. For each , there exists a ball in such that
[TABLE]
Denote by and the centre and radius of respectively. Consider the affine map given by
[TABLE]
It is clear that maps balls to balls and is mapped to a ball centred at with radius . We claim is contained in . Indeed, take , then there exists such that . Denote by
[TABLE]
then
[TABLE]
so and
[TABLE]
so that .
Therefore, we have found a ball with radius , thus we can infer that
[TABLE]
Since the proof holds for all , we see the desired result.
Given 7.4a and 7.4b, the final result 7.4c follows directly from Rem. 6.10. ∎
We can also use to define an inverse lift. Given an element domain , lifted element domain and lifting map . We know that is invertible onto its image, namely . So for , we denote its inverse lift by defined by
[TABLE]
Lemma 7.6**.**
If the assumptions of Lem. 7.2 and 7.3 hold then, for and , we have that there exists such that for all , we have
[TABLE]
Proof.
The same proof can be applied to the case of the inverse lift as well as the usual lift. ∎
7.2. Lifted surface finite element space
Let be a smooth hypersurface and for and let be a triangulated hypersurface (Def. 6.14(a)) equipped with a conforming subdivision (Def. 6.14(d)) and a surface finite element space (Def. 6.19(a)). Let each be associated with a lifted finite element with lifted map .
Definition 7.7*.*
- (a)
We denote by the set of all lifted element domains
[TABLE]
If the global map is single valued and forms a conforming subdivision of , we say that is an exact subdivision of the surface . In this case we set to be the set of facets between adjacent elements. 2. (b)
We define a global lifting map by . We define the inverse lift in a similar element-wise fashion by . 3. (c)
We denote by the global lift given for by
[TABLE]
and by the global inverse lift given for by
[TABLE]
We will also use the notation and . 4. (d)
If for each surface finite element there is an associated lifted surface finite element , then, we define a lifted surface finite element space (c.f. 3.17) by
[TABLE]
Proposition 7.8**.**
Assume additionally that the family of triangulations is regular (Def. 6.22(a)), satisfies the assumptions of Lem. 7.2 and there exists such that
[TABLE]
for all for all . Then is regular and is a surface finite element space and there exists constants , which are independent of such that
[TABLE]
Furthermore, if satisfies the assumptions of Lem. 7.3 for all and all , then for , there exists independent of such that
[TABLE]
Proof.
The regularity of the family of triangulations follows from 7.4a and 7.4b from assumption 7.5.
The results follow by combining the previous results for each element in . We achieve bounds independently of since the regularity of the subdivisions implies that is independent of . ∎
Lemma 7.9**.**
Let be an exact decomposition of then .
Proof.
First, let . Then we have that and it is left to show that has a weak derivative in (Dziuk and Elliott, 2013b, Def. 2.11). We have a candidate given element-wise by . It is clear that and for with compact support and , we have
[TABLE]
We note that we can write
[TABLE]
where is the set of facets between adjacent elements in , since the traces of and to coincide for any adjacent pair . We note that the sum over edges is zero since the conormals from adjacent elements are equal and opposite in the exact triangulation and we see that is the weak derivative of .
Second, let then it is clear that , and by the trace theorem has common trace between adjacent elements. ∎
We can also show a special interpolation estimate which interpolates smooth functions over the continuous surface into the lifted surface finite element space. We denote by an interpolation operator defined by
[TABLE]
Theorem 7.10** (Global lifted interpolation theorem).**
For , let be a -regular (Def. 6.22(b)), quasi-uniform (Def. 6.22(c)) family of triangulations of equipped with a surface finite element space (Def. 6.19(a)) consisting of -surface finite elements (Def. 6.4) over a reference element which satisfies the assumptions of Lem. 4.3 for some , . Let each element be equipped with a lifting map such that and 7.5 and 7.1 hold each uniformly for , then is a -regular, quasi-uniform family of triangulations of and is a surface finite element space consisting of -surface finite elements. Let be a continuous function, then is well defined. Furthermore, if the assumptions of Thm. 6.13 hold for the reference element , there exists a constant such that for all functions ,
[TABLE]
Proof.
The result follows by combining previous lemmas in the appropriate way. We see the lifted triangulation is quasi-uniform by applying the results in Lem. 7.5 and the assumption 7.5. The interpolation result then follows by applying Thm. 6.24. ∎
Thm. 7.10 will be used to show the abstract approximation properties I1 and I2.
Corollary 7.11**.**
Let , then the interpolant of into , denoted by , is given by
[TABLE]
Furthermore, there exists a constant such that for all functions , and all , we have
[TABLE]
Proof.
We apply the inverse lift result to the estimates in the theorem. ∎
7.3. Evolving lifted finite element spaces
For , let be a smoothly evolving hypersurface with flow map : i.e. and . For , let be an evolving triangulated hypersurface (Def. 6.27(b)) with global discrete flow (Def. 6.27(c)) and equipped with an evolving conforming subdivision (Def. 6.27(a)) and an evolving surface finite element space (Def. 6.28(a)). We assume that we are given a global lifting map (Def. 7.7(b)) which gives an exact subdivision of (Def. 7.7(a)).
Definition 7.12* (Lifted discrete flow map, material velocity and pushed forward map).*
- (a)
The lifted flow map of the smooth hypersurface is defined by
[TABLE]
We note that . 2. (b)
The lifted discrete material velocity on is defined by
[TABLE] 3. (c)
The family of lifted push forward maps is denoted by for and is defined by
[TABLE]
See also 3.13.
Remark 7.13*.*
Note that in general is different to , but each describes a different parametrisation of the same evolving surface. Also , the material velocity of , and define the same surface evolving from , so have the same normal components whilst the tangential components may not agree.
Proposition 7.14**.**
For , let be an evolving surface finite element space over a -regular (Def. 6.22(b)), uniformly quasi-uniform (Def. 6.29(c)), evolving conforming subdivision consisting of -evolving surface finite elements (Def. 6.25(d)) and the global push-forward map (Def. 6.27(e)). For each and each , let each element be equipped with a lifting map such that and 7.5 and 7.1 hold each uniformly for and . Then is also an evolving surface finite element space over a -regular, uniformly quasi-uniform, evolving conforming subdivision consisting of -evolving surface finite elements. Furthermore, for , , there exists independent of and such that
[TABLE]
In particular, the pair is compatible (Def. 2.2). Furthermore the pairs , equipped with the norm , and are also compatible.
Proof.
The initial properties of follow for the same reasoning as Prop. 7.8 since the constants are bounded uniformly in time. The bounds 7.12 and the compatibility of follow since the assumptions imply that are compatible and we can simply lift this result with uniform bounds. The compatibility of and follow since they are closed subspaces of (Lem. 6.20, 6.18 and 7.9) and Rem. 2.3. ∎
Part III Application to parabolic equations on evolving domains
8. Application I: Parabolic equation on an evolving bulk domain
In this section, we will formulate and analyse a finite element method for a parabolic problem posed in an evolving bulk domain 1.4. We will begin with notation and definitions for the initial value problem. The numerical method is based on the general theory of Sec. 4 applied here with isoparametric bulk finite element spaces of order . The analysis is based on the abstract theory presented in Sec. 2 and 3.
8.1. The domain and function spaces
We set , and . We will also make use of the spaces and . We see that for each and the inclusions are uniformly continuous. We define the push-forward operator by
[TABLE]
and are compatible pairs (Def. 2.2), the spaces , and are well defined (c.f. 2.1) and we have a well defined strong material derivative, denoted 2.4. form an evolving Hilbert triple (Def. 2.4). We also have that and are compatible pairs. For details see Alphonse et al. (2015b).
Remark 8.1*.*
For the well posedness of the partial differential equation Prob. 8.3 we require that the boundary is a -hypersurface and that the flow map . See Alphonse et al. (2015b) for more details. For the approximation properties we derive we require that is a -hypersurface and that with and both of class .
We introduce a signed distance function for the boundary surface . The oriented signed distance function for is given by
[TABLE]
For each , we orient by choosing the unit normal for . Our assumptions on imply that there exists a neighbourhood of and normal projection operator given as the unique solution of
[TABLE]
See Gilbarg and Trudinger (1983, Lem. 14.16); Foote (1984) for more details.
The final result we show in this section will be useful in the error analysis of our methods.
Lemma 8.2** (Narrow band trace inequality).**
For , let be the band given by
[TABLE]
Then there exists a constant such that for all ,
[TABLE]
Proof.
The proof for stationary domains in given by Elliott and Ranner (2013, Lem. 4.10) which can be easily extended to the evolving case. ∎
8.2. The initial value problem
We assume that is an symmetric diffusion tensor, a smooth vector field and is a smooth scalar field. We assume that for each , is uniformly positive definite: There exists such that for all
[TABLE]
We consider the initial value problem
Problem 8.3**.**
Find such that
[TABLE]
Remark 8.4*.*
The problem (1.4) is recovered by setting and .
8.2.1. Transport formulae
The following transport formulae hold on portions
| of the domain , which follow the flow for , and have Lipschitz boundaries at each time. |
|---|
- •
For by
[TABLE]
- •
For , we have the identity
[TABLE]
where is given by
[TABLE]
and is the rate of deformation tensor
[TABLE]
- •
For , , we have
[TABLE]
where is given by
[TABLE]
8.3. The bilinear forms and transport formulae
We define
[TABLE]
We can apply 8.5, 8.6 and 8.8 to see that we have the transport laws
[TABLE]
with the new forms
[TABLE]
where
[TABLE]
and is the rate of deformation tensor:
[TABLE]
We define to be
[TABLE]
The smoothness assumptions on the coefficients and the velocity imply that , and are uniformly bounded. There exists a constant such that
[TABLE]
8.4. Variational formulation
We consider the following variational formulation of 1.4:
Problem 8.5**.**
Given \emph{\mathrm{u}}_{0}\in L^{2}(\Omega_{0}), find \emph{\mathrm{u}}\in\mathcal{W}(\mathcal{V},\mathcal{V}^{*})such that for almost every we have
[TABLE]
Theorem 8.6**.**
There exists a unique solution to Prob. 8.5 which satisfies the stability bound:
[TABLE]
and if
[TABLE]
Proof.
We simply apply the abstract theory of Thm. 2.15. For Ass. 2.6 and 2.9 we refer to Alphonse et al. (2015b, Sec. 4 and Sec. 5). Ass. 2.13 is a consequence of Rem. 2.14. Also Assumptions M1, M2, G1, G2, A1, A2, A3, As4, Bs1, Bn1, Bs2 and Bn2 hold. It is clear that M1 and M2 hold since is equal to the -inner product. Similarly assumptions G1 and G2 follow from the transport formula 8.9 and the boundedness of the velocity. We know that the map is differentiable hence measurable which shows A1. The coercivity A3 and boundedness As4 of and the boundedness of An4 follow from standard arguments. The existence of the bilinear forms Bs1 and Bn1 has been shown in 8.10 and 8.11 and the estimates Bs2 and Bn2 are shown in 8.13 and 8.14. ∎
8.5. Discretisation of the domain and finite element spaces
The first stage in constructing our finite element method is to define an approximate computation domain . Our construction satisfies that the boundary Lagrange points of lie on the boundary of and all Lagrange points evolve with the prescribed velocity . We will consider as an interpolant of . We recall is the order of isoparametric bulk finite element spaces we wish to use. Throughout the remainder of this section we will denote global discrete quantities with a subscript , which is related to element size. We assume implicitly that these structures exist for each in this range (see also Rem. 4.11).
We will use the simplical, Lagrangian reference element over order from Ex. 4.2. We start by constructing a family of time dependent element reference maps (Def. 4.30(b)) which will define an evolving conforming subdivision (Def. 4.32(a)) of an evolving triangulated open domain (Def. 4.32(b)).
Let be a polyhedral approximation of equipped with a quasi-uniform, conforming subdivision into simplicies (see Sec. 4.3 for details). We denote by . We restrict that the vertices of lie on the surface . We assume that the normal projection operator 8.2, is a homomorphism from onto . More precisely, for each , there exists an affine map which satisfies the assumption of Def. 4.4 so that we can define a bulk finite element using 4.3. We assume the Lagrange points satisfy 4.19.
We extend to construct a bijection which we will define element-wise. A similar construction is used by Elliott and Ranner (2013). We first decompose into boundary elements, which have more than one vertex on the boundary, and interior elements. For an interior element , we define
[TABLE]
Otherwise, let be a boundary element and consider . Denote by the vertices of ordered so that lie on (recall that ). First, decompose into barycentric coordinates:
[TABLE]
We introduce the function and the singular set by
[TABLE]
The scalar represents how far we are from with on . Note that we have . The set is the set of points in furthest from . If , we denote the projection onto by . We see that is given by
[TABLE]
Then, we define by
[TABLE]
Remark 8.7*.*
This mapping takes points on onto , since here and and we recover . Points in remain unchanged by this mapping, since here and furthermore the power on ensures that the mapping is on the closed domain .
We write for interpolation over 4.12 and define an initial element reference map by
[TABLE]
We call the union of all domains constructed in this way , which is a conforming subdivision (Def. 4.14(d)), and call the union of element domains , which is a triangulated bulk domain (Def. 4.14(a)). Finally, we call the Lagrange nodes of which satisfy 4.19 since Lagrange points on the boundaries of each element are not moved by . An example of domains constructed using the above are shown in Fig. 8.1.
To complete the construction, for , we consider the discrete domain given by the discrete flow defined by
[TABLE]
which is a bijection onto its image for small enough. We denote its inverse by . Using 4.22 defines an evolving reference element map and with 4.3 we have a bulk finite element (Def. 4.30(a)). We call the set of such domains at each time and we define by
[TABLE]
and write a global discrete flow map defined element-wise by (Def. 4.32(c)). By we denote the maximum mesh diameter over time 4.24:
[TABLE]
We introduce the Hilbert spaces and (c.f. Lem. 4.20). Since is a globally continuous function, it is clear that the Lagrange points in each element satisfy 4.25. We define a global discrete function space (Def. 4.34(a)) by
[TABLE]
Using Lem. 4.23 we can identify elements in as continuous functions and .
We will assume that this construction results in a uniformly quasi-uniform evolving subdivision (Def. 4.35(c)) - that is that the velocity is such that the simplicies in do not become too distorted. In order to show the result of this construction satisfies the other assumptions we require we first state a result shown by Elliott and Ranner (2013):
Lemma 8.8** (Prop. 4.4, Elliott and Ranner 2013).**
The mapping is of class when restricted to each element and the composed map satisfies
[TABLE]
Proposition 8.9**.**
The space defined by the above construction is an evolving bulk finite element space consisting of -surface finite elements over a uniformly -regular evolving subdivision.
Proof.
Let be a single element in . We can write the element parametrisation as
[TABLE]
where is the flow map from defined as the piecewise linear interpolant of . We note that is linear so we define the splitting
[TABLE]
Since and agree at the vertices of , we have that is a linear interpolant of over . Thus, we infer that
[TABLE]
Here, we have used the notation for piecewise linear interpolation over and applied the Bramble-Hilbert Lemma (Lem. 4.3) and the rescaling 4.10.
Hence we have
[TABLE]
Since we have assumed that is uniformly quasi-uniform we see that is uniformly bounded, it remains to show that is uniformly bounded. However, this follows from the definition of and the smoothness of 8.19 and the smooth flow map . Therefore, for sufficiently small, .
To show the higher order bounds, we compute similarly that for
[TABLE]
Hence, we see that
[TABLE]
In the final inequality, we use that . The bounds on then follow since we have assumed that the mesh is uniformly quasi-uniform and that is sufficiently smooth. ∎
The element flow map defines a velocity on each element (Def. 4.30(c)) by
[TABLE]
This can be combined into a global velocity (4.32(d)). We note that the global velocity is determined purely by the velocity of the vertices :
[TABLE]
We also have a discrete push forward map (Def. 4.32(e)) for given by
[TABLE]
Since we have constructed a uniformly -regular mesh, we infer that is a compatible pair (Lem. 4.36) so we may define the space 2.2 and a discrete material derivative by 3.3:
[TABLE]
We note that the pairs and are also compatible (Lem. 4.36 and 4.37) so we may also define the spaces and .
Lemma 8.10**.**
For and for each , let be a smooth, positive-definite, diffusion tensor on and be a smooth vector field on for each and . For we have the transport formula:
[TABLE]
Furthermore for we have
[TABLE]
and for
[TABLE]
where and are given by
[TABLE]
and is the rate of deformation tensor
[TABLE]
Proof.
We note that the left hand side maybe decomposed into individual elements then apply 8.5, 8.6 and 8.8 on each element.∎
8.6. Construction of the lifted finite element space
We construct a bijection between the computation domain and the continuous problem domain which we will call the lifting operator. We do this using a similar construction to used to define at the start of Sec. 8.5. It will again be based on using an extension of the normal projection operator used as a lifting operator in Sec. 9.6.
Fix . We wish to construct a bijection which we will define element-wise. We decompose into boundary elements, which have more than one vertex on the boundary, and interior elements. For an interior element , we define
[TABLE]
Otherwise, let be a boundary element and consider . Denote by the vertices of ordered so that lie on . We recall that the element is given by a parametrisation over a reference element so that, fixing a point , we can define points and vertices in by
[TABLE]
We decompose into barycentric coordinates on :
[TABLE]
We introduce the function and the singular set by
[TABLE]
Again, we note that , and for . If , we denote the projection onto by given by
[TABLE]
We then define by
[TABLE]
The lifting map induces a lifted subdivision given as the set of all elements
[TABLE]
Remark 8.11*.*
The ideas behind this construction are similar to the construction of . See Rem. 8.7 for more discussion and also Lem. 8.12.
We next follow a sequence of calculations to show the properties of . These estimates are based on previous work by Bernardi (1989) and Elliott and Ranner (2013). It is useful to recall the following Faá di Bruno formula (Bernardi, 1989, Eq. 2.9) for two smooth functions
[TABLE]
where are constants and is the set given by
[TABLE]
We wish to calculate derivatives of of order for . A direct calculation shows that
[TABLE]
for a constant independent of and . Then, we have that
[TABLE]
where in the second line, we have used that is a -surface finite element so
[TABLE]
Next, applying 8.26 and 8.27, we have
[TABLE]
Using a similar geometric construction to Lem. 9.10, we have
[TABLE]
Hence, we infer that
[TABLE]
Finally, using the Leibniz formula and the properties of , we have
[TABLE]
Lemma 8.12**.**
The lifting function is an element-wise -diffeomorphism and satisfies
[TABLE]
Proof.
The first result is clear for all internal elements. Consider a time and a single boundary element . The smoothness away from follows from the fact that restricted to each element is the composition of smooth functions. The above calculations, combined with 4.11 and 4.9b, show that for ,
[TABLE]
The mapping is of class on with derivatives of order less than or equal to tending to zero when tends to a point in . Hence, it can be extended to a mapping on . ∎
Lemma 8.13**.**
Furthermore, we have that the lifted triangulation is a uniformly -regular evolving subdivision (Def. 4.35(b)).
Proof.
We use the splitting
[TABLE]
and write . We start by taking in 8.29 to give
[TABLE]
Clearly the right hand side is less than for small enough for a uniformly regular subdivision. This inequality can clearly be translated to the required estimate on . ∎
For and a function , we define its lift by
[TABLE]
We will also make use of an inverse lift for functions on . For , we define the inverse lift of , denoted by by
[TABLE]
Lemma 8.14**.**
Let and denote its lift by . Then there exists constants such that
[TABLE]
Proof.
We apply Prop. 5.8 and 8.13. ∎
For each , we use Def. 5.1 to construct an associated lifted bulk finite element . We assume that the domains are such that the set of lifted element domains defines an exact decomposition of (Def. 5.7(a)) for each . For , we define the space of lifted functions to be given by (c.f. 3.17)
[TABLE]
Using the lift and the discrete flow , we can define the lifted flow map (Def. 5.11(a)), lifted discrete velocity (Def. 5.11(b)) and lifted push forward maps (Def. 5.11(c)). The results of Lem. 8.13 imply with Prop. 5.13 that and and form compatible pairs (Def. 2.2). We will use the notations and for the spaces of functions smoothly evolving in time 2.2 with respect to the push-forward map in and respectively. We may define a lifted strong material derivative 3.15:
[TABLE]
Lemma 8.15**.**
The flow map induces a new transport formula on . For we have
[TABLE]
Let be a smooth, positive-definite, diffusion tensor on and be a smooth vector field on for all . Furthermore, we have for
[TABLE]
and for we have
[TABLE]
where and are defined as in Lem. 8.10.
Proof.
We note that the left hand side may be decomposed into individual elements then apply 8.5, 8.6 and 8.8 on each element.∎
For the later analysis, we will also require bounds on the time derivative of . We consider an element and the trajectory of a point which follows the velocity field . From the definition of , we have that
[TABLE]
In particular, the barycentric coordinate representation of does not depend on time. Therefore, writing and using from 8.24, we have
[TABLE]
Then we can compute that if is a boundary element, recalling the definition of 8.25, we have
[TABLE]
A similar calculation to those prior to Lem. 8.12 show that for
[TABLE]
This follows by using the smoothness of the surface along with the fact that is an interpolant of (Cor. 5.10).
We conclude this section by showing some geometric estimates arising from the use of the lifting function .
Lemma 8.16**.**
Under the above assumptions, we have the estimates that
[TABLE]
Additionally, writing J_{h}(\cdot,t)=\sqrt{\det\big{(}(\nabla\Lambda_{h}(\cdot,t))^{t}(\nabla\Lambda_{h}(\cdot,t))\big{)}}, we have
[TABLE]
Proof.
We have already shown 8.36 as part of the proof of Lem. 8.12. Inequality 8.38 then follows by a result of (Ipsen and Rehman, 2008, Cor. 2.11).
To show the time derivative bounds, we note that for each we can write
[TABLE]
Then we can compute, for
[TABLE]
Hence we have that (up to a set of measure zero)
[TABLE]
Then applying 8.35, 8.36 together with Lem. 8.18, we have
[TABLE]
This shows 8.37.
For 8.39 we have, applying 8.36 and 8.37, that
[TABLE]
8.7. The discrete problem and stability
The practical finite element method is based on the variational formulation 2.12 of Prob. 8.5. We introduce a element-wise smooth -diffusion tensor , an element-wise smooth dimensional vector field and an element-wise smooth scalar field . We will use the notation and , for all , which we assume satisfy:
[TABLE]
We assume that is uniformly positive definite: There exists such that for all and
[TABLE]
Problem 8.17**.**
Given , find such that for every
[TABLE]
where we have
[TABLE]
To show the properties of these bilinear forms we require one further lemma:
Lemma 8.18**.**
The discrete velocity of the discrete evolving domain is uniformly bounded in . That is, there exists a constant such that for all
[TABLE]
Proof.
The bound follows using the characterisation 8.20 by using the interpolation bound shown in Cor. 5.10. ∎
We have a transport formula for the domain .
Lemma 8.19**.**
There exists bilinear forms , and such that
[TABLE]
where
[TABLE]
and
[TABLE]
Furthermore, there exists a constants such that for all and all we have
[TABLE]
Proof.
The bilinear forms exist due to the more general Lem. 8.10. The estimates follow from Lem. 8.18. ∎
Theorem 8.20**.**
There exists a unique solution of the finite element scheme 8.40. The solution satisfies the stability bound:
[TABLE]
Proof.
We apply the abstract result of Thm. 3.3. It is left to check the required assumptions. The assumptions on , Mh1 and Mh2, follow directly since is equal to the inner-product. That is measurable Ah1 and the estimates on , Ah2 and Ah3 follow in the same manner as Thm. 8.6. The transport formulae and estimates for and , Gh1, Gh2 Bh1 and Bh2, are shown in Lem. 8.19. ∎
8.8. Error analysis
The space is equipped with the following approximation property.
Lemma 8.21**.**
For there exists a Lagrangian interpolation operator that is well defined. Furthermore, the following bounds hold for constants independent of and time:
[TABLE]
Proof.
We apply Thm. 5.9. The second result applies the theorem in the obvious way. The first result applies the theorem with noting that and the inclusions for Lem. 4.3 still hold. ∎
We can also use the lift to define an evolving lifted triangulation. For each and , we define
[TABLE]
The edges of these curvilinear-simplicies evolve with a velocity which can be characterised as follows. Let be the trajectory of a point on according to the flow . Then we have that
[TABLE]
Now consider a point . The trajectory of defines the velocity field by
[TABLE]
Equivalently this means the flow is given by
[TABLE]
Lemma 8.22**.**
The pairs and are compatible and we may define a material derivative for and transport formula: There exists a bilinear for given by
[TABLE]
such that
[TABLE]
and there exists a constant such that for all and we have
[TABLE]
Furthermore, we have a new transport formula for the bilinear form. There exists a bilinear form given by
[TABLE]
such that
[TABLE]
and there exists a constant such that for all and we have
[TABLE]
Proof.
We simply apply Lem. 8.15. ∎
The bounds in Lem. 8.16 allow to show some of the abstract error bounds required use the result of Thm. 3.11.
Lemma 8.23**.**
We have the estimate:
[TABLE]
If then and if then . Moreover
[TABLE]
Proof.
We write for that
[TABLE]
This allows us to apply the interpolation theorem (Cor. 5.10) to the embedded velocity with the estimate 8.35 to achieve the estimate 8.51. The inclusions and bounds 8.52 follow from a simple calculation finding for appropriate and 8.51. ∎
For the remainder of this section, we will take , and . and assume that is of class in space.
Lemma 8.24**.**
There exists a constant such that for all and all the following holds for all with lifts :
[TABLE]
For with inverse lifts , we have
[TABLE]
Proof.
We use the notation J_{h}=\sqrt{\det\big{(}(\nabla\Lambda_{h})^{t}(\nabla\Lambda_{h})\big{)}}. We note that is the identity on elements away from the boundary and hence on these elements. We denote by the union of boundary elements where and note that we have
[TABLE]
For 8.53, we have
[TABLE]
Hence, we have, applying 8.38 and Lem. 8.14, that
[TABLE]
Applying the narrow band trace inequality Lem. 8.2 this can be improved to
[TABLE]
Similarly we have
[TABLE]
Applying 8.36, 8.38 and Lem. 8.14 we can see 8.55. Again, by applying Lem. 8.2 we show the improved bound in 8.59.
We apply a similar process to the proof of (Ranner, 2013, Lemma 3.3.14) combined with the results of Lem. 8.16 and the narrow band trace inequality (Lem. 8.2) to show the estimates 8.54, 8.56 and 8.60.
Finally, 8.57 and 8.58 follow from the estimate 8.51. The bound 8.61 follows from 8.59, the fact that and the estimate 8.52. Indeed we can compute that
[TABLE]
The remaining assumtion to verify is the estimate B3.
Lemma 8.25**.**
For any , let be as in 3.20 and then for all , we have
[TABLE]
Proof.
The proof is based on a duality argument from Douglas and Dupont (1973). Fix . We use the Hilbert triple (at each time ) and identify
[TABLE]
For any , let be the solution of the variational problem
[TABLE]
By the trace theorem the trace of lies in so that interpreting the duality product on the right-hand side as an inner-product 8.63, we see that
[TABLE]
This is a simple consequence of the trace theorem and the Riesz representation theorem.
As in Sec. 3.3, we introduce such that there exists a constant such that
[TABLE]
We wish to estimate in the -norm. We consider the dual problem: Given , find such that
[TABLE]
The problem is a weak form of an elliptic problem with inhomogenous Neumann boundary data and has a unique solution which satisfies the regularity estimate (Ladyzhenskaya and Uraltseva, 1968; Gilbarg and Trudinger, 1983)
[TABLE]
for a constant independent of time .
We see using 8.64b and in 8.65, that
[TABLE]
For the first term on the right hand side of 8.67, we apply the boundedness of and the interpolation estimate 8.45 to see
[TABLE]
For the second term on the right hand side, we apply the geometric estimates 8.55, 8.53 and 8.59 and the interpolation estimate 8.45 to see that
[TABLE]
Hence, combining the previous estimates and applying the dual regularity 8.66, we infer that
[TABLE]
and we have shown that
[TABLE]
Returning to B3, we now see that for and that
[TABLE]
For , using the smoothness of and the divergence theorem, we have
[TABLE]
We interpret the first term, , on the right hand side as the duality pairing between and its dual 8.63 and apply 8.68 and the trace theorem to see
[TABLE]
For , we use the smoothness of and to see
[TABLE]
Similarly for and we have
[TABLE]
Finally, we have collected all the estimates we require to show the error bound.
Theorem 8.26**.**
Let , , and and let be the solution of 8.15 which we assume satisfies
[TABLE]
Let be the solution of the finite element scheme 8.40 and denote its lift by . Then we have the following error estimate
[TABLE]
Proof.
The proof is performed by applying the abstract result from Thm. 3.11. We know the lift is stable from Lem. 8.14. The existence and boundedness of and are dealt with in Lem. 8.22. The interpolation properties I1 and I2 are shown in Lem. 8.21. The geometric perturbation estimates P1, P2, P3, P4, P5, P6, P4’, P5’, P8 and P9 are shown in Lem. 8.24 and 8.23 and P7 follow from P1, P4 and P4’ and that fact that lifting and taking material derivatives commute (Lem. 3.5). We have shown B3 in Lem. 8.25. ∎
9. Application II: Parabolic equation on an evolving surface
In this section, we will formulate and analyse a finite element method for a parabolic equation posed on an evolving surface 1.5. We begin with some notation and a definition of the initial value problem. Our numerical approach will be to first discretise the domain and construct isoparametric surface finite element spaces based on the general theory in Sec. 6 and 7. We will consider isoparametric finite elements of order , which will be fixed throughout this section. We will analyse this method using both techniques from the general surface finite element constructions in Sec. 6 and the abstract theory from Sec. 3.
9.1. The domain and function spaces
Let and be a compact sufficiently smooth hypersurface without boundary. Let be a family of evolving hypersurfaces such that there exists a sufficiently smooth mapping (called the flow map) such that:
- (1)
is a diffeomorphism of onto for each ; 2. (2)
.
We call
[TABLE]
the material velocity field of which also satisfies
[TABLE]
We denote , and . We will also make use of the spaces and . We see that for each and the inclusions are uniformly continuous. We define the push-forward operator by
[TABLE]
and are compatible pairs (Def. 2.2), the spaces , and are well defined (c.f. 2.1) and we have a well defined strong material derivative denoted by 2.4 (see Vierling 2014, Lem. 3.2,3.3). form a Hilbert triple (Def. 2.4). We also have that and are compatible pairs. For details see Alphonse et al. (2015b).
Remark 9.1*.*
For the well posedness of the partial differential equation Prob. 9.2 we require that the boundary is a -hypersurface and that the flow map . See Alphonse et al. (2015b) for more details. For the approximation properties we derive we require that is a -hypersurface and that with and both of class .
We introduce a signed distance function for a closed surface . We assume that is the boundary of an open bounded domain. The oriented distance function for is defined by
[TABLE]
We can orient by choosing the unit normal as
[TABLE]
This allows us to define the (extended) Weingarten map by and the mean curvature by . Note that, for example when , this definition of is the sum of the principal curvatures rather than the mean. For each time , there exists a narrow band such that the distance function is smooth and, for each there exists a unique point such that (see Lemma 14.16 of Gilbarg and Trudinger (1983) and Foote (1984))
[TABLE]
We call the operator the normal projection operator and note that is also smooth. We use this projection to extend the unit normal and Weingarten map to be defined in by and . See Gilbarg and Trudinger (1983, Lem. 14.16); Foote (1984) for more details.
9.2. The initial value problem
We assume that is a symmetric diffusion tensor which maps the tangent space of at a point into itself, and there exists such that for all
[TABLE]
is a smooth tangential vector field and is a smooth scalar field.
We consider the initial value problem
Problem 9.2**.**
Find such that
[TABLE]
Remark 9.3*.*
Writing for a decomposition of into tangential and normal components, the problem 1.5 is recovered by setting and .
9.2.1. Transport formulae
The following transport formulae hold on and on portions of the domain ,
| which follow the flow for , and have Lipschitz boundaries at each time |
|---|
- •
For
[TABLE]
where
[TABLE]
More precisely, we can show that Ass. 2.6 holds for (Alphonse et al., 2015b, Sec. 4.1).
- •
For , we have the identity
[TABLE]
where is given by
[TABLE]
and is the rate of deformation tensor
[TABLE]
- •
For , , we have
[TABLE]
where is given by
[TABLE]
Remark 9.4*.*
The identity 9.9 is equivalent to (Elliott and Venkataraman, 2015, Lem. A.1). The proof of 9.7 and 9.9 follows from applying 9.6 with the identity (in which is the unit normal to )
[TABLE]
9.3. The bilinear forms and transport formulae
We define
[TABLE]
We can apply 9.6 9.7 and 9.9 to see that:
[TABLE]
with
[TABLE]
We define to be
[TABLE]
9.4. Variational formulation
The weak formulation of Prob. 9.2 becomes
Problem 9.5**.**
Given \emph{\mathrm{u}}_{0}\in L^{2}(\Gamma_{0}), find such that for almost every we have
[TABLE]
We have the following well-posedness result.
Theorem 9.6**.**
There exists a unique solution to 9.14 which satisfies
[TABLE]
and if then
[TABLE]
Proof.
We simply check the assumptions required for Thm. 2.15. For Ass. 2.6 and 2.9 we refer to Alphonse et al. (2015b, Sec. 4 and Sec. 5). Ass. 2.13 is a consequence of Rem. 2.14. The assumptions M1 and M2 follow simply since is the -inner product. G1 holds from 9.6 and 2.6 from the assumption that is uniformly bounded in space and time. The bilinear form is differentiable in time, hence measurable A1. The coercivity of A3 follows from a standard calculation since is positive definite and is bounded.. The smoothness of and imply that and are bounded As4, An4. The existence of the bilinear forms Bs1 and Bn1 follow from 9.12 and 9.13 and the bounds Bs2 and Bn2 from the smoothness of and . ∎
9.5. Discretisation of the domain and finite element spaces
The first stage in developing our finite element method is to define the approximate computational domain . We do this by constructing an isoparametric approximation of which is then pushed-forwards under an approximation to the flow . The result will be that the Langrange points of lie on the surface for all times and evolving according to the velocity . In this sense, can be considered as an interpolation of . Recall that denotes the normal projection operator 9.4 and that is the order of isoparametric finite elements we will use. Throughout the remainder of this section we will denote global discrete quantities with a subscript , which is related to element size. We assume implicitly that these structures exist for each in this range (see also Rem. 6.11).
We will use the simplical, Lagrange reference element from Ex. 4.2. We start by constructing a family of time dependent element reference maps (Def. 6.1) which defines an evolving conforming subdivision (Def. 6.27(a)) of an evolving triangulated hypersurface (Def. 6.27(b)).
Let be a polyhedral approximation of equipped with a quasi-uniform, conforming subdivision (Def. 6.14(d)). We restrict the vertices of to lie on the surface and denote by the maximum element domain diameter on (Def. 6.14(b)). We assume that is sufficiently small and is such that is a smooth bijection from onto . More precisely, for each , there exists an affine map which satisfies the assumptions of Def. 6.1 so that we can define a surface finite element using 6.2. Note that the vertices of lie on but the other Lagrange points may not. We assume the collection of all Lagrange points satisfy 6.20. We write for the local interpolation operator over 6.13 and define an initial element reference map by
[TABLE]
An example of element domains 6.2a defined by is given for in Fig. 9.1. We denote the union of all elements constructed in this way which is a triangulated hypersurface (Def. 6.14(a)) equipped with a conforming subdivision (Def. 6.14(d)), the set of all element domains .
To complete the construction of , for , we consider the element flow map (Def. 6.25(b)) defined by
[TABLE]
which is a bijection onto its image and we denote its inverse by . An example is shown in Fig. 9.2. Using 6.23, defines an evolving reference element map:
[TABLE]
Using 6.2, this defines an evolving surface finite element (Def. 6.25(a)).
We call the set of such elements and the union and write a global discrete flow map given element-wise by (Def. 6.27(c)). By we denote the maximum mesh diameter over time 6.25:
[TABLE]
Our construction implies that 6.26 holds since each of the transformation from are continuous.
We will use two Hilbert spaces defined over . First we will denote by and by 6.19. We equip each space with a norm:
[TABLE]
We note that each element reference map is an element of so that at each time the triple is an isoparametric surface finite element. Furthermore, we can use the basis functions of to decompose :
[TABLE]
In particular, this implies that the evolving triangulated surface only depends on the evolving of the Lagrange points which we can infer satisfy
[TABLE]
We define a global evolving surface finite element space (Def. 6.28(a)) by
[TABLE]
Using Lem. 6.20 we can identify elements in as continuous functions and .
Remark 9.7*.*
This construction is a generalisation of the construction of Dziuk and Elliott (2007). Indeed, in the case that we wish to consider affine finite elements, it is worth noting that for and . A different view of the same construction is given by Kovács (2018) and at each time our construction coincides with the work of Demlow (2009).
We will assume that the evolving subdivision is uniformly quasi-uniform (Def. 6.29(c)). It is clear that our construction maintains the conformity of the initial base triangulation .
Proposition 9.8**.**
The evolving surface finite element space defined by the above construction consists of evolving -surface finite elements (Def. 6.25(d)) over a uniformly -regular subdivision (Def. 6.29(b)).
Proof.
The proof follows in the same way as Prop. 8.9 and we do not give full details here. The only part to check is that the discrete flow map is uniformly bounded in . However, this follows directly from the definition of as an interpolation of which is a smooth function. ∎
The element flow map defines a velocity on each element (Def. 6.25(c)) by
[TABLE]
This can be combined into a global velocity (Def. 6.27(d)). We note that the global velocity is determined purely by the velocity of the vertices :
[TABLE]
We also have a discrete push forward map (Def. 6.27(e)) on by
[TABLE]
Since we have constructed a uniformly -regular mesh, we infer that is a compatible pair (Lem. 6.30) and we may define the space 2.2 and a discrete material derivative by 3.3:
[TABLE]
We also have that and are compatible pairs (Rem. 6.31) so we may define the spaces and .
Lemma 9.9**.**
For and for each , let be a smooth, positive-definite, diffusion tensor on , which maps the tangent space of to itself, and be a smooth tangential vector field on for each and . Then for :
[TABLE]
For all , we have
[TABLE]
and for all , we have
[TABLE]
where and are given by
[TABLE]
and is the rate of deformation tensor
[TABLE]
Proof.
We note that the left hand side may be decomposed into individual elements then apply 9.6, 9.7 and 9.9 on each element.∎
9.6. Construction of the lifted finite element space
Recalling the normal projection operator 9.4, we define the global lifting map (Def. 7.7(b)) by
[TABLE]
and denote the restriction to each element by for each . For , we denote its lift by given by
[TABLE]
and for , we denote its inverse lift by given by
[TABLE]
For each and each , we use Def. 7.1 to construct an associated lifted surface finite element . We assume the domains are such that the set of lifted element domains defines an exact decomposition of (Def. 7.7(a)) at each .
We use the lift to define the space of lifted finite element functions (Def. 7.7(d)) by
[TABLE]
We continue by showing some basic geometric estimates. These geometric estimates have been given by (Kovács, 2018, Lemma 5.2).
Lemma 9.10**.**
Under the above smoothness assumptions, we have
[TABLE]
[TABLE]
where . Writing for the quotient between discrete and continuous surface measures so that , we have
[TABLE]
Lemma 9.11**.**
Let be the evolving subdivision defined by the lifting map and assume that is a uniformly -regular, uniformly quasi-uniform, evolving conforming subdivision. Then is also a uniformly -regular, uniformly quasi-uniform, evolving conforming exact subdivision of .
Proof.
We use the decomposition that
[TABLE]
and write for . From Rem. 7.4 and Prop. 7.14, we have to show
[TABLE]
and
[TABLE]
The second estimate follows from the smoothness of 9.4 (which follows from the smoothness of (Foote, 1984)). Next we compute directly that
[TABLE]
so that using the notations and for , we have
[TABLE]
Applying 9.24 and 9.25, we see that
[TABLE]
Clearly the right hand side of this equation is less that for sufficiently small. ∎
Lemma 9.12**.**
Let and denote its lift by . Then there exists constants such that
[TABLE]
Proof.
We apply Prop. 7.8 using Lem. 9.11. ∎
Using the lift and the discrete flow , we can define the lifted flow map (Def. 7.12(a)), lifted discrete velocity (Def. 7.12(b)) and lifted push forward maps (Def. 7.12(c) and 3.13). Lem. 9.11 implies, with Prop. 7.14, that , and are compatible pairs uniformly for . We will use the notations and for the spaces of functions smoothly evolving in time 2.2 with respect to the push-forward map in and respectively. We may also define a lifted material derivative for functions using 3.15 by
[TABLE]
Lemma 9.13**.**
The push forward map induces a new transport formula on . For we have
[TABLE]
Let be a smooth, positive-definite, diffusion tensor on , which maps the tangent space of to itself, and be a smooth tangential vector field on for all . Furthermore, we have for
[TABLE]
and for we have
[TABLE]
where and are defined as in Lem. 9.9.
Proof.
We note that the left hand side may be decomposed into individual elements then apply 9.6, 9.7 and 9.9 on each element.∎
Lemma 9.14**.**
[TABLE]
Proof.
See (Kovács, 2018, Lem. 5.2). ∎
We can characterise the lifted discrete material derivatives 9.32 using the difference between smooth, , and lifted, , velocities. The result also shows the abstract inclusions L3.
Lemma 9.15**.**
If , then and conversely if then . In either case, we have the identity
[TABLE]
Furthermore, if then .
Proof.
We use a chart and write . Note that . Then
[TABLE]
Using the notation and , we have
[TABLE]
Using the substitutions , , , and , we see
[TABLE]
Next, we take a time derivative of to see that
[TABLE]
By multiplying by , summing over and multiplying we see that .
Next, we take a time derivative of , to see that
[TABLE]
Multiplying by we see that
[TABLE]
from which we infer that
[TABLE]
Combining the previous expressions we see that
[TABLE]
Using the substitutions as above and the parametric definition of tangential gradient gives the result.
The second result can be shown by applying the tangential gradient to the basic result. ∎
9.7. The discrete problem and stability
For each , and , we assume that is an element-wise smooth symmetric diffusion tensor defined element-wise with for each . We assume that maps the tangent space of at a point into itself and is uniformly positive definite on the tangent space: There exists such that for all , , and
[TABLE]
We assume we are also given a element-wise smooth tangential vector field (with ) and element-wise smooth scalar field (with ). We assume that
[TABLE]
Example 9.16*.*
Here we are thinking of the case that , and .
We consider the following semi-discrete problem (c.f. Prob. 3.2):
Problem 9.17**.**
Given , find with and such that for every ,
[TABLE]
where
[TABLE]
We note that the assumption that is uniformly quasi-uniform regular implies that is a compatible pair when equipped with the or -norms (Lem. 6.30).
To show the properties of these bilinear forms we require one further lemma:
Lemma 9.18**.**
The discrete velocity of the discrete evolving surface 9.20 is uniformly bounded in . That is, there exists a constant such that for all
[TABLE]
Proof.
The result follows from the interpolation bound (Cor. 7.11) and the stability of the lift (Prop. 7.8). ∎
We have transport formulae on the surface .
Lemma 9.19**.**
There exists a bilinear forms and such that
[TABLE]
where
[TABLE]
and
[TABLE]
Furthermore, there exist a constant such that for all and all we have
[TABLE]
Proof.
The transport theorem 9.21 directly gives 9.45 and additionally 9.22 and 9.23 give 9.46. To see the boundedness properties, we directly apply Lem. 9.18. ∎
Theorem 9.20**.**
There exists a unique solution of the finite element scheme 9.41. The solution satisfies the stability bound:
[TABLE]
Proof.
The result is shown in the abstract setting in Thm. 3.3 so we are left to check the assumptions. The assumptions Mh1 and Mh2 follow since is simply the inner product. For Gh1, we use 9.21 and the product rule . The bound Gh2 is shown in 9.47. The map is differentiable and hence measurable Ah1. The bounds Ah2 and Ah3 follow from standard calculations and our assumptions on . Finally we have shown Bh1 and Bh2 in 9.46 and 9.48. ∎
9.8. Error analysis
The space is equipped with the following approximation property:
Lemma 9.21**.**
The interpolation operator is well defined and satisfies
[TABLE]
Proof.
We simply apply Thm. 7.10. The second result applies the theorem in the obvious way. The first result applies the theorem with noting that and the inclusions for Lem. 4.3 still hold. ∎
We can further relate the lifted material velocity with the discrete material velocity. Let evolve with velocity and . Then
[TABLE]
Then from the above calculation of we have
[TABLE]
Lemma 9.22**.**
We have the estimate:
[TABLE]
Proof.
See (Kovács, 2018, Lem 5.4). ∎
The lifting operator also defines transport formulae:
Lemma 9.23**.**
There exists bilinear forms and given by
[TABLE]
These bilinear forms satisfy the following transport formulae on :
[TABLE]
Furthermore, the two new bilinear forms are uniformly bounded in the sense that there exists a constant such that for all and all ,
[TABLE]
Proof.
The transport formulae are direct translations of Lem. 9.13. The bounds follow from the fact that is bounded uniformly from Lem. 9.22. ∎
For the remainder of this section, we will additionally assume that is uniformly smooth in space. We set .
Lemma 9.24**.**
There exists a constant such that for all and all with lifts we have
[TABLE]
Proof.
The results 9.61a, 9.61b, 9.61c and 9.61d easily follow using ideas from (Kovács, 2018, Lem 5.6) and Lem. 8.24. 9.61e and 9.61f follow directly from Lem. 9.22. ∎
Lemma 9.25**.**
[TABLE]
Furthermore for all and , we have
[TABLE]
Proof.
We recall 9.39:
[TABLE]
We combine this calculation with 9.54 to see 9.62.
We may apply the tangential gradient to the above equation and use 9.54 again to obtain
[TABLE]
9.64 follows from 9.61b and 9.63. ∎
Theorem 9.26**.**
Let be the solution of 9.14 which we assume satisfies the further regularity requirement
[TABLE]
Let be the solution of 9.41 and denote its lift by . Then we have the following error estimate
[TABLE]
Proof.
We simply check the assumptions of Thm. 3.11. We know the lift is stable from Lem. 9.12. The existence and boundedness of and are dealt with in Lem. 9.23. The interpolation properties I1 and I2 are shown in Lem. 9.21. The geometric perturbation estimates P1, P2, P3, P4, P5, P6, P4’, P5’, P7, P8 and P9 are shown in the sequence of Lem. 9.24, 9.22 and 9.25. Finally, to show B3 we follow a calculation given in the proof of (Dziuk and Elliott, 2013a, Thm. 6.2). In this setting the simpler version B3’ holds; see Rem. 3.7. Observe that for any , , that
[TABLE]
For , using the additional smoothness assumptions on , we can apply integration by parts to see that
[TABLE]
The bounds for and are obvious. This implies that
[TABLE]
10. Application III: A coupled bulk-surface parabolic system
In this section we will consider a finite element method for the coupled bulk-surface problem 1.6. The method is based on combining the isoparametric approaches from the problem in a bulk domain (Sec. 8) and the problem on a surface (Sec. 9). The discretisation will be posed in the product of a bulk isoparametric finite element space of order and a surface isoparametric finite element space of order . We take our notation from the previous sections (Sec. 9 and 8).
10.1. The domain and function spaces
We set , and We will also make use of the spaces and . We see that for each and the inclusions are uniformly continuous. We define the push forward operator by
[TABLE]
and are compatible pairs (Def. 2.2) and the spaces and (c.f. 2.1) and and (c.f. 2.2) are well defined. For , we define the strong material derivative, denoted , using 2.4. The product material derivative coincides with the product of surface and bulk material derivatives:
[TABLE]
form an evolving Hilbert triple (Def. 2.4). We also have that and are compatible pairs. For further information on this functional analytic setting see Alphonse et al. (2015b, Sec. 5.3).
Remark 10.1*.*
For the well posedness of the partial differential equation Prob. 10.2 we require that the boundary is a -hypersurface and that the flow map . See Alphonse et al. (2015b) for more details. For the approximation properties we derive we require that is a -hypersurface and that with and both of class .
We introduce a signed distance function for the boundary surface . The oriented signed distance function for is given by
[TABLE]
For each , we orient by choosing the unit normal for . Our assumptions on imply that there exists a neighbourhood of and normal projection operator given as the unique solution of
[TABLE]
See Gilbarg and Trudinger (1983, Lem. 14.16); Foote (1984) for more details.
10.2. The initial value problem
Given as in Sec. 8.2 and as in Sec. 9.2, we consider the initial value problem
Problem 10.2**.**
Find and such that
[TABLE]
Remark 10.3*.*
The problem 1.6 is recovered by setting and and and , where is a decomposition of into tangential and normal components on .
10.3. The bilinear forms and transport formulae
We set
[TABLE]
We can combine the transport formula for the bulk and surface only cases 8.9 and 9.11 for , 8.10 and 9.12 for and 8.11 and 9.13 to derive transport laws for these coupled bilinear forms First, for we have
[TABLE]
for , we have
[TABLE]
and for , we have
[TABLE]
where is given for by
[TABLE]
and is given for by
[TABLE]
We define to be
[TABLE]
We also have the estimates that there exists a constant such that for all we have
[TABLE]
10.4. Variational formulation
We consider a weak form of 1.6.
Problem 10.4**.**
Given (\emph{\mathrm{u}}_{0},v_{0})\in\mathcal{H}_{0}, find (\emph{\mathrm{u}},\emph{\mathrm{v}})\in\mathcal{W}(\mathcal{V},\mathcal{V}^{*}) such that for almost every we have
[TABLE]
Theorem 10.5**.**
There exists a unique solution pair which satisfies the stability bound:
[TABLE]
and if then
[TABLE]
Proof.
We again apply the abstract theory of Thm. 2.15 and check the assumptions. Ass. 2.9 and 2.13 are shown in Alphonse et al. (2015b, Sec 5.3). It is clear that M1 and M2 hold since is equal to the -inner product. The assumptions G1 and G2 are shown in 10.4 and 10.7. We know that the map is differentiable hence measurable which shows A1. The coercivity A3 and boundedness As4 of and boundedness of An4 follow from standard arguments since the extra cross term is clearly positive (see also Elliott and Ranner (2013, Thm. 3.2)). The existence of the bilinear forms Bs1 and Bn1 has been shown in 10.5 and 10.6 and the estimates Bs2 and Bn2 are shown in 10.8 and 10.9. ∎
10.5. Discretisation of the domain and finite element spaces
In order to define our computational method we use the construction of the isoparametric domain of order used in Sec. 8.5. This defines an evolving triangulated bulk domain (Def. 4.32(b)) equipped with an evolving conforming subdivision (Def. 4.32(a)). We will assume that is a uniformly quasi-uniform subdivision (Def. 4.35(a)). We will make the same assumptions on the domain as in Sec. 8.5 which allow us to show Prop. 8.9. Namely, we conclude that we can define an evolving bulk finite element space (8.18 and 4.34(a)) consisting of Lagrange finite elements of order (Ex. 4.7(b)) over a uniformly -regular evolving subdivision (Def. 4.35(b)). By we denote the maximum mesh diameter over time 4.24:
[TABLE]
Throughout the remainder of this section we will denote global discrete quantities with a subscript . We assume implicitly that these structures exist for each in this range (see also Rem. 4.11 and 6.11).
For , we write and for the boundary faces of :
[TABLE]
We note that is also an evolving conforming subdivision (Def. 6.27(a)) which we assume is uniformly quasi-uniform (Def. 6.29(c)). In fact, this is the construction we have previously used for an evolving surface in Sec. 9 so that we may make the same conclusions as Prop. 9.8. That is that we can define an evolving surface finite element space (9.19 and 6.28(a)) consisting of Lagrange finite elements of order (Ex. 6.7(b)) over a uniformly -regular, evolving subdivision (Def. 6.29(b)).
Our analysis will make use of the product Hilbert spaces and and the product finite element space . Using Lem. 4.23 and 6.20 we can identify elements of as a product of continuous functions on and and that . We equip with the norms:
[TABLE]
The previous constructions define a global flow map (Def. 4.32(c)) and discrete velocity (Def. 4.32(d)), with well defined surface restrictions (Def. 6.27(c)) and (Def. 6.27(d)). For each , we define the discrete push forward map by
[TABLE]
Since we have shown that and are both uniformly -regular and uniformly quasi-uniform, the spaces , and form a compatible pair (Def. 2.2). Further, we can define the spaces 2.1 and 2.2 and we can define a material derivative for functions which can be identified as 3.3:
[TABLE]
10.6. Construction of lifted finite element space
We have already constructed a bijection between the computational domain and the continuous domain . In Sec. 8.6, for each , we constructed element-wise a bijection . Furthermore, we note that the restriction of the lifting operator to , , is simply the normal projection operator which is the lifting operator used in Sec. 9.6.
For each bulk finite element , we can use Def. 5.1 to construct an associated lifted bulk finite element . We will denote the set of lifted bulk finite elements by . For each surface finite element , we can use Def. 7.1 to construct an associated lifted surface finite element . We will denote the set of lifted surface finite elements by .
For and a function pair , we define the lift by
[TABLE]
We will often write to signify that the lifting process is simply a combination the previous lifts for the surface and bulk components.
We will also make use of an inverse lift for functions on . For , we define he inverse lift of , denoted by by
[TABLE]
Lemma 10.6**.**
Let and denote their lift by . Then there exists constants such that
[TABLE]
Proof.
We simply combine the results of Lem. 8.14 and Lem. 9.12. ∎
We use the lifts to define a product space of lifted finite element functions (Def. 5.7(d) and 7.7(d)) by
[TABLE]
Using the lift and the discrete flow , we can defined a lifted flow map (Def. 5.11(a)) and lifted discrete velocity (Def. 5.11(b)), with well defined restrictions (Def. 7.12(a)) and (Def. 7.12(b)). This allows us to define a lifted push forward map (3.13, 7.12(c) and 5.11(c)):
[TABLE]
Using the previous constructions, applying Lem. 9.11 and 8.13, we see that , and are compatible pairs uniformly for . We will use the notations and for the spaces of functions smoothly evolving in time 2.2 with respect to the push-forward map in and respectively. This implies we have a well defined strong material derivative 3.15:
[TABLE]
10.7. The discrete problem and stability
The finite element method is based on the variation form 2.12 of Prob. 10.4. We assume we have and for as in Sec. 8.7 and and for as in Sec. 9.7.
Problem 10.7**.**
Given , find such that for all
[TABLE]
where for , we define
[TABLE]
and, for , we define
[TABLE]
We also have discrete transport formula from the bulk and surface cases:
Lemma 10.8**.**
There exists bilinear forms and such that for all we have
[TABLE]
and for all we have
[TABLE]
where
[TABLE]
and
[TABLE]
Further, there exists a constant such that, for all ,
[TABLE]
Proof.
We simply combine Lem. 9.19 and 8.19. ∎
Theorem 10.9**.**
There exists a unique solution pair of the finite element scheme (Prob. 10.7) which satisfies the stability bound
[TABLE]
Proof.
We apply the abstract result of Thm. 3.3 and check the assumptions. The assumptions on , Mh1 and Mh2, follow directly since is equal to the inner-product. The estimates on , Ah2 and Ah3 follow in the same manner as Thm. 10.5. The transport formulae and estimates for and , Gh1, Gh2 Bh1 and Bh2, are shown in Lem. 10.8. ∎
10.8. Error analysis
We assume in this section that , , , , , . The space of lifted finite element functions is equipped with the follow approximation property:
Lemma 10.10** (Approximation property).**
For the Lagrangian interpolation operator is well defined. Furthermore, the following bounds hold for a constant for all and :
[TABLE]
Proof.
We define the interpolation operator to be for . The proof follows by combining the result of Lem. 9.21 and Lem. 8.21. ∎
Lemma 10.11**.**
There exists a bilinear forms and given by
[TABLE]
such that
[TABLE]
Furthermore, there exists a constant such that for all and all we have
[TABLE]
Proof.
We combine the results of Lem. 9.23 and 8.22. ∎
The geometric perturbation results now follow directly by combining the appropriate results from Sec. 9.8 and 8.8.
Lemma 10.12**.**
We have the estimates
[TABLE]
In particular, this implies
[TABLE]
Proof.
We combine the results of Lem. 9.22, 9.25 and 8.23. ∎
Lemma 10.13**.**
There exists a constant such that for all and all the following holds for all with lifts :
[TABLE]
For with inverse lifts , we have
[TABLE]
Proof.
We combine the results of Lem. 9.24 and 8.24. ∎
We require one final result in order to show the error bound.
Lemma 10.14**.**
For any , let be as in 3.20 and , then for all , we have
[TABLE]
Proof.
The proof follows a similar path to Lem. 8.25. We start by fixing . We recall that for any , there exists such that
[TABLE]
and satisfies
[TABLE]
As in Sec. 3.3, we introduce such that there exists a constant such that
[TABLE]
We wish to estimate the trace of , the bulk component of , in the -norm. We consider the dual problem: Given , find such that
[TABLE]
This problem is the weak form of coupled bulk surface elliptic problems and has a unique weak solution in (see Elliott and Ranner 2013, Thm. 3.2) which satisfies the estimate
[TABLE]
The problem with is a weak form of the problem:
[TABLE]
This is a Robin problem which satisfies the regularity estimate (Ladyzhenskaya and Uraltseva, 1968; Gilbarg and Trudinger, 1983)
[TABLE]
The problem with is a weak form of a surface elliptic problem with right hand side which satisfies the regularity estimate (Aubin, 1982).
[TABLE]
Combining the two regularity estimates we see that
[TABLE]
We note that the constant here is independent of time .
We see using 10.37b and in 10.38 that
[TABLE]
For the first term on the right hand side of 10.42, we apply the boundedness of and the interpolation estimate 10.19 to see
[TABLE]
For the second term on the right hand side of 10.42, we apply the geometric estimates 10.29, 10.27 and 10.33 and the interpolation estimate 10.19 to see that
[TABLE]
Hence, combining the previous estimates and applying the dual regularity 10.41, we infer that
[TABLE]
and we have shown that
[TABLE]
Returning to 10.36, we now see that for and that
[TABLE]
The estimates for and are clear
[TABLE]
For , we apply similar reasoning to Lem. 8.25 with our estimate for 10.43 and for , we apply similar reasoning to 9.67. Combining these estimates shows the desired bound. ∎
Theorem 10.15**.**
Let be the solution of 10.10 which we assume satisfies the regularity bound
[TABLE]
Let , be the solution of the finite element scheme 10.15 and write . Then we have the following error estimate
[TABLE]
Proof.
We apply abstract Thm. 3.11 and check the assumptions. We know the lift is stable from Lem. 10.6. The existence and boundedness of and are dealt with in Lem. 10.11. The interpolation properties I1 and I2 are shown in Lem. 10.10. The geometric perturbation estimates P1, P2, P3, P4, P5, P6, P4’, P5’, P7, P8 and P9 are shown in Lem. 10.13 and 10.12. Finally, B3 is shown in Lem. 10.14. ∎
11. Numerical results
11.1. Implementation
The finite element methods were implemented using DUNE. In our numerical examples, we integrate the coefficients , , , and using a sufficiently accurate quadrature which does not affect the order of convergence of the schemes.We discretise in time using an implicit Euler time stepping scheme. The time step is scaled so that the optimal error scales are recovered. At each time step we solve the full system using the generalised minimal residual method,
The code produced to run these computations is available at
https://github.com/tranner/dune-evolving-domains
Let , and . For , we define via a parametrisation , for the unit ball in three-dimensions. The parametrisation is given by
[TABLE]
with velocity field given by
[TABLE]
The geometry is the same for each problem, which corresponds to an ellipsoidal domain growing along a single axis, but we solve in and on different parts of the domain.
For each test problem, for each iteration we complete an appropriate number of bisectional refinements (projecting boundary nodes on to the exact surface using the normal projection operator 9.4) from a macro triangulation in order to approximately half the mesh size and scale the time step to recover the optimal order of convergence – i.e. . We show the error in an -norm at the final time. The experimental order of convergence (eoc) at level is computed by
[TABLE]
Errors in an -norm demonstrate an order of convergence less and are not listed here.
11.2. Problem on a bulk domain 1.4
We set the parameters in the equation as , , and compute additional right hand sides in 1.4a and 1.4b and take appropriate initial data so that the solution is given by
[TABLE]
We compute with . The results are shown in Tab. 11.1 and 11.2.
11.3. Problem on a closed surface 1.5
We set the parameters in the equation as , , and compute an additional right hand side in 1.5a and take appropriate initial data so that the solution is given by
[TABLE]
We compute with . The results are shown in Tab. 11.3 and 11.4.
11.4. Problem on a coupled bulk-surface domain 1.6
We set the parameters in the equation as , , , for and , and , and compute additional right hand sides in 1.6a, 1.6b and 1.6c and take appropriate initial data so that the solution is given by
[TABLE]
We compute with . The results are shown in Tab. 11.5 and 11.6.
Appendix A A Faá di Bruno formula for parametric surfaces
A partition of the set is a collection of non-empty subsets such that if and . We call the order of the partition and denote by the number of elements in . We say that two sets and are ordered and write if .
We denote by the set of ordered non-empty partitions of of order :
[TABLE]
We note that we have and . That is that and each contain one partition.
For a subset , given , we write for smooth
[TABLE]
Note that it is possible to have repeated indexes which are not distinguished in this formula.
Theorem A.1**.**
Let be a surface with parametrisation over . We denote the components of by for and the components of by . Let be in and write for the function defined by for . Then and
[TABLE]
Before we show the result we will include some examples to show how this result can be interpreted and related to previous results.
Example A.2*.*
- (1)
Consider the case and is a flat hypersurface. In this case . Then the A.1 translates to
[TABLE]
We note that if and [math] otherwise and that . Finally we note that the terms involving can be reordered as
[TABLE]
where is given by
[TABLE]
The set rather than being a set of partitions maps each partitions to a vector such that is the number of sets in the partition such that for . We can make the simplification in this case since, here, the order of derivatives is not important, whereas we wish to track this in A.1. Furthermore, the mapping from partitions to vectors is not one-to-one (in particular because does not care about ordering) which results in the constant which counts how many partitions map to the same .
We recover the result of (Bernardi, 1989, Eq. 2.9) in the scalar case:
[TABLE] 2. (2)
We return to the general case, but to low numbers of derivatives. It is clear that
[TABLE]
To start to understand the inductive procedure we will apply in the proof of Thm. A.1, we construct a second derivative. We apply the product rule to see
[TABLE]
We see that the extra derivative either applies to the terms from the lower order derivative involving or else apply to the terms involving , in which case an extra first derivative of is included. The result is a sum of these two terms.
To understand further how the terms involving arise, we compute a third order derivative. Applying the same procedure as above we see
[TABLE]
Again, we see that the derivative can either be applied to the terms involving or the terms involving and the result is a sum of these terms. Computing further we see
[TABLE]
We end up with terms involving one, two and three derivatives of times sums of products of one, two and three terms involving derivatives of respectively. We notice that the derivatives on are not all arbitrary combinations of ’s and ’s. The derivative always is against , the derivative is always against or and the derivative is against , or . Moreover, we see that the minimum index against is increasing in . This property is written more formally in the definition of . Making this association rigorous is the route we take to proving the theorem.
We will show the result through the following two lemmas.
Lemma A.3**.**
The derivatives of can be written as
[TABLE]
where satisfies the recurrence relationships
[TABLE]
Proof.
For the base case A.3a, we see that
[TABLE]
Then, given A.2 holds up to derivatives, we see that
[TABLE]
Reading off coefficients gives the result. ∎
Lemma A.4**.**
It holds that
[TABLE]
Proof.
For the base case, , we see that and clearly the set of ordered partitions is .
Suppose the identity holds for derivatives up to order .
First, we consider , then
[TABLE]
We also have that there is only one partition of into subset: .
Next, we consider , then
[TABLE]
We also have that there is only partition of into ordered sets: so this case is complete.
Next, we consider . We see
[TABLE]
where if and otherwise and is given by
[TABLE]
The proof will be complete if we show .
First, let . It is clear that is a partition. There are two cases to check the ordering (i) or for (ii). For case (i), we have that and . The first assertion shows that for and the second shows that for . So we see that . For case (ii), we have that and for some . Since is an ordered partition, we see that is non empty and the smallest index in is the same as the smallest index in so that the ordering property is preserved.
Second, let . Let be such that . First, suppose that removing from results in a non-empty set. Then is a partition in . In this case, . Otherwise suppose is empty. Then, by the ordering of partitions in , we must have that and is an ordered partition of . In this case .
Thus we have shown the desired form of . ∎
Proof of Thm. A.1.
Simply combine the two previous lemmas. ∎
We could also apply a similar result using an inverse parametrisation and recover coefficients involving derivatives of . In the applications we consider, higher derivatives of are hard to estimate. As an alternative we give the following result which is based on rearranging terms in A.1.
Corollary A.5**.**
Under the same assumptions as Thm. A.1 we have:
[TABLE]
where is projection onto the tangent space of at .
Proof.
We start by noting that we have two equivalent ways of computing surface derivatives: either using the projection onto the tangent space of or using the parametrisation of . We apply each of these formulae to the function and take the derivative at to see
[TABLE]
In particular, we see that
[TABLE]
The next step of the proof is to split the right hand side of A.1 into terms involving th order derivatives of and the rest. We multiply by and sum of each in turn to see that
[TABLE]
For the first term on the right hand side we apply A.6 and that the tangential gradient is already in the tangent space to see
[TABLE]
The result then follows by simply rearranging the terms. ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Alphonse et al. (2015 a) A. Alphonse, C. M. Elliott and B. Stinner. An abstract framework for parabolic PD Es on evolving spaces. Port. Math. , 72(1):1–46, 2015 a.
- 2Alphonse et al. (2015 b) A. Alphonse, C. M. Elliott and B. Stinner. On some linear parabolic PD Es on moving hypersurfaces. Interfaces Free Bound. , 17:157–187, 2015 b.
- 3Alphonse et al. (2018) A. Alphonse, C. M. Elliott and J. Terra. A coupled ligand-receptor bulk surface system, on a moving domain: Regularity and convergence to equilibrium. SIAM Journal on Mathematical Analysis , 50(2):1544–1592, 2018.
- 4Arnold et al. (2002) D. Arnold, F. Brezzi, B. Cockburn and L. Marini. Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM Journal on Numerical Analysis , 39(5):1749–1779, 2002. doi: 10.1137/S 0036142901384162 .
- 5Aubin (1982) T. Aubin. Nonlinear analysis on manifolds, Monge-Ampère equations , volume 252 of Grundlehren der mathematischen Wissenschaften . Springer Verlag, Berlin New Yoork, 1982.
- 6Badia and Codina (2006) S. Badia and R. Codina. Analysis of a stabilized finite element approximation of the transient convection diffusion equation using an ALE framework. SIAM Journal on Numerical Analysis , 44(5):2159–2197, 2006. doi: 10.1137/050643532 .
- 7Barrett et al. (2015) J. W. Barrett, H. Garcke and R. Nürnberg. On the stable numerical approximation of two-phase flow with insoluble surfactant. ESAIM :Mathematical modelling and numerical analysis , 49(2):421–458, 2015.
- 8Barrett et al. (2017) J. W. Barrett, K. P. Deckelnick and V. M. Styles. Numerical analysis for a system coupling curve evolution to reaction diffusion on the curve. SIAM Journal on Numerical Analysis , 55(2):1080–1100, 2017. doi: 10.1137/16M 1083682 .
