Implementing a Convolutional Perfectly Matched Layer in a finite-difference code for the simulation of seismic wave propagation in a 2D elastic medium
|
|
- Roy Johnson
- 5 years ago
- Views:
Transcription
1 Implemenng a Conoluonal Perfecl Mached Laer n a fne-dfference code for he smulaon of sesmc wae propagaon n a D elasc medum Progress Repor BRGM/RP-559-FR December 7
2
3 Implemenng a Conoluonal Perfecl Mached Laer n a fne-dfference code for he smulaon of sesmc wae propagaon n a D elasc medum Progress Repor BRGM/RP-559-FR December 7 Sud carred ou as par of Research aces - BRGM 7 PDR7ARN4 Duceller A., Aoch H. Checked b: Name: Phlppe Jousse Dae: Sgnaure: Approed b: Name: Hormoz Modaress Dae: Sgnaure: BRGM's qual managemen ssem s cerfed ISO 9: b AFAQ IM 3 ANG Aprl 5
4 Kewords: Sesmc wae propagaon, Absorbng condon, Perfecl Mached Laer, Conoluonal Perfecl Mached Laer, Non-reflecng condons, Fne-dfference. In bblograph, hs repor should be ced as follows: Duceller A., Aoch H. (7) Implemenng a Conoluonal Perfecl Mached Laer n a fne-dfference code for he smulaon of sesmc wae propagaon n a D elasc medum, BRGM echncal repor, BRGM/RP-559-FR. BRGM, 5. No par of hs documen ma be reproduced whou he pror permsson of BRGM.
5 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon Snopss One of he mos popular mehods o smulae numercall he sesmc wae propagaon n an elasc meda s he fne-dfference mehod. In he cone of numercal modellng n unbounded meda, he wae needs o be absorbed a he arfcal boundares of he compuaonal doman and herefore, s necessar o defne non-reflecng condons a hese boundares o mmc an unbounded medum. The Perfecl Mached Laer (PML) condon has he remarkable proper of hang a zero reflecon coeffcen for all angles of ncdence and all frequences before dscrezaon and has become wdel used (e.g. Collno and Tsogka, ). Howeer, hs reflecon coeffcen s no zero anmore afer he dscrezaon and becomes een er large a grazng ncdence. Therefore, an mproed erson of he PML condon has been deeloped: he Conoluonal Perfecl Mached Laer (CPML) condon (e.g. Komasch and Marn, 7). The am of hs projec s o mplemen hese wo non-reflecng condons n a D fne-dfference compuaonal code for an elasc medum. We wll frs eplan he heorecal formulaon of he PML and CPML condons. Then, we wll epose he srucure of our compuaonal code. Fnall, we wll presen he resuls of some numercal ess performed o erf he effcenc of he mehod. BRGM/RP-559-FR Progress repor 3
6
7 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon Conens. Inroducon...7. The CPML formulaon CLASSICAL PML FORMULATION CPML FORMULATION Implemenaon n he D code GRID CELL BOUNDARY CONDITIONS SOURCE FONCTION STRUCTURE OF THE CODE Numercal ess PML CPML Concluson References...3 BRGM/RP-559-FR Progress repor 5
8 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon Ls of llusraons Fgure - Elemenar grd cell of he wo dmensonnal saggered spaal fnedfference mehod... 7 Fgure - Source me funcon mposed a pon source... Fgure 3 Tme eoluon of he and componens of he wo-dmensonal eloc ecor a he frs, second and hrd saons of he numercal soluon wh classcal PML boundares for he hn slce and he large slce Fgure 4 Snapshos of he and componens of he wo-dmensonal eloc ecor of he numercal soluon wh classcal PML boundares for he hn slce, a me sep 3, 5, 7, 9, and Fgure 5 Tme eoluon of he and componens of he wo-dmensonal eloc ecor a he frs, second and hrd saons of he numercal soluons wh CPML boundares for he hn slce and he large slce Fgure 6 Snapshos of he and componens of he wo-dmensonnal eloc ecor of he numercal soluon wh CPML boundares for he hn slce, a me sep 3, 5, 7, 9, and BRGM/RP-559-FR Progress repor
9 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon. Inroducon One of he mos popular mehods o smulae numercall he sesmc wae propagaon n an elasc meda s he fne-dfference mehod. In he cone of numercal modellng n unbounded meda, he wae needs o be absorbed a he arfcal boundares of he compuaonal doman and herefore, s necessar o defne non-reflecng condons a hese boundares o mmc an unbounded medum. The Perfecl Mached Laer (PML) condon has he remarkable proper of hang a zero reflecon coeffcen for all angles of ncdence and all frequences before dscrezaon and has become wdel used (e.g. Collno and Tsogka, ). Howeer, hs reflecon coeffcen s no zero anmore afer he dscrezaon and becomes een er large a grazng ncdence. Therefore, an mproed erson of he PML condon has been deeloped: he Conoluonal Perfecl Mached Laer (CPML) condon (e.g. Komasch and Marn, 7). The am of hs projec s o mplemen hese wo non-reflecng condons n a D fne-dfference compuaonal code for an elasc medum. We wll frs eplan he heorecal formulaon of he PML and CPML condons. Then, we wll epose he srucure of our compuaonal code. Fnall, we wll presen he resuls of some numercal ess performed o erf he effcenc of he mehod. BRGM/RP-559-FR Progress repor 7
10
11 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon. The CPML formulaon One of he mos popular mehods o smulae numercall he sesmc wae propagaon n an elasc meda s he fne-dfference mehod. In he cone of numercal modellng n unbounded meda, he wae needs o be absorbed a he arfcal boundares of he compuaonal doman and herefore, s necessar o defne non-reflecng condons a hese boundares o mmc an unbounded medum. The Perfecl Mached Laer (PML) condon has he remarkable proper of hang a zero reflecon coeffcen for all angles of ncdence and all frequences before dscrezaon and has become wdel used (e.g. Collno and Tsogka, ). Howeer, hs reflecon coeffcen s no zero anmore afer he dscrezaon and becomes een er large a grazng ncdence. Therefore, an mproed erson of he PML condon has been deeloped: he Conoluonal Perfecl Mached Laer (CPML) condon (Komasch and Marn, 7). We wll frs eplan he classcal PML formulaon, and hen we wll nroduce he CPML formulaon... CLASSICAL PML FORMULATION The dfferenal form of he D elasc wae equaon can be wren as: ρ ρ for he eloces ( and ) and: λ µ λ µ BRGM/RP-559-FR Progress repor 9
12 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon µ for he sress (, and ), where ρ s dens and λ and µ are elasc coeffcens. In he parcular case of a homogenous medum, hs equaon has plane wae r r r Aep k. ω. soluons of he form ( ( ) Le us consder a regular doman locaed n. In order o aod reflecons a, we need o defne an absorbng laer locaed n >. Ths laer wll be an absorbng laer f he amplude of he wae decreases eponenall n he drecon for >. To defne hs laer, we wll defne a new wae equaon, such ha s soluon s a plane wae n he regular doman and an eponenall-decang plane wae n he PML. We nroduce a new comple coordnae ~. Thus, he PML can be ewed as an analcal connuaon of he real coordnaes n he comple space (Collno and Tsogka, ). The coordnae ~ s defned b: ~ ω ( ) d ( s) ds Equaon wh d () beng a dumpng profle n he PML regon, such ha d nsde he regular doman and d > n he PML. Equaon ges us upon dfferenang: ω ~ ϖ d s Equaon wheres ω d ω d ω We wll now change he orgnal wae equaon wren n erms of arable and no a new wae equaon wren n erms of arables ~ and. We frs rewre he elasc wae equaon n he frequenc doman: ωρ ωρ BRGM/RP-559-FR Progress repor
13 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon BRGM/RP-559-FR Progress repor for he eloces and: µ λ ω µ λ ω µ ω for he sress. We hen replace he wae equaon wren n erms of wh a generalzed wae equaon wren n erms of ~. Insde he regular doman, boh equaons are dencal because d. ~ ωρ ~ ωρ for he eloces and: ~ ~ µ λ ω µ λ ω ~ ~ µ ω for he sress. In he PML, he new wae equaon has eponenall-decang plane wae soluons of he form: ( ) ( ) ( ) ( ) ( ) ds s d k k A k k A ep. ep ~ ep ω ω ω r r r r
14 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon BRGM/RP-559-FR Progress repor n he drecon. Usng Equaon : s ~, we rewre he wae equaon n he whole doman n erms of raher han ~ (n he regular doman, s ): s ωρ s ωρ for he eloces and: s s µ λ ω s µ λ ω s µ ω for he sress. The eloc and sress felds are hen spl no wo pars, and, such ha: s ωρ ωρ s ωρ ϖρ
15 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon BRGM/RP-559-FR Progress repor 3 for he eloces and: s s µ λ ω λ ω s λ ω µ λ ω s µ ω µ ϖ for he sress. Conerng back o he me doman, we fnall ge: d ρ ρ d ρ ρ for he eloces and: d λ µ
16 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon λ d λ λ µ d µ µ for he sress. Ths new wae equaon perms o hae an eponenall-decang plane wae soluon n he PML laer ( > ). Therefore he boundar a appears o be a nonreflecng boundar. We hen use he same echnque for he hree oher boundares of he regular doman. Ths formulaon perms o hae numercal soluons wh reduced arfcal reflecons bu s effcenc becomes poor a grazng ncdences. The CPML formulaon was deeloped o address hs ssue... CPML FORMULATION The man dea of he CPML formulaon s o choose a more general defnon of s b nroducng new arables α and κ : s d κ α ω Equaon 3 As hs epresson depends on frequenc, when we go back n he me doman, we wll ge a me conoluon on each modfed spaal derae. Denong s () he nerse Fourer ransform of s, we can wre: 4 BRGM/RP-559-FR Progress repor
17 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon BRGM/RP-559-FR Progress repor 5 ( ) s * ~ Equaon 4 * beng a conoluon operaor. We hae from Equaon 3: ( ) ( ) ( ) ( ) d H d s α κ κ κ δ ep ) ( Equaon 5 where δ() and H() are he Drac and Heasde funcons. If we denoe: ( ) ( ) ( ) ( ) d H d α κ κ ζ ep Equaon 6 we can see from Equaons 4, 5 and 6 ha: ( ) * ~ ζ κ Le us denoe n ϕ he conoluon erm compued a me n, wh beng he me sep: ( ) ( ) n n n d * τ τ ς ζ ϕ τ Equaon 7 Equaon 7 ges us b dscrezaon: ( ) ( ) ( ) ( ) ( ) n m m n n n m m m m n n n m m m n n m Z d d ϕ τ τ ς ϕ τ τ ς ϕ τ Equaon 8 wh: ( ) ( ) ( ) m m d m Z τ τ ς Equaon 9 From Equaons 6 and 9, we oban:
18 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon d ( m ) Z ( m) ep ( ( d ) ) m κ α τdτ a κ ( b ) m wh a d ( ) ( b ) and b ep ( ( d κ α ) ) d κ α κ We consder hen φ as a memor arable whose me eoluon s goerned a each sep b: ϕ n b ϕ n a n From a praccal pon of ew, o mplemen he CPML echnque n a compuaonal code, we jus need o replace he spaal derae b: ~ κ ϕ and updae φ a each me sep. 6 BRGM/RP-559-FR Progress repor
19 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon 3. Implemenaon n he D code 3.. GRID CELL We use a classcal second-order saggered grd n space and me (Madaraga, 976, and Vreu, 986): Fgure - Elemenar grd cell of he wo dmensonnal saggered spaal fne-dfference mehod of Madaraga (976) and Vreu (986) used classcall o dscreze he equaons of elasodnamcs. 3.. BOUNDARY CONDITIONS We pu non-reflecng boundares on he four sdes of he compuaonal doman (Komasch and Marn, 7) so as o smulae wae propagaon n nfne elasc medum. In hs case, he boundar condons are: (,. ) ( mn ma,.) (., jmn ) (., j ) ma and (,. ) ( mn ma,.) (., jmn ) (., j ) ma BRGM/RP-559-FR Progress repor 7
20 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon mn, ma, j mn and j ma beng he lms of he regular doman and beng he number of cells n he (C)PML laer SOURCE FONCTION The source s gen b a eloc ecor, he same as Komasch and Marn (7) for aldaon. A he pon source, we jus add: ( ) d ρ ( α ) a ( ) ep a ( ) sn f * Equaon o ( source,j source ) and: ( α ) a ( ) ep a ( ) ( ) d ρ cos f * Equaon o ( source ½, j source ½), beng he me, α, f, a and beng source parameers gen n he code and ρ beng he dens STRUCTURE OF THE CODE The srucure of he compuaonal code s he followng: Inalzaons Calculaon of he alues of d (), α (), κ (), d (j), α (j) and κ (j) defned n Equaon 3 Begnnng of he me loop Tme (n/) (*) 8 BRGM/RP-559-FR Progress repor
21 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon for j { for { f (,j n he CPML) { Calculaon of and n he CPML } else { Calculaon of and n he regular doman } } } for j { for { f (,j n he CPML) { Calculaon of n he CPML } else { Calculaon of n he regular doman } } } Tme (n) for j { for { f (,j n he CPML) { Calculaon of n he CPML } else { Calculaon of n he regular doman } } } for j { BRGM/RP-559-FR Progress repor 9
22 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon for { f (,j n he CPML) { Calculaon of n he CPML } else { Calculaon of n he regular doman } } } Addng a eloc ecor a he source Imposng he boundar condons Sang he alues n fles for he sesmograms and he snapshos Incremen of n and reurn o (*) BRGM/RP-559-FR Progress repor
23 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon 4. Numercal ess 4.. PML To es our compuaonal code, we use he same approach as Komasch and Marn (7). We consder wo models: he sze of he frs model s 8*64 meers, represenng a doman much longer han wde (.e., a hn slce). In hs model, he wae wll propagae a grazng ncdence along he edge, whch ma cause spurous reflecons. To hghlgh hese reflecons, we wll consder a second, larger model (34*64 meers), whch boundares are far enough from he source o aod an spurous reflecons, and we wll compare he resuls of he wo models. Ths hn model s dscrezed usng a grd comprsng 8 pons * 64 pons. The large model s dscrezed usng a grd comprsng 34 * 64 pons. The sde of a grd cell s m. The compressonal wae speed s Vp 33 m.s -, he shear wae speed s Vs Vp/ m.s - and he dens s ρ 8 kg.m -3. The me sep s. s and we perform he smulaon for me seps. The pon source s a eloc ecor locaed n source 6 m and source 48 m. The me funcon s gen b Equaons and. The parameers for he source are: α 35, f 7, a π*f d,. / f d and f d 7 Hz Source funcon 5 source (m/s) 5-5,8,36,54,7,9,8,6,44,6,8,98,6,34,5,7,88 3,6 3,4 3,4 3,6 3,78 3, me (s) Fgure - Source me funcon mposed a pon source The saons were we compue sesmograms are locaed a he followng pons: BRGM/RP-559-FR Progress repor
24 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon Saon number X coordnae Y coordnae Saon s locaed a he mddle of he model where here s no spurous reflecon. Saon and 3 are locaed along he edge a grazng and er grazng ncdence, where spurous reflecons could dsor he sesmograms. For he dumpng coeffcen n he PML, we use he followng profle: d ( ) d L N where s he dsance beween he pon where d s compued and he begnnng of he PML laer, and L s he hckness of he PML: L *, N and: d ( N ) log( R ) p c wh Rc.. L Absorbng laers of hckness grd pons are mplemened on he four sdes of he model. The resuls are shown hereafer n Fgures 3 and 4. Saon - V - PML Saon - V - PML,5,8,5,6,4 (m/s) -,5 - Thn slce Large slce (m/s), Thn slce Large slce -,5 -, - -,4 -,5 -,6,,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s),,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s) BRGM/RP-559-FR Progress repor
25 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon Saon - V - PML Saon - V - PML,5,8,5,6,4 (m/s) -,5 - Thn slce Large slce (m/s), c Thn slce Large slce -,5 -, - -,4 -,5 -,6,,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s),,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s) Saon 3 - V -PML Saon 3 - V - PML,8,5,6,4 (m/s) -,5 Thn slce Large slce (m/s), Thn slce Large slce -, - -,4 -,5 -,6,,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s),,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s) Fgure 3 Tme eoluon of he (lef) and (rgh) componens of he wo-dmensonal eloc ecor a he frs (op), second (mddle) and hrd (boom) saons of he numercal soluon wh classcal PML boundares for he hn slce (blue) and he large slce (pnk). A he frs saon, relael far from he begnnng of he PML laer and wh non-grazng ncdence, he agreemen s almos perfec. A he second saon, spurous oscllaons, n parcular for he S-wae of he componen, sar o appear. A he hrd receer, he oscllaons become large, he P-wae s no correcl calculaed and he shape of he S-wae s compleel dsored. BRGM/RP-559-FR Progress repor 3
26 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon Fgure 4 Snapshos of he and componens of he wo-dmensonal eloc ecor of he numercal soluon wh classcal PML boundares for he hn slce, a me sep 3, 5, 7 (lef, op o boom), 9, and 3 (rgh, op o boom). The hree dos ndcae he posons of he saons a whch he sesmograms represened n Fgure are recorded. Spurous waes appear a grazng ncdence along he edge of he model and send spurous energ back no he man doman. We fnd ha some large reflecon waes reman especall a saons and 3. In fac, he snapshos show some srong reflecon a he boundar closes o he source n he begnnng. We also obsere he oscllaon socked n he PML laer, whch propagaes parallel o he boundar. Ths s he problem of he PML preousl poned ou (e.g. Komasch and Marn, 7). 4 BRGM/RP-559-FR Progress repor
27 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon 4.. CPML To es he CPML, we use he same approach as for he PML. We choose o make α ar n a lnear fashon beween mamal alue α ma a he begnnng of he CPML and zero a s op. We ake α ma π*f d. We choose κ consan, equal o. The resuls are hereafer n Fgures 5 and 6. Saon - V - CPML Saon - V - CPML,5,8,5,6,4 (m/s) -,5 - Thn slce Large slce (m/s), Thn slce Large slce -,5 -, - -,4 -,5 -,6,,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s),,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s) Saon - V - PML Saon - V - PML,5,8,5,6,4 (m/s) -,5 - Thn slce Large slce (m/s), c Thn slce Large slce -,5 -, - -,4 -,5 -,6,,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s),,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s) Saon 3 - V - CPML Saon 3 - V - CPML,8,5,6,4 (m/s) -,5 Thn slce Large slce (m/s), Thn slce Large slce -, - -,4 -,5 -,6,,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s),,4,6,8,3,4,44,65,85,6,7,47,68,88 3,9 3,3 3,5 3,7 3,9 me (s) BRGM/RP-559-FR Progress repor 5
28 Implemenng a CPML n a D fne-dfference code for he smulaon of sesmc wae propagaon Fgure 5 Tme eoluon of he (lef) and (rgh) componens of he wo-dmensonal eloc ecor a he frs (op), second (mddle) and hrd (boom) saons of he numercal soluons wh CPML boundares for he hn slce (blue) and he large slce (pnk). A he frs saon, relael far from he begnnng of he CPML laer and wh non-grazng ncdence, he agreemen s almos prefec. A he second saon, wh grazng ncdence and close o he begnnng of he CPML laer, he agreemen remans ecellen. A he hrd receer, n he dffcul case of er grazng ncdence, of a long dsance of propagaon and of a saon locaed close o he begnnng of he CPML laer, he agreemen remans er sasfacor. 6 BRGM/RP-559-FR Progress repor
29 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon Fgure 6 Snapshos of he and componens of he wo-dmensonnal eloc ecor of he numercal soluon wh CPML boundares for he hn slce, a me sep 3, 5, 7 (lef, op o boom), 9, and 3 (rgh, op o boom). The hree dos ndcae he posons of he saons a whch he sesmograms represened n Fgure 4 are recorded. No spurous wae of sgnfcan amplude s sble, een a grazng ncdence. In he CPML case, we are able o ge he dencal resuls as Komasch and Marn (7) well mproed wh respec o he preous ones (PML). We do no hae an more he arfcal wae reflecon een n he snapshos. We confrm here ha he CPML absorbng condon works much beer. BRGM/RP-559-FR Progress repor 7
30
31 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon 5. Concluson The behaour of he Perfecl mached Laer a grazng ncdence s well mproed wh he new Conoluonnal PML approach. I can be useful when sudng wae propagaon n hn slces or when he source s locaed near he edge of he model. The cos of CPML n memor sorage s smlar o he classcal PML. BRGM/RP-559-FR Progress repor 9
32
33 Implemenng a CPML n D fne-dfference code for he smulaon of sesmc wae propagaon 6. References Collno F., Tsogka C. () - Applcaon of he PML absorbng laer model o he lnear elasodnamc problem n ansoropc heerogeneous meda. Geophscs, 66(), Graes R.W. (996) - Smulang Sesmc Wae propagaon n 3D Elasc Meda Usng Saggered-Grd Fne Dfferences. Bull. Sesm. Soc. Am., 86, 9-6. Komasch D., Marn R. (7) - An unspl conoluonal Perfecl Mached Laer mproed a grazng ncdence for he sesmc wae equaon. Geophscs, 7, SM55- SM67. Madaraga R. (976) - Dnamcs of an epandng crcular faul. Bull. Sesm. Soc. Am., 65, Vreu J. (986) - P-SV wae propagaon n heerogeneous meda: eloc-sress fne-dfference mehod. Geophscs, 5, BRGM/RP-559-FR Progress repor 3
34
35
36 Scenfc and Techncal Cenre Deelopmen Plannng and Naural Rsks Dson 3, aenue Claude-Gullemn - BP Orléans Cede France Tel.: 33 ()
Implementing a Convolutional Perfectly Matched Layer in a finite-difference code for the simulation of seismic wave propagation in a 3D elastic medium
Implemenng a Conoluonal Perfecl Mache Laer n a fne-fference coe for he mulaon of emc wae propagaon n a 3D elac meum Progre Repor BRGM/RP-559-FR December, 7 Implemenng a Conoluonal Perfecl Mache Laer n
More informationSolution in semi infinite diffusion couples (error function analysis)
Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of
More information[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5
TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres
More informationHEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD
Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,
More informationJ i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.
umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal
More informationVariants of Pegasos. December 11, 2009
Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on
More information( ) () we define the interaction representation by the unitary transformation () = ()
Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger
More informationTHE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS
THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he
More informationApproximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy
Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae
More informationChapter Lagrangian Interpolation
Chaper 5.4 agrangan Inerpolaon Afer readng hs chaper you should be able o:. dere agrangan mehod of nerpolaon. sole problems usng agrangan mehod of nerpolaon and. use agrangan nerpolans o fnd deraes and
More information2/20/2013. EE 101 Midterm 2 Review
//3 EE Mderm eew //3 Volage-mplfer Model The npu ressance s he equalen ressance see when lookng no he npu ermnals of he amplfer. o s he oupu ressance. I causes he oupu olage o decrease as he load ressance
More informationPHYS 1443 Section 001 Lecture #4
PHYS 1443 Secon 001 Lecure #4 Monda, June 5, 006 Moon n Two Dmensons Moon under consan acceleraon Projecle Moon Mamum ranges and heghs Reerence Frames and relae moon Newon s Laws o Moon Force Newon s Law
More informationP R = P 0. The system is shown on the next figure:
TPG460 Reservor Smulaon 08 page of INTRODUCTION TO RESERVOIR SIMULATION Analycal and numercal soluons of smple one-dmensonal, one-phase flow equaons As an nroducon o reservor smulaon, we wll revew he smples
More informationNotes on the stability of dynamic systems and the use of Eigen Values.
Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon
More informationNearly perfectly matched layer method for seismic wave propagation in poroelastic media
ρ u + ρ w = λ + µ u + µ u + α M w c ρ u + m w + b w = α M u + µ u + M w CANADIAN JOURNAL o EXPLORATION GEOPHYSlCS VOL 37, NO 1 June 01, P -7 Nearly perecly mached layer mehod or sesmc wave propagaon n
More informationOn One Analytic Method of. Constructing Program Controls
Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure
More informationCH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal
More informationLecture 18: The Laplace Transform (See Sections and 14.7 in Boas)
Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on
More informationScattering at an Interface: Oblique Incidence
Course Insrucor Dr. Raymond C. Rumpf Offce: A 337 Phone: (915) 747 6958 E Mal: rcrumpf@uep.edu EE 4347 Appled Elecromagnecs Topc 3g Scaerng a an Inerface: Oblque Incdence Scaerng These Oblque noes may
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ
More informationFTCS Solution to the Heat Equation
FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence
More informationLet s treat the problem of the response of a system to an applied external force. Again,
Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem
More informationVolatility Interpolation
Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4
CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu
More informationII. Light is a Ray (Geometrical Optics)
II Lgh s a Ray (Geomercal Opcs) IIB Reflecon and Refracon Hero s Prncple of Leas Dsance Law of Reflecon Hero of Aleandra, who lved n he 2 nd cenury BC, posulaed he followng prncple: Prncple of Leas Dsance:
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as
More informationCubic Bezier Homotopy Function for Solving Exponential Equations
Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.
More informationAdvanced time-series analysis (University of Lund, Economic History Department)
Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng
More informationCHAPTER 10: LINEAR DISCRIMINATION
CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g
More informationA NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION
S19 A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION by Xaojun YANG a,b, Yugu YANG a*, Carlo CATTANI c, and Mngzheng ZHU b a Sae Key Laboraory for Geomechancs and Deep Underground Engneerng, Chna Unversy
More informationResponse of MDOF systems
Response of MDOF syses Degree of freedo DOF: he nu nuber of ndependen coordnaes requred o deerne copleely he posons of all pars of a syse a any nsan of e. wo DOF syses hree DOF syses he noral ode analyss
More informationELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION
THE 19 TH INTERNATIONAL ONFERENE ON OMPOSITE MATERIALS ELASTI MODULUS ESTIMATION OF HOPPED ARBON FIBER TAPE REINFORED THERMOPLASTIS USING THE MONTE ARLO SIMULATION Y. Sao 1*, J. Takahash 1, T. Masuo 1,
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems
Swss Federal Insue of Page 1 The Fne Elemen Mehod for he Analyss of Non-Lnear and Dynamc Sysems Prof. Dr. Mchael Havbro Faber Dr. Nebojsa Mojslovc Swss Federal Insue of ETH Zurch, Swzerland Mehod of Fne
More informationComprehensive Integrated Simulation and Optimization of LPP for EUV Lithography Devices
Comprehense Inegraed Smulaon and Opmaon of LPP for EUV Lhograph Deces A. Hassanen V. Su V. Moroo T. Su B. Rce (Inel) Fourh Inernaonal EUVL Smposum San Dego CA Noember 7-9 2005 Argonne Naonal Laboraor Offce
More informationChapters 2 Kinematics. Position, Distance, Displacement
Chapers Knemacs Poson, Dsance, Dsplacemen Mechancs: Knemacs and Dynamcs. Knemacs deals wh moon, bu s no concerned wh he cause o moon. Dynamcs deals wh he relaonshp beween orce and moon. The word dsplacemen
More informationLi An-Ping. Beijing , P.R.China
A New Type of Cpher: DICING_csb L An-Png Bejng 100085, P.R.Chna apl0001@sna.com Absrac: In hs paper, we wll propose a new ype of cpher named DICING_csb, whch s derved from our prevous sream cpher DICING.
More informationDual Approximate Dynamic Programming for Large Scale Hydro Valleys
Dual Approxmae Dynamc Programmng for Large Scale Hydro Valleys Perre Carpener and Jean-Phlppe Chanceler 1 ENSTA ParsTech and ENPC ParsTech CMM Workshop, January 2016 1 Jon work wh J.-C. Alas, suppored
More informationDynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005
Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc
More informationMotion in Two Dimensions
Phys 1 Chaper 4 Moon n Two Dmensons adzyubenko@csub.edu hp://www.csub.edu/~adzyubenko 005, 014 A. Dzyubenko 004 Brooks/Cole 1 Dsplacemen as a Vecor The poson of an objec s descrbed by s poson ecor, r The
More informationMulti-Fuel and Mixed-Mode IC Engine Combustion Simulation with a Detailed Chemistry Based Progress Variable Library Approach
Mul-Fuel and Med-Mode IC Engne Combuson Smulaon wh a Dealed Chemsry Based Progress Varable Lbrary Approach Conens Inroducon Approach Resuls Conclusons 2 Inroducon New Combuson Model- PVM-MF New Legslaons
More informationStochastic Maxwell Equations in Photonic Crystal Modeling and Simulations
Sochasc Maxwell Equaons n Phoonc Crsal Modelng and Smulaons Hao-Mn Zhou School of Mah Georga Insue of Technolog Jon work wh: Al Adb ECE Majd Bade ECE Shu-Nee Chow Mah IPAM UCLA Aprl 14-18 2008 Parall suppored
More informationTHEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that
THEORETICAL AUTOCORRELATIONS Cov( y, y ) E( y E( y))( y E( y)) ρ = = Var( y) E( y E( y)) =,, L ρ = and Cov( y, y ) s ofen denoed by whle Var( y ) f ofen denoed by γ. Noe ha γ = γ and ρ = ρ and because
More informationIn the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!
ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal
More informationAdvanced Machine Learning & Perception
Advanced Machne Learnng & Percepon Insrucor: Tony Jebara SVM Feaure & Kernel Selecon SVM Eensons Feaure Selecon (Flerng and Wrappng) SVM Feaure Selecon SVM Kernel Selecon SVM Eensons Classfcaon Feaure/Kernel
More informationPart II CONTINUOUS TIME STOCHASTIC PROCESSES
Par II CONTINUOUS TIME STOCHASTIC PROCESSES 4 Chaper 4 For an advanced analyss of he properes of he Wener process, see: Revus D and Yor M: Connuous marngales and Brownan Moon Karazas I and Shreve S E:
More informationChapter 6: AC Circuits
Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.
More informationOrdinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s
Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class
More informationDEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL
DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL Sco Wsdom, John Hershey 2, Jonahan Le Roux 2, and Shnj Waanabe 2 Deparmen o Elecrcal Engneerng, Unversy o Washngon, Seale, WA, USA
More informationTesting a new idea to solve the P = NP problem with mathematical induction
Tesng a new dea o solve he P = NP problem wh mahemacal nducon Bacground P and NP are wo classes (ses) of languages n Compuer Scence An open problem s wheher P = NP Ths paper ess a new dea o compare he
More informationDisplacement, Velocity, and Acceleration. (WHERE and WHEN?)
Dsplacemen, Velocy, and Acceleraon (WHERE and WHEN?) Mah resources Append A n your book! Symbols and meanng Algebra Geomery (olumes, ec.) Trgonomery Append A Logarhms Remnder You wll do well n hs class
More informationFI 3103 Quantum Physics
/9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon
More informationGENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim
Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran
More information. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.
Lnear Algebra Lecure # Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons
More informationF-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction
ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or
More information. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.
Mah E-b Lecure #0 Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons are
More informationExistence and Uniqueness Results for Random Impulsive Integro-Differential Equation
Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal
More informationMotion of Wavepackets in Non-Hermitian. Quantum Mechanics
Moon of Wavepaces n Non-Herman Quanum Mechancs Nmrod Moseyev Deparmen of Chemsry and Mnerva Cener for Non-lnear Physcs of Complex Sysems, Technon-Israel Insue of Technology www.echnon echnon.ac..ac.l\~nmrod
More informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system Y X (2) and write EOM (1) as two first-order Eqs.
Handou # 6 (MEEN 67) Numercal Inegraon o Fnd Tme Response of SDOF mechancal sysem Sae Space Mehod The EOM for a lnear sysem s M X DX K X F() () X X X X V wh nal condons, a 0 0 ; 0 Defne he followng varables,
More informationPerformance Analysis for a Network having Standby Redundant Unit with Waiting in Repair
TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen
More informationThe ray paths and travel times for multiple layers can be computed using ray-tracing, as demonstrated in Lab 3.
C. Trael me cures for mulple reflecors The ray pahs ad rael mes for mulple layers ca be compued usg ray-racg, as demosraed Lab. MATLAB scrp reflec_layers_.m performs smple ray racg. (m) ref(ms) ref(ms)
More information10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :
. A. IUITS Synopss : GOWTH OF UNT IN IUIT : d. When swch S s closed a =; = d. A me, curren = e 3. The consan / has dmensons of me and s called he nducve me consan ( τ ) of he crcu. 4. = τ; =.63, n one
More informationSpurious oscillations and conservation errors in interface-capturing schemes
Cener for Turbulence Research Annual Research Brefs 8 115 Spurous oscllaons and conservaon errors n nerface-capurng schemes By E. Johnsen Movaon and objecves When shock-capurng schemes are appled o flows
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours
NATONAL UNVERSTY OF SNGAPORE PC5 ADVANCED STATSTCAL MECHANCS (Semeser : AY 1-13) Tme Allowed: Hours NSTRUCTONS TO CANDDATES 1. Ths examnaon paper conans 5 quesons and comprses 4 prned pages.. Answer all
More informationAppendix H: Rarefaction and extrapolation of Hill numbers for incidence data
Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs
More informationA Paper presentation on. Department of Hydrology, Indian Institute of Technology, Roorkee
A Paper presenaon on EXPERIMENTAL INVESTIGATION OF RAINFALL RUNOFF PROCESS by Ank Cakravar M.K.Jan Kapl Rola Deparmen of Hydrology, Indan Insue of Tecnology, Roorkee-247667 Inroducon Ranfall-runoff processes
More informationCS286.2 Lecture 14: Quantum de Finetti Theorems II
CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2
More informationRobust and Accurate Cancer Classification with Gene Expression Profiling
Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem
More informationGORDON AND NEWELL QUEUEING NETWORKS AND COPULAS
Yugoslav Journal of Operaons Research Vol 9 (009) Number 0- DOI:0.98/YUJOR0900C GORDON AND NEWELL QUEUEING NETWORKS AND COPULAS Danel CIUIU Facul of Cvl Indusral and Agrculural Buldngs Techncal Unvers
More informationComb Filters. Comb Filters
The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of
More informationImplementation of Quantized State Systems in MATLAB/Simulink
SNE T ECHNICAL N OTE Implemenaon of Quanzed Sae Sysems n MATLAB/Smulnk Parck Grabher, Mahas Rößler 2*, Bernhard Henzl 3 Ins. of Analyss and Scenfc Compung, Venna Unversy of Technology, Wedner Haupsraße
More informationVEHICLE DYNAMIC MODELING & SIMULATION: COMPARING A FINITE- ELEMENT SOLUTION TO A MULTI-BODY DYNAMIC SOLUTION
21 NDIA GROUND VEHICLE SYSTEMS ENGINEERING AND TECHNOLOGY SYMPOSIUM MODELING & SIMULATION, TESTING AND VALIDATION (MSTV) MINI-SYMPOSIUM AUGUST 17-19 DEARBORN, MICHIGAN VEHICLE DYNAMIC MODELING & SIMULATION:
More informationEVALUATION OF FORCE COEFFICIENTS FOR A 2-D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL
The Sevenh Asa-Pacfc Conference on Wnd Engneerng, November 8-, 009, Tape, Tawan EVALUATION OF FORCE COEFFICIENTS FOR A -D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL S. Chra Ganapah, P. Harkrshna,
More informationAnisotropic Behaviors and Its Application on Sheet Metal Stamping Processes
Ansoropc Behavors and Is Applcaon on Shee Meal Sampng Processes Welong Hu ETA-Engneerng Technology Assocaes, Inc. 33 E. Maple oad, Sue 00 Troy, MI 48083 USA 48-79-300 whu@ea.com Jeanne He ETA-Engneerng
More informationOn computing differential transform of nonlinear non-autonomous functions and its applications
On compung dfferenal ransform of nonlnear non-auonomous funcons and s applcaons Essam. R. El-Zahar, and Abdelhalm Ebad Deparmen of Mahemacs, Faculy of Scences and Humanes, Prnce Saam Bn Abdulazz Unversy,
More informationComparison of Differences between Power Means 1
In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,
More informationTrack Properities of Normal Chain
In. J. Conemp. Mah. Scences, Vol. 8, 213, no. 4, 163-171 HIKARI Ld, www.m-har.com rac Propes of Normal Chan L Chen School of Mahemacs and Sascs, Zhengzhou Normal Unversy Zhengzhou Cy, Hennan Provnce, 4544,
More informationFall 2010 Graduate Course on Dynamic Learning
Fall 200 Graduae Course on Dynamc Learnng Chaper 4: Parcle Flers Sepember 27, 200 Byoung-Tak Zhang School of Compuer Scence and Engneerng & Cognve Scence and Bran Scence Programs Seoul aonal Unversy hp://b.snu.ac.kr/~bzhang/
More informationV.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS
R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon
More informationCoupled Method for Solving Time-Fractional Navier-Stokes Equation
INTERNATIONAL JOURNAL OF CIRCUITS, SYSTEMS AND SIGNAL PROCESSING Volume, 8 Coupled Mehod for Solng Tme-Fraconal Naer-Sokes Equaon S. O. Edek, and G. O. Aknlab Absrac Ths paper wnesses he couplng of wo
More informationOnline Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading
Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng
More informationIncluding the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.
Soluons o Ordnary Derenal Equaons An ordnary derenal equaon has only one ndependen varable. A sysem o ordnary derenal equaons consss o several derenal equaons each wh he same ndependen varable. An eample
More informationLecture 2 M/G/1 queues. M/G/1-queue
Lecure M/G/ queues M/G/-queue Posson arrval process Arbrary servce me dsrbuon Sngle server To deermne he sae of he sysem a me, we mus now The number of cusomers n he sysems N() Tme ha he cusomer currenly
More informationDepartment of Economics University of Toronto
Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of
More informationFirst-order piecewise-linear dynamic circuits
Frs-order pecewse-lnear dynamc crcus. Fndng he soluon We wll sudy rs-order dynamc crcus composed o a nonlnear resse one-por, ermnaed eher by a lnear capacor or a lnear nducor (see Fg.. Nonlnear resse one-por
More informationTime-interval analysis of β decay. V. Horvat and J. C. Hardy
Tme-nerval analyss of β decay V. Horva and J. C. Hardy Work on he even analyss of β decay [1] connued and resuled n he developmen of a novel mehod of bea-decay me-nerval analyss ha produces hghly accurae
More informationTSS = SST + SSE An orthogonal partition of the total SS
ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
Ths documen s downloaded from DR-NTU, Nanyang Technologcal Unversy Lbrary, Sngapore. Tle A smplfed verb machng algorhm for word paron n vsual speech processng( Acceped verson ) Auhor(s) Foo, Say We; Yong,
More informationTransient Numerical of Piston Wind in Subway Station. Haitao Bao
Appled Mechancs and Maerals Submed: 2014-07-20 ISSN: 1662-7482, Vols. 644-650, pp 467-470 Acceped: 2014-07-21 do:10.4028/www.scenfc.ne/amm.644-650.467 Onlne: 2014-09-22 2014 Trans Tech Publcaons, Swzerland
More informationABSTRACT KEYWORDS. Bonus-malus systems, frequency component, severity component. 1. INTRODUCTION
EERAIED BU-MAU YTEM ITH A FREQUECY AD A EVERITY CMET A IDIVIDUA BAI I AUTMBIE IURACE* BY RAHIM MAHMUDVAD AD HEI HAAI ABTRACT Frangos and Vronos (2001) proposed an opmal bonus-malus sysems wh a frequency
More informationLearning Objectives. Self Organization Map. Hamming Distance(1/5) Introduction. Hamming Distance(3/5) Hamming Distance(2/5) 15/04/2015
/4/ Learnng Objecves Self Organzaon Map Learnng whou Exaples. Inroducon. MAXNET 3. Cluserng 4. Feaure Map. Self-organzng Feaure Map 6. Concluson 38 Inroducon. Learnng whou exaples. Daa are npu o he syse
More informationUNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION
INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he
More informationby Lauren DeDieu Advisor: George Chen
b Laren DeDe Advsor: George Chen Are one of he mos powerfl mehods o nmercall solve me dependen paral dfferenal eqaons PDE wh some knd of snglar shock waves & blow-p problems. Fed nmber of mesh pons Moves
More informationShould Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth
Should Exac Index umbers have Sandard Errors? Theory and Applcaon o Asan Growh Rober C. Feensra Marshall B. Rensdorf ovember 003 Proof of Proposon APPEDIX () Frs, we wll derve he convenonal Sao-Vara prce
More informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationDensity Matrix Description of NMR BCMB/CHEM 8190
Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon If we say wh one bass se, properes vary only because of changes n he coeffcens weghng each bass se funcon x = h< Ix > - hs s how we calculae
More informationNormal Random Variable and its discriminant functions
Noral Rando Varable and s dscrnan funcons Oulne Noral Rando Varable Properes Dscrnan funcons Why Noral Rando Varables? Analycally racable Works well when observaon coes for a corruped snle prooype 3 The
More informationAnisotropy and oblique total transmission at a planar negative-index interface
Ansoropy and oblque oal ransmsson a a planar negave-ndex nerface Le Zhou, C.T. Chan and P. Sheng Deparmen of Physcs, The Hong Kong Unversy of Scence and Technology Clear Waer Bay, Kowloon, Hong Kong, Chna
More informationM. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria
IOSR Journal of Mahemacs (IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume 0, Issue 4 Ver. IV (Jul-Aug. 04, PP 40-44 Mulple SolonSoluons for a (+-dmensonalhroa-sasuma shallow waer wave equaon UsngPanlevé-Bӓclund
More information