cGFEM for hyperbolic integro-differential eqs.

# Continuous Galerkin finite element methods for hyperbolic integro-differential equations

Fardin Saedpanah Department of Mathematics, University of Kurdistan, P. O. Box 416, Sanandaj, Iran
March 8, 2013
###### Abstract.

A hyperbolic integro-differential equation is considered, as a model problem, where the convolution kernel is assumed to be either smooth or no worse than weakly singular. Well-posedness of the problem is studied in the context of semigroup of linear operators, and regularity of any order is proved for smooth kernels. Energy method is used to prove optimal order a priori error estimates for the finite element spatial semidiscrete problem. A continuous space-time finite element method of order one is formulated for the problem. Stability of the discrete dual problem is proved, that is used to obtain optimal order a priori estimates via duality arguments. The theory is illustrated by an example.

###### Key words and phrases:
integro-differential equation, linear semigroup theory, continuous Galerkin finite element method, convolution kernel, stability, a priori estimate.
65M60, 45K05

## 1. Introduction

We consider, for any fixed , a hyperbolic type integro-differential equation of the form

 (1.1) ¨u+Au−∫t0K(t−s)Au(s)ds=f,t∈(0,T),with u(0)=u0, ˙u(0)=u1,

(we use ‘’ to denote ‘’) where is a self-adjoint, positive definite, uniformly elliptic second order operator on a Hilbert space. The kernel is considered to be either smooth (exponential), or no worse than weakly singular, and in both cases with the properties that

 (1.2)

This kind of problems arise e.g., in the thoery of linear and fractional order viscoelasticity. For examples and applications of this type of problems see, e.g., , , and references therein.

For our analysis, we define a function by

 (1.3) ξ(t)=κ−∫t0K(s)ds=∫∞tK(s)ds,

and, having (1.2), it is easy to see that

 (1.4) Dtξ(t)=−K(t)<0,ξ(0)=κ,limt→∞ξ(t)=0,0<ξ(t)≤κ.

Hence, is a completely monotone function, since

 (−1)jDjtξ(t)≥0,t∈(0,∞),j=0,1,2,

and consequently is a positive type kernel, that is, for any and ,

 (1.5) ∫T0∫t0ξ(t−s)ϕ(t)ϕ(s)dsdt≥0.

From the extensive literature on theoritical and numerical analysis for partial differential equations with memory, we mention , , , , , and their references.

The fractional order kernels, such as Mittag-Leffler type kernels in fractional viscoelasticity, interpolate between smooth (exponential) kernels and weakly singular kernels, that are singular at origin but integrable on finite time intervals , for any , see  and references therein. This is the reason for considering problem (1.1) with convolution kernels satisfying (1.2).

In  well-posedness of a problem, similar to (1.1) with a Mittag-Leffler type kernel, was studied in the framework of the linear semigroup theory. Here we first extend the theory to prove higher regularity of the solution for more smooth kernels, such that a priori error estimates are fulfilled. We prove optimal order a priori error estimate, by energy methods, for finite element spatial semidiscrete approximate solution. This provides an alternative proof to what we presented in , and is straightforward. The continuous space-time finite element method of order one, cG(1)cG(1), is used to formulate the fully dicrete problem. A similar method has been applied to the wave equation in , where adaptive methods based on dual weighted residual (DWR) method has been studied. An energy identity is proved for the discrete dual problem, using the positive type auxiliary function . This is then used to prove and optimal order a priori error estimates by duality. This and , where a posteriori error analysis of this method has been studied via duality, complete the error analysis of this method for model problems similar to (1.1).

The present work also extend previous works, e.g., , , , on quasi-static fractional order viscoelsticity to the dynamic case. Spatial finite element approximation of integro-differential equations similar to (1.1) have been studied in  and , however, for optimal order a priori error estimate for the solution , they require one extra time derivative regularity of the solution. A dynamic model for viscoelasticity based on internal variables is studied in . The memory term generates a growing amount of data that has to be stored and used in each time step. This can be dealt with by introducing “sparse quadrature” in the convolution term . For a different approach based on “convolution quadrature”, see . However, we should note that this is not an issue for exponentially decaying memory kernels, in linear viscoelasticity, that are represented as a Prony series. In this case recurrence relationships can be derived which means recurrence formula are used for history updating, see  and  for more details. In practice, the global regularity needed for a priori error analysis is not present, e.g., due to the mixed boundary conditions, that calls for adaptive methods based on a posteriori error analysis. We plan to address these issues in future work.

