arxiv: v1 [math.na] 3 Mar 2016

Transcription

1 MATHEMATICS OF COMPUTATION Volum 00, Numbr 0, Pags S (XX) arxiv: v1 [math.na] 3 Mar 2016 RESIDUAL-BASED A POSTERIORI ERROR ESTIMATE FOR INTERFACE PROBLEMS: NONCONFORMING LINEAR ELEMENTS ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG Abstract. In this papr, w study a modifid rsidual-basd a postriori rror stimator for th nonconforming linar finit lmnt approximation to th intrfac problm. Th rliability of th stimator is analyzd by a nw and dirct approach without using th Hlmholtz dcomposition. It is provd that th stimator is rliabl with constant indpndnt of th jump of diffusion cofficints across th intrfacs, without th assumption that th diffusion cofficint is quasi-monoton. Numrical rsults for on tst problm with intrscting intrfacs ar also prsntd. 1. Introduction During th past dcad, th construction, analysis, and implmntation of robust a postriori rror stimators for various finit lmnt approximations to partial diffrntial quations with paramtrs hav bn on of th focuss of rsarch in th fild of th a postriori rror stimation. For th lliptic intrfac problm, various robust stimators hav bn constructd, analyzd, and implmntd (s,.g., [4, 23, 22, 8, 9, 11, 26, 12] for conforminglmnts, [1, 20, 10] for nonconforminglmnts, [10] for mixd lmnts, and [7] for discontinuous lmnts). Th robustnss for rsidual basd stimators in th rliability bound is stablishd thortically undr th assumption of th quasi-monoton distribution of th diffusion cofficints, s [4] for mor dtails. Howvr, numrical rsults by many rsarchrs including ours strongly suggst that thos stimators ar robust vn whn th diffusion cofficints ar not quasi-monoton. In this papr, w provid a thortical vidnc for th nonconforming linar lmnt without th quasi-monoton assumption. On of th ky stps in obtaining th robust rliability bound of classical rsidual basd stimator is to construct a modifid Clémnt-typ intrpolation oprator satisfying spcific approximation and stability proprtis in th nrgy norm (s [4] for dtails). For th conforming linar lmnt, th dgrs of frdom ar th nodal valus at vrtics of triangls. Th nodal valu of th modifid Clémnttyp intrpolation is dfind by th avrag valu of th function ovr connctd lmnts whos corrsponding diffusion cofficints ar th gratst. Undr th quasi-monoton assumption, Brnardi and Vrfürth [4] wr abl to stablish th rquird proprtis of th intrpolation oprator to guarant th robust rliability 2010 Mathmatics Subjct Classification. Primary 65N30. This work was supportd in part by th National Scinc Foundation undr grants DMS and DMS , th Purdu Rsarch Foundation, and th Rsarch Grants Council of th Hong ong SAR, China undr th GRF Projct No , CityU c XXXX Amrican Mathmatical Socity

2 2 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG bound. A ky advantag for th nonconforming linar lmnt is that its dgrs of frdom ar nodal valus at th middl points of dgs of triangls and that ach middl point is shard by at most two triangls. Hnc, w ar abl to construct a modifid Clémnt-typ intrpolation satisfying th dsird proprtis without th quasi-monotonicity assumption (s Sction 4). Th a postriori rror stimation for th nonconforming lmnts has bn studid by many rsarchrs. Du to th lack of th rrorquation, Dari, Duran, Padra, and Vampa [14] stablishd th rliability bound of th rsidual-basd rror stimator for th Poisson quation through th Hlmholtz dcomposition of th tru rror. Thir analysis is widly usd by othr rsarchrs (s,.g., [13, 5, 1, 6]), and th Hlmholtz dcomposition bcoms a ncssary tool for obtaining th rliability bound for th nonconforming lmnts. This approach has also bn applid to th mixd finit lmnt mthod [21] and discontinuous Garlrkin finit lmnt mthod [3, 2, 7]. It is obvious that application of thir analysis to th intrfac problm will lad to th sam distribution assumption as th conforming lmnts in [4]. Ainsworth [1] constructd an quilibratd stimator without using th Clémnt typ intrpolation but th rror bounds dpnd on th jump of diffusion constants. Dspit th main trnd of using Hlmholtz dcomposition in th nonconforming finit lmnt analysis, thr ar svral othr intrsting paprs that approachd diffrntly. Hopp and Wohlmuth [18] constructd two a postriori rror stimators by using th hirarchical basis undr th saturation assumption. Schiwck [24] constructd a two-sidd bound of th nrgy rror using th analysis of conforming cas with som simpl additional argumnts. Nvrthlss, conforming Clémnt typ intrpolation was applid in that papr hnc again impos th assumption of quasi-monotonicity. Th purpos of this papr is to prsnt a nw and dirct analysis, which dos not involv th Hlmholtz dcomposition, for stimating th rliability bound with th aim of rmoving th quasi-monoton assumption. To do so, our analysis maks us of (a) our nwly dvlopd rror quation for th nonconforming finit lmnt approximation in [7] and (b) th structur of th nonconforming lmnts. Combining with our obsrvation on th modifid Clémnt-typ intrpolation for th nonconforming lmnts, w ar abl to bound both th lmnt rsiduals and th numrical flux jumps uniformly without th quasi-monotonicity assumption. Unfortunatly, w ar unabl to do th sam for th numrical dg solution jump. As an altrnativ, w modify th dg solution jump at lmnts whr th quasimonotonicity assumption is not satisfid. Th modifid stimator is provd to b rliabl with constant indpndnt of th jump of th diffusion cofficints across intrfacs without th quasi-monotonicity assumption. By using th standard argumnt (s,.g., [4]), w also stablish local fficincy bounds uniformly with rspct to th jump of th diffusion cofficint. This robustnss is obtaind for th standard (modifid) indicators without (with) th quasi-monotonicity assumption. Nvrthlss, numrical rsults prsntd in Sction 7 for a bnchmark tst problm sms to suggst that th modifid indicator gnrats a bttr msh than th standard indicator. Existing rsidual basd stimators consist of th lmnt rsidual, th dg flux jump, and th dg tangntial drivativ jump du to th Hlmholtz dcomposition. As a by-product of our dirct approach (s (2.9)), th rsidual basd stimators

