Discontinuous Galerkin approximation of flows in fractured porous media
|
|
- Toby Joseph
- 5 years ago
- Views:
Transcription
1 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 Milano (Italy)
2 Discontinuous Galrkin approximation of flows in fracturd porous mdia Paola F. Antonitti, Chiara Facciolà, Alssandro Russo # and Marco Vrani May 15, 2016 MOX- Laboratory for Modling and Scintific Computing Dipartimnto di Matmatica Politcnico di Milano Piazza Lonardo da Vinci 32, Milano, Italy paola.antonitti@polimi.it, chiara.facciola@polimi.it, marco.vrani@polimi.it # Dipartimnto di Matmatica Applicazioni Univrsità dgli Studi di Milano-Bicocca Via Cozzi 55, Milano, Italy alssandro.russo@unimib.it Abstract W prsnt a numrical approximation of Darcy s flow through a fracturd porous mdium which mploys discontinuous Galrkin mthods. For simplicity, w considr th cas of a singl fractur rprsntd by a (d 1)-dimnsional intrfac btwn two d-dimnsional subdomains, d = 2, 3. W propos a discontinuous Galrkin finit lmnt approximation for th flow in th porous matrix which is coupld with a conforming finit lmnt schm for th flow in th fractur. Suitabl (physically consistnt) coupling conditions complt th modl. W thortically analys th rsulting formulation and prov its wll-posdnss. Morovr, w driv optimal a priori rror stimats in a suitabl (msh-dpndnt) nrgy norm and w prsnt two-dimnsional numrical xprimnts assssing thir validity. Introduction Modlling flows in fracturd porous mdia has rcivd incrasing attntion in th past dcads, bing fundamntal for addrssing many 1
3 2 nvironmntal and nrgy problms, such as watr rsourcs managmnt, oil migration tracmnt, isolation of radioactiv wast, ground watr contamination. In ths applications th flow is strongly influncd by th prsnc of fracturs, which can act as prfrntial paths (whn thir prmability is highr than that of th surrounding mdium), or as barrirs for th flow (whn thy ar filld with low prmabl matrial). A fractur is typically dfind as a rgion charactrizd by a small aprtur compard to its lngth and th siz of th domain and with a diffrnt porous structur than th surrounding mdium. Th task of ffctivly modlling th intraction btwn th systm of fracturs and th porous matrix is particularly challnging. In th following, lt us brifly commnt on a popular modlling choic to handl such a problm, s.g [27, 5, 20]. To rduc th complxity of th problm, a common modlling choic consists in trating fracturs as (d 1)-dimnsional intrfacs btwn d-dimnsional porous matrics, d = 2, 3. Th dvlopmnt of this kind of rducd modls, which can b justifid in cas of fracturs with vry small width, has bn addrssd for singl-phas flows in svral works, s.g. [2, 1, 27, 23]. In this papr w adopt th prspctiv of th singl fractur modl dscribd in [27]. A first vrsion of this modl has bn introducd in [2] and [1] undr th assumption of larg prmability in th fractur. In [27] th modl has bn furthr gnralisd to handl also fracturs with low prmability. On th othr hand, th flow in th porous mdium is assumd to b govrnd by Darcy s law and a suitabl rducd vrsion of this law is formulatd on th surfac modlling th fractur. Physically consistnt coupling conditions ar addd to account for th xchang of fluid btwn th fractur and th porous mdium. Th xtnsion of such a coupld modl to th cas of two-phas flow has bn addrssd in [24] and [26], whil a totally immrsd fractur has bn considrd in [3]. Various numrical mthods hav bn mployd in th litratur for th approximation of th rsulting coupld bulk-fractur flow. Roughly spaking, thy can b classifid dpnding on th condition (matching or non-matching) btwn th bulk and th fractur mshs. A traditional approach combining mixd finit lmnts and bulk mshs conforming to th fractur msh was adoptd for xampl in [2, 23, 27]. Th us of non-matching grids coupld with th Xtndd Finit Elmnt Mthod (XFEM) has bn proposd in [24, 20]. Mor rcntly, an approximation basd on th us of conforming polygonal mshs and Mimtic Finit Diffrncs (MFD) has bn xplord in [5]. W also mntion a promising framwork to trat flows in systms of fractur ntworks introducd in [12, 13, 14, 15, 11]. Th aim of this papr is to mploy discontinuous Garlrkin (DG) finit
4 3 lmnts to discrtiz th coupld bulk-fractur problm stmming from th modlling of flows in fracturd porous mdia. Th inhritd flxibility of DG mthods in handling arbitrarily shapd, non-ncssarly matching, grids rprsnts th idal stting to handl such kind of problms that typically fatur a high-lvl of gomtrical complxity. Discontinuous Galrkin mthods wr first introducd in th arly 1970s (s for xampl [28, 22, 8, 31, 6]) as a tchniqu to numrically solv partial diffrntial quations. Thy hav bn succssfully dvlopd and applid to hyprbolic, lliptic and parabolic problms arising from a wid rang of applications: various xampls can b found, for xampl, in [7, 9, 18, 19, 16, 29, 25, 21]. Th primary motivations for th us of DG mthods ar th nhancd flxibility affordd by discontinuous lmnts [6] and th possibility of handling mshs mad of arbitrarily shapd lmnts and with hanging nods. Morovr, th local natur of th trial spac allows lmntwis variabl polynomial ordrs which nabls mor accurat approximation of solutions which vary in charactr from on part of th domain to anothr. W rfr to [7] for a unifid prsntation and analysis on DG mthods for lliptic problms. Mor spcifically, th choic of DG mthods for addrssing th problm of th flow in a fracturd porous mdium ariss quit spontanously in viw of th discontinuous natur of th solution at th matrix-fractur intrfac. Howvr, this is not th only motivation to mploy DG mthods in this spcific contxt. Indd, our diffrntial modl is basd on th primal form of th Darcy s quations for both th bulk and fractur flows, which ar coupld with suitabl conditions at th intrfac. Ths coupling conditions can b naturally formulatd using jump and avrag oprators, so that DG mthods turn out to b a vry natural and powrful tool for fficintly handling th coupling of th two problms, which is indd naturally mbddd, in wak form, in th variational formulation. In this papr w propos a discrtization which combins a DG approximation for th problm in th bulk with a conforming finit lmnt approximation in th fractur. Th us of conforming finit lmnts to discrtiz th quations in th fractur is mad just for th sak of simplicity, othr discrtization tchniqus can b mployd and our approach is as gnral to tak into account straightforwardly also such cass. W analys th rsulting mthod and prov a priori rror stimats, which w numrically tst in a two-dimnsional stting. Th papr is structurd as follows. In Sction 1 w introduc th govrning quations for th coupld problm. Th problm is thn writtn in a wak form in Sction 2, whr w also prov its wll-posdnss. In Sction 3 w
5 4 introduc th DG discrtization of th coupld problm and w show som tchnical rsults ndd to prov its wll-posdnss. In Sction 4 w driv a priori rror stimats in a suitabl (msh-dpndnt) norm, whil in Sction 5 w prsnt two-dimnsional numrical xprimnts assssing th validity of th thortical rror stimats. 1 Modl problm Throughout th papr w will mploy th following notation. For an opn, boundd domain D R d, d = 2, 3, w dnot by H s (D) th standard Sobolv spac of ordr s, for a ral numbr s 0. Th usual norm on H s (D) is dnotd by H s (D) and th usual sminorm by H s (D). Furthrmor, w will dnot by P k (D) th spac of polynomials of dgr lss than or qual to k 1 on D. Throughout th papr th symbol (and ) will signify that th inqualitis hold up to multiplicativ constants which ar indpndnt of th discrtization paramtr. In th following w prsnt th govrning quations for our modl, which is a variant of th modl drivd in [27]. Th flow of an incomprssibl fluid through a fracturd d-dimnsional porous mdium, d = 2, 3, can b dscribd by th following thr ingrdints: 1. th govrning quations for th flow in th porous mdium; 2. th govrning quations for th flow in th fracturs; 3. a st of physically consistnt conditions which coupl th problms in th bulk and fracturs along thir intrfacs. For simplicity, w will assum that thr is a uniqu fractur in th porous mdium and that th fractur cuts th domain xactly into two disjoint connctd subrgions (s Figur 1 for a two-dimnsional xampl), following th approach of [5] and [20]. Th xtnsion to multipl fracturs can b tratd analogously, whil th cas of an immrsd fractur is mor complx to b analysd [3] and will b th subjct of futur rsarch. Mor prcisly, lt Ω R d, d = 2, 3, b an opn, boundd, convx polygonal/polyhdral domain rprsnting th porous matrix. W suppos that th fractur is a (d 1)-dimnsional C manifold Γ R d 1, d = 2, 3, whos masur is uniformly boundd (i.., Γ = O(1)), and assum that Γ sparats Ω into two connctd subdomains, which ar disjoint, i.., Ω \ Γ = Ω 1 Ω 2 with Ω 1 Ω 2 =. For i = 1, 2, w dnot by γ i th part of boundary of Ω i shard
6 5 Γ γ 2 Ω 2 n Γ γ 1 Ω 1 Figur 1: Th subdomains Ω 1 and Ω 2 sparatd by th fractur Γ considrd as an intrfac. with th boundary of Ω, i.., γ i = Ω i Ω. W dnot by n i, i = 1, 2 th unit normal vctor to Γ pointing outwards from Ω i and, for a (rgular nough) scalar-valud function v and a (rgular nough) vctor-valud function τ, w dfin th standard jump and avrag oprators across Γ as {v} = 1 2 (v 1 + v 2 ) v = v 1 n 1 + v 2 n 2, {τ } = 1 2 (τ 1 + τ 2 ) τ = τ 1 n 1 + τ 2 n 2, (1) whr th subscript i = 1, 2 dnots th rstriction to th subdomain Ω i. Morovr w dnot by n Γ th normal unit vctor on Γ with a fixd orintation from Ω 1 to Ω 2, so that w hav n Γ = n 1 = n Govrning quations According to th abov discussion, w suppos that th flow in th bulk is govrnd by Darcy s law. Lt K = K(x) R d d b th bulk prmability tnsor, which satisfis th following rgularity assumptions: (i) K is a symmtric, positiv dfinit tnsor whos ntris ar boundd, picwis continuous ral-valud functions; (ii) K is uniformly boundd by blow and abov, i.., x T x x T Kx x T x x R d. (2)
7 6 Givn a function f L 2 (Ω) rprsnting a sourc trm and g H 1/2 ( Ω), th motion of a incomprssibl fluid in ach domain Ω i, i = 1, 2, with prssur p i is dscribd by: (K i p i ) = f i in Ω i, i = 1, 2, (3) p i = g i on γ i, i = 1, 2. (4) Hr w hav dnotd by K i and f i, th rstrictions of K and f to Ω i, i = 1, 2, rspctivly, and by g i th rstriction of g to γ i, i = 1, 2 (for simplicity, w hav imposd Dirichlt boundary conditions on both γ 1 and γ 2 ). Th scond ingrdint for th modl is rprsntd by th govrning quations for th fractur flow. In our modl th fractur is tratd as a (d 1)-dimnsional manifold immrsd in a d-dimnsional objct. If w assum that th fracturs ar filld by a porous mdium with diffrnt porosity and prmability than th surroundings, Darcy s law can b usd also for modlling th flow along th fracturs [10]. Th rducd modl is thn obtaind through a procss of avraging across th fractur: in th bginning th fractur is assumd to b a d-dimnsional subdomain of Ω, that sparats it into two disjoint subdomains. Thn Darcy s quations ar writtn on th fractur in th normal and tangntial componnts and th tangntial componnt is intgratd along th thicknss l Γ = l Γ (x) > 0 of th fractur domain, which is typically som ordrs of magnitud smallr than th siz of th domain. W rfr to [27] for a rigourous drivation of th rducd mathmatical modl. Not that in [27] this avraging procss is carrid out for th flow quations writtn in mixd form. Hr, w rwrit th modl in primal form. Th fractur flow is thn charactrizd by th fractur prmability tnsor K Γ, which is assumd to satisfy th sam rgularity assumptions as thos satisfid by th bulk prmability K and to hav a block-diagonal structur of th form [ ] K n K Γ = Γ 0 0 KΓ τ, (5) whn writtn in its normal and tangntial componnts. Hr K τ Γ R(d 1) (d 1) is a positiv dfinit tnsor (it rducs to a positiv numbr for d = 2) that rprsnts th tangntial componnt of th prmability of th fractur. Lt us assum that f Γ L 2 (Γ). Stting Γ = Γ Ω, and dnoting by
8 7 p Γ th fractur prssur, th govrning quations for th fractur flow rad τ (K τ Γl Γ τ p Γ ) = f + K p in Γ, (6) p Γ = g Γ on Γ, (7) whr g Γ H 1/2 ( Γ) and τ and τ dnot th tangntial gradint and divrgnc oprators, rspctivly. Equation (6) rprsnts Darcy s law in th dirction tangntial to th fractur, whr a sourc trm K p is introducd to tak into account th contribution of th subdomain flows to th fractur flow [27]. For th sak of simplicity, w impos Dirichlt boundary conditions at th boundary Γ of th fractur Γ. Finally, following [27], w provid th intrfac conditions to coupl problms (3)-(4) and (6)-(7). Lt ξ b a positiv ral numbr, ξ 1 2, that will b chosn latr on. Th coupling conditions ar givn by whr 2{K p} n Γ = β Γ (p 1 p 2 ) on Γ, (8) K p = α Γ ({p} p Γ ) on Γ, (9) β Γ = 1 2η Γ, α Γ = 2 η Γ (2ξ 1) (10) and η Γ = l Γ KΓ n, KΓ n bing th normal componnt of th fractur prmability tnsor, s (5). Not that th coupling conditions ar formulatd mploying jump and avrag oprators. This turns out to b convnint for mploying DG mthods in th discrtization. In conclusion, th coupld modl problm rads: (K i p i ) = f i in Ω i, i = 1, 2, p i = g i on γ i, i = 1, 2, τ (KΓl τ Γ τ p Γ ) = f + K p in Γ, p Γ = g Γ on Γ, 2{K p} n Γ = β Γ (p 1 p 2 ) on Γ, K p = α Γ ({p} p Γ ) on Γ. (11) Not that th introduction of th paramtr ξ yilds a family of modls, s [27] for mor dtails.
9 8 2 Wak formulation and its wll-posdnss In this sction w prsnt a wak formulation of our modl problm (11) and prov its wll-posdnss. For th sak of simplicity w will assum that homognous Dirichlt boundary conditions ar imposd for both th bulk and fractur problms. Th xtnsion to th gnral cas is straightforward. W introduc th following spacs V b = {p = (p 1, p 2 ) V b 1 V b 2 }, V Γ = H 1 0 (Γ) H s (Γ), (12) whr w dfin, for i = 1, 2 and s 1, Vi b H0,γ 1 i (Ω i ) = {q H 1 (Ω i ) s.t. q γi = 0}. = H s (Ω i ) H 1 0,γ i (Ω i ), with Nxt w introduc th bilinar forms A b : V b V b R, A Γ : V Γ V Γ R and I : (V b V Γ ) (V b V Γ ) R, dfind as follows 2 A b (p, q) = K p i q i, i=1 Ω i A Γ (p Γ, q Γ ) = KΓl τ Γ p Γ q Γ, Γ I((p, p Γ ), (q, q Γ )) = β Γ p q + α Γ ({p} p Γ )({q} q Γ ). Γ Clarly, th bilinar forms A b (, ) and A Γ (, ) tak into account th problms in th bulk and in th fractur, rspctivly, whil I(, ) taks into account th intrfac conditions (8)-(9). W also introduc th linar functionals L b : V b R and L Γ : V Γ R dfind as 2 L b (q) = fq i, i=1 Ω i L Γ (q Γ ) = fq Γ, which rprsnt th sourc trms in th bulk and fractur, rspctivly. With th abov notation, th wak formulation of our modl problm rads as follows: Find (p, p Γ ) V b V Γ such that, for all (q, q Γ ) V b V Γ Γ Γ A ((p, p Γ ), (q, q Γ )) = L(q, q Γ ), (13) whr A : (V b V Γ ) (V b V Γ ) R is dfind as th sum of th bilinar forms just introducd: A ((p, p Γ ), (q, q Γ )) = A b (p, q) + A Γ (p Γ, q Γ ) + I((p, p Γ ), (q, q Γ )), (14)
10 9 and th linar oprator L : V b V Γ R is dfind as L(q, q Γ ) = L b (q) + L Γ (q Γ ). (15) Nxt, w show that formulation (13) is wll-posd. introduc th following norm on V b V Γ : To this aim w (q, q Γ ) 2 E = 2 i=1 K 1/2 i q i 2 L 2 (Ω i ) + (Kτ Γl Γ ) 1/2 q Γ 2 L 2 (Γ) + β Γ q 2 L 2 (Γ) + α Γ {q} q Γ 2 L 2 (Γ). (16) 2 This is clarly a norm if α Γ 0. Sinc α Γ = η Γ (2ξ 1), from now on, w will assum that ξ > 1/2. W rmark that th sam condition on th paramtr ξ has bn found also in [27] and [5]. Thorm 2.1. Lt ξ > 1/2. Thn, problm (13) is wll-posd. Proof. W show that A(, ) is continuous and corciv on V b V Γ quippd with th norm (16), as wll as L( ) is continuous on V b V Γ with rspct to th sam norm. Thn, xistnc and uniqunss of th solution, as wll as linar dpndnc on th data, follow dirctly from Lax-Milgram s lmma. Corcivity is straightforward, as w clarly hav that A((q, q Γ ), (q, q Γ )) = (q, q Γ ) 2 E for any (q, q Γ ) V b V Γ. On th othr hand, continuity is a dirct consqunc of Cauchy-Schwarz inquality, whil continuity of L( ) on V b V Γ is guarantd by th rgularity of th forcing trm f. 3 Numrical discrtization In this sction w prsnt a numrical discrtization of our problm which combins a Discontinuous Galrkin approximation for th problm in th bulk with a conforming finit lmnt approximation in th fractur (s Rmark 2 blow). DG mthods rsult to b vry convnint for handling th discontinuity of th bulk prssur across th fractur, as wll as th coupling of th bulk-fractur problms, which has bn formulatd using jump and avrag oprators. As a rsult, w can mploy th standard tools offrd by DG mthods to prov th wll-posdnss of our discrt mthod (Proposition 3.3). W start with th introduction of som usful notation. W considr a family of shap-rgular, not-ncssarily matching, simplicial mshs T h which ar alignd with th fractur Γ, so that any triangl E T h cannot b cut by Γ.
11 10 Not that, sinc Ω 1 and Ω 2 ar disjoint, ach lmnt E blongs xactly to on of th two subdomains. Clarly ach msh T h inducs a subdivision of th fractur Γ into dgs (dg mans fac for d = 3), that w will dnot by Γ h. Morovr, th st of all dgs of th dcomposition T h is dnotd by E h and w hav E h = E I h EB h Γ h, whr Eh B is th st of boundary dgs and EI h is th st of intrior dgs not blonging to th fractur. Not that, sinc th prsnc of hanging nods is prmittd, an intrior dg is dfind as th (non-mpty) intrsction of th boundaris of two nighbouring lmnts of T h. For ach lmnt E T h, w dnot by h E its diamtr and w st h = max E Th h E. W mak th assumption that, for all Eh I, if = E+ E, it holds h E + h and h E h, whr h is th diamtr of. Finally, givn an lmnt E T h, for any dg E w dfin n as th unit normal vctor on that points outsid E. W can thn dfin th standard jump and avrag oprators across an dg E h for (rgular nough) scalar and vctor-valud functions similarly to (1). Givn a partition T h of th domain, w dnot by H s (T h ), s 0, th standard brokn Sobolv spac. For s = 0 w writ L 2 (Ω) in plac of H 0 (T h ). With th aim of building a DG-conforming finit lmnt approximation, w choos to st th discrt problm in th finit-dimnsional spacs V b h = {q h L 2 (Ω) : q h E P k (E) E T h }, V Γ h = {qγ h C0 (Γ) : q Γ h P k () Γ h }. W rmark that th polynomial dgr in th bulk and fractur discrt spacs just dfind can b chosn indpndntly. Hr, for th sak of simplicity, w dnot both by k. Nxt, w introduc th bilinar forms A DG b : Vh b V h b R and I DG : (V b h V Γ h ) (V h b A DG b (p h, q h ) = E T E h V h Γ ) R, dfind as follows K p h q h {K p h } q h E h \Γ h E h \Γ h {K q h } p h + E h \Γ h σ p h q h, I DG ((p h, p Γ h ), (q h, qh Γ )) = β Γ p h q h + α Γ ({p h } p Γ h )({q h} qh Γ ). Γ h Γ h
12 11 Th cofficint σ is dfind by σ = σ0 for any β 0 EI h EB h. Th pnalty paramtr σ 0 > 0 and th positiv numbr β 0 will b dfind latr on. Finally w dfin th linar functional L DG b : Vh b R as L DG b (q h ) = fq h. E T E h Rmark 1. Sinc w ar imposing homognous boundary conditions L DG b has th sam structur of th linar functional L b prviously dfind. In gnral, for g 0, L DG b contains som additional trms: L DG b (q h ) = E T h E fq h + Eh B ( K q h n + σ q h )g. Th DG discrtization of problm (13) rads as follows: Find (p h, p Γ h ) V h b V h Γ such that A h ( (ph, p Γ h ), (q h, q Γ h )) = L h (q h, q Γ h ) (q h, q Γ h ) V b h V Γ h, (17) whr A h : (Vh b V h Γ) (V h b V h Γ ) R is dfind as ( A h (ph, p Γ h ), (q h, qh Γ )) = A DG b (p h, q h ) + A Γ (p Γ h, qγ h ) + IDG ((p h, p Γ h ), (q h, qh Γ )), (18) and L h : Vh b V h Γ R is dfind as L h (q h, qh Γ ) = LDG b (q h ) + L Γ (qh Γ ). (19) Not that th discrt bilinar form A h has th sam structur as th bilinar form A b,γ prviously dfind, bing th sum of thr diffrnt componnts, ach rprsnting a spcific part of th problm. Bfor proving that formulation (17) is wll-posd, w stat (and prov) som auxiliary rsults, s Lmma 3.1 and 3.2 blow. W primarily introduc th following norm on V b h V Γ h (q h, q Γ h ) 2 E h = q h 2 DG + q Γ h 2 Γ + (q h, q Γ h ) 2 I, (20) whr q h 2 DG = K 1/2 h q h 2 L 2 (Ω) + E h \Γ h σ 1/2 q h 2 L 2 (), qh Γ 2 Γ = (KΓl τ Γ ) 1/2 h qh Γ 2 L 2 (Γ), (q h, qh Γ ) 2 I = β Γ q h 2 L 2 () + α Γ {q h } qh Γ 2 L 2 (). Γ h Γ h
13 12 It is asy to show that DG is a norm if σ 0 > 0 for all and that I is a norm if α Γ 0 (that is ξ > 1/2). Nxt w rcall two wll-known trac inqualitis (s for xampl [6]) for a polynomial q P k (E) ovr an dg E, whr E is a shap-rgular lmnt: q L 2 () h 1/2 E q L 2 (E), (21) q n L 2 () h 1/2 E q L 2 (E). (22) Ths inqualitis ar instrumntal to prov th subsqunt rsults. Lmma 3.1. Th bilinar form A DG b (, ) is corciv and continuous on Vh b V h b with rspct to th norm DG. Proof. Th proof is rathr standard in th analysis of Discontinuous Galrkin mthods (s for xampl [29]). W rcall that in ordr to prov corcivity, w us th trac inqualitis (21)-(22) and th bounddnss of th prmability tnsor K, s (2), to show that, for any q h V b h, E T h E K q h q h 2 E h \Γ h E T h {K q h } q h + E K q h q h + E h \Γ h σ E h \Γ h σ q h q h q h q h. Not that th inquality holds if th stabilization cofficints σ 0 ar chosn larg nough and β 0 (d 1) 1. Lmma 3.2. Th bilinar form A Γ (, ) is corciv and continuous on Vh Γ V h Γ with rspct to th norm Γ Proof. Sinc A Γ (q Γ h, qγ h ) = qγ h 2 Γ q Γ h V Γ h, A Γ (, ) is clarly corciv. Continuity follows dirctly from Cauchy-Schwarz inquality. Employing Lmma 3.1 and Lmma 3.2, w can asily prov th wllposdnss of th discrt problm (17). Proposition 3.3. Problm (17) is wll-posd.
14 13 Proof. In ordr to apply Lax-Milgram s lmma, w show that A h (, ) and L h ( ) ar continuous and corciv on Vh b V h Γ quippd with th norm (20). In ordr to prov corcivity w considr for (q h, qh Γ) V h b V h Γ, th trm I DG ((q h, qh Γ), (q h, qh Γ )). Clarly w hav I DG ((q h, q Γ h ), (q h, q Γ h )) = (q h, q Γ h ) 2 I. Morovr from Lmma 3.1 and Lmma 3.2 w know that A DG b (q h, q h ) q h 2 DG and A Γ(qh Γ, qγ h ) = qγ h 2 Γ, rspctivly. Thrfor ( A h (qh, qh Γ ), (q h, qh Γ )) (q h, qh Γ ) 2 E h (q h, qh Γ ) V h b V h Γ. (23) Nxt w prov continuity. Lt (q h, qh Γ), (w h, wh Γ) V h b V h Γ. Thn, from Lmma 3.1, w know that A DG b (q h, w h ) q h DG w h DG (q h, q Γ h ) E h (w h, w Γ h ) E h. Similarly w hav from Lmma 3.2 that A Γ (q Γ h, wγ h ) qγ h Γ w Γ h Γ (q h, q Γ h ) E h (w h, w Γ h ) E h. Finally, from Cauchy-Schwarz inquality, w gt I DG ((q h, qh Γ ), (w h, wh Γ )) β Γ q h 2 L 2 () w h 2 L 2 () Γ h + α Γ In conclusion w hav provd that Γ h {q h } q Γ h 2 L 2 () {w h} w Γ h 2 L 2 () (q h, q Γ h ) E h (w h, w Γ h ) E h. A h ( (qh, q Γ h ), (w h, w Γ h )) (q h, q Γ h ) E h (w h, w Γ h ) E h. Th continuity of L h ( ) on Vh b V h Γ can b asily provd using Cauchy- Schwarz inquality, thanks to th rgularity assumptions on th forcing trm f. Rmark 2. Th choic of mploying a conforming finit lmnt approximation for th flow in th fractur has bn mad in ordr to kp th analysis of th numrical mthod as clar as possibl. W rmark that DG mthods could b mployd for th fractur problm, too. Howvr, th analysis in th lattr cas would rquir unncssary tchnical dtails without giving any dpr undrstanding of th problm.
15 14 4 Error stimats In ordr to prov rror stimats, w assum that th xact solution (p, p Γ ) blongs to th spac V b V Γ dfind in (12), and w assum that in th dfinition of Vi b, i = 1, 2, w hav s 2. W will show that th discrt solution (p h, p Γ h ) to problm (17) convrgs to th xact solution (p, p Γ ), driving an a priori stimat for th rror in th norm (20). To this aim w introduc th spac V b (h) Vh Γ = (V h b + V b ) Vh Γ ndowd with th norm (p, p Γ h ) 2 = p 2 DG + p Γ h 2 Γ + (p, p Γ h ) 2 I, (24) whr DG : V b (h) R is dfind as q DG = q 2 DG + h 2 E q 2 H 2 (E) E Th 1/2. (25) Using classical intrpolation stimats (s for xampl [7]), it can b shown that, for any p H k+1 (T h ), k 1, thr xists p I b V h b such that p p I b DG E T h h k E p H k+1 (E). (26) Analogously, for any p Γ H k+1 (Γ h ), k 1, thr xists p I Γ V h Γ such that p Γ p I Γ Γ h k p Γ H k+1 (). Γ h (27) Th nxt continuity rsult will b crucial for th proof of our a priori rror stimat. Lmma 4.1. Th bilinar form A DG b (, ) is continuous on V b (h) V b (h) quippd with th norm DG. Proof. Lt q, w V b (h), thn w hav A DG b (v, w) = K q w E T E h {K w} q + E h \Γ h = T 1 + T 2 + T 3 + T 4. E h \Γ h {K q} w E h \Γ h σ q w
16 15 W bound ach of th trms sparatly. Using Cauchy-Schwarz inquality, w hav T 1 E T h K 1/2 q L 2 (E) K 1/2 w L 2 (E) q DG w DG, T 4 E h \Γ h σ q L 2 () w L 2 () q DG w DG q DG w DG. In ordr to bound T 2 and T 3 w will mak us of th following wll-known trac inquality (s [6]), for a function v H s (E), s 2, and E: v n 2 L 2 () ( 1 v 2 H 1 (E) + v 2 H 2 (E) ). In fact combining th abov rsult with (2) it follows that, for vry q L 2 (), K w n q w n L 2 () q L 2 () and this lattr implis that T 3 = {K w} q E h \Γ h [ ] 1/2 w 2 H 1 (E) + h2 E w 2 1/2 H 2 (E) h E q L 2 (), w 2 H 1 (E) + h2 E w 2 H 2 (E) E Th 1/2 w DG q DG w DG q DG. E h \Γ h σ q 2 Similarly T 2 q DG w DG w DG q DG. Summing all th contributions concluds th proof. 1/2 W now hav all th ingrdints to prov th following rror stimat: Thorm 4.2. Lt (p, p Γ ) b th solution of problm (13) and lt (p h, p Γ h ) V h b V h Γ b its approximation obtaind with th mthod (17). Thn, if (p, p Γ ) H k+1 (T h ) H k+1 (Γ h ), k 1, it holds (p, p Γ ) (p h, p Γ h ) E h h k ( p H k+1 (T h ) + p Γ H k+1 (Γ)). (28)
17 16 Proof. First w obsrv that th bilinar form A h (, ) is continuous on (V b (h) Vh Γ) (V b (h) Vh Γ ) with rspct to th norm (24). Rcalling that A h ( (q, q Γ h ), (w, w Γ h )) = A DG b (q, w) + A Γ (q Γ h, wγ h ) + IDG ((q, q Γ h ), (w, wγ h )), w considr ach of th trms sparatly. Lmma 4.1 implis A DG b (q, w) q DG w DG (q, q Γ h ) (w, wγ h ). Using Lmma 3.2 w also know that for all (q, q Γ h ), (w, wγ h ) V b (h) V Γ h A Γ (q Γ h, wγ h ) (q, qγ h ) (w, wγ h ). Finally w considr th trm I((q, qh Γ), (w, wγ h )). W hav, using Cauchy- Schwarz inquality, that I((q, q Γ h ), (w, wγ h )) (q, qγ h ) I (w, w Γ h ) I (q, q Γ h ) (w, wγ h ). Combining all th prvious bounds w obtain that A h ( (q, q Γ h ), (w, w Γ h )) (q, q Γ h ) (w, wγ h ). Nxt, lt p I b and pi Γ b th two intrpolant functions of stimats (26) and (27), rspctivly. Thus, w hav (p p I b, p Γ p I Γ) 2 I β Γ p p I b 2 L 2 () Γ h + α Γ {p p I b } 2 L 2 () + α Γ p Γ p I Γ 2 L 2 (). Γ h Γ h Using wll known intrpolation rror stimats, w hav β Γ p p I b 2 L 2 () β Γ p p I b 2 L 2 ( E) Γ h Similarly, α Γ Γ h {p p I b } 2 L 2 () E Γ E Γ p p I b 2 L 2 ( E) Morovr, from (27) w know that α Γ Γ h p Γ p I Γ 2 L 2 () Γ h h 2k E Γ E Γ p Γ 2 H k+1 (). h 2(k+1/2) E p 2 H k+1 (E). h 2(k+1/2) E p 2 H k+1 (E).
18 17 If w assum that for any E T h and for any E it holds h h E, combining all prvious bounds with (26) and (27), w obtain that ( ) (p, p Γ ) (p I b, pi Γ) h k p H k+1 (T h ) + p Γ H k+1 (Γ). In conclusion, sinc w hav alrady provd that A h (, ) is corciv on th discrt spac Vh b V h Γ (s (23)), and sinc Galrkin orthogonality trivially holds, th thsis follows (s [7]). 5 Numrical rsults In this sction w prsnt som two-dimnsional numrical xprimnts to confirm th validity of th a priori rror stimats that w hav drivd for our mthod. Th numrical rsults hav bn obtaind in Matlab R. Throughout this sction w st th fractur thicknss (apparing in th coupling conditions (8)-(9)) qual to l Γ = = η Γ and KΓ τ = 1. Morovr in our tsts w tak th dgr k of polynomials in th bulk and fractur finit dimnsional spacs qual to 1. W hav also implmntd a vrsion of our mthod which is abl to dal with mshs mad of gnral polygonal lmnts. Th analysis of th aformntiond mthod, which is basd on th framwork dvlopd in [17, 4], will b th objct of a futur publication. For th gnration of polygonal mshs conforming to th fracturs w hav usd a modification of th Matlab R cod calld PolyMshr implmntd by G.C. Paulino and collaborators [30]. In ach of our xampls w also rport th rsults of th numrical xprimnts prformd on polygonal mshs. 5.1 Exampl 1 In this first tst cas w tak Ω = (0, 1) 2, and choos as xact solutions in th bulk and in th fractur Γ = {(x, y) Ω : x + y = 1} p = { x+y in Ω 1, x+y + 4η Γ 2 in Ω 2, p Γ = + 2η Γ 2. It is asy to prov that p and p Γ satisfy th coupling conditions (8)-(9) with ξ = 1 and K = I. Not that w nd to choos f = 0 on Γ sinc th solution is constant and K p = 0.
19 18 K1/2 (p ph ) L2 (Ω) p ph L2 (Ω) 10 1 Γ q Γ qh L2 (Γ) h h Figur 2: Exampl 1: Computd rrors as a function of th invrs of th msh siz (loglog scal). In Figur 2 w plot th computd rrors K1/2 (p ph ) L2 (Ω) and q Γ qhγ L2 (Γ) as a function of th msh siz (loglog scal). In both cass th numrical rsults ar in agrmnt with th thortical stimats. In Figur 2 w also rport th bhaviour of th rror p ph L2 (Ω). On ordr of convrgnc is clarly gaind. (a) msh 1 (b) msh 2 (c) msh 3 Figur 3: Exampl 1: Thr rfinmnts of th polygonal msh grid conforming to th fractur. Finally, in Figur 4 w rport th computd rrors for th mthod implmntd on polygonal msh grids. Th sam convrgnc rats as in th cas of triangular mshs ar achivd. In Figur 3 w also rport thr rfinmnts of a polygonal msh conforming to th fractur.
20 K 1/2 (p p h ) L 2 (Ω) p p h L 2 (Ω) q Γ qh Γ L 2 (Γ) h h Figur 4: Exampl 1: Computd rrors as a function of th invrs of th msh siz (loglog scal) on polygonal grids. 5.2 Exampl 2 Hr w considr again Ω = (0, 1) 2 and Γ = {(x, y) Ω : x+y = 1}. W tak as xact solutions in th bulk and in th fractur th following functions p = { x+y in Ω 1, x+y 2 + ( η Γ 2 ) in Ω 2, p Γ = (1 + 2η Γ ). W choos ξ = 1 and tak again K = I. In this cas w st th sourc trm as 2 x+y in Ω 1, f = 2 in Γ, x+y in Ω 2. Notic that on th fractur th sourc trm satisfis f Γ = τ (K τ Γl Γ τ p Γ ) + K p, and, sinc p Γ is constant, it must b f Γ = K p. Figur 5 shows th computd rrors K 1/2 (p p h ) L 2 (Ω) for th bulk problm and th computd rrors q Γ q Γ h L 2 (Γ) in th fractur. Again th validity of th rror stimats is confirmd. Morovr from Figur 5 on can clarly s that also in this tst cas on ordr of convrgnc is gaind if w comput th rror p p h L 2 (Ω) for th bulk problm.
21 20 1/ K (p ph ) L 2 (Ω) p p h L 2 (Ω) 10 2 q Γ q Γ h L 2 (Γ) h h Figur 5: Exampl 2: Computd rrors as a function of th invrs of th msh siz (loglog scal). 1/2 K (p ph ) L 2 (Ω) 10 1 p p h L 2 (Ω) q Γ qh Γ L 2 (Γ) h h Figur 6: Exampl 2: Computd rrors as a function of th invrs of th msh siz (loglog scal) on polygonal grids. In Figur 6 w plot th computd rrors on polygonal mshs, which ar in agrmnt with th sam convrgnc rats as thos xpctd for simplicial mshs.
22 Exampl 3 In this last xampl w considr th circular fractur Γ = {(x, y) Ω : x2 + y 2 = R}, with R = 0.7 includd in th domain Ω = (0, 1)2. W choos th xact solutions in th bulk and in th fractur as follows ( 2 2 x +y in Ω1, 7 ηγ 2 pγ = 1 + p = x2r+y2, 1 3 4R + R ηγ + 2 in Ω2, 2R2 so that thy satisfy coupling conditions sourc trm is chosn as 4 R2 f = R1 2 R2 (a) msh 1 (8)-(9) with ξ = (b) msh and K = I. Th in Ω1, in Γ, in Ω2. (c) msh 3 Figur 7: Exampl 3: Thr rfinmnts of th msh grid. Figur 7 rports thr rfinmnts of th triangular msh grid. On can s that th fractur is approximatd by a polygonal lin. Analogously, Figur 8 shows thr rfinmnts of th polygonal msh grid mployd in our xprimnts. In Figur 9 w hav plottd th computd rrors on triangular mshs for th problms in th bulk and in th fractur, which clarly validat th thortical stimats. Figur 10 shows th computd rrors for th mthod implmntd on polygonal msh grids.
23 22 (a) msh 1 (b) msh 2 (c) msh 3 Figur 8: Exampl 3: Thr rfinmnts of th polygonal msh grid. K1/2 (p ph ) L2 (Ω) p ph L2 (Ω) 10 1 Γ L2 (Γ) q Γ qh h h Figur 9: Exampl 3: Computd rrors as a function of th invrs of th msh siz (loglog scal). 6 Conclusions In this work w proposd a numrical approximation of Darcy s flow in a fracturd porous mdium basd on DG mthods. In particular w combind a DG approximation for th flow in th porous matrix with a conforming finit lmnt schm for th fractur flow. Th mthod was formulatd in th cas of a non-immrsd singl fractur, taking as a rfrnc th modl proposd in [27], rwrittn in primal form. W analysd th proprtis of our discrt mthod and provd its wll-posdnss and convrgnc. In particular w obtaind optimal a priori rror stimats for both th problm in th bulk and in th fractur. Finally, w carrid out two dimnsional numrical
24 23 1/2 K (p ph ) L 2 (Ω) 10 1 p p h L 2 (Ω) q Γ qh Γ L 2 (Γ) h 10 2 h Figur 10: Exampl 3: Computd rrors as a function of invrs of th msh siz (loglog scal) on polygonal grids. xprimnts assssing th validity of th thortical stimats. Sinc th gomtric conformity of th grid to th fractur can lad to low-quality lmnts, a natural dvlopmnt of th prsnt work will b th xtnsion to mor gnral mshs of polygonal/polyhdral lmnts, including th cas of dgnrating dgs/facs. This will b th objct of a futur publication mploying th framwork of [17, 4]. Rfrncs [1] C. Alboin, J. Jaffré, J. E. Robrts, and C. Srrs. Modling fracturs as intrfacs for flow and transport in porous mdia. In Fluid flow and transport in porous mdia: mathmatical and numrical tratmnt (South Hadly, MA, 2001), volum 295 of Contmp. Math., pags Amr. Math. Soc., Providnc, RI, [2] C. Alboin, J. Jaffré, J. E. Robrts, X. Wang, and C. Srrs. Domain dcomposition for som transmission problms in flow in porous mdia. In Numrical tratmnt of multiphas flows in porous mdia (Bijing, 1999), volum 552 of Lctur Nots in Phys., pags Springr, Brlin, 2000.
25 24 [3] P. Angot, F. Boyr, and F. Hubrt. Asymptotic and numrical modlling of flows in fracturd porous mdia. M2AN Math. Modl. Numr. Anal., 43(2): , [4] P. F. Antonitti, A. Cangiani, J. Collis, Z. Dong, E. H. Gorgoulis, S. Giani, and P. Houston. Rviw of discontinuous galrkin finit lmnt mthods for partial diffrntial quations on complicatd domains [5] P. F. Antonitti, L. Formaggia, A. Scotti, M. Vrani, and N. Vrzotti. Mimtic finit diffrnc approximation of flows in fracturd porous mdia. Accptd for publication on M2AN Math. Modl. Numr. Anal., [6] D. N. Arnold. An intrior pnalty finit lmnt mthod with discontinuous lmnts. SIAM J. Numr. Anal., 19(4): , [7] 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(5): , 2001/02. [8] G. A. Bakr. Finit lmnt mthods for lliptic quations using nonconforming lmnts. Math. Comp., 31(137):45 59, [9] 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(2): , [10] J. Bar, C. F. Tsang, and G. d. Marsily. Flow and contaminant transport in fracturd rocks [11] M. F. Bndtto, S. Brron, and S. Scialò. A globally conforming mthod for solving flow in discrt fractur ntworks using th virtual lmnt mthod. Finit Elmnts in Analysis and Dsign, 109:23 36, [12] S. Brron, S. Piraccini, and S. Scialò. A PDE-constraind optimization formulation for discrt fractur ntwork flows. SIAM J. Sci. Comput., 35(2):B487 B510, [13] S. Brron, S. Piraccini, and S. Scialò. An optimization approach for larg scal simulations of discrt fractur ntwork flows. J. Comput. Phys., 256: , 2014.
26 25 [14] S. Brron, S. Piraccini, S. Scialò, and F. Vicini. A paralll solvr for larg scal DFN flow simulations. SIAM J. Sci. Comput., 37(3):C285 C306, [15] S. Brron, S. Piraccini, S. Scialò, and F. Vicini. A paralll solvr for larg scal DFN flow simulations. SIAM J. Sci. Comput., 37(3):C285 C306, [16] F. Brzzi, B. Cockburn, L. D. Marini, and E. Süli. Stabilization mchanisms in discontinuous Galrkin finit lmnt mthods. Comput. Mthods Appl. Mch. Engrg., 195(25-28): , [17] A. Cangiani, E. H. Gorgoulis, and P. Houston. hp-vrsion discontinuous Galrkin mthods on polygonal and polyhdral mshs. Math. Modls Mthods Appl. Sci., 24(10): , [18] P. Castillo, B. Cockburn, I. Prugia, and D. Schötzau. An a priori rror analysis of th local discontinuous Galrkin mthod for lliptic problms. SIAM J. Numr. Anal., 38(5): (lctronic), [19] B. Cockburn and C. Dawson. Som xtnsions of th local discontinuous Galrkin mthod for convction-diffusion quations in multidimnsions. In Th mathmatics of finit lmnts and applications, X, MAFELAP 1999 (Uxbridg), pags Elsvir, Oxford, [20] C. D Anglo and A. Scotti. A mixd finit lmnt mthod for darcy flow in fracturd porous mdia with non-matching grids. ESAIM: Mathmatical Modlling and Numrical Analysis, 46(02): , [21] D. A. Di Pitro and A. Ern. Mathmatical aspcts of discontinuous Galrkin mthods, volum 69. Springr Scinc & Businss Mdia, [22] J. Douglas, Jr. and T. Dupont. Intrior pnalty procdurs for lliptic and parabolic Galrkin mthods. In Computing mthods in applid scincs (Scond Intrnat. Sympos., Vrsaills, 1975), pags Lctur Nots in Phys., Vol. 58. Springr, Brlin, [23] N. Frih, J. E. Robrts, and A. Saada. Modling fracturs as intrfacs: a modl for Forchhimr fracturs. Comput. Gosci., 12(1):91 104, [24] A. Fumagalli and A. Scotti. A numrical mthod for two-phas flow in fracturd porous mdia with non-matching grids. Advancs in Watr Rsourcs, 62, Part C: , 2013.
27 26 [25] J. S. Hsthavn and T. Warburton. Nodal discontinuous Galrkin mthods: algorithms, analysis, and applications. Springr Scinc & Businss Mdia, [26] J. Jaffrré, M. Mnjja, and J. Robrts. A discrt fractur modl for two-phas flow with matrix-fractur intraction. Procdia Computr Scinc, 4: , [27] V. Martin, J. Jaffré, and J. E. Robrts. Modling fracturs and barrirs as intrfacs for flow in porous mdia. SIAM J. Sci. Comput., 26(5): (lctronic), [28] W. Rd and T. Hill. Triangular msh mthods for th nutron transport quation. Los Alamos Rport LA-UR , [29] B. Rivièr. Discontinuous Galrkin mthods for solving lliptic and parabolic quations, volum 35 of Frontirs in Applid Mathmatics. Socity for Industrial and Applid Mathmatics (SIAM), Philadlphia, PA, Thory and implmntation. [30] C. Talischi, G. H. Paulino, A. Prira, and I. F. Mnzs. Polymshr: a gnral-purpos msh gnrator for polygonal lmnts writtn in matlab. Structural and Multidisciplinary Optimization, 45(3): , [31] M. F. Whlr. An lliptic collocation-finit lmnt mthod with intrior pnaltis. SIAM J. Numr. Anal., 15(1): , 1978.
28 MOX Tchnical Rports, last issus Dipartimnto di Matmatica Politcnico di Milano, Via Bonardi Milano (Italy) 21/2016 Ambrosi, D.; Zanzottra, A. Mchanics and polarity in cll motility 19/2016 Gurciotti, B.; Vrgara, C. Computational comparison btwn Nwtonian and non-nwtonian blood rhologis in stnotic vssls 20/2016 Wilhlm, M.; Sangalli, L.M. Gnralizd Spatial Rgrssion with Diffrntial Rgularization Gurciotti, B.; Vrgara, C. Computational comparison btwn Nwtonian and non-nwtonian blood rhologis in stnotic vssls 18/2016 Frroni, A.; Antonitti, P.F.; Mazziri, I.; Quartroni, A. Disprsion-dissipation analysis of 3D continuous and discontinuous spctral lmnt mthods for th lastodynamics quation 17/2016 Pnati, M.; Miglio, E. A nw mixd mthod for th Stoks quations basd on strss-vlocity-vorticity formulation 15/2016 Iva, F.; Paganoni, A.M. A taxonomy of outlir dtction mthods for robust classification in multivariat functional data 16/2016 Agosti, A.; Antonitti, P.F.; Ciarltta, P.; Grasslli, M.; Vrani, M. A Cahn-Hilliard typ quation with dgnrat mobility and singl-wll potntial. Part I: convrgnc analysis of a continuous Galrkin finit lmnt discrtization. 14/2016 Bonomi, D.; Manzoni, A.; Quartroni, A. A matrix discrt mpirical intrpolation mthod for th fficint modl rduction of paramtrizd nonlinar PDEs: application to nonlinar lasticity problms 13/2016 Gurciotti, B; Vrgara, C; Ippolito, S; Quartroni, A; Antona, C; Scrofani, R. Computational study of th risk of rstnosis in coronary bypasss
Discontinuous 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 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 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 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 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 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 informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 45, No. 3, pp. 1269 1286 c 2007 Socity for Industrial and Applid Mathmatics NEW FINITE ELEMENT METHODS IN COMPUTATIONAL FLUID DYNAMICS BY H (DIV) ELEMENTS JUNPING WANG AND XIU
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationDivision of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.
Division of Mchanics Lund Univrsity MULTIBODY DYNMICS Examination 7033 Nam (writ in block lttrs):. Id.-numbr: Writtn xamination with fiv tasks. Plas chck that all tasks ar includd. clan copy of th solutions
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 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 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 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 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 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 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 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 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 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 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 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 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 information22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.
Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M
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 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 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 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 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 informationBackground: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.
Chaptr 7 Th Hydrogn Atom Background: W hav discussd th PIB HO and th nrgy of th RR modl. In this chaptr th H-atom and atomic orbitals. * A singl particl moving undr a cntral forc adoptd from Scott Kirby
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 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 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 information16. Electromagnetics and vector elements (draft, under construction)
16. Elctromagntics (draft)... 1 16.1 Introduction... 1 16.2 Paramtric coordinats... 2 16.3 Edg Basd (Vctor) Finit Elmnts... 4 16.4 Whitny vctor lmnts... 5 16.5 Wak Form... 8 16.6 Vctor lmnt matrics...
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 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 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 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 informationPrinciples of Humidity Dalton s law
Principls of Humidity Dalton s law Air is a mixtur of diffrnt gass. Th main gas componnts ar: Gas componnt volum [%] wight [%] Nitrogn N 2 78,03 75,47 Oxygn O 2 20,99 23,20 Argon Ar 0,93 1,28 Carbon dioxid
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 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 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 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 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 informationCOHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.
MTH rviw part b Lucian Mitroiu Th LOG and EXP functions Th ponntial function p : R, dfind as Proprtis: lim > lim p Eponntial function Y 8 6 - -8-6 - - X Th natural logarithm function ln in US- log: function
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 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 informationConstruction of Mimetic Numerical Methods
Construction of Mimtic Numrical Mthods Blair Prot Thortical and Computational Fluid Dynamics Laboratory Dltars July 17, 013 Numrical Mthods Th Foundation on which CFD rsts. Rvolution Math: Accuracy Stability
More informationEXST Regression Techniques Page 1
EXST704 - Rgrssion Tchniqus Pag 1 Masurmnt rrors in X W hav assumd that all variation is in Y. Masurmnt rror in this variabl will not ffct th rsults, as long as thy ar uncorrlatd and unbiasd, sinc thy
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 informationMA 262, Spring 2018, Final exam Version 01 (Green)
MA 262, Spring 218, Final xam Vrsion 1 (Grn) INSTRUCTIONS 1. Switch off your phon upon ntring th xam room. 2. Do not opn th xam booklt until you ar instructd to do so. 3. Bfor you opn th booklt, fill in
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 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 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 informationSliding Mode Flow Rate Observer Design
Sliding Mod Flow Rat Obsrvr Dsign Song Liu and Bin Yao School of Mchanical Enginring, Purdu Univrsity, Wst Lafaytt, IN797, USA liu(byao)@purdudu Abstract Dynamic flow rat information is ndd in a lot of
More informationPropositional Logic. Combinatorial Problem Solving (CPS) Albert Oliveras Enric Rodríguez-Carbonell. May 17, 2018
Propositional Logic Combinatorial Problm Solving (CPS) Albrt Olivras Enric Rodríguz-Carbonll May 17, 2018 Ovrviw of th sssion Dfinition of Propositional Logic Gnral Concpts in Logic Rduction to SAT CNFs
More informationMapping properties of the elliptic maximal function
Rv. Mat. Ibroamricana 19 (2003), 221 234 Mapping proprtis of th lliptic maximal function M. Burak Erdoğan Abstract W prov that th lliptic maximal function maps th Sobolv spac W 4,η (R 2 )intol 4 (R 2 )
More informationMCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)
MCB37: Physical Biology of th Cll Spring 207 Homwork 6: Ligand binding and th MWC modl of allostry (Du 3/23/7) Hrnan G. Garcia March 2, 207 Simpl rprssion In class, w drivd a mathmatical modl of how simpl
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 informationBINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES. 1. Statement of results
BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES DONALD M. DAVIS Abstract. If p is a prim and n a positiv intgr, lt ν p (n dnot th xponnt of p in n, and u p (n n/p νp(n th unit part of n. If α
More informationCHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle
CHPTER 1 Introductory Concpts Elmnts of Vctor nalysis Nwton s Laws Units Th basis of Nwtonian Mchanics D lmbrt s Principl 1 Scinc of Mchanics: It is concrnd with th motion of matrial bodis. odis hav diffrnt
More informationGEOMETRICAL PHENOMENA IN THE PHYSICS OF SUBATOMIC PARTICLES. Eduard N. Klenov* Rostov-on-Don, Russia
GEOMETRICAL PHENOMENA IN THE PHYSICS OF SUBATOMIC PARTICLES Eduard N. Klnov* Rostov-on-Don, Russia Th articl considrs phnomnal gomtry figurs bing th carrirs of valu spctra for th pairs of th rmaining additiv
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 informationData Assimilation 1. Alan O Neill National Centre for Earth Observation UK
Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal
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 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 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 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 informationA Sub-Optimal Log-Domain Decoding Algorithm for Non-Binary LDPC Codes
Procdings of th 9th WSEAS Intrnational Confrnc on APPLICATIONS of COMPUTER ENGINEERING A Sub-Optimal Log-Domain Dcoding Algorithm for Non-Binary LDPC Cods CHIRAG DADLANI and RANJAN BOSE Dpartmnt of Elctrical
More informationHigh Energy Physics. Lecture 5 The Passage of Particles through Matter
High Enrgy Physics Lctur 5 Th Passag of Particls through Mattr 1 Introduction In prvious lcturs w hav sn xampls of tracks lft by chargd particls in passing through mattr. Such tracks provid som of th most
More informationCOMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.
C:\Dallas\0_Courss\03A_OpSci_67\0 Cgh_Book\0_athmaticalPrliminaris\0_0 Combath.doc of 8 COPUTER GENERATED HOLOGRAS Optical Scincs 67 W.J. Dallas (onday, April 04, 005, 8:35 A) PART I: CHAPTER TWO COB ATH
More informationChemical Physics II. More Stat. Thermo Kinetics Protein Folding...
Chmical Physics II Mor Stat. Thrmo Kintics Protin Folding... http://www.nmc.ctc.com/imags/projct/proj15thumb.jpg http://nuclarwaponarchiv.org/usa/tsts/ukgrabl2.jpg http://www.photolib.noaa.gov/corps/imags/big/corp1417.jpg
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 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 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 informationLarge Scale Topology Optimization Using Preconditioned Krylov Subspace Recycling and Continuous Approximation of Material Distribution
Larg Scal Topology Optimization Using Prconditiond Krylov Subspac Rcycling and Continuous Approximation of Matrial Distribution Eric d Sturlr*, Chau L**, Shun Wang***, Glaucio Paulino** * Dpartmnt of Mathmatics,
More informationCE 530 Molecular Simulation
CE 53 Molcular Simulation Lctur 8 Fr-nrgy calculations David A. Kofk Dpartmnt of Chmical Enginring SUNY Buffalo kofk@ng.buffalo.du 2 Fr-Enrgy Calculations Uss of fr nrgy Phas quilibria Raction quilibria
More informationDirect Approach for Discrete Systems One-Dimensional Elements
CONTINUUM & FINITE ELEMENT METHOD Dirct Approach or Discrt Systms On-Dimnsional Elmnts Pro. Song Jin Par Mchanical Enginring, POSTECH Dirct Approach or Discrt Systms Dirct approach has th ollowing aturs:
More informationGeneral Notes About 2007 AP Physics Scoring Guidelines
AP PHYSICS C: ELECTRICITY AND MAGNETISM 2007 SCORING GUIDELINES Gnral Nots About 2007 AP Physics Scoring Guidlins 1. Th solutions contain th most common mthod of solving th fr-rspons qustions and th allocation
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 information