BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
|
|
- George Roberts
- 5 years ago
- Views:
Transcription
1 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 discretization of the Laplace operator is obtained by enriching the space element wise with a non-conforming quadratic bubble. This enriched space allows for a wider range of pressure spaces. We prove optimal convergence estimates and local conservation of both mass and linear momentum independent of numerical parameters. 1. Introduction Discontinuous Galerkin (DG methods for incompressible flow has been studied in Hansbo and Larson [9] in the framework of Nitsche s method and inf-sup stable velocity pressure pairs. The analysis was extended by Toselli to the hp-framework using mixed or equal order stabilized formulations in [12]. Local discontinuous Galerkin methods with equal order velocity and pressure spaces stabilized using penalty on the interelement pressure jumps was proposed by Cockburn et al. [6]. The Navier-Stokes equations discretized using DG has recently been given a full analysis in the framework of domain decomposition on non-matching meshes using DG techniques by Girault et al. [8]. In this paper we extend our previous work on low order discontinuous Galerkin methods for scalar second order elliptic problems to the case of incompressible flow problems [2, 3]. Using piecewise affine discontinuous finite elements enriched with non-conforming bubbles we can eliminate all stabilization terms from the formulation without compromising the adjoint consistency. The upshot is that linear momentum is conserved locally and optimal convergence in the L 2 -norm may be proven using a duality argument. This is a consequence of the fact that the vectors of the velocity gradient matrix are functions in the lowest order Raviart- Thomas space for our choice of velocity finite element space (see also [1]. Several choices for the pressure space are possible without introducing a penalty term. Indeed we can use either globally continuous, piecewise affine functions, piecewise constant, discontinuous functions or functions from a direct sum of these two sets of functions and still satisfy the inf-sup condition uniformly with respect to the mesh size. Depending on the choice of pressure space slightly different results may be obtained. When the pressure space consists of continuous functions we prove that the stresses are continuous. Using spaces with discontinuous functions for the pressure on the other hand leads to a method that also enjoys local mass conservation and Date: January 7, Key words and phrases. discontinuous Galerkin methods, Stokes problem, bubble stabilization. 1
2 2 E. BURMAN AND B. STAMM for which the divergence of the non-pressure stress has optimal convergence. We give a unified analysis for these three choices of pressure space and we then discuss the differences of these approaches from numerical and analytical point of view. 2. Notation Let Ω be a convex polygon (polyhedron in three space dimensions in R d, d =2, 3, with outer normal n. Let K be a subdivision of Ω R d into non-overlapping d- simplices κ and denote by N K the number of simplices of the mesh. Suppose that each κ K is an affine image of the reference element κ, i.e. for each element κ there exists an affine transformation T κ : κ κ. Let F i denote the set of interior faces ((d 1-manifolds of the mesh, i.e. the set of faces that are not included in the boundary Ω. The set F e denotes the faces that are included in Ω and define F = F i F e. Define by N F = card(f and N Fi = card(f i the number of faces resp. interior faces of the mesh. Assume that K is shape-regular, does not contain any hanging node and covers Ω exactly. For an element κ K, h κ denotes its diameter and for a face F F, h F denotes the diameter of F. Set h = max κ K h κ and let h be the function such that h κ = h κ and h F = h F for all κ K and F F. For a subset R Ω or R F, (, R denotes the L 2 (R scalar product, R = (, 1/2 R the corresponding norm, and s,r the H s (R norm. The element-wise counterparts will be distinguished using the discrete partition as subscript, for example (, K = κ K (, K. For s 1, let H s (K be the space of piecewise Sobolev H s functions and denote its norm by s,k. Further let us denote H s (Ω = [H s (Ω] d and H s (K =[H s (K] d. Moreover the following space is defined { } L 2 0(Ω = v L 2 (Ω v dx =0. Let v =(v 1,...,v d H 1 (K, then we define v κ [L 2 (κ] d d by ( v i,j = xj v i,1 i, j d, for each κ K. Based on the scalar L 2 -product we define (v, w K = d (v i,w i K and ( v, w K = i=1 Ω d ( xj v i, xj w i K i,j=1 for all v, w H 1 (K. Further let us define the jump and average operators. Fix F F i and thus F = κ 1 κ 2 with κ 1,κ 2 K. Let v H 1 (K, v H 1 (K and denote by v 1, v 2 resp. v 1, v 2 the restriction of v, v to the element κ 1, κ 2, i.e. v i = v κi resp. v i = v κi, i =1, 2, and denote by n 1, n 2 the exterior normal of κ 1 resp. κ 2. We then define the average and jump operators by and {v} = 1 2 (v 1 + v 2, {v} = 1 2 (v 1 + v 2, {v } = 1 2 (v 1 + v 2, { v } = 1 2 ( v 1 + v 2, [v] = v 1 n 1 + v 2 n 2, [v] = v 1 n 1 + v 2 n 2, [v ] = v 1 n 1 + v 2 n 2, [ v ] = v 1 n 1 + v 2 n 2,
3 BUBBLE STABILIZED DG METHOD 3 where a b R d d is defined by (a b i,j = a i b j,1 i, j d, for all a, b R d. On outer faces F F e we define them by {v} = v, [v] = vn, {v} = v, [v] = v n, {v } = v, [v ] = v n, { v } = v, [ v ] = vn, where n is the outer normal of the domain Ω. Further we introduce some additional notation of the jump and average operators. Let n F {n 1, n 2 } be arbitrarly chosen but fixed and introduce the vectorial quantities, indexed by v, { v } v = { v }n F and [v ] v = [v ]n F. This notation is necessary to define some projection, see Lemma 3.8. Observe that using these definitions that (1 [v ] : { w } = [v ] v { w } v for all v H 1 (K, v, w H 1 (K and where a : b = d i,j=1 a ijb ij for all a, b R d d. Additionally remark that and thus [v ] = [v ] v n F, [v] = [v ]n F n F (2 [v ] F = [v ] v F, [v] F c [v ] F for some mesh independent constant c>0. Lemma 2.1 (Integration by parts. Let v H 1 (K and v, w H 1 (K, then (v, w K = ( v, w K +({v}, [w] F + ([v], {w} Fi, ( v, w K = ( v, w K +({ v }, [w ] F + ([ v ], {w } Fi. Proof. Element-wise integration by parts and applying the definitions of the jump and average operators leads to the result. Respecting the equalities (1 leads to the following corollary. Corollary 2.2 (Integration by parts. Let v H 1 (K and v, w H 1 (K, then ( v, w K = ( v, w K +({ v } v, [w ] v F + ([ v ], {w } Fi. Let us denote by V p h p 0 defined by 3. Bubble stabilized finite element space the standard discontinuous finite element space of degree V p h = { v h L 2 (Ω vh κ P p (κ, κ K }, where P p (κ denotes the set of polynomials of maximum degree p on κ. Consider then the enriched finite element space V bs = V 1 h { v h L 2 (Ω vh (x =α x x, α V 0 h }, where x =(x 1,...,x d denotes the physical variables. Let us additionally define some functional spaces that consists of functions only defined on the skeleton of the mesh: W 0 h = { v h L 2 (F vh F P 0 (F, F F }. Define also the vectorial versions V p h =[V p h ]d, V bs =[V bs ] d and W 0 h =[W 0 h ]d.
4 4 E. BURMAN AND B. STAMM Let v [L 2 (F] m, m {1, d, d 2 }, and define by v the L 2 -projection of v onto [W 0 h ]m, i.e. ( v, w h F =(v, w h F, w h [W 0 h ]m Properties of the bubble stabilized finite element space. Let us discuss some important properties of the space V bs. Lemma 3.1. For v h V bs we have that v h V 0 h. Moreover the application : V bs /V 1 h V 0 h is bijective. Proof. Observe that w h = 0 for all w h Vh 1 0 and that (αx x = 2dα Vh where d is the dimension of Ω. Let us denote by RT 0 (κ the local lowest order Raviart-Thomas space on κ. The following Lemma holds. Lemma 3.2. For all v h V bs with v h =(v h,1,..., v h,d there holds v h,i κ RT 0 (κ, κ K, and for all κ K and r h =(r h,1,...,r h,d with r h,i RT 0 (κ, there exists v h V bs such that v h,i κ = r h,i for all 1 i d. Proof. For the scalar case we refer to the proof of Lemma 3.4 in [2] and the vectorial version is constructed componentwise. Corollary 3.3. For v h V bs we have that { v h } v W 0 h, [ v h ] W 0 h. Moreover the applications { } v : V bs W 0 h and [ ] : V bs W 0 h are surjective. Lemma 3.4. There is a constant c>0 independent of h such that for all v h V bs there holds h 1 2 [vh ] 2 F ( h c 1 2 [vh ] 2 F + v h 2 K. Proof. We refer to the proof of Lemma 4.1 in [2] for the scalar case. Its vectorial conterpart follows from the componentwise construction. Lemma 3.5 (Poincaré inequality. There is a constant c>0 independent of h such that for all v h V bs there holds ( v h 2 K c P h 1 2 [vh ] 2 F + v h 2 K. Proof. We refer to the proof of Corollary 4.2 in [2] for the scalar case. Its vectorial counterpart follows from the componentwise construction. Observe that Lemma 3.4 and 3.5 are only valid on discrete spaces and thus does not hold for functions in H 1 (K. Thus we define a norm for H 1 (K by v 2 = v 2 K + h 1 2 [v] 2 F for all v H 1 (K. Nevertheless observe that there exists a constant c d > 0 such that (3 c d v h 2 v h 2 K + h 1 2 [vh ] 2 F v h 2 for all v h V bs.
5 BUBBLE STABILIZED DG METHOD Technical lemmas. Let us cite some well known results. For the proofs we refer to [4]. Lemma 3.6 (Inverse inequality. Let v h [V bs ] m, m {1,d}, then there exists a constant c I > 0 independent of h such that c 1 I h 2 v h K h v h K c I v h K. Lemma 3.7 (Trace inequality. Let v [H 1 (K1] m1 and v h [V bs ] m1, m 1 {1,d} resp. w [H 1 (K1] m2 and w h [V bs ] m2, m 2 {d, d 2 }, then there holds ( {v} F + [v] F c T h 1 2 v K + h 1 2 v K, {v h } F + [v h ] F c T h 1 2 vh K and {w } F + [w ] F c T ( h 1 2 w K + h 1 2 w K, {w h } F + [w h ] F c T h 1 2 wh K where c T > 0 is a constant independent of h Projections. Moreover denote by π p : L 2 (Ω V p h the L2 -projection onto V p h defined by π k (v h w h dx = v h w h dx w h V p h Ω Ω and by π p : L 2 (Ω V p h its vectorial counterpart defined by (π pv i = π p v i,1 i d. Then π p and π p satisfy the following approximation result: Let v H p+1 (K resp. v H p+1 (K, then (4 (5 v π p v K + h (v π p v K ch p+1 v p+1,k, v π p v K + h (v π p v K ch p+1 v p+1,k. Denote the Crouzeix-Raviart space by CR and its vectorial counterpart by CR. Additionally let us denote by i c : H 1 (Ω CR the vectorial Crouzeix-Raviart interpolant satisfying the following approximation result: if v H 2 (Ω, then (6 v i c v K + h (v i c v K ch 2 v 2,K. Further denote by C h : L 2 (Ω V 1 h,c the Clément interpolant [5] where V 1 h,c = { vh C 0 ( Ω vh κ P p (κ, κ K } is the continuous finite element space of degree 1. Remind that the Clément interpolant satisfies (7 v C h v K + h (v C h v K ch γ+1 v γ+1,k for all v H γ+1 (K, γ {0, 1}. Note that one can prove that the Clément interpolant conserves the mean of a function over Ω, i.e., Ω C hv(x dx = v(x dx, Ω and thus C h (L 2 0(Ω L 2 0(Ω. Note also that C h is H 1 -stable, i.e. (C h v K c v K. We denote by C h the vectorial version of C h sharing all properties. Furthermore, we present the following projection which will be used in the analysis.
6 6 E. BURMAN AND B. STAMM Lemma 3.8. Let a h Vh 0 and b h, c h Wh 0 be fixed. Then, there exists a unique function φ h V bs such that π 0 φ h = a h, (8 { φ h } v F = b h F F F, {φ h } F = c h F F F i. Moreover φ h satisfies the following stability result (9 h 1 φ h 2 K + φ h 2 c ( h 1 a h 2 K + h 1 2 bh 2 F + h 1 2 ch 2 Fi. Proof. Let us first establish the a priori estimate. Firstly by the trace inequality observe that (10 h 1 2 [φh ] 2 F c h 1 φ h 2 K and using (5 that (11 h 1 φ h 2 K h 1 π 0 φ h 2 K + h 1 (φ h π 0 φ h 2 K h 1 a h 2 K + c φ h 2 K for some constant c > 0. Secondly integrate by parts, use Lemma 3.1 and Corollary 3.3 φ h 2 K = ( φ h,π 0 φ h K +({ φ h } v, [φ h ] v F + ([ φ h ], {φ h } Fi = ( φ h, a h K +([φ h ] }{{} v, b h F + ([ φ h ], c h Fi. }{{}}{{} I II III Applying (11 and the Cauchy-Schwarz, the inverse or the trace and Young s inequality for each term yields respectively I c I φ h K h 1 a h K 1 4 φ h 2 K + c2 I h 1 a h 2 K II c T h 1 φ h K h 1 2 bh F 1 4 φ h 2 K + 1 4c h 1 a h 2 K + c c 2 T h 1 2 bh 2 F, III c T φ h K h 1 2 ch Fi 1 4 φ h 2 K + c 2 T h 1 2 ch 2 F i, and thus φ h 2 K c ( h 1 a h 2 K + h 1 2 bh 2 F + h 1 2 ch 2 F i, which, combined with (3, (10 and (11, completes the a priori estimate. To conclude the proof, it now suffices to observe that (8 is a square linear system of size d (N K + N F + N Fi. Hence, existence and uniqueness of a solution of the linear system are equivalent. Let us denote by Aw = b the square linear system and assume that there is a vector w 1 and w 2 such that Aw i = b, i =1, 2. Further let us denote the difference between them by e = w 1 w 2 and therefore Ae = 0. The a priori estimate (9 implies that e = 0 and thus the solution is unique and hence the matrix is regular. 4. Bubble stabilized discontinuous Galerkin method for Stokes problem Consider the steady Stokes problem: find u H 1 (Ω and p L 2 0(Ω such that u + p = f in Ω, (12 u = 0 in Ω, u = g on Ω, where f H 1 (Ω and g H 1/2 ( Ω such that Ω g n ds = 0. This setting ensures a unique solution to the model problem (12, see [7].
7 BUBBLE STABILIZED DG METHOD Bubble stabilized discontinuous Galerkin method. Let Q h L 2 0 (Ω be some scalar finite element space that will be precised later. We introduce the bubble stabilized Galerkin method by: find (u h,p h V bs Q h such that (13 a(u h, v h +b(p h, v h b(q h, u h =F (v h,q h (v h,q h V bs Q h where the linear form F (, and the bilinear forms a(, and b(, are defined by a(u h, v h = ( u h, v h K ( { u h }, [v h ] F ([u h ], { v h } F, b(p h, v h = (p h, v h K +({p h }, [v h ] F, F (v h,q h = (f, v h K (g n,q h Fe (g, v h n Fe. Remark 4.1. The discrete solution u h and p h of (13 satisfies the following local linear momentum conservation property { u h }n κ ds + {p h }n κ ds = f dx κ for all κ K and where n κ denotes the outer normal of κ. κ Lemma 4.2 (Consistency. Let the couple (u,p H 2 (Ω H 1 (Ω be the exact solution of (12 and let (u h,p h V bs Q h be the approximation defined by (13. Then, there holds that for all (v h,q h V bs Q h. a(u u h, v h +b(p p h, v h b(q h, u u h = 0 Proof. By integration by parts, Lemma 2.1. Details are left to the reader. 5. Convergence analysis Assume for simplicity homogeneous boundary conditions, i.e. g = 0, in this section. We will specify the choice for the pressure space Q h to get inf-sup stable approximations. Let us first introduce three possible choices for Q h. The first alternative to define the pressure approximation space is as the continuous piecewise linear finite element space defined by Q 1 h,c = { v h C 0 ( Ω L 2 0(Ω vh κ P 1 (κ, κ K } and the second one is the space of piecewise constant functions Q 0 h,d = { v h L 2 0(Ω vh κ P 0 (κ, κ K }. Further let us also consider the direct sum of the above defined spaces Q 1 h,c Q0 h,d. Remark that h q h K and h 1 2 [q h ] Fi is a norm for q h Q 1 h,c resp. q h Q 0 h,d and thus q h 2 q = h q h 2 K + h 1 2 [qh ] 2 F i is a norm for q h Q 1 h,c Q0 h,d. For the following we chose Q h = Q 1 h,c Q0 h,d and observe that the cases Q h = Q 1 h,c and Q h = Q 0 h,d are covered by the analysis as well. Thus let us define the following triple norm for all v h q h V bs Q h. v h,q h 2 = v h 2 + q h 2 q κ
8 8 E. BURMAN AND B. STAMM Proposition 5.1 (Inf-sup condition. There exists a constant c>0 independent of h such that there holds (u h,p h V bs Q h a(u h, v h +b(p h, v h b(q h, u h c u h,p h sup. 0 (v h,q h V bs Q h v h,q h Proof. For the proof of Proposition 5.1 we introduce the following two lemmas. Lemma 5.2. There exists a constant c>0 independent of h such that there exists for each fixed couple (u h,p h V bs Q h a couple (v h,q h V bs Q h with c u h,p h 2 a(u h, v h +b(p h, v h b(q h, u h. Lemma 5.3. There exists a constant c>0 independent of h such that for each (u h,p h and (v h,q h defined by Lemma 5.2 there holds v h,q h c u h,p h. Indeed combining Lemma 5.2 and 5.3 yields a(u h, v h +b(p h, v h b(q h, u h c u h,p h 2 c u h,p h v h,q h. Proof of Lemma 5.2. Firstly fix u h V bs and p h Q h. Observe that chosing v h = u h and q h = p h in (13 yields a(u h, u h (14 = u h 2 K 2( { u h }, [u h ] F u h 2 K 2 h 1 2 { uh } F h 1 2 [uh ] F u h 2 K 2c T u h K h 1 2 [uh ] F 1 2 u h 2 K c 1 h 1 2 [uh ] 2 F using Corollary 3.3, the Cauchy-Schwarz, the trace and Young s inequality. Further let w h V bs be the projection defined by Lemma 3.8 with a h = 0, b h = h 1 [u h ] v and c h = 0. By its properties, integration by parts, Lemma 3.1, Corollary 3.3 and (2 we get (15 a(u h, w h = ( u h, w h K + ([ u h ], {w h } Fi ([u h ] v, { w h } v F = ( u h, π 0 w h K + ([ u h ], {w h } Fi ([u h ] v, { w h } v F = h 1 2 [uh ] 2 F and again by integration by parts b(p h, w h = ( p h, π 0 w h K ([p h ], {w h } Fi = ([p h ], {w h } Fi since p h Vh 0 and by the property of the projection w h. Further since p h Q h we write p h = p h,c + p h,d with p h,c Q 1 h,c and p h,d Q 0 h,d and thus b(p h, w h = ([p h,d ], {w h } Fi =0 since [p h,c ] = 0 and by the property of the projection w h. Let z h be the projection defined by Lemma 3.8 with a h = h 2 p h, b h = 0 and c h = h[p h ] and observe that b(p h, z h = ( p h, π 0 z h K ([p h,d ], {z h } Fi = h p h 2 K + h 1 2 [ph ] 2 F i = p h 2 q.
9 BUBBLE STABILIZED DG METHOD 9 Also, applying integration by parts combined with the Cauchy-Schwarz and Young s inequality, yields a(u h, z h = ( u h, π 0 z h K + ([ u h ], {z h } Fi ([u h ] v, { z h } v F = ( u h, h 2 p h K ([ u h ], h[p h ] Fi c u h K ( h p h K + h 1 2 [ph ] Fi c u h 2 K 1 2 p h 2 q which implies that (16 a(u h, z h +b(p h, z h 1 2 p h 2 q c 2 u h 2 K. Defining v h = u h +(c w h + 1 4c 2 z h, q h = p h and respecting (14-(16 yields a(u h, v h +b(p h, v h b(q h, u h c u h,p h 2 using (3. Proof of Lemma 5.3. Observe that by (9 we get w h 2 c u h 2 and z h 2 c p h 2 q, and therefore using the definition of v h implies v h,q h 2 = v h 2 + p h 2 q 4 u h 2 +4 w h 2 +4 z h 2 + p h 2 q c ( u h 2 + p h 2 q=c u h,p h 2. Thus the numerical scheme is inf-sup stable for all three choices Q h = Q 1 h,c, Q h = Q 0 h,d and Q h = Q 1 h,c Q0 h,d. Let us define the following auxilliary norm v,q 2 a = h 1 v 2 K + v 2 K + h 1 2 {v} 2 Fi + h 1 2 { v } 2 F + q 2 K + h 1 2 {q} 2 F for all v H 2 (K and q L 2 0(Ω which will be used for the continuity result. We denote further by A + B the (in general not direct sum of the functional spaces A and B. Then we prove the following continuity result. Proposition 5.4 (Continuity. Let v H 2 (K +CR, v h V bs, q H 1 (K and q h Q h. There exists a constant c>0 independent of h such that a(v, v h +b(p, v h b(q h, v c v,q a v h,q h.
10 10 E. BURMAN AND B. STAMM Proof. Applying the Cauchy-Schwarz inequality, Lemma 3.4 and (2 yields a(v, v h = ( v, v h K ( { v }, [v h ] F ( v 2 K + h 1 2 { v } 2 F v,q a v h,q h, b(q, v h = (q, v h K +({q}, [v h ] F 1 2 ( v h 2 K + h 1 2 [vh ] 2 F ( q 2 K + h 1 2 {q} 2 F 1 2 ( v h 2 K + h 1 2 [vh ] 2 F ( q 2 K + h 1 1 ( 2 {q} 2 2 F v h 2 K + h [vh ] 2 2 F v,q a v h,q h, b(q h, v = ( q h, v K ([q h ], {v} Fi ( h q h 2 K + h 1 1 ( 2 [qh ] 2 2 Fi h 1 v 2 K + h {v} 2 2 Fi c v,q a v h,q h. Assume that Q h = αq 1 h,c βq0 h,d for (α, β {(1, 0, (0, 1, (1, 1} and define in a standard manner the following quantities (17 η u = u i c u, ξ u = u h i c u, and η p = p P α p, ξ p = p h P α p, where the projection P α : L 2 (Ω Q h is defined by P α = (1 α π 0 + αc h. Observe that ξ u V bs and ξ p Q h since the projection P α preserves the property of zero mean. Further P α satisfies the error estimate (18 v P α v K + h (v P α v K ch αγ+1 v αγ+1,k for all v H γ+1 (K, γ {0, 1}. Proposition 5.5 (Approximability. Let η u H 2 (Ω + CR and η p H γ+1 (Ω + Q h, with γ {0, 1} and Q h = αq 1 h,c βq0 h,d for (α, β {(1, 0, (0, 1, (1, 1}, be defined by (17. Then, there exists a constant c>0 independent of h such that η u,η p + η u,η p a ch u 2,K + ch 1+αγ p 1+αγ,K. Proof. This is a direct consequence of the trace inquality and the error estimates (6 and (18. Theorem 5.6 (Convergence in Energy norm. Let u H 2 (Ω, p H γ+1 (Ω, γ {0, 1}, be the exact solution of problem (12 and let u h V bs, p h Q h, with Q h = αq 1 h,c βq0 h,d for (α, β {(1, 0, (0, 1, (1, 1}, be the approximation defined by (13. Then, there exists a constant c>0 independent of h such that u u h + p p h K + p p h q ch u 2,K + ch 1+αγ p 1+αγ,K. Proof. Let us first establish the a priori estimate for the triple norm. Split the error in a standard manner in two parts (19 u u h,p p h η u,η p + ξ u,ξ p. By Proposition 5.5 it follows that η u,η p ch u 2,K + ch 1+αγ p 1+αγ,K.
11 BUBBLE STABILIZED DG METHOD 11 For the second term of the right hand side of (19 observe that applying the infsup condition of Proposition 5.1, the consistency of Lemma 4.2, the continuity of Proposition 5.4 and the approximability of Proposition 5.5 yields a(ξ u, v h +b(ξ p, v h b(q h, ξ u ξ u,ξ p sup 0 (v h,q h V bs Q h v h,q h Therefore we conclude that a(η u, v h +b(η p, v h b(q h, η u = sup 0 (v h,q h V bs Q h v h,q h c η u,η p a ch u 2,K + ch 1+αγ p 1+αγ,K. (20 u u h + p p h q ch u 2,K + ch 1+αγ p 1+αγ,K. In order to prove the error estimate for the quantity p p h K we follow the standard technique, see [7, 11]. Let H 1 0 (Ω = { v H 1 (Ω v Ω = 0 }. Since p p h L 2 0 (Ω there exists v p H 1 0 (Ω such that v p = p p h and v p K c p p h K. Applying integration by parts yields p p h 2 K =(p p h, v p K = ( (p p h, v p K + ([p p h ], {v p } Fi. Further split this equation (21 p p h 2 K = I + II with I = ( (p p h, v p C h v p K + ([p p h ], {v p C h v p } Fi, II = ( (p p h, C h v p K + ([p p h ], {C h v p } Fi and where C h denotes the vectorial Clément interpolation operator. Firstly observe using the Cauchy-Schwarz inequality, the trace inequality and the approximation properties of C h that (22 I c ( h (p p h K + h 1 2 [p ph ] Fi v p K c p p h q p p h K. Secondly using the consistency of Lemma 4.2 implies II = b(p p h, C h v p =a(u u h, C h v p (23 = ( (u u h, C h v p K ([u u h ], { C h v p } F ( (u u h 2 K + h 1 2 [u uh ] 2 F c ( (u u h K + h 1 2 [u uh ] F C h v p K c ( (u u h K + h 1 2 [u uh ] F p p h K 1 2 ( C h v p 2 K + h 1 2 { Ch v p } 2 F using that C h v p H0 1 (Ω, the Cauchy-Schwarz inequality, the trace inquality, the H 1 (Ω-stability of C h and the stability of v p. Combining (20-(23 leads to the result. Remark 5.7 (Optimal L 2 -convergence. Optimal convergence in the L 2 -norm can be shown using Nitsche s trick. The details are left to the reader. 1 2
12 12 E. BURMAN AND B. STAMM Proposition 5.8. Let (u h,p h V bs Q h be the solution of (13 with Q h = αq 1 h,c βq0 h,d for (α, β {(1, 0, (0, 1, (1, 1}. Then, there holds (24 h f + u h p h K + (1 β h 1 2 [ uh ] Fi + h 1 2 [uh ] F ch f π 0 f K. Additionally, if f H 1 (K there holds (25 h f + u h p h K + (1 β h 1 2 [ uh ] Fi + h 1 2 [uh ] F ch 2 f K. Proof. Let w h be the function defined by Lemma 3.8 with a h = h 2 ( p h π 0 f u h, b h = h 1 [u h ] v and c h = (1 β h[ u h ]. By its properties, integration by parts, Lemma 3.1, Corollary 3.3 and (2 we get a(u h, w h = ( u h, w h K + ([ u h ], {w h } Fi ([u h ] v, { w h } v F = ( u h, π 0 w h K + ([ u h ], {w h } Fi ([u h ] v, { w h } v F = ( u h, a h K + (1 β h 1 2 [ uh ] 2 F i + h 1 2 [uh ] 2 F. On the other hand since p h = αp h,c + βp h,d, using again integration by parts and the properties of the projections we get b(p h, w h = ( p h, w h K ([p h ], {w h } Fi =( p h, π 0 w h K β([p d,h ], {w h } Fi = ( p h, a h K β(1 β([p d,h ], h[ u h ] Fi =( p h, a h K since β(1 β = 0 for the considered set of values of β. Moreover for this set we have that (1 β = (1 β 2 and since a(u h, w h +b(p h, w h = (f, w h K we get h 2 π 0 f + u h p h 2 K + (1 β2 h 1 2 [ uh ] 2 F i + h 1 2 [uh ] 2 F = a(u h, w h +b(p h, w h (π 0 f, w h K =(f π 0 f, w h K f π 0 f K w h K. Using further the stability estimate (9 of w h yields w h K ch (h 2 π 0 f + u h p h 2 K + (1 β2 h 1 2 [uh ] 2 F + h [ uh ] 2 2 Fi and therefore we get h π 0 f + u h p h K + (1 β h 1 2 [ uh ] Fi + h 1 2 [uh ] F ch f π 0 f K. Finally we observe that f + u h p h K f π 0 f K + π 0 f + u h p h K which yields (24 and additionally noting that for f H 1 (K there holds f π 0 f K ch f K proves (25. Corollary 5.9. If f is piecewise constant, i.e. f V 0 h, then [u h ] F =0, [ u h ] Fi =0, if α =1,β =0, [u h ] F =0, u u h K =0, if α =0,β =1, [u h ] F =0, if α =1,β =1. Corollary If f is piecewise constant, i.e. f Vh 0, and Q h = Q 1 h,c then, the discrete solution u h and p h of (13 satisfies the following local linear momentum conservation property u h n κ ds + p h n κ ds = f dx κ κ κ for all κ K and where n κ denotes the outer normal of κ.
13 BUBBLE STABILIZED DG METHOD 13 Q h h u u h K u u h p p h K (u u h K Q 1 h,c E E E E E-04 ( E-03 ( E-04 ( E-03 ( E-05 ( E-03 ( E-04 ( E-03 ( E-06 ( E-03 ( E-05 ( E-03 (0.98 Q 0 h,d E E E E E-04 ( E-02 ( E-02 ( E-02 ( E-04 ( E-02 ( E-02 ( E-03 ( E-05 ( E-03 ( E-03 ( E-03 (1.00 Q 1 h,c Q0 h,d E E E E E-05 ( E-03 ( E-03 ( E-03 ( E-05 ( E-03 ( E-03 ( E-03 ( E-06 ( E-03 ( E-04 ( E-04 (1.00 Table 1. Smooth problem: Different error quantities of the numerical solution for all three choices of the pressure approximation space with respect to the mesh size h. The quantities in the brackets correspond to the convergence rates. Q h h h 1 2 [u h ] F h 1 2 [[u h ]] F h 1 2 [[ u h ]] Fi Q 1 h,c E E E E-03 ( E-04 ( E-04 ( E-03 ( E-05 ( E-05 ( E-03 ( E-05 ( E-06 (2.94 Q 0 h,d E E E E-02 ( E-04 ( E-02 ( E-02 ( E-05 ( E-02 ( E-03 ( E-05 ( E-02 (0.99 Q 1 h,c Q0 h,d E E E E-03 ( E-04 ( E-03 ( E-03 ( E-05 ( E-03 ( E-03 ( E-05 ( E-03 (0.98 Table 2. Smooth problem: Different error quantities of the numerical solution for all three choices of the pressure approximation space with respect to the mesh size h. The quantities in the brackets correspond to the convergence rates. 6. Numerical tests Let us introduce the numerical examples tested in this section. We consider two numerical tests proposed in [11]. i Problem with smooth solution Consider problem (12 with Ω = (0, 1 2 and f(x imposed such that the exact solution is given by ( (x 4 u(x = 1 2x x2 1 (4x3 2 6x2 2 +2x 2 (4x 3 1 6x2 1 +2x 1(x 4 2 2x3 2 + x2 2, p(x =x 1 + x 2 1. Observe that (u,p satisfies the regularity assumption and thus Theorem 5.6 and Proposition 5.8 are valid. A sequence of structured meshes is considered.
14 14 E. BURMAN AND B. STAMM (0,1 x 2 (0,0.5 Ω (2,0 (12,0 x 1 Figure 1. Geometry of the backstep channel problem. ii Backstep channel problem Consider problem (12 with f = 0 and (8(1 x 2 (x 2 0.5, 0 if x 1 =0, g(x = ((1 x 2 x 2, 0 if x 1 = 12, 0 otherwise, on the domain defined by Figure 1. The domain is not convex and thus the solution does not lie in H 2 (Ω. In consequence Theorem 5.6 is no longer valid. But observe that Proposition 5.8 is independent of the geometry resp. regularity assumption and therefore Proposition 5.8 still holds. A sequence of unstructured and globally uniform meshes is considered. We consider the approximations defined by (13 using a pressure approximation space Q h = αq 1 h,c βq0 h,d for (α, β {(1, 0, (0, 1, (1, 1}. For the computations we use FreeFem++ [10] Smooth problem. The convergence results for test problem i with all three choices of the pressure approximation space are illustrated in Table 1 and 2. We observe the optimal convergence as predicted by Theorem 5.6 and the superconvergence of Proposition Backstep channel problem. The convergence results for test problem ii with all three choices of the pressure approximation space is illustrated in Table 3. Since the exact solution is not knows only the non-conforming error-quantities are given. Observe that Proposition 5.8 is still valid, in contrast to Theorem 5.6 since Ω is non-convex, and since f = 0 the quantity h 1 2 [u h ] F is for any choice of pressure space equal to zero (machine precision. Acknowledgements The authors acknowledge the financial support provided through the Swiss National Science Foundation under grant References [1] F. Brezzi and L. D. Marini. Bubble stabilization of discontinuous Galerkin methods. In W. Fitzgibbon, R. Hoppe, J. Periaux, O. Pironneau, and Y. Vassilevski, editors, Advances in numerical mathematics, Proc. International Conference on the occasion of the 60th birthday of Y.A. Kuznetsov, pages Institute of Numerical Mathematics of The Russian Academy of Sciences, Moscow, [2] E. Burman and B. Stamm. Low order discontinuous Galerkin methods for second order elliptic problems. SIAM J. Numer. Anal., 47(1: , [3] E. Burman and B. Stamm. Symmetric and non-symmetric discontinuous Galerkin methods stabilized using bubble enrichment. C. R. Math. Acad. Sci. Paris, 346(1-2: , 2008.
15 BUBBLE STABILIZED DG METHOD 15 Q h h (u u h K h 2 1 [u h ] F h 2 1 [[u h ]] F h 1 2 [[ u h ]] Fi Q 1 h,c E E E E-01 ( E-01 ( E-01 ( E-01 ( E-01 ( E-01 ( E-01 ( E-01 ( E-01 ( Q 0 h,d E E E E-01 ( E-01 ( E+00 ( E-01 ( E-01 ( E+00 ( E-01 ( E-01 ( E-01 (0.87 Q 1 h,c Q0 h,d E E E E E-01 ( E-01 ( E-01 ( E+00 ( E-01 ( E-01 ( E-01 ( E-01 ( E-02 ( E-02 ( E-01 ( E-01 (0.86 Table 3. Backstep channel problem: Different error quantities of the numerical solution for all three choices of the pressure approximation space with respect to the mesh size h. The quantities in the brackets correspond to the convergence rates and 0 corresponds to zero in machine precision. [4] P. G. Ciarlet. The finite element method for elliptic problems, volume 40 of Classics in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM, Philadelphia, PA, Reprint of the 1978 original [North-Holland, Amsterdam; MR (58 #25001]. [5] P. Clément. Approximation by finite element functions using local regularization. Rev. Française Automat. Informat. Recherche Opérationnelle Sér. RAIRO Analyse Numérique, 9(R-2:77 84, [6] B. Cockburn, G. Kanschat, D. Schötzau, and C. Schwab. Local discontinuous Galerkin methods for the Stokes system. SIAM J. Numer. Anal., 40(1: (electronic, [7] V. Girault and P.-A. Raviart. Finite element methods for Navier-Stokes equations, volume 5 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin, Theory and algorithms. [8] V. Girault, B. Rivière, and M. F. Wheeler. A discontinuous Galerkin method with nonoverlapping domain decomposition for the Stokes and Navier-Stokes problems. Math. Comp., 74(249:53 84 (electronic, [9] P. Hansbo and M. G. Larson. Discontinuous Galerkin methods for incompressible and nearly incompressible elasticity by Nitsche s method. Comput. Methods Appl. Mech. Engrg., 191(17-18: , [10] F. Hecht, O. Pironneau, A. L. Hyaric, and K. Ohtsuka. FreeFEM++ Manual, [11] B. Rivière and V. Girault. Discontinuous finite element methods for incompressible flows on subdomains with non-matching interfaces. Comput. Methods Appl. Mech. Engrg., 195(25-28: , [12] A. Toselli. hp discontinuous Galerkin approximations for the Stokes problem. Math. Models Methods Appl. Sci., 12(11: , Erik Burman, Department of Mathematics, Mantell Building, University of Sussex, Falmer, Brighton, BN1 9RF, United Kingdom address: E.N.Burman@sussex.ac.uk Benjamin Stamm, EPFL, IACS-CMCS, Station 8, 1015 Lausanne, Switzerland address: benjamin.stamm@epfl.ch
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS ERIK BURMAN AND BENJAMIN STAMM Abstract. In this paper we give an analysis of a bubble stabilized discontinuous Galerkin
More informationBubble stabilized discontinuous Galerkin methods on conforming and non-conforming meshes
Calcolo 2011) 48: 189 209 DOI 10.1007/s10092-010-0036-3 Bubble stabilized discontinuous Galerkin methods on conforming and non-conforming meshes Erik Burman Benjamin Stamm Received: 4 January 2010 / Accepted:
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
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 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 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 information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
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 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 informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA PENALTY FREE NON-SYMMETRIC NITSCHE TYPE METHOD FOR THE WEAK IMPOSITION OF BOUNDARY CONDITIONS
A PENALTY FREE NON-SYMMETRIC NITSCHE TYPE METHOD FOR THE WEAK IMPOSITION OF BOUNDARY CONDITIONS ERIK BURMAN Abstract. In this note we show that the non-symmetric version of the classical Nitsche s method
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
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 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 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 informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
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 informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
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 informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
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 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 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 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 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 informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
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 informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
More informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More informationarxiv: v1 [math.na] 13 Sep 2016
THE PENALTY FREE NITSCHE METHOD AND NONCONFORMING FINITE ELEMENTS FOR THE SIGNORINI PROBLEM ERIK BURMAN, PETER HANSBO, AND MATS G. LARSON arxiv:1609.03745v1 [math.na] 13 Sep 2016 Abstract. We design and
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 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 informationA Stabilized Finite Element Method for the Darcy Problem on Surfaces
arxiv:1511.03747v1 [math.na] 12 Nov 2015 A Stabilized Finite Element Method for the Darcy Problem on Surfaces Peter Hansbo Mats G. Larson October 8, 2018 Abstract We consider a stabilized finite element
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
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 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 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 informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
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 informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
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 informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
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 informationMixed Discontinuous Galerkin Methods for Darcy Flow
Journal of Scientific Computing, Volumes 22 and 23, June 2005 ( 2005) DOI: 10.1007/s10915-004-4150-8 Mixed Discontinuous Galerkin Methods for Darcy Flow F. Brezzi, 1,2 T. J. R. Hughes, 3 L. D. Marini,
More informationAn hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity
An hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity Paul Houston 1 Department of Mathematics, University of Leicester, Leicester LE1 7RH, UK (email: Paul.Houston@mcs.le.ac.uk)
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 informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationAn a posteriori error indicator for discontinuous Galerkin discretizations of H(curl) elliptic partial differential equations
IMA Journal of Numerical Analysis 005) Page of 7 doi: 0.093/imanum/ An a posteriori error indicator for discontinuous Galerkin discretizations of Hcurl) elliptic partial differential equations PAUL HOUSTON
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
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 informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
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 informationSTABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW. Math. Models Methods Appl. Sci., Vol. 13, pp , 2003
Mathematical Models and Methods in Applied Sciences c World Scientiic Publishing Company STABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW D. SCHÖTZAU1 C. SCHWAB 2 A. TOSELLI 2 1 Department o Mathematics, University
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationEnergy Norm A Posteriori Error Estimation for Mixed Discontinuous Galerkin Approximations of the Stokes Problem
idgenössische Technische Hochschule Zürich cole polytechnique fédérale de Zurich Politecnico federale di Zurigo Swiss Federal Institute of Technology Zurich nergy Norm A Posteriori rror stimation for Mixed
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
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 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 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 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 information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationDISCONTINUOUS GALERKIN FINITE ELEMENT DISCRETIZATION FOR STEADY STOKES FLOWS WITH THRESHOLD SLIP BOUNDARY CONDITION
DISCONTINUOUS GALRKIN FINIT LMNT DISCRTIZATION FOR STADY STOKS FLOWS WITH THRSHOLD SLIP BOUNDARY CONDITION J.K. Djoko Department of Mathematics and Applied Mathematics, University of Pretoria, Private
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
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 informationA Petrov-Galerkin Enriched Method: A Mass Conservative Finite Element Method For The Darcy Equation
University of Colorado at Denver and Health Sciences Center A Petrov-Galerkin Enriched Method: A Mass Conservative Finite Element Method For The Darcy Equation Gabriel R. Barrenechea, Leopoldo P. Franca,
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
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 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 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 informationStabilized Finite Element Approximation of the Mean Curvature Vector on Closed Surfaces
arxiv:1407.3043v1 [math.na] 11 Jul 2014 Stabilized Finite Element Approximation of the Mean Curvature Vector on Closed Surfaces Peter Hansbo Mats G. Larson Sara Zahedi October 17, 2018 Abstract We develop
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationA MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp. 2078 2095 c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG
More informationAn Equal-Order DG Method for the Incompressible Navier-Stokes Equations
J Sci Comput 009) 40: 188 10 DOI 10.1007/s10915-008-961-1 An qual-order DG Method for the Incompressible Navier-Stokes quations Bernardo Cockburn Guido anschat Dominik Schötzau Received: 4 September 008
More information