3 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS 3 to b studid in this papr rplac th dg tangntial drivativ jump by th dg solution jump. Evn though thy ar quivalnt in two dimnsions, numrical rsult shows that our stimator is mor accurat than th xisting stimators (s Figur 6). Th outlin of th papr is as follows. Th intrfac problm and its nonconforming finit lmnt approximation ar introducd in Sction 2 as wll as th L 2 rprsntation of th tru rror in th (brokn) nrgy norm. Th standard and modifid indicators and stimators ar prsntd in Sction 3. Th modifid Clémnt-typ intrpolation oprator is dfind and its approximation proprtis ar provd in Sction 4. Robust local fficincy and global rliability bounds ar stablishd in Sctions 5 and 6, rspctivly. Finally, w provid som numrical rsults in Sction Nonconforming Linar Finit Elmnt Approximation to Intrfac Problm 2.1. Intrfac Problm. For simplicity of th prsntation, w considr only two dimnsions. Extnsion of th rsults in this papr to thr dimnsions is straightforward. Lt Ω b a boundd, opn, connctd subst of R 2 with a Lipschitz continuous boundary Ω. Dnot by n = (n 1, n 2 ) t th outward unit vctor normal to th boundary. W partition th boundary of th domain Ω into two opn substs Γ D and Γ N such that Ω = Γ D Γ N and that Γ D Γ N =. For simplicity, w assum that Γ D is not mpty (i.., ms(γ D ) 0). Considr th following lliptic intrfac problm (2.1) (α(x) u) = f in Ω with boundary conditions (2.2) u = g D on Γ D and n (α u) = g N on Γ N, whr and ar th divrgnc and gradint oprators, rspctivly; f, g D, and g N ar givn scalar-valud functions; and th diffusion cofficint α > 0 is picwis constant with rspct to a partition of th domain Ω = n Ω i. Hr th subdomain Ω i is opn and polygonal. Th jump of th α across intrfacs (subdomain boundaris) ar possibly vry larg. For simplicity, assum that f, g D, and g N ar picwis linar functions. W us th standard notations and dfinitions for th Sobolv spacs H s (Ω) and H s ( Ω) for s 0. Th standard associatd innr products ar dnotd by (, ) s,ω and (, ) s, Ω, and thir rspctiv norms ar dnotd by s,ω and s, Ω. (W omit th subscript Ω from th innr product and norm dsignation whn thr is no risk of confusion.) For s = 0, H s (Ω) coincids with L 2 (Ω). In this cas, th innr product and norm will b dnotd by and (, ), rspctivly. Lt H 1 g,d (Ω) := {v H1 (Ω) : v = g D on Γ D }. Th corrsponding variational formulation of problm (2.1)-(2.2) is to find u Hg,D 1 (Ω) such that (2.3) a(u,v) = f(v), v H0,D(Ω), 1 whr th bilinar and linar forms ar dfind by a(u,v) = (α(x) u, v) Ω and f(v) = (f,v) Ω +(g N,v) ΓN. i=1

4 4 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG 2.2. Nonconforming Linar Finit Elmnt Approximation. Lt T h b a triangulation of th domain Ω. Assum that T h is rgular; i.., for all T h, thr xist a positiv constant κ such that h κρ, whr h dnots th diamtr ofth lmnt and ρ th diamtr ofth largst circl that may b inscribd in. Not that th assumption of th msh rgularity dos not xclud highly, locally rfind mshs. Lt N h = N h I N h D N h N and E h = E h I E h D E h N, whr N h I (Eh I ) is th st of all intrior vrtics (dgs) in T h, and N h D (Eh D ) and Nh N (E h N ) ar th rspctiv sts of all vrtics (dgs) on Γ D and Γ N. For ach E h, dnot by m th mid-point of th dg. Furthrmor, assum that intrfacs F = { Ω i Ω j : i,j = 1, n} do not cut through any lmnt T h. Lt P k () b th spac of polynomials of dgr lss than or qual to k on th lmnt. Dnot th conforming picwis linar finit lmnt spac associatd with th triangulation T h by and its subst by U c = {v H 1 (Ω) : v P 1 (), T h } U c g,d = {v U c : v ΓD = g D }. Dnot th nonconforming picwis linar finit lmnt spac, i.., th Crouzix- Raviart lmnt [17], associatd with th triangulation T h by U nc = {v L 2 (Ω) : v P 1 (), T h, and v is continuous at m, E h I } and its subst by Lt U nc g,d = {v U nc : v(m ) = g D (m ), E h D}. H 1 (T h ) = {v L 2 (Ω) : v H 1 (), T h }. For any v, w H 1 (T h ), dnot th (brokn) bilinar form by a h (v, w) = T h (α v, w) and th (brokn) nrgy norm by v Ω = a h (v, v) = α 1/2 v 2 0, T h Th nonconforming finit lmnt approximation is to find u h Ug,D nc such that (2.4) a h (u h,v) = f(v), v U nc 0,D. 1/2.

