hp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems
|
|
- Alice Hancock
- 5 years ago
- Views:
Transcription
1 Report no. 3/11 hp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems Andris Lasis Endre Süli We consider the hp version interior penalty discontinuous Galerkin finite element method hp DGFEM for semilinear parabolic equations with mixed Dirichlet and Neumann boundary conditions. Our main concern is the error analysis of the hp DGFEM on shape regular spatial meshes. We derive error bounds under various hypotheses on the regularity of the solution, for both the symmetric and non symmetric versions of DGFEM. Subject classifications: AMSMOS: 65N1, 65N15, 65N3 Key words and phrases: hp finite element methods, discontinuous Galerkin methods, semilinear parabolic PDEs Oxford University Computing Laboratory Numerical Analysis Group Wolfson Building Parks Road Oxford, England OX1 3QD October, 3
2 Contents 1 Introduction 3 Model Problem 4 3 Function Spaces 4 4 Broken Weak Formulation 6 5 Finite Element Space 7 6 Discontinuous Galerkin Finite Element Method 8 7 hp Error Estimates 9 8 Error Analysis The Non Symmetric Version of DGFEM The Symmetric Version of DGFEM The Broken Elliptic Projector A priori Error Bounds Conclusions 33
3 3 1 Introduction Discontinuous Galerkin finite element methods DGFEMs were introduced in the early 197 s for the numerical solution of first order hyperbolic problems see [3, 6, 4, 3, 16, 17, 18, 19, 31, 3]. Simultaneously, but independently, they were proposed as non standard schemes for the numerical approximation of second order elliptic equations [9, 36, 4]. In recent years there has been renewed interest in discontinuous Galerkin methods due to their favourable properties, such as a high degree of locality, stability in the absence of streamline diffusion stabilisation for convection dominated diffusion problems [1], and the flexibility of locally varying the polynomial degree in hp version approximations, since no pointwise continuity requirements are imposed at the element interfaces. Much attention has been paid to the analysis of DG methods applied to non linear hyperbolic equations and hyperbolic systems [,13, 14], several other types of non linear equations including the Hamilton Jacobi equation [], the non linear Schrödinger equation [5], and other non linear problems [15]. The analysis of the spatial discretisation of non linear parabolic problems by the Interior Penalty type of the DGFEM see [4] has been pursued by Rivière & Wheeler in [33], where the non linearities were assumed to be uniformly Lipschitz continuous with respect to the unknown solution. The resulting error bounds were based on the projection operator described in [34], and were not p optimal in the H 1 norm. In this work we shall be concerned with the error analysis of the hp version interior penalty discontinuous Galerkin finite element method hp DGFEM, for an initial boundary value problem for a semilinear PDE of parabolic type in n spatial dimensions on shape regular quadrilateral meshes see.1 below. Here, we consider only the spatial discretisation of the problem, leaving the choice of time stepping techniques and their analysis for a future work. We shall suppose that the non linearity satisfies the local Lipschitz condition.. The paper is structured as follows. In Section we state the model problem, followed by the definition of function spaces used throughout our work Section 3. Next, we state the broken weak formulation Section 4. After selecting the finite element space that will be used for the discretisation of the model problem in space Section 5, we state the hp DGFEM Section 6. Section 7 contains the approximation theory, required in the subsequent error analysis. The error analysis of the hp DGFEM for semilinear parabolic equations is discussed in Section 8. We begin by establishing the local Lipschitz continuity of the mapping f : L q Ω L Ω. Section 8.1 contains the error analysis of the non symmetric version of the interior penalty hp DGFEM: we prove an h optimal and p suboptimal by half an order of p a priori error bound. The bound indicates that the presence of the non linearity, obeying condition., does not degrade the accuracy of the hp DGFEM in the H 1 norm. Section 8. is concerned with the derivation of the L norm error bounds in the case of the symmetric version of the interior penalty hp DGFEM. For this purpose, we first derive error bounds on the broken elliptic projector Section 8..1 defined by the symmetric version of the hp DGFEM. Section 8.. is concerned with the error analysis and derivation of the a priori error bound for the L norm, and is largely based on the techniques developed in the analysis of the non
4 4 symmetric version of the hp DGFEM. Section 9 contains some final comments on the results in this work. Model Problem Let Ω be a bounded open domain in R n, n, with a sufficiently smooth boundary Ω. We consider the semilinear partial differential equation of parabolic type u u = fu in Ω,T],.1 where u u/ t, T >, and f C 1 R. We also assume the following growth condition on the function f: fv fw C g 1 + v + w α v w for all v, w R,. where C g > and α >. Upon decomposing the boundary Ω into two parts, Γ D and Γ N, so that Γ D Γ N = Ω, we impose Dirichlet and Neumann boundary conditions respectively: u = g D u n = g N on Γ D [,T], on Γ N [,T], where n = nx denotes the unit outward normal vector to Ω at x Ω. Finally, we impose the initial condition.3 u = u on Ω {},.4 where u = u x. As the solution of this problem may exhibit blow up in finite time, we shall assume that, for the potential blow up time T, ], the time interval [,T] on which the problem is defined is bounded by the blow up time, i.e., T < T. 3 Function Spaces Since hp-dgfem is a non conforming method, it is necessary to introduce Sobolev spaces defined on a subdivision T of the domain Ω; we call such piecewise Sobolev spaces broken Sobolev spaces. A subdivision T of the domain Ω R n, n, is a family of disjoint open sets elements such that Ω =. Before we define broken Sobolev spaces, we shall introduce the basic principles of constructing a subdivision T. Let T be a subdivision of the polyhedral domain Ω R n, n, into disjoint open polyhedra elements such that Ω =, where T is regular or 1 irregular, i.e., each face of has at most one hanging node. We assume that the family of subdivisions T is shape regular see pages 61, 113, and Remark. on page 114 in [1], and require each T to be an affine image of a fixed master element ˆ, i.e., = F ˆ for all T, where ˆ is either the open unit simplex or the open unit hypercube in R n, n.
5 5 Definition 3.1 The broken Sobolev space of composite order s = {s : T } on a subdivision T of Ω is defined as W s pω, T := { u L p Ω : u W s p for all T }, s being the local Sobolev index on the element. The associated broken norm and seminorm are defined as u W s p Ω,T := u p W s p 1/p, u W s pω,t := u p Wp s 1/p. When s = s, we write W s pω, T, and for p = we denote H s W s. As our main concern are time dependent problems, we need to introduce Sobolev spaces comprising functions that map a closed bounded subinterval of R, with the interval in question thought of as a time interval, into Banach spaces. For further reference, let X denote a Banach space, with the norm, and let the time interval of interest be [,T] with T >. Definition 3. The space L p,t;x consists of all strongly measurable functions u : [,T] X with the norm and T 1/p u L p,t;x := ut dt p < for 1 p <, u L,T;X := ess.sup ut <. t T In order to move to Banach space valued Sobolev spaces, we shall define the weak derivative of a function belonging to L 1,T;X Definition 3.3 The function v L 1,T;X is the weak derivative of u L 1,T;X, written u = v, provided that, for all scalar test functions ϕ C,T, we have T T ϕtut dt = ϕtvt dt. Definition 3.4 The Sobolev space W 1 p,t;x
6 6 consists of all functions u L p,t;x such that u exists in the weak sense and belongs to L p,t;x, with the associated norm and T 1/p u W 1 p,t;x := { ut p + ut p } dt < for 1 p <, u W 1,T;X := ess.sup ut + ut. t T Further, for simplicity, we shall write H 1,T;X W 1,T;X. 4 Broken Weak Formulation Before presenting the broken weak formulation of the problem described in Section, we shall introduce some notation. Let T be a subdivision of Ω R n, n, into disjoint open polyhedra as in Section 3. By E we denote the set of all open n 1-dimensional faces of the subdivision T, containing the smallest common n 1-dimensional interfaces e of neighbouring elements. We define E int := e and E := e. e E\ Ω e E Ω Numbering the elements of the subdivision T, and choosing any internal interface e E int, there exist positive integers i, j such that i > j and elements i and j which share this interface e. We define the jump of a function u H s Ω, T across the face e and the mean value of u on e by [u] e := u e u e and u e := 1 u e + u e, respectively, with denoting the union of all open faces of the element. With each face e we associate the unit normal vector ν pointing from the element i to j when i > j; when the face belongs to E, we choose ν to be the unit outward normal vector n. Finally, we decompose the set of all faces on the boundary E into two sets, E D and E N, such that Γ D = e ED e and Γ N = e EN e. Now we are ready to introduce the broken weak formulation of the problem.1.4. We define the bilinear form B, by Bu,v := u v dx + {θ v ν [u] u ν [v]} ds + σ [u] [v] ds 4.1 Γ int Γ int + {θ v nu u nv} ds + σuv ds, Γ D Γ D
7 7 and the linear functional l by lv := g N v ds + θ Γ N v ng D ds + Γ D σg D v ds. Γ D 4. Here σ is called the discontinuity penalisation parameter and is defined by σ e = σ e for e E int E, where σ e is a non negative constant on the face e. The precise choice of σ e will be discussed in Section 8. The subscript e in these definitions will be suppressed when no confusion is likely to occur. The parameter θ here takes the values ±1. The choice of θ = 1 leads to a symmetric bilinear form B, ; we call this method a Symmetric Interior Penalty, or SIP, method. On the other hand, the choice of θ = 1 leads to a non symmetric, but coercive bilinear form B, ; we call such method a Nonsymmetric Interior Penalty, or NSIP, method. Further we shall label the bilinear form 4.1 and the linear functional 4. with indices S and NS in the symmetric and non symmetric cases respectively. Then, the broken weak formulation of the problem.1.4 reads: find u H 1,T;A such that uv dx + Bu,v fuv dx = lv, for all v H Ω, T, u = u, 4.3 where by A we denote the function space A = { w H Ω, T : w, w ν are continuous across each e E int }. 5 Finite Element Space Here we define the finite dimensional subspace of H 1 Ω, T on which the finite element method will be posed. It makes sense to construct this space in such a way that the degree of piecewise polynomials contained in the space can be different on every element of the subdivision T. This will allow us to vary the approximation order according to the local regularity of the solution on the element by changing the degree of the polynomial on elements. As we are concerned with the discontinuous Galerkin method here, we do not need to make any additional assumptions to ensure continuity of the approximation across element interfaces. Henceforth, this method will be referred to as hp DGFEM see [35] for a description of hp FEM. For a non negative integer p, we denote by P p ˆ the set of polynomials of total degree p on a bounded open set ˆ. When ˆ is the unit hypercube, we also consider Q p ˆ, the set of all tensor product polynomials on ˆ of degree p in each coordinate
8 8 direction. To each T we assign a non negative integer p the local polynomial degree and a non negative integer s the local Sobolev index. Recalling the construction of the subdivision T see Section 3, we collect the p and the F into vectors p = {p : T } and F = {F : T }, and consider the finite element space where R is either P or Q. S p Ω, T,F := { u L Ω : u F R p ˆ, T }, Discontinuous Galerkin Finite Element Method Using the finite element space S p Ω, T,F, defined in the previous section, and the broken weak formulation of the problem 4.3, the approximation u DG to the solution u of the problem.1.4, discretised by the discontinuous Galerkin finite element method in space, is defined as follows: find u DG H 1,T;S p Ω, T,F such that u DG v dx + Bu DG,v fu DG v dx = lv, for all v S p Ω, T,F, u DG = u DG, 6.1 where u DG denotes the approximation of the function u from the finite element space S p Ω, T,F, and the parameter σ in 4.1 and 4. is to be defined in the error analysis. The equation 6.1 can be interpreted as a system of ordinary differential equations in t for the coefficients in the expansion of u DG t in terms of basis functions of the finite dimensional space S p Ω, T,F. Thus, 6.1 defines an autonomous system of ordinary differential equations with C 1 and, therefore, locally Lipschitz continuous right hand side, given that f C 1 R and the other terms are linear. By the Cauchy Picard theorem this, in turn, implies the existence of a unique local solution to 6.1. Since no pointwise continuity requirement is imposed on the elements of the finite element space, the approximation u DG in S p Ω, T,F to the solution u will be, in general, discontinuous. Remark 6.1 If the continuity assumptions made in the construction of the space A are violated, i.e., u and u ν are discontinuous across the element interfaces, we have to modify the DGFEM accordingly. This could be done, for example, by performing a DGFEM discretisation on every subdomain of Ω where the continuity requirements are satisfied, and incorporating into the definition of the method transmission conditions on interfaces where discontinuities in the solutions occur. Such situations include, for example, heat transfer problems in heterogeneous or layered media or problems that contain different phases of material. There the solution u and/or the diffusive fluxes u n can be discontinuous across element interfaces. This information has to be
9 9 incorporated into the definition of the method and, in particular, into the choice of the discontinuity penalisation parameter σ, to avoid penalising physical discontinuities. 7 hp Error Estimates The first analysis of the p version of FEM for Poisson s equation was given by Babuška et al. [9], and was subsequently refined by Babuška & Suri in [7] and [8]. The analysis relied on the use of the Babuška Suri projection operator. For the special case of n =, the analysis in W s q norms was carried out by Ainsworth & Kay in [] and [3], where the approximation bounds were used for deriving a priori error bounds for p and hp version FEMs for the r Laplacian, using approximation by continuous piecewise polynomials on both quadrilateral and triangular elements. The error bounds obtained in these works contain logarithmic terms in p, and thus are only optimal in p up to a logarithmic factor. We shall proceed with the derivation of local approximation error bounds, avoiding such suboptimal logarithmic terms by using some very recent results due to Melenk [8]. From Proposition A. and Theorem A.3 in [8], we conclude the following result concerning polynomial approximation of functions defined on hypercubes. Lemma 7.1 Let Q := 1, 1 n, n 1, and let u W k q Q, where q [1, ]; then there exists a sequence of algebraic polynomials z p u R p Q, p N, such that, for any l k, u z p u W l q Q Cp k l u W k q Q, 1 q, 7.1 where C > is a constant, independent of u and p, but dependent on q and k. To derive the general hp estimates for the projection operator u z p u, recall from Sections 3 and 5 the construction of the subdivision T of the computational domain Ω. Let ˆ be the n dimensional open unit hypercube, which we shall call the reference element. We construct each element T via an affine mapping from the reference element = F ˆ, based on scaling each coordinate of the reference element by the factor h. We shall also need the following result. Lemma 7. Suppose that T is an n dimensional parallelepiped of diameter h, and that u Wq k for some k and T. Define û Wq k ˆ by the rule ûˆx ˆ = uf ˆx ; Then where s = minp + 1,k. inf û ˆv W ˆv R p ˆ q k ˆ Chs n/q u W k q,
10 1 Proof. See [7], Lemma 4.4, and [3], Lemma 1. Assume that k is an integer. If k =, then the result follows by bounding the left hand side of the inequality by û L q ˆ and scaling to u L q. Suppose, therefore, that k 1. For any ˆv R p ˆ, we have k û ˆv W k q ˆ û ˆv Wq s ˆ + û W l q ˆ, l =s +1 with the convention that if s = k then the summation is over an empty index set of l. Using Theorem in [1], we obtain k inf û ˆv W k ˆv R p ˆ q ˆ û W l q ˆ. l =s Scaling back to the element T, we obtain the result for integer k. The result for general k follows by a standard function space interpolation argument. Now we are ready to state our main result concerning the approximation properties of the projection operator u z p u. Theorem 7.3 Suppose that T is an n dimensional parallelepiped of diameter h, and that u Wq k for some k and T ; then, there exists a sequence of algebraic polynomials zp h u R p, p 1, such that for any l, with l k, u zp h u hs l C Wq l p k l u W k q, 1 q, 7. and, for q =, u z h p u C hs L e 1 p k 1 u H k, 7.3 u z h p u C hs L e 3 p k 3 u H k, 7.4 where e is any face edge e, s = minp + 1,k, and C is a constant independent of u, h, and p, but dependent on k = max k and q. Proof. See also [7]. Let u W k q and define û W k q ˆ by the rule ûˆx ˆ = uf ˆx. First, we note that, by Lemma 4.1 in [7], for any ˆv R p ˆ, we have the property that ẑh p ˆv = ˆv. By Lemma 7.1, 7.1, we have, for l k, û ẑh p û W l q ˆ Cp k l û W k q ˆ.
11 11 Noting that ẑh p ûˆx = zp h uf ˆx, and applying Lemma 7. with ˆv R p ˆ, we obtain û ẑh p û = W l q ˆ û ˆv ẑh Cp k l p û ˆv inf ˆv R p ˆ W l q ˆ û ˆv W k q Thus, by a scaling argument, for m l k, we have u zp h u h s m u W k q, and therefore u z h p u W m q Cp k l Wq l Cp k l h s l and hence 7.. By setting q = in 7. and using the trace inequality ˆ Cp k l h s n q u W k q. u W k q, u L C h 1 u L + u 1 L u 1 L, we obtain 7.3 and Error Analysis This section is be concerned with the derivation of a priori error bounds for the initial boundary value problem for the semilinear parabolic equation described in Section. We shall assume that the polynomial degree vector p, with p 1 for each T, has bounded local variation, i.e., there exists a constant ρ 1 such that, for any pair of elements and which share an n 1 dimensional face, ρ 1 p p ρ. 8.1 We also recall our regularity assumptions on the subdivision T : namely, T is shape regular, and regular or 1 irregular. We shall consider the error analysis of the hp version of the discontinuous Galerkin finite element method on shape regular meshes. In particular, we shall derive a priori error bounds for both the symmetric and the non symmetric version of DGFEM. Let us begin with the following lemma which establishes the local Lipschitz continuity of the non linearity f, required in the a priori error analysis of the hp version of DGFEM 6.1 for the model problem.1.4. Lemma 8.1 Let f C 1 R satisfy the growth condition. with < α < when n = and < α /n when n 3, and suppose that < q <. Let αq ˆq = max q,. q
12 1 Then, there exists a positive constant C = Cα,C g,q, Ω such that fu fv L Ω C u v L q Ω 1 + u α L q αq + Ω v α L q αq Ω 8. for all u, v Lˆq Ω. Suppose that q = α + 1; then ˆq = α + 1. Moreover, if n =, < α < then < ˆq <, and if n 3, < α /n then < ˆq n/n. Proof. From., by Hölder s inequality, we have fu fv L Ω = Ω fu fv dx C g As 1 /q = q /q and q >, we have Ω u v 1 + u + v α dx Cg u v q q dx 1 + u + v α 1 q 1 dx Ω Ω q fu fv L Ω C g u v q q dx + u + v Ω Ω1 αq q q dx q = Cg u v L q Ω 1 + u + v αq q q dx Ω q = Cg u v L q Ω 1 + u + v αq αq α q dx Ω Cg u v L q Ω Ω q αq + u L q αq + v Ω L q αq Ω C u v L q Ω 1 + u α L q αq + Ω v α L q αq, Ω 1 q. and hence 8. for all u,v Lˆq Ω. The statement in the final sentence of the lemma follows from our hypothesis on the range of α and the fact that q = α + 1. Hypothesis A. Let f C 1 R satisfy the growth condition. with < α < when n =, and < α /n when n 3. We define q = α + 1. With this hypothesis in mind, we can remove the dependence on q in the constant C in Lemma 8.1 in terms of α. For the sake of clarity of the exposition, in the rest of the paper we shall confine ourselves to the case of n 3. Our proofs can be easily adjusted to cover the case of n = with < α <. α
13 The Non Symmetric Version of DGFEM Let the bilinear form B be as in 4.1. Here we shall be concerned with the non symmetric version of DGFEM corresponding to θ = 1 in 4.1, so we write B NS, in place of B,. We begin our error analysis with the following definition. Definition 8. We define the quantity DG on H 1 Ω, T, associated with the DGFEM, as follows: w DG := w L Γ + σw ds + σ [w] ds D Γ int where σ is a non negative function on Γ. Remark 8.3 Let us observe some properties of DG defined above. 1. If σ > on Γ, then DG defines a norm in H 1 Ω, T.. If σ = on Γ, then DG defines a seminorm in H 1 Ω, T. 1, Clearly, w DG = B NSw,w, for all w H 1 Ω, T. The first step in the error analysis is to decompose the error u u DG, where u denotes the analytical solution, as u u DG = ξ + η, where ξ Πu u DG, η u Πu, with Π defined element wise by Πu := Πu, and Π denoting an appropriate projection operator on the element. Thus, using the triangle inequality for the H 1 norm, we have u u DG H 1 Ω,T η H 1 Ω,T + ξ H 1 Ω,T. 8.4 We assume for simplicity that the initial value is chosen as u DG = Πu, 8.5 and thus ξ =. We shall proceed by deriving a bound on ξ H 1 Ω,T in terms of norms of η. Then, by choosing a suitable projection operator Π, we shall be able to use the bounds on various norms of the projection error η derived in Section 7 to deduce an a priori error bound for the method. Let us prove the continuity of the bilinear form B NS, which will also provide the necessary bound for our error analysis. Lemma 8.4 Let T be a shape regular subdivision of Ω and assume that the parameter σ is positive on Γ int Γ D ; then, the following inequality holds for all v H 1 Ω, T and
14 14 w S p Ω, T,F, with C a positive constant that depends only on the dimension n and the shape regularity of T : { B NS v,w C w DG σ v ds + σ [v] ds + v L Γ D Γ int + τv + 1 L Γ D σ v L Γ D + τ [v] + 1 L Γ int σ v L Γ int }1, 8.6 where τ e = p e /h e and h e is the diameter of a face e E int E D ; when e E D the contribution from outside Ω in the definition of τ e is set to. Proof. See also [1]. Let us decompose B NS v,w I + II + III + IV, where I v w dx, II {v w n v nw} ds, Γ D III {[v] w ν v ν [w]} ds, Γ int IV σvw ds + σ [v] [w] ds. Γ D Γ int For the term I we have and for the term IV we have that I w DG 1 v L, IV w DG σ v Γ ds + σ [v] ds. 8.8 D Γ int To deal with the term II, we first note that II 1 1 v L γ Γ D γ w L Γ D + 1 v σ L Γ D 1 1 σ w L Γ D 1
15 15 for any set of positive numbers {γ : T }. Here we can apply the inverse inequality w L Γ D K p h w L, 8.9 where K depends only on the shape regularity of T see Schwab [35], Theorem 4.76, On letting γ = h /p and defining τ e = p /h e for an n 1 dimensional face e Γ D, we obtain II C w DG τv + 1 L Γ D σ v Similarly, we have III C w DG τ [v] + 1 L Γ int σ v By collecting the results we have the desired bound. L Γ D L Γ int Now, let us derive a bound on the H 1 norm of the error u u DG Lemma 8.5 Let T be a shape regular subdivision of Ω and assume that f C 1 R satisfies Hypothesis A. Suppose further that the positive parameter σ is defined on Γ int Γ D and σ e = σ e h 1 e on each face e E int E D. In addition, suppose that a the local polynomial degree p on each T ; b the local Sobolev smoothness k 3.5 on each T ; c the hp mesh is quasi uniform in the sense that there exists a positive constant C such that h h max C p min. p Then, for all t [,T], there exists h > such that for all h,h ], h = max h, the following inequality holds, with C a positive constant that depends only on the domain Ω, the quasi uniformity of T, on the final time T, the exponent α in the growth condition for the function f, and the Lebesgue and Sobolev norms of u over the time interval [,T]: t u u DG s H 1 Ω,T ds C t { ηs L + ηs L α+1 + ηs H 1 + σηs + σ [ηs] + L Γ D L Γ int τηs L Γ D + 1 ηs σ + τ [ηs] } + 1 L Γ int σ ηs ds 8.1 L Γ D L Γ int
16 16 where τ e = p e /h e and h e is the diameter of a face e E int E D, in which for e E D the contribution from outside Ω is set to. Proof. From the formulation of the hp DGFEM 6.1, for all v S p Ω, T,F, we have u DG v dx + B NS u DG,v = fu DG v dx + l NS v On the other hand, the broken weak formulation 4.3 of the problem can be rewritten as uv dx + B NS Πu,v = Π fuv dx + l NS v Π u uv dx + B NS Πu u,v. Upon subtracting 8.13 from 8.14 and choosing v = ξ, we obtain ξξ dx + B NS ξ,ξ = {fu fu DG } ξ dx ηξ dx B NS η,ξ. By noting that ξξ dx = 1 d dt we can rewrite the above expression as 1 d dt ξ L Ω + ξ DG ξ L = 1 d dt ξ L Ω, {fu fπu}ξ dx + {fπu fu DG }ξ dx + ηξ dx + B NSη,ξ By the Cauchy Schwarz and Young inequalities, with ε 1 >, we have ηξ dx η L 1 ξ L 1 ε 1 η L Ω + 1 ε 1 ξ L Ω, and, by the same argument, with ε, ε 3 >, {fu fπu}ξ dx ε fu fπu L Ω + 1 ξ L ε Ω, {fπu fu DG }ξ dx ε 3 fπu fu DG L Ω + 1 ξ L ε Ω. 3
17 17 Further, by Lemma 8.1, upon absorbing all constants into C and noting the definition of q in Hypothesis A, we have fu fπu L Ω C η L q Ω 1 + u α L q αq + Ω Πu α L q αq Ω C η L q Ω 1 + u α L q αq + Ω Πu α L q αq Ω = C η L q Ω 1 + u α L q αq + u η α Ω L q αq Ω C η L q Ω 1 + u α L q αq + Ω η α C η L q Ω = C η L α+1 Ω, L q αq Ω where the constant C > depends only on the domain Ω, the growth condition for the function f, and on Lebesgue norms of u over the time interval [,T]. By Lemma 8.4 and Young inequality, with ε 4 >, we have the bound B NS η,ξ ε 4 ξ DG + 1 ε 4 F 1 η, where F 1 η := C ηs L + σηs L Γ D + σ [ηs] L Γ int + τηs + 1 L Γ D σ ηs L Γ D + τ [ηs] + 1 L Γ int σ ηs L Γ int. Applying these bounds on the right hand side of 8.15 and absorbing all constants into C 1 and C, we obtain d dt ξ L Ω + ε 4 ξ DG C 1 Fη + C ξ L Ω + ε 3 fπu fu DG L Ω, 8.16 where Fη := F 1 η + η L α+1 Ω + η L Ω. To bound fπu fu DG L Ω, we first note that, by the same argument as above, fπu fu DG L Ω C ξ L α+1 Ω 1 + u DG α L q αq, Ω
18 18 where the constant C > depends only on the domain Ω, the growth condition for the function f, and on Lebesgue norms of u over the time interval [,T]. Let us choose u DG = Πu, thus giving ξ =, and let < t T be the largest time such that u DG exists for all t [,t ] and ξ H 1 Ω,T 1 for all t [,t ]; existence of such a t is guaranteed by the Cauchy Picard theorem. Since, by Hypothesis A, αq/q n/n, this implies that u DG α L αq q Ω Const. for all t [,t ] by the broken Sobolev Poincaré inequality see Theorem 3.7 in [7] 1 ; here Const. is a constant that is independent of the discretisation parameters and t, and only depends on Sobolev norms of u over the time interval [,t ]. This implies that fπu fu DG L Ω C ξ DG, where the constant C > depends only on the domain Ω, the growth condition for the function f, and on Lebesgue and Sobolev norms of u over the time interval [,t ]. On choosing ε 4 + ε 3 C 1, 8.16 takes the form d dt ξ L Ω + ξ DG C 1 Fη + C ξ L Ω, 8.17 with the constant C 1 > depending only on the domain Ω, the growth condition for the function f, and on Lebesgue and Sobolev norms of u over the time interval [,t ]. Upon integrating from to t t and noting that ξ =, this yields t t t ξt L Ω + ξs DG ds C 1 Fηs ds + C ξs L Ω ds, 8.18 with the constant C 1 as above. According to this inequality, if Fη were zero, we would have ξ L Ω = for all t [,t ]. More generally, by choosing an appropriate projection operator Π, we can make Fη as small as we like for example, by fixing the local polynomial degree p on each element T and reducing h = max h. 1 Using the notation of the cited paper, we define Ψ as in Example 3.6 of that paper, with ψ L Ω, ψ on Γ N and Ψξ C h 1 e e E D ξ ds. e
19 19 Let us choose C 3 = C α and h > so small that, for all h h and t [,t ], the following inequality holds: t C 1 Fηs ds < T e C 3T C 1 h inv C max, p where C inv is the constant from the inverse inequality ξ H 1 Ω,T C inv max p h ξ L Ω We note in passing that in order to be able to extract the factor max h /p from Fη, we need hypotheses a and b above. Hence 8.18 becomes t ξt L Ω + ξs DG ds < T e C 3T C 1 h inv C max p t + C ξs L Ω ds, which, by the Gronwall Bellman inequality, implies that ξt L Ω < C 1 inv C max By the inverse inequality 8.19 we have that, ξt H 1 Ω,T < C max h max p h for all t [,t p ]. p h = C max h min p h, p for all t [,t ], which, by the quasi uniformity hypothesis c above, is 1. Hence, then, for h h, we have ξ H 1 Ω,T < 1 for all t [,t ]. By continuity of the mapping t ξt H 1 Ω,T, the assumption t < T implies that either ξt H 1 Ω,T 1 for all t [,T], or that there exists a time t t,t] such that ξt H 1 Ω,T = 1. In either case, we arrive at a contradiction with the fact that t is the largest time in the interval [,T] such that, for all t [,t ], we have ξt H 1 Ω,T 1. Thus we deduce that t = T, for < h h. From 8.18 by the Gronwall Bellman inequality we obtain t t ξt L Ω + ξs H 1 Ω,T ds C Fηs ds, t T,
20 and hence t t ξs H 1 Ω,T ds C Fηs ds, t T, with the constant C > depending only on the domain Ω, the quasi uniformity of T, on the time T, the growth condition for the function f, and on Lebesgue and Sobolev norms of u over the time interval [,T]. Employing the triangle inequality yields t t } u u DG s H 1 Ω,T { η ds C H 1Ω,T + Fηs ds, t T, and hence 8.1. Our next result concerns the accuracy of the hp version NSIP DGFEM 6.1. Theorem 8.6 Let Ω R n, n, be a bounded polyhedral domain, T = {} a shape regular and quasi uniform subdivision of Ω into n parallelepipeds, and p a polynomial degree vector of bounded local variation. Let each face e E int E D be assigned a positive real number σ e = p e h e, 8. where h e is the diameter of e, with the convention that for e E D the contributions from outside Ω in the definition of σ e are set to. Suppose that the function f C 1 R, that f satisfies the growth condition. for some positive constant C g, and that Hypothesis A holds. Then, if u,t H k with k 3.5 on each T, there exists h > such that for all h,h ], h = max h, and all t [,T], the solution u DG,t S p Ω, T,F of the NSIP DGFEM 6.1 satisfies the following error bound: u u DG L,T;H 1 Ω,T C h s u p k 3 X ; 8.1 with 1 s minp + 1,k, p on each T, where C is a positive constant depending only on the domain Ω, the shape regularity and quasi uniformity of T, the time T, the growth condition on the function f, the parameter ρ in 8.1, on k = max k, and on the Lebesgue and Sobolev norms of u over the time interval [,T]; the norm u X signifies the collection of norms u L,T;H k + u L,T;H k 1. Proof. Let us choose the projector Π to be the projection operator u z h p u, defined in Section 7. From Theorem 7.3, inequalities , we have the estimates η L C hs 1 u p k 1 H k, η H 1 C hs u p k H k, η L C hs 3 u p k 3 H k, η L C hs u p k H k.
21 1 Let us collect all the terms on the right hand side of 8.1, except η L α+1 : t I C { ηs L + ηs H 1 + σηs L Γ D + σ [ηs] + L Γ int τηs + 1 L Γ D σ ηs + τ [ηs] + 1 L Γ int σ ηs L Γ D L Γ int From the above approximation results, by choosing σ e as in 8., noting 8.1 and the shape regularity of T to relate h e to h, and taking the maximum over t [,T], we obtain I C { h s p k h u s L,T;H k 1 + p k h s 1 h p k 1 + p C } } u L,T;H k h s p k 3 ds. u X, 8. with 1 s minp + 1,k, p on each T, where C is a positive constant depending only on the domain Ω, the shape regularity and quasi uniformity of T, the time T, the growth condition for the function f, the parameter ρ in 8.1, on k = max k, and on the Lebesgue and Sobolev norms of u over the time interval [,T]. Further, by the broken Sobolev Poincaré inequality [7] we have the bound η L α+1 Ω C η L + h 1 e [η] ds + h 1 e η ds, e E e int e E e D and thus by the above approximation bounds, by noting the shape regularity of T to relate h e to h, we obtain η L α+1 C h s u p k H k, with the constant C as above. Applying this bound to the right hand side of 8.1, noting 8., and taking the maximum over t T, we obtain the desired bound. Remark The estimate 8.1 is optimal in h and p suboptimal by p 1.. By the broken Sobolev Poincaré inequality, the same bound holds for the L norm of the error. The bound in this case is not hp optimal.
22 3. From the error bound we conclude that the presence of the non linearity f, satisfying the conditions stated in Section, does not diminish the rate of hp convergence rate in the H 1 norm compared to the linear case. 8. The Symmetric Version of DGFEM The symmetric version of the interior penalty discontinuous Galerkin finite element method appeared in the literature much earlier than the non symmetric formulation, see Wheeler [36]. It was not widely accepted as an effective numerical method until very recently, due to an additional condition on the size of the penalty parameter which is required in order to ensure the coercivity of the bilinear form of the method; this will be discussed in the next section. The renewed interest in the symmetric formulation of the IP DGFEM is due to the optimality of its convergence rate in the L norm and for linear functionals of the solution. The non symmetric formulation of the IP method suffers from lack of adjoint consistency see [6, 5], and results in suboptimal a priori error bounds in the L norm and in linear functionals of the solutions. The symmetric version, due to its adjoint consistency, does not suffer from these drawbacks. We start our a priori error analysis in the L norm by deriving the error bounds on the broken elliptic projector defined by the symmetric version of the interior penalty discontinuous Galerkin finite element method. This part of the L norm error analysis is crucial, as it will allow us to remove the terms in the error bound, containing the H 1 seminorm, which would otherwise result in a suboptimal convergence rate in the L norm The Broken Elliptic Projector Consider the boundary value problem for the elliptic equation in the form u = in Ω, u = g D on Γ D, u n = g N on Γ N, 8.3 with Γ D Γ N = Ω, Γ D having positive measure, and g D H 1 Γ D, g N L Γ N. We shall also assume that the solution u exists, that it is unique, and that u A. In view of Section 4, the SIP formulation of the DGFEM for this problem is find u DG S p Ω, T,F such that B S u DG,v = l S v for all v S p Ω, T,F, 8.4 where the symmetric bilinear form B S is defined by 4.1, and the linear functional l S is defined by 4., with θ = 1. Let us check whether and under what conditions the solution u DG to 8.4 exists and is unique.
23 3 The proof of continuity of the symmetric bilinear form B S u,v is essentially the same as in the non symmetric case see Lemma 8.4. The coercivity, though, requires further investigation. For the symmetric bilinear form 4.1 with θ = 1, we have, for any w S p Ω, T,F, B S w,w = w L + σw Γ w w n ds+ σ [w] [w] w ν ds. D Γ int Clearly the integrands in the last two terms need not be positive unless σ is chosen sufficiently large: the purpose of the analysis that now follows is to explore just how large σ needs to be to ensure coercivity of B S, over S p Ω, T,F S p Ω, T,F. For any positive number τ e we have w w n ds Γ D e E D τ e w ds + e e τe 1 w n ds. Omitting the summations, the second term on the right hand side can be further bounded by using the inverse inequality 8.9, the shape regularity condition to relate h to h e, where is the element whose face is e E D and the bounded local variation condition to relate p to p e, by absorbing all constants into C τ, we obtain e τe 1 w n ds e τ 1 e w n ds τe 1 p C e τ w L h. e Similarly, for the term involving interior faces, we have [w] w ν ds τ e [w] ds + τ 1 p e C e τ w Γ int h L w L, e e E int e using the restriction imposed by the bounded local variation condition 8.1: here and are the two elements that have e as their common face. Now letting we get τ 1 e := 1 p 1 C e n n 1 τ for e E D E int, h e B S w,w 1 w L Γ + σ τw ds + σ τ [w] ds. D Γ int Thus the symmetric bilinear form B S u,v is coercive if σ e τ e = n n 1 C τ p e h e. The factor n n 1 stems for the fact that in n dimensions the summation over e E int may count, any one element, n n 1 times, as we allow one hanging node per interface.
24 4 Choosing σ e appropriately, i.e., σ e = C σ p e h e 8.5 with the constant C σ > large enough, σ e τ e will be ensured and by the Lax Milgram theorem the solution to 8.4 then exists and is unique. In view of the above arguments, we conclude that the SIP DGFEM solution of the problem 8.3 uniquely determines the projection operator Π e on A onto the finite element space S p Ω, T,F with the property for u A B S u Π e u,v = for all v S p Ω, T,F. 8.6 Next, we state the approximation error bounds in the H 1 and L norms for the broken elliptic projector Π e. Lemma 8.8 Let Ω R n, n, be a bounded polyhedral domain, T = {} a shape regular subdivision of Ω into n parallelepipeds, and suppose that u H k for some Sobolev index k and T. Let Π e u be the projection of u A onto S p Ω, T,F, defined by 8.6, with p for T, and σ e chosen as in 8.5. Then, the following error estimate holds: Furthermore, if Ω is convex, then u Π e u H 1 Ω,T C u Π e u L Ω C max h s u p k 3 H k. 8.7 h p h s u p k 3 H k, 8.8 where s = minp + 1,k, and the constant C is independent of u, p and h, but dependent on k = max k and C σ. Proof. By recalling the definition of the DG norm 8.3, we have, from the assumption on σ e, that u DG B S u,u for all u A, and thus by writing u Π e u = u Πu Πu Π e u = η + ξ, where the projection operator Π will be chosen later, taking v ξ in the definition of the broken elliptic projector 8.6, we deduce that ξ DG B S ξ + η η,ξ B S ξ + η,ξ + B S η,ξ = B S η,ξ. By continuity of the bilinear form B S η,ξ see Lemma 8.4 and the comments above, after applying Young inequality we have ξ DG C ση + σ [η] + L Γ D L Γ int η L + τη L Γ D + 1 η σ + τ [η] + 1 L Γ int σ η, 8.9 L Γ D L Γ int
25 5 where τ e = n n 1 C τ p e /h e, h e is the diameter of a face e E int E D, and for e E D the contribution from outside Ω is set to. Further, by noting that ξ L ξ DG, and employing the triangle inequality 8.4, we obtain the bound u Π e u H 1 C ση + σ [η] + L Γ D L Γ int η H 1 + τη + 1 L Γ D σ η L Γ D + τ [η] + 1 L Γ int σ η, 8.3 L Γ int Let us choose Π to be the u z h p u see Section 7. From Theorem 7.3, inequalities , we have the estimates η L C hs 1 u p k 1 H k, η L C hs 3 u p k 3 H k, η H 1 C hs u p k H k. Applying these inequalities to the right hand side of 8.3, choosing σ e as in 8.5, noting the bounded local variation condition 8.1 and the shape regularity of T to relate h e to h, we obtain u Π e u H 1 Ω,T C h s p k h s 1 h p k 1 + p u H k, and hence 8.7. Let us note that the same bound 8.7 is also valid for the DG norm u Π e u DG ; this follows from 8.9 and the fact that u Π e u DG C h s p k 3 u H k. To estimate u Π e u L Ω, we shall use the Aubin Nitsche duality argument see [11]. Let, signify the L inner product. Then, for every g L Ω, by the Cauchy Schwarz inequality we have u Π e u,g u Π e u L Ω g L Ω, and therefore u Π e u L Ω = sup g L Ω g u Π e u,g g L Ω
26 6 Further, let the function w H Ω be the solution of the problem w = g in Ω, w = on Γ D, w n = on Γ N, 8.3 with g L Ω, and Γ D, Γ N as in 8.3. Then the SIP DGFEM formulation of this problem is find w A such that B S w,v = l g v for all v H Ω, T, where B S w,v is defined by 4.1 with θ = 1, and l g v = g,v + l S v, with l S v defined by 4. with θ = 1 and g D =, g N = : clearly, then, l S v = for all v in H Ω, T. Consider the SIP DGFEM approximation of 8.3 in the form find w DG S p Ω, T,F such that B S w DG,v = l g v for all v S p Ω, T,F. By Galerkin orthogonality, we have and thence B S u Π e u,w DG =, u Π e u,g = g,u Π e u = l g u Π e u = B S w,u Π e u = B S u Π e u,w = B S u Π e u,w Πw, where Π is the projection operator u z h p u. Further, by Lemma 8.4, 8.6, and by noting that the bilinear form B S, is symmetric, we have u Π e u,g B S u Π e u,w Πw C u Π e u DG { σ w Πw ds + σ [w Πw] ds + w Πw L Γ D Γ int + τw Πw + 1 L Γ D σ w Πw L Γ D + }1 τ [w Πw] + 1 L Γ int σ w Πw L Γ int 8.33 with τ e = n n 1 C τ p e /h e.
27 7 By the previous argument, we have the estimate u Π e u DG C h s p k 3 and from Theorem 7.3, inequalities , we have the estimates u H k, 8.34 w Πw L C h3 w p 3 H, w Πw L C h p w H, w Πw L C h w p H. Applying these inequalities and the estimate 8.34 to the right hand side of 8.33, choosing σ e as in 8.5 and noting the bounded local variation condition 8.1 and the shape regularity of T to relate h e to h, we obtain u Π e u,g C h s p k 3 u H k Further, by noting that for a suitable constant C > we have h w H p C max h w H p max = C h p w H 1. h w H p Ω, and, recalling that Ω is convex, on employing elliptic regularity, we obtain and therefore u Π e u,g C max h p h s p k 3 u H k 1 g L Ω, u Π e u,g g L Ω C max h p h s p k 3 u H k 1. Noting 8.31, taking the supremum over g L Ω, g, and squaring the resulting expression yields A priori Error Bounds Having defined the broken elliptic projector and obtained the respective approximation error bounds, we are ready to state our main result about the accuracy of the symmetric version of the hp DGFEM.
28 8 Theorem 8.9 Let Ω R n, n, be a bounded convex polyhedral domain, T = {} a shape regular subdivision of Ω into n parallelepipeds, and p a polynomial degree vector of bounded local variation. Let each face e E int E D be assigned a real positive number σ e = C σ p e h e, 8.35 where h e is the diameter of e, with the convention that for e E D the contributions from outside Ω in the definition of σ e are set to, and C σ is sufficiently large. Suppose that the function f C 1 R and obeys the growth condition. for some positive constant C g, and that Hypothesis A holds. Then, if u,t H k, k, T, for t T there exists h > such that for all < h h, h = max h, the solution u DG,t S p Ω, T,F of the SIP DGFEM 6.1 obeys the following error bounds: and u u DG L,T;L Ω C max ess.sup u u DG DG C t T h p + max h s p k 3 h s u p k 3 X αn h α+1 u X h s u p k 3 L,T;H k, 8.37 with 1 s minp + 1,k, p 1, for T, where C is a positive constant depending only on the domain Ω, the shape regularity of T, the final time T, the growth condition for the function f, the parameter ρ in 8.1, the Lebesgue and Sobolev norms of u, and k = max k ; the norms u X 1, signify the collection of norms u L,T;H k Ω + u L,T;H k Ω + u L,T;H k Ω and u L,T;H k + u L,T;H k, respectively. Proof. By the same argument as in the proof of Lemma 8.5, upon subtracting 8.13 from 8.14 and choosing v = ξ, we obtain ξ L Ω + B Sξ, ξ = {fu fu DG } ξ dx η ξ dx B S η, ξ Let us choose the projection operator Π to be the broken elliptic projector Π e. Then, by definition 8.5, B S η, ξ =. With the constant C σ in 8.35 chosen large enough, the symmetric bilinear form B S, is coercive, and therefore defines an inner product on H 1 Ω, T, which induces the norm DG on this space. Hence we deduce that B S ξ, ξ = 1 d dt ξ DG.
29 9 Thus, we can rewrite 8.38 in the form ξ L Ω + 1 d dt ξ DG = {fu fu DG } ξ dx η ξ dx η ξ dx + {fu fπ e u} ξ dx + {fπ e u fu DG } ξ dx By the Cauchy Schwarz and Young inequalities, we have 1 η ξ 1 dx η L ξ L ε 1 η L Ω + 1 ξ L ε Ω, 1 and {fu fπ e u} ξ dx ε fu fπ eu L Ω + 1 ξ L ε Ω, {fπ e u fu DG } ξ dx ε 3 fπ eu fu DG L Ω + 1 ξ L ε Ω, 3 with ε 1,ε,ε 3 >. Next, by the result of Lemma 8.1, we have, upon absorbing the constants into C, and noting Hypothesis A, fu fπ e u L Ω C η L q Ω 1 + u α L q αq + Π eu α Ω L q αq Ω C η L q Ω 1 + u α L q αq + Π eu α Ω = C η L q Ω 1 + u α L q αq + u η α Ω L αq C η L q Ω 1 + u α L q αq + Ω η α C η L q Ω = C η L α+1 Ω, L αq q Ω L q αq Ω q Ω where the constant C > depends only on the growth condition for the function f, on Lebesgue norms of u over the time interval [,T]. Choosing ε 1,ε,ε 3 such that ε ε 1 + ε 1 3, and inserting the above bounds into 8.39, we obtain d dt ξ DG C 1 η L Ω + η L α+1 Ω + C fπ e u fu DG L Ω. 8.4
30 3 To bound fπ e u fu DG L Ω we note that, by the same argument as above, we have fπ e u fu DG L Ω C ξ L α+1 Ω 1 + ξ α L q αq, Ω where the constant C > depends only on the growth condition for the function f, on Lebesgue norms of u over the time interval [,T]. Let us choose u DG = Π e u, thus having ξ =, and let < t T be the largest time such that the solution ξt DG of 8.38 and thus u DGt exists and ξ DG 1 for t [,t ]; the existence of such t is guaranteed by the Cauchy Picard theorem from the theory of ODEs. By Hypothesis A, we have α +1 n/n, and hence by the broken Sobolev Poincaré inequality, fπ e u fu DG L Ω C ξ DG. Inserting this bound into 8.4, we obtain the differential inequality d dt ξ DG C 1 η L Ω + η L α+1 Ω + C ξ DG, 8.41 which, upon integrating from to t t and noting that ξ =, yields ξt DG C 1 t { } ηs L Ω + ηs L α+1 Ω t ds + C ξs DG ds. 8.4 By Lemma 8.8, the first argument on the right hand side can be bounded in terms of h and p. Fixing the polynomial degree p for all T and denoting < h = max h, let us define C 3 = C α, and let h > be small enough so that for all h h we have C 1 t { } ηs L Ω + ηs L α+1 Ω ds T e C 3T. Thus, for h h and t [,t ], from 8.4 we have ξt DG < T e C 3T + C 3 t ξs DG ds; using the Gronwall Bellmann inequality, we deduce that ξt DG < 1 for all t [,t ] with h h. By continuity of the mapping t ξt DG, the assumption t < T implies that either ξt DG 1 for all t [,T], or that there exists a time t t,t] such that ξt DG = 1. In either case, we have a contradiction with the fact that t is the largest time in the interval [,T] such that, for all t [,t ], we have ξt DG 1. Thus we deduce that t = T for < h h.
31 31 Taking into account this fact, setting h h, and applying the Gronwall Bellman inequality to 8.4 gives us the following bound: ξt DG C t { } ηs L Ω + ηs L α+1 Ω ds, t T, 8.43 where the constant C > depends only on the domain Ω, the growth condition for the function f, the time T, on Lebesgue and Sobolev norms of u over the time interval [,T]. Further, by the broken Sobolev Poincaré inequality, we have the bound η L α+1 Ω C η DG, and, employing the triangle inequality, we thus obtain t { } u u DG t DG C ηt DG + ηs L Ω + ηs DG ds, t T, with the constant C as above. By the results of Lemma 8.8 we have that η L Ω max C h p h s u p k 3 H k and η DG C h s u p k 3 H k, with 1 s minp + 1,k, p 1, for T. Inserting these bounds in the above error bound, denoting u X 1 := u L,T;H k Ω + u L,T;H k Ω + u L,T;H k Ω, and taking the maximum over t [,T] yields From 8.43, by the broken Sobolev Poincaré inequality, we deduce that ξt L Ω C t { } ηs L Ω + ηs L α+1 Ω ds, t T. Employing the triangle inequality yields, for all t T, t { } ut u DG t L Ω C ηt L Ω + ηs L Ω + ηs L α+1 Ω ds, with the constant C as above Further, by the Sobolev inequality see [1], we have, for 1 α + 1 n/n, n 3, and 1 α + 1 <, n =, η L α+1 ˆ C η H 1 ˆ, where ˆ is the unit reference element the unit hypercube. By scaling back from the reference element, we obtain η L α+1 C h n 1 α η L + h1+n α+1 1, η H 1
32 3 and thus, upon squaring and summing over T, taking the square root and noting that 1 1 q a i q a i, q, we obtain η L α+1 Ω C max i 1 hn α+1 1 Inserting this inequality into 8.44 gives us and t u u DG t L Ω C ηt L Ω + + max hn 1 α+1 ηs L Ω + max i η L Ω + max h1+n 1 α+1 1 { ηs L Ω h+n 1 α+1 η H 1 Ω,T ηs H 1 Ω,T }. ds From Lemma 8.8, error bound 8.8, for k and s = minp + 1,k, we have ηt L Ω max C h h s ut p p k 3 H k max hn 1 α+1 ηt L Ω C Similarly, from 8.7 we have max h+n 1 α+1 max ηt H 1 Ω,T C max h +n 1 α+1 p h+n 1 α+1 h s ut p k 3 H k. h s ut p k 3 H k, Inserting these error bounds into 8.45 and taking the maximum over t [,T], we obtain u u DG L,T;L Ω C max h p h s u p k 3 X + max αn h α+1 h s u p k 3 L,T;H k, with 1 s minp + 1,k, p 1, for T, where the constant C > depends only on the domain Ω, the shape regularity of T, the time T, the parameter ρ in 8.1 the growth condition for the function f, on k = max k, and on Lebesgue norms of u over the time interval [,T]; here we denote u X := u L,T;H k Ω + u L,T;H k Ω, and hence 8.37.
33 33 9 Conclusions This work was concerned with the spatial discretisation of initial boundary value problems with mixed Dirichlet and Neumann boundary conditions for second order semilinear equations of parabolic type by the hp version interior penalty discontinuous Galerkin finite element method. Our goal was to derive hp version a priori error bounds. For this purpose, we derived hp version error bounds in the L and broken H 1 norms for the non local broken elliptic projection operator. We also developed the techniques of handling the non linearity in the error analysis of the hp version interior penalty discontinuous Galerkin finite element method, which allows for the proofs to be conducted on the entire time interval of existence of the solution. These enabled us to prove general error bounds for hp version discontinuous Galerkin finite element methods symmetric and non symmetric variants on shape regular meshes. The bounds, in the H 1 norm at least, are optimal in h and slightly suboptimal in p. To the best of our knowledge, these are the first error bounds of this kind for semilinear parabolic equations with a non linearity of such general type. With these bounds, we have shown that the presence of the non linearity, satisfying certain growth conditions, does not degrade the convergence rate in the H 1 norm compared to the rates obtained in the linear case. In the case of the symmetric version of the DGFEM, an attempt of the L analysis has been made; here, the impact of the non linearity on the optimality of the convergence rate is clearly seen, as the presence of the non linear term introduces a non optimal term into the error bound. Acknowledgement Part of this work was carried out while visiting the Isaac Newton Institute for Mathematical Sciences, Cambridge. We wish to express our sincere gratitude to the Institute for its support.
Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method
Report no. 03/06 Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method Emmanuil H. Georgoulis and Endre Süli Oxford University Computing Laboratory Numerical
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationarxiv: v1 [math.na] 19 Nov 2018
QUASINORMS IN SEMILINEAR ELLIPTIC PROBLEMS JAMES JACKAMAN AND TRISTAN PRYER arxiv:1811.07816v1 [math.na] 19 Nov 2018 (1.1) Abstract. In this note we examine the a priori and a posteriori analysis of discontinuous
More informationA posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model
A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationA-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems
Report no. NA-0/1 A-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems Paul Houston 1, Endre Süli, and Thomas P. Wihler
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationThe Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations
The Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations Iain Smears INRIA Paris Linz, November 2016 joint work with Endre Süli, University of Oxford Overview
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
More informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationANALYSIS OF AN INTERFACE STABILISED FINITE ELEMENT METHOD: THE ADVECTION-DIFFUSION-REACTION EQUATION
ANALYSIS OF AN INTERFACE STABILISED FINITE ELEMENT METHOD: THE ADVECTION-DIFFUSION-REACTION EQUATION GARTH N. WELLS Abstract. Analysis of an interface stabilised finite element method for the scalar advectiondiffusion-reaction
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationFinite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem
Journal manuscript No. (will be inserted by the editor Finite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem Constantin Christof Christof Haubner Received:
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationA Posteriori Existence in Adaptive Computations
Report no. 06/11 A Posteriori Existence in Adaptive Computations Christoph Ortner This short note demonstrates that it is not necessary to assume the existence of exact solutions in an a posteriori error
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationCONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS
J. Korean Math. Soc. 34 (1997), No. 3, pp. 515 531 CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS S. K. CHUNG, A.K.PANI AND M. G. PARK ABSTRACT. In this paper,
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationhp-dgfem FOR SECOND-ORDER MIXED ELLIPTIC PROBLEMS IN POLYHEDRA
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718(XX)0000-0 hp-dgfem FOR SECOND-ORDER MIXED ELLIPTIC PROBLEMS IN POLYHEDRA DOMINIK SCHÖTZAU, CHRISTOPH SCHWAB, AND THOMAS P. WIHLER
More informationNUMERICAL PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2011/12 Christian Clason January 31, 2012 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationKey words. optimal control, heat equation, control constraints, state constraints, finite elements, a priori error estimates
A PRIORI ERROR ESTIMATES FOR FINITE ELEMENT DISCRETIZATIONS OF PARABOLIC OPTIMIZATION PROBLEMS WITH POINTWISE STATE CONSTRAINTS IN TIME DOMINIK MEIDNER, ROLF RANNACHER, AND BORIS VEXLER Abstract. In this
More informationu = f in Ω, u = q on Γ. (1.2)
ERROR ANALYSIS FOR A FINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET BOUNDARY CONTROL PROBLEMS S. MAY, R. RANNACHER, AND B. VEXLER Abstract. We consider the Galerkin finite element approximation of
More informationarxiv: v1 [math.oc] 12 Nov 2018
EXTERNAL OPTIMAL CONTROL OF NONLOCAL PDES HARBIR ANTIL, RATNA KHATRI, AND MAHAMADI WARMA arxiv:1811.04515v1 [math.oc] 12 Nov 2018 Abstract. Very recently Warma [35] has shown that for nonlocal PDEs associated
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationThomas Apel 1, Ariel L. Lombardi 2 and Max Winkler 1
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANISOTROPIC MESH REFINEMENT IN POLYHEDRAL DOMAINS: ERROR ESTIMATES WITH DATA IN
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationDIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2013/14 Christian Clason July 15, 2014 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationThe continuity method
The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationAnalysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems
Analysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 340 April 2012 Key words: transport problem, Nitsche method, XFEM,
More informationPartially Penalized Immersed Finite Element Methods for Parabolic Interface Problems
Partially Penalized Immersed Finite Element Methods for Parabolic Interface Problems Tao Lin, Qing Yang and Xu Zhang Abstract We present partially penalized immersed finite element methods for solving
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationA Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 UNCORRECTED PROOF
1 A Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 Herbert Egger 1 and Christian Waluga 2 4 1 Center of Mathematics, Technische Universität München, Boltzmannstraße 3, 85748 5 Garching
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More informationFINITE ELEMENT METHODS
FINITE ELEMENT METHODS Lecture notes arxiv:1709.08618v1 [math.na] 25 Sep 2017 Christian Clason September 25, 2017 christian.clason@uni-due.de https://udue.de/clason CONTENTS I BACKGROUND 1 overview of
More informationA GAGLIARDO NIRENBERG INEQUALITY, WITH APPLICATION TO DUALITY-BASED A POSTERIORI ESTIMATION IN THE L 1 NORM
27 Kragujevac J. Math. 30 (2007 27 43. A GAGLIARDO NIRENBERG INEQUALITY, WITH APPLICATION TO DUALITY-BASED A POSTERIORI ESTIMATION IN THE L 1 NORM Endre Süli Oxford University Computing Laboratory, Wolfson
More informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
More informationA metric space X is a non-empty set endowed with a metric ρ (x, y):
Chapter 1 Preliminaries References: Troianiello, G.M., 1987, Elliptic differential equations and obstacle problems, Plenum Press, New York. Friedman, A., 1982, Variational principles and free-boundary
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More informationEquations paraboliques: comportement qualitatif
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Equations paraboliques: comportement qualitatif par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 25/6 1 Contents
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION
Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE
More informationUniversité de Metz. Master 2 Recherche de Mathématiques 2ème semestre. par Ralph Chill Laboratoire de Mathématiques et Applications de Metz
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Systèmes gradients par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 26/7 1 Contents Chapter 1. Introduction
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationA u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
More informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
More information