c 2007 Society for Industrial and Applied Mathematics
|
|
- Brittany Austin
- 5 years ago
- Views:
Transcription
1 SIAM J. NUMER. ANAL. Vol. 45, No. 3, pp c 2007 Socity for Industrial and Applid Mathmatics NEW FINITE ELEMENT METHODS IN COMPUTATIONAL FLUID DYNAMICS BY H (DIV) ELEMENTS JUNPING WANG AND XIU YE Abstract. In this papr, th authors prsnt two formulations for th Stoks problm which mak us of th xisting H(div) lmnts of th Raviart Thomas typ originally dvlopd for th scond-ordr lliptic problms. In addition, two nw H(div) lmnts ar constructd and analyzd particularly for th nw formulations. Optimal-ordr rror stimats ar stablishd for th corrsponding finit lmnt solutions in various Sobolv norms. Th finit lmnt solutions fatur a full satisfaction of th continuity quation whn xisting Raviart Thomas-typ lmnts ar mployd in th numrical schm. y words. finit lmnt mthods, Stoks problm AMS subjct classifications. Primary, 65N15, 65N30, 76D07; Scondary, 35B45, 35J50 DOI / Introduction. This papr is concrnd with numrical solutions of incomprssibl fluid flow problms by finit lmnt mthods. Our objctiv is to introduc a finit lmnt schm with attntion paid to th discrtization of th mass continuity quation. For illustrativ purposs, w show how th mthod works for th Stoks problm, which sks a pair of unknown functions (u; p) satisfying (1.1) (1.2) (1.3) ν u + p = f in Ω, u = 0 in Ω, u = 0 on Ω, whr ν dnots th fluid viscosity;,, and dnot th Laplacian, gradint, and divrgnc oprators, rspctivly; Ω R d is th rgion occupid by th fluid; f = f(x) (L 2 (Ω)) d is th unit xtrnal volumtric forc acting on th fluid at x Ω. Th commonly usd finit lmnt mthods for th Stoks problm (1.1) (1.3) ar basd on a variational quation which is obtaind by tsting th momntum quation (1.1) by functions in (H0 1 (Ω)) d and th continuity quation (1.2) by functions in L 2 (Ω) (s sction 2 for thir dfinition). Th corrsponding finit lmnt mthod rquirs a pair of finit lmnt spacs which ar conforming in (H0 1 (Ω)) d L 2 (Ω) and satisfy th inf-sup condition of Babu ska [2] and Brzzi [3]. Ths constraints rsult in finit lmnt approximations, dnotd by (u h ; p h ), which hardly satisfy th continuity quation (1.4) u h (x) = 0 x Ω. Radrs ar rfrrd to [8, 19, 21] for mor dtails rgarding th approximation mthods and thir proprtis. Rcivd by th ditors January 5, 2006; accptd for publication (in rvisd form) Fbruary 1, 2007; publishd lctronically May 22, Division of Mathmatical Scincs, National Scinc Foundation, Arlington, VA (jwang@ nsf.gov). This author s rsarch was supportd in part by th NSF IR/D program. This rsarch was initiatd whil th author was visiting NASA Cntr for Computational Scincs (NCCS) at GSFC through th NASA/ASEE Faculty Fllowship Program. Dpartmnt of Mathmatics, Univrsity of Arkansas at Littl Rock, Littl Rock, AR (xxy@ualr.du). 1269
2 1270 JUNPING WANG AND XIU YE Th rcnt dvlopmnt in discontinuous Galrkin mthods [1, 4, 5, 6, 10, 11, 13] provids nw mans of solving th Stoks quations numrically. Howvr, th corrsponding finit lmnt solutions ar usually totally discontinuous and fail to satisfy th continuity quation (1.4) in th classical sns [12, 22, 24, 27]. Th continuity quation (1.4) rquirs th numrical solution u h to b a mmbr of th Sobolv spac H(div; Ω). Thrfor, th discontinuous Galrkin mthods [12, 22, 24, 27] appar to b noncomptitiv whn (1.4) nds to b satisfid. On th othr hand, th (H0 1 ) d L 2 conforming finit lmnt mthods rquir th total continuity of u h, which is too much to satisfy for (1.4). Thrfor, it sms that th H(div) lmnts of Raviart Thomas typ [25, 7, 8, 17] might b good candidats for producing nw numrical schms that satisfy (1.4). Th goal of this papr is to prsnt a mthod that dmonstrats th us of H(div) lmnts in solving th Stoks problm. Our main contribution is on th dvlopmnt of a nw formulation for th Stoks problm which maks us of th xisting H(div) lmnts in numrical schms. Optimal-ordr rror stimats ar drivd for th rsulting H(div) finit lmnt approximations. In addition, two nw familis of H(div) lmnts ar proposd and analyzd in this articl. This papr is organizd as follows. In sction 2, w introduc som prliminaris and notations for Sobolv spacs. A nw variational formula is prsntd in sction 3 for th Stoks problm. In sction 4, w prsnt a H(div) finit lmnt mthod by using th variational formula dvlopd in sction 3. In sction 5, w stablish som optimal-ordr rror stimats for th nw finit lmnt approximations in H 1 and L 2 norms. Finally, in sction 6, w rviw som rprsntativs of H(div) lmnts, followd with a dtaild dscription of two nw H(div) lmnts. 2. Prliminaris and notations. Lt D b any domain in R d, d = 2, 3. For simplicity, th mthod will b prsntd for two-dimnsional problms only. An xtnsion to highr-dimnsional problms can b mad formally for gnral polyhdral domains. W us standard dfinitions for th Sobolv spacs H s (D) and thir associatd innr products (, ) s,d, norms s,d, and sminorms s,d for s 0. For xampl, for any intgr s 0, th sminorm s,d is givn by with th usual notation v s,d = Th Sobolv norm m,d is givn by α =s D α v 2 dd α = (α 1, α 2 ), α = α 1 + α 2, α = α 1 x 1 α 2 x 2. m v m,d = j=0 v 2 j,d Th spac H 0 (D) coincids with L 2 (D), for which th norm and th innr product ar dnotd by D and (, ) D, rspctivly. Whn D = Ω, w shall drop th subscript D in th norm and innr product notation. W also us L 2 0(Ω) to dnot th subspac of L 2 (Ω) consisting of functions with man valu zro ,
3 H(DIV) FINITE ELEMENT METHOD 1271 Th spac H(div; Ω) is dfind as th st of vctor-valud functions on Ω which, togthr with thir divrgnc, ar squar intgrabl; i.., H(div; Ω) = { v : v (L 2 (Ω)) 2, v L 2 (Ω) }. Th norm in H(div; Ω) is dfind by v H(div;Ω) = ( v 2 + v 2) 1 2. Lt Ω b a triangl or quadrilatral. For any smooth vctor-valud functions w and v, it follows from th divrgnc thorm that w (2.1) ( w) vd = ( w, v) v ds, n whr ds rprsnts th boundary lmnt, n is th outward normal dirction on, and ( w, v) = 2 i,j=1 w i x j v i x j d. Lt τ b th tangntial dirction to so that n and τ form a right-hand coordinat systm. It follows from th rprsntation that (2.2) v = (v n )n + (v τ )τ w v = (w n ) (v n ) + (w τ ) (v τ ). n n n 3. A variational formula. For simplicity, w lt ν = 1 for th fluid viscosity in th Stoks quation (1.1). Furthrmor, w assum that Ω is a plan polygonal domain without cracks. Lt T h b a finit lmnt partition of th domain Ω with msh siz h. Assum that th partition T h is shap rgular so that th routin invrs inquality in finit lmnts holds tru (s [9]). Dfin th finit lmnt spacs V h and W h for th vlocity and prssur variabls, rspctivly, by V h = {v H(div; Ω) : v V r () T h ; v n Ω = 0} T h W h = {q L 2 0(Ω) : q W m () T h }, whr n is th outward normal dirction on th boundary of Ω, V r () is a spac of vctor-valud polynomials on th lmnt with indx r 1, and W m () is a st of polynomials on th lmnt with indx m 0. Exampls of V r () and W m () will b givn in sction 6. To driv a wak formulation, w multiply th quation (1.1) by any v V h and us (2.1) to obtain ( ) u (3.1) ( u, v) v ds (p, v) = (f, v), n
4 1272 JUNPING WANG AND XIU YE whr w hav also usd th intgration by parts to dduc p vdω = (p, v). T h Ω Th fact that v V h implis that v n is continuous across ach intrior boundary. Thus, it follows from (2.2) that u v ds = (u τ ) (3.2) v τ ds. n n Introduc th following notation: T h ( h u, h v) = By substituting (3.2) into (3.1) w obtain ( h u, h v) (p, v) (3.3) T h T h ( u, v). (u τ ) n v τ ds = (f, v), which is th basis of our first quation in th nw variational form. Our scond quation can b drivd from tsting (1.2) against any q W h, yilding (3.4) ( u, q) = 0. W now rformulat th boundary intgrals in (3.3). Lt b an intrior dg shard by two lmnts 1 and 2, and lt n 1 and n 2 b unit normal vctors on pointing xtrior to 1 and 2, rspctivly. Dnot by τ 1 and τ 2 th two tangntial dirctions which mak th right-hand coordinat systms with n 1 and n 2, rspctivly. W dfin th avrag { } and jump [[ ]] on for vctor-valud functions w as follows: {ε(w) } = 1 2 (n 1 (w τ 1 ) 1 + n 2 (w τ 2 ) 2 ), [[w ] = w 1 τ 1 + w 2 τ 2. For boundary dg = 1 Ω, th abov two oprations must b modifid by {ε(w) } = n 1 (w τ 1 ) 1, [[w]] = w 1 τ 1. Lt E h dnot th union of th boundaris of all lmnts in T h. For sufficintly smooth u (.g., u H 3 2 +ɛ (Ω) for som ɛ > 0), it is not hard to s that (u τ ) v τ ds = {ε(u) }[[v]]ds. n E h T h Substituting th abov into (3.3) w obtain ( h u, h v) ( v, p) (3.5) E h {ε(u) }[[v]]ds = (f, v). Lt V (h) = V h + (H s (Ω) H0 1 (Ω)) 2, with s > 3 2. Dnot by a o (u, v) = ( h u, h v) {ε(u) }[v]]ds E h
5 H(DIV) FINITE ELEMENT METHOD 1273 and b(v, q) = ( v, q) two bilinar forms on V (h) V (h) and V (h) L 2 0(Ω). With th conditions spcifid in this papr, it can b provd that th xact solution (u; p) of th Stoks problm in 2D blongs to V (h) for som s > 3 2. Radrs ar rfrrd to [20, 15, 14, 23] for dtails. As a rsult, it follows from (3.5) and (3.4) that th xact solution of th 2D Stoks problm satisfis th following variational quations: (3.6) (3.7) a o (u, v) b(v, p) = (f, v) v V h, b(u, q) = 0 q W h. Howvr, it is not clar if th sam statmnt can b mad for th Stoks problm in thr-dimnsional spac without assuming a smooth boundary Ω or a convx polyhdral domain Ω [16, 18]. 4. Finit lmnt schms. Our goal of this sction is to propos two finit lmnt schms basd on two modifications of th wak formulation (3.6) (3.7) for th Stoks problm (1.1) (1.3). To this nd, lt us introduc a symmtric bilinar form on V (h) V (h) as follows: a s (w, v) = a o (w, v) + ( αh 1 [[w]][[v]] {ε(v) }[[w]] ) ds, E h whr α > 0 is a paramtr to b dtrmind latr, and h is th lngth of th dg. For th xact solution (u; p) of th Stoks problm, w clarly hav a s (u, v) = a o (u, v) v V h. Thrfor, it follows from (3.6) and (3.7) that (4.1) (4.2) a s (u, v) b(v, p) = (f, v) v V h, b(u, q) = 0 q W h. Th corrsponding finit lmnt schm for (1.1) (1.3) sks (u h ; p h ) V h W h such that (4.3) (4.4) a s (u h, v) b(v, p h ) = (f, v) v V h, b(u h, q) = 0 q W h. To invstigat th proprtis of th abov numrical schm, w introduc two norms 1 and for th st V (h) as follows: (4.5) (4.6) v 2 1 = v 2 1,h + E h h 1 [[v]] 2, v 2 = v E h h {ε(v) } 2, whr v 2 1,h = T h v 2 1, and v 2 = v vds.
6 1274 JUNPING WANG AND XIU YE Lt b an lmnt with as an dg. It is wll known that thr xists a constant C such that for any function g H 1 () (4.7) In particular, for any v V h, w hav g 2 C ( h 1 g 2 + h g 2 ). (4.8) h {ε(v) } 2 C ( v 2 + h 2 2 v 2 ). Th standard invrs inquality can b mployd to th last trm of th abov inquality, yilding (4.9) h {ε(v) } 2 C v 2 for som constant C indpndnt of th msh siz h. Consquntly, thr is a constant C such that (4.10) v C 0 v 1 v V h. Th following rsult is concrnd with th llipticity of th bilinar form a s (, ) in V h V h. Lmma 4.1. Thr xists a constant α 0 indpndnt of h such that for any v V h w hav (4.11) a s (v, v) α 0 v 2, providd that α is sufficintly larg. Proof. It follows from th Cauchy Schwarz inquality that thr is a constant C such that E h ( {ε(w) }[[v]]ds C h {ε(w) } 2 E h ( C w 1,h h 1 [[v]] 2 E h ) 1 2 ( ) 1 2 h 1 E h [[v ] 2 ) w 2 1,h + C E h h 1 [[v]] 2, whr w hav usd th stimat (4.9) in th scond lin. Using th abov inquality and (4.10), w obtain a s (v, v) = ( h v, h v) + α h 1 [v]] 2 ds 2 {ε(v) }[[v ]ds E h E h v 2 1,h + α E h h 1 [[v]] v 2 1,h C h 1 [[v ] 2 E h = 1 2 v 2 1,h + (α C) E h h 1 [[v ] 2 α 1 v 2 1 α 0 v 2, with α 1 = min( 1 2, α C) and α 0 = α 1 /C 0. For xampl, on may hav α 0 = 1/(2C 0 ) if th paramtr α is chosn so that α C
7 H(DIV) FINITE ELEMENT METHOD 1275 In th rst of th papr, w assum that th paramtr α is chosn so that (4.11) holds tru for th symmtric bilinar from a s (, ). Th proof of Lmma 4.1 indicats that th valu of α dpnds upon th constant in th invrs inquality for finit lmnt functions. Thrfor, th valu of α for which a s (, ) is corciv is mshdpndnt. Existing rsults for saddl-point problms indicat that it is thortically and computationally important to hav th corcivity (4.11). Thrfor, th msh dpndnc of th paramtr α maks th finit lmnt schm (4.3) (4.4) conditionally intrsting in practical computation. To ovrcom th difficulty on th paramtr slction, w introduc a scond finit lmnt schm which is paramtr-insnsitiv. To this nd, w dfin a nonsymmtric bilinar form on V (h) V (h) as follows: a ns (w, v) = a o (w, v) + E h ( αh 1 [[w ][[v ] + {ε(v) }[[w]] ) ds. Similar to th bilinar form a s (, ), for th xact solution (u; p) of th Stoks problm w hav a ns (u, v) = a o (u, v) v V h. Consquntly, th solution of th Stoks problm satisfis th following variational quations: (4.12) (4.13) a ns (u, v) b(v, p) = (f, v) v V h, b(u, q) = 0 q W h. Our scond finit lmnt schm for (1.1) (1.3) sks (u h ; p h ) V h W h such that (4.14) (4.15) a ns (u h, v) b(v, p h ) = (f, v) v V h, b(u h, q) = 0 q W h. To s th corcivity of th bilinar form a ns (, ), w us its dfinition and (4.10) to obtain a ns (v, v) = ( h v, h v) + α h 1 [[v]] 2 ds E h min(1, α) v 2 1 min(1, α)c 1 0 v 2, whr v V h. Thus, th corcivity (4.11) holds tru for th bilinar form a ns (, ) with any valu of α > 0. Th following is a rsult on th bounddnss of th bilinar forms a s (, ) and a ns (, ). Lmma 4.2. Thr xists a constant C indpndnt of h such that (4.16) a i (w, v) C w v w, v V (h), whr i = s, ns.
8 1276 JUNPING WANG AND XIU YE Proof. For simplicity, w shall prsnt th analysis for a s (, ) only. By th dfinition of a s (w, v) and th Schwarz inquality, thr xists a constant C such that ( ) 1 ( ) a s (w, v) C w 1,h v 1,h + h {ε(w) } 2 h 1 [[v]] 2 E h E h ( ) 1 ( ) h {ε(v) } 2 h 1 [w]] 2 E h E h ( ) 1 ( ) α h 1 [[w[] 2 h 1 [v[] 2 E h E h C w v, which provs th dsird bounddnss. 5. Error stimats. Th first goal of this sction is to driv an optimal-ordr rror stimat for th prssur in L 2 (Ω) and th vlocity in th norm givn by (4.6). Th scond goal is to driv an optimal-ordr rror stimat for th vlocity approximation in th L 2 -norm for th symmtric schm (4.3) (4.4). Assumption 1. Thr xists an oprator Π h : (H 1 0 (Ω)) 2 V h such that (5.1) b(v Π h v, q) = 0 q W h. In addition, th oprator Π h is assumd to satisfy th following: (5.2) v Π h v s, Ch t s v t, T h, s = 0, 1, whr th constant C dpnds only on th shap of and 1 t r + 1. From (5.2) and th inquality (4.7) it is not hard to s that v Π h v 1 C v 1 v (H 1 0 (Ω)) 2. Thus, it follows from v 1 = v 1 v 1 and th triangl inquality that (5.3) Π h v 1 C v 1. For our finit lmnt formulations, th inf-sup condition givn in Brzzi s framwork would rad as follows: Thr xists a positiv constant β, indpndnt of h, such that b(v, q) (5.4) sup β q q W h. v V h v To vrify (5.4), w first us th oprator Π h to obtain (5.5) b(v, q) b(π h v, q) b(v, q) sup sup = sup v V h v v (H0 1 Π (Ω))2 h v v (H0 1 Π (Ω))2 h v. Obsrv that by using (5.3), and (4.10), w hav for all v (H 1 0 (Ω)) 2 (5.6) Π h v C Π h v 1 C v 1. Thus, substituting (5.6) into th inquality (5.5) givs b(v, q) sup v V h v C 1 b(v, q) sup β q, v (H0 1 v (Ω))2 1 whr w hav usd th inf-sup condition for th continuous cas [19, 8].
9 H(DIV) FINITE ELEMENT METHOD Error stimats in H 1 L 2. Th rror analysis rquirs th L 2 projction from L 2 0(Ω) to th finit lmnt spac W h, which is dnotd by Q h. In addition, th following rror quations turn out to b usful: (5.7) (5.8) a s (u u h, v) b(v, p p h ) = 0 v V h, b(u u h, q) = 0 q W h. Ths rror quations can b obtaind by subtracting (4.3) (4.4) from (4.1) (4.2). Similar rror quations hold tru for th nonsymmtric schm (4.14) (4.15) with a s (, ) bing rplacd by a ns (, ). Now w ar in a position to prsnt an rror stimat for th nw finit lmnt approximations. Thorm 5.1. Lt (u; p) b th solution of (1.1) (1.3) and (u h ; p h ) V h W h b obtaind from ithr (4.3) (4.4) or (4.14) (4.15). Assum that Assumption 1 holds tru. Thn, thr xists a constant C indpndnt of h such that (5.9) u u h + p p h C ( u Π h u + p Q h p ). Proof. Lt ξ h = u h Π h u, η h = p h Q h p b th rror btwn th finit lmnt solution (u h ; p h ) and th projction (Π h u; Q h p) of th xact solution. Dnot by ξ = u Π h u, η = p Q h p th rror btwn th xact solution (u; p) and its projction. It follows from th rror quations (5.7) and (5.8) that (5.10) (5.11) a(ξ h, v) b(v, η h ) = a(ξ, v) b(v, η), b(ξ h, q) = b(ξ, q) = 0 for any v V h and q W h. Hr and in what follows of this sction, a(, ) dnots ithr a s (, ) or a ns (, ). By ltting v = ξ h in (5.10) and q = η h in (5.11), th sum of (5.10) and (5.11) givs a(ξ h, ξ h ) = a(ξ, ξ h ) b(ξ h, η). Thus, it follows from th corcivity (4.11) and th bounddnss (4.16) that which implis th following: Th abov stimat can b rwrittn as α 0 ξ h 2 C( ξ ξ h + η ξ h ), ξ h C ( ξ + η ). (5.12) u h Π h u C ( u Π h u + p Q h p ). Now using th triangl inquality and th rror stimat (5.12) w gt (5.13) u u h C ( u Π h u + p Q h p ), which complts th stimat for th vlocity approximation.
10 1278 JUNPING WANG AND XIU YE It rmains to stimat th prssur approximation p h. To this nd, w us th discrt inf-sup condition (5.4) to obtain p h Q h p 1 β sup b(v, p h Q h p) v V h v which, togthr with (5.13), implis that = 1 β sup b(v, p h p) + b(v, p Q h p) v V h v = 1 β sup a(u u h, v) + b(v, p Q h p) v V h v 1 C sup v V h v v ( u u h + p Q h p ) C( u u h + p Q h p ), p h Q h p C ( u Π h u + p Q h p ). Th rror stimat for th prssur approximation is thn compltd by combining th abov inquality with th standard triangl inquality An L 2 -rror stimat for th vlocity approximation. Considr only th finit lmnt approximat solutions arising from th symmtric finit lmnt schm. To driv an L 2 -rror stimat for th vlocity approximation, w sk (w; λ) (H0 1 (Ω)) 2 L 2 0(Ω) satisfying w + λ = u u h in Ω, w = 0 in Ω, w = 0 on Ω. Not that for any (v; q) V (h) L 2 0(Ω) th solution (w; λ) satisfis (5.14) (5.15) a s (w, v) b(v, λ) = (u u h, v), b(w, q) = 0. Assum that th Stoks problm has th H 2 (Ω) H 1 (Ω)-rgularity proprty in th sns that th solution (w; λ) (H 2 (Ω)) 2 H 1 (Ω) and th following a priori stimat holds tru: (5.16) w 2 + λ 1 C u u h. In addition, w assum that th finit lmnt spac V h and th projction oprator Π h hav th following proprty: (5.17) w Π h w Ch w 2. With ths assumptions, it is not hard to s that thr xists a constant C indpndnt of h such that (5.18) w Π h w + λ Q h λ Ch u u h. It must b pointd out that th H 2 H 1 -rgularity proprty assumption statd as abov rquirs that th polygonal domain Ω b convx. For nonconvx but smooth
11 H(DIV) FINITE ELEMENT METHOD 1279 domains, th rgularity (5.16) can b provd to b valid. Howvr, isoparamtric lmnts would nd to b mployd in th finit lmnt schm in ordr to maintain optimal-ordr rror stimats in ithr H 1 - or L 2 -norms. Thorm 5.2. Lt (u h ; p h ) V h W h and (u; p) b th solutions of (4.3) (4.4) and (1.1) (1.3), rspctivly. Assum that Assumption 1 and th stimat (5.17) hold tru and that th Stoks problm (1.1) (1.3) has th H 2 (Ω) H 1 (Ω)-rgularity proprty. Thn thr xists a constant C indpndnt of h such that (5.19) u u h Ch( u Π h u + p Q h p ). Proof. By ltting v = u u h in (5.14) w arriv at (5.20) a s (u u h, w) b(u u h, λ) = u u h 2. Notic that (5.21) b(u u h, λ) = b(u u h, λ Q h λ) and (5.22) a s (u u h, w) = a s (u u h, w Π h w) + a s (u u h, Π h w) = a s (u u h, w Π h w) + b(π h w, p p h ) = a s (u u h, w Π h w) + b(π h w w, p p h ). Substituting (5.21) and (5.22) into (5.20) w obtain u u h 2 = a s (u u h, w Π h w) + b(π h w w, p p h ) b(u u h, λ Q h λ). Thus, u u h 2 C ( u u h + p p h ) ( w Π h w + λ Q h λ ). Substituting (5.18) into th abov stimat w obtain which implis that u u h 2 Ch ( u u h + p p h ) u u h, u u h Ch ( u u h + p p h ). Th abov inquality and th rror stimat (5.9) imply This complts th proof of th thorm. u u h Ch( u Π h u + p Q h p ). 6. Exampls of H(div) lmnts. Lt us rcall that th rror stimats stablishd in sction 5 ar basd on th following thr proprtis: B1. V h H(div; Ω), B2. Assumption 1 and th stimat (5.17) as dscribd in sction 5, and B3. th H 2 H 1 -rgularity proprty assumption for th Stoks problm.
12 1280 JUNPING WANG AND XIU YE Th last proprty (B3) is rquird only for th L 2 -rror stimat for th vlocity approximation. This mans that any finit lmnt pair V h W h satisfying proprtis B1 B2 is applicabl for th formulations prsntd arlir in this manuscript. Dnot by P k () th spac of polynomials of dgr k and P k1,k 2 () = p(x 1, x 2 ) : p(x 1, x 2 ) = a ij x i 1x j 2. 0 i k 1,0 j k 2 P k1,k 2,k 3 () is dfind similarly in thr-dimnsional spacs. Dfin Q k () as follows: { Pk,k () for d = 2, Q k () = P k,k,k () for d = 3. Obsrv that th finit lmnt pair V h W h is constructd from local lmnts V r () and W m () as dscribd in sction 3. Thrfor, it suffics to spcify th local pair V r () W m () for ach xampl to b prsntd Existing lmnts. All of th xisting H(div) lmnts dsignd for th scond-ordr lliptic problms (.g., s [25, 8, 7, 17, 19]) satisfy proprtis B1 B2, xcpt th stimat (5.17) for th lowst-ordr Raviart Thomas lmnt on triangls and quadrilatrals. Thrfor, thr ar plnty of finit lmnt spacs applicabl to th nw formulation of th Stoks problm. For illustrativ purposs, w mntion thr xampls. Radrs ar rfrrd to th book by Brzzi and Fortin [8] for mor xampls of th H(div) lmnt Raviart Thomas lmnts on triangls or ttrahdra: RT k (). Lt k 1 b any intgr. For any triangular or ttrahdral lmnt, th local lmnt V r () W m () is dfind by V k () = (P k ()) d xp k (), W k () = P k (), whr d = 2 if is a triangl and d = 3 if is a ttrahdron. Th projction oprator Π h satisfying all of th rquird proprtis is givn locally on ach lmnt. For xampl, th rstriction of Π h on th lmnt, dnotd by Π, is dfind as follows: (v Π v) nqds = 0 q P k ( ), (v Π v) qd = 0 q (P k 1 ()) d, k BDM lmnts on triangls or ttrahdra: BDM k () [8]. Lt k 1 b any intgr. For any triangular or ttrahdral lmnt, th local lmnt V r () W m () is dfind by stting r = m + 1 = k and V k () = (P k ()) d, W k 1 () = P k 1 (). On a triangular lmnt, th local projction oprator Π : (H 1 ()) 2 V k () is dfind by (v Π v) nqds = 0 q P k ( ), (v Π v) qd = 0 q P k 1 (), (v Π v) curl(b q)d = 0 q P k 2 (), k 2,
13 H(DIV) FINITE ELEMENT METHOD 1281 whr b is th bubbl function dfind on. On a ttrahdral lmnt, th corrsponding local projction Π is givn by (v Π v) nqds = 0 q P k ( ), (v Π v) qd = 0 q P k 1 (), (v Π v) qd = 0 q Φ k (), whr Φ() = {φ (P k ()) 3 : φ = 0, φ n = 0 on } BDM lmnts on quadrilatrals: BDM [k] (). It is sufficint to dscrib th lmnt on th unit squar. Lt k 1 b any intgr. Th local lmnt V r () W m () is dfind by V k () = (P k ()) 2 curl(x k+1 1 x 2 ) curl(x 1 x k+1 2 ), W k 1 () = P k 1 (). On th unit squar lmnt, th local projction oprator Π : (H 1 ()) 2 V k () is dfind by (v Π v) nqds = 0 q P k ( ), (v Π v) wd = 0 w (P k 2 ()) 2, k Error stimats for th xisting lmnts. Rcall that th vlocity V h and th prssur spac W h ar dfind, rspctivly, by (6.1) V h = {v H(div; Ω) : v V r () T h ; v n Ω = 0} and (6.2) W h = {q L 2 0(Ω) : q W m () T h }. For th xisting H(div) lmnts listd abov, w hav V r () = RT k (), BDM k (), or BDM [k] () and W m () = P k (), P k 1 (), or P k 1 (), rspctivly. Th projction oprator Π h is givn by (6.3) (Π h v) = Π (v ). Th dfinition of Π h implis that (6.4) b(v Π h v, q) = 0 q W h. Furthrmor, it has bn provd in [8] that (5.2) and (5.17) hold tru for Π h dfind in (6.3). Thrfor, proprtis B1 B2 ar wll justifid. Lt Q h b th L 2 projction from L 2 0(Ω) to W h. It is not hard to s that W h has th following local approximation proprtis: For BDM k () and BDM [k] () (6.5) p Q h p s, Ch k s p k, T h, s = 0, 1,
14 1282 JUNPING WANG AND XIU YE and for RT k () (6.6) p Q h p s, Ch k+1 s p k+1, T h, s = 0, 1. Th constant C in (6.5) (6.6) dpnds only on k and th shap of. Th following rsult follows from (5.17), (6.5) (6.6), and Thorms 5.1 and 5.2. Proposition 6.1. Lt (u; p) b th solution of (1.1) (1.3) and (u h ; p h ) V h W h b obtaind from ithr (4.3) (4.4) or (4.14) (4.15). Assum that (u; p) (H t+1 (Ω)) 2 H t (Ω) for som 1 t k. Thn thr xists a constant C indpndnt of h such that for BDM k (), BDM [k] (), and RT k () (6.7) u u h + p p h Ch t ( u t+1 + p t ). Furthrmor, if th H 2 H 1 -rgularity proprty holds tru for th Stoks problm, thn thr is a constant C such that th finit lmnt approximation (u h ; p h ) from th symmtric formulation has th following rror stimat: u u h Ch t+1 ( u t+1 + p t ). W commnt that th abov rror stimats hold tru for all of th H(div) lmnts listd in [8] Nw lmnts. Stability and accuracy ar two main factors in th construction of nw finit lmnts. For th variational schms prsntd in this papr, th stability part is ralizd by a combination of th inf-sup condition and th corcivity for th corrsponding bilinar forms. Th accuracy part is charactrizd by Assumption 1 and a balancd prssur spac W h. For xampl, Proposition 6.1 indicats that th Raviart Thomas lmnt can b usd to approximat th solution of th Stoks quations, which is convrgnt as th msh siz dcrass (not that w do not know any convrgnc whn RT 0 () is usd). Howvr, th Raviart Thomas lmnts do not appar to b wll balancd, bcaus both th vlocity and th prssur unknowns ar approximatd by polynomials of ordr k. In contrast, th BDM lmnts offr a bttr/optimal combination for th solution of th Stoks problm. But th BDM [k] () lmnt on rctangls is constructd in an awkward way by involving th curl of som polynomials. W fl that bttr constructd lmnts should b xplord on rctangls and cubs for solving th Stoks problm. For this purpos, w would lik to propos som altrnativs on rctangls and cubs which ar suitabl for approximating th solution of th Stoks problm. Ths lmnts can b usd on quadrilatrals through local transformations as dscribd in [26] A nw lmnt on rctangls: NE1 k (). W illustrat th construction of th nw NE1 k () lmnt on th unit squar = [0, 1] [0, 1]. Lt k 1 b any intgr. W dfin local lmnts V r () W m () by V k () = (Q k ()) 2, W k 1 () = Q k 1 (). For th first componnt of v = (v 1, v 2 ), w dfin an oprator Π,1 : H 1 () Q k () as follows: (6.8) (v 1 Π,1 v 1 )φds = 0 φ P k (), = wst, ast, (6.9) (v 1 Π,1 v 1 )ψd = 0 ψ P k 2,k (),
15 H(DIV) FINITE ELEMENT METHOD 1283 whr = wst mans that is th wst dg (i.., = {(0, x 2 ) : x 2 [0, 1]}) of th unit squar; th ast dg is dfind accordingly. Th systm (6.8) involvs xactly 2(k + 1) linar quations and (6.9) involvs (k 1)(k + 1) linar quations. Th total numbr of quations is givn by 2(k + 1) + (k 1)(k + 1) = (k + 1) 2, which is th sam as th total numbr of dgrs of frdom for a polynomial in Q k (). Th following proposition shows that th linar systms (6.8) and (6.9) uniquly dtrmin th projction Π,1 v 1. Proposition 6.2. Lt v Q k () b such that (6.10) vφds = 0 φ P k (), = wst, ast, (6.11) vψd = 0 ψ P k 2,k (). Thn w must hav v 0. Proof. Th condition (6.10) implis that v = 0 at th ast and wst dgs of th unit squar. Thus, thr is a polynomial g = g(x 1, x 2 ) P k 2,k () such that v = x 1 (1 x 1 )g. Substitut v = x 1 (1 x 1 ) into (6.11), and thn lt ψ = g. It follows that g 0. This shows that v 0. Th projction of th scond componnt of v, dnotd by Π,2 v 2, can b dfind in a similar fashion. Th local projction oprator is thn givn by Π v = (Π,1 v 1, Π,2 v 2 ). It is not hard to show that such a dfind projction satisfis all of th conditions rquird in th prvious sctions. As a rsult, th lmnt NE1 k () can b usd to approximat th Stoks problm A nw lmnt on cubs: NE2 k (). Again, w shall dscrib dtails only on th unit cub = [0, 1] 3. Lt k 1 b an intgr. A straightforward xtnsion of th NE1 k lmnt to thr-dimnsional spac is givn by V k () = (Q k ()) 3, W k 1 () = Q k 1 (). Our goal hr is to show that th abov xtnsion actually works. To this nd, it suffics to construct a projction oprator Π which satisfis th rquird proprtis. Lt v = (v 1, v 2, v 3 ) ( H 1 (Ω) ) 3 b a vctor-valud function. For ach componnt v i, w dfin its projction to Q k () as follows: (6.12) (v i Π,i v i )φds = 0 φ Q k ( i ), i (6.13) (v i Π,i v i )ψd = 0 ψ P k1,k 2,k 3 (), whr i = {(x 1, x 2, x 3 ) : x j [0, 1], j i; x i = 0 or 1} ar th two facs of th cub which ar orthogonal to th x i -axis, and k i = k 2, k j = k for j i. Thr ar 2(k + 1) 2 linar quations from th condition (6.12) and (k 1)(k + 1) 2 linar quations from th condition (6.13). Th total numbr of linar quations is thn givn by 2(k + 1) 2 + (k 1)(k + 1) 2 = (k + 1) 3,
16 1284 JUNPING WANG AND XIU YE which is th sam as th total numbr of dgrs of frdom for a polynomial in th spac V k (). A similar argumnt as in th prvious subsction for NE1 k () can b applid to show that th projction Π,i v i is uniquly dtrmind by (6.12) and (6.13). Furthrmor, th local projction givn by Π v = (Π,1 v 1, Π,2 v 2, Π,3 v 3 ) can b vrifid to satisfy all of th proprtis rquird by th convrgnc thory dvlopd in prvious sctions for th nw finit lmnt mthods Anothr nw lmnt on rctangls: NE3 k (). Again for simplicity, w shall dscrib th nw lmnt on th unit squar = [0, 1] 2. This lmnt will b a simplifid vrsion of NE1 k () but with th sam ordr of accuracy. Lt k 1 b any intgr, and dfin V k () = (P k () {x 1 x k 2}) (P k () {x 2 x k 1}), W k 1 () = P k 1 (). For th first componnt of v = (v 1, v 2 ), w dfin its projction Π,1 v 1 P k () {x 1 x k 2} by using th following quations: (6.14) (v 1 Π,1 v 1 )φds = 0 φ P k (), = wst, ast, (6.15) (v 1 Π,1 v 1 )ψd = 0 ψ P k 2 (). Thr ar 2(k + 1) quations from th condition (6.14) and 1 2 (k 1)k quations from th condition (6.15). Th total numbr of linar quations is givn by 2(k + 1) (k 1)k = 1 (k + 1)(k + 2) + 1, 2 which is th sam as th total numbr of dgrs of frdom for functions in th spac P k () {x 1 x k 2}. Using th sam tchniqu as in th analysis for NE1 k (), it can b provd that Π,1 v 1 is uniquly dtrmind by (6.14) and (6.15). Th projction of th scond componnt of v can b dtrmind in a similar way. Th rsulting local projction Π v = (Π,1 v 1, Π,2 v 2 ) satisfis all of th proprtis rquird in th convrgnc thory Error stimats for th nw lmnts. First, w dfin th vlocity spac V h by (6.16) V h = {v H(div; Ω) : v V r () T h ; v n Ω = 0} and th prssur spac W h by (6.17) W h = {q L 2 0(Ω) : q W m () T h }, whr V r () = NE1 k (), NE2 k (), or NE3 k () and W m () = Q k 1 (), Q k 1 (), or P k 1 (), rspctivly. For any v (H 1 0 (Ω)) d, with d = 2, 3, dfin Π h v V h by (6.18) (Π h v) = Π v T h,
17 H(DIV) FINITE ELEMENT METHOD 1285 whr Π is th corrsponding local projction oprator on ach lmnt. From th construction of Π, it is asy to s that it is indd tru that Π h v V h, and, morovr, on has (6.19) b(v Π h v, q) = 0 q W h and that proprtis B1 B2 ar satisfid for th thr nw lmnts NE1 k (), NE2 k (), and NE3 k (). Similar to Proposition 6.1, w hav th following convrgnc stimats. Proposition 6.3. Lt (u; p) b th solution of (1.1) (1.3) and (u h ; p h ) V h W h b obtaind from ithr (4.3) (4.4) or (4.14) (4.15) by using th nw lmnts dscribd in this subsction. Assum that (u; p) (H t+1 (Ω)) d H t (Ω) for som < t k. Thn thr xists a constant C indpndnt of h such that 1 2 u u h + p p h Ch t ( u t+1 + p t ), and for th symmtric formulation w also hav u u h Ch t+1 ( u t+1 + p t ), providd that th H 2 H 1 -rgularity proprty holds tru for th Stoks problm. W point out that, unlik th xisting H(div) lmnts, th nw H(div) lmnts dscribd in this sction do not yild numrical vlocitis that satisfy th continuity quation (1.4) in th classical sns. Howvr, th numrical approximations arising from th nw lmnts indd consrv mass locally on ach lmnt. REFERENCES [1] D. N. Arnold, F. Brzzi, B. Cockburn, and L. D. Marini, Unifid analysis of discontinuous Galrkin mthods for lliptic problms, SIAM J. Numr. Anal., 39 (2002), pp [2] I. Babu ska, Th finit lmnt mthod with Lagrangian multiplir, Numr. Math., 20 (1973), pp [3] F. Brzzi, On th xistnc, uniqunss, and approximation of saddl point problms arising from Lagrangian multiplirs, RAIRO, Anal. Numér., 2 (1974), pp [4] I. Babu ska and M. Zlámal, Nonconforming lmnts in th finit lmnt mthod with pnalty, SIAM J. Numr. Anal., 10 (1973), pp [5] F. Bassi and S. Rbay, A high-ordr accurat discontinuous finit lmnt mthod for th numrical solution of th comprssibl Navir-Stoks quations, J. Comput. Phys., 131 (1997), pp [6] C.E. Baumann and J.T. Odn, A discontinuous hp finit lmnt mthod for convctiondiffusion problms, Comput. Mthods Appl. Mch. Engrg., 175 (1999), pp [7] F. Brzzi, J. Douglas, and L. Marini, Two familis of mixd finit lmnts for scond ordr lliptic problms, Numr. Math., 47 (1985), pp [8] F. Brzzi and M. Fortin, Mixd and Hybrid Finit Elmnts, Springr-Vrlag, Nw York, [9] P.G. Ciarlt, Th Finit Elmnt Mthod for Elliptic Problms, North-Holland, Nw York, [10] B. Cockburn, S. Hou, and C.-W. Shu, TVB Rung-utta local projction discontinuous Galrkin finit lmnt mthod for consrvation laws IV: Th multidimnsional cas, Math. Comp., 54 (1990), pp [11] B. Cockburn, G. anschat, and D. Schotzau, A locally consrvativ LDG mthod for th incomprssibl Navir-Stoks quations, Math. Comp., 74 (2005), pp [12] B. Cockburn, G. anschat, D. Schötzau, and C. Schwab, Local discontinuous Galrkin mthods for th Stoks systm, SIAM J. Numr. Anal., 40 (2002), pp [13] B. Cockburn and C.-W. Shu, Th local discontinuous Galrkin mthod for tim-dpndnt convction-diffusion systms, SIAM J. Numr. Anal., 35 (1998), pp [14] M. Costabl and M. Daug, Crack singularitis for gnral lliptic systms, Math. Nachr., 235 (2002), pp
18 1286 JUNPING WANG AND XIU YE [15] M. Daug, Elliptic Boundary Valu Problms in Cornr Domains - Smoothnss and Asymptotics of Solutions, Lctur Nots in Math. 1341, Springr-Vrlag, Brlin, [16] M. Daug, Stationary Stoks and Navir-Stoks systms on two- or thr-dimnsional domains with cornrs. Part I. Linarizd Equations, SIAM J. Math. Anal., 20 (1989), pp [17] J. Douglas and J. Wang, A nw family of spacs in mixd finit lmnt mthods for rctangular lmnts, Comput. Appl. Math., 12 (1993), pp [18] V. Girault and J.-L. Lions, Two-grid finit-lmnt schms for th stady Navir-Stoks problm in polyhdra, Port. Math. (N.S.), 58 (2001), pp [19] V. Girault and P.A. Raviart, Finit Elmnt Mthods for th Navir-Stoks Equations: Thory and Algorithms, Springr-Vrlag, Brlin, [20] P. Grisvard, Boundary Valu Problms in Non-Smooth Domains, Pitman, London, [21] M. D. Gunzburgr, Finit Elmnt Mthods for Viscous Incomprssibl Flows, A Guid to Thory, Practic and Algorithms, Acadmic, San Digo, [22] P. Hansbo and M. Larson, Discontinuous Galrkin mthods for narly incomprssibl lasticity by Nitsch s mthod, Comput. Mthods Appl. Mch. Engrg., 191 (2002), pp [23] J. R. won, Cornr singularity for comprssibl Navir-Stoks problm, in Procdings of th Intrnational Congrss of Mathmatics, 2002, Highr Education Prss, China, [24] J. Liu and C. Shu, A high ordr discontinuous Galrkin mthod for 2D incomprssibl flows, J. Comput. Phys., 160 (2000), pp [25] P. Raviart and J. Thomas, A mixd finit lmnt mthod for scond ordr lliptic problms, Mathmatical Aspcts of th Finit Elmnt Mthod, I. Galligani, E. Magns, ds., Lcturs Nots in Math. 606, Springr-Vrlag, Nw York, [26] J. Wang and T. Mathw, Mixd finit lmnt mthod ovr quadrilatrals, in Procdings of th Third Intrnational Confrnc on Advancs in Numrical Mthods and Applications, I. T. Dimov, Bl. Sndov, and P. Vassilvski, ds., World Scintific, Rivr Edg, NJ, 1994, pp [27] X. Y, Discontinuous stabl lmnts for incomprssibl flow, Adv. Comput. Math., 20 (2004) pp
A Weakly Over-Penalized Non-Symmetric Interior Penalty Method
Europan Socity of Computational Mthods in Scincs and Enginring (ESCMSE Journal of Numrical Analysis, Industrial and Applid Mathmatics (JNAIAM vol.?, no.?, 2007, pp.?-? ISSN 1790 8140 A Wakly Ovr-Pnalizd
More informationA LOCALLY CONSERVATIVE LDG METHOD FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
A LOCALLY CONSERVATIVE LDG METHOD FOR THE INCOMPRESSIBLE NAVIER-STOES EQUATIONS BERNARDO COCBURN, GUIDO ANSCHAT, AND DOMINI SCHÖTZAU Abstract. In this papr, a nw local discontinuous Galrkin mthod for th
More informationSymmetric Interior Penalty Galerkin Method for Elliptic Problems
Symmtric Intrior Pnalty Galrkin Mthod for Elliptic Problms Ykatrina Epshtyn and Béatric Rivièr Dpartmnt of Mathmatics, Univrsity of Pittsburgh, 3 Thackray, Pittsburgh, PA 56, U.S.A. Abstract This papr
More informationA ROBUST NUMERICAL METHOD FOR THE STOKES EQUATIONS BASED ON DIVERGENCE-FREE H(DIV) FINITE ELEMENT METHODS
A ROBUST NUMERICAL METHOD FOR THE STOKES EQUATIONS BASED ON DIVERGENCE-FREE H(DIV) FINITE ELEMENT METHODS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. A computational procdur basd on a divrgnc-fr H(div)
More informationAnother view for a posteriori error estimates for variational inequalities of the second kind
Accptd by Applid Numrical Mathmatics in July 2013 Anothr viw for a postriori rror stimats for variational inqualitis of th scond kind Fi Wang 1 and Wimin Han 2 Abstract. In this papr, w giv anothr viw
More informationA NEW FINITE VOLUME METHOD FOR THE STOKES PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volum -, Numbr -, Pags 22 c - Institut for Scintific Computing and Information A NEW FINITE VOLUME METHOD FOR THE STOKES PROBLEMS Abstract. JUNPING
More informationA Family of Discontinuous Galerkin Finite Elements for the Reissner Mindlin plate
A Family of Discontinuous Galrkin Finit Elmnts for th Rissnr Mindlin plat Douglas N. Arnold Franco Brzzi L. Donatlla Marini Sptmbr 28, 2003 Abstract W dvlop a family of locking-fr lmnts for th Rissnr Mindlin
More informationA POSTERIORI ERROR ESTIMATION FOR AN INTERIOR PENALTY TYPE METHOD EMPLOYING H(DIV) ELEMENTS FOR THE STOKES EQUATIONS
A POSTERIORI ERROR ESTIMATION FOR AN INTERIOR PENALTY TYPE METHOD EMPLOYING H(DIV) ELEMENTS FOR THE STOKES EQUATIONS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This papr stablishs a postriori rror
More informationUNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This papr is concrnd with rsidual typ a postriori rror stimators
More informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by D. Klain Vrsion 207.0.05 Corrctions and commnts ar wlcom. Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix A A k I + A + k!
More informationDiscontinuous Galerkin and Mimetic Finite Difference Methods for Coupled Stokes-Darcy Flows on Polygonal and Polyhedral Grids
Discontinuous Galrkin and Mimtic Finit Diffrnc Mthods for Coupld Stoks-Darcy Flows on Polygonal and Polyhdral Grids Konstantin Lipnikov Danail Vassilv Ivan Yotov Abstract W study locally mass consrvativ
More informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by Dan Klain Vrsion 28928 Corrctions and commnts ar wlcom Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix () A A k I + A + k!
More informationA Propagating Wave Packet Group Velocity Dispersion
Lctur 8 Phys 375 A Propagating Wav Packt Group Vlocity Disprsion Ovrviw and Motivation: In th last lctur w lookd at a localizd solution t) to th 1D fr-particl Schrödingr quation (SE) that corrsponds to
More informationAn interior penalty method for a two dimensional curl-curl and grad-div problem
ANZIAM J. 50 (CTAC2008) pp.c947 C975, 2009 C947 An intrior pnalty mthod for a two dimnsional curl-curl and grad-div problm S. C. Brnnr 1 J. Cui 2 L.-Y. Sung 3 (Rcivd 30 Octobr 2008; rvisd 30 April 2009)
More informationAnalysis of a Discontinuous Finite Element Method for the Coupled Stokes and Darcy Problems
Journal of Scintific Computing, Volums and 3, Jun 005 ( 005) DOI: 0.007/s095-004-447-3 Analysis of a Discontinuous Finit Elmnt Mthod for th Coupld Stoks and Darcy Problms Béatric Rivièr Rcivd July 8, 003;
More information3 Finite Element Parametric Geometry
3 Finit Elmnt Paramtric Gomtry 3. Introduction Th intgral of a matrix is th matrix containing th intgral of ach and vry on of its original componnts. Practical finit lmnt analysis rquirs intgrating matrics,
More informationHomotopy perturbation technique
Comput. Mthods Appl. Mch. Engrg. 178 (1999) 257±262 www.lsvir.com/locat/cma Homotopy prturbation tchniqu Ji-Huan H 1 Shanghai Univrsity, Shanghai Institut of Applid Mathmatics and Mchanics, Shanghai 272,
More informationIntroduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013
18.782 Introduction to Arithmtic Gomtry Fall 2013 Lctur #20 11/14/2013 20.1 Dgr thorm for morphisms of curvs Lt us rstat th thorm givn at th nd of th last lctur, which w will now prov. Thorm 20.1. Lt φ:
More informationA Family of Discontinuous Galerkin Finite Elements for the Reissner Mindlin Plate
Journal of Scintific Computing, Volums 22 and 23, Jun 2005 ( 2005) DOI: 10.1007/s10915-004-4134-8 A Family of Discontinuous Galrkin Finit Elmnts for th Rissnr Mindlin Plat Douglas N. Arnold, 1 Franco Brzzi,
More informationA LOCALLY DIVERGENCE-FREE INTERIOR PENALTY METHOD FOR TWO-DIMENSIONAL CURL-CURL PROBLEMS
A LOCALLY DIVERGENCE-FREE INTERIOR PENALTY METHOD FOR TWO-DIMENSIONAL CURL-CURL PROBLEMS SUSANNE C. BRENNER, FENGYAN LI, AND LI-YENG SUNG Abstract. An intrior pnalty mthod for crtain two-dimnsional curl-curl
More informationFinite element discretization of Laplace and Poisson equations
Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization
More informationRELIABLE AND EFFICIENT A POSTERIORI ERROR ESTIMATES OF DG METHODS FOR A SIMPLIFIED FRICTIONAL CONTACT PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volum 16, Numbr 1, Pags 49 62 c 2019 Institut for Scintific Computing and Information RELIABLE AND EFFICIENT A POSTERIORI ERROR ESTIMATES OF DG
More informationAPPROXIMATION THEORY, II ACADEMIC PRfSS. INC. New York San Francioco London. J. T. aden
Rprintd trom; APPROXIMATION THORY, II 1976 ACADMIC PRfSS. INC. Nw York San Francioco London \st. I IITBRID FINIT L}ffiNT }ffithods J. T. adn Som nw tchniqus for dtrmining rror stimats for so-calld hybrid
More informationNumerische Mathematik
Numr. Math. 04 6:3 360 DOI 0.007/s00-03-0563-3 Numrisch Mathmatik Discontinuous Galrkin and mimtic finit diffrnc mthods for coupld Stoks Darcy flows on polygonal and polyhdral grids Konstantin Lipnikov
More informationDerangements and Applications
2 3 47 6 23 Journal of Intgr Squncs, Vol. 6 (2003), Articl 03..2 Drangmnts and Applications Mhdi Hassani Dpartmnt of Mathmatics Institut for Advancd Studis in Basic Scincs Zanjan, Iran mhassani@iasbs.ac.ir
More informationLINEAR DELAY DIFFERENTIAL EQUATION WITH A POSITIVE AND A NEGATIVE TERM
Elctronic Journal of Diffrntial Equations, Vol. 2003(2003), No. 92, pp. 1 6. ISSN: 1072-6691. URL: http://jd.math.swt.du or http://jd.math.unt.du ftp jd.math.swt.du (login: ftp) LINEAR DELAY DIFFERENTIAL
More informationu 3 = u 3 (x 1, x 2, x 3 )
Lctur 23: Curvilinar Coordinats (RHB 8.0 It is oftn convnint to work with variabls othr than th Cartsian coordinats x i ( = x, y, z. For xampl in Lctur 5 w mt sphrical polar and cylindrical polar coordinats.
More informationOn spanning trees and cycles of multicolored point sets with few intersections
On spanning trs and cycls of multicolord point sts with fw intrsctions M. Kano, C. Mrino, and J. Urrutia April, 00 Abstract Lt P 1,..., P k b a collction of disjoint point sts in R in gnral position. W
More informationH(curl; Ω) : n v = n
A LOCALLY DIVERGENCE-FREE INTERIOR PENALTY METHOD FOR TWO-DIMENSIONAL CURL-CURL PROBLEMS SUSANNE C. BRENNER, FENGYAN LI, AND LI-YENG SUNG Abstract. An intrior pnalty mthod for crtain two-dimnsional curl-curl
More informationConstruction of asymmetric orthogonal arrays of strength three via a replacement method
isid/ms/26/2 Fbruary, 26 http://www.isid.ac.in/ statmath/indx.php?modul=prprint Construction of asymmtric orthogonal arrays of strngth thr via a rplacmnt mthod Tian-fang Zhang, Qiaoling Dng and Alok Dy
More informationOn the irreducibility of some polynomials in two variables
ACTA ARITHMETICA LXXXII.3 (1997) On th irrducibility of som polynomials in two variabls by B. Brindza and Á. Pintér (Dbrcn) To th mmory of Paul Erdős Lt f(x) and g(y ) b polynomials with intgral cofficints
More information(Upside-Down o Direct Rotation) β - Numbers
Amrican Journal of Mathmatics and Statistics 014, 4(): 58-64 DOI: 10593/jajms0140400 (Upsid-Down o Dirct Rotation) β - Numbrs Ammar Sddiq Mahmood 1, Shukriyah Sabir Ali,* 1 Dpartmnt of Mathmatics, Collg
More informationAddition of angular momentum
Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th
More informationAddition of angular momentum
Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat
More informationSymmetric Interior penalty discontinuous Galerkin methods for elliptic problems in polygons
Symmtric Intrior pnalty discontinuous Galrkin mthods for lliptic problms in polygons F. Müllr and D. Schötzau and Ch. Schwab Rsarch Rport No. 2017-15 March 2017 Sminar für Angwandt Mathmatik Eidgnössisch
More informationDeift/Zhou Steepest descent, Part I
Lctur 9 Dift/Zhou Stpst dscnt, Part I W now focus on th cas of orthogonal polynomials for th wight w(x) = NV (x), V (x) = t x2 2 + x4 4. Sinc th wight dpnds on th paramtr N N w will writ π n,n, a n,n,
More informationHigher-Order Discrete Calculus Methods
Highr-Ordr Discrt Calculus Mthods J. Blair Prot V. Subramanian Ralistic Practical, Cost-ctiv, Physically Accurat Paralll, Moving Msh, Complx Gomtry, Slid 1 Contxt Discrt Calculus Mthods Finit Dirnc Mimtic
More informationDirect Discontinuous Galerkin Method and Its Variations for Second Order Elliptic Equations
DOI 10.1007/s10915-016-0264-z Dirct Discontinuous Galrkin Mthod and Its Variations for Scond Ordr Elliptic Equations Hongying Huang 1,2 Zhng Chn 3 Jin Li 4 Ju Yan 3 Rcivd: 15 Dcmbr 2015 / Rvisd: 10 Jun
More informationEinstein Equations for Tetrad Fields
Apiron, Vol 13, No, Octobr 006 6 Einstin Equations for Ttrad Filds Ali Rıza ŞAHİN, R T L Istanbul (Turky) Evry mtric tnsor can b xprssd by th innr product of ttrad filds W prov that Einstin quations for
More informationarxiv: v1 [math.na] 3 Mar 2016
MATHEMATICS OF COMPUTATION Volum 00, Numbr 0, Pags 000 000 S 0025-5718(XX)0000-0 arxiv:1603.01024v1 [math.na] 3 Mar 2016 RESIDUAL-BASED A POSTERIORI ERROR ESTIMATE FOR INTERFACE PROBLEMS: NONCONFORMING
More informationDG Methods for Elliptic Equations
DG Mthods for Elliptic Equations Part I: Introduction A Prsntation in Profssor C-W Shu s DG Sminar Andras löcknr Tabl of contnts Tabl of contnts 1 Sourcs 1 1 Elliptic Equations 1 11
More informationDynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *
17 nd Intrnational Confrnc on Mchanical Control and Automation (ICMCA 17) ISBN: 978-1-6595-46-8 Dynamic Modlling of Hoisting Stl Wir Rop Da-zhi CAO, Wn-zhng DU, Bao-zhu MA * and Su-bing LIU Xi an High
More informationCOUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM
COUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM Jim Brown Dpartmnt of Mathmatical Scincs, Clmson Univrsity, Clmson, SC 9634, USA jimlb@g.clmson.du Robrt Cass Dpartmnt of Mathmatics,
More informationNONCONFORMING FINITE ELEMENTS FOR REISSNER-MINDLIN PLATES
NONCONFORMING FINITE ELEMENTS FOR REISSNER-MINDLIN PLATES C. CHINOSI Dipartimnto di Scinz Tcnologi Avanzat, Univrsità dl Pimont Orintal, Via Bllini 5/G, 5 Alssandria, Italy E-mail: chinosi@mfn.unipmn.it
More informationAS 5850 Finite Element Analysis
AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form
More informationBasic Polyhedral theory
Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist
More informationA SPLITTING METHOD USING DISCONTINUOUS GALERKIN FOR THE TRANSIENT INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
ESAIM: MAN Vol. 39, N o 6, 005, pp. 5 47 DOI: 0.05/man:005048 ESAIM: Mathmatical Modlling and Numrical Analysis A SPLITTING METHOD USING DISCONTINUOUS GALERKIN FOR THE TRANSIENT INCOMPRESSIBLE NAVIER-STOKES
More informationME469A Numerical Methods for Fluid Mechanics
ME469A Numrical Mthods for Fluid Mchanics Handout #5 Gianluca Iaccarino Finit Volum Mthods Last tim w introducd th FV mthod as a discrtization tchniqu applid to th intgral form of th govrning quations
More informationMiddle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis
Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 43 Introduction to Finit Elmnt Analysis Chaptr 3 Computr Implmntation of D FEM Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt
More informationDiscontinuous Galerkin Approximations for Elliptic Problems
Discontinuous Galrkin Approximations for lliptic Problms F. Brzzi, 1,2 G. Manzini, 2 D. Marini, 1,2 P. Pitra, 2 A. Russo 2 1 Dipartimnto di Matmatica Univrsità di Pavia via Frrata 1 27100 Pavia, Italy
More informationSupplementary Materials
6 Supplmntary Matrials APPENDIX A PHYSICAL INTERPRETATION OF FUEL-RATE-SPEED FUNCTION A truck running on a road with grad/slop θ positiv if moving up and ngativ if moving down facs thr rsistancs: arodynamic
More informationRecall that by Theorems 10.3 and 10.4 together provide us the estimate o(n2 ), S(q) q 9, q=1
Chaptr 11 Th singular sris Rcall that by Thorms 10 and 104 togthr provid us th stimat 9 4 n 2 111 Rn = SnΓ 2 + on2, whr th singular sris Sn was dfind in Chaptr 10 as Sn = q=1 Sq q 9, with Sq = 1 a q gcda,q=1
More informationA ROBUST NONCONFORMING H 2 -ELEMENT
MAHEMAICS OF COMPUAION Volum 70, Numbr 234, Pags 489 505 S 0025-5718(00)01230-8 Articl lctronically publishd on Fbruary 23, 2000 A ROBUS NONCONFORMING H 2 -ELEMEN RYGVE K. NILSSEN, XUE-CHENG AI, AND RAGNAR
More informationCPSC 665 : An Algorithmist s Toolkit Lecture 4 : 21 Jan Linear Programming
CPSC 665 : An Algorithmist s Toolkit Lctur 4 : 21 Jan 2015 Lcturr: Sushant Sachdva Linar Programming Scrib: Rasmus Kyng 1. Introduction An optimization problm rquirs us to find th minimum or maximum) of
More information1 Isoparametric Concept
UNIVERSITY OF CALIFORNIA BERKELEY Dpartmnt of Civil Enginring Spring 06 Structural Enginring, Mchanics and Matrials Profssor: S. Govindj Nots on D isoparamtric lmnts Isoparamtric Concpt Th isoparamtric
More informationcycle that does not cross any edges (including its own), then it has at least
W prov th following thorm: Thorm If a K n is drawn in th plan in such a way that it has a hamiltonian cycl that dos not cross any dgs (including its own, thn it has at last n ( 4 48 π + O(n crossings Th
More informationA UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE VOLUME METHODS FOR THE STOKES EQUATIONS
A UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE VOLUME METHODS FOR THE STOKES EQUATIONS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. In tis papr, t autors stablisd a unifid framwork for driving and
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 43, No. 4, pp. 1728 1749 c 25 Socity for Industrial and Applid Mathmatics SUPERCONVERGENCE OF THE VELOCITY IN MIMETIC FINITE DIFFERENCE METHODS ON QUADRILATERALS M. BERNDT, K.
More informationABEL TYPE THEOREMS FOR THE WAVELET TRANSFORM THROUGH THE QUASIASYMPTOTIC BOUNDEDNESS
Novi Sad J. Math. Vol. 45, No. 1, 2015, 201-206 ABEL TYPE THEOREMS FOR THE WAVELET TRANSFORM THROUGH THE QUASIASYMPTOTIC BOUNDEDNESS Mirjana Vuković 1 and Ivana Zubac 2 Ddicatd to Acadmician Bogoljub Stanković
More informationDifference -Analytical Method of The One-Dimensional Convection-Diffusion Equation
Diffrnc -Analytical Mthod of Th On-Dimnsional Convction-Diffusion Equation Dalabav Umurdin Dpartmnt mathmatic modlling, Univrsity of orld Economy and Diplomacy, Uzbistan Abstract. An analytical diffrncing
More informationAdrian Lew, Patrizio Neff, Deborah Sulsky, and Michael Ortiz
AMRX Applid Mathmatics Rsarch Xprss 004, No. 3 Optimal V stimats for a Discontinuous Galrkin Mthod for Linar lasticity Adrian Lw, Patrizio Nff, Dborah Sulsky, and Michal Ortiz 1 Introduction Discontinuous
More informationQuasi-Classical States of the Simple Harmonic Oscillator
Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats
More informationON THE DISTRIBUTION OF THE ELLIPTIC SUBSET SUM GENERATOR OF PSEUDORANDOM NUMBERS
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 8 2008, #A3 ON THE DISTRIBUTION OF THE ELLIPTIC SUBSET SUM GENERATOR OF PSEUDORANDOM NUMBERS Edwin D. El-Mahassni Dpartmnt of Computing, Macquari
More informationProperties of Phase Space Wavefunctions and Eigenvalue Equation of Momentum Dispersion Operator
Proprtis of Phas Spac Wavfunctions and Eignvalu Equation of Momntum Disprsion Oprator Ravo Tokiniaina Ranaivoson 1, Raolina Andriambololona 2, Hanitriarivo Rakotoson 3 raolinasp@yahoo.fr 1 ;jacqulinraolina@hotmail.com
More information2.3 Matrix Formulation
23 Matrix Formulation 43 A mor complicatd xampl ariss for a nonlinar systm of diffrntial quations Considr th following xampl Exampl 23 x y + x( x 2 y 2 y x + y( x 2 y 2 (233 Transforming to polar coordinats,
More informationThat is, we start with a general matrix: And end with a simpler matrix:
DIAGON ALIZATION OF THE STR ESS TEN SOR INTRO DUCTIO N By th us of Cauchy s thorm w ar abl to rduc th numbr of strss componnts in th strss tnsor to only nin valus. An additional simplification of th strss
More informationRECOVERY-BASED ERROR ESTIMATORS FOR INTERFACE PROBLEMS: MIXED AND NONCONFORMING FINITE ELEMENTS
SIAM J. NUMER. ANAL. Vol. 48 No. 1 pp. 30 52 c 2010 Socity for Industrial Applid Mathmatics RECOVERY-BASED ERROR ESTIMATORS FOR INTERFACE PROBLEMS: MIXED AND NONCONFORMING FINITE ELEMENTS ZHIQIANG CAI
More informationBINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES
BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES DONALD M. DAVIS Abstract. If p is a prim (implicit in notation and n a positiv intgr, lt ν(n dnot th xponnt of p in n, and U(n n/p ν(n, th unit
More informationApplication of Vague Soft Sets in students evaluation
Availabl onlin at www.plagiarsarchlibrary.com Advancs in Applid Scinc Rsarch, 0, (6):48-43 ISSN: 0976-860 CODEN (USA): AASRFC Application of Vagu Soft Sts in studnts valuation B. Chtia*and P. K. Das Dpartmnt
More informationBoundary layers for cellular flows at high Péclet numbers
Boundary layrs for cllular flows at high Péclt numbrs Alxi Novikov Dpartmnt of Mathmatics, Pnnsylvania Stat Univrsity, Univrsity Park, PA 16802 Gorg Papanicolaou Dpartmnt of Mathmatics, Stanford Univrsity,
More informationCONVERGENCE RESULTS FOR THE VECTOR PENALTY-PROJECTION AND TWO-STEP ARTIFICIAL COMPRESSIBILITY METHODS. Philippe Angot.
DISCRETE AND CONTINUOUS doi:.3934/dcdsb..7.383 DYNAMICAL SYSTEMS SERIES B Volum 7, Numbr 5, July pp. 383 45 CONVERGENCE RESULTS FOR THE VECTOR PENALTY-PROJECTION AND TWO-STEP ARTIFICIAL COMPRESSIBILITY
More informationInference Methods for Stochastic Volatility Models
Intrnational Mathmatical Forum, Vol 8, 03, no 8, 369-375 Infrnc Mthods for Stochastic Volatility Modls Maddalna Cavicchioli Cá Foscari Univrsity of Vnic Advancd School of Economics Cannargio 3, Vnic, Italy
More informationCramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter
WHEN THE CRAMÉR-RAO INEQUALITY PROVIDES NO INFORMATION STEVEN J. MILLER Abstract. W invstigat a on-paramtr family of probability dnsitis (rlatd to th Parto distribution, which dscribs many natural phnomna)
More informationSECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.
SETION 6. 57 6. Evaluation of Dfinit Intgrals Exampl 6.6 W hav usd dfinit intgrals to valuat contour intgrals. It may com as a surpris to larn that contour intgrals and rsidus can b usd to valuat crtain
More informationChapter 10. The singular integral Introducing S(n) and J(n)
Chaptr Th singular intgral Our aim in this chaptr is to rplac th functions S (n) and J (n) by mor convnint xprssions; ths will b calld th singular sris S(n) and th singular intgral J(n). This will b don
More informationSymmetric centrosymmetric matrix vector multiplication
Linar Algbra and its Applications 320 (2000) 193 198 www.lsvir.com/locat/laa Symmtric cntrosymmtric matrix vctor multiplication A. Mlman 1 Dpartmnt of Mathmatics, Univrsity of San Francisco, San Francisco,
More informationA SIMPLE FINITE ELEMENT METHOD OF THE CAUCHY PROBLEM FOR POISSON EQUATION. We consider the Cauchy problem for Poisson equation
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volum 14, Numbr 4-5, Pags 591 603 c 2017 Institut for Scintific Computing and Information A SIMPLE FINITE ELEMENT METHOD OF THE CAUCHY PROBLEM FOR
More informationFinding low cost TSP and 2-matching solutions using certain half integer subtour vertices
Finding low cost TSP and 2-matching solutions using crtain half intgr subtour vrtics Sylvia Boyd and Robrt Carr Novmbr 996 Introduction Givn th complt graph K n = (V, E) on n nods with dg costs c R E,
More informationMiddle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis
Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME Introduction to Finit Elmnt Analysis Chaptr 5 Two-Dimnsional Formulation Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt
More informationThe Equitable Dominating Graph
Intrnational Journal of Enginring Rsarch and Tchnology. ISSN 0974-3154 Volum 8, Numbr 1 (015), pp. 35-4 Intrnational Rsarch Publication Hous http://www.irphous.com Th Equitabl Dominating Graph P.N. Vinay
More informationMutually Independent Hamiltonian Cycles of Pancake Networks
Mutually Indpndnt Hamiltonian Cycls of Pancak Ntworks Chng-Kuan Lin Dpartmnt of Mathmatics National Cntral Univrsity, Chung-Li, Taiwan 00, R O C discipl@ms0urlcomtw Hua-Min Huang Dpartmnt of Mathmatics
More informationHigher order derivatives
Robrto s Nots on Diffrntial Calculus Chaptr 4: Basic diffrntiation ruls Sction 7 Highr ordr drivativs What you nd to know alrady: Basic diffrntiation ruls. What you can larn hr: How to rpat th procss of
More informationBoundary layers for cellular flows at high Péclet numbers
Boundary layrs for cllular flows at high Péclt numbrs Alxi Novikov Gorg Papanicolaou Lnya Ryzhik January 0, 004 Abstract W analyz bhavior of solutions of th stady advction-diffusion problms in boundd domains
More informationThe van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012
Th van dr Waals intraction D. E. Sopr 2 Univrsity of Orgon 20 pril 202 Th van dr Waals intraction is discussd in Chaptr 5 of J. J. Sakurai, Modrn Quantum Mchanics. Hr I tak a look at it in a littl mor
More information1 Minimum Cut Problem
CS 6 Lctur 6 Min Cut and argr s Algorithm Scribs: Png Hui How (05), Virginia Dat: May 4, 06 Minimum Cut Problm Today, w introduc th minimum cut problm. This problm has many motivations, on of which coms
More informationFourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.
Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform
More informationSearch sequence databases 3 10/25/2016
Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an
More informationUNIFIED ERROR ANALYSIS
UNIFIED ERROR ANALYSIS LONG CHEN CONTENTS 1. Lax Equivalnc Torm 1 2. Abstract Error Analysis 2 3. Application: Finit Diffrnc Mtods 4 4. Application: Finit Elmnt Mtods 4 5. Application: Conforming Discrtization
More informationHydrogen Atom and One Electron Ions
Hydrogn Atom and On Elctron Ions Th Schrödingr quation for this two-body problm starts out th sam as th gnral two-body Schrödingr quation. First w sparat out th motion of th cntr of mass. Th intrnal potntial
More informationDiscontinuous Galerkin approximation of flows in fractured porous media
MOX-Rport No. 22/2016 Discontinuous Galrkin approximation of flows in fracturd porous mdia Antonitti, P.F.; Facciola', C.; Russo, A.;Vrani, M. MOX, Dipartimnto di Matmatica Politcnico di Milano, Via Bonardi
More informationPROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS
Intrnational Journal Of Advanc Rsarch In Scinc And Enginring http://www.ijars.com IJARSE, Vol. No., Issu No., Fbruary, 013 ISSN-319-8354(E) PROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS 1 Dharmndra
More informationMORTAR ADAPTIVITY IN MIXED METHODS FOR FLOW IN POROUS MEDIA
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volum 2, Numbr 3, Pags 241 282 c 25 Institut for Scintific Computing and Information MORTAR ADAPTIVITY IN MIXED METHODS FOR FLOW IN POROUS MEDIA
More informationSquare of Hamilton cycle in a random graph
Squar of Hamilton cycl in a random graph Andrzj Dudk Alan Friz Jun 28, 2016 Abstract W show that p = n is a sharp thrshold for th random graph G n,p to contain th squar of a Hamilton cycl. This improvs
More informationBrief Introduction to Statistical Mechanics
Brif Introduction to Statistical Mchanics. Purpos: Ths nots ar intndd to provid a vry quick introduction to Statistical Mchanics. Th fild is of cours far mor vast than could b containd in ths fw pags.
More informationCS 361 Meeting 12 10/3/18
CS 36 Mting 2 /3/8 Announcmnts. Homwork 4 is du Friday. If Friday is Mountain Day, homwork should b turnd in at my offic or th dpartmnt offic bfor 4. 2. Homwork 5 will b availabl ovr th wknd. 3. Our midtrm
More informationReparameterization and Adaptive Quadrature for the Isogeometric Discontinuous Galerkin Method
Rparamtrization and Adaptiv Quadratur for th Isogomtric Discontinuous Galrkin Mthod Agns Silr, Brt Jüttlr 2 Doctoral Program Computational Mathmatics 2 Institut of Applid Gomtry Johanns Kplr Univrsity
More informationEXISTENCE OF POSITIVE ENTIRE RADIAL SOLUTIONS TO A (k 1, k 2 )-HESSIAN SYSTEMS WITH CONVECTION TERMS
Elctronic Journal of Diffrntial Equations, Vol. 26 (26, No. 272, pp. 8. ISSN: 72-669. URL: http://jd.math.txstat.du or http://jd.math.unt.du EXISTENCE OF POSITIVE ENTIRE RADIAL SOLUTIONS TO A (k, k 2 -HESSIAN
More informationSection 6.1. Question: 2. Let H be a subgroup of a group G. Then H operates on G by left multiplication. Describe the orbits for this operation.
MAT 444 H Barclo Spring 004 Homwork 6 Solutions Sction 6 Lt H b a subgroup of a group G Thn H oprats on G by lft multiplication Dscrib th orbits for this opration Th orbits of G ar th right costs of H
More informationVSMN30 FINITA ELEMENTMETODEN - DUGGA
VSMN3 FINITA ELEMENTMETODEN - DUGGA 1-11-6 kl. 8.-1. Maximum points: 4, Rquird points to pass: Assistanc: CALFEM manual and calculator Problm 1 ( 8p ) 8 7 6 5 y 4 1. m x 1 3 1. m Th isotropic two-dimnsional
More informationNumerical considerations regarding the simulation of an aircraft in the approaching phase for landing
INCAS BULLETIN, Volum, Numbr 1/ 1 Numrical considrations rgarding th simulation of an aircraft in th approaching phas for landing Ionl Cristinl IORGA ionliorga@yahoo.com Univrsity of Craiova, Alxandru
More information