5 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS L 2 Rprsntation of th Error. For ach dg E h, dnot by h th lngth of ; dnot by n a unit vctor normal to. Whn ED h Eh N, dnot by + th boundary lmnt with th dg, and assum that n is th unit outward normal vctor of. For any EI h, lt + and b th two lmnts sharing th common dg assuming that α + α + α α, and that n coincids with th unit outward normal vctor of +. Dnot by v + and v, rspctivly, th tracs of th doubl valud function v ovr rstrictd on + and. For any v H1 (T h ), dnot th normal flux jump ovr dg E h by (α v n ) + (α v n ), Eh I, α v n := 0, E D h, (α v n ) g N, EN h, and th valu jump ovr dg E h by v + v, Eh I, v := v g D, ED h, 0, EN h. Th arithmtic avrag ovr dg E h is dnotd by 1 ( v + {v} := 2 +v ), E h I, v, ED h Eh N. A simpl calculation lads to th following idntity: (2.5) uv = {u} v + u {v}, E h I. For any v U0,D nc, it is wll known that th following orthogonality proprty holds (2.6) v ds = 0, EI h EN h and vds = 0, ED. h Lt u and u h b th solutions of (2.3) and (2.4), rspctivly. It is shown in [7] that (2.7) a h (u,v h ) = f(v h )+ (α u n ) v h ds+ (α u n )v h ds, v h U nc E h I Dnot th tru rror by E h D E = u u h. Diffrnc of (2.7) and (2.4) yilds th following rror quation: (2.8) a h (E,v h ) = (α u n ) v h ds+ (α u n )v h ds, v h U nc E h I E h D Introducing th lmnt rsidual, th numrical flux jump, and th numrical solution jump r = ( f + (α u h ) ), T h, j σ, = α u h n and j u, = u h, E h, 0,D. 0,D.

6 6 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG rspctivly, thn th tru rror in th (brokn) nrgy norm may b xprssd in trms of thos quantitis. Lmma 2.1. Lt E h U0,D nc b an intrpolation of E, thn w hav th following L 2 rprsntation of th rror E in th (brokn) nrgy norm: (2.9) a h (E,E) = (r,e E h ) j σ, {E E h }ds {α E n }j u, ds. T h E h I Eh N E h I Eh D Proof. Firstnotthat{α u h n } isaconstantforvry E h. Thorthogonality in (2.6) lads to (2.10) {α u h n } E h ds = 0, EI h and (α u h n )E h ds = 0, ED. h It follows from intgration by parts, (2.5), th continuitis of th normal componnt of th flux α u and th solution u, and (2.10) that a h (E, E E h ) = (α E, (E E h )) T h = (r, E E h ) + (α E n )(E E h ) ds T h + E h D Eh N E h I (α E n )(E E h )ds = (r, E E h ) + α E n {E E h }ds T h EI h + {α E n } ( E E h ) ds+ (α E n )(E E h )ds E h I = (r, E E h ) T h E h I Eh N E h D Eh N j σ, {E E h }ds E h I Eh D (α u n ) E h ds (α u n )E h ds, E h I E h D which, togthr with th rror quation in (2.8) with v h = E h, yilds a h (E,E) = a h (E, E E h )+a h (E,E h ) = (r, E E h ) j σ, {E E h }ds T h E h I Eh N E h I Eh D {α E n }j u, ds {α E n }j u, ds. This complts th proof of th lmma.

7 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS 7 3. Indicator and Estimator In this sction, basd on th L 2 rprsntation of th tru rror in th nrgy norm in Lmma 2.1, w first introduc th standard indicator and th corrsponding stimator. Our standard stimator consists of th usual lmnt rsidual and dg flux jump plus th dg solution jump that rplacs th dg tangntial drivativ jump of xisting rsidual basd stimators. Sinc th robustnss of th rliability bound of stimators was stablishd undr th quasi-monotonicity assumption on th distribution of th diffusion cofficint, to avoid such an assumption w introduc a nw indicator and th corrsponding stimator by modifying th dg solution jump at lmnts whr th quasi-monotonicity assumption fails. For any T h, dnot by N h and Eh, rspctivly, th sts of thr vrtics and thr dgs of. Dnot th rspctiv indicators of lmnt rsidual, dg flux jump, and dg solution jump by η 2 R f, = h2 α r 2 0,, η 2 J σ, = E h Eh I and η 2 J u, = E h Eh I h 2α + j σ, 2 0, + E h Eh N α j u, 2 0, + 2h E h Eh D h α j σ, 2 0,, α h j u, 2 0,. Thn th standard indicator associatd with T h is dfind by ( 1/2, (3.1) η = ηr 2 f, +η2 J σ, +η2 J u,) and th standard stimator by (3.2) η = T h η 2 Rmark 3.1. Instad of using th dg tangntial drivativ jump as xisting rsidualbasdrrorstimators, th indicatorη and th rsulting stimatorη abovmploy th dg solution jump η Ju,. In two dimnsions, th dg solution jump is quivalnt to th tangntial drivativ jump (s (5.4)). Nvrthlss, numrical rsults for a tst problm show that our stimator is mor accurat than th xisting stimators using th dg tangntial drivativ jump. By th standard argumnt [4], it is shown in Sction 5 that th indicator η is fficint uniformly with rspct to th jump of th diffusion cofficint. By using th Hlmholtz dcomposition and th modifid Clémnt-typ intrpolation, on can also prov that th stimator η is rliabl. Morovr, th rliability constant is indpndnt of th jump of α(x) providd that th distribution of α(x) is quasimonoton [23]. In ordr to rmov this assumption, w prsnt a nw analysis for stimating th rliability bound without using th Hlmholtz dcomposition. Th analysis will mak us of th structur of th nonconforming lmnt in twodimnsions, and it nabls us to bound both th lmnt rsidual and th numrical flux jump uniformly without th quasi-monotonicity. Unfortunatly, w ar unabl to do th sam for th numrical solution jump. As an altrnativ, w modify th 1/2.

