Running head: Weighted ENO Schemes Prfs sent t: Xu-Dng Liu Curant Institute f Mathematical Sciences 5 Mercer Street New Yrk, NY 00 AMS(MOS) subject cl
|
|
- Denis Phillips
- 6 years ago
- Views:
Transcription
1 Weighted Essentially Nn-Oscillatry Schemes Xu-Dng Liu Stanley Osher y Tny Chan z Department f Mathematics, UCLA, Ls Angeles, Califrnia, Research supprted by NSF grant DMS y Department f Mathematics, UCLA, Ls Angeles, Califrnia, Research supprted by NSF DMS and ONR NOO04-9-J-890 z Department f Mathematics,UCLA, Ls Angeles, Califrnia, Research supprted by NSF ASC , ARO DAAL03-9-G-50 and ONR NOO04-90-J-695
2 Running head: Weighted ENO Schemes Prfs sent t: Xu-Dng Liu Curant Institute f Mathematical Sciences 5 Mercer Street New Yrk, NY 00 AMS(MOS) subject classicatins: Primary 65M0; Secndary 65M05 Keywrds: Hyperblic Cnservatin Laws, ENO
3 Abstract In this paper we intrduce a new versin f ENO (Essentially Nn- Oscillatry) shck-capturing schemes which we call Weighted ENO. The main new idea is that, instead f chsing the \smthest" stencil t pick ne interplating plynmial fr the ENO recnstructin, we use a cnvex cmbinatin f all candidates t achieve the essentially nn-scillatry prperty, while additinally btaining ne rder f imprvement in accuracy. The resulting Weighted ENO schemes are based n cell-averages and a TVD Runge-Kutta time discretizatin. Preliminary encuraging numerical experiments are given. 3
4 Intrductin In this paper we present a new versin f ENO (Essentially Nn-Oscillatry) schemes. The cell-average versin f ENO schemes riginally was intrduced and develped by Harten and Osher in [] and Harten, Engquist, Osher and Chakravarthy in []. Later Shu and Osher develped the ux versin f ENO schemes and intrduced the TVD Runge-Kutta time discretizatin in [3] and [4]. The ENO schemes wrk well in many numerical experiments. The new ENO schemes which we call the Weighted ENO schemes are based n cell-averages and the TVD Runge-Kutta time discretizatin. The nly dierence between these schemes and the standard cell-average versin f ENO is hw we dene a recnstructin prcedure which prduces a high-rder accurate glbal apprximatin t the slutin frm its given cellaverages. The cell-average versin f ENO schemes attempts t avid grwth f spurius scillatins by an adaptive-stencil apprach, in which each cell is assigned its wn stencil f cells fr the purpses f recnstructin. Fr each cell the cell-average versin f ENO schemes selects an interplating stencil in which the slutin is smthest in sme sense. Thus a cell near a discntinuity is assigned a stencil frm the smth part f the slutin and a Gibbs-like phenmenn is s avided (see [5]). The Weighted ENO schemes develped here fllw this basic idea by using a cnvex cmbinatin apprach, in which each cell is assigned all crrespnding stencils and a cnvex cmbinatin f all crrespnding interplating plynmials n the stencils is cmputed t be the apprximating plynmial. This is dne by assigning prper weights t the cnvex cmbinatin. T achieve the essentially nn-scillatry prperty as the cell-average versin f ENO, the Weighted ENO schemes require that the cnvex cmbinatin be essentially a cnvex cmbinatin f the interplating plynmials n the smth stencils and that the interplating plynmials n the discntinuus stencils have essentially n cntributin t the cnvex cmbinatin. Thus, as in the cell-average versin f ENO schemes, a cell near a discntinuity is essentially assigned stencils frm the smth part f the slutin and a Gibbs-like phenmenn is als avided. In additin t this, the cnvex cmbinatin apprach results in cancellatin f truncatin errrs f crrespnding interplating plynmials and imprves the rder f accuracy by ne. Anther pssible advantage f Weighted ENO is smther dependence n data which may lessen sme f ENO's scillatry behavir near cnvergence and may help in getting a cnvergence prf. 4
5 In x we intrduce sme ntatins and basic ntins and give the TVD Runge-Kutta time discretizatin. In x3 we describe the prcedure f recnstructin frm given cell averages. In x4 we present sme preliminary numerical experiments. Basic Frmulatin and TVD Runge-Kutta Time Discretizatin We cnsider a hyperblic cnservatin law u t + f(u) x = 0; u(x; 0) = u (x): (:) Let fi j g be a partitin f R, where I j = [x j? ; x j+ ] is the j-th cell, x j+? = h. Dente fu j (; t)g t be the sliding averages f the weak slutin x j? u(x; t) f (.) i.e. R u j (; t) = h I j u(x; t) dx: (:) Integrating (.) ver each cell I j, we btain that the sliding averages fu j (; t)g j(; t) =? [f(u(x h j+ ; t))? f(u(x j? ; t))]: (:3) T evaluate j(; t), we need t evaluate f(u(x; t)) at each interface x j+. First f all, frm given cell-averages u = fu j g in which u j apprximates u j (; t), we recnstruct the slutin t btain R(x) = fr j (x)g which is a piecewise plynmial with unifrm plynmial degree r?, and in which each R j (x) is a plynmial apprximating u(x; t) n I j. We shall shw hw t btain R(x) frm u = fu j g in x3. Next at each interface x j+ have tw apprximating values R j (x j+ ) and R j+ (x j+ ) fr u(x j+, R(x) may ; t). We need a tw-pint Lipschitz mntne ux ~ h(; ) which is nndecreasing fr the rst argument and nnincreasing fr the secnd argument. Sme pssible chices are (i) Engquist-Osher h EO R (a; b) = b 0 R min(f 0 (s); 0) ds + a max(f 0 (s); 0) ds + f(0); (:4) 0 5
6 (ii) Gdunv h G (a; b) = ( minaub f(u) if a b; max aub f(u) if a > b; (:5) (iii) Re with entrpy x 8 >< f(a) if f 0 (u) 0 fr u [min(a; b); max(a; b)]; h RF (a; b) = f(b) if f 0 (u) 0 fr u [min(a; b); max(a; b)]; >: h LLF (a; b) therwise; (:6a) where h LLF (a; b) is dened as h LLF (a; b) = [f(a) + f(b)? (b? a)]; = max j f 0 (u) j : min(a;b)umax(a;b) (:6b) We apprximate f(u(x j+ and f(u(x j? where ; t)) by ~ h(r j? (x j? ; t)) by ~ h(r j (x j+ ); R j (x j(; t) L j (u); ); R j+ (x j+ )) )). Therefre (:7a) L j (u) =? h [~ h(r j (x j+ ); R j+ (x j+ ))? ~ h(r j? (x j? ); R j (x j? ))]: (:7b) In sectin x3, in which we intrduce the recnstructin prcedure, we shall btain that, in each cell I j, u(x; t) = R j (x) + O(h r ) 8x I j ; (:8) and at ne chsen pint f tw end pints f I j, u(x j; t) = R j (x j) + O(h r+ ) x j = x j? r x j = x j+ : (:9) Here and belw we always cnsider smth slutins when we discuss accuracy. Fr general upwind schemes, away frm snic pints (where f 0 (u) = 0), ( 0 f(a) in the regins f f > 0 ~h(a; b) = f(b) in the regins f f 0 < 0 6
7 In the regins f f 0 > 0, frm (.7b), and if we chse x j = x j+ hence L j (u) =? [f(r h j(x j+ ))? f(r j? (x j? ))]; u(x j+ u(x j? and x j? = x j? ; t) = R j (x j+ ; t) = R j? (x j? in (.9) i.e. ) + O(h r+ ) ) + O(h r+ j(; t) = L j (u) + O(h r+ ): (:0) Similarly, we shall have the abut frmula (.0) in the regins f f 0 < 0 by chsing x j = x j? and x j+ = x j+ in (.9). This will be detailed in x3:4. As usual, in the regins arund f 0 = 0 (snic pints), we j(; t) = L j (u) + O(h r ): Fr high rder time discretizatin, because f (.0), we need (r + )-th rder TVD Runge-Kutta time discretizatins intrduced by Shu and Osher in [3]. We need nly t spell ut the 3rd and 4th rder methds, which will be implemented in ur numerical experiments. Fr 3rd rder, 8j, Fr 4th rder, 8j, u (0) j = u n j ; u () j = u (0) u () j u (3) j u n+ j = u(0) j = 9 u(0) j = 3 u() j u (0) j = u n j ; u () j = u (0) u () j u n+ j j + L j (u (0) ) + u() + 9 u() j + 3 u() j = 3 4 u(0) j = 3 u(0) j j + L j (u (0) ) + 4 u() + 3 u() j + L 4 j(u () ) j + L 3 j(u () ): j? 4 L j(u (0) ) + L j(u () ) + 3 u() + 3 u(3) j? L 9 j(u (0) )? L 3 j(u () ) + L j (u () ) j + L 6 j(u () ) + L 6 j(u (3) ): T cmplete the cnstructin f ur schemes we frm ur nvel recnstructin prcedure. 7
8 3 Recnstructin Prcedure 3. Purpses f Recnstructin In this sectin we present the recnstructin prcedure. The R(x) is required t satisfy (i) In each cell I j, 8x I j and ne chsen pint x j I j, we have and R j (x) = u(x; t) + O(h r ); R j (x j) = u(x j; t) + O(h r+ ); (3:a) (3:b) where (3.b) will lead t ne rder f imprvement in accuracy, see x3:4 in this paper. (ii) R(x) has cnservatin frm i.e. 8j h R I j R j (x) dx = u j : (3:) (iii) Every R j (x) achieves the \ENO prperty" which will be specied later. 3. Interplatin n Each Stencil Fllwing the recnstructin prcedure in [], given the cell averages fu j g, we can immediately evaluate the pint values f the slutin's primitive functin W (x) at interfaces fw (x j+ )g, where the primitive functin is dened as W (x) = R x x u(x; t) dx; 0 j? (3:3) where x j 0? and bviusly culd be any interface, hence u(x; t) = W 0 (x) = d W (x); (3:4) dx P W (x j+ ) = j u i h: (3:5) i=j 0 T recnstruct the slutin, we interplate W (x) n each stencil S j = (x j?r+ ; x j?r+ 3 ; ; x j+ ) t btain a plynmial p j (x) i.e. p j (x l+ ) = W (x l+ ); l = j? r; ; j: 8
9 Obviusly the crrespnding plynmial p 0 j(x) (with degree r? ) apprximates the slutin u(x; t) i.e. see []. u(x; t) = p 0 j(x) + O(h r ) 8x (x j?r+ ; x j+ ); Als fr each stencil S j = (x j?r+ ; x j?r+ 3 ; ; x j+ ), we dene an indi- catr f the smthness IS j f u(x; t) n S j as fllwing: First we cmpute a table f dierences f fu j g n S j, where [u j?r+ ]; [u j?r+ ]; ; [u j? ]; [u j?r+ ]; [u j?r+ ]; ; [u j? ];. r? [u j?r+ ]; [u l ] = u l+? u l k [u l ] = k? [u l+ ]? k? [u l ]: Next we dene IS j t be the summatin f all averages f square values f the same rder dierences, That is, fr r =, and, fr r = 3, IS j = r? P P ( l ( r?l [u j?r+k ]) )=l: l= k= IS j = ([u j? ]) ; IS j = (([u j? ]) + ([u j? ]) )= + ( [u j? ]) : We bserve that if u(x; t) is discntinuus n S j, IS j O(), and if u(x; t) is cntinuus n S j, IS j O(h ). Hence fr each stencil S j, we btain p 0 j(x) apprximating u(x; t) n S j and IS j indicating the smthness f u(x; t) n S j. In the fllwing subsectin, t recnstruct the slutin in I j, we shall use r interplating plynmials fp 0 j+k(x)g r? k=0 n the stencils fs j+k g r? k=0, in which all S j+k cver the I j, t btain a cnvex cmbinatin f them, and we shall explre fis j+k g r? k=0 t assign a prper weight fr each f fp 0 j+k(x)g r? k=0 in the cnvex cmbinatin fr the purpses f recnstructin. 9
10 0 3.3 Cnvex Cmbinatin f fp (x)gr? j+k k=0 fr Each Cell I j Fr each cell I j we have r stencils fs j+k g r? k=0 = f(x j+k?r+ ; x j+k?r+ 3 ; ; x j+k+ )g r? k=0 which all include tw end pints x j? and x j+ f I j. We als have r interplating plynmials fp 0 j+k(x)g r? k=0 n the crrespnding stencils fs j+k g r? k=0. The main idea f the cell-average versin f ENO is t chse the \smthest" ne frm these r interplating plynmials. Fr Weighted ENO, instead f chsing ne, we use all r interplating plynmials and cmpute a cnvex cmbinatin f them t btain a plynmial R j (x) as fllws R j (x) = r? P j k r? P k=0 j l l=0 p 0 j+k(x); (3:6) where the j k > 0 (k = 0; ; ; ; r? ). Obviusly u(x; t) = R j (x) + O(h r ) in the smth regins f u(x; t) which is the purpse f (3.a). Nte that fr any k = 0; ; ; r?, p j+k (x j? ) = W (x j? ) and p j+k (x j+ ) = W (x j+ ), hence we achieve the purpse f (3.) R h I j R j (x) dx = r? P j k h r? P k=0 j l l=0 (p j+k (x j+ = h fw (x j+ )? W (x j? )? p j+k (x j? )) r? P )g k=0 j k r? P j l l=0 = u j : (3:7) Nte that n matter hw we dene f j k gr? k=0, R j (x) satises the purpses f (3.a) and (3.). We specify the \ENO prperty" f R j (x) by the crrespnding f j k gr? k=0. Denitin : The R j (x) has the \ENO prperty" if the crrespnding f j k gr? k=0 satisfy that (i) If the stencil S j+k is in the smth regins, the crrespnding j k satisfy j k P r? l=0 j l = O(): (3:8a) (ii) If the stencil S j+k is in a discntinuus regin f the slutin u(x; t), the crrespnding j k satisfy j k P r? l=0 j l O(h r ): (3:8b) 0
11 r? P Nte that, if f j k gr? k=0 satisfy the \ENO prperty" (3.8), the R j (x) = j k r? P k=0 j l l=0 p 0 j+k(x) will be a cnvex cmbinatin f the interplating plynmials n the smth stencils (3.8a), and the interplating plynmials n the discntinuus stencils have essentially n cntributin t R j (x) (3.8b). Dene, j k = C j k=( + IS j+k ) r ; k = 0; ; r? ; (3:9) where C j k = O() and C j k > 0 will be dened later fr imprvement f accuracy. Nte that because IS j+k culd be zer and =x is t sensitive as x is near zer, we add a small psitive number = 0?5 in the denminatr. Nte that if the stencil S j+k is in the smth regins j k P r? l=0 j l = O(); and if the stencil S j+k is in the discntinuus regins f u(x; t) j k P r? l=0 j l max(o( r ); O(h r )): Hence these f j k gr? k=0 (3.9) satisfy the \ENO prperty" (3.8)(O( r ) O(0?0 )). Here we assume there is at least ne stencil f fs j+k g r? k=0 in the smth regins. N matter hw we dene the cnstants fc j k gr? k=0, we have achieved the purpses f (3.a), (3.) and the \ENO prperty" (3.8). Hwever we shall specify fc j k gr? k=0 fr (3.b) which will lead ut ne rder imprvement in accuracy in sectin x3.4, ur last purpse f the recnstructin. Fr analysis we assume that in [x j?r+ ; x j+r+ ]. u(x; t) C r+ ; (3:0)
12 Fr each p 0 j+k(x), we express its truncatin errr as e j+k (x) where a j k(x) = r P = u(x; t)? p 0 j+k(x) = W 0 (x)? p 0 j+k(x) = d fw [x; x rq dx j+k?r+ ; ; x j+k+ ] (x? x j+k?l+ )g l=0 = d W [x; x rq dx j+k?r+ ; ; x j+k+ ] (x? x j+k?l+ ) l=0 rp rq +W [x; x j+k?r+ ; ; x j+k+ ] f (x? x j+k?l+ s=0 l=0;l6=s = W [x; x j+k?r+ ; ; x j+k+ ] a j k(x) + O(h r+ ); s=0 rq f l=0;l6=s (x? x j+k?l+ )g. We express the truncatin errr fr R j (x) E j (x) = u(x; t n )? R j (x) = W 0 (x)? R j (x) = r? P j k r? P (W 0 (x)? p 0 j+k(x)) = r? P P k=0 j l l=0 j k r? k=0 j l l=0 Because f the assumptin (3.0), 8k = 0; ; ; r?, j IS j+k j O(h ) j IS j+k? IS j j O(h ) j a j k(x) j O(h r ) j W [x; x j+k?r+ ; ; x j+k+ We have, frm (3.a) and (3.b), E j (x) = r? P j k r? P k=0 j l l=0 = r? P j k r? P k=0 j l l=0 r? C j k r? P = f P k=0 C j k l=0 e j+k (x) W [x; x j+k?r+ ]? W [x; x j?r+ a j k(x)g W [x; x j?r+ e j+k (x): ; ; x j+ ] j O(h): ; ; x j+k+ ] a j k(x) + O(h r+ ) ; ; x j+ ] + O(h r+ ): )g (3:a) (3:b) (3:c) The idea is that fr ne chsen pint x j [x j? ; x j+ ], we dene C j k t make the rst term in (3.c) equal t zer and btain E j (x j) = O(h r+ ):
13 Fr x j [x j? ; x j+ ], we dente p be the number f psitive terms in fa j k(x j)g r? and k=0 n be the number f negative terms in fa j k(x j)g r?, then we k=0 dene 8 if a j >< k(x j) = 0, C j h r k = if a j pja j k (x)j k (x j) > 0, (3:) j h >: r if a j nja k(x j) < 0. j k (x)j j Obviusly the C j k are independent f grid size h. E j (x j) r? P = f k=0 = f C j k r? P C j l l=0 P a j k(x j)gw [x; x j?r+ ; ; x j+ ] + O(h r+ ) p r? P a j k (x)>0 C j j l l=0 = 0 + O(h r+ ) = O(h r+ ):? P a j k (x j )<0 n r? P C j l l=0 gw [x; x j?r+ ; ; x j+ ] + O(h r+ ) (3:3) Remark : We have t have p and n t guarantee (3.3). Thus we btain that, fr ne chsen pint x j and any ther pint x [x j? ; x j+ ], dening C j k by (3.) gives us E j (x) = O(h r ); (3:4a) and E j (x j) = O(h r+ ): (3:4b) Up t nw, we have achieved all purpses f recnstructin (3.a), (3.b), (3.) and (3.8). 3.4 One Order Imprvement in Accuracy using (3.b) In this subsectin, we shall see hw (3.b) r (3.4b) gives us ne rder f imprvement in accuracy by chsing x j prperly in each cell. Let us cnsider the numerical spatial apprximatin (.7b) L j (u) =? h [ ~ h(rj (x j+ ~h(r j? (x j? 3 ); R j+ (x j+ ))? ))]: ); R j (x j?
14 Cnsider three cells in a smth regin, say cells I j?, I j and I j+, which are away frm snic pints. If f 0 (R(x)) > 0 in the cells, we have In (3.b), we chse x j = x j+ that Thus L j (u) =? [f(r h j(x j+ ))? f(r j? (x j? ))]: R j (x j+ R j? (x j? )? u(x j+ )? u(x j? L j (u) =? h [f(u(x j+ If f 0 (R(x)) < 0 in the cells, we have In (3.b), we chse x j = x j? that and x j? = x j?, then by (3.4b) we btain ; t) = E j (x j+ ) = O(h r+ ); ; t) = E j? (x j? ) = O(h r+ ): ; t))? f(u(x j? ; t))] + O(h r+ ): L j (u) =? [f(r h j+(x j+ ))? f(r j (x j? ))]: R j+ (x j+ R j (x j? )? u(x j+ )? u(x j? and x j+ = x j+, then by (3.4b) we btain ; t) = E j+ (x j+ ) = O(h r+ ); ; t) = E j (x j? ) = O(h r+ ): Thus L j (u) =? [f(u(x h j+ ; t))? f(u(x j? ; t))] + O(h r+ ): Hence in the smth regins and away frm snic pints, the numerical spatial peratrs fl j (u)g apprximate j(; t)g t the rder O(h r+ ). We specify x j in each I j in the fllwing way: First we cmpute f 0 (u j ). Then (i) if f 0 (u j ) > 0 we chse x j = x j+, (ii) if f 0 (u j ) < 0 we chse x j = x j?, (iii) if f 0 (u j ) = 0 we chse x j = x j+ r x j = x j?. If the cell I j is in the smth regins and away frm snic pints, then in general f 0 (R(x)) f 0 (u j ) > 0 arund the cell I j, hence accrding t the u j(; t) = L j (u) + O(h r+ ): (3:5) 4
15 Because snic pints are islated, in general, we btain (3.5) in mst f the cells u j(; t) = L j (u) + O(h r ) in a bunded, in fact small, number f cells near which there are snic pints as h decreases t zer. Remark 3: We have achieved ne rder imprvement in accuracy. Fr r = and r = 3, the cst f cmputing f the Weighted ENO schemes is cmparable t (f curse a little mre expensive than) that f standard ENO schemes (with the same rder accuracy) n sequential cmputers. Hwever n parallel cmputers, t achieve the same rder accuracy, the frmer schemes are much less expensive than the latter because the latter need mre expensive data transprt between cells. 3.5 Schemes fr r = The purpse f the fllwing tw subsectins x3:5 and x3:6 is t spell ut the details f the general schemes fr tw specic values f r, perhaps t aid the reader in implementin. In this subsectin, we cnsider ur schemes when r =. In this case we use linear interplatin t achieve the \ENO prperty" and 3rd rder accuracy (in ur numerical experiments, we achieved 4th rder accuracy) with cnservatin frm. Here we give the recnstructin prcedure fr r =. Fr each cell I j, we have tw stencils S j = (x j? 3 crrespnding t I j = [x j? linear interplatins and ; x j+ ; x j? ; x j+ ) and S j+ = (x j? ; x j+ ; x j+ 3 ) ]. On these tw stencils, we btain tw p 0 j(x) = u j + u j?u j? h (x? x j ) p 0 j+(x) = u j + u j+?u j h (x? x j ); and tw indicatrs f smthness IS j = (u j?u j? ) and IS j+ = (u j+?u j ). The recnstructed slutin R j (x) will be a cnvex cmbinatin f p 0 j(x) and p 0 j+(x) i.e. R j (x) = j 0 p 0 j(x) + j p 0 j 0 +j j+(x); (3:6) j 0 +j 5
16 where j 0 = C0=( j + IS j ), j = C=( j + IS j+ ). We shall specify C j 0 and C j in the fllwing tw cases. Case : If f 0 (u j ) > 0, we chse x j = x j+. We cmpute a j 0(x j+ ) = h and a j (x j+ ) =?h, and btain p = and n =, hence C j 0 = = and C j =. Thus in (3.6). j 0 = (+IS j ) j = (+IS j+ ) Case : If f 0 (u j ) 0, we chse x j = x j? and a j (x j? C j = =. Thus in (3.6). (3:7a). We cmpute a j 0(x j? ) =?h ) = h, and btain p = and n =, hence C j 0 = and j 0 = (+IS j ) j = (+IS j+ ) (3:7b) 3.6 Schemes fr r = 3 In this subsectin, we cnsider ur schemes when r = 3. In this case we use quadratic interplatin t achieve the \ENO prperty" and 4th rder accuracy (in ut numerical experiments, we achieved 5th rder accuracy) with cnservatin frm. Here we give ut the recnstructin prcedure fr r = 3. Fr each I j, we have three stencils S j = (x j? 5 ; x j? 3 ; x j? ; x j+ ), S j+ = (x j? 3 ; x j? ; x j+ ; x j+ 3 ), and S j+ = (x j? ; x j+ ; x 3 j+ ; x 5 j+ ) crrespnding t I j = [x j? three stencils, we btain three quadratic interplatins and p 0 j(x) = u j?u j? +u j? h (x? x j? ) + u j?u j? h (x? x j? )+ u j?? u j?u j? +u j? 4 p 0 j+(x) = u j+?u j +u j? h (x? x j ) + u j+?u j? h (x? x j )+ u j? u j+?u j +u j? 4 p 0 j+(x) = u j+?u j+ +u j h (x? x j+ ) + u j+?u j (x? x h j+ )+ u j+? u j+?u j+ +u j 4 ; 6 ; x j+ ]. On these
17 and three indicatrs f smthness IS j = ((u j?? u j? ) + (u j? u j? ) )= + (u j?u j? +u j? ), IS j+ = ((u j?u j? ) +(u j+?u j ) )=+(u j+?u j +u j? ) and IS j+ = ((u j+? u j ) + (u j+? u j+ ) )= + (u j+? u j+ + u j ). The recnstructed slutin R j (x) will be a cnvex cmbinatin f p 0 j(x), p 0 j+(x) and p 0 j+(x) i.e. R j (x) = j 0 j 0 +j +j p 0 j(x) + j j 0 +j +j p 0 j+(x) + j j 0 +j +j p 0 j+(x); (3:8) where j 0 = C0=( j + IS j ) 3, j = C=( j + IS j+ ) 3, j = C=( j + IS j+ ) 3. We shall specify C0, j C j and C j in the fllwing tw cases. Case : If f 0 (u j ) > 0, we chse x j = x j+. We cmpute a j 0(x j+ ) = 6h 3, a j (x j+ ) =?h 3 and a j (x j+ ) = h 3, and btain p = and n =, hence C j 0 = =, C j = = and C j = =4. Thus j 0 = (+IS j ) 3 j = (+IS j+ ) 3 j = 4(+IS j+ ) 3 (3:9a) in (3.8). Case : If f 0 (u j ) 0, we chse x j = x j??h 3,a j (x j? ) = h 3, and a j (x j? hence C j 0 = =4, C j = = and C j = =. Thus. We cmpute a j 0(x j? ) = ) =?6h 3 and btain p = and n =, j 0 = 4(+IS j ) 3 j = (+IS j+ ) 3 j = (+IS j+ ) 3 (3:9b) in (3.8). 4 Numerical Experiments 4. Scalar Cnservatin Laws In this subsectin we use sme mdel prblems t numerically test ur schemes. We use the Re ux with entrpy x as numerical ux and chse 7
18 r = which means we use a linear plynmial t recnstruct the slutin, and/r r = 3 which means we use a quadratic plynmial t recnstruct the slutin, and we expect t achieve 3rd and 4th rder accuracy respectively (at least away frm snic pints) accrding t ur analysis in the previus sectin. Example. We slve the mdel equatin u t + u x = 0? x u(x; 0) = u 0 (x) u 0 (x) peridic with perid : (4:) Five dierent initial data u 0 (x) are used. The rst ne is u 0 (x) = sin(x) and we list the errrs at time t = in Table. The secnd ne is u 0 (x) = sin 4 (x) and we list the errrs at time t = in Table. TABLE (=h = 0:8; t = ) l L errr L rder L errr L rder r = 80.77D-03.D D D D D r = D-05.03D D D D D D D TABLE (=h = 0:8; t = ) 8
19 l L errr L rder L errr L rder r = 80.77D-0 7.3D D D D D D D r = D D D D D D D D Here and belw l is the ttal number f cells and the step size h = =l in all scalar examples. Fr the rst tw initial data, we btain abut 4th (fr r = ) and 5th (fr r = 3) rder f accuracy respectively in the smth regin in bth L and L nrms which is surprisingly better than the 3rd and 4-th rder, the theretical results. We nte that standard ENO schemes applied t the example with the secnd initial data experienced an (easily xed) lss f accuracy, see [6], [7]. N such degeneracy was fund with ur present methds. The third initial functin is (? u 0 (x) = therwise, the furth is and the fth is u 0 (x) = ( (? ( 0 3 x) )? 3 0 x 3 0 ; 0 therwise; u 0 (x) = e?300x : We see the gd reslutin f the slutins in Figures -3 which are btained by ur scheme with r = 3. Linear discntinuities are smeared a bit. We expect t x this in the future using either the subcell reslutin technique f Harten [0] r the articial cmpressin technique f Yang [] tgether with the present technique. Figure (=h = 0:8) 9
20 The slutin by WENO at T = Number f Pints = true slu ++ apprx. slu r = 3 Figure (=h = 0:8) The slutin by WENO at T = Number f Pints = true slu ++ apprx. slu r = 3 0
21 Figure 3 (=h = 0:8) The slutin by WENO at T = Number f Pints = true slu ++ apprx. slu r = 3 Example. We slve Burgers' equatin with a peridic bundary cnditin u t + ( u ) x = 0? x (4:) u(x; 0) = u 0 (x) u 0 (x) peridic with perid : Fr the initial data u 0 (x) = + sin(x), the exact slutin is smth up t t =, then it develps a mving shck which interacts with a rarefactin wave. Observe that there is a snic pint. At t = 0:5 the slutin is still smth. We list the errrs in Table 3. Nte we als have abut 5th (fr r = 3) rder f accuracy respectively bth in L and L nrms.
22 TABLE 3 (=h = 0:6; t = 0:5) l L errr L rder L errr L rder r = D-05.84D D D D D At t = the shck just begins t frm, at t = 0:55 the interactin between the shck and the rarefactin waves is ver, and the slutin becmes mntne between shcks. In Figures 4-5 which are btain by ur scheme with r = 3 we can see the excellent behavir f the schemes in bth cases. The errrs at a distance 0: away frm the shck (i.e. j x? shck lcatin j 0:) are listed in Table 4 at t = 0:55. These errrs are f same magnitude as the nes in the smth case f Table 3 and shw abut 5th (fr r = 3) rder f accuracy respectively bth in L and L in the smth regins 0: away frm the shck. This shws that the errr prpagatin f the scheme is still very lcal. TABLE 4 (=h = 0:6; t = 0:55) l L errr L rder L errr L rder r = D D D D D D
23 Figure 4 (=h = 0:6).5 The slutin by WENO at T = Number f Pints = true slu ++ apprx. slu r = 3 Figure 5 (=h = 0:6).5 The slutin by WENO at T = 0.55 Number f Pints = true slu ++ apprx. slu r = 3 Example 3. we use tw nncnvex uxes t test the cnvergence t the physically crrect slutins. The true slutins are btained frm the Lax- Friedrichs scheme n a very ne grid. We use ur scheme with r = 3 in this example. The rst ne is a Riemann prblem with the ux f(u) = 3
24 4 (u? )(u? 4), and the initial data ( ul x < 0 u 0 (x) = u r x > 0: The tw cases we test are (i) u l =, u r =?, Figure 6; (ii) u l =?3, u r = 3, Figure 7. Fr mre details cncluding this prblem see [] Figure 6 (=h = 0:3) The slutin by WENO at T =.5 h = true slu apprx. slu l=80 Figure 7 (=h = 0:04) The slutin by WENO at T = 0.04 h = true slu apprx. slu l=80 4
25 The secnd ux is the Buckley-Leverett ux used t mdel il recvery [], f(u) = 4u =(4u + (? u) ), with initial data u = in [? ; 0] and u = 0 elsewhere. The result is displayed in Figure 8. Figure 8 (=h = 0:4) The slutin by WENO at T = h = true slu apprx. slu l=80 In this example, we bserve cnvergence with gd reslutin t the entrpy slutins in bth cases. In all the examples that we have illustrated abve, we bserve that the schemes are f abut 4th (fr r = ) and 5th (fr r = 3) rder f accuracy respectively and cnvergent with gd reslutin t the entrpy slutins. 4. Euler Equatins f Gas Dynamics In this subsectin we apply ur schemes t the Euler equatin f gas dynamics fr a plytrpic gas, u t + f(u) x = 0 u = (; m; E) T f(u) = qu + (0; P; qp ) T P = (? )(E? q ) m = q; (4:) where = :4 in the fllwing cmputatin. Fr details f the Jacbian, its eigenvalues, eigenvectrs, etc., see []. 5
26 Example 4. We cnsider the fllwing Riemann prblems: ( ul x < 0 u 0 (x) = u r x > 0: Tw sets f initial data are used. One is prpsed by Sd in [8]: ( l ; q l ; P l ) = (; 0; ); ( r ; q r ; P r ) = (0:5; 0; 0:): The ther is used by Lax [9]: ( l ; q l ; P l ) = (0:445; 0:698; 3:58); ( r ; q r ; P r ) = (0:5; 0; 0:57): We test ur schemes with r = 3. We use the characteristic recnstructin and Re ux with entrpy x frmed by Re's average as numerical ux. Fr details see []. The results are displayed in Figure 9-0. Figure 9a (=h = 0:4; t = ) DENSITY at time T = The Number f Pints = 00 Figure 9b (=h = 0:4; t = ) VELOCITY at time T = The Number f Pints = 00 6
27 Figure 9c (=h = 0:4; t = ) PRESSURE at time T = The Number f Pints = 00 Figure 0a (=h = 0:; t = :5).4 DENSITY at time T = The Number f Pints = 00 Figure 0b (=h = 0:; t = :5).6 VELOCITY at time T = The Number f Pints = 00 7
28 Figure 0c (=h = 0:; t = :5) 4 PRESSURE at time T = The Number f Pints = 00 Example 5. In this example we shall test the accuracy f ur schemes (r = 3) fr the Euler equatin f gas dynamics fr a plytrpic gas. We chse initial data as = +sin(x), m = +sin(x) and E = +sin(x), and peridic bundary cnditin. The true slutin was btained by applying the schemes t a very ne grid. Fr time t = when shcks haven't frmed, ur schemes achieve 5th (r = 3) rder accuracy in all three cmpnents, see Table 5. We can als see the slutin fr time t = in Figure. 8
29 TABLE 5 (=h = 0:6; t = ) l L errr L rder L errr L rder DENSITY 80.99D-04.9D D D D D D D MOMENTUM 80.7D-04.76D D D D D D D ENERGY 80.0D-04.9D D D D D D D
30 Figure a (=h = 0:6; t = ).3 Density The Number f Pints = 00 Figure b (=h = 0:6; t = ).6 Mmentum The Number f Pints = 00 Figure c (=h = 0:6; t = ).4 Energy The Number f Pints = 00 Acknwledgment We are grateful t Prfessr Chi-Wang Shu fr his 30
31 suggestin f the smth indicatr functin i.e. instead f IS j = r? P P ( l ( r?l [u j?r+k ]) )=l; l= l= k= IS j = r? P P ( l j r?l [u j?r+k ] j)=l; k= which we used riginally. Bth functins wrk well, hwever the latter ne leads t a smther (C vs. Lipschitz) numerical ux which may be helpful fr steady state cnvergence r cnvergence prf. References [] A. Harten and S. Osher,\Unifrmly High-Order Accurate Nn- Oscillatry Schemes I," SIAM Jurnal n Numerical Analysis, V4, pp , 987; als MRC Technical Summary Reprt N. 83, May 985. [] A. Harten, B. Engquist, S. Osher and S. Chakravarthy, \Unifrmly High Order Accurate Essentially Nn-Oscillatry Schemes III," Jurnal f Cmputatinal Physics, V7, pp , 987; als ICASE Reprt N. 86-, April 986. [3] C.-W. Shu and S. Osher, \Ecient Implementatin f Essentially Nn-scillatry Shck-Capturing Schemes," Jurnal f Cmputatinal Physics, V77, 988, pp [4] Chi-Wang Shu, Stanley Osher, \Ecient Implementatin f Essentially Nn-scillatry Shck-Capturing Schemes, II," J. Cmput. Phys., V83, 989, pp [5] A. Harten and S. Chakraverthy, \Multi-Dimensinal ENO Schemes fr General Gemetries," UCLA CAM reprt N. 9-6, August 99. [6] A. Rgersn and E. Meiburg, \A Numerical Study f the Cnvergence Prperties f ENO Schemes," J. Scientic Cmputing, V5, N., 990, pp
32 [7] Chi-Wang Shu, \Numerical Experiments n the Accuracy f ENO and Mdied ENO Schemes," J. Scientic Cmputing, V5, N., 990, pp [8] G. Sd, J. Cmput. Phys. 7, (978). [9] P. Lax. Cmmun. Pure Appl. Math. 46, (986). [0] A. Harten, \ENO Schemes with Subcell Reslutin," J. Cmput. Phys., V83 (989), pp [] H. Yang, \ An Articial Cmpressin Methd fr ENO Schemes, the Slpe Mdicatin Methd," J. Cmput. Phys., V89 (990) pp
Lyapunov Stability Stability of Equilibrium Points
Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),
More information1996 Engineering Systems Design and Analysis Conference, Montpellier, France, July 1-4, 1996, Vol. 7, pp
THE POWER AND LIMIT OF NEURAL NETWORKS T. Y. Lin Department f Mathematics and Cmputer Science San Jse State University San Jse, Califrnia 959-003 tylin@cs.ssu.edu and Bereley Initiative in Sft Cmputing*
More informationBootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >
Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);
More informationCOMP 551 Applied Machine Learning Lecture 11: Support Vector Machines
COMP 551 Applied Machine Learning Lecture 11: Supprt Vectr Machines Instructr: (jpineau@cs.mcgill.ca) Class web page: www.cs.mcgill.ca/~jpineau/cmp551 Unless therwise nted, all material psted fr this curse
More informationCOMP 551 Applied Machine Learning Lecture 9: Support Vector Machines (cont d)
COMP 551 Applied Machine Learning Lecture 9: Supprt Vectr Machines (cnt d) Instructr: Herke van Hf (herke.vanhf@mail.mcgill.ca) Slides mstly by: Class web page: www.cs.mcgill.ca/~hvanh2/cmp551 Unless therwise
More informationPattern Recognition 2014 Support Vector Machines
Pattern Recgnitin 2014 Supprt Vectr Machines Ad Feelders Universiteit Utrecht Ad Feelders ( Universiteit Utrecht ) Pattern Recgnitin 1 / 55 Overview 1 Separable Case 2 Kernel Functins 3 Allwing Errrs (Sft
More informationChapter 3 Kinematics in Two Dimensions; Vectors
Chapter 3 Kinematics in Tw Dimensins; Vectrs Vectrs and Scalars Additin f Vectrs Graphical Methds (One and Tw- Dimensin) Multiplicatin f a Vectr b a Scalar Subtractin f Vectrs Graphical Methds Adding Vectrs
More informationHomology groups of disks with holes
Hmlgy grups f disks with hles THEOREM. Let p 1,, p k } be a sequence f distinct pints in the interir unit disk D n where n 2, and suppse that fr all j the sets E j Int D n are clsed, pairwise disjint subdisks.
More informationAdmissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs
Admissibility Cnditins and Asympttic Behavir f Strngly Regular Graphs VASCO MOÇO MANO Department f Mathematics University f Prt Oprt PORTUGAL vascmcman@gmailcm LUÍS ANTÓNIO DE ALMEIDA VIEIRA Department
More informationV. Balakrishnan and S. Boyd. (To Appear in Systems and Control Letters, 1992) Abstract
On Cmputing the WrstCase Peak Gain f Linear Systems V Balakrishnan and S Byd (T Appear in Systems and Cntrl Letters, 99) Abstract Based n the bunds due t Dyle and Byd, we present simple upper and lwer
More informationDetermining the Accuracy of Modal Parameter Estimation Methods
Determining the Accuracy f Mdal Parameter Estimatin Methds by Michael Lee Ph.D., P.E. & Mar Richardsn Ph.D. Structural Measurement Systems Milpitas, CA Abstract The mst cmmn type f mdal testing system
More informationQuantum Harmonic Oscillator, a computational approach
IOSR Jurnal f Applied Physics (IOSR-JAP) e-issn: 78-4861.Vlume 7, Issue 5 Ver. II (Sep. - Oct. 015), PP 33-38 www.isrjurnals Quantum Harmnic Oscillatr, a cmputatinal apprach Sarmistha Sahu, Maharani Lakshmi
More informationOF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
U. S. FOREST SERVICE RESEARCH PAPER FPL 50 DECEMBER U. S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
More informationKinetic Model Completeness
5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins
More informationLecture 17: Free Energy of Multi-phase Solutions at Equilibrium
Lecture 17: 11.07.05 Free Energy f Multi-phase Slutins at Equilibrium Tday: LAST TIME...2 FREE ENERGY DIAGRAMS OF MULTI-PHASE SOLUTIONS 1...3 The cmmn tangent cnstructin and the lever rule...3 Practical
More informationThe blessing of dimensionality for kernel methods
fr kernel methds Building classifiers in high dimensinal space Pierre Dupnt Pierre.Dupnt@ucluvain.be Classifiers define decisin surfaces in sme feature space where the data is either initially represented
More informationinitially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur
Cdewrd Distributin fr Frequency Sensitive Cmpetitive Learning with One Dimensinal Input Data Aristides S. Galanpuls and Stanley C. Ahalt Department f Electrical Engineering The Ohi State University Abstract
More information, which yields. where z1. and z2
The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin
More informationFall 2013 Physics 172 Recitation 3 Momentum and Springs
Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.
More informationThermodynamics and Equilibrium
Thermdynamics and Equilibrium Thermdynamics Thermdynamics is the study f the relatinship between heat and ther frms f energy in a chemical r physical prcess. We intrduced the thermdynamic prperty f enthalpy,
More informationSupport-Vector Machines
Supprt-Vectr Machines Intrductin Supprt vectr machine is a linear machine with sme very nice prperties. Haykin chapter 6. See Alpaydin chapter 13 fr similar cntent. Nte: Part f this lecture drew material
More informationLeast Squares Optimal Filtering with Multirate Observations
Prc. 36th Asilmar Cnf. n Signals, Systems, and Cmputers, Pacific Grve, CA, Nvember 2002 Least Squares Optimal Filtering with Multirate Observatins Charles W. herrien and Anthny H. Hawes Department f Electrical
More informationSmoothing, penalized least squares and splines
Smthing, penalized least squares and splines Duglas Nychka, www.image.ucar.edu/~nychka Lcally weighted averages Penalized least squares smthers Prperties f smthers Splines and Reprducing Kernels The interplatin
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 2: Mdeling change. In Petre Department f IT, Åb Akademi http://users.ab.fi/ipetre/cmpmd/ Cntent f the lecture Basic paradigm f mdeling change Examples Linear dynamical
More informationCS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007
CS 477/677 Analysis f Algrithms Fall 2007 Dr. Gerge Bebis Curse Prject Due Date: 11/29/2007 Part1: Cmparisn f Srting Algrithms (70% f the prject grade) The bjective f the first part f the assignment is
More informationLead/Lag Compensator Frequency Domain Properties and Design Methods
Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin
More informationFebruary 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA
February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCurseWare http://cw.mit.edu 2.161 Signal Prcessing: Cntinuus and Discrete Fall 2008 Fr infrmatin abut citing these materials r ur Terms f Use, visit: http://cw.mit.edu/terms. Massachusetts Institute
More informationChapter 3: Cluster Analysis
Chapter 3: Cluster Analysis } 3.1 Basic Cncepts f Clustering 3.1.1 Cluster Analysis 3.1. Clustering Categries } 3. Partitining Methds 3..1 The principle 3.. K-Means Methd 3..3 K-Medids Methd 3..4 CLARA
More informationDead-beat controller design
J. Hetthéssy, A. Barta, R. Bars: Dead beat cntrller design Nvember, 4 Dead-beat cntrller design In sampled data cntrl systems the cntrller is realised by an intelligent device, typically by a PLC (Prgrammable
More informationPressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects
Pressure And Entrpy Variatins Acrss The Weak Shck Wave Due T Viscsity Effects OSTAFA A. A. AHOUD Department f athematics Faculty f Science Benha University 13518 Benha EGYPT Abstract:-The nnlinear differential
More informationx x
Mdeling the Dynamics f Life: Calculus and Prbability fr Life Scientists Frederick R. Adler cfrederick R. Adler, Department f Mathematics and Department f Bilgy, University f Utah, Salt Lake City, Utah
More informationPerfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart
Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Sandy D. Balkin Dennis K. J. Lin y Pennsylvania State University, University Park, PA 16802 Sandy Balkin is a graduate student
More informationDifferentiation Applications 1: Related Rates
Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm
More informationMath 302 Learning Objectives
Multivariable Calculus (Part I) 13.1 Vectrs in Three-Dimensinal Space Math 302 Learning Objectives Plt pints in three-dimensinal space. Find the distance between tw pints in three-dimensinal space. Write
More informationENGI 4430 Parametric Vector Functions Page 2-01
ENGI 4430 Parametric Vectr Functins Page -01. Parametric Vectr Functins (cntinued) Any nn-zer vectr r can be decmpsed int its magnitude r and its directin: r rrˆ, where r r 0 Tangent Vectr: dx dy dz dr
More informationDepartment of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets
Department f Ecnmics, University f alifrnia, Davis Ecn 200 Micr Thery Prfessr Giacm Bnann Insurance Markets nsider an individual wh has an initial wealth f. ith sme prbability p he faces a lss f x (0
More informationEQUADIFF 6. Erich Martensen The ROTHE method for nonlinear hyperbolic problems. Terms of use:
EQUADIFF 6 Erich Martensen The ROTHE methd fr nnlinear hyperblic prblems In: Jarmír Vsmanský and Milš Zlámal (eds.): Equadiff 6, Prceedings f the Internatinal Cnference n Differential Equatins and Their
More informationAPPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL
JP2.11 APPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL Xingang Fan * and Jeffrey S. Tilley University f Alaska Fairbanks, Fairbanks,
More informationSections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.
Tpic : AC Fundamentals, Sinusidal Wavefrm, and Phasrs Sectins 5. t 5., 6. and 6. f the textbk (Rbbins-Miller) cver the materials required fr this tpic.. Wavefrms in electrical systems are current r vltage
More informationLHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers
LHS Mathematics Department Hnrs Pre-alculus Final Eam nswers Part Shrt Prblems The table at the right gives the ppulatin f Massachusetts ver the past several decades Using an epnential mdel, predict the
More informationRevision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax
.7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins
More informationIAML: Support Vector Machines
1 / 22 IAML: Supprt Vectr Machines Charles Suttn and Victr Lavrenk Schl f Infrmatics Semester 1 2 / 22 Outline Separating hyperplane with maimum margin Nn-separable training data Epanding the input int
More informationRevisiting the Socrates Example
Sectin 1.6 Sectin Summary Valid Arguments Inference Rules fr Prpsitinal Lgic Using Rules f Inference t Build Arguments Rules f Inference fr Quantified Statements Building Arguments fr Quantified Statements
More informationALE 21. Gibbs Free Energy. At what temperature does the spontaneity of a reaction change?
Name Chem 163 Sectin: Team Number: ALE 21. Gibbs Free Energy (Reference: 20.3 Silberberg 5 th editin) At what temperature des the spntaneity f a reactin change? The Mdel: The Definitin f Free Energy S
More informationOptimization Programming Problems For Control And Management Of Bacterial Disease With Two Stage Growth/Spread Among Plants
Internatinal Jurnal f Engineering Science Inventin ISSN (Online): 9 67, ISSN (Print): 9 676 www.ijesi.rg Vlume 5 Issue 8 ugust 06 PP.0-07 Optimizatin Prgramming Prblems Fr Cntrl nd Management Of Bacterial
More informationTime-domain lifted wavelet collocation method for modeling nonlinear wave propagation
Lee et al.: Acustics Research Letters Online [DOI./.] Published Online 8 August Time-dmain lifted wavelet cllcatin methd fr mdeling nnlinear wave prpagatin Kelvin Chee-Mun Lee and Wn-Seng Gan Digital Signal
More informationCOASTAL ENGINEERING Chapter 2
CASTAL ENGINEERING Chapter 2 GENERALIZED WAVE DIFFRACTIN DIAGRAMS J. W. Jhnsn Assciate Prfessr f Mechanical Engineering University f Califrnia Berkeley, Califrnia INTRDUCTIN Wave diffractin is the phenmenn
More informationModeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function
www.ccsenet.rg/mer Mechanical Engineering Research Vl. 1, N. 1; December 011 Mdeling the Nnlinear Rhelgical Behavir f Materials with a Hyper-Expnential Type Functin Marc Delphin Mnsia Département de Physique,
More informationLecture 5: Equilibrium and Oscillations
Lecture 5: Equilibrium and Oscillatins Energy and Mtin Last time, we fund that fr a system with energy cnserved, v = ± E U m ( ) ( ) One result we see immediately is that there is n slutin fr velcity if
More informationInterference is when two (or more) sets of waves meet and combine to produce a new pattern.
Interference Interference is when tw (r mre) sets f waves meet and cmbine t prduce a new pattern. This pattern can vary depending n the riginal wave directin, wavelength, amplitude, etc. The tw mst extreme
More informationROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS
Particle Acceleratrs, 1986, Vl. 19, pp. 99-105 0031-2460/86/1904-0099/$15.00/0 1986 Grdn and Breach, Science Publishers, S.A. Printed in the United States f America ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS
More informationENSC Discrete Time Systems. Project Outline. Semester
ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding
More informationAP Statistics Notes Unit Two: The Normal Distributions
AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).
More informationA solution of certain Diophantine problems
A slutin f certain Diphantine prblems Authr L. Euler* E7 Nvi Cmmentarii academiae scientiarum Petrplitanae 0, 1776, pp. 8-58 Opera Omnia: Series 1, Vlume 3, pp. 05-17 Reprinted in Cmmentat. arithm. 1,
More informationA new Type of Fuzzy Functions in Fuzzy Topological Spaces
IOSR Jurnal f Mathematics (IOSR-JM e-issn: 78-578, p-issn: 39-765X Vlume, Issue 5 Ver I (Sep - Oct06, PP 8-4 wwwisrjurnalsrg A new Type f Fuzzy Functins in Fuzzy Tplgical Spaces Assist Prf Dr Munir Abdul
More informationMath Foundations 20 Work Plan
Math Fundatins 20 Wrk Plan Units / Tpics 20.8 Demnstrate understanding f systems f linear inequalities in tw variables. Time Frame December 1-3 weeks 6-10 Majr Learning Indicatrs Identify situatins relevant
More information39th International Physics Olympiad - Hanoi - Vietnam Theoretical Problem No. 1 /Solution. Solution
39th Internatinal Physics Olympiad - Hani - Vietnam - 8 Theretical Prblem N. /Slutin Slutin. The structure f the mrtar.. Calculating the distance TG The vlume f water in the bucket is V = = 3 3 3 cm m.
More informationModule 4: General Formulation of Electric Circuit Theory
Mdule 4: General Frmulatin f Electric Circuit Thery 4. General Frmulatin f Electric Circuit Thery All electrmagnetic phenmena are described at a fundamental level by Maxwell's equatins and the assciated
More informationBOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS. Christopher Costello, Andrew Solow, Michael Neubert, and Stephen Polasky
BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS Christpher Cstell, Andrew Slw, Michael Neubert, and Stephen Plasky Intrductin The central questin in the ecnmic analysis f climate change plicy cncerns
More informationREPRESENTATIONS OF sp(2n; C ) SVATOPLUK KR YSL. Abstract. In this paper we have shown how a tensor product of an innite dimensional
ON A DISTINGUISHED CLASS OF INFINITE DIMENSIONAL REPRESENTATIONS OF sp(2n; C ) SVATOPLUK KR YSL Abstract. In this paper we have shwn hw a tensr prduct f an innite dimensinal representatin within a certain
More informationPreparation work for A2 Mathematics [2017]
Preparatin wrk fr A2 Mathematics [2017] The wrk studied in Y12 after the return frm study leave is frm the Cre 3 mdule f the A2 Mathematics curse. This wrk will nly be reviewed during Year 13, it will
More informationWe can see from the graph above that the intersection is, i.e., [ ).
MTH 111 Cllege Algebra Lecture Ntes July 2, 2014 Functin Arithmetic: With nt t much difficulty, we ntice that inputs f functins are numbers, and utputs f functins are numbers. S whatever we can d with
More informationPreparation work for A2 Mathematics [2018]
Preparatin wrk fr A Mathematics [018] The wrk studied in Y1 will frm the fundatins n which will build upn in Year 13. It will nly be reviewed during Year 13, it will nt be retaught. This is t allw time
More informationLim f (x) e. Find the largest possible domain and its discontinuity points. Why is it discontinuous at those points (if any)?
THESE ARE SAMPLE QUESTIONS FOR EACH OF THE STUDENT LEARNING OUTCOMES (SLO) SET FOR THIS COURSE. SLO 1: Understand and use the cncept f the limit f a functin i. Use prperties f limits and ther techniques,
More informationSlide04 (supplemental) Haykin Chapter 4 (both 2nd and 3rd ed): Multi-Layer Perceptrons
Slide04 supplemental) Haykin Chapter 4 bth 2nd and 3rd ed): Multi-Layer Perceptrns CPSC 636-600 Instructr: Ynsuck Che Heuristic fr Making Backprp Perfrm Better 1. Sequential vs. batch update: fr large
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationCHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India
CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 11: Mdeling with systems f ODEs In Petre Department f IT, Ab Akademi http://www.users.ab.fi/ipetre/cmpmd/ Mdeling with differential equatins Mdeling strategy Fcus
More informationParticle Size Distributions from SANS Data Using the Maximum Entropy Method. By J. A. POTTON, G. J. DANIELL AND B. D. RAINFORD
3 J. Appl. Cryst. (1988). 21,3-8 Particle Size Distributins frm SANS Data Using the Maximum Entrpy Methd By J. A. PTTN, G. J. DANIELL AND B. D. RAINFRD Physics Department, The University, Suthamptn S9
More informationMargin Distribution and Learning Algorithms
ICML 03 Margin Distributin and Learning Algrithms Ashutsh Garg IBM Almaden Research Center, San Jse, CA 9513 USA Dan Rth Department f Cmputer Science, University f Illinis, Urbana, IL 61801 USA ASHUTOSH@US.IBM.COM
More informationRigid Body Dynamics (continued)
Last time: Rigid dy Dynamics (cntinued) Discussin f pint mass, rigid bdy as useful abstractins f reality Many-particle apprach t rigid bdy mdeling: Newtn s Secnd Law, Euler s Law Cntinuus bdy apprach t
More informationSOLUTION OF THREE-CONSTRAINT ENTROPY-BASED VELOCITY DISTRIBUTION
SOLUTION OF THREECONSTRAINT ENTROPYBASED VELOCITY DISTRIBUTION By D. E. Barbe,' J. F. Cruise, 2 and V. P. Singh, 3 Members, ASCE ABSTRACT: A twdimensinal velcity prfile based upn the principle f maximum
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationBuilding to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.
Building t Transfrmatins n Crdinate Axis Grade 5: Gemetry Graph pints n the crdinate plane t slve real-wrld and mathematical prblems. 5.G.1. Use a pair f perpendicular number lines, called axes, t define
More informationo o IMPORTANT REMINDERS Reports will be graded largely on their ability to clearly communicate results and important conclusions.
BASD High Schl Frmal Lab Reprt GENERAL INFORMATION 12 pt Times New Rman fnt Duble-spaced, if required by yur teacher 1 inch margins n all sides (tp, bttm, left, and right) Always write in third persn (avid
More informationPhysics 2010 Motion with Constant Acceleration Experiment 1
. Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin
More informationA New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation
III-l III. A New Evaluatin Measure J. Jiner and L. Werner Abstract The prblems f evaluatin and the needed criteria f evaluatin measures in the SMART system f infrmatin retrieval are reviewed and discussed.
More informationIntroduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem
A Generalized apprach fr cmputing the trajectries assciated with the Newtnian N Bdy Prblem AbuBar Mehmd, Syed Umer Abbas Shah and Ghulam Shabbir Faculty f Engineering Sciences, GIK Institute f Engineering
More informationDistributions, spatial statistics and a Bayesian perspective
Distributins, spatial statistics and a Bayesian perspective Dug Nychka Natinal Center fr Atmspheric Research Distributins and densities Cnditinal distributins and Bayes Thm Bivariate nrmal Spatial statistics
More information1 The limitations of Hartree Fock approximation
Chapter: Pst-Hartree Fck Methds - I The limitatins f Hartree Fck apprximatin The n electrn single determinant Hartree Fck wave functin is the variatinal best amng all pssible n electrn single determinants
More informationBASD HIGH SCHOOL FORMAL LAB REPORT
BASD HIGH SCHOOL FORMAL LAB REPORT *WARNING: After an explanatin f what t include in each sectin, there is an example f hw the sectin might lk using a sample experiment Keep in mind, the sample lab used
More informationDepartment: MATHEMATICS
Cde: MATH 022 Title: ALGEBRA SKILLS Institute: STEM Department: MATHEMATICS Curse Descriptin: This curse prvides students wh have cmpleted MATH 021 with the necessary skills and cncepts t cntinue the study
More informationTree Structured Classifier
Tree Structured Classifier Reference: Classificatin and Regressin Trees by L. Breiman, J. H. Friedman, R. A. Olshen, and C. J. Stne, Chapman & Hall, 98. A Medical Eample (CART): Predict high risk patients
More informationSection 6-2: Simplex Method: Maximization with Problem Constraints of the Form ~
Sectin 6-2: Simplex Methd: Maximizatin with Prblem Cnstraints f the Frm ~ Nte: This methd was develped by Gerge B. Dantzig in 1947 while n assignment t the U.S. Department f the Air Frce. Definitin: Standard
More informationYou need to be able to define the following terms and answer basic questions about them:
CS440/ECE448 Sectin Q Fall 2017 Midterm Review Yu need t be able t define the fllwing terms and answer basic questins abut them: Intr t AI, agents and envirnments Pssible definitins f AI, prs and cns f
More informationNAME: Prof. Ruiz. 1. [5 points] What is the difference between simple random sampling and stratified random sampling?
CS4445 ata Mining and Kwledge iscery in atabases. B Term 2014 Exam 1 Nember 24, 2014 Prf. Carlina Ruiz epartment f Cmputer Science Wrcester Plytechnic Institute NAME: Prf. Ruiz Prblem I: Prblem II: Prblem
More informationFIELD QUALITY IN ACCELERATOR MAGNETS
FIELD QUALITY IN ACCELERATOR MAGNETS S. Russenschuck CERN, 1211 Geneva 23, Switzerland Abstract The field quality in the supercnducting magnets is expressed in terms f the cefficients f the Furier series
More information^YawataR&D Laboratory, Nippon Steel Corporation, Tobata, Kitakyushu, Japan
Detectin f fatigue crack initiatin frm a ntch under a randm lad C. Makabe," S. Nishida^C. Urashima,' H. Kaneshir* "Department f Mechanical Systems Engineering, University f the Ryukyus, Nishihara, kinawa,
More informationPhys. 344 Ch 7 Lecture 8 Fri., April. 10 th,
Phys. 344 Ch 7 Lecture 8 Fri., April. 0 th, 009 Fri. 4/0 8. Ising Mdel f Ferrmagnets HW30 66, 74 Mn. 4/3 Review Sat. 4/8 3pm Exam 3 HW Mnday: Review fr est 3. See n-line practice test lecture-prep is t
More informationKinetics of Particles. Chapter 3
Kinetics f Particles Chapter 3 1 Kinetics f Particles It is the study f the relatins existing between the frces acting n bdy, the mass f the bdy, and the mtin f the bdy. It is the study f the relatin between
More information(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f
1.0 Review f Electrmagnetic Field Thery Selected aspects f electrmagnetic thery are reviewed in this sectin, with emphasis n cncepts which are useful in understanding magnet design. Detailed, rigrus treatments
More informationOn Boussinesq's problem
Internatinal Jurnal f Engineering Science 39 (2001) 317±322 www.elsevier.cm/lcate/ijengsci On Bussinesq's prblem A.P.S. Selvadurai * Department f Civil Engineering and Applied Mechanics, McGill University,
More informationAerodynamic Separability in Tip Speed Ratio and Separability in Wind Speed- a Comparison
Jurnal f Physics: Cnference Series OPEN ACCESS Aerdynamic Separability in Tip Speed Rati and Separability in Wind Speed- a Cmparisn T cite this article: M L Gala Sants et al 14 J. Phys.: Cnf. Ser. 555
More informationSequential Allocation with Minimal Switching
In Cmputing Science and Statistics 28 (1996), pp. 567 572 Sequential Allcatin with Minimal Switching Quentin F. Stut 1 Janis Hardwick 1 EECS Dept., University f Michigan Statistics Dept., Purdue University
More informationCHM112 Lab Graphing with Excel Grading Rubric
Name CHM112 Lab Graphing with Excel Grading Rubric Criteria Pints pssible Pints earned Graphs crrectly pltted and adhere t all guidelines (including descriptive title, prperly frmatted axes, trendline
More informationStudy Group Report: Plate-fin Heat Exchangers: AEA Technology
Study Grup Reprt: Plate-fin Heat Exchangers: AEA Technlgy The prblem under study cncerned the apparent discrepancy between a series f experiments using a plate fin heat exchanger and the classical thery
More informationA Simple Set of Test Matrices for Eigenvalue Programs*
Simple Set f Test Matrices fr Eigenvalue Prgrams* By C. W. Gear** bstract. Sets f simple matrices f rder N are given, tgether with all f their eigenvalues and right eigenvectrs, and simple rules fr generating
More information