arxiv: v1 [math.na] 22 Oct 2010

Transcription

1 arxiv: v1 [matna] 22 Oct 2010 ABSOLUTELY STABLE LOCAL DSCONTNUOUS GALERKN METHODS FOR THE HELMHOLTZ EQUATON WTH LARGE WAVE NUMBER XAOBNG FENG AND YULONG XNG Abstract Two local discontinuous Galrkin (LDG mtods using som nonstandard numrical fluxs ar dvlopd for t Hlmoltz quation wit t first ordr absorbing boundary condition in t ig frquncy rgim t is sown tat t proposd LDG mtods ar absolutly stabl (nc wllposd wit rspct to bot t wav numbr and t ms siz Optimal ordr (wit rspct to t ms siz rror stimats ar provd for all wav numbrs in t prasymptotic rgim To analyz t proposd LDG mtods, ty ar rcastd and tratd as (non-conforming mixd finit lmnt mtods T crux of t analysis is to stablis a gnralizd inf-sup condition, wic olds witout any ms constraint, for ac LDG mtod T gnralizd inf-sup conditions tn asily infr t dsird absolut stability of t proposd LDG mtods n rturn, t stability rsults not only guarant t wll-posdnss of t LDG mtods but also play a crucial rol in t drivation of t rror stimats Numrical xprimnts, wic confirm t tortical rsults and compar t proposd two LDG mtods, ar also prsntd in t papr 1 ntroduction Tis papr is t tird installmnt in a sris [10, 11] wic dvot to dvloping absolutly stabl discontinuous Galrkin (DG mtods for t following prototypical Hlmoltz problm wit larg wav numbr: (11 (12 u k 2 u = f in Ω R d, d = 2, 3, u n Ω + ıku = g on Γ = Ω, wr ı = 1 dnots t imaginary unit k R + is a givn positiv (larg numbr and known as t wav numbr (12 is t so-calld first ordr absorbing boundary condition [9] W rcall tat [10, 11] focusd on dsigning and analyzing - and p-intrior pnalty discontinuous Galrkin (PDG mtods wic ar absolutly stabl (wit rspct to wav numbr k and ms siz and optimally convrgnt (wit rspct to T main idas of [10, 11] ar to introduc som novl intrior pnalty trms 1991 Matmatics Subjct Classification 65N12, 65N15, 65N30, 78A40 Ky words and prass Hlmoltz quation, tim armonic wavs, local discontinuous Galrkin mtods, stability, rror stimats T work of t first autor was partially supportd by t NSF grants DMS and DMS T rsarc of t scond autor was partially sponsord by t Offic of Advancd Scintific Computing Rsarc; US Dpartmnt of Enrgy T work of t scond autor was prformd at t ORNL, wic is managd by UT-Battll, LLC undr Contract No DE-AC05-00OR

2 2 XAOBNG FENG AND YULONG XNG in t ssquilinar forms of t proposd PDG mtods and to us a non-standard analytical tool, wic is basd on a Rllic idntity tcniqu, to prov t dsird stability and rror stimats T numrical xprimnt rsults sown tat t absolutly stabl PDG mtods significantly outprform t standard finit lmnt and finit diffrnc mtods, wic ar known only to b stabl undr stringnt ms constraints k 1 or k 2 1 (cf [8, 16], for t Hlmoltz problm Morovr, t numrical xprimnt rsults also sown tat ts PDG mtods ar capabl to corrctly track t pass of t igly oscillatory wavs vn wn t ms violats t rul-of-tumb condition (i, 6 10 grid points must b usd in a wav lngt T main difficulty of analyzing t Hlmoltz typ problms is causd by t strong indfinitnss of t Hlmoltz quation wic in turn maks it ard to stablis stability stimats for its numrical approximations T loss of stability in t cas of larg wav numbrs rsults in an additional pollution rror (bsids t intrpolation rror in t global rror bounds Extnsiv rsarc as bn don to addrss t qustion wtr it is possibl to rduc t pollution ffct, w rfr t radr to Captr 4 of [15] and t rfrncs trin for an dtaild xposition in tis dirction Motivatd by t succss of [10, 11], t primary objctiv of tis papr is to xtnd t work of [10, 11] to t local discontinuous Galrkin (LDG formulation, wic is known to b mor pysical and flxibl tan t PDG formulation on dsigning DG scms [1, 5] As it is wll-known now, t ky stp for constructing LDG mtods is to dsign t numrical fluxs As soon as t numrical fluxs ar slctd, for a larg class of corciv lliptic and parabolic scond ordr problms, tr is a gnral framwork for carrying out convrgnc analysis of LDG mtods [1] Unfortunatly, tis gnral framwork dos not apply to t Hlmoltz typ problms wic is xtrmly noncorciv/indfinit for larg wav numbr k Nvrtlss, wn dsigning t numrical fluxs for our LDG mtods, w borrow t ida of [1] by stablising t connction btwn our LDG mtods and t PDG mtods of [10, 11] altoug it turns out tat t PDG mtods of [10, 11] do not av xactly quivalnt LDG formulations du to t non-standard pnalty trms usd in [10, 11] Tis tn lads to t construction of our first LDG mtod t is provd and numrically vrifid tat tis LDG mtod is absolutly stabl and optimally convrgnt for t scalar variabl Howvr, it is sub-optimal for t vctor/flux variabl To improv t approximation accuracy for t vctor/flux variabl, w dsign anotr st of numrical fluxs wic rsult in t construction of our scond LDG mtod t is provd tat t scond LDG mtod is also absolutly stabl and givs a bttr approximation for t vctor/flux variabl tan t first mtod On t otr and, it is computationally mor xpnsiv tan t first LDG mtod, wic is xpctd To analyz t proposd LDG mtods, w tak an opposit approac to tat advocatd in [1], tat is, instad of convrting LDG mtods to tir quivalnt PDG mtods in t primal form, w rcast and trat our LDG mtods as nonconforming mixd finit lmnt mtods To avoid using t standard tcniqus suc as Scatz argumnt (cf [2, 8] or Babuška s inf-sup condition argumnt [16] to driv rror stimats (and to prov stability, bot approacs would crtainly lad to stringnt ms constraints, our main ida is to stablis a gnralizd inf-sup condition, wic olds witout any ms constraint, for ac LDG mtod T gnralizd inf-sup conditions tn immdiatly infr t dsird absolut stability

3 LDG METHODS FOR THE HELMHOLTZ EQUATON 3 of t proposd LDG mtods n rturn, t stability rsults not only guarant t wll-posdnss of t LDG mtods but also play a crucial rol in t drivation of t (optimal rror stimats t sould b pointd out tat a lot of work as rcntly bn don on dvloping DG mtods using picwis plan wav functions, oppos to simplr picwis polynomial functions as don in tis papr, for t Hlmoltz typ problms Howvr, to t bst of our knowldg, non of ts plan wav DG mtod is provd to b absolutly stabl wit rspct to wav numbr k and ms siz W rfr t radr to [12, 14, 17] and t rfrncs trin for mor discussions in tis dirction W also rfr to [10, 11] for mor discussions and rfrncs on otr discrtization tcniqus for t Hlmoltz typ problms Tis papr consists of four additional sctions n Sction 2, w introduc t notations usd in tis papr and prsnt t drivations of our two LDG mtods n Sction 3, w prsnt a dtaild stability analysis for bot LDG mtods T main task of t sction is to prov a gnralizd inf-sup condition for ac proposd LDG mtod Similar to [10, 11], a nonortodox tst function trick is t ky to gt t job don n Sction 4, a non-standard two-stp rror stimat procdur is usd to driv rror stimats for t proposd LDG mtods Onc again, t stability stimats stablisd in Sction 3 play a crucial rol Finally, Sction 5 contains som numrical xprimnts wic ar dsignd to vrify t tortical rror bounds provd in Sction 4 and to compar t prformanc of t proposd two LDG mtods 2 Formulation of local discontinuous Galrkin mtods T standard spac, norm and innr product notation ar adoptd in tis papr Tir dfinitions can b found in [1, 2, 4, 19] n particular, (, Q and, Σ for Σ Q dnot t L 2 -innr product on complx-valud L 2 (Q and L 2 (Σ spacs, rspctivly (, := (, Ω and, :=, Ω Trougout t papr, C is usd to dnot a gnric positiv constant wic is indpndnt of and k W also us t sortand notation A B and B A for t inquality A CB and B CA A B is a sortand notation for t statmnt A B and B A Assum tat Ω R d (d = 2, 3 is a boundd and strictly star-sapd domain wit rspct to a point x Ω Ω W now rcall t dfinition of star-sapd domains Dfinition 21 Q R d is said to b a star-sapd domain wit rspct to x Q Q if tr xists a nonngativ constant c Q suc tat (21 (x x Q n Q c Q x Q Q R d is said to b strictly star-sapd if c Q is positiv Lt T b a family of partitions of Ω paramtrizd by > 0 For any triangl/ttradron K T, w dfin K := diam(k and := max K T K Similarly, for ac dg/fac of K T, dfin := diam( W assum tat t lmnts of T satisfy t minimal angl condition W also dfin E := st of all intrior dgs/facs of T, E B := st of all boundary dgs/facs of T on Γ = Ω, E := E E B

4 4 XAOBNG FENG AND YULONG XNG Lt b an intrior dg sard by two lmnts K 1 and K 2 wos unit outward normal vctors ar dnotd by n 1 and n 2 For a scalar function v, lt v i = v Ki, and dfin {v} = 1 2 (v 1 + v 2, [v] = v K v K, [[v]] = v 1 n 1 + v 2 n 2 on E, wr K is K 1 or K 2, wicvr as t biggr global labling and K is t otr For a vctor fild v, lt v i = v Ki and dfin {v} = 1 2 (v 1 + v 2, [v] = v K v K, [[v]] = v 1 n 1 + v 2 n 2 on E As it is wll-known now (cf [1] tat t first stp for formulating an LDG mtod is to rwrit t givn PDE as a first ordr systm by introducing an auxiliary variabl For t Hlmoltz problm (11 (12 w av (22 (23 (24 σ = u in Ω, divσ k 2 u = f in Ω, u n Ω + ıku = g on Γ, Clarly, t vctor-valud function (oftn calld t flux variabl σ is t auxiliary variabl Tn, multiplying (22 and (23 by tst functions τ and v, rspctivly, and intgrating bot quations ovr an lmnt K T yilds (25 σ τ dx = u divτ dx + u n K τ ds, K K K (26 σ v dx k 2 uv dx = fv dx + σ n K v ds, K K K wr n K dnots t unit outward normal vctor to K T abov quations form t wak formulation on uss to dfin LDG mtods for t Hlmoltz problm (11 (12 Nxt, w dfin LDG spacs as follows V := {v L 2 (Ω; R(v K, m(v P r (K K T }, Σ := {τ (L 2 (Ω d ; R(τ K, m(τ (P l (K d K T }, wr P r (K (r 1 stands for t st of all polynomials of dgr lss tan or qual to r on K Finally, w ar rady to dfin t following gnral LDG formulation: Find (u, σ V Σ suc tat for all K Γ tr old (27 σ τ dx = u divτ dx + û K n K τ ds, K K K (28 σ v dx k 2 u v dx = fv dx + ˆσ K n K v ds K K for any (v, τ V Σ Wr t quantitis û K and ˆσ K, wic ar calld numrical fluxs, ar rspctivly approximations to σ = u and u on t boundary K of K As it is wll-known now tat t most important issu for all LDG mtods is ow to coos t numrical fluxs T diffrnt coics of t numrical fluxs obviously lad to diffrnt LDG mtods t is asy to undrstand tat K K K

5 LDG METHODS FOR THE HELMHOLTZ EQUATON 5 ts numrical fluxs must b cosn carfully in ordr to nsur t stability and accuracy of t rsultd LDG mtods n tis papr, w sall only considr t linar lmnt cas (i, r = l = 1 and propos two sts of numrical fluxs (û K, ˆσ K, wic lad to two LDG mtods Our coics of numrical fluxs ar inspird by t intrior pnalty discontinuous Galrkin (PDG mtods proposd by Fng and Wu [10] and ar idntifid wit t lp of t unifid DG framwork of [1] wic bridgs t primal DG formulations (g PDG mtods and t flux DG formulations (i, LDG mtods (1 LDG mtod #1: St ˆσ K = { u } ıβ[[u ]], û K = {u } + ıδ[[ u ]] on E, ˆσ K = ıku n K + gn K, û K = u on E B (2 LDG mtod #2: St ˆσ K = {σ } ıβ[[u ]], û K = {u } + ıδ[[σ ]] on E, ˆσ K = ıku n K + gn K, û K = u on E B Wr β and δ ar positiv constants to b spcifid latr and dnots t picwisly dfind gradint oprator ovr T, tat is, K = K K T For t radr s convninc, w now brifly sktc t drivation of t primal DG formulation corrsponding to our LDG mtod #1 by adapting t drivation givn in t gnral framwork of [1] Substituting t numrical fluxs of LDG mtod #1 into (27 and (28, summing t rsulting quations ovr all lmnt K T and using t following intgration by parts idntity w gt (u, divτ Ω = ( u, τ Ω + E ( {u }, [[τ ]] + [[u ]], {τ } + u, n K τ Γ, (29 (210 (σ, τ Ω ( u, τ Ω (ıδ [[ u ]], [[τ ]] [[u ]], {τ } = 0, E (σ, v Ω k 2 (u, v Ω { u } ıβ[[u ]], [[v ]] E + ık u, v Γ = (f, v Ω + g, v Γ Stting τ = v in (29 and subtracting t rsulting quation from (210 tn lads to t following formulation: (211 A (u, v k 2 (u, v Ω = F (v v V,

6 6 XAOBNG FENG AND YULONG XNG wr A (u, v := ( u, v Ω + ık u, v Γ + ı (δ [[ u ]], [[ v ]] + β [[u ]], [[v ]] E ( [[u ]], { v } + { u }, [[v ]] E F (v := (f, v Ω + g, v Γ Hnc, (211 is t corrsponding primal (i, PDG formulation of our LDG mtod #1 Comparing (211 wit t PDG formulation of [10], w s tat all t intrior pnalty trms of (211 also appar in t PDG formulation of [10] n fact, it is xactly by rvrsing t ordr of t abov drivation tat lads to t discovry of t numrical fluxs of LDG mtod #1 3 Stability analysis From tir constructions, it is asy to s tat bot LDG mtods proposd in t prvious sction ar consistnt scms for t Hlmoltz problm (11 (12 For corciv lliptic (and parabolic problms, t stability of suc a numrical scm can b provd asily as dmonstratd in [1] (t sam statmnt is tru for tir corrsponding PDE stability analysis Howvr, t Hlmoltz problm (11 (12 is an indfinit problm and it bcoms notoriously non-corciv for larg wav numbr k Driving its stability stimats (i, a priori stimats of its PDE solution, particularly wav-numbr-dpndnt stimats, as bn provd not to b an asy job (cf [6, 7, 13, 18] and t rfrncs trin Numrically, suc a qust as bn known to b vn ardr bcaus of t lowr ordr of t smootnss and t inflxibility of (picwis approximation functions (cf [6, 16, 18] and t rfrncs trin T stability of t numrical mtods in all t abov quotd rfrncs was provd undr som vry rstrictiv ms constraints An opn qustion was tn raisd by Zinkiwicz [20] wic asks wtr it is possibl to construct an absolutly stabl (and optimally convrgnt numrical mtod (i, no rstriction on t ms siz and t wav numbr k for t Hlmoltz quation Almost a dcad latr Fng and Wu [10, 11] wr abl to dsign for t first tim suc numrical mtods, wic appn to blong to t PDG family, for t Hlmoltz problm (11 (12 T goal of tis sction is to sow tat in t cas of t linar lmnt (i, r = l = 1 t LDG mtod #1 and #2 proposd in t prvious sction av t sam stability proprty as t PDG mtods of [10, 11] do, tat is, t LDG mtod #1 and #2 ar absolutly stabl for all ms siz > 0 and all wav numbr k > 0 witout imposing any constraint on tm To stablis tis rsult, unlik t approac usd and advocatd in t gnral framwork of [1] wic convrts an LDG mtod into its quivalnt primal mtod and tn analyzs t lattr using t standard finit lmnt Galrkin tcniqus, w sall dirctly fit t LDG mtod #1 and #2 into t (nonconforming mixd mtod framwork as don in [16] and adapting t mixd mtod tcniqus to prov t dsird stability rsult t is wll-known tat t ky ingrdint for t mixd mtod approac is to stablis t inf-sup or Babuška-Brzzi condition for t (augmntd ssquilinar forms for ac mtod Howvr, w ar not abl to prov suc an inf-sup condition for itr

7 LDG METHODS FOR THE HELMHOLTZ EQUATON 7 mtod witout imposing ms constraints (wic w bliv is not possibl To ovrcom t difficulty, our main ida is to prov a gnralizd (and wakr inf-sup condition wic olds for all, k > 0 t turns out tat tis gnralizd inf-sup condition is sufficint for us to stablis t dsird absolut stability for t LDG mtod #1 and #2 To prov t gnralizd inf-sup condition, t ky tcniqu w us is a spcial tst function tcniqu wic was first introducd and dvlopd in [10] 31 Absolut stability of LDG mtod #1 W first rcast t LDG mtod #1 as t following nonconforming mixd mtod, wic is asily obtaind by adding (29 and (210 Find (u, σ V Σ suc tat (31 A (u, σ ; v, τ = F (v, τ (v, τ V Σ, wr (32 A (w, χ ; v, τ = (χ, v Ω k 2 (w, v Ω + ık w, v Γ { w } ıβ[[w ]], [[v ]] E (ıδ [[ w ]], [[τ ]] [[w ]], {τ } (33 F (v, τ := (f, v Ω + g, v Γ E + (χ, τ Ω ( w, τ Ω 311 A gnralizd inf-sup condition T goal of tis subsction is to sow tat t ssquilinar form A dfind in (32 satisfis a gnralizd inf-sup condition, wic will play a vital rol for us to stablis t absolut stability of t LDG mtod #1 in t nxt subsction Proposition 31 Tr xists an - and k-indpndnt constant c 1 > 0 suc tat for any (w, χ V Σ (34 wr (35 (36 sup (v,τ V Σ v 0 γ 1 := 1 + k + R A (w, χ ; v, τ v DG + sup (v,τ V Σ v 0 β ( k 2 δ + max E β β 3, w DG := ( w 1, := w 2 L 2 (K K T m A (w, χ ; v, τ v DG c 1 γ 1 w DG, ( k 2 w 2 L 2 (Ω + k2 w 2 L 2 (Γ + c Ω w 2 L 2 (Γ + w 2 1, 1 2 Proof T main ida of t proof is tat for a fixd (w, χ V Σ w nd to pick up two sts of spcial functions (v, τ V Σ, wic, as xpctd, must dpnd on (w, χ V Σ, suc tat bot quotints in (34 can b boundd 1 2,

8 8 XAOBNG FENG AND YULONG XNG from blow by w DG Wn tat is don, t inf-sup constant c 1 /γ 1 will b rvald in t procss Sinc t proof is vry long, w divid it into four stps Stp 1: Taking t first tst function W first coos t tst function (v, τ = (w, w to gt A (w, χ ; v, τ = A (w, χ ; w, w = ( w, w Ω k 2 (w, w Ω + ık w, w Γ + (ıδ [[ w ]], [[ w ]] [[w ]], { w } E E Taking t ral and imaginary parts yilds (37 (38 R A (w, χ ; w, w { w } ıβ[[w ]], [[w ]] = w 2 1, k 2 w 2 L 2 (Ω 2 R [[w ]], { w }, E m A (w, χ ; w, w = ( δ [[ w ]] 2 L 2 ( + β [[w ]] 2 L 2 ( + k w 2 L 2 (Γ E Stp 2: Taking t scond tst function nspird by t spcial tst function tcniqu of [10], w now coos anotr tst function (v, τ = (α w, w wit α := x x Ω (s Dfinition 21 and us t fact tat v = w to gt A (w, χ ; v, τ = A (w, χ ; α w, w (39 = ( w, w Ω k 2 (w, α w Ω + ık w, α w Γ + (ıδ [[ w ]], [[ w ]] [[w ]], { w } E E { w } ıβ[[w ]], [[α w ]] Taking t ral part immdiatly givs (not tat v = α w R A (w, χ ; v, w = w 2 1, k 2 (w, v Ω k m w, v Γ β m [[w ]], [[v ]] E ( R [[w ]], { w } + { w }, [[v ]] E Stp 3: Driving an uppr bound for k 2 w 2 L 2 (Ω To bound (37 from blow, w nd to gt an uppr bound for t trm k 2 w 2 L 2 (Ω on t rigt and sid of (37 Tis will b don by carfully and judicially combining (38 and (39 wit som otr diffrntial idntitis, wic w now xplain

9 LDG METHODS FOR THE HELMHOLTZ EQUATON 9 Using t intgral idntity (310 2 w 2 L 2 (K = α n K w 2 ds 2 R(w, v K (d 2 w 2 L 2 (K, K and (37, (39 w gt (311 2k 2 w 2 L 2 (Ω = 2 R A (w, χ, v, w + (d 2 R A (w, χ ; w, w + 2k m w, v Γ 2 w 2 1, (d 2 w 2 1, + 2 R ( [[w ]], { w } + { w }, [[v ]] E + 2(d 2 R E + E [[w ]], { w } (2β m [[w ]], [[v ]] + k 2 K T K α n K w 2 ds By t lmntary idntity a 2 b 2 = R(a + b(ā b for any two complx numbrs a and b w av (312 K T K α n K w 2 ds = 2 R E α n {w }, [w ] + α n Ω, w 2 Γ Nxt, using t local Rllic idntity (s [10, Lmma 41] (d 2 w 2 L 2 (K + 2 R( w, v K = α n K v 2 ds K and t fact tat v = w, w gt (313 d w 2 1, = (d 2 w 2 1, + 2 R = R K T = 2 R E K α n K v 2 ds K T ( w, v K α n { w }, [ w ] + α n, w 2 E B

10 10 XAOBNG FENG AND YULONG XNG Substituting (312 and (313 into (311 w gt (314 2k 2 w 2 L 2 (Ω = 2 R A (w, χ ; v, w + (d 2 R A (w, χ ; w, w + 2k 2 R α n {w }, [w ] + k 2 α n Ω, w 2 Γ E + 2k m w, v Γ E B α n, w 2 ( α n { w }, [ w ] { w }, [[v ]] 2 R E + 2(d 1 R E [[w ]], { w } + 2 m E β [[w ]], [[v ]] To gt an uppr bound for k 2 w 2 L 2 (Ω, w nd to bound t trms on t rigt-and sid of (314, wic w now bound as follows (315 (316 2k 2 R E α n {w }, [w ] Ck 2 E k 2 α n Ω, w 2 Γ Ck2 w 2 L 2 (Γ 1 2 w L2 (K K [w ] L2 ( k2 2 w 2 L 2 (Ω + C k 2 β [w ] 2 L β 2 ( E t follows from t star-sapd assumption on Ω tat (317 2k m w, v Γ α n, w 2 E B Ck w L 2 ( w L 2 ( c Ω w 2 L 2 ( E B Ck 2 w 2 L 2 (Γ c Ω 2 w 2 L 2 (Γ E B By t trac inquality [2], w also av (318 2d R E [[w ]], { w } 2d E 1 2 K=K,K 1 4 w 2 1, + C E w L2 (K [w ] L2 ( 1 β β [w ] 2 L 2 (

11 (319 2 R E LDG METHODS FOR THE HELMHOLTZ EQUATON 11 ( α n { w }, [ w ] { w }, [[v ]] = 2 R E d 1 ( (α τ j { w n } j=1 ] (α n { w τ j } [w ] τ j E d 1 j=1 1 2 K=K,K 1 4 w 2 1, + C E w L 2 (K [ w τ j ] L 2 ( 1 d 1 β [ w τ j β ] 2 L 2 ( j=1 By t dfinition of v := α w, w gt (320 2 m E = 2 m E β [[w ]], [[v ]] = 2 m E [ β [w ], β [w ], [v ] d 1 ] (α n w n + (α τ j w τ j j=1 C E C + C E β δ β [w ] L 2 ( [ w n ] L 2 ( d 1 β [w ] L 2 ( [ w τ j ] L 2 ( E j=1 ( β [w ] 2 L 2 ( + δ [ w n ] 2 L 2 ( + Cβ ( [w ] 2 L 2 ( + d 1 [ w τ j ] 2 L 2 (, E j=1 wr {τ j } d 1 j=1 dnots an ortogonal tangntial fram on t dg/fac, and w av usd t dcomposition α = (α n n + d 1 j=1 (α τ j τ j

12 12 XAOBNG FENG AND YULONG XNG Now substituting stimats (315 (319 into (314 w obtain (321 2k 2 w 2 L 2 (Ω 2 R A (w, χ ; v, w + (d 2 R A (w, χ ; w, w + Ck 2 w 2 L 2 (Γ c Ω 2 w 2 L 2 (Γ + k2 2 w 2 L 2 (Ω + C E k 2 β [w ] 2 L β 2 ( 2 R [[w ]], { w } E w 2 1, + C 1 d 1 β [ τ j β ] 2 L 2 ( E j= w 2 1, + C 1 β [w ] 2 L β 2 ( E + C E β ( δ β [w ] 2 L 2 ( + δ [ w n ] 2 L 2 ( + C (β [w ] 2 L 2 ( + d 1 β [ w τ j ] 2 L 2 ( E j=1 On noting tat (38 provids uppr bounds for trms [ w n ] 2 L 2 (, [w ] 2 L 2 ( and k 2 w 2 L 2 (Γ in trms of m A (w, χ ; w, w, using ts bounds in (321 w gt (322 2k 2 w 2 L 2 (Ω + k2 2 w 2 L 2 (Γ + c Ω 2 w 2 L 2 (Γ 2 R A (w, χ ; v, w + (d 2 R A (w, χ ; w, w w 2 1, + k2 2 w 2 L 2 (Ω 2 R { w }, [[w ]] E d 1 + C ( β [ w τ j β ] 2 L 2 ( E j=1 + M 1 m A (w, χ ; w, w, wr (323 ( β M 1 := C 1 + k + δ + max k E β To bound t jumps of t tangntial drivativs [ w τ j ] 2 L 2 ( w appal to t invrs inquality in (321, (324 [ w τ j ] 2 L 2 ( C 2 [w ] 2 L 2 (,

13 LDG METHODS FOR THE HELMHOLTZ EQUATON 13 and using (38 and (37 to gt 2k 2 w 2 L 2 (Ω + k2 2 w 2 L 2 (Γ + c Ω 2 w 2 L 2 (Γ 2 R A (w, χ ; v, w + (d 2 R A (w, χ ; w, w + k2 2 w 2 L 2 (Ω + M 2 m A (w, χ ; w, w 2 R { w }, [[w ]] w 2 1, E 2 R A (w, χ ; v, w + (d 1 R A (w, χ ; w, w + 3k2 2 w 2 L 2 (Ω + M 2 m A (w, χ ; w, w 1 2 w 2 1,, wr (325 C M 2 = M 1 + max E( 2 + C β 3 ( β ( k 2 = C 1 + k + δ + max + 1 E β β 3 Using t linarity of t ssquilinar form A w av (326 k 2 w 2 L 2 (Ω + k2 w 2 L 2 (Γ + c Ω w 2 L 2 (Γ + w 2 1, 4 R A (w, χ ; v, w + 2(d 1 R A (w, χ ; w, w + 2M 2 m A (w, χ ; w, w = R A (w, χ ; w, w + 2M 2 m A (w, χ ; w, w wit w = 4v + 2(d 1w Stp 4: Finising up By t dfinition of DG in (36 and t fact tat v = w w av (327 v 2 DG = k 2 α w 2 L 2 (Ω + k2 α w 2 L 2 (Γ + c Ω w 2 L 2 (Γ + w 2 1, Ck 2 w 2 L 2 (Ω + Ck2 w 2 L 2 (Γ ( + c Ω w 2 L 2 (Γ + w 2 1, C(1 + k 2 w 2 1, + c Ω w 2 L 2 (Γ C(1 + k 2 w 2 DG Tus, it follows from t triangl inquality tat (328 w DG 4 v DG + 2(d 1 w DG C(1 + k w DG

14 14 XAOBNG FENG AND YULONG XNG Now from (326 and (328 w av (329 R A (w, χ ; w, w + m A (w, χ ; w, w w DG w DG R A (w, χ ; w, w + m A (w, χ ; w, w C(1 + k w DG w DG 1 R A (w, χ ; w, w + 2M 2 m A (w, χ ; w, w 2M 2 w DG c 1 γ 1 w DG for som constant c 1 > 0 and γ 1 is dfind by (35 Hnc, (34 olds T proof is complt Rmark 31 (a W not tat γ 1 dpnds on bot and k (b T gnralizd inf-sup condition is a wak stimat bcaus it dos not provid a control for t variabl χ As a comparison, w rcall tat t standard inf-sup condition for t ssquilinar form A sould b sup (v,τ V Σ A (w, χ ; v, τ (v, τ c 1 (w, χ for som positiv constant c 1 = c 1 (k, β, δ, Ω Wr (w, χ := ( k 2 w 2 L 2 (Ω + χ L 2 (Ω (w, χ V Σ Howvr, t abov standard inf-sup condition can b provd only undr t ms constraint = O(k 2 and w bliv tat it dos not old witout a ms constraint 312 Stability stimats T goal of tis subsction is to stablis t absolut stability for t LDG mtod #1 using t gnralizd inf-sup condition provd in t prvious subsction Torm 31 Lt (u, σ V Σ solv (31 Dfin (330 M(f, g := f Ω + g L2 (Γ Tn tr old t following stability stimats: (331 (332 u DG γ 1 k 1 M(f, g σ L 2 (Ω γ 1 k 1( 1 + (δ + k 1 ( max 1 K M(f, g K T Proof By Scwarz inquality w av (333 F (v, χ f L 2 (Ω v L 2 (Ω + g L 2 (Γ v L 2 (Γ Ck 1 M(f, g ( k 2 v 2 L 2 (Ω + k2 v 2 L 2 (Γ Ck 1 M(f, g v DG 1 2

15 LDG METHODS FOR THE HELMHOLTZ EQUATON 15 Lt (w, χ = (u, σ in (34 By quation (31 and (333 w gt c 1 R A (u, σ ; v, τ u DG sup γ 1 (v,τ V Σ v DG m A (u, σ ; v, τ + sup (v,τ V Σ v DG R F (v, τ m F (v, τ = sup + sup (v,τ V Σ v DG (v,τ V Σ v DG F (v, τ 2 sup v DG (v,τ V Σ Ck 1 M(f, g Hnc (331 olds To sow (332, stting (v, τ = (0, σ in (31 and using t trac and Scwarz inqualitis yilds σ 2 L 2 (Ω = ( u, σ Ω + E u 2 L 2 (Ω σ 2 L 2 (Ω + C E (ıδ [[ u ]], [[σ ]] [[u ]], {σ } K=K,K + u L2 (K ( 1 K δ u L 2 (K σ L2 (K Tus, u 2 L 2 (Ω σ 2 L 2 (Ω + C E K=K,K 2 K ( δ 2 u 2 L 2 (K + u 2 L 2 (K σ 2 L 2 (Ω (1 + δ 2( max 2 K u 2 1, + ( max 2 K u 2 L K T K T 2 (Ω T dsird stimat (332 follows from combining t abov inquality wit (331 T proof is complt An immdiat consqunc of t stability stimats is t following uniqu solvability torm Torm 32 Tr xists a uniqu solution to t LDG mtod (31 for all k,, δ, β > 0 Proof Sinc problm (31 is quivalnt to a linar systm, nc, it suffics to sow t uniqunss But t uniqunss follows immdiatly from t stability stimats as t zro sourcs imply tat any solution must b a trivial solution 32 Absolut stability of LDG mtod #2 n tis subsction, w considr t LDG mtod #2 By adding (27 and (28 wit t givn numrical fluxs, w tn rcast t LDG mtod #2 as t following nonconforming mixd mtod: Find (u, σ V Σ suc tat (334 B (u, σ ; v, τ = F (v, τ (v, τ V Σ,

16 16 XAOBNG FENG AND YULONG XNG wr F is dfind in (33 and (335 B (w, χ ; v, τ = (χ, v Ω k 2 (w, v Ω + ık w, v Γ {χ } ıβ[[w ]], [[v ]] E (ıδ [[χ ]], [[τ ]] [[w ]], {τ } E + (χ, τ Ω ( w, τ Ω 321 A gnralizd inf-sup condition T goal of tis subsction is to sow tat t ssquilinar form B dfind in (335 for t LDG mtod #2 satisfis anotr gnralizd inf-sup condition To t nd, w introduc t following spac notation: (336 S := {(w, χ V Σ ; (w, χ satisfis (337}, wr (337 (χ, τ Ω ( w, τ Ω (ıδ [[χ ]], [[τ ]] [[w ]], {τ } = 0 (v, τ V Σ E W not tat it is asy to cck tat (w, χ S implis tat it satisfis (27 wit û K bing dfind by t LDG mtod #2 Lmma 31 For any (w, χ S, tr olds t following stimats: (338 ( w 1, 16 χ L 2 (Ω ( + C E ( χ w L 2 (Ω C ( 1 [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 ( E 1 2, ( 1 [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 ( 1 2

17 LDG METHODS FOR THE HELMHOLTZ EQUATON 17 Proof On noting tat (w, χ satisfis (337, stting τ = w in (337, w gt w 2 1, = R(χ, w Ω R (ıδ [[χ ]], [[ w ]] [[w ]], { w } E 1 2 w 2 1, χ 2 L 2 (Ω + δ 1 2 E + E 1 2 K=K,K K=K,K 1 2 w 2 1, χ 2 L 2 (Ω w 2 1, + C E [[χ ]] L2 ( w L2 (K [[w ]] L 2 ( w L 2 (K ( 1 [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 ( Trfor, w 2 1, χ 2 L 2 (Ω + C ( 1 [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 (, E wic givs (338 T stimat (339 follows from t sam drivation by stting τ = χ w in (337 T proof is complt W now ar rady to stat a gnralizd inf-sup condition for t ssquilinar form B Proposition tr olds (340 wr (341 ( Tr xits constant c 2 > 0 suc tat for any (w, χ S R B (w, χ sup ; v, τ (v,τ V Σ (v, τ DG m B (w, χ + sup ; v, τ c 2 (w, χ (v,τ V Σ (v, τ DG γ DG, 2 ( k γ 2 := k + max + β + δ + δ E β β 3, ( (w, χ DG := k 2 w 2 L 2 (Ω + k2 w 2 L 2 (Γ + χ 2 L 2 (Ω + c Ω w 2 L 2 (Γ Proof Sinc t proof follows t sam lins as tat of Proposition 31, w sall only igligt t main stps and point out t diffrncs Stp 1: Taking t first tst function 1 2

18 18 XAOBNG FENG AND YULONG XNG W first coos t tst function (v, τ = (w, χ to gt B (w, χ ; v, τ = B (w, χ ; w, χ = (χ, χ Ω k 2 (w, w Ω + ık w, w Γ {χ } ıβ[[w ]], [[w ]] E (ıδ [[χ ]], [[χ ]] [[w ]], {χ } Taking t ral and imaginary parts yilds E + (χ, w Ω ( w, χ Ω (343 R B (w, χ ; w, χ = χ 2 L 2 (Ω k2 w 2 L 2 (Ω, (344 m B (w, χ ; w, χ = k w 2 L 2 (Γ + ( δ [[χ ]] 2 L 2 ( + β [[w ]] 2 L 2 ( E + 2 m(χ, w Ω + E 2 m [[w ]], {χ } On noting tat (w, χ satisfis (337, stting τ = χ in (337, w gt χ 2 L 2 (Ω ( w, χ Ω = (ıδ [[χ ]] 2 L 2 ( [[w ]], {χ }, E Taking t imaginary part yilds m(χ, w Ω + E Hnc, (344 bcoms (345 m [[w ]], {χ } = E δ [[χ ]] 2 L 2 ( m B (w, χ ; w, χ = k w 2 L 2 (Γ + ( δ [[χ ]] 2 L 2 ( + β [[w ]] 2 L 2 ( E Stp 2: Taking t scond tst function Nxt, w coos anotr tst function (v, τ = (α w, χ, wic is diffrnt from t on usd in t proof of Proposition 31, and us t fact tat v = w to gt B (w, χ ; v, τ = B (w, χ, α w, χ = (χ, χ Ω k 2 (w, α w Ω + ık w, α w Γ + E ıβ [[w ]], [[α w ]] E ıδ [[χ ]], [[χ ]] + E E + 2ı m(χ, w {χ }, [[α w ]] [[w ]], {χ }

19 LDG METHODS FOR THE HELMHOLTZ EQUATON 19 Taking t ral part immdiatly givs (v = α w (346 R B (w, χ ; v, χ = χ 2 L 2 (Ω k2 R(w, v Ω k m w, v Γ + (R {χ }, [[w ]] [[v ]] β m [[w ]], [[v ]] E Stp 3: Driving an uppr bound for k 2 w 2 L 2 (Ω To bound (343 from blow, w again nd to gt an uppr bound for t trm k 2 w 2 L 2 (Ω on t rigt and sid of (343 Using t intgral idntity (310, (312, (343 and (346, w av (347 2k 2 w 2 L 2 (Ω = 2 R B (w, χ ; v, χ + (d 2 R B (w, χ ; w, χ d χ 2 L 2 (Ω + 2k m w, v Γ + k 2 α n Ω, w 2 Γ + 2k 2 R E α n {w }, [w ] + 2 R E + 2(d 1 R {χ }, [[w ]] + 2 m E E + d w 2 1, E B α n, w 2 {χ }, [[v ]] β [[w ]], [[v ]] 2 R E ± 2 R E α n { w }, [ w ] α n {χ }, [ w ] W not tat by (313 t sum of t scond and tird lins to t last is zro, and t contribution of t last lin is obviously zro Ts trms ar purposly addd in ordr to gt sarpr uppr bounds wn ty ar combind wit t trms prcding tm W now nd to bound t trms on t rigt-and sid of (347 Som of ts av bn obtaind in t proof of t Proposition 31, and t otrs ar drivd as follows (348 2(d 1 R E {χ }, [[w ]] 2(d 1 E 1 2 K=K,K 1 16 χ 2 L 2 (Ω + C E χ L2 (K [w ] L2 ( 1 β β [w ] 2 L 2 (

20 20 XAOBNG FENG AND YULONG XNG (349 2 R {χ }, [[v ]] 2 R α n {χ }, [ w ] E E = 2 R d 1 ( (α τ j {χ n } E j=1 ] (α n {χ τ j } [w ] τ j E d 1 j=1 1 2 K=K,K 1 16 χ 2 L 2 (Ω + C E χ L2 (K [ w τ j ] L2 ( 1 d 1 β [ w τ j β ] 2 L 2 ( j=1 t follows from (338 and (339 tat (350 2 R E E α n {χ w }, [ w ] 1 C E K=K,K χ w L2 (K w L2 (K C max 2 E χ w 2 L 2 (Ω w L 2 (Ω ( 1 3 [[w ]] 2 L 2 ( + δ2 3 [[χ ]] 2 L 2 ( χ 2 L 2 (Ω + C E 1 16 χ 2 L 2 (Ω + C 1 E( 3 ( 1 [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 ( [[w ]] 2 L 2 ( + δ2 3 [[χ ]] 2 L 2 ( Similar to t drivation of (320, w av (351 2 m E C E C E β [[w ]], [[v ]] = 2 m E β 1 2 K=K,K β [w ], [α w ] [w ] L 2 ( w L 2 (K β 2 1 [w ] 2 L 2 ( w 2 L 2 (Ω 1 16 χ 2 L 2 (Ω + C E ( β [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 (

21 LDG METHODS FOR THE HELMHOLTZ EQUATON 21 Now substituting stimats (348 (351, togtr wit (315, (316 and (317, into (347 w obtain (352 2k 2 w 2 L 2 (Ω 2 R B (w, χ ; v, χ + (d 2 R B (w, χ ; w, χ d χ 2 L 2 (Ω + Ck 2 w 2 L 2 (Γ c Ω w 2 L 2 2 ( + k2 2 w 2 L 2 (Ω + C E E B k 2 β [w ] 2 L β 2 ( χ 2 L 2 (Ω + C E ( β d χ 2 L 2 (Ω + Cd E χ 2 L 2 (Ω + C E 1 β β [w ] 2 L 2 ( [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 ( 1 d 1 β [ w τ j β ] 2 L 2 ( j= χ 2 L 2 (Ω + C 1 E( 3 [[w ]] 2 L 2 ( + δ2 3 [[χ ]] 2 L 2 ( χ 2 L 2 (Ω + C E ( β [[w ]] 2 L 2 ( + δ2 [[χ ]] 2 L 2 ( On noting tat (345 provids uppr bounds for trms [[χ ]] 2 L 2 (, [w ] 2 L 2 ( and k 2 w 2 L 2 (Γ in trms of m B (w, χ ; w, χ, and t jumps of t tangntial drivativs [ w τ j ] 2 L 2 ( in (352 can b boundd by t invrs inquality (324, w gt (353 wr (354 2k 2 w 2 L 2 (Ω + k2 2 w 2 L 2 (Γ + c Ω 2 w 2 L 2 ( E B 2 R B (w, χ ; v, χ + (d 2 R B (w, χ ; w, χ ( + M 3 m B (w, χ ; w, χ + k2 1 2 w 2 L 2 (Ω d χ 16 2 L 2 (Ω 2 R B (w, χ ; v, χ + (d 1 R B (w, χ ; w, χ + M 3 m B (w, χ ; w, χ + 3k2 2 w 2 L 2 (Ω 1 2 χ 2 L 2 (Ω, ( ( k M 3 = C k + max E β + β + δ + δ 3 Using t linarity of t ssquilinar form B, w av + 1 β 3 k 2 w 2 L 2 (Ω + k2 w 2 L 2 (Γ + χ 2 L 2 (Ω + c Ω w 2 L 2 (Γ 4 R B (w, χ ; v, χ + 2(d 1 R B (w, χ ; w, χ + 2M 3 m B (w, χ ; w, χ = R B (w, χ ; w, χ + 2M 3 m B (w, χ ; w, χ,

22 22 XAOBNG FENG AND YULONG XNG wr w = 4v + 2(d 1w and v = α w Stp 4: Finising up t follows from t invrs inquality (324 tat k 2 v 2 L 2 (Ω + k2 v 2 L 2 (Γ Ck2( max 2 K w 2 L K T 2 (Ω + c Ω w 2 L 2 (Γ Tn by t dfinition of t norm (w, χ DG in (342 and t fact tat v = w, w gt Trfor, (355 ( w, χ DG 4 (v, χ DG + 2(d 1 (w, χ DG C ( max 1 K + k + 1 (w, χ K T DG R B (w, χ ; w, χ + m B (w, χ ; w, χ ( w, χ DG (w, χ DG R B (w, χ ; w, χ C ( max K T 1 K + k m B (w, χ ; w, χ (w, χ DG (w, χ DG 1 R B (w, χ ; w, χ + M 3 m B (w, χ ; w, χ 2M 3 (w, χ DG c 2 γ 2 (w, χ DG, for som constant c 2 > 0 and γ 2 dfind by (341 Hnc, (340 olds and t proof is complt 322 Stability stimats T gnralizd inf-sup condition provd in t last subsction immdiatly infrs t following (absolut stability and wll-posdnss torms for t LDG mtod #2 Torm 33 Lt (u, σ V Σ solv (334 Tn tr olds (356 (u, σ DG γ 2 k 1 M(f, g Proof On noting tat any solution (u, σ of (334 blongs to t st S, t dsird stimat (33 follows radily from (340 wit (w, χ = (u, σ, (334 and (333 T proof is complt Torm 34 T LDG mtod (334 as a uniqu solution for all k,, δ, β > 0 Sinc t proof of t abov torm is a vrbatim copy of tat of Torm 32, w omit it 4 Error stimats T goal of tis sction is to driv rror stimats for t LDG mtod #1 and #2 Following t ida of [10], tis will b don in two stps First, w introduc an lliptic projction of t solution (u, σ using a corrsponding corciv ssquilinar form of A (rsp B and driv rror bounds for t projction W not tat t rror analysis for t lliptic projctions as an indpndnt intrst in itslf (cf [3] Scond, w bound t rror btwn t projction and t LDG solution using t stability stimats obtaind in Sction 3 Sinc t rror analysis for t two LDG mtods ar similar, w sall giv mor dtails of t rror analysis for

23 LDG METHODS FOR THE HELMHOLTZ EQUATON 23 t LDG mtod #1 but sall b brif for t LDG mtod #2 Trougout tis sction, w lt for j = 1, 2 H j (T = H j (K, = max K max, β = β 0 1, δ = δ 0 K T E K T for som positiv constants β 0 and δ 0 41 Error stimats for t LDG mtod #1 411 Elliptic projction and its rror stimats For any (w, χ H 2 (T H 1 (T d, w dfin t lliptic projction ( w, χ V Σ of (w, χ by (41 wr (42 a ( w, χ ; v, τ = a (w, χ; v, τ (v, τ V Σ, a (w, χ ; v, τ : = A (w, χ ; v, τ + k 2 (w, v Ω = (χ, v Ω + ık w, v Γ { w } ıβ 0 1 [[w ]], [[v ]] E (ıδ [[ w ]], [[τ ]] [[w ]], {τ } E + (χ, τ Ω ( w, τ Ω To driv rror bounds for t abov lliptic projction, w first notic tat a (w, χ ; v, v = ( w, v Ω + ık w, v Γ { w } ıβ 0 1 [[w ]], [[v ]] As a rsult, w V satisfis (43 Morovr, sinc E + (ıδ [[ w ]], [[ v ]] [[w ]], { v } E =: A (w, v A ( w, v = A (w, v v V a (w, χ ; 0, τ = (χ, τ Ω ( w, τ Ω (ıδ [[ w ]], [[τ ]] [[w ]], {τ }, E w av tat χ Σ satisfis (44 ( χ, τ Ω = ( w, τ Ω + (ıδ [[ w w]], [[τ ]] [[ w w]], {τ } E + (χ, τ Ω ( w, τ Ω τ Σ

24 24 XAOBNG FENG AND YULONG XNG Lmma 41 For any w, v H 2 (T, tr xists a k- and -indpndnt constant C suc tat (45 A (w, v C w 1, v 1, Morovr, for any ɛ (0, 1, tr xists a constant c ɛ > 0 suc tat (46 R A (v, v + (1 ɛ + c ɛ m A (v, v (1 ɛ v 2 1,, wr (47 (48 w 1, := w 1, := ( w 2 L 2 (Ω + k w 2 L 2 (Γ + E ( β [w] 2 L 2 ( + δ [[ w]] 2 L 2 ( ( w 2 1, + β 1 { w n } 2 L 2 ( E , Sinc t proof of t abov lmma is lmntary, w omit it W now rcall t following stability stimat for u (cf [7, 10]: u H2 (Ω (k 1 + km(f, g, wic is ndd to prov t nxt lmma and will b usd svral tims in t rst of tis sction Proposition 41 Lt u H 2 (Ω b t solution to problm (11-(12 and σ = u Lt (ũ, σ V Σ dnot t lliptic projction of (u, σ dfind by (41 Tn tr old t following rror stimats: (49 u ũ 1, + k 1 2 u ũ L2 (Γ (1 + k 1 2 k, (410 u ũ L2(Ω (1 + kk 2, (411 σ σ L2 (Ω (1 + k 1 2 k Proof Sinc t proof of (49 and (410 is ssntially sam as tat of [10, Lmma 52], w omit it to sav t spac and rfr t radr to [10] for t dtails To sow (411, on noting tat (44 and t idntity (σ, τ = ( u, τ imply (412 (σ σ, τ Ω = ( u ũ, τ Ω E [[u ũ ]], {τ } ( ıδ 0 [[ (u ũ ]][[τ ]] τ Σ

25 LDG METHODS FOR THE HELMHOLTZ EQUATON 25 For any χ Σ, w st τ = χ σ Tn by (412, t trac inquality, Scwarz inquality w gt σ σ 2 L 2 (Ω = (σ σ, σ χ Ω + (σ σ, τ Ω = (σ σ, σ χ Ω + ( (u ũ, τ Ω (ıδ 0 [[ (u ũ ]], [[τ ]] [[u ũ ]], {τ } E σ σ L2 (Ω σ χ L2 (Ω + (u ũ L2 (Ω τ L2 (Ω ( + C δ 0 [[ (u ũ ]] 2 L 2 ( + β 0 1 [u ũ ] 2 L 2 ( E + ε E τ 2 L 2 ( 1 4 σ σ 2 L 2 (Ω τ 2 L 2 (Ω + σ χ 2 L 2 (Ω + C u ũ 2 1, 1 2 σ σ 2 L 2 (Ω σ χ 2 L 2 (Ω + C u ũ 2 1, Hnc, it follows from t abov inquality, (49, and t polynomial approximation tory (cf [2] tat σ σ L 2 (Ω C u ũ 1, + 2 inf σ χ L χ Σ 2 (Ω (1 + k 1 2 k + (k + k 1 (1 + k 1 2 k, wic givs (411 T proof is complt 412 Global rror stimats for t LDG mtod #1 n t prcding subsction w av drivd t rror bounds for (u ũ, σ σ By t dcomposition u u = (u ũ + (ũ u and σ σ = (σ σ + ( σ σ and t triangl inquality, it suffics to gt rror bounds for (ũ u, σ σ W sall accomplis tis task by xploiting t linarity of t Hlmoltz quation and using t stability stimat for t LDG mtod #1 obtaind in Sction 31 First, on noting tat (u, σ satisfis (413 A (u, σ; v, τ = F (v, τ (v, τ V Σ Subtracting (31 from (413 yilds t following rror quation (or Galrkin ortogonality: (414 A (u u, σ σ ; v, τ = 0 (v, τ V Σ Nxt, to procd w introduc t notation u u = + q, := u ũ, q := ũ u, σ σ = ψ + φ, ψ := σ σ, φ := σ σ

26 26 XAOBNG FENG AND YULONG XNG Tn by (414 and t dfinitions of t ssquilinar form a and t lliptic projction w av (415 A (q, φ ; v, τ = A (, ψ ; v, τ = a (, ψ ; v, τ + k 2 (, v = k 2 (, v, (v, τ V Σ T abov quation implis tat (q, φ V Σ is t LDG solution to t Hlmoltz problm wit sourc trms f = k 2 and g = 0 Tn an application of t stability stimats of Torm 31 immdiatly yilds t following lmma Proposition 42 Tr old t following stimats for (q, φ : (416 (417 q DG γ 1 (1 + kk 2 2, φ L 2 (Ω γ 1 (1 + k + δ 0 k 2 (1 + kk Combining Proposition 41 and 42, using t triangl inquality and t standard duality argumnt giv t following rror stimats for (u, σ Torm 41 Lt u H 2 (Ω b t solution to problm (11 (12 and σ := u, and (u, σ b t solution to problm (31 Tn tr old t following rror stimats for (u, σ : (418 (419 (420 u u 1, + k 1 2 u u L2 (Γ ( (1 + k γ1 (1 + kk k, u u L 2 (Ω (1 + γ 1 (1 + kk 2, σ σ L2 (Ω ( (1 + k γ1 (1 + k(1 + k + δ 0 k 2 k 42 Error stimats for t LDG mtod #2 T rror analysis for t LDG mtod #2 ssntially follows t sam lins as tat for t LDG mtod #1 givn in t prvious subsction Howvr, tr ar tr main diffrncs wic w now xplain First, t ssquilinar form a nds to b rplacd by anotr ssquilinar form b in t dfinition of t lliptic projction (41, wr b is dfind by (421 b (w, χ ; v, τ : = B (w, χ ; v, τ + k 2 (w, v Ω = (χ, v Ω + ık w, v Γ {χ } ıβ[[w ]], [[v ]] E (ıδ [[χ ]], [[τ ]] [[w ]], {τ } E + (χ, τ Ω ( w, τ Ω Scond, du to strong coupling btwn ũ and σ, t rror stimats for t nw lliptic projction (ũ, σ must b drivd diffrntly To t nd, w nd t following lmma, wic rplacs Lmma 41 Lmma 42 Lt β = β 0 1 and δ = δ 0 for som positiv constants β 0 and δ 0 (i Tr xists an - and k-indpndnt constant c 3 > 0 suc tat t ssquilinar form b satisfis t following gnralizd inf-sup condition: for any fixd

27 LDG METHODS FOR THE HELMHOLTZ EQUATON 27 (w, χ V Σ (422 R b (w, χ sup ; v, τ (v,τ V Σ (v, τ DG + sup (v,τ V Σ m b (w, χ ; v, τ (v, τ DG c 3 (w, χ DG (ii Tr xists an - and k-indpndnt constant C > 0 suc tat for any (w, χ, (v, τ H 2 (T H 1 (T d, tr olds (423 wr (424 (425 b (w, χ; v, τ C (w, χ 1, (v, τ 1,, (w, χ DG := (w, χ 1, := ( w 2 1, + χ 2 L 2 (Ω 1 2, ( (w, χ 2 DG + β 1 {χ} 2 L 2 ( E T proof of (i is basd on valuating t first quotint on t lft-and sid of (422 at (v, τ = ( (1+C 1 w, C 1 χ w and valuating t scond quotint at (v, τ = (C 2 w, C 2 χ for som sufficintly larg positiv constants C 1 and C 2 T proof of (ii is a straigtforward application of Scwarz and trac inqualitis W skip t rst of t drivation to sav spac T abov gnralizd inf-sup condition, t bounddnss of t ssquilinar form b, and t duality argumnt (cf [2] radily infr t following rror stimats for t nw lliptic projction (ũ, σ W omit t proof sinc it is standard Proposition 43 Undr t assumptions of Proposition 41, tr old t following stimats: (426 (427 u ũ 1, + σ σ L2 (Ω k, u ũ L 2 (Ω k T tird diffrnc is tat t nw rror function (q, φ now satisfis (428 B (q, φ ; v, τ = k 2 (, v (v, τ V Σ As a rsult, by Torm 33 and (338 w gt (429 q 1, + (q, φ DG γ 2 (1 + kk 2 2, wic rplacs stimats (416 and (417 Aftr aving stablisd Proposition 43 and (429, onc again, by t triangl inquality w arriv at t following rror stimats for t solution (u, σ to t LDG mtod #2 Torm 42 Lt u H 2 (Ω b t solution to problm (11 (12 and σ := u, and (u, σ b t solution to problm (334 Tn tr old t following rror stimats for (u, σ : (430 (431 u u 1, + k 1 2 u u L 2 (Γ + σ σ L 2 (Ω ( 1 + γ 2 (1 + kk k, u u L2 (Ω (1 + γ 2 (1 + kk 2 2

28 28 XAOBNG FENG AND YULONG XNG Rmark 41 (429 sows tat φ := σ σ as an optimal ordr (in rror bound for t LDG mtod #2, wil (417 sows tat φ only as a sub-optimal ordr rror bound for t LDG mtod #1 W bliv tat tis is t main rason wy in practic t LDG mtod #2 givs a bttr approximation to t flux variabl σ tan t LDG mtod #1 dos altoug bot mtods av t sam asymptotic rat of convrgnc in 5 Numrical xprimnts n tis sction w sall provid som numrical rsults of t two proposd LDG mtods Our tsts ar don for t following 2-d Hlmoltz problm: (51 u k 2 u = f := sin(kr r in Ω, (52 u + ıku = g n Ω on Γ R := Ω Hr Ω is t unit squar [ 05, 05] [ 05, 05], and g is cosn so tat t xact solution is givn by (53 u = cos(kr k cos k + ı sin k k ( J 0 (k + ıj 1 (k J 0(kr in polar coordinats, wr J ν (z ar Bssl functions of t first kind Assum T 1/m b t rgular triangulation tat consists of 2m 2 rigt-angld quicrural triangls of siz = 1/m, for any positiv intgr m S Figur 1 for t sampl triangulation T 1/4 and T 1/10 Figur 1 T computational domain and sampl mss Lft: T 1/4 tat consists of rigt-angld quicrural triangls of siz = 1 4 ; Rigt: T 1/10 wit = Snsitivity wit rspct to t paramtrs δ and β n tis subsction, w xamin t snsitivity of t rror of t LDG solutions in H 1 -sminorm wit rspct to t paramtrs δ and β T LDG mtod #1 is considrd first W start by fixing δ = 01 and tsting t snsitivity in t paramtr β Wit two wav numbrs k = 5 and 50, w comput t solutions of t LDG mtod #1 wit diffrnt valus of β:

29 LDG METHODS FOR THE HELMHOLTZ EQUATON , 001 1, 1 and 1 T rlativ rrors, dfind by t rrors in t H 1 -sminorm dividd by t xact solution in t H 1 -sminorm, ar sown in t lft grap of Figur 2 W obsrv tat t rlativ rrors av similar baviors and dcay as ms siz bcoms smallr Tis sows tat t rrors ar not snsitiv to t paramtr β Nxt, w fix β = , and rpat t tst wit diffrnt δ T rigt grap of Figur 2 sows t rlativ rrors wit paramtrs δ = 0001, 01, 10 and 01, and wav numbrs k = 5 and 50 W obsrv tat t rrors av similar baviors for small valus δ = 0001 and 01 Largr δ rsults in largr rror T snsitivity tsts of t LDG mtod #2 ar sown in Figur 3, similar baviors ar also obsrvd Figur 2 Rlativ rror in t H 1 -sminorm of t LDG mtod #1 wit diffrnt paramtrs for two wav numbrs k = 5 and 50 Lft: δ = 01 is fixd, β = , 001 1, 1 and 1; Rigt: β = is fixd, δ = 0001, 01, 10 and Errors of t LDG solutions n tis subsction, w fix t paramtrs and invstigat t cangs of t numrical rrors as functions of t ms siz W start from t LDG mtod #1 As suggstd by t snsitivity tsts in t prvious subsction, w pick (54 δ = 01, β = T rlativ rror of t LDG mtod, and t finit lmnt intrpolation ar sown in t lft grap of Figur 4, wit four diffrnt wav numbrs k = 5, 10, 50 and 100 T rlativ rror of t LDG solution stays around 100% bfor a critical ms siz is racd, tn dcays at a rat gratr tan 1 in t log-log scal but convrgs as fast as t finit lmnt intrpolation (wit slop 1 for small T critical ms siz dcrass as k incrass T rigt grap of Figur 4 contains t rlativ rror wn w fix k = 1 and k = 05 t indicats tat unlik t rror of t finit lmnt intrpolation t rror of t LDG is not controlld by t magnitud of k, wic suggsts tat tr is a pollution contribution in t total rror T lft grap of Figur 5 contains

30 30 XAOBNG FENG AND YULONG XNG Figur 3 Rlativ rror in t H 1 -sminorm of t LDG mtod #2 wit diffrnt paramtrs for two wav numbrs k = 5 and 50 Lft: δ = 01 is fixd, β = , 001 1, 1 and 1; Rigt: β = is fixd, δ = 0001, 01, 10 and 01 t rlativ rror of t LDG mtod wit t ms siz satisfying k 3 2 = 1 for diffrnt valus of T rror dos not incras wit rspct to k T LDG mtod #2 as also bn tstd using t sam paramtrs in (54 Similar baviors av bn obsrvd as sown in Figur 5 and 6 Figur 4 Lft: rlativ rror of t LDG mtod #1 (solid lin and t finit lmnt intrpolation (dottd lin in t H 1 - sminorm for k = 5, 10, 50 and 100; Rigt: rlativ rror of t LDG mtod #1 in t H 1 -sminorm for k = 1,, 200, k = 1 and k = 05 At t nd, w look closly at t situation wit a larg rlativ rror wn k > 1 T LDG mtod #1 solution wit paramtrs δ = 01, β = , k = 100 and = 1/45 as a larg rlativ rror of siz T surfac plots of t finit lmnt intrpolation and t LDG solution ar givn in Figur 7 t

31 LDG METHODS FOR THE HELMHOLTZ EQUATON 31 Figur 5 Lft: rlativ rror of t LDG mtod #1 in t H 1 - sminorm wit k 3 2 = 1; Rigt: rlativ rror of t LDG mtod #2 (solid lin and t finit lmnt intrpolation (dottd lin in t H 1 -sminorm for k = 5, 10, 50 and 100 Figur 6 Rlativ rror of t LDG mtod #2 in t H 1 - sminorm Lft: k = 1 and k = 05 Rigt: k 3 2 = 1 sows tat t LDG solution as t corrct sap/pas altoug its amplitud is smallr 53 Comparison btwn t two LDG mtods Two diffrnt LDG mtods ar proposd in tis papr T first on is drivd following t PDG mtod proposd in [10], and t scond on as a mor standard numrical flux formulation and is supposd to av a bttr approximation for t vctor/flux variabl n tis subsction, w provid a comparison btwn ts two mtods, in trms of t rror and computational cost W start by rvisiting t tst xampls of Subsction 51 nstad of computing t rlativ rror of u in t H 1 -sminorm, w comput t rlativ rror of σ in t L 2 -norm for t LDG mtod #2 T numrical rsults ar prsntd in Figur

32 32 XAOBNG FENG AND YULONG XNG Figur 7 Lft: surfac plots of t finit lmnt intrpolation (lft and t LDG mtod #1 solution (rigt wit paramtrs δ = 01, β = , k = 100 and = 1/45 8, wic sow tat altoug t solution is still not snsitiv to t paramtr β, bttr approximation to σ is acivd for largr δ t confirms our prdiction tat t LDG mtod #2 givs a bttr approximation for t vctor/flux variabl Figur 8 Rlativ rror of σ in L 2 -norm of t LDG mtod #2 wit diffrnt paramtrs for two wav numbrs k = 5 and 50 Lft: δ = 01 is fixd, β = , 001 1, 1 and 1; Rigt: β = is fixd, δ = 0001, 01, 10 and 01 Tabl 1 provids a dtaild comparison of ts two mtods for diffrnt ms sizs, wit t paramtrs δ = 01, β = and k = 10 t sows tat t computational cost of t LDG mtod #2 is about twic largr tan tat of t LDG mtod #1 Also, as xpctd, t rror of t vctor/flux variabl of t LDG mtod #2 is smallr tan tat of t LDG mtod #1, and bot mtods dmonstrat a first ordr rat of convrgnc 54 Comparison btwn LDG and finit lmnt solutions W av sown t prformanc and comparison of t two LDG mtods in prvious subsctions

33 LDG METHODS FOR THE HELMHOLTZ EQUATON 33 Tabl 1 Comparison of t two LDG mtods wit paramtrs δ = 01 and β = / u u H 1 ordr σ σ L 2 ordr CPU tim (s E E E E LDG # E E E E E E E E E E E E LDG # E E E E E E E E Figur 9 T tracs of t LDG #1 solution (lft and t finit lmnt solution (rigt in t xz-plan, for k = 100 and = 1/50 (top, 1/120 (middl and 1/200 (bottom, rspctivly T dottd lins ar t tracs of t xact solution n tis subsction, w provid a brif comparison btwn t LDG solution and t P 1 conforming finit lmnt solution

34 34 XAOBNG FENG AND YULONG XNG W considr t Hlmoltz problm (51-(52 wit wav numbr k = 100 Wit ms siz = 1/50, 1/120 and 1/200, w plot t tracs of t LDG mtod #1 solution wit paramtrs (54 in xz-plan in t lft column of Figur 9 T xact solution is also providd as a rfrnc T tracs of t finit lmnt solution ar sown in t rigt column of Figur 9 t is clar tat t LDG mtod #1 as a bttr approximation to t xact solution Largr pas rror in t finit lmnt solution is obsrvd in all tr cass Also, t LDG solution as a bttr approximation for t amplitud of t xact solution Rfrncs [1] D Arnold, F Brzzi, B Cockburn, and D Marini Unifid analysis of discontinuous Galrkin mtods for lliptic problms SAM J Numr Anal, 39: , 2002 [2] S C Brnnr and L R Scott T matmatical tory of finit lmnt mtods 3rd d Txts in Applid Matmatics, vol 15, Springr, Nw York, 2008 [3] P Castillo, B Cockburn, Prugia and D Scötzau Local discontinuous Galrkin mtod for lliptic problms Commun Numr Mt Engrg, 18:69-75, 2002 [4] B Cockburn, G E Karniadakis, C-W Su Discontinuous Galrkin Mtods, Tory, Computation, and Applications Springr lctur Nots in Computational Scinc and Enginring, vol 11, Springr Vrlag, 2000 [5] B Cockburn and C -W Su T local discontinuous Galrkin mtod for convctiondiffusion systms SAM J Numr Anal, 35: , 1998 [6] P Cummings Analysis of Finit Elmnt Basd Numrical Mtods for Acoustic Wavs, Elastic Wavs and Fluid-Solid ntractions in t Frquncy Domain PD tsis, T Univrsity Tnnss, 2001 [7] P Cummings and X Fng Sarp rgularity cofficint stimats for complx-valud acoustic and lastic Hlmoltz quations M 3 AS, 16: , 2006 [8] J Douglas, Jr, J E Santos, D Sn, and L S Bnntum Frquncy domain tratmnt of on-dimnsional scalar wavs Mat Modls Mtods Appl Sci, 3(2: , 1993 [9] B Engquist and A Majda Radiation boundary conditions for acoustic and lastic wav calculations Comm Pur Appl Mat, 32(3: , 1979 [10] X Fng and H Wu Discontinuous Galrkin mtods for t Hlmoltz quation wit larg wav numbrs SAM J Numr Anal, 47: , 2009 [11] X Fng and H Wu p-discontinuous Galrkin mtods for t Hlmoltz quation wit larg wav numbrs Mat Comp (to appar [12] R Grismair and P Monk Error analysis for a ybridizabl discontinuous Galrkin mtod for t Hlmoltz quation prprint [13] U Htmaniuk Stability stimats for a class of Hlmoltz problms Commun Mat Sci, 5(3: , 2007 [14] R Hiptmair and Prugia Mixd plan wav discontinuous Galrkin mtods Lctur Nots in Computational Scinc and Enginring, 70:51-62, 2009 [15] F lnburg Finit Elmnt Analysis of Acoustic Scattring Springr-Vrlag, Nw York, 1998 [16] F lnburg and Babuška Finit lmnt solution of t Hlmoltz quation wit ig wav numbr T -vrsion of t FEM Comput Mat Appl, 30(9:9 37, 1995 [17] T Luostari, T Huttunn, and P Monk Plan wav mtods for approximating t tim armonic wav quation Higly oscillatory problms London Mat Soc Lctur Nots, 366: , 2009 [18] J M Mlnk On gnralizd finit lmnt mtods PD tsis, Univrsity of Maryland at Collg Park, 1995 [19] B Rivièr Discontinuous Galrkin mtods for solving lliptic and parabolic quations SAM, Piladlpia, PA, 2008 [20] O C Zinkiwicz Acivmnts and som unsolvd problms of t finit lmnt mtod ntrnat J Numr Mtods Engrg, 47:9 28, 2000

35 LDG METHODS FOR THE HELMHOLTZ EQUATON 35 Dpartmnt of Matmatics, T Univrsity of Tnnss, Knoxvill, TN 37996, USA addrss: xfng@matutkdu Dpartmnt of Matmatics, T Univrsity of Tnnss, Knoxvill, TN 37996, Computr Scinc and Matmatics Division, Oak Ridg National Laboratory, Oak Ridg, TN 37830, USA addrss: xingy@matutkdu