8 8 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG indicator of th numrical solution jump at lmnts whr th quasi-monotonicity is not satisfid. To this nd, for ach vrtx z N h, dnot by ωz h and Eh z, rspctivly, th sts of all lmnts T h and all dgs E h having z as a common vrtx. Lt ω z h = { ωz h : α = max α } ω h ωz h z b th st of all lmnts in ωz h such that th corrsponding diffusion cofficints ar th gratst. For any intrfac intrscting point z N h, th vrtx patch ωz h is calld quasi-monoton (s [23]) if for ach ωz, h thr xists a subst ŵz, h of ωz h such that th union of lmnts in ŵz, h is a Lipschitz domain and that Dnot by if z N h \N h D, thn {} ωh z ŵh z, and α α, ŵ h z, ; ifz N h D, thn ŵh z,, ( ŵ h z, ) Γ D andα α, ŵ h z,. N M = {z N h : ω h z is not quasi-monoton } th st of all intrfac intrscting points whos vrtx patchs ar not quasimonoton. For ach lmnt T h, subdivid it into four sub-triangls by conncting thr mid-points of dgs of, and dnot by T h/2 th rfind triangulation. Lt N h/2 b th sts of all vrtics basd on T h/2. N h/2 = N h/2 I N h/2 D Nh/2 N and E h/2 = E h/2 I E h/2 D Eh/2 N whr N h/2 I (E h/2 I ), N h/2 D (Eh/2 D ), and Nh/2 N (Nh/2 N ) ar th sts of all intrior vrtics (dgs) of T h/2, all boundary vrtics (dgs) on Γ D and Γ N, rspctivly. Lt { } U h/2,c g,d = v H 1 (Ω) : v T P 1 (T), T T h/2 and v ΓD = g D, which is th continuous picwis linar finit lmnt spac associatd with th triangulation T h/2. Nxt, w introduc an intrpolation oprator, I h/2 : Ug,D nc Uh/2,c g,d, from th nonconforming picwis linar finit lmnt spac on T h to th conforming picwis linar finit lmnt spac on T h/2. For a givn v Ug,D nc, th nodal valus of I h/2 v U h/2,c g,d ar dfind as follows: (i) st (ii) st (iii) st (I h/2 v)(z) = g D (z), z N h D ; (I h/2 v)(m ) = v(m ), E h ; (I h/2 v)(z) = v z (z), z (N h I N h N), whr z is chosn to b on lmnt in ω h z. For ach vrtx z N h, dnot by ωz h/2 th sts of all lmnts T T h/2 having z as a common vrtx. For lmnt T h with at last on vrtx in N M, th indicator of th numrical solution jump η Ju, is modifid as follows:

