A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne France-Comté, 16 route de Gray, 5030 Besançon Cedex, France. e-mail: alexei.lozinski@univ-fcomte.fr Marc 0, 018 Abstract We propose an efficient variant of a primal Discontinuous Galerkin metod wit interior penalty for te second order elliptic equations on very general meses (polytopes wit eventually curved boundaries). Efficiency, especially wen iger order polynomials are used, is acieved by static condensation, i.e. a local elimination of certain degrees of freedom element by element. is alters te original metod in a way tat preserves te optimal error estimates. Numerical experiments confirm tat te solutions produced by te new metod are indeed very close to tat produced by te classical one. Keywords: Discontinuous Galerkin, static condensation, polyedral (polygonal) meses. 1 Introduction e recent years ave seen te emergence (or te revival) of several numerical metods capable to solve approximately elliptic partial differential equations using general polygonal/polyedral meses. is is witnessed for example by te book [4]. e metods reviewed in tis book (if we restrict our attention only to finite element type metods using piecewise polynomial approximation spaces in one form or anoter) include interior penalty discontinuous Galerkin (DG) metods [], ybridizable discontinuous Galerkin (HDG) metods ([6], introduced in [8]), te Virtual Element (VE) metod ([18], introduced in [17, 5]), te Hybrid Hig-Order (HHO) metod ([1], introduced in [10, 11]). One can add to tis list te weak Galerkin finite element [19], wic is similar to HDG. e relations between HHO and HDG metods were exibited in [7]. Among te above, te primal interior penalty DG metods are te most classical. In te symmetric form, also referred to as SIP symmetric interior penalty, tis metod dates back to [0, 3] and is now presented and torougly studied in several monograps, for example [16, 9]. It is well suited to te discretization on very general meses because its approximation space is populated by polynomials of degree, say k, on eac mes element witout any constraints linking te polynomials on two adjacent elements. It leaves tus a lot of freedom on te coice of te mes elements wic can be not only polytopes but also virtually any geometrical sapes. It is generally admitted owever tat te SIP metod is too expensive especially wen iger order polynomials are employed. Indeed, its cost, i.e. te dimension of te approximation space, is te product of te number of mes elements and te dimension of te space of polynomials of degree k. e cost on a given mes is tus proportional to k in D (resp. k 3 in 3D). is sould be contrasted wit te cost of HDG or HHO metods wic is proportional to k in D (resp. k in 3D). e goal of te present article is to modify te SIP metod so tat its cost is reduced to tat of HDG or HHO metods. In doing so, we inspire ourselves from te static condensation procedure for te standard continuous Galerkin (CG) finite element metods. It is indeed well known tat te dimension of te approximation space in CG is proportional to k in D on a given mes, but te degrees of freedom interior 1
to eac mes element can be locally eliminated wic leaves a global problem of te size proportional to k (tese numbers are canged to, respectively, k 3 and k in 3D). Altoug te notion of interior degrees of freedom does not make sense in te DG context, we sall be able to select, on eac mes element, a subspace of te approximating polynomials tat can be used to construct a local approximation troug te solution of a local problem. e remaining degrees of freedom will ten be used in a global problem. We sall tus acieve a significant reduction of te global problem size in te DG SIP-like metod, similarly to tat acieved in CG by te static condensation. e resulting DG metod, wic can be refereed to as scsip (static condensation SIP), will not produce exactly te same approximation as te original SIP metod. We sall prove owever tat tese two solutions satisfy te same optimal a priori error bounds in H 1 and L norms. Moreover, tey turn out to be very close to eac oter in our numerical experiments. We treat ere only te diffusion equation wit variable, but sufficiently smoot, coefficients. e extension to oter problems, suc as convection-reaction-diffusion, linear elasticity, Stokes, as well as to oter DG variants (IIP, NIP) seems relatively straigt-forward. Our assumptions on te mes allow for elements of general sape, not necessarily te polytopes. e article is organized as follows: in te next section, we present te idea of our metod starting by te description of te governing equations. We te recall te static condensation for te classical CG FEM. Our variants of DG FEM (SIP and scsip) are first introduced in Subsection.. e convergence proofs are in Section 3. ey are done assuming some properties of te discontinuous FE spaces and te underlying mes. In Section 4, we give an example of te ypoteses on te mes under wic te necessary properties of te FE spaces can be establised. Finally, a numerical illustration is presented in Section 5. Description of te problem and static condensation for FEM (CG and DG cases) We consider te second-order elliptic problem L u = f in Ω, u = g on Ω (1) were Ω R d, d = or 3, is a bounded Lipscitz domain, f L (Ω) and g H 1/ ( Ω) are given functions. e differential operator L is defined by L u = i (A i j (x) j u) wit i denoting te partial derivative in te direction x i, i = 1,...,d. e coefficients A i j are supposed to form a positive definite matrix A = A(x) for any x Ω wic is sufficiently smoot wit respect to x so tat and wit some constants β α > 0, M > 0..1 Static condensation for CG FEM α ξ ξ A(x)ξ β ξ, ξ R d,x Ω () A i j (x) M, x Ω,i, j = 1,...,d (3) o present our idea, we start by recalling te idea of static condensation, going back to [14], as applied to te usual CG finite element metod for problem (1). Let us assume for te moment (in tis subsection only) tat Ω is a polygon (polyedron) and introduce a conforming mes on Ω consisting of triangles (tetraedrons). Assuming for simplicity g = 0, te usual continuous P k finite element discretization of (1) is ten written: find u V suc tat a(u,v ) := A u v = f v, v V (4) Ω Ω
were V is te space of continuous piecewise polynomial functions (polynomials of degree k on eac mes element for some k 1) vanising on Ω. e size of tis problem, i.e. te dimension of V, is of order k on a given mes in D (resp. k 3 in 3D). o reduce tis cost, one can decompose te space V as follows V = V loc V were te subspace V loc consists of functions of V tat vanis on te boundaries of all te mes elements, and V is te complement of V loc, ortogonal wit respect to te bilinear form a. Decomposing u = u loc + u wit uloc V loc and u V we see tat (4) is split into two independent problems u loc V loc : a(u loc,vloc ) = f v loc, vloc V loc (5) Ω u V : a(u,v ) = f v, v V (6) e first problem above is furter split into a collection of mutually independent problems on every mes element : Find u loc, suc tat A u loc, v loc, = f v loc,, v loc, V loc, (7) Ω were V loc, is te restriction of V loc on, i.e. te set of all polynomials of degree k vanising on. e cost of solution of tese local problems is negligible and we tus get very ceaply u loc = u loc,. Note also tat Problem (7) can be recast as Ω π L (u loc ) = π f, (8) were π is te projection to V loc,, ortogonal in L ( ). On te oter and, Problem (6) remains global but its size is only proportional to k in D (resp. k in 3D) wic is muc smaller tan tat of te original problem (4). Indeed, te degrees of freedom are associated to te standard interpolation points of P k finite elements on te edges of te mes. Note also tat a basis for V can be constructed solving ceap local problems of te type π L (v ) = 0, (9) wit appropriate boundary conditions on insuring te continuity of functions in V.. DG FEM: SIP and scsip metods We turn now to te main subject of tis paper: te DG metods. We now let Ω be a bounded domain of general sape, and be a splitting of Ω into a collection of non-overlapping subdomains (again of general sape, te precise definitions and assumptions on te mes will be given Sections 3 and 4). Let V denote te space of discontinuous piecewise polynomial functions (polynomials of degree k on eac mes element for some k 1). e SIP metod is ten written as: find u V suc tat a (u,v ) = L (v ), v V (10) wit te bilinear form a and te linear form L defined by a (u,v) = A u v ({A u n}[v] + {A v n}[u]) E E E + E E γ E E [u][v] (11) and ( ) γ L (v) = f v + g v A v n E E b E E (1) 3
were E is te set of all te edges/faces of te mes, E b E regroups te edges/faces on te boundary Ω, n, [ ] and { } denote te unit normal, te jump and te mean over E E. More precisely, for any internal edge E sared by two mes elements 1 E and E, we coose n as te unit vector, normal to E and looking from 1 E to E. We ten define for any function v wic is H1 on bot 1 E and E but discontinuous across E [v] E := v E v 1 E, {A v n} E := 1 ( ) A v E + v 1 E n On a boundary edge E E b, n is te unit normal looking outward Ω and [v] = v, {A v n} = A v n. Unlike te case of continuous finite elements, Problem (10) does not allow directly for a static condensation. However, we can construct a modification of (10) tat mimics te caracterization of local and global components of te solution by te projectors on local polynomial spaces (8) (9). ese spaces are now defined simply as V loc, = P k We also let π,k to be te projection to V loc,, ortogonal in L ( ), and propose te following DG sceme: Compute v loc V by solving π,k L (u loc ) = π,k f, (13) We sall sow tat tis problem admits an infinity of solutions. We coose arbitrarily a solution on eac mes element. Define te subspace of V V = { v V : π,k L (v ) = 0, } (14) and compute u V suc tat a (u,v ) = L (v ) a (u loc,v ), v V (15) Set u = u loc + u. (16) e local projection step (13) is not necessarily consistent wit te original formulation (10) so tat te solution u given by (13) (16) is different from tat of (10). We sall prove owever tat SIP and scsip approximations satisfy te same a priori error bounds. Moreover, tey turn out to be very close to eac oter in numerical experiments. Note tat te dimension of te global space V on a given mes is of order k in D (k in 3D) so tat global problem (15) is muc ceaper tan (10) for large k. 3 Well posedness of SIP and scsip metods and a priori error estimates Let us now be more precise about te ypoteses on te mes. Recall tat Ω R d, d = or 3, is a Lipscitz bounded domain and is a general (not necessarily polygonal or polyedral) mes on Ω. We mean by tis tat is a decomposition of Ω into mutually disjoint cells Ω = so tat eac cell is a Lipscitz subdomain of Ω and for every 1, we ave eiter 1 = or 1 = (te cells 1 are 4
treated ere as open sets). We also introduce te sets of internal and boundary edges/faces as respectively E i = {E = 1 1 for some 1, } E b = {E = Ω for some } and denote by E := E i E b te union of all te edges/faces. Let B, for any, denote te smallest ball containing, and B in denote te largest ball inscribed in. Set = diam( ) and = max. From now on, we assume tat mes is Sape regular: tere is a mes-independent parameter ρ 1 > 1 suc tat, for all, R ρ 1 r (17) were r is te radius of B in and R is te radius of B. is also implies ρ 1 r and R ρ 1. Locally quasi-uniform: tere is a mes-independent parameter ρ > 1 suc tat for all 1, wit one as ρ 1 (18) Coose an integer k 1 and introduce te FE space V = {v L (Ω) : v P k ( ), } (19) In te sequel, we assume tat V as two following properties (and we sall prove in Section 4 tat tese properties old under some additional assumptions on te mes): Optimal interpolation: tere exists an operator I : H k+1 (Ω) V suc tat for any v H k+1 (Ω) ( ( v I v H 1 ( ) + 1 v I v L ( ) + v I v H ( ) (0) + v I v L ( ) + 1 )) 1 ( v I v L ( ) C k v H k+1 ( ) Inverse inequalities: for any v V and any v L ( ) C v L ( ) v L ( ) C v L ( ) (1) ) 1 and v H ( ) C v H 1 ( ) () We can now study te well posedness and establis optimal a priori error estimates for te classical DG metod (10). ese results are well known but we sall sketc teir proofs for te sake of completeness. Lemma 1. Under te above assumptions on te mes and on V, setting ( 1 E = + 1 ) 1 for any E E i 1 wit E = 1, (3) E = for any E E b wit E = Ω, and coosing γ large enoug, γ γ 0, te bilinear form a defined by (11) is coercive, i.e. wit some c > 0 and te triple norm defined by a (v,v ) c v, v V (4) u = u H 1 ( ) + 1 [u] E E L (E) (5) E e constants c,γ 0 depend only on te parameters in te assumptions on te mes and on V, as well as on α,β in (). 5
Proof. A direct calculation sows a (v,v ) = A v v {A v n}[v ] + E E E ( α v H 1 ( ) A v n[v ] + Ω ( γ + [v ] + γ ) v Ω Ω α v H 1 ( ) + γ E E E Ω γ C β /(α) [v ] L ( ) E [v ] (A v n)v ) wit C te constant from te inverse inequality (1). Rewriting te sum over te boundaries of mes cells back to te sum over te edges and recalling te local quasi-uniformity assumption we get (4). Lemma 1 implies tat problem (10) of te SIP metod is well posed. Moreover, we ave te error estimate: eorem 1. Assume tat te solution u to (1) is in H k+1 (Ω). Under te assumptions of Lemma 1, tere exists te unique solution u to (10) and it satisfies u u H 1 ( ) Ck u H k+1 (Ω) were H 1 ( ) is te broken H 1 space on te mes and H 1 ( ) := ( H 1 ( )) 1. If, moreover, te elliptic regularity property olds for (1), ten u u L (Ω) C u H k+1 (Ω) k+1 We skip te proof of tis well known result. Actually, it goes along te same lines as tat of our fortcoming eorem. We turn now to te study of te scsip metod (13) (15) and start by te following tecnical lemma. Lemma. ere exists 0 > 0 suc tat for all wit 0 and for all q P k ( ) one can find u P k ( ) suc tat q (L u ) 1 q L ( ) (6) and u H 1 ( ) + 1 u L ( ) C q L ( ) (7) e constants 0 and C depend only on te regularity of te mes and on α, β and M in () and (3). Proof. Let χ be te polynomial of degree vanising on B in, i.e. ) χ (x) = ( d (x i xi 0 ) r i=1 were x 0 = (x 0 1,...,x0 d ) is te center of Bin and r is its radius. Set A 0 i j = A i j(x 0 ) and L 0 = i A 0 i j j. Consider te linear map Q : P k ( ) P k ( ) defined by Q(v) = L 0 (χ v) e kernel of Q is {0}. Indeed, if Q(v) = 0 ten w := χ v is te solution to L 0 w = 0 in B in, w = 0 on B in 6
so tat w = 0 as a solution to an elliptic problem wit vanising rigt-and side and boundary conditions. Since Q is a linear map on te finite dimensional space P k ( ), tis means tat Q is one-to-one. ake any q P k ( ) and let u = χ v wit v P k ( ) suc tat Q(v ) = q. We ave tus constructed u P k ( ) suc tat L 0 u = q. is immediately proves (6) in te case of an operator L = L 0 wit constant coefficients. Moreover, by scaling, u W, (B ) + 1 u W 1, (B ) + 1 u L (B ) wit a constant C depending only on α, β and te ratio R /r. us, C d/ u H 1 ( ) 1/ u W 1, (B ) C q L ( ) q L (B in ) (8) wic proves te estimate in H 1 ( ) norm in (7). Similarly, u L ( ) C q L ( ) and te estimate in L ( ) norm in (7) follows by te trace inverse inequality. It remains to prove (6) in te case of operator L wit variable coefficients. o tis end, we use te estimates in (8) as follows q L u = q L 0 u + q i ((A i j A 0 i j) j u ) q L ( ) q L ( ) 1/ [max A(x) A 0 u W x, ( ) + max A(x) u W x 1, ( ) ] for sufficiently small. q L ( ) q L ( ) 1/ max x A(x) C q L ( ) C q L ( ) 1 q L ( ) Corollary 1. Introduce te bilinear form and te space b (q,v) = d/ ql u q L (B in ) M = {v L (Ω) : v P k ( ), } Equip te space V wit te triple norm (5) and te space M wit ( q = q L ( ) e bilinear form b satisfies te inf-sup condition ) 1/ inf q M sup v V b(q,v ) q v δ (9) wit a mes-independent constant δ > 0. Moreover, b is continuous on M V wit a mes-independent continuity bound. Proof. ake any q M, denote q = q, construct u as in Lemma and introduce u V by u = u on all. is yields using (6) and (7), b (q,v ) v q L ( ( ) 1/ = q C q C L ( ) wic is equivalent to (9) wit δ = / C. Finally, te continuity of b is easily seen from inverse inequality (). 7
Lemma implies tat operator π,k appearing in (13) is surjective from P k to P k so tat (13) as indeed a solution at least on sufficiently refined meses. e existence of a solution to (15) follows from te coercivity of a. us, sceme (13) (16) produced indeed a solution u V. In order to establis te error estimates for tis solution, we reinterpret its definition as a saddle point problem. Lemma 3. e problem of finding u V and p M suc tat a (u,v ) + b (p,v ) = L (v ), v V (30) b (q,u ) = f q, Ω q M (31) as te unique solution. Moreover, u given by (30) (31) coincides wit u given by (13) (16). is implies tat u produced by te sceme (13) (16) is unique. Proof. e existence and uniqueness of te solution to (30) (31) follows from te standard teory of saddle point problems, cf. for example Proposition.36 from [13], tanks to te coercivity of a (Lemma 1) and to te inf-sup property on b (Corollary 1). In order to explore its relation wit u from (13) (16), we note b (q,u loc ) = Ω f q for all q M and b (q,v ) = 0 for all q M, v V. We obtain tus b (q,u loc + u ) = f q, q M Ω wic is equivalent to (31). Eq. (15) can be rewritten as a (u loc + u,v ) = L (v ), q uad v V is is equivalent to (30) after eliminating p since V is precisely te kernel of te bilinear form b. We conclude tat u = u loc + u coincides wit te solution to (30) (31).. eorem. Assume tat te solution u to (1) is in H k+1 (Ω). Under te assumptions of Lemma 1 and sufficiently small, te scsip metod (13) (16) produces te unique solution u V, wic satisfies If, moreover, te elliptic regularity property olds for (1), ten u u H 1 ( ) Ck u H k+1 (Ω) (3) u u L (Ω) C u H k+1 (Ω) k+1 (33) Proof. We sall use te saddle point reformulation (30) (31). is discretization is consistent. Indeed setting p = 0 we ave a (u,v ) + b (p,v ) = L (v ), b (q,u) = f q, v V q M us, by te standard teory of saddle point problems, cf. for example Proposition.36 from [13], recalling te coercivity of a (Lemma 1) and te inf-sup property on b (Corollary 1), we get u I u + p C sup (v,q ) V M v + q = 1 = C sup (v,q ) V M v + q = 1 ( C (a (u I u,v ) + b (p,v ) + b (q,u I u)) (a (u I u,v ) + b (q,u I u)) u I u a + L (u I u) L ( ) ) 1 8
wit te augmented triple norm a defined by v a := v + E E E {A v n} L (E) Applying te interpolation estimates (0) and te triangle inequality gives u I u a + p C k u H k+1 (34) is implies in particular (3). o prove te L error estimate, we consider te auxiliary problem for z H (Ω) L z = u u in Ω, z = 0 on Ω en, for all v H 1 ( ) and q L (Ω), a (v,z) + b (q,z) = (u u )v + (u u )q Ω Setting v = u u and q = p p (wit p = 0) and using Galerkin ortogonality yields u u L (Ω) + (u u )(p p ) = a (u u,z) + b (p p,z) = a (u u,z z ) + b (p p,z z ) for any z V. us, taking z = I z and applying te interpolation estimates, ( ) 1/ u u L (Ω) a (u u,z z ) + p p L z L z L ( ) + p p u u L (Ω) C u I u a z H (Ω) +C p ( z H (Ω) + u u L (Ω) ) Recalling z H (Ω) C u u L (Ω) and (34) yields (33). 4 An example of assumptions on te mes tat guarantee te interpolation and inverse estimates In tis section, we adopt te following assumptions on te mes. M1: is sape regular in te sense (17) wit a parameter ρ 1 > 1. M: is locally quasi-uniform in a sense stronger tan (18): one can assign a local messize x > 0 to any point x Ω so tat for all and any x B wit a parameter ρ > 1. 1 ρ x ρ x M3: e cell boundaries are not too wiggly: for all wit a parameter ρ 3 > 0. ρ 3 d 1 9
We sall sow tat tese assumptions allow us to construct an interpolation operator I to te discontinuous finite element space (19) at least for k and to prove te optimal error (0) te inverse estimates (1) (). First of all, te assumptions tat te mes is sape regular and locally quasi-uniform entail te following Lemma 4. Let N ball (x) be te number of balls B, covering a point x, x R d. Under assumptions M1 and M, we ave N ball (x) N int for all x R d wit a constant N int depending only on ρ 1 and ρ. Proof. ake any x R d. If x B for some ten ρ x. Hence, is inside te ball B x of radius ρ x centered at x. Recall tat contains an inscribed ball of radius r R ρ 1 ρ 1 ρ 1 ρ x. If tere are several suc mes cells covering x, ten teir inscribed balls do not intersect eac oter and tey are all inside B x. us, teir number is bounded by B x min x B B in (ρ x ) d ( ρ1 ρ ) d ( ) d = ρ 1 ρ x Recall tat V is te discontinuous FE space on of degree k 1, cf. (19). ake any v H k+1 (Ω). Lemma 5. (Local interpolation estimate) ake any. Let π denote te L (B )-ortogonal projection to te space of polynomials, i.e. given v L (B ), v = π v is a polynomial from P k suc tat v ϕ = vϕ B B ϕ P k Under Assumptions M1 and M3, we ave ten for any v H k+1 (B ) v v H 1 ( ) + 1 v v L ( ) + v v H ( ) wit a constant C > 0 depending only on ρ 3. + (v v ) L ( ) + 1 v v L ( ) Ck v H k+1 (B ) Proof. Since te embedding of H k+1 (B ) into L (B ) is compact, Petree-artar lemma togeter wit a scaling argument entail v v L (B ) C k+1 d/ v H k+1 (B ) Hence, and, in view of te ypotesis ρ 3 d 1, v v L ( ) 1 v v L (B ) C k+1 v H k+1 (B ) v v L ( ) (ρ 3 d 1 ) 1 v v L (B ) C k+1/ v H k+1 (B ) e estimates for v v H 1 ( ) and v v L ( ) are proven in te same way starting from v v L (B ) C k d/ v H k+1 (B ) is is valid since te embedding of H k+1 (B ) into W 1, (B ) is compact for k. Finally, te estimate for v v H ( ) olds by te compactness of te embedding of Hk+1 (B ) into H (B ) (k ). 10
Lemma 6. (Global interpolation estimate) Let I : H k+1 (Ω) V denote te operator obtained by first extending a function v H k+1 (Ω) by a function ṽ H k+1 (R d ) and ten applying te local operator π from Lemma 5 to ṽ on every, i.e. I v := (π ṽ) on any. en, under Auumptions M1 M3, (0) olds for any v H k+1 (Ω) Proof. First, extension teorem for Sobolev spaces [1] insure tat tere exists ṽ H k+1 (R d ) suc tat ṽ = v on Ω and ṽ H k+1 (R d ) C ṽ H k+1 (Ω o prove (0), we sum te local interpolation estimates of Lemma 5 over all te mes elements and ten use te estimates of Lemma 4: ( v v H 1 ( ) + 1 v v L ( ) + v v L ( ) + 1 ) v v L ( ) C k k+1 v dx Cρ k B Ω k x N ball (x) k+1 v dx CN int ρ 4k k v H k+1 ( ) Lemma 7. (Inverse inequalities) Under assumptionsone M1 and M3, (1) and () old for any v V and any. Proof. Bot bounds in (1) follow immediately from te following one: for any polynomial v of degree l (wit l = k or l = k 1) one as wic in turn follows from v L ( ) v L (B ) C d/ C d/ v L ( ) v L (B in ) were B in is te largest ball inscribed in. Scaling te ball B to a ball of radius 1 B 1 and considering all te possible positions of te inscribed ball, te last inequality can be rewritten as v L (B 1 ) C min B in B 1,B in a ball of radius ρ 1 1 is is valid for any polynomial of degree l by equivalence of norms. e remaining inverse inequality () can be proven similarly: v L (B in ) v H ( ) Cd/ v W, (B ) Cd/ 1 v W 1, (B ) C v H 1 (B in ) C v H 1 ( ) 5 Numerical results We sall illustrate te convergence of SIP and scsip metods on polygonal meses obtained by agglomerating te elements of a background triangular mes. Bot te mes construction and te following calculations are done in FreeFEM++ [15]. An example of suc a mes is given in Fig. 1. o construct it, we take a positive integer n (n = 4 in te Figure), let FreeFEM++ to construct a Delaunay triangulation of Ω = (0,1) wit 4n boundary nodes on eac side of te square, and finally agglomerate te triangles of tis mes into n n cells as follows. We start by attributing te triangle containing te point ( i 1/ n, j 1/ n ), i, j = 1,...,n 11
Figure 1: On te left, a polygonal mes consisting of 4 4 cells, wic are obtained by te agglomeration of te triangles of a finer mes seen on te rigt. Figure : e error in L norm and H 1 semi-norm vs. mes-size. e solid lines wit squares represent te SIP metod. e dased lines wit circles represent te scsip metod. 1
to te cell number i + jn. en, iteratively, we run over all te cells and attac yet unattributed triangles neigboring a triangle from a cell to te same cell, until all te triangles are attributed. We ave considered Poisson equation, i.e. (1) wit A = I, on Ω = (0,1) wit omogeneous Diriclet boundary conditions g = 0 and te exact solution u = sin(πx) sin(πy). We ave applied SIP metod (10) and scsip metod (13) (15) to tis problem on te agglomerated meses as described above. e results are presented in Fig.. In a sligt deviation from te general notations, we set ere te mes-size as = 1/n on te n n mes, and E = on all te edges in (11). ree coices for te polynomial space degree k were investigated, namely k =,3,4 and te penalty parameter γ in (11) was set, respectively, to 10,30,50. e numerical results confirm te teoretically expected order of convergence in bot L norm and H 1 semi-norm. ey also demonstrate tat te approximation produced by SIP and scsip metods are very close to eac oter. Acknowledgements I am grateful to Simon Lemaire for interesting discussions about HHO and mshho metods, wic were te starting point of conceiving te present article. References [1] Adams, R.A.: Sobolev spaces. Academic Press [A subsidiary of Harcourt Brace Jovanovic, Publisers], New York-London (1975). Pure and Applied Matematics, Vol. 65 [] Antonietti, P.F., Cangiani, A., Collis, J., Dong, Z., Georgoulis, E.H., Giani, S., Houston, P.: Review of discontinuous galerkin finite element metods for partial differential equations on complicated domains. In [4], pp. 81 310. Springer, [Cam] (016) [3] Arnold, D.N.: An interior penalty finite element metod wit discontinuous elements. SIAM J. Numer. Anal. 19(4), 74 760 (198). [4] Barrenecea, G.R., Brezzi, F., Cangiani, A., Georgoulis, E.H. (eds.): Building bridges: connections and callenges in modern approaces to numerical partial differential equations, Lecture Notes in Computational Science and Engineering, vol. 114. Springer, [Cam] (016). Selected papers from te 101st LMS-EPSRC Symposium eld at Duram University, Duram, July 8 16, 014 [5] Brezzi, F., Marini, L.D.: Virtual element metods for plate bending problems. Comput. Metods Appl. Mec. Engrg. 53, 455 46 (013). [6] Cockburn, B.: Static condensation, ybridization, and te devising of te HDG metods. In [4], pp. 19 177. Springer, [Cam] (016) [7] Cockburn, B., Di Pietro, D.A., Ern, A.: Bridging te ybrid ig-order and ybridizable discontinuous Galerkin metods. ESAIM Mat. Model. Numer. Anal. 50(3), 635 650 (016). [8] Cockburn, B., Gopalakrisnan, J., Lazarov, R.: Unified ybridization of discontinuous Galerkin, mixed, and continuous Galerkin metods for second order elliptic problems. SIAM J. Numer. Anal. 47(), 1319 1365 (009). [9] Di Pietro, D.A., Ern, A.: Matematical aspects of discontinuous Galerkin metods, Matématiques & Applications (Berlin) [Matematics & Applications], vol. 69. Springer, Heidelberg (01). [10] Di Pietro, D.A., Ern, A.: A ybrid ig-order locking-free metod for linear elasticity on general meses. Comput. Metods Appl. Mec. Engrg. 83, 1 1 (015). [11] Di Pietro, D.A., Ern, A., Lemaire, S.: An arbitrary-order and compact-stencil discretization of diffusion on general meses based on local reconstruction operators. Comput. Metods Appl. Mat. 14(4), 461 47 (014). 13
[1] Di Pietro, D.A., Ern, A., Lemaire, S.: A review of ybrid ig-order metods: formulations, computational aspects, comparison wit oter metods. In [4], pp. 05 36. Springer, [Cam] (016) [13] Ern, A., Guermond, J.L.: eory and practice of finite elements, Applied Matematical Sciences, vol. 159. Springer-Verlag, New York (004) [14] Guyan, R.J.: Reduction of stiffness and mass matrices. AIAA journal 3(), 380 (1965) [15] Hect, F.: New development in freefem++. J. Numer. Mat. 0(3-4), 51 65 (01) [16] Riviere, B.: Discontinuous Galerkin metods for solving elliptic and parabolic equations: teory and implementation. SIAM (008) [17] Beirão da Veiga, L., Brezzi, F., Marini, L.D.: Virtual elements for linear elasticity problems. SIAM J. Numer. Anal. 51(), 794 81 (013). [18] Beirão da Veiga, L., Brezzi, F., Marini, L.D., Russo, A.: Virtual element implementation for general elliptic equations. In [4], pp. 39 71. Springer, [Cam] (016) [19] Wang, J., Ye, X.: A weak Galerkin finite element metod for second-order elliptic problems. J. Comput. Appl. Mat. 41, 103 115 (013). [0] Weeler, M.F.: An elliptic collocation-finite element metod wit interior penalties. SIAM Journal on Numerical Analysis 15(1), 15 161 (1978) 14