In the sequel, in , well-posedness of the problem is proved and high regularity of the solution of the problem with smooth kernels is verified. In , the spatial finite element discretization is studied and, using energy method, optimal order a priori error estimates are proved. The continuous space-time finite element method of order one is applied to the problem in , and stability estimates for the discrete dual problem are obtained. These are then used to prove optimal order a priori error estimates in by duality. Finally, in , we illustrate the theory by a simple example.

## 2. Well-posedness and regularity

We use the semigroup theory of linear operators to show that there is a unique solution of (1.1), and we prove that under appropriate assumptions on the data we get higher regularity of the solution. In we quote the main framework from , to prove existence and uniqueness, to be complete. Here we restrict to pure homogeneous Dirichlet boundary condition, though the presented framework applies also to mixed homogeneous Dirichlet-Neumann boundary conditions. But it does not admit mixed homogeneuos Dirichlet nonhomogeneous Neumann boundary conditions, and this case has been studied in  for a more general problem, by means of Galerkin approximation method. Then in we extend the semigroup framework to prove regularity of any order for models with smooth kernels. To this end, we specialize to the homogeneous Dirichlet boundary condition.

### 2.1. Existence and uniqueness

We let , be a bounded convex domain with smooth boundary . In order to describe the spatial regularity of functions, we recall the usual Sobolev spaces with the corresponding norms and inner products, and we denote . We equip with the energy inner product and norm . We recall that is a selfadjoint, positive definite, unbounded linear operator, with , and we use the norms . We note that with mixed homogeneous Dirichlet-Neumann boundary conditions, we have

 V={v∈H1:v=0\ on Dirichlet boundary}.

We extend by for with to be chosen. By adding to both sides of (1.1), changing the variables in the convolution terms and defining , we get

 (2.1) ¨u(t)+(1−κ)Au(t)+∫∞0K(s)Aw(t,s)ds=f(t)−∫∞tK(s)Ah(t−s)ds,

where, we recall that . For latter use, we note that equation (1.1) can be retained from (2.1) by backward calculations.

For a given integer number , we use the Taylor expansion of order of the solution at to define the extension for . That is, we set

 (2.2) u(t)=hr(t)=r∑n=0tnn!un(0),t<0,

where we use the notation , with .

Now we reformulate the model problem (1.1) to an abstract Cauchy problem. First, we choose in (2.2), that is , and for the initial data we assume that and . Therefore, from (2.1), we have

 (2.3) ¨u(t)+(1−κ)Au(t)+∫∞0K(s)Aw(t,s)ds=f(t)−Au0∫∞tK(s) ds,