9 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS 9 η 2 J u, = z N h \NM Ez h/2 E h/2 Eh/2 I α 4h j u, 2 0, + Ez h/2 E h/2 Eh/2 D α j u, 2 0, 2h + z N h NM α 2h Ih/2 u h u h 2 0, T,z, whr T,z = ωz h/2. Thn th modifid indicator is dfind as follows: (3.3) η = ( η 2 R f, +η2 J σ, + η2 J u,) 1/2, if N h N M η, othrwis. Th corrsponding modifid stimator is thn givn by ( ) 1/2 (3.4) η = η 2. T h Rmark 3.2. In th cas that N M = ; i.., th distribution of th diffusion cofficint is quasi-monoton, thn η Ju, = η Ju, for all T h and hnc η = η and η = η. 4. Th Modifid Clémnt-typ Intrpolation In this sction, following th ida in [4, 15], w introduc th modifid Clémnttyp intrpolation oprator for th non-conforming linar lmnt and stablish its approximation proprtis. Dnot by vdx = ω 1 vdx mas(ω) ω th manvalu ofagivnfunction v on agivnmasurablst ω in R 2 with positiv 2-dimnsional Lbsgu masur mas(ω). With this convntion, st { ffl π (v) = vdx, + Eh I, ffl vdx, ED h Eh N. Th modifid Clémnt intrpolation oprator I h : H 1 (T h ) U nc is dfind by (4.1) I h (v) = (π v)φ, E h whr φ is th nodal basis function of U nc which taks valu 1 at m and taks 0 at mid-points of othr dgs. For any T h, lt b th union of lmnts in T h sharing an dg with. For any E h, lt b th union of lmnts in T h having th common dg.

10 10 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG Lmma 4.1. For any function v H 1 (T h ), thn th modifid Clémnt intrpolation satisfis th following approximation proprtis: (4.2) v I h v 0, h v α 1/2 + ( ) α 1/2 v 0,, T h h and (4.3) E h v + π v ( ) 1/2 h 0, α + v 0, +, Eh. Hr and thraftr, w us th a b notation to indicat that a cb for a furthr not spficifid constant c, which dpnds only on th shap rgularity of T h but not on th data of th undrlying problms, in particular, th jump of th diffusion cofficint. Unlik th modifid Clémnt intrpolation for th conforming lmnts, thr is an xtra jump trm in th approximation proprty in (4.2) which is du to th discontinuity of th function v across th dgs of. Proof. For any T h, sinc th nodal basis functions form a partition of th unity, th triangl inquality givs v I h v 0, = φ (v π v) 0, v π v 0,. E h E h Hnc, to show th validity of (4.2), it suffics to prov that ( (4.4) v π v 0, h ( ) α 1/2 ) v α 1/2 + v 0,, E h. h Sinc th st contains only two lmnts for all EI h, it is obvious that = + or. If = +, thn (4.4) is a dirct consqunc of th Poincaré inquality: v π v 0, = v vdx 0, h α 1/2 v. In th cas that =, th triangl and th Poincaré inqualitis imply v π v 0, v vdx 0, + ( h α 1/2 v +h 1/2 ( h α 1/2 v +h 1/2 vdx v vdx v ) + 0, + vdx + vdx 0, ( v + vdx + + vdx v + ) v 0, 0, + ) v 0,. Nxt, w bound th thr trms abov. It follows from th trac thorm and th Poincaré inquality that h 1/2 vdx v 0, vdx v 0, +h vdx v 1, h α 1/2 v. Similarly, w hav (4.5) h 1/2 + vdx v + 0, h α 1/2 v + +.

11 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS 11 Not that α = α α +, combining abov thr inqualitis givs ( v π v 0, h ( ) α 1/2 ) v α 1/2 + v 0,, h which provs th validity of (4.4) whn =. Whn ED h Eh N, (4.4) is a dirct consqunc of th Poincaréinquality. This complts th proof of (4.4) and hnc (4.2). (4.3) is a dirct consqunc of (4.5). This complts th proof of th lmma. 5. Local Efficincy Bound In this sction, w stablish local fficincy bounds for th indicators η and η dfind, rspctivly, in (3.1) and (3.3). Both th bounds ar indpndnt of th jump of th diffusion cofficint, this robustnss is obtaind without (with) th quasi-monotonicity assumption (s [4]) for th indicator η ( η ). Nvrthlss, numrical rsults prsntd in Sction 7 for a bnchmark tst problm sms to suggst that th modifid indicator η gnrats a bttr msh than th standard indicator η. By using local dg and lmnt bubbl functions, ψ and ψ (s [25] for thir dfinitions and proprtis), it is a common practic to obtain th local fficincy bound for th rsidual-basd a postriori rror stimator. By proprly wighting trms in th indicator by th diffusion cofficint, on can show that th local fficincy bound is robust without th quasi-monotonicity assumption. For th convninc of radrs, w only sktch th proof in th following thorm. Thorm 5.1. (Local Efficincy) Assuming that u H 1+ǫ (Ω) and u h ar th solutions of (2.3) and (2.4), rspctivly, thn th indicator η satisfis th following local fficincy bound: (5.1) η E, T h. Proof. For any T h, it follows from th proprtis of ψ, intgration by parts, and th Cauchy-Schwarz inquality that r 2 0, (f + (α u h ))r ψ dx = α (u u h ) (r ψ )dx α 1/2 u u h r ψ 1, α 1/2 h 1 u u h r 0,, which implis (5.2) r 0, α1/2 u u h, T h. h For any EI h, by using th proprtis of ψ, intgration by parts, th Cauchy- Schwartz inquality, and (5.2), w hav j σ, 2 0, α u h n j σ, ψ ds = (α E n)j σ, ψ ds = ( (α E) (j σ, ψ )dx ) r j σ, ψ dx ( ) α + 1/2 E j σ, 0,. h

12 12 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG Togthr with a similar bound for ED h Eh N, it implis ( ) α + 1/2 (5.3) j σ, 0, E, E h. h For any E h, lt n = (n 1, n 2 ) t, thn τ = ( n 2, n 1 ) t is th unit vctor tangnt to th dg. Dnot by j τ, = u h τ th jump of th tangntial drivativ of th numrical solution u h along th dg. By th continuity of u h at th midpoint m, w hav (5.4) j u, 0, = 1 h j τ, 0,, EI. h 12 ( ) t v For a scalar-valud function v, dnot by v = y, v th formal adjoint x oprator of th curl oprator in two-dimnsions. For any EI h, it follows from th proprtis of ψ, intgration by parts, and th Cauchy-Schwartz inquality that j τ, 2 0, uh τ jτ, ψ ds = ( E τ)j τ, ψ ds = α 1/2 E (j τ, ψ )dx α 1/2 E j τ, ψ 1, E h 1/2 j τ, 0, ( α h 1/2 ) E j τ, 0,, which, togthr with (5.4) and a similar bound for ED h Eh N, yilds ( ) 1/2 h (5.5) j u, 0, E, E h. α Now, th fficincy bound in (5.1) is a dirct consqunc of th bounds in (5.2), (5.3), and (5.5). This complts th proof of th thorm. Nxt, w stablish th local fficincy bound for th modifid indicator η. To this nd, lt T h b an lmnt having at last on vrtx z N M. From th lmnt ωz h to th lmnt z ω z h ωz, h whr z is chosn in th dfinition of th oprator I h/2 in Sction 3, thr ar at most two possibl paths (clockwis or countr-clockwis) in ωz. h Dnot by ω h/2,z th union of lmnts of ωh/2 z on on of th paths such that th maximum of th ratio btwn α and th diffusion cofficints ovr lmnts on that path is th smallst. Lt α. α T C,z max T ω h/2,z Rmark 5.2. If ω h z is quasi-monoton, thn C,z = 1 for all ω h z. Lmma 5.3. Lt u h b th solution of (2.4), and lt T h b an lmnt having at last on vrtx z N M and lt T,z = ω h/2 z. Lt I h/2 b th intrpolation oprator dfind in Sction 3 with z ω h/2,z dscribd abov. Thn w hav (5.6) α I h/2 u h u h 2 0, T h,z 2C,z E h/2 z α h uh 2 0,.

13 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS 13 Proof. Without loss of gnrality, w prov (5.6) only whn z NI h. In th cas whn z ND h Nh N, (5.6) may b provd in a similar fashion. To this nd, assum that thr ar k lmnts in ω h/2,z and dnot ths lmnts by T i (i = 1,..., k) starting from T 1 = T,z along th path and nding at T k = z ωz h/2. Dnot by u i and α i th rstrictions of u h (z) and α on T i, rspctivly, and by i th dg btwn T i 1 and T i. A dirct calculation givs that α h I h/2 u h u h 2 0, T,z = 2 3 α 1(u k u 1 ) 2 and that α i u h 2 0, h i = 1 i 3 α ui i u i 1 2, i = 2,...,k. Togthr with th triangl inquality, w hav α h I h/2 u h u h 2 0, T,z = 2 3 α 1(u k u 1 ) k i=2 α 1 α i α i (u i u i 1 ) 2 2C,z k i=2 α i h i u h 2 0, i 2C,z E h/2 z E h/2 I This complts th proof of (5.6) and, hnc, th lmma. α h u h 2 0,. Thorm 5.4. (Local Efficincy) Assuming that u H 1+ǫ (Ω) and u h ar th solutions of (2.3) and (2.4), rspctivly. Thn th modifid indicator η satisfis th following local fficincy bound: (5.7) η max C,z E, T h. z N h Proof. (5.7) is a dirct consqunc of (5.1) and (5.6). Lt η Ju, = ( α h 6. Global Rliability Bound ) 1/2 I h/2u h u h and η Ju = T h α I h/2 u h u h 2 h Lmma 6.1. Lt u h b th solution of (2.4) and I h/2 b th intrpolation oprator dfind in Sction 3, thn th jump of th numrical solution has th following uppr bound: (6.1) {α E n }j u, ds η Ju E Ω. E h I Eh D Proof. Sinc th intgral ovr dg sgmnt on th lft-hand sid of inquality (6.1) may b only rgardd as th duality pair btwn H δ 1/2 () and H 1/2 δ () for an arbitrarily small δ > 0, w ar not abl to bound this intgral dirctly. To ovrcom this difficulty, w xprss thm in trms of intgrals along th boundary of lmnts. To this nd, first not that I h/2 u h = 0 and α u n = 0, E h I. 1/2.

14 14 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG By (2.5) and th fact that I h/2 u h = g D on Γ D, w hav {α E n }j u, ds EI h Eh D = {α E n } I h/2 u h u h ds+ (α E n )(g D u h )ds EI h ED h = (α E n)(i h/2 u h u h )ds α E n {I h/2 u h u h }ds T h EI h (α E n )(I h/2 u h u h )ds+ (α E n )(g D u h )ds = ED h Eh N (α E n)(i h/2 u h u h )ds+ E h D T h EI h Eh N I 1 +I 2. j σ, {I h/2 u h u h }ds Th I 1 may b boundd abov by using th dfinition of th dual norm, th trac thorm (s,.g., [7]), th invrs inquality, and (5.2) as follows: I 1 T h α E n 1/2, Ih/2 u h u h 1/2, α 1/2 T h ( α E 0, +h r 0, ) η Ju, (6.2) η Ju E Ω. To bound th I 2, first not that j σ, I h/2 u h u h ds = 0, EI h, which is a consqunc of th orthogonality proprty in (2.6) and th facts that j σ, is a constant and that I h/2 u h = 0 for all E h I. Hnc, j σ, {I h/2 u h u h }ds = = j σ, {I h/2 u h u h }ds+ 1 j σ, I h/2 u h u h ds 2 j σ, (I h/2 u h u h + )ds, Eh I.

15 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS 15 Now, it follows from th Cauchy-Schwartz inquality and (5.3) that I 2 = j σ, (I h/2 u h u h + )ds EI h Eh N E h I Eh N E h I Eh N α + I h/2 u h u h + h 2 0, α + I h/2 u h u h + h 2 0, 1/2 1/2 h α + EI h Eh N E Ω j σ, 2 0, 1/2 (6.3) η Ju E Ω. (6.1) is thn a consqunc of (6.2) and (6.3). This complts th proof of th lmma. Dnot by η Ju = ( T h η 2 J u, with th solution jump. Thn η 2 J u = z N h \NM E h/2 z α 2h j u, 2 0, + ) 1/2 th part ofth modifid stimator associatd z N h NM T ω h/2 z Lmma 6.2. Th η Ju is boundd abov by th η Ju ; i.., (6.4) η Ju η Ju. Proof. Sinc u h I h/2 u h vanishs on all boundary dgs of wz h/2 hav η Ju = T h α h I h/2 u h u h 2 0, = z N h T w h/2 z α T 4h T Ih/2 u h u h 2 0, T. for all z N h, w α T 2h T I h/2 u h u h 2 0, T. To prov th validity of (6.4), it suffics to show that for all z N h \N M, α T I h/2 u h u h 2 0, T h T E h/2 z α h u h 2 0,, T ωh/2 z. This may b provd in a similar fashion as (5.6) with th fact that C,z = 1 for all z N h \N M. This complts th proof of th lmma. Thorm 6.3. (Global Rliability) Lt u and u h b th solutions of (2.3) and (2.4), rspctivly. Thn th stimator η satisfis th following global rliability bound: (6.5) E Ω η.

16 16 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG Proof. Lt I h b th modifid Clémnt intrpolation oprator dfind in Sction 4. Thn (2.9) with E h = I h E bcoms a h (E,E) = (r,e I h E) j σ, {E I h E}ds T h (6.6) E h I Eh D I 1 +I 2 +I 3. E h I Eh N {α E n }j u, ds Th first trm in (6.6) may b boundd by th Cauchy-Schwartz inquality, Lmma 4.1, and (5.5) as follows: I 1 r 0, E I h E 0, T h (6.7) T h η Rf, T h η 2 R f, E 0, + ( ) α 1/2 u h 0, 1/2 E Ω. E h To bound th scond trm in (6.6), first notic that h E I h E = u h +I h E, E h I. Sinc u h +I h E U nc and th th fact that j σ, is a constant for all E h, (2.6) yilds j σ, E E h ds = 0, EI h. Hnc, (6.8) h E} ds+ {E I 1 2 E I h E ds = (E I h E) + ds = (E + π E)ds, for all EI h. Th last quality coms from th proprty of th nonconforming nodal basis functions: ffl i φ j = δ ij. It thn follows from (6.8), th Cauchy-Schwartz inquality, and Lmma 4.1 that I 2 = j σ, (E + π E)ds j σ, 0, E + π E 0, (6.9) E h I Eh N E h I Eh N ( h α + E h I Eh N ) 1/2 j σ, 0, E 0, + T h η 2 J σ, 1/2 E Ω. Now, (6.5) is a dirct consqunc of (6.6), (6.7), (6.9), and Lmmas 6.1 and 6.2. This complts th proof of th thorm.

17 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS Numrical Exprimnts In this sction, w rport som numrical rsults for an intrfac problm with intrscting intrfacs usd by many authors (s,.g., [19, 7]), which is considrd as a bnchmark tst problm. Lt Ω = ( 1,1) 2 and u(r,θ) = r β µ(θ) in th polar coordinats at th origin with cos((π/2 σ)β) cos((θ π/2+ρ)β) if 0 θ π/2, cos(ρβ) cos((θ π +σ)β) if π/2 θ π, µ(θ) = cos(σβ) cos((θ π ρ)β) if π θ 3π/2, cos((π/2 ρ)β) cos((θ 3π/2 σ)β) if 3π/2 θ 2π, whr σ and ρ ar numbrs. Th function u(r,θ) satisfis th intrfac problm in (2.1) with Γ N =, f = 0, and { R in (0,1) 2 ( 1,0) 2, α(x) = 1 in Ω\([0,1] 2 [ 1,0] 2 ). Th numbrs β, R, σ, and ρ satisfy som nonlinar rlations. For xampl, whn β = 0.1, thn R , ρ = π/4, and σ Not that whn β = 0.1, this is a difficult problm for computation. Rmark 7.1. This problm dos not satisfy Hypothsis 2.7 in [4] as th quasimonotonicity is not satisfid about th origin. Startd with a coars triangulation, a squnc of mshs is gnratd by using a standard adaptiv mshing algorithm that adopts th L 2 stratgy: (i) mark lmnts whos indicators ar among th first 20 prcnt of th nrgy norm of th total rror, and (ii) rfin th markd triangls by bisction. Th stopping critria rl-rr := u u u Ω tol u Ω is usd, and numrical rsults with tol = 0.1 ar rportd. Figur 1. msh gnratd by η. Figur 2. msh gnratd by η.

18 18 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG Enrgy rror Enrgy rror u u h A N 0.5 ta tild ta u u h A N 0.5 ta tild ta Figur 3. rror and stimators on msh gnratd by η Figur 4. rror and stimators on msh gnratd by η. Mshs gnratd by th standard and modifid indicators, η and η, ar dpictd rspctivly in Figurs 1 and 2. Both th rfinmnts ar cntrd at th origin. Thr ar11974and 5524lmnts in th rspctivfigurs1and 2. Hnc, this tst problm suggsts that th modifid indicator gnrats a much bttr msh than th standard indicator vn though th local fficincy bound of th modifid indicator dpnds on th jump of th diffusion cofficint. Th comparisons of th tru rror and th stimators η and η on th mshs gnratd by th standard and modifid indicators ar shown in Figurs 3 and 4, rspctivly. Th slop of th log(dof)- log(rror) for both th stimators on both th mshs ar vry clos to 1/2, which indicats th optimal dcay of th rror with rspct to th numbr of unknowns. Th fficincy indx is dfind by ff-indx := stimator u u h Ω. Th fficincy indics for th η and η ar about and on th msh gnratd by η and about and on th msh gnratd by η, rspctivly. Enrgy rror 10 1 u u h A 10 2 N 0.5 ta with tangntial jump ta with solution jump Figur 5. msh gnratd by η with solution jump Figur 6. rror and stimators on msh gnratd by η

19 RESIDUAL-BASED A POSTERIORI ESTIMATE FOR NONCONFORMING ELEMENTS 19 Our numrical rsults also show that our standard stimator using th dg solution jump is mor accurat than th xisting stimators using th dg tangntial jump. To illustrat this fact, w prsnt numrical rsults for a tst problm [16]: a Poisson quation dfind on th L-shapd domain Ω = ( 1, 1) 2 \[0, 1] [ 1, 0] with th following xact solution ( ) 2θ +π u(x, θ) = r 2/3 sin, θ [0, 3π/2]. 3 Th stopping critria is st as tol Th fficincy indics in th final stp ar and for th rspctiv stimators with th dg solution and tangntial drivativ jumps. This indicats that our standard stimator is mor accurat than th xisting rsidual stimator (s Figur 6). Rfrncs 1. M. Ainsworth, Robust a postriori rror stimation for nonconforming finit lmnt approximation, SIAM J. Numr. Anal., 42:6 (2005), M. Ainsworth, A postriori rror stimation for discontinuous Galrkin finit lmnt approximation, SIAM J. Numr. Anal., 45:4 (2007), P. Bckr, P. Hansbo, and M. G. Larson, Enrgy norm a postriori rror stimation for discontinuous Galrkin mthod, Comput. Mth. Appl. Mch. Engrg., 192 (2003), C. Brnardi and R. Vrfürth, Adaptiv finit lmnt mthods for lliptic quations with non-smooth cofficints, Numr. Math., 85 (2000), C. Carstnsn, S. Bartls, and S. Jansch, A postriori rror stimats for nonconforming finit lmnt mthods, Numr. Math., 92 (2002), C. Carstnsn, J. Hu, and A. Orlando, Framwork for th a postriori rror analysis of th nonconforming finit lmnts, SIAM J. Numr. Anal., 45:1 (2007), Z. Cai, X. Y, and S. Zhang, Discontinuous Galrkin finit lmnt mthods for intrfac problms: a priori and a postriori rror stimations, SIAM J. Numr. Anal., 49:5 (2011), Z. Cai and S. Zhang, Rcovry-basd rror stimator for intrfac problms: conforming linar lmnts, SIAM J. Numr. Anal., 47:3 (2009), Z. Cai and S. Zhang, Flux rcovry and a postriori rror stimators: conforming lmnts for scalar lliptic quations, SIAM J. Numr. Anal., 48:2 (2010), Z. Cai and S. Zhang, Rcovry-basd rror stimator for intrfac problms: mixd and nonconforming lmnts, SIAM J. Numr. Anal., 48:1 (2010), Z. Cai and S. Zhang, Robust rsidual- and rcovry a postriori rror stimators for intrfac problms with flux jumps, Numr. Mthods for PDEs, 28:2 (2012), Z. Cai and S. Zhang, Robust quilibratd rsidual rror stimator for diffusion problms: conforming lmnts, SIAM J. Numr. Anal., 50:1 (2012), E. Dari, R. Duran, and C. Padra, Error stimations for nonconforming finit lmnt approximations of th Stoks problm, Math. Comp., 64 (1995), E. Dari, R. Duran, C. Padra, and V. Vampa, A postriori rror stimators for nonconforming finit lmnt mthods, RAIRO Modél Anal. Numér., 30:4 (1996), M. Dryja, M. V. Sarkis, and O. B. Widlund, Multilvl Schwartz mthod for lliptic problms with discontinuous in thr dimnsions, Numr. Math., 72 (1996), H. Elman, D. Silvstr, and A. Wathn, Finit Elmnts and Fast Itrativ Solvrs: With Applications in Incomprssibl Fluid Dynamics, Numr. Math. Sci. Comput., Oxford Univrsity Prss, Oxford, U, V. Girault and P. A. Raviart, Finit Elmnt Mthods for Navir-Stoks Equations, Springr-Vrlag, Brlin, R. H. W. Hopp and B. Wohlmuth, Elmnt-orintd and dg-orintd local rror stimators for nonconforming finit lmnt mthods, RAIRO Modél. Math. Anal. Numér., 30:2 (1996), R. B. llogg, On th Poisson quation with intrscting intrfacs, Appl. Anal., 4 (1975),

20 20 ZHIQIANG CAI, CUIYU HE, AND SHUN ZHANG 20.. Y. im, A postriori rror analysis for locally consrvativ mixd mthods, Math. Comp., 76 (2007), C. Lovadina and R. Stnbrg, Enrgy norm a postriori rror stimats for mixd finit lmnt mthods, Math. Comp., 75 (2006), R. Luc and B. I. Wohlmuth, A local a postriori rror stimator basd on quilibratd fluxs, SIAM J. Numr. Anal., 42:4 (2004), M. Ptzoldt, A postriori rror stimators for lliptic quations with discontinuous cofficints, Adv. Comput. Math., 16 (2002), F. Schiwck, A postriori rror stimats with post-procssing for nonconforming finit lmnts, ESAIM Math. Mod. Numr. Anal., 36:3 (2002), R. Vrfürth, A Postriori Error Estimation Tchniqus for Finit Elmnt Mthods, Oxford Univrsity Prss, Oxford, Unitd ingdom, M. Vohralík, Guarantd and fully robust a postriori rror stimats for conforming discrtizations of diffusion problms with discontinuous cofficints, J. Sci. Comput., 46:3 (2011), Dpartmnt of Mathmatics, Purdu Univrsity, 150 N. Univrsity Strt, Wst Lafaytt, IN addrss: caiz@purdu.du Dpartmnt of Mathmatics, Purdu Univrsity, 150 N. Univrsity Strt, Wst Lafaytt, IN addrss: h75@purdu.du Dpartmnt of Mathmatics, City Univrsity of Hong ong, Hong ong addrss: shun.zhang@cityu.du.hk