Convergence results for the vector penalty-projection and two-step artificial compressibility methods
|
|
- Silvia Sanders
- 5 years ago
- Views:
Transcription
1 Convrgnc rsults for th vctor pnalty-projction and two-stp artificial comprssibility mthods Philipp Angot, Pirr Fabri To cit this vrsion: Philipp Angot, Pirr Fabri. Convrgnc rsults for th vctor pnalty-projction and twostp artificial comprssibility mthods. Discrt and Continuous Dynamical Systms - Sris B, Amrican Institut of Mathmatical Scincs,, 7 5), pp , communicatd by Rogr Tmam. <.3934/dcdsb >. <hal-6533v> HAL Id: hal Submittd on 9 Fb HAL is a multi-disciplinary opn accss archiv for th dposit and dissmination of scintific rsarch documnts, whthr thy ar publishd or not. Th documnts may com from taching and rsarch institutions in Franc or abroad, or from public or privat rsarch cntrs. L archiv ouvrt pluridisciplinair HAL, st dstiné au dépôt t à la diffusion d documnts scintifiqus d nivau rchrch, publiés ou non, émanant ds établissmnts d nsignmnt t d rchrch français ou étrangrs, ds laboratoirs publics ou privés.
2 Manuscript submittd to AIMS Journals Volum X, Numbr X, XX X Wbsit: pp. X XX CONVERGENCE RESULTS FOR THE VECTOR PENALTY-PROJECTION AND TWO-STEP ARTIFICIAL COMPRESSIBILITY METHODS Philipp Angot Aix-Marsill Univrsité Laboratoir d Analys, Topologi, Probabilités - CNRS UMR7353 Cntr d Mathématiqus t Informatiqu 3453 Marsill cdx 3 - Franc Pirr Fabri Univrsité d Bordaux & IPB Institut Mathématiqus d Bordaux - CNRS UMR55 ENSEIRB-MATMECA, Talnc - Franc Communicatd by Rogr Tmam) Abstract. In this papr, w propos and analyz a nw artificial comprssibility splitting mthod which is issud from th rcnt vctor pnalty-projction mthod for th numrical solution of unstady incomprssibl viscous flows introducd in [], [] and [3]. This mthod may b viwd as an hybrid two-stp prdiction-corrction mthod combining an artificial comprssibility mthod and an augmntd Lagrangian mthod without innr itration. Th prturbd systm can b viwd as a nw approximation to th incomprssibl Navir- Stoks quations. In th main rsult, w stablish th convrgnc of solutions to th wak solutions of th Navir-Stoks quations whn th pnalty paramtr tnds to zro.. Introduction and stting of th problm. Th artificial comprssibility mthod was introducd by Chorin [6] and Tmam [7] for th solution of th unstady incomprssibl Stoks or Navir-Stoks quations; s also [] for th thortical analysis. Thn, som othr numrical schms to fficintly comput th solutions of Navir-Stoks problms can b viwd as discrtizations of prturbd systms of th typ of pnalization [4] or psudo-comprssibility. This is th cas with th famous projction mthods from Chorin [7] and Tmam [8, 9] and thir variants [], s.g. [5]. Hr, w prsnt a nw approximation mthod for th Navir-Stoks quations modling incomprssibl viscous flows in a boundd rgular opn st ndowd with Dirichlt boundary conditions on Γ = Lipschitz-continuous). With a givn sourc trm f, th Navir-Stoks systm rads: Mathmatics Subjct Classification. Primary: 35Q3, 76D5, 76N, 35A35; Scondary: 65M, 65N. Ky words and phrass. Artificial comprssibility, Navir-Stoks quations, vctor pnaltyprojction, psudo-comprssibility, pnalty mthod.
3 PHILIPPE ANGOT AND PIERRE FABRIE v t +v )v v + p = f, div v =, v) = v,v Γ =, whr dnots th Rynolds numbr. Accordingtothidntity, ϕ = curl curl ϕ div ϕ, wconsidrthfollowing approximat mthod to obtain a solution of th abov Navir-Stoks systm, with th paramtrs r,γ > and,ε > ṽ ε t +v ε )ṽ ε + div v ε)ṽ ε + curl curl ṽ ε div ṽ ε r div ṽ ε + p ε = f v ε t +v ε ) v ε + div v ε) v ε + curl curl v ε div v ε r div v ε div v ε +div ṽ ε ) = ε v ε = ṽ ε + v ε γ p ε t +γp ε + ε div v ε +rdiv ṽ ε =. W associat to th prvious systm th following boundary conditions and initial data, ṽ ε ) = v, v ε ) =, p ε ) = p, ṽ ε ν Γ =, v ε ν Γ =, curl ṽ ε ) ν Γ =, curl v ε ) ν Γ =, whr ν dnots th outward unit normal vctor on Γ. To vanish, at th limit procss, th two tangntial componnt of th vlocity filds, ṽ ε ν and v ε ν, w us a pnalization mthod which will b dtaild blow. This mthod is clos to th artificial comprssibility mthod of Chorin [6] and Tmam[7], but prsnts on important diffrnc. It is a two-stp splitting mthod. Thfirstquationofthprvioussystmgivsaprdictdvlocityṽ ε andthscond on is th approximat projction of ṽ ε on th fr-divrgnc vctor filds. This quation may b sn as an approximat mthod to solv th wll-posd problm s appndix A) : div v ε = div ṽ ε, curl v ε =, v ε.ν Γ =. Rmark. This approximat mthod is issud from th Vctor Pnalty-Projction VPP r,ε ) mthods for th numrical solution of unstady incomprssibl viscous flows introducd in [] and [3]. A fast vrsion of ths mthods, th so-calld VPP ε ) mthod, is rcntly proposd also for th numrical solution of th nonhomognous Navir-Stoks quations in [, 4]. It is shown to b vry fficint to comput multiphas flows, i.. fast, chap, and robust whatvr th dnsity, viscosity or prmability jumps.
4 THE VECTOR PENALTY-PROJECTION METHOD 3 Evn for r =, th rsulting mthod which corrsponds to a two-stp psudocomprssibility mthod, is diffrnt from th original artificial comprssibility mthod of Chorin [6] and Tmam [7, ]. Th nw important point is th pnalty trm ε div v ε + div ṽ ε ) that appars in th vlocity corrction stp which allows us a dirct stimat on th divrgnc of th vlocity. Morovr, this systm is quit asy to solv and prsnts good stability proprtis, s [,, 3]. Th vlocity v ε and th prssur p ε satisfy th quations: v ε t + v ε )v ε + div v ) ε)v ε + curl curl v ε div v ε r div v ε div v ε + p ε = f ε γ p ε t +γp ε + ε div v ε +rdiv ṽ ε =, v ε ) = v, p ε ) = p, v ε ν =, curl v ε ) ν Γ =. Γ Th vanishing of th tangntial componnt at th limit procss, is fullfilld by a pnalization mthod, which implis that this boundary condition is satisfid at th ordr ε for th approximat solution... Notations. Lt b a rgular boundd and connctd opn st of R d, for d = or 3. W not H s ) th classical Sobolv spac, and H s th associatd norm. Th norm of a function in L p ) is dnotd by L p, and if B is a Banach spac, w dnot by. L p,b th norm in L p ],T[;B). L p ) = L p )) d H div ) = {v L )) d,div v L )} H = {v L )) d,div v =,v ν Γ = } H ν) = {v H )) d,v ν Γ = } G = {v L )) d, q H ),v = q}.. Mathmatical rcalls. Proposition. Undr th prvious hypothsis, on has th following proprtis: L ) = H G, Kr curl ) = G. Morovr, thr xists on constant C > dpnding only on such that: u H = u L + u L C u L + div u L + curl u L ), u H ν ). ) Bsids, if w suppos that th opn st is simply-connctd, thr xist two constants λ and λ dpnding only on such that: and w hav: u L λ div u L + curl u L ), u H ν ), u L + u L λ div u L + curl u L ), u H ν ) Kr curl ) H = {}. Proof. All ths rsults may b found in [9] and [8].
5 4 PHILIPPE ANGOT AND PIERRE FABRIE For a Banach spac E w introduc th Nikolskii spac dfind for q < +, < σ < : N σ q ],T[;E) = { } f L q f +h) f ) Lq ],T h[;e) ],T[;E), sup <h<t h σ < +, ndowd with th following norm: f N σ q = ) q ) f q q L q ],T[;E) + sup <h<t h σ f +h) f ) L q ],T h[;e). Lt us rcall th following proprty s for xampl [5] pag 5, Proposition. Lt H b an Hilbrt spac and f a function givn in L ],T[;H) such that, for som < σ <, R τ σ F f)τ) H dτ C, whr f dnots th xtnsion by of th function f outsid [,T]. Thn f N σ ],T[;H) and w hav f N σ M σ +C), whr M σ is a constant dpnding only on σ. W now rcall th important compactnss thorm, s for xampl [6] Thorm.. Aubin-Lions-Simon Lt B, B, B thr Banach spacs with B B B with continuous imbdding. Suppos morovr that th injction of B in B is compact. Thn, for all q + and < σ <, th imbdding L q ],T[;B ) N σ q ],T[;B ) L q ],T[;B ) is compact.. Main rsult. W associat to th prvious approximat systm, th variational problm whr th tangntial componnts of th vlocitis ṽ ε and v ε ar pnalizd. This problm is studid in th nxt sction.
6 Find ṽ ε, v ε,p ε ) in THE VECTOR PENALTY-PROJECTION METHOD 5 L ],T[; L )) L ],T[; H ν))) L ],T[; L )) satisfying in D ],T[), ṽ ε v t ϕdω + ε )ṽ ε + ) div v ε)ṽ ε ϕdω + curl ṽ ε curl ϕdω + div ṽ ε div ϕdω +r div ṽ ε div ϕdω p ε div ϕdω + ṽ ε ν) ϕ ν)dσ = f ϕdx, ε Γ v ε v t ψdω + ε ) v ε + ) div v ε) v ε ψdω + curl v ε curl ψdω + div v ε div ψdω Γ +r div v ε div ψdω + ε + v ε ν) ψ ν)dσ =, ε div ṽ ε +div v ε ) div ψdω ) v ε = ṽ ε + v ε, p ε γ πdω +γ ε πdω + t p πdiv v ε dω +r πdiv ṽ ε dω =, ε ϕ,ψ,π) H ν)) L ), ṽ ε ) = v, v ε ) =, p ε ) = p. Thn th vlocity v ε and th prssur p ε satisfy in D ],T[), v ε v t ϕdω + ε )v ε + ) div v ε)v ε ϕdω + curl v ε curl ϕdω + div v ε div ϕdω +r div v ε div ϕdω + div v ε div ϕdω ε ε div ϕdω + p v ε ν) ϕ ν)dσ ε Γ = f ϕdω, ϕ H ν), v ε ) = v 3) Rmark. In ordr to stablish th strong convrgnc of th squnc v ε ) ε> whn ε, w us in Sction 4 th Lray s orthogonal dcomposition in th
7 6 PHILIPPE ANGOT AND PIERRE FABRIE boundd domain. Th curl-fr componnt vanishs with th pnalty trm introducd by our mthod, whras th divrgnc-fr componnt strongly convrgs thanks to an stimat of a fractional drivativ in tim, s []. Howvr, this rquirs to considr vlocity filds having only thir normal componnt which is zro on th boundary. Sinc at th limit procss, w aim at solving th Navir- Stoks problm with homognous Dirichlt boundary condition, w also pnaliz th tangntial part of th vlocity filds. W prov in sction 3 th following rsults. Lmma.. Lt us suppos that f blongs to L ],T[;L )). Thn, thr xists at last a solution to th systm ). This solution is uniqu in two spac dimnsion. For th dimnsion d 3, this solution satisfis th following nrgy inquality: d rε ṽε L dt +ε v ε L + v ε L +γε p ε ) L +γε p ε L + rε curl ṽ ε L R + ε curl v ε L R + curl v ε L R + div v ε L + rε div ṽ ε L + ε div v ε L +εr div ṽ ε L +εr div v ε L + ε div v ε L +r ṽ ε ν) L Γ) + v ε ν) L Γ) + ε v ε ν) L Γ) λ +rε) f L For th dimnsion d =, on has th following nrgy quality: d rε ṽε L dt + v ε L +γε p ε ) rε L + curl ṽ ε L R + curl v ε L +γε p ε R L + rε div ṽ ε L R + div v ε L R +εr div ṽ ε L + div v ε L εr + ṽ ε ν) L Γ + ε v ε ν) L Γ) = rε f ṽ ε dω + f v ε dω. This rsult is quit classical and w only giv th sktch of proof in th sction 3. In fact, w can prcis th prvious nrgy inquality if w suppos that th data f blongs to L R + ;L )). This shows th absolut stability of th approximat mthod. Thorm.. Suppos that th data f satisfis f L R +,L )), thn, thr xists a constant α indpndnt of th data, such that rε ṽε t) L +ε v εt) L + v εt) L +γε p εt) L ) αt +rε) v L +γε p L ) + λ α +rε) f L,L, t R +. Th main goal of this papr is to prov th following convrgnc thorm:
8 THE VECTOR PENALTY-PROJECTION METHOD 7 Thorm.3. For d 3, thr xists a subsqunc v εk,p εk ) k solution of 3) that convrgs to a solution v, p) to th systm of Navir-Stoks quations with homognous Dirichlt boundary conditions. For d =, th solution v,p) is uniqu, and for all squncs ε k, v εk,p εk ) εk convrgs to v,p). Morovr, for all squncs ε k, v εk ) k convrgs strongly to v in L,T;H ν)). W now giv an intrprtation of th prssur and prcis its convrgnc. Lt us dfin q ε = p ε +ε ε +r ) div v ε. Th scalar function q ε appars to b th ffctiv approximat prssur, and w hav Thorm.4. Th function q εk satisfis if d = 3, q εk convrgs wakly to p in H,T) ) ) 3 if d =, q εk convrgs strongly to p in H,T) ) ) Ths convrgnc rsults for both vlocity and prssur ar provd in Sction Enrgy stimats. W first stablish th following xistnc rsult. Proposition 3. For v ε, p givn in L ) L ), thr xists at last a solution of th systm ) satisfying for d = 3: ṽ ε L ],T[;L ) L ],T[;H ν)), ṽ ε t L4 3 ],T[;H ν )) ) v ε L ],T[;L ) L ],T[;H ν)), v ε t L4 3 ],T[;H ν )) ) p ε L ],T[;L )), p ε t L ],T[;L )) ṽ ε ) = v εo in H ν)), v ε ) = in H ν)), p ε ) = p in L ). If d =, th uniqu solution of ) satisfis th following rgularity rsults: ṽ ε L ],T[;L ) L ],T[;H ν)), ṽ ε t L ],T[;H ν)) ))+L ],T[;L 3 )) v ε L ],T[;L ) L ],T[;H ν)), v ε t L ],T[;H ν)) ))+L 4 3 ],T[;L 4 3 )) p ε L ],T[;L )), p ε t L ],T[;L )) ṽ ε ) = v εo in L ), v ε ) = in L ), p ε ) = p in L ). Rmark 3. In th thr-dimnsional cas, th qualitis ar valid in th trac sns. ṽ ε ) = v εo in H ν)), v ε ) = in H ν)),
9 8 PHILIPPE ANGOT AND PIERRE FABRIE Proof. For fixd paramtrs ε >, r and γ >, w build approximat solutions by a classical Galrkin procss. Lt us introduc th slf-adjoint oprator A = curl curl div dfind on th domain H ν) H )) d. Thn, for th approximation of th two filds of vlocity, w us as spcial basis th ignfunctions of this oprator associatd with th following boundary conditions: u ν Γ =, curl u) ν Γ =. For th prssur, on can us as spcial basis th ignfunctions of th slf-adjoint oprator A = with domain H ) associatd to th Numann boundary conditions. This approximat finit dimnsional systm is thn a classical ordinary diffrntial quation which has a uniqu solution. Nxt, to prform th limit w us th sam stratgy as for th classical Navir-Stoks quations i.. a priori stimats and compactnss rsults using an stimat on th tmporal drivativ, s for xampl [3], [],[5]. Now, w will focus our attntion on th stimats on th tim drivativ according to th dimnsion d. Lt us bgin with th thr-dimnsional cas. W hav to stimat th two nonlinar trms v ε w, wdiv v ε, with ithr w = ṽ ε or w = v ε and th prssur p ε. Suppos first that d = 3. By Sobolv mbdding, th two nonlinar trms of th form v ε w and wdiv v ε blong to L 4 3,T;H ν) ) sinc w hav, for xampl for all ϕ H ν) : v ε )w ϕdω + div v ε )w ϕdω C ) v L 3 w H ν ϕ L 6 + w L 3 v ε H ν ϕ L 6 ) C v ε L v ε H w ν H + w ν L w H v ε H ϕ ν ν H ν. Th bounds of th linar trms ar straightforward and w hav curl ṽ ε curl ϕdω + div ṽ ε div ϕdω +r div ṽ ε div ϕdω C ṽ ε H ν ϕ H ν p ε div ϕdω C p ε L ϕ H ν ) div ṽε +div v ε div ϕdω C ṽ ε H + v ) ν ε H ϕ H ν ν. and, with standard trac thorms v ε ν) ϕ ν)dω C v ε H ν ϕ H ν. Thus it follows, from quation ) that Γ ṽ ε t L4 3 ],T[;H ν )) ), p ε t L ],T[;L )) v ε t L4 3 ],T[;H ν )) )
10 THE VECTOR PENALTY-PROJECTION METHOD 9 Ths stimats show that th vlocitis ṽ ε, v ε ) ar qual almost vrywhr to continuous functions with valus in H ν)). Bsids, th prssur p ε is qual almost vrywhr to a continuous function with valu in L ). For th two-dimnsional cas, th situation is quit diffrnt. W obsrv first that th vlocity filds v ε and ṽ ε blong to L 4 ],T[;L 4 )), so that: vε )ṽ ε + div vε )ṽε L 4 3 ],T[;L 4 3 )), vε ) v ε + div vε ) vε L 4 3 ],T[;L 4 3 )). So it follows from th quation ) that ṽ ε t L ],T[;H ν)) )+L ],T[;L 3 )) v ε t L ],T[;H ν)) )+L ],T[;L 3 )) p ε t L ],T[;L )) W now obsrv that th two vlocity filds ṽ ε and v ε blong to which is th dual spac of L 4 ],T[;L 4 )) L ],T[;H ν)) ) L ],T[;H ν)) )+L 4 3 ],T[;L 4 3 )). Thus th functions ṽ ε, v ε ) ar qual almost vrywhr to continuous functions with valus in L ). This nds th proof of proposition Stability. In th cas of thr-dimnsional vctor spacs, w do not hav an quality for th consrvation of th nrgy, w hav only an inquality. Nvrthlss for two-dimnsional vctor spacs, th wak solutions satisfy th nrgy quality. Proof. Through classical computations on obtains, with quations ) and 3): d dt ṽ ε L + ṽ R curl ε L + div ṽ ε L +r div ṽ ε R L + ε ṽ ε ν) L Γ) 4) p ε div ṽ ε dω = f ṽ ε dω, d dt v ε L + v R curl ε L + div v ε L +r div v ε R L + ε v ε ν) L Γ) + 5) div v ε div v ε dω =, ε d dt v ε L + v R curl ε L + div v ε L + div v R ε L εr +r div v ε L + 6) ε v ε ν) L Γ) p ε div v ε dω = f v ε dω, γε d dt p ε L +γε p ε L + p ε div v ε dω +rε p ε div ṽ ε dω =. 7)
11 PHILIPPE ANGOT AND PIERRE FABRIE Multiplying 4) by rε and 5) by ε and summing with 6) and 7), on obtains: d dt rε ṽε L +ε v ε L + v ε L +γε p ε L ) + rε curl ṽ ε L + ε curl v ε L + curl v ε L +γε p ε R L + rε div ṽ ε L R + ε div v ε L R + div v ε L R +r ṽ ε ν) L Γ) + v ε ν) L Γ) + ε v ε ν) L Γ) +εr div ṽ ε L +rε div v ε L + div v ε L εr = div v ε div v ε dω +rε f ṽ ε dω + f v ε dω, R Lt us now giv som bounds of th right-hand sid trms. Th trm Rdiv v ε,div v ε ) is boundd by div v ε L εr + ε div v ε L R. According to th stimat of th L norm in H ν) givn by quation ), w bound th sourc trms in th following way: f vdω f L v L λ f L div v L + curl v L ) div v L + curl v L ) + λ f L. Using ths bounds, w gt from th prvious quation th following fundamntal stimat: d rε ṽε L dt +ε v ε L + v ε L +γε p ε ) L +γε p ε L + rε curl ṽ ε L R + ε curl v ε L R + curl v ε L R + div v ε L + rε div ṽ ε L + ε div v ε L +εr div ṽ ε L +εr div v ε L + ε div v ε L +r ṽ ε ν) L Γ) + v ε ν) L Γ) + ε v ε ν) L Γ) λ +εr) f L. 8) Aftr intgration in tim, w dduc from th prvious stimat that thr xists a continuous function g dfind on [,T] such that: for all t >,
12 THE VECTOR PENALTY-PROJECTION METHOD ) rε ṽ ε t) +ε v ε t) + v ε t) L +γε p εt) L t +γε p ε s) rε t Lds+ curl ṽ ε s) L R ds+ ε t with + + ε +r + ε +εr t t t t t curl v ε s) rε Lds+ div v ε s) L ds+ t ṽ ε s) ν) L Γ) ds+ t t div ṽ ε s) L ds div v ε s) L ds v ε s) ν) L Γ) ds t v ε s) ν) L Γ) ds+εr div ṽ ε s) L ds div v ε s) L ds+ ε t div v ε s) Lds gt), gt) = rε ṽ ε ) +ε v ε ) + v ε ) L +γε p ε) L ) + λ +rε) This inquality is th ky point to stablish th convrgnc rsult. curl v ε s) L ds t 9) fs) L ds. To improv th convrgnc rsult in th two-dimnsional cas, w us th following nrgy quality drivd as abov without using 5). d dt rε ṽε L + v ε L +γε p ε ) rε L + curl ṽ ε L R + curl v ε L R +γε p ε L + rε div ṽ ε L R + div v ε L R +εr div ṽ ε L + div v ε L εr +r ṽ ε ν) L Γ) + ε v ε ν) L Γ) = rε f ṽ ε dω + f v ε dω. This concluds th proof of lmma.. ) 3.. Uniform stability for th approximat solution. In this sction, w dal with th stability of th proposd approximation mthod. W notic that this proprty is valid for all solutions satisfying th nrgy inquality 8), as it is th cas whn thy ar built by a finit dimnsional approximation mthod such as th Galrkin mthod for xampl. Proof. Lt us writ χ ε t) = εr ṽ ε t) L +ε v εt) L + v εt) L +γε p εt) L. W not λ > th smallst ignvalu of th slf-adjoint oprator A = curl curl div with th domain D = H ν) H )) d,
13 PHILIPPE ANGOT AND PIERRE FABRIE ) λ and w introduc α = min,. Classically, th inquality 8) lads to th diffrntial inquality d dt χ εt)+αχ ε t) λ +rε) ft) L, which implis th following uniform bound χ ε t) αt χ ε )+ λ α +rε) f L,L. This concluds th proof of thorm.. 4. Convrgnc analysis and compactnss rsults. 4.. Compactnss rsults for th vlocity. Lt us introduc th Lray projction w ε of a vlocity fild v ε t) H ν) dfind as follows v ε = w ε + q ε, div w ε =, w ε ν Γ =, q ε ν Γ =, q ε dω =. By th stimat 9), w s that th irrotational part of v ε gos to zro with ε. Thus it rmains to bring to th for th bhavior of th fr divrgnc part w ε and to obtain an stimat on a fractional tim drivativ of this trm. W dtail th diffrnt stps of this stratgy. From th rgularity of th Lray projctor s R. Tmam [] pag 8), on has: Morovr, w can asily prov th following lmma. w ε L,L c v ε L,L, w ε L,H c v ε L,H. ) Lmma 4.. Thr xists two constants dpnding only on T and such that: q ε L,H c ε, q ε L,L c ) Proof: Th function q ε blongs to H ) and satisfis This implis using th stimat 9) q ε t) = div v ε t), q ε t) ν Γ =. q ε t) L,L = div v ε L,L C ε. Bsids, w hav q ε q ε dω = q ε q ε dω, so that, with Poincaré-Numann inquality, w gt q ε t) L,L C q εt) L,L q ε L,L, C ε q ε t) L,L.
14 THE VECTOR PENALTY-PROJECTION METHOD 3 Th rgularity proprtis of th Numann problm giv q ε L,H C q ε L,H C q ε L,L C ε. 3) Morovr, by orthogonality of th Lray projctor in L, on has This concluds th proof of th lmma 4.. q ε L,L v ε L,L C. 4) So by intrpolation and using stimats 3)-4), w hav provd th rsult blow. Corollary. Th function q ε satisfis: q ε strongly convrgs to in L p ],T[;L )) ), p, p < +. Now w hav to writ th quation satisfid by w ε. As th Lray projction is orthogonal in L ), this quation rads ϕ H ν), div ϕ =, w ε ϕdω + t + = v ε )v ε + ) div v ε)v ε ϕdω v ε curl ϕdω + curl v ε ν) ϕ ν)dγ ε Γ f ϕdω in L,T). Now w introduc th xtnsion by of w ε rsp. v ε ) outsid [,T] dnotd, only in this part, by w ε rsp. ṽ ε ) and w tak th Fourir transform in tim of th quation 5) to obtain ϕ H ν), div ϕ =, iτ F w ε )τ) ϕdω + + = F ṽ ε )ṽ ε + ) div ṽ ε)ṽ ε τ) ϕdω Fṽ ε ν)τ) ϕ ν)dγ Fṽ ε )τ) curl ϕdω + curl ε Γ F f)τ) ϕdω + v ε ) ϕdω iτt π π v ε T) ϕdω. Following Boyr-Fabri [5, pag 53], w tak ϕ = F w ε )τ) as tst function in th prvious quation to obtain for all τ R: iτ F w ε )τ) ) dω = Fṽε )ṽ ε τ) F wε )τ)dω ) Fdiv ṽε )ṽ ε τ) F wε )τ)dω curl Fṽ ε )τ) curl F w ε )τ)dω R Fṽ ε ν)τ) F w ε )τ) ν)dγ ε Γ + F f)τ) F w ε )τ)dω + π v ε ) F w ε )τ)dω iτt π 5) v ε T) F w ε )τ)dω.
15 4 PHILIPPE ANGOT AND PIERRE FABRIE As w look for an stimat indpndnt of ε, w hav to pay a spcial attntion to th imaginary part of th pnalty trm: A ε = Fṽ ε ν)τ) F w ε )τ) ν)dγ. 6) ε Γ By writing w ε = v ε q ε, w hav: ε ν)τ) F w ε )τ) ν)dγ = ε ΓFṽ Fṽ ε )τ) ν dγ ε Γ F qε )τ) ν ) Fṽ ε )τ) ν) ) dγ. ε So, th imaginary part of A ε is boundd as follows: ε F qε )τ) ν ) Fṽ ε )τ) ν) ) dγ Γ From stimats ) and 9), w hav Γ ε F q ε)τ) L Γ) Fṽ ε)τ) ν L Γ), C ε F q ε)τ) H Fṽ ε )τ) ν L Γ). q ε L,H C ε, v ε ν L,L Γ) C ε. So thr xists a function fετ) 4 L R) boundd indpndntly on ε such that: ε F qε )τ) ν ) Fṽ ε )τ) ν) ) dγ f4 ετ) 8) Γ Now w can driv th stimat of τ F w ε )τ) dω, and w hav τ F w ε )τ) ) dω Fṽε )ṽ ε τ) F wε )τ)dω + ) Fdiv ṽε )ṽ ε τ) F wε )τ)dω + curl Fṽ ε )τ) curl F w ε )τ)dω R + ε F qε )τ) ν ) Fṽ ε )τ) ν) ) dγ Γ + Ff)τ) F w ε )τ)dω + v ε ) F w ε )τ)dω π + v ε T) F w ε )τ)dω π f ετ)+f ετ)+f 3 ετ)+f 4 ετ)+f 5 ετ)+f 6 ετ)+f 7 ετ). W now stimat ach trm of th right-hand sid of th prvious inquality for d 3. ) Trm fε = Fṽε )ṽ ε τ) F wε )τ)dω 7) 9)
16 THE VECTOR PENALTY-PROJECTION METHOD 5 Accordingtothnrgystimat9),thfunctionv ε isbounddinl,t;h ν)) and hnc, by Sobolv injction, it is boundd in L 6 5,T;L 6 )). So by Hausdorff- Young thorm, F w ε )) ε is boundd in L 6 R;L 6 )). ) W also hav th inquality: ṽ ε ṽ ε L 6 5 v ε L v ε L 3 C v ε 3 L v ε L, which implis, according to 9), that v ε v ε is boundd in L 4 3,T;L 6 5)), and ncssarily in L 6 5,T;L 6 5)). So, by Hausdorff-Young thorm, th family of functions Fṽ ε )ṽ ε ) is boundd in L 6 R;L 6 5)). Thn with Höldr inquality, Trm fε = Th sam argumnts show that Trm fε 3 = R fε) ε is boundd in L 3 R). ) ) Fdiv ṽε )ṽ ε τ) F wε )τ)dω fε) ε is boundd in L 3 R). ) curl Fṽ ε )τ) curl F w ε )τ)dω According to th rgularity of th Lray projction rcalld abov and stimat 9), on has: Trms fε 4 = ε Γ fε) 3 ε is boundd in L R). 3) F qε )τ) ν ) Fṽ ε )τ) ν) ) dγ As w hav sn by th stimat 8), Trm fε 5 = fε) 4 ε is boundd in L R). 4) F f)τ) F w ε )τ)dω By hypothsis, f is a givn function in L R;L )) and from stimat 9), w gt that ṽ ε is boundd in L R;L )), so as w ε is th Lray projction of ṽ ε, th function w ε is also boundd in L R;L )). Finally, w obtain Trms f 6 ε = π fε) 5 ε is boundd in L R). 5) v ε ) F w ε )τ)dω andf7 ε = π v εt) F w ε )τ)dω Ths two trms com from th Dirac masur whn w driv discontinuous functions. Lt us considr f 7 ε. f 7 ετ) π v ε T) L F w ε )τ) L. According to 9), this trm is boundd in L R). Morovr, th st of functions F wε )τ) ) ε is boundd in L R;L )) so, f 7 ε) ε is boundd in L R).
17 6 PHILIPPE ANGOT AND PIERRE FABRIE W trat in th sam way th trm f 6 ε, and thus: f 6 ε) ε and f 7 ε) ε ar boundd in L R) 6) Warnowabltoshowthatthstoffunctions w ε ) ε isbounddinanappropriat Nikolskii spac. For all γ <, thr xists a constant d such that: τ γ d + τ ) + τ γ, so that, τ γ F w ε )τ) L d F w ε )τ) L + τ ) + τ γ F w ε)τ) L Lt us dnot f ετ) = F w ε )τ) L, which blongs to L R), th prvious inquality rads with 9): τ γ F w ε )τ) L dfετ)+ d ) + τ γ fετ)+f ετ)+f ετ)+f 3 ετ)+f 4 ετ)+f 5 ετ)+f 6 ετ) 7 hτ) If w suppos that th function τ + τ γ blongs to L R) L R), thn th function τ hτ) blongs to L R). This condition is satisfid for γ ] 3,[. So w hav provd : Lmma 4.. Lt us suppos that σ ], 6 [, thn thr xists a constant C such that τ σ Fw ε )τ) Ldτ C. 7) R Thn, from lmma 4. and 4. w dduc th following ky rsult: Thorm 4.3. Thr xists a squnc ε k ) k which convrgs to zro and a function v L ],T[;L )) satisfying div v = such that: v εk ) k v in L ],T[;L )) strongly. Proof. Th function v ε is th sum of two trms q ε and w ε. From corollary th first trm convrgs strongly to in L ],T[;L )). Now, from Aubin-Lions- Simon Thorm, it follows from lmma 4., that thr xists a squnc ε k ) k such that: w εk ) k v in L ],T[;L )) strongly. Morovr, sinc div w εk =, w hav div v =. 4.. Convrgnc of th mthod. W first giv a gnral convrgnc thorm for a subsqunc solution of th approximat schm 3), to a wak solution of th initial Navir-Stoks problm, in th cas d 3. For th two-dimnsional cas, sinc th wak solution of th Navir-Stoks quation is uniqu, th whol squnc of approximat solution v ε convrgs to v. Morovr, in this cas, w prov that th convrgnc is strong.
18 THE VECTOR PENALTY-PROJECTION METHOD Th gnral cas d 3. Lt θ an lmnt of C,T), satisfying θt) =, and ϕ a fr-divrgnc vctor fild in H) ) d H )) d. An intgration by parts givs from th quation 3) T T v ε ϕdωθ τ)dτ + v ε )v ε + ) div v ε)v ε ϕdωθτ)dτ + T = T curl v ε curl ϕdωθτ)dτ f ϕdωθτ)dτ + v ϕdω θ). According to stimat 9), thr xists a squnc ε k such that v εk v in L,T;L )) strongly, div v εk div v = in L,T;L )) strongly, curl v εk curl v in L,T;L )) wakly. Sinc u L + div u L + curl u L ) is a norm quivalnt to th H -norm on H ν), w hav v εk v in L,T;L )) wakly. and so, w can tak th limit on th trm v ε )v ε as v εk v εk v v in L,T;L )) div vεk ) vεk in L,T;L )). 8) From stimat 9), w hav for th tangntial tracs t and sinc for any function v in H ) ) d, v ε ν)τ) L Γ) dτ εgt), v L Γ) C v L v H, w obtain that v εk v in L ],T[;L Γ)) d) strongly. This implis v ν) Γ =, and so, sinc by construction v ν) Γ =, v blongs to H ) ) d Finally, at th limit procss w obtain T T v ϕdωθ τ)dτ + v )v ϕdωθτ)dτ + T T curl v curl ϕdωθτ)dτ = f ϕdωθτ)dτ + v ϕdω θ), div v =, v Γ =. From th idntity v : ϕdω = curl v curl ϕdω + div v div ϕdω,
19 8 PHILIPPE ANGOT AND PIERRE FABRIE which is valid for all functions v,ϕ) H )) d H )) d, th prvious quality shows that th limit function v satisfis th classical Navir-Stoks quations in a wak sns Th spcial cas d =. Th ky point to stablish th strong convrgnc in th two-dimnsional cas, lis on an ida of R. Tmam []. It is basd on th fact that in this cas, th solution of th approximat problm and th solution of th Navir-Stoks quation vrify th quality of nrgy. Th ida is to bring to th for an nrgy quation satisfid by th diffrnc btwn th approximat solution and th xact solution. W first obsrv that, according to th quality = curl curl div th classical wak solution v of th Navir-Stoks quation satisfis for any tst function in H)) d, with fr-divrgnc: v t ϕdω + v )v ) ϕdω + curl v curl ϕdω = f ϕdω. Thquationsatisfidbythrrorv ε v = u ε withafr-divrgnctstfunction ϕ in H) ) d H ν ). u ε t ϕdω + + u ε )v + div u ε)v u ε )u ε + ) div u ε)u ε ϕdω + ) ϕdω + Aftr intgration in tim, this quation givs t u ε t) ϕdω + t t + v )u ε ϕdωdτ + + t v )u ε ϕdω curl u ε curl ϕdω =. u ε )u ε + div u ε)u ε ) ϕdωdτ curl u ε curl ϕdωdτ = u ε )v + ) div u ε)v ϕdωdτ u ε ) ϕdω. Taking th limit whn ε gos to, on obtains with th convrgncs proprtis statd in th prvious sction: Lmma 4.4. ϕ H)) d, div ϕ =, lim u ε t) ϕdω =. 9) ε Following R. Tmam [], w introduc χ ε t) = εr ṽε t) L +γε p εt) L + v ε v)t) ) L + rε t t curl ṽ ε τ) L R dτ +γε p ε τ) L dτ + t curl v ε v)τ) L R dτ + t div v ε v)τ) L R dτ +r + εr t t ṽ ε ν)τ) L Γ) dτ + ε t div ṽ ε τ) L dτ + ε t v ε ν)τ) L Γ) dτ div v ε τ) L dτ.
20 THE VECTOR PENALTY-PROJECTION METHOD 9 Hr v is th uniqu solution of th two-dimnsional Navir-Stoks quations. Du to th nrgy quality ), χ ε t) satisfis χ ε t) = εr ṽε ) L + v ε) L +γε p ε) ) t L +rε fτ) ṽ ε τ)dωdτ t + fτ) v ε τ)dωdτ v ε t) vt)dω t curl v ε τ) curl vτ)dωdτ t div v ε τ)div vτ)dωdτ + vt) L + t curl vτ) L R dτ + t div vτ) L R dτ. By wak convrgnc in L ],t[;h)) of th squnc v ε ) ε and from th lmma 4.4, w obsrv that : t lim ε v ε) L v ε t) vt)dω + fτ) v ε τ)dωdτ t curl v ε τ) curl vτ)dωdτ t ) div v ε τ) div vτ)dωdτ = t v) L + t div vτ) L R dτ. fτ) vτ)dωdτ vt) L t curl vτ) L R dτ In th two-dimnsional cas, th uniqu solution of th Navir-Stoks quation satisfis th following nrgy quality vt) L + t curl uτ) L dτ + t div uτ) L R dτ = t v) L + fτ) vτ)dωdτ. Morovr, from stimat ), t lim ε fτ) ṽ ε τ)dωdτ =. ε So, w hav provd that lim χ εt) =. ε In othr words, w hav stablishd th following rsult v ε v in C [,T];L )) L ],T[;H ν)) v C [,T];H) L ],T[;H )) This concluds th proof of thorm.3. W ar now abl to prcis th convrgnc for th ffctiv prssur and stablish th thorm.4 Lt us writ th quation satisfid by th vlocity v ε and th prssur p ε in th distribution sns. W hav v ε t +v ε )v ε + div v ε)v ε + curl curl v ε + p ε +ε +r ) ) div v ε = f. ε
21 PHILIPPE ANGOT AND PIERRE FABRIE Introducing th ffctiv prssur this quation rads q ε = f q ε = p ε +ε ε +r ) div v ε, vε t +v ε )v ε + div v ε)v ε + ) curl curl v ε and th proof follows, from th prvious stps. 5. Appndix. Lt us considr th following problm: Proposition 4. For a fixd ε > and a coupl of functions f,g) givn in L ) L ), thr xists a uniqu solution v ε {w H div ),w ν Γ = } solution of: εv ε div v ε = f +εg. 3) Morovr, if u,u ) H div ) H is solution of div u = f,curl u =, u ν Γ =, 3) div u =,curl u = curl g, u ν Γ =, and w hav th following stimat v ε u u H ε u+u g L. Proof. Stp : Existnc of v ε Lt us not H div, ) = {v L ), div v L, v ν Γ = }. Th xistnc of a uniqu solution to th quation 3) is obtaind by a straightforward application of th Lax-Milgram thorm with th bilinar form dfind on H div, ) by εu,v)+div u,div v), and th right-hand sid: f,div v)+εg,v). Stp : xistnc of u and u Th xistnc of u satisfying 3) coms from th rsolution of th following Numann problm q = f, q ν Γ =, and w st u = q, with q H )/R ) H ). Th xistnc of u H is th consqunc of th Lray projction applid to g by writing g = u + p, div u =, u ν Γ =. Now, writting v ε = u+u +u ε, w gt that u ε H div, satisfis εu ε div u ε = εu εu +εg, curl u ε =,
22 THE VECTOR PENALTY-PROJECTION METHOD and w hav th stimat ε u ε L + div u ε L ε u+u g L u ε L. W obsrv that, according to [8], H ν) ndowd with th norm div v + curl v ) is qual to {w H )) d,w ν γ }. So th prvious stimat givs ε u ε L + div u ε L + curl u ε L ε u+u g L u ε H ν Cε u+u g L div uε L + curl u ε ) L, which implis using Young inquality and th proof of proposition 4 follows. ε u ε L + div u ε L ε u+u g L, Rmark 4.. Th function u blongs to H ν), which is not th cas for v ε or u, without som additional rgularity hypothss on th function g. Nvrthlss, th function u ε = v ε u u blongs to H ν).. In th cas whr g =, th function v ε blongs to H ν). REFERENCES [] Ph. Angot, J.-P. Caltagiron and P. Fabri, Vctor pnalty-projction mthod for th solution of unstady incomprssibl flows, in Finit Volums for Complx Applications V - Problms & Prspctivs ds R. Eymard and J.-M. Hérard), ISTE Ltd and J. Wily & Sons, 8), [] Ph. Angot, J.-P. Caltagiron and P. Fabri, A spctacular vctor pnalty-projction mthod for Darcy and Navir-Stoks problms, in Finit Volums for Complx Applications VI - Problms & Prspctivs ds J. Fořt t al.), Intrnational Symposium FVCA6 in Pragu, Jun 6-, Springr Procdings in Mathmatics 4, Vol., Springr-Vrlag Brlin), ), [3] Ph. Angot, J.-P. Caltagiron and P. Fabri, A nw fast mthod to comput saddl-points in constraind optimization and applications, Applid Mathmatics Lttrs, 53) ), [4] Ph. Angot, J.-P. Caltagiron and P. Fabri, A fast vctor pnalty-projction mthod for incomprssibl non-homognous or multiphas Navir-stoks problms, Applid Mathmatics Lttrs, in prss), DOI:.6/j.aml [5] F. Boyr and P. Fabri, Élémnts d analys pour l étud d qulqus modèls d écoulmnts d fluids visquux incomprssibls, Mathématiqus & Applications 5, Springr-Vrlag, 6. [6] A.J. Chorin, A numrical mthod for solving incomprssibl viscous flow problms, J. Comput. Phys., 967), 6. [7] A.J. Chorin, Numrical solution of th Navir-Stoks quations, Math. Comput., 968), [8] C. Foias and R. Tmam, Rmarqus sur ls équations d Navir-Stoks stationnairs t ls phénomèns succssifs d bifurcation, Annali dlla Scuolo Normal Suprior di Pisa, Class di Scinz 4 séri, 5) 978), [9] V. Girault and P.A. Raviart, Finit Elmnt Mthods for th Navir-Stoks Equations, Springr Sris in Comput. Math. 5, Springr-Vrlag Nw York), 986. [] J.-L. Gurmond, P.D. Minv and J. Shn, An ovrviw of projction mthods for incomprssibl flows, Comput. Mth. Appl. Mch. Engrg., 95 6), [] O.A. Ladyzhnskaya, Th mathmatical thory of viscous incomprssibl flow, nd dition, Gordon and Brach Nw York), 969.
23 PHILIPPE ANGOT AND PIERRE FABRIE [] J. Lray, Essai sur ls mouvmnts plans d un liquid visquux qu limitnt ds parois, J. Math. Purs Appl., 3 934), [3] J.-L. Lions, Qulqus méthods d résolution ds problèms aux limits non linéairs, Dunod & Gauthir-Villars Paris), 969. [4] J. Shn, On rror stimats of th pnalty mthod for unstady Navir-Stoks quations, SIAM J. Numr. Anal., 3) 995), [5] J. Shn, On a nw psudocomprssibility mthod for th incomprssibl Navir-Stoks quations, Appl. Numr. Math., 996), 7 9. [6] J. Simon, Compact sts in th spac L p,t;b), Ann. Mat. Pura Appl., 464) 987), [7] R. Tmam, Un méthod d approximation d la solution ds équations d Navir-Stoks, Bull. Soc. Math. Franc, ), 5 5. [8] R. Tmam, Sur l approximation d la solution ds équations d Navir-Stoks par la méthod ds pas fractionnairs I, Arch. Ration. Mch. Anal., 3) 969), [9] R. Tmam, Sur l approximation d la solution ds équations d Navir-Stoks par la méthod ds pas fractionnairs II, Arch. Ration. Mch. Anal., 335) 969), [] R. Tmam, Navir-Stoks Equations; Thory and Numrical Analysis, 3 rd dition, North- Holland Amstrdam), 984. Rcivd 3/7 ; rvisd 7/. addrss: angot@cmi.univ-mrs.fr addrss: pirr.fabri@math.u-bordaux.fr
CONVERGENCE 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 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 informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by Dan Klain Vrsion 28928 Corrctions and commnts ar wlcom Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix () A A k I + A + k!
More informationA 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 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 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 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 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 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 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 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 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 informationLimiting value of higher Mahler measure
Limiting valu of highr Mahlr masur Arunabha Biswas a, Chris Monico a, a Dpartmnt of Mathmatics & Statistics, Txas Tch Univrsity, Lubbock, TX 7949, USA Abstract W considr th k-highr Mahlr masur m k P )
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 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 informationDiscontinuous Galerkin and Mimetic Finite Difference Methods for Coupled Stokes-Darcy Flows on Polygonal and Polyhedral Grids
Discontinuous Galrkin and Mimtic Finit Diffrnc Mthods for Coupld Stoks-Darcy Flows on Polygonal and Polyhdral Grids Konstantin Lipnikov Danail Vassilv Ivan Yotov Abstract W study locally mass consrvativ
More informationRecall that by Theorems 10.3 and 10.4 together provide us the estimate o(n2 ), S(q) q 9, q=1
Chaptr 11 Th singular sris Rcall that by Thorms 10 and 104 togthr provid us th stimat 9 4 n 2 111 Rn = SnΓ 2 + on2, whr th singular sris Sn was dfind in Chaptr 10 as Sn = q=1 Sq q 9, with Sq = 1 a q gcda,q=1
More 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 informationBifurcation Theory. , a stationary point, depends on the value of α. At certain values
Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local
More informationChapter 10. The singular integral Introducing S(n) and J(n)
Chaptr Th singular intgral Our aim in this chaptr is to rplac th functions S (n) and J (n) by mor convnint xprssions; ths will b calld th singular sris S(n) and th singular intgral J(n). This will b don
More 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 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 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 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 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 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 information10. The Discrete-Time Fourier Transform (DTFT)
Th Discrt-Tim Fourir Transform (DTFT Dfinition of th discrt-tim Fourir transform Th Fourir rprsntation of signals plays an important rol in both continuous and discrt signal procssing In this sction w
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 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 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 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 informationProperties of Phase Space Wavefunctions and Eigenvalue Equation of Momentum Dispersion Operator
Proprtis of Phas Spac Wavfunctions and Eignvalu Equation of Momntum Disprsion Oprator Ravo Tokiniaina Ranaivoson 1, Raolina Andriambololona 2, Hanitriarivo Rakotoson 3 raolinasp@yahoo.fr 1 ;jacqulinraolina@hotmail.com
More informationA Simple Formula for the Hilbert Metric with Respect to a Sub-Gaussian Cone
mathmatics Articl A Simpl Formula for th Hilbrt Mtric with Rspct to a Sub-Gaussian Con Stéphan Chrétin 1, * and Juan-Pablo Ortga 2 1 National Physical Laboratory, Hampton Road, Tddinton TW11 0LW, UK 2
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 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 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 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 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 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 informationSCHUR S THEOREM REU SUMMER 2005
SCHUR S THEOREM REU SUMMER 2005 1. Combinatorial aroach Prhas th first rsult in th subjct blongs to I. Schur and dats back to 1916. On of his motivation was to study th local vrsion of th famous quation
More informationCOUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM
COUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM Jim Brown Dpartmnt of Mathmatical Scincs, Clmson Univrsity, Clmson, SC 9634, USA jimlb@g.clmson.du Robrt Cass Dpartmnt of Mathmatics,
More 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 informationResearch Article Norm and Essential Norm of an Integral-Type Operator from the Dirichlet Space to the Bloch-Type Space on the Unit Ball
Hindawi Publishing Corporation Abstract and Applid Analysis Volum 2010, Articl ID 134969, 9 pags doi:10.1155/2010/134969 Rsarch Articl Norm and Essntial Norm of an Intgral-Typ Oprator from th Dirichlt
More informationSome remarks on Kurepa s left factorial
Som rmarks on Kurpa s lft factorial arxiv:math/0410477v1 [math.nt] 21 Oct 2004 Brnd C. Kllnr Abstract W stablish a connction btwn th subfactorial function S(n) and th lft factorial function of Kurpa K(n).
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 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 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 informationQuasi-Classical States of the Simple Harmonic Oscillator
Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats
More informationOn 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 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 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 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 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 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 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 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 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 informationCalculus II (MAC )
Calculus II (MAC232-2) Tst 2 (25/6/25) Nam (PRINT): Plas show your work. An answr with no work rcivs no crdit. You may us th back of a pag if you nd mor spac for a problm. You may not us any calculators.
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 informationECE602 Exam 1 April 5, You must show ALL of your work for full credit.
ECE62 Exam April 5, 27 Nam: Solution Scor: / This xam is closd-book. You must show ALL of your work for full crdit. Plas rad th qustions carfully. Plas chck your answrs carfully. Calculators may NOT b
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 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 informationSpectral Synthesis in the Heisenberg Group
Intrnational Journal of Mathmatical Analysis Vol. 13, 19, no. 1, 1-5 HIKARI Ltd, www.m-hikari.com https://doi.org/1.1988/ijma.19.81179 Spctral Synthsis in th Hisnbrg Group Yitzhak Wit Dpartmnt of Mathmatics,
More informationarxiv: v1 [math.na] 3 Mar 2016
MATHEMATICS OF COMPUTATION Volum 00, Numbr 0, Pags 000 000 S 0025-5718(XX)0000-0 arxiv:1603.01024v1 [math.na] 3 Mar 2016 RESIDUAL-BASED A POSTERIORI ERROR ESTIMATE FOR INTERFACE PROBLEMS: NONCONFORMING
More 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 informationMath 102. Rumbos Spring Solutions to Assignment #8. Solution: The matrix, A, corresponding to the system in (1) is
Math 12. Rumbos Spring 218 1 Solutions to Assignmnt #8 1. Construct a fundamntal matrix for th systm { ẋ 2y ẏ x + y. (1 Solution: Th matrix, A, corrsponding to th systm in (1 is 2 A. (2 1 1 Th charactristic
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 information2008 AP Calculus BC Multiple Choice Exam
008 AP Multipl Choic Eam Nam 008 AP Calculus BC Multipl Choic Eam Sction No Calculator Activ AP Calculus 008 BC Multipl Choic. At tim t 0, a particl moving in th -plan is th acclration vctor of th particl
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 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 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 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 informationNEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA
NE APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA Mirca I CÎRNU Ph Dp o Mathmatics III Faculty o Applid Scincs Univrsity Polithnica o Bucharst Cirnumirca @yahoocom Abstract In a rcnt papr [] 5 th indinit intgrals
More informationEwald s Method Revisited: Rapidly Convergent Series Representations of Certain Green s Functions. Vassilis G. Papanicolaou 1
wald s Mthod Rvisitd: Rapidly Convrgnt Sris Rprsntations of Crtain Grn s Functions Vassilis G. Papanicolaou 1 Suggstd Running Had: wald s Mthod Rvisitd Complt Mailing Addrss of Contact Author for offic
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 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 information10. Limits involving infinity
. Limits involving infinity It is known from th it ruls for fundamntal arithmtic oprations (+,-,, ) that if two functions hav finit its at a (finit or infinit) point, that is, thy ar convrgnt, th it of
More informationWeek 3: Connected Subgraphs
Wk 3: Connctd Subgraphs Sptmbr 19, 2016 1 Connctd Graphs Path, Distanc: A path from a vrtx x to a vrtx y in a graph G is rfrrd to an xy-path. Lt X, Y V (G). An (X, Y )-path is an xy-path with x X and y
More informationMATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations
MATH 39, WEEK 5: Th Fundamntal Matrix, Non-Homognous Systms of Diffrntial Equations Fundamntal Matrics Considr th problm of dtrmining th particular solution for an nsmbl of initial conditions For instanc,
More informationFinite Element Models for Steady Flows of Viscous Incompressible Fluids
Finit Elmnt Modls for Stad Flows of Viscous Incomprssibl Fluids Rad: Chaptr 10 JN Rdd CONTENTS Govrning Equations of Flows of Incomprssibl Fluids Mid (Vlocit-Prssur) Finit Elmnt Modl Pnalt Function Mthod
More informationExercise 1. Sketch the graph of the following function. (x 2
Writtn tst: Fbruary 9th, 06 Exrcis. Sktch th graph of th following function fx = x + x, spcifying: domain, possibl asymptots, monotonicity, continuity, local and global maxima or minima, and non-drivability
More informationBSc Engineering Sciences A. Y. 2017/18 Written exam of the course Mathematical Analysis 2 August 30, x n, ) n 2
BSc Enginring Scincs A. Y. 27/8 Writtn xam of th cours Mathmatical Analysis 2 August, 28. Givn th powr sris + n + n 2 x n, n n dtrmin its radius of convrgnc r, and study th convrgnc for x ±r. By th root
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 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 informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 401 Digital Signal Procssing Prof. Mark Fowlr Dtails of th ot St #19 Rading Assignmnt: Sct. 7.1.2, 7.1.3, & 7.2 of Proakis & Manolakis Dfinition of th So Givn signal data points x[n] for n = 0,, -1
More informationA C 0 INTERIOR PENALTY METHOD FOR A FOURTH ORDER ELLIPTIC SINGULAR PERTURBATION PROBLEM
A C 0 INERIOR PENALY MEHOD FOR A FOURH ORDER ELLIPIC SINGULAR PERURBAION PROBLEM SUSANNE C. BRENNER AND MICHAEL NEILAN Abstract. In tis papr, w dvlop a C 0 intrior pnalty mtod for a fourt ordr singular
More information(1) Then we could wave our hands over this and it would become:
MAT* K285 Spring 28 Anthony Bnoit 4/17/28 Wk 12: Laplac Tranform Rading: Kohlr & Johnon, Chaptr 5 to p. 35 HW: 5.1: 3, 7, 1*, 19 5.2: 1, 5*, 13*, 19, 45* 5.3: 1, 11*, 19 * Pla writ-up th problm natly and
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 information2F1120 Spektrala transformer för Media Solutions to Steiglitz, Chapter 1
F110 Spktrala transformr för Mdia Solutions to Stiglitz, Chaptr 1 Prfac This documnt contains solutions to slctd problms from Kn Stiglitz s book: A Digital Signal Procssing Primr publishd by Addison-Wsly.
More informationFourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.
Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform
More 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 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 informationSelf-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016
Ronald I. Frank 06 Adjoint https://n.wikipdia.org/wiki/adjoint In gnral thr is an oprator and a procss that dfin its adjoint *. It is thn slf-adjoint if *. Innr product spac https://n.wikipdia.org/wiki/innr_product_spac
More informationSection 11.6: Directional Derivatives and the Gradient Vector
Sction.6: Dirctional Drivativs and th Gradint Vctor Practic HW rom Stwart Ttbook not to hand in p. 778 # -4 p. 799 # 4-5 7 9 9 35 37 odd Th Dirctional Drivativ Rcall that a b Slop o th tangnt lin to th
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 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 informationLecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields
Lctur 37 (Schrödingr Equation) Physics 6-01 Spring 018 Douglas Filds Rducd Mass OK, so th Bohr modl of th atom givs nrgy lvls: E n 1 k m n 4 But, this has on problm it was dvlopd assuming th acclration
More informationPROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS
Intrnational Journal Of Advanc Rsarch In Scinc And Enginring http://www.ijars.com IJARSE, Vol. No., Issu No., Fbruary, 013 ISSN-319-8354(E) PROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS 1 Dharmndra
More information1 Minimum Cut Problem
CS 6 Lctur 6 Min Cut and argr s Algorithm Scribs: Png Hui How (05), Virginia Dat: May 4, 06 Minimum Cut Problm Today, w introduc th minimum cut problm. This problm has many motivations, on of which coms
More informationAn Extensive Study of Approximating the Periodic. Solutions of the Prey Predator System
pplid athmatical Scincs Vol. 00 no. 5 5 - n xtnsiv Study of pproximating th Priodic Solutions of th Pry Prdator Systm D. Vnu Gopala Rao * ailing addrss: Plot No.59 Sctor-.V.P.Colony Visahapatnam 50 07
More information