where,

 w(t,s)={u(t)−u(t−s),s∈[0,t],u(t)−u0,s∈[t,∞).

Then we write (2.3), together with the initial conditions, as an abstract Cauchy problem and prove well-posedness.

We set and define the Hilbert spaces

 W=L2,K(R+;V)={w:∥w∥2W=∫∞0K(s)∥w(s)∥2Vds<∞},Z=V×H×W={z=(u,v,w):∥z∥2Z=(1−κ)∥u∥2V+∥v∥2+∥w∥2W<∞}.

We also define the linear operator on such that, for

 Az=(v,−A((1−κ)u+∫∞0K(s)w(s)ds),v−Dw),

with domain of definition

 D(A)={(u,v,w)∈Z:v∈V,(1−κ)u+∫∞0K(s)w(s)ds∈D(A),w∈D(D)}.

Here with

Therefore, a solution of (1.1) satisfies the system of delay differential equations, for ,

This can be writen as the abstract Cauchy problem

 (2.4) ˙z(t)=Az(t)+F(t),t∈(0,T),z(0)=z0,

where and , since

 w(0,s)=u(0)−u(−s)=u(0)−h(−s)=u0−u0=0.

We note that , so that .

We quote from [7, Theorem 2.2], that generates a -semigroup of cotractions on .

###### Corollary 1.

The linear operator is an infinitesimal generator of a -semigroup of contractions on the Hilbert space .

Now, we look for a strong solution of the initial value problem (2.4), that is, a function which is differentiable a.e. on with , if , , and a.e. on .

Recalling the assumptions and , we know that if be a strong solution of the abstract Cauchy problem (2.4) with , then is a solution of (1.1) by [7, Lemma 2.1]. Hence, to prove that there is a unique solution for (1.1), we need to prove that there is a unique strong solution for (2.4). This has been proved in [7, Theorem 2.2], if is Lipschitz continuous, using the fact that the linear operator generates a -semigroup of contractions on . Moreover, for some , we have the regularity estimate, for ,

 (2.5) ∥u(t)∥V+∥˙u(t)∥≤C(∥Au0∥+∥u1∥+∫t0∥f∥ds).

### 2.2. High order regularity

In order to prove higher regularity of order (), we assume that the bounded domain is convex, and we specialize to the homogeneous Dirichlet boundary condition. Hence, the elliptic regularity estimate holds, that is

 (2.6) ∥u∥2≤C∥Au∥,u∈D(A)=H2∩V.

We note that the case is the choice for (2.5). We substitute from (2.2), with , in (2.1). Then, differentiating and using the notation , we have

 (2.7) ¨ur(t)+(1−κ)Aur(t)+∫∞0K(s)Awr(t,s)ds=fr(t)−A∂r∂tr∫∞tK(s)r∑n=0(t−s)nn!un(0) ds=fr(t)+Ar−1∑n=0un(0)Kr−n−1(t)−Aur(0)ξ(t)=:~fr(t),

with the initial data .

Recalling the initial data and , from (1.1), we have . To obtain , we differentiate of equation (1.1), and we have

 (2.8) ¨um−2(t)+Aum−2(t) −∫t0Km−2(t−s)Au(s)ds =fm−2(t)+Am−3∑n=0un(t)Km−n−3(0),t∈(0,T),

that, with , implies the initial condition

 (2.9) um(0)=fm−2(0)−Aum−2(0)+m−3∑n=0Aun(0)Km−n−3(0), m≥3.

Throughout, obviously any sum is supposed to be suppressed from the formulas, when .

###### Remark 1.

We note that, if we assume , then by Sobolev inequality, and therefore is well-defined.

###### Remark 2.

One can show, by induction and the fact that by (2.6)

 (2.10) ∥v∥≤∥v∥H2≤C∥Av∥,for v∈D(A),

we have ,

 (2.11) |Amu2k(0)|≤C(|Am+ku0|+|Am+k−1u1|+k−1∑j=0|Am+jf2k−2j−2(0)|+k−2∑j=1|Am+jf2k−2j−3(0)|),
 (2.12) |Amu2k+1(0)|≤C(|Am+ku0|+|Am+ku1|+k−1∑j=1|Am+jf2k−2j−2(0)|+k−1∑j=0|Am+jf2k−2j−1(0)|).

Now we note that, in (2.7), we have

 wr(t,s)={ur(t)−ur(t−s),s∈[0,t],ur(t)−ur(0),s∈[t,∞),

so that . Therefore, considering continuty of , we have .

Then, in the same way as in the previous section, with , we can reformulate (2.7), with , as the abstract Cauchy problem

 (2.13) ˙zr(t)=Azr(t)+Fr(t),0

where and , since .

In particular, for , we have

 F1(t)=(0,f1(t)+Au0K(t)−Au1ξ(t),0),

with initial data .

Now, we need to show that from a strong solution of the abstract Cauchy problem (2.13), for , we get a solution of the main problem (1.1). Therfore we should prove that the abstract Cauchy problem (2.13) has a unique strong solution, under certain conditions on the data. The proof is by induction, and therefore we recall some facts from , for .

###### Lemma 1.

Let be a strong solution of the abstract Cauchy problem (2.13) with . Then is a solution of (1.1).

###### Theorem 1.

There is a unique solution of (1.1) if with , , and is Lipschitz continuous. Moreover, for some , we have the regularity estimate, for ,

 (2.14) ∥˙u(t)∥V+∥¨u(t)∥≤C(∥Au0∥+∥Au1∥+∥f(0)∥+∫t0∥˙f∥ds).
###### Proof.

There exists a unique strong solution for (2.13), with , by [7, Theorem 2.4]. Hence, the proof is complete by Lemma 1. ∎

###### Lemma 2.

Let , for , be a strong solution of the abstract Cauchy problem (2.13) with . Then is a solution of (1.1).

###### Proof.

The proof is by induction. The case follows from Theorem 1.

Now, we assume that the lemma valids for some , and we prove that it holds also for . To this end, we show that if be a strong solution of (2.13) (for ) with , then is a strong solution of (2.13) with , that completes the proof by induction assumption.

Since a.e. on , we have, for ,

 ˙ur+1(t)=vr+1(t),˙vr+1(t)=−A((1−κ)ur+1(t)+∫∞0K(s)wr+1(t,s)ds)+~fr+1(t),˙wr+1(t,s)=vr+1(t)−Dwr+1(t,s), s∈(0,∞).

The first and the third equation implies that satisfies the first order partial differential equation

 wr+1t+wr+1s=ur+1t.

This, with , has the unique solution , that implies, by integration with respect to ,

 wr(t,s)=ur(t)−ur(t−s)=∫t0wr+1(τ,s)dτ.

From the first and the second equations we obtain equation (2.7) with , that is obtained from equation (2.1) by differentiating . We recall that equations (1.1) and (2.1) are equivalent, that implies equivalence of equations (2.7) and (2.8). Therefore also satisfies (2.8) with . Then, integrating with respect to , we have, for ,

 (2.15) ¨ur(t)−¨ur (0)+Aur(t)−Aur(0)−∫t0∫τ0Kr+1(τ−s)Au(s) ds dτ =fr(t)−fr(0)+Ar∑n=0∫t0un(τ) dτKr−n(0).

Now, we need to show that (2.15) implies (2.8). We note that

 ∫t0∫τ0Kr+1(τ−s)Au(s) ds dτ =∫t0∫tsKr+1(τ−s)Au(s) dτ ds =∫t0Kr(t−s)Au(s) ds−AKr(0)∫t0u(s) ds,

and

 Ar∑n=0 ∫t0un(τ) dτKr−n(0) =A(∫t0u(τ) dτKr(0)+r∑n=1(un−1(t)−un−1(0))Kr−n(0)) =AKr(0)∫t0u(τ) dτ+Ar−1∑n=0un(t)Kr−n−1(0)−Ar−1∑n=0un(0)Kr−n−1(0).

Using these and (2.9) in (2.15) we conclude (2.8), that is equivalent to (2.7). This means that, is a strong solution of (2.13) with . Hence, by induction assumption, is a solution of (1.1), and this completes the proof. ∎

In the next theorem we find the circumstances under which there is a unique strong solution of the abstract Cauchy problem (2.13), that by Lemma 2 implies existence of a unique solution of (1.1) with higher regularity. We also obtain regularity estimates, which are extensions of (2.5) and (2.14).

We note that, recalling Remark 1 and having the assumptions from the next theorem, the calculations in the proof of Lemma 2 make sense.

###### Theorem 2.

For a given integer number , let be Lipschitz continuous and with . We also, recalling , assume the following compatibility conditions:

for ,

 (2.16) Aku0∈D(A),Aku1∈V,fr−2j(0)∈H2j∩V,j=1,⋯,k, k≥1,fr−2j+1(0)∈H2(j−1)∩V,j=1,⋯,k, k≥1,

and for ,

 (2.17) Aku1∈D(A),Ak+1u0∈V,fr−2j(0)∈H2j∩V,j=1,⋯,k, k≥1,fr−2j+1(0)∈H2(j−1)∩V,j=1,⋯,k+1, k≥0.

Then there is a unique solution of (1.1).

Moreover, for some :

for , we have the regularity estimate

 (2.18)

and, for , we have the estimate

 (2.19)
###### Proof.

1. The case follows from Theorem 1. Then, for a given , we show that

• ,

• is differentiable almost everywhere on and .

These imply existence of a unique strong solution of the abstract Cauchy problem (2.13), by [11, Corollary 4.2.10], that yields existence of a unique solusion of (1.1), by Lemma 2.

2. First we note that (i) holds, if and , by the definition of . This can be verified by applying the compatibility conditions (2.16)–(2.17) in (2.11)–(2.12), using (2.10).

3. Now we prove (ii). By assumption is Lipschitz continuous. Therefore, by a classical result from functional analysis, is differentiable almost everywhere on and , since is a Hilbert space. Then, recalling the assumption and the fact that

 ˙~fr(t)=˙fr(t)+Ar∑n=0un(0)Kr−n(t),

we conclude that is differentiable almost everywhere on and , that completes the proof of (ii).

4. Hence, since generates a -semigroup of contractions on by Corollary 1, we conclude, by [11, Corollary 4.2.10], that there exists a unique strong solution for the abstract Cauchy problem (2.13). This, by Lemma 2, proves that there is a unique solution of (1.1), that completes the first part of the theorem.

5. Finally, we prove the regularity estimates (2.18) and (2.19) for , since the case follows from Theorem 1.

The unique strong solution of (2.13), is given by

 zr(t)=etAzr,0+∫t0e(t−s)AFr(s)ds,

and we recall the fact that , since is an infinitesimal generator of a semigroup of contractions on . Therefore

 ∥zr(t)∥Z≤∥zr,0∥Z+∫t0∥Fr(s)∥Zds.

Since , , and

 ∥Fr(s)∥Z=∥~fr(s)∥=∥fr(s)∥+r−1∑n=0∥Aun(0)∥|Kr−n−1(s)|+∥Aur(0)∥|ξ(s)|,

therefore we have

 ((1−κ)∥ur(t)∥2V+∥ur+1(t)∥2)1/2≤((1−κ)∥ur(0)∥2V+∥ur+1(0)∥2)1/2+∫t0∥fr(s)∥ ds+r−1∑n=0∥Aun(0)∥∫t0|Kr−n−1(s)| ds+∥Aur(0)∥∫t0|ξ(s)| ds.

Hence, considering the assumption that , we have, for some ,

 (2.20) ∥ur(t)∥V+∥ur+1(t)∥≤C(∥ur(0)∥V+∥ur+1(0)∥+∫t0∥fr(s)∥ ds+r∑n=0∥Aun(0)∥).

Since, by elliptic regularity estimate (2.6),

 ∥ur(0)∥V≤∥ur(0)∥H2≤C∥Aur(0)∥,

so we have

 ∥ur(t)∥V+∥ur+1(t)∥≤C(∥ur+1(0)∥+r∑n=0∥Aun(0)∥+∫t0∥fr(s)∥ ds),

that, by (2.10)–(2.12), implies the regularity estimates (2.18)–(2.19). Now, the proof is complete. ∎

## 3. The spatial finite elment discretization

The variational form of (1.1) is to find , such that , , and for ,

 (3.1) (¨u,v)+a(u,v)−∫t0K(t−s)a(u(s),v) ds=(f,v),∀v∈V.

Let be a convex polygonal domain and be a regular family of triangulations of with corresponding family of finite element spaces , consisting of continuous piecewise polynomials of degree at most , that vanish on (so the mesh is required to fit ). Here is an integer number. We define piecewise constant mesh function for , and for our error analysis we denote . We note that the finite element spaces have the property that

 (3.2) minχ∈Vlh{∥v−χ∥+h∥v−χ∥1}≤Chi∥v∥i,for v∈Hi∩V, 1≤i≤l.

We recall the -projection and the Ritz projection defined by

 a(Rhv,χ)=a(v,χ)and(Phv,χ)=(v,χ),∀χ∈Vlh.

We also recall the elliptic regularity estimate (2.6), such that the error estimates (3.2) hold true for the Ritz projection , see , i.e.,

 (3.3) ∥(Rh−I)v∥+h∥(Rh−I)v∥1≤Chi∥v∥i,for v∈Hi∩V, 1≤i≤l.

Then, the spatial finite element discretization of (3.1) is to find such that , and for ,

 (3.4) (¨uh,vh)+a(uh,vh)−∫t0K(t−s)a(uh(s),vh) ds=(f,vh),∀vh∈Vlh,

where and are suitable approximations to be chosen, respectively, for and in .

###### Theorem 3.

Assume that is a convex polygonal domain. Let and be, respectively, the solutions of (3.1) and (3.4). Then

 (3.5) ∥uh(T)−u(T)∥≤C∥u0h−Rhu0∥+Chl(∥u(T)∥l+∫T0∥