arxiv: v1 [math.na] 19 Mar 2018
|
|
- Juliet Richards
- 5 years ago
- Views:
Transcription
1 A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv: v1 [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. 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
2 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) Ω Ω
3 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
4 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
5 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
6 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
7 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
8 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
9 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
10 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
11 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
12 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
13 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 Springer, [Cam] (016) [3] Arnold, D.N.: An interior penalty finite element metod wit discontinuous elements. SIAM J. Numer. Anal. 19(4), (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 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, (013). [6] Cockburn, B.: Static condensation, ybridization, and te devising of te HDG metods. In [4], pp 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), (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(), (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), (014). 13
14 [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 Springer, [Cam] (016) [13] Ern, A., Guermond, J.L.: eory and practice of finite elements, Applied Matematical Sciences, vol 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), (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(), (013). [18] Beirão da Veiga, L., Brezzi, F., Marini, L.D., Russo, A.: Virtual element implementation for general elliptic equations. In [4], pp Springer, [Cam] (016) [19] Wang, J., Ye, X.: A weak Galerkin finite element metod for second-order elliptic problems. J. Comput. Appl. Mat. 41, (013). [0] Weeler, M.F.: An elliptic collocation-finite element metod wit interior penalties. SIAM Journal on Numerical Analysis 15(1), (1978) 14
1. 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 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 informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More informationPart VIII, Chapter 39. Fluctuation-based stabilization Model problem
Part VIII, Capter 39 Fluctuation-based stabilization Tis capter presents a unified analysis of recent stabilization tecniques for te standard Galerkin approximation of first-order PDEs using H 1 - conforming
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationarxiv: v1 [math.na] 20 Nov 2018
An HDG Metod for Tangential Boundary Control of Stokes Equations I: Hig Regularity Wei Gong Weiwei Hu Mariano Mateos Jon R. Singler Yangwen Zang arxiv:1811.08522v1 [mat.na] 20 Nov 2018 November 22, 2018
More informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationarxiv: v1 [math.na] 11 May 2018
Nitsce s metod for unilateral contact problems arxiv:1805.04283v1 [mat.na] 11 May 2018 Tom Gustafsson, Rolf Stenberg and Jua Videman May 14, 2018 Abstract We derive optimal a priori and a posteriori error
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationA MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD
A MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD KOLJA BRIX, MARTIN CAMPOS PINTO, AND WOLFGANG DAHMEN Abstract. In tis article we present a multilevel preconditioner for
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
More informationCONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH UNITY INF SUP AND CONTINUITY CONSTANTS FOR PARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICS OF COMPUTATION Volume 82, Number 284, October 203, Pages 923 952 S 0025-578(203)02697-X Article electronically publised on April 23, 203 CONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH
More informationarxiv: v1 [math.na] 6 Dec 2010
MULTILEVEL PRECONDITIONERS FOR DISCONTINUOUS GALERKIN APPROXIMATIONS OF ELLIPTIC PROBLEMS WITH JUMP COEFFICIENTS BLANCA AYUSO DE DIOS, MICHAEL HOLST, YUNRONG ZHU, AND LUDMIL ZIKATANOV arxiv:1012.1287v1
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
More informationarxiv: v1 [math.na] 18 Jul 2015
ENERGY NORM ERROR ESTIMATES FOR FINITE ELEMENT DISCRETIZATION OF PARABOLIC PROBLEMS HERBERT EGGER arxiv:150705183v1 [matna] 18 Jul 2015 Department of Matematics, TU Darmstadt, Germany Abstract We consider
More informationarxiv:math/ v1 [math.ca] 1 Oct 2003
arxiv:mat/0310017v1 [mat.ca] 1 Oct 2003 Cange of Variable for Multi-dimensional Integral 4 Marc 2003 Isidore Fleiscer Abstract Te cange of variable teorem is proved under te sole ypotesis of differentiability
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationA New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions
Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical
More informationHYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE STOKES PROBLEM
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 533 563 S 0025-5718(05)01804-1 Article electronically publised on December 16, 2005 HYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationA Hybridizable Discontinuous Galerkin Method for the Compressible Euler and Navier-Stokes Equations
A Hybridizable Discontinuous Galerkin Metod for te Compressible Euler and Navier-Stokes Equations J. Peraire and N. C. Nguyen Massacusetts Institute of Tecnology, Cambridge, MA 02139, USA B. Cockburn University
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationChemnitz Scientific Computing Preprints
G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Cemnitz Scientific Computing Preprints Impressum:
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationPosteriori Analysis of a Finite Element Discretization for a Penalized Naghdi Shell
International Journal of Difference Equations ISSN 0973-6069, Volume 8, Number 1, pp. 111 124 (2013) ttp://campus.mst.edu/ijde Posteriori Analysis of a Finite Element Discretization for a Penalized Nagdi
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationWeierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stocastik im Forscungsverbund Berlin e.v. Preprint ISSN 0946 8633 Stability of infinite dimensional control problems wit pointwise state constraints Micael
More informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2811 2834 c 2017 Society for Industrial and Applied Matematics GUARANTEED, LOCALLY SPACE-TIME EFFICIENT, AND POLYNOMIAL-DEGREE ROBUST A POSTERIORI ERROR ESTIMATES
More informationH(div) conforming and DG methods for incompressible Euler s equations
H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationKey words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness
PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent
More informationA UNIFORM INF SUP CONDITION WITH APPLICATIONS TO PRECONDITIONING
A UNIFORM INF SUP CONDIION WIH APPLICAIONS O PRECONDIIONING KEN ANDRE MARDAL, JOACHIM SCHÖBERL, AND RAGNAR WINHER Abstract. A uniform inf sup condition related to a parameter dependent Stokes problem is
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationFINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES
FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More informationHigh-Order Extended Finite Element Methods for Solving Interface Problems
Hig-Order Extended Finite Element Metods for Solving Interface Problems arxiv:1604.06171v1 [mat.na] 21 Apr 2016 Fei Wang Yuanming Xiao Jincao Xu ey words. Elliptic interface problems, unfitted mes, extended
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 information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS. 1. Introduction. We consider the parabolic PDE of the form,
ERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS K. CHRYSAFINOS AND NOEL. J. WALKINGTON Abstract. We analyze te classical discontinuous Galerkin metod for a general parabolic
More informationKey words. Sixth order problem, higher order partial differential equations, biharmonic problem, mixed finite elements, error estimates.
A MIXED FINITE ELEMENT METHOD FOR A SIXTH ORDER ELLIPTIC PROBLEM JÉRÔME DRONIOU, MUHAMMAD ILYAS, BISHNU P. LAMICHHANE, AND GLEN E. WHEELER Abstract. We consider a saddle point formulation for a sixt order
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 informationarxiv: v1 [physics.comp-ph] 12 Jan 2014
arxiv:4.259v [pysics.comp-p] 2 Jan 24 Caracteristic local discontinuous Galerkin metods for solving time-dependent convection-dominated Navier-Stokes equations Suqin Wang Scool of Matematics and Statistics,
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 informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationMixed Finite Element Methods for Incompressible Flow: Stationary Stokes Equations
Mixed Finite Element Metods for Incompressible Flow: Stationary Stoes Equations Ziqiang Cai, Carles Tong, 2 Panayot S. Vassilevsi, 2 Cunbo Wang Department of Matematics, Purdue University, West Lafayette,
More informationCELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7, Number 1, Pages 1 9 c 010 Institute for Scientific Computing and Information CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More informationERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS
ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis
More informationA trace finite element method for a class of coupled bulk-interface transport problems
Numerical Analysis and Scientific Computing Preprint Seria A trace finite element metod for a class of coupled bulk-interface transport problems S. Gross M.A. Olsanskii A. Reusken Preprint #28 Department
More informationDownloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
More informationThird order Approximation on Icosahedral Great Circle Grids on the Sphere. J. Steppeler, P. Ripodas DWD Langen 2006
Tird order Approximation on Icosaedral Great Circle Grids on te Spere J. Steppeler, P. Ripodas DWD Langen 2006 Deasirable features of discretisation metods on te spere Great circle divisions of te spere:
More informationGuaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: a unified framework
Guaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: a unified framework Eric Cancès, Geneviève Dusson, Yvon Maday, Benjamin Stamm, Martin Voralík To cite tis version: Eric
More informationInstitut für Numerische und Angewandte Mathematik
Institut für Numerisce und Angewandte Matematik Hybrid discontinuous Galerkin metods wit relaxed H(div-conformity for incompressible flows. Part I Lederer, P. L., Lerenfeld, C., Scöberl, J. Nr. 4 Preprint-Serie
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More information