Implementing a Convolutional Perfectly Matched Layer in a finite-difference code for the simulation of seismic wave propagation in a 3D elastic medium
|
|
- Gilbert Hutchinson
- 5 years ago
- Views:
Transcription
1 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
2
3 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 Su carre ou a par of Reearch ace - BRGM 7 PDR7ARN4 Duceller A., Aoch H. Checke b: Name: Phlppe Joue Dae: Sgnaure: Approe b: Name: Hormo Moare Dae: Sgnaure: BRGM' qual managemen em cerfe ISO 9: b AFAQ IM 3 ANG Aprl 5
4 Kewor: Semc wae propagaon, Aborbng conon, Perfecl Mache Laer, Conoluonal Perfecl Mache Laer, Non-reflecng conon, Fne-fference. In bblograph, h repor houl be ce a follow: Duceller A., Aoch H. (7) Implemenng a Conoluonal Perfecl Mache Laer n a fne-fference coe for he mulaon of emc wae propagaon n a 3D elac meum, BRGM echncal repor, BRGM/RP-559-FR. BRGM, 5. No par of h ocumen ma be reprouce whou he pror permon of BRGM.
5 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Snop One of he mo popular meho o mulae numercall he emc wae propagaon n an elac mea he fne-fference meho. In he cone of numercal moellng n unboune mea, he wae nee o be aborbe a he arfcal bounare of he compuaonal oman an herefore, necear o efne non-reflecng conon a hee bounare o mmc an unboune meum. The Perfecl Mache Laer (PML) conon ha he remarkable proper of hang a ero reflecon coeffcen for all angle of ncence an all frequence before creaon an ha become wel ue (e.g. Collno an Togka, ). Howeer, h reflecon coeffcen no ero anmore afer he creaon an become een er large a grang ncence. Therefore, an mproe eron of he PML conon ha been eelope: he Conoluonal Perfecl Mache Laer (CPML) conon (e.g. Komach an Marn, 7). The am of h repor o ecrbe how are mplemene hee wo non-reflecng conon n a 3D fne-fference compuaonal coe for an elac meum. We wll fr eplan he heorecal formulaon of he PML an CPML conon. Then, we wll epoe he rucure of our compuaonal coe. Fnall, we wll preen he reul of ome numercal e performe o erf he effcenc of he meho. BRGM/RP-559-FR Progre repor 3
6
7 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Conen. Inroucon...7. The CPML formulaon CLASSICAL PML FORMULATION CPML FORMULATION Implemenaon n he 3D coe GRID CELL BOUNDARY CONDITIONS SOURCE FONCTION STRUCTURE OF THE CODE Numercal e PML CPML Concluon Reference...35 BRGM/RP-559-FR Progre repor 5
8 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon L of lluraon Fgure Elemenar gr cell of he hree-menonal aggere paal fnefference meho... 9 Fgure - Source me funcon mpoe a pon ource... 3 Fgure 3 Dampng profle ()... 4 Fgure 4 Tme eoluon of he, an componen of he hree-menonal eloc ecor a he fr aon of he numercal oluon wh clacal PML bounare for he hn lce an he large lce Fgure 5 Tme eoluon of he, an componen of he hree-menonal eloc ecor a he econ aon of he numercal oluon wh clacal PML bounare for he hn lce an he large lce Fgure 6 Tme eoluon of he, an componen of he hree-menonal eloc ecor a he hr aon of he numercal oluon wh clacal PML bounare for he hn lce an he large lce Fgure 7 Snapho of he, an componen of he hree-menonal eloc ecor of he numercal oluon wh clacal PML bounare for he hn lce, a me ep 375, 65, 875, 5, 375 an 65, a 39 m... 8 Fgure 8 Tme eoluon of he, an componen of he hree-menonal eloc ecor a he fr aon of he numercal oluon wh CPML bounare for he hn lce an he large lce Fgure 9 Tme eoluon of he, an componen of he hree-menonal eloc ecor a he econ aon of he numercal oluon wh CPML bounare for he hn lce an he large lce Fgure Tme eoluon of he, an componen of he hree-menonal eloc ecor a he hr aon of he numercal oluon wh CPML bounare for he hn lce an he large lce Fgure Snapho of he, an componen of he hree-menonal eloc ecor of he numercal oluon wh CPML bounare for he hn lce, a me ep 375, 65, 875, 5, 375 an 65, a 39 m BRGM/RP-559-FR Progre repor
9 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon. Inroucon One of he mo popular meho o mulae numercall he emc wae propagaon n an elac mea he fne-fference meho. In he cone of numercal moellng n unboune mea, he wae nee o be aborbe a he arfcal bounare of he compuaonal oman an herefore, necear o efne non-reflecng conon a hee bounare o mmc an unboune meum. The Perfecl Mache Laer (PML) conon ha he remarkable proper of hang a ero reflecon coeffcen for all angle of ncence an all frequence before creaon an ha become wel ue (e.g. Collno an Togka, ). Howeer, h reflecon coeffcen no ero anmore afer he creaon an become een er large a grang ncence. Therefore, an mproe eron of he PML conon ha been eelope: he Conoluonal Perfecl Mache Laer (CPML) conon (e.g. Komach an Marn, 7). The am of h repor o ecrbe how are mplemene hee wo non-reflecng conon n a 3D fne-fference compuaonal coe for an elac meum. We wll fr eplan he heorecal formulaon of he PML an CPML conon. Then, we wll epoe he rucure of our compuaonal coe. Fnall, we wll preen he reul of ome numercal e performe o erf he effcenc of he meho. BRGM/RP-559-FR Progre repor 7
10
11 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon. The CPML formulaon One of he mo popular meho o mulae numercall he emc wae propagaon n an elac mea he fne-fference meho (Grae, 996, Leaner, 988). In he cone of numercal moellng n unboune mea, he wae nee o be aborbe a he arfcal bounare of he compuaonal oman an herefore, necear o efne non-reflecng conon a hee bounare o mmc an unboune meum. The Perfecl Mache Laer (PML) conon ha he remarkable proper of hang a ero reflecon coeffcen for all angle of ncence an all frequence before creaon an ha become wel ue (e.g. Collno an Togka, ). Howeer, h reflecon coeffcen no ero anmore afer he creaon an become een er large a grang ncence. Therefore, an mproe eron of he PML conon ha been eelope: he Conoluonal Perfecl Mache Laer (CPML) conon (Komach an Marn, 7). We wll fr eplan he clacal PML formulaon, an hen we wll nrouce he CPML formulaon... CLASSICAL PML FORMULATION The fferenal form of he 3D elac wae equaon can be wren a: ρ ρ ρ for he eloce (, an ) an: µ BRGM/RP-559-FR Progre repor 9
12 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon BRGM/RP-559-FR Progre repor µ µ µ µ µ for he re (,,,, an ), where ρ en an an µ are elac coeffcen. In he parcular cae of a homogenou meum, h equaon ha plane wae oluon of he form ( ) ( ) k A r r r. ep. Le u coner a regular oman locae n. In orer o ao reflecon a, we nee o efne an aborbng laer locae n >. Th laer wll be an aborbng laer f he amplue of he wae ecreae eponenall n he recon for >. To efne h laer, we wll efne a new wae equaon, uch ha oluon a plane wae n he regular oman an an eponenall-ecang wae n he PML. We nrouce hen a new comple coornae ~. Thu, he PML can be ewe a an analcal connuaon of he real coornae n he cople pace (Collno an Togka, ). The coornae ~ efne b: ( ) ( ) ~ Equaon Equaon ge u upon fferenang: ~ ϖ Equaon where
13 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon BRGM/RP-559-FR Progre repor We wll now change he orgnal wae equaon wren n erm of arable, an no a new wae equaon wren n erm of arable ~, an. We fr rewre he elac wae equaon n he frequenc oman: ρ ρ ρ for he eloce an: µ µ µ µ µ µ for he re. We hen replace he wae equaon wren n erm of wh a generale wae equaon wren n erm of ~. Ine he regular oman, boh equaon are encal becaue.
14 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon BRGM/RP-559-FR Progre repor ~ ρ ~ ρ ~ ρ for he eloce an: ~ ~ µ µ ~ µ ~ ~ µ ~ µ µ for he re. In he PML, he new wae equaon ha eponenall-ecang plane wae oluon of he form: ( ) ( ) ( ) ( ) ( ) k k A k k k A ep. ep ~ ep r r r r n he recon.
15 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon BRGM/RP-559-FR Progre repor 3 Ung Equaon : ~, we rewre he wae equaon n erm of raher han ~ (n he regular oman, ): ρ ρ ρ for he eloce an: µ µ µ µ µ µ for he re. The eloc an re fel are hen pl no wo par, an, uch ha:
16 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 4 BRGM/RP-559-FR Progre repor ρ ρ ρ ϖρ ρ ϖρ for he eloce an: µ µ µ
17 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon BRGM/RP-559-FR Progre repor 5 µ µ ϖ µ µ ϖ µ for he re. Conerng back o he me oman, we fnall ge: ρ ρ ρ ρ ρ ρ for he eloce an:
18 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 6 BRGM/RP-559-FR Progre repor µ µ µ µ µ µ µ µ for he re. Th new wae equaon perm o hae an eponenall-ecang plane wae oluon n he PML laer ( > ). Therefore he bounar a appear o be a non-
19 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon reflecng bounar. We hen ue he ame echnque for he fe oher bounare of he regular oman. In each ege of he oman, we appl h echnque n wo recon. In each corner, we appl h echnque n he hree recon. Th formulaon perm o hae numercal oluon wh reuce arfcal reflecon bu effcenc become poor a grang ncence. The CPML formulaon wa eelope o are h ue... CPML FORMULATION The man ea of he CPML formulaon o chooe a more general efnon of b nroucng new arable α en κ : κ α Equaon 3 A h epreon epen on frequenc, when we go back n he me oman, we wll ge a me conoluon on each mofe paal erae. Denong () he nere Fourer ranform of ( ) *, we can wre: ~ Equaon 4 * beng a conoluon operaor. We hae from Equaon 3: ( ) δ ) H ep κ κ ( ( ) ( ( κ α ) ) Equaon 5 where δ() an H() are he Drac an Heae rbuon. If we enoe: ζ κ ( ) H ( ) ep( ( κ α ) ) Equaon 6 we can ee from Equaon 4, 5 an 6 ha: ~ ζ κ ( ) * BRGM/RP-559-FR Progre repor 7
20 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 8 BRGM/RP-559-FR Progre repor Le u enoe n ϕ he conoluon erm compue a me ep n, wh beng he me ep: ( ) ( ) τ τ ζ ζ ϕ τ n n n * Equaon 7 Equaon 7 ge u b creaon: ( ) ( ) ( ) ( ) ( ) n m m n n m m m n n m n n n m m m n m Z ϕ τ τ ζ ϕ τ τ ζ ϕ τ Equaon 8 wh: ( ) ( ) ( ) m m m Z τ τ ζ Equaon 9 From Equaon 6 an 9, we oban: ( ) ( ) ( ) ( ) ( ) ( ) m a m Z m m α κ τ τ α κ κ ep ep, wh ( ) ( ) b a α κ κ an ( ) ( ) b α κ ep We coner hen φ a a memor arable whoe me eoluon goerne a each ep b: n n n a b ϕ ϕ From a praccal pon of ew, o mplemen he CPML echnque n a compuaonal coe, we ju nee o replace he paal erae b: ϕ κ ~ an upae φ a each me ep.
21 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 3. Implemenaon n he 3D coe 3.. GRID CELL We ue a clacal econ-orer aggere gr n pace an me (Maaraga, 976, an Vreu, 986): Fgure Elemenar gr cell of he hree-menonal aggere paal fne-fference meho of Maaraga (976) an Vreu (986) ue clacall o cree he equaon of elaonamc. 3.. BOUNDARY CONDITIONS We pu non-reflecng bounare on he e of he compuaonal oman (Komach an Marn, 7) o a o mulae wae propagaon n nfne elac meum. In h cae, he bounar conon are: BRGM/RP-559-FR Progre repor 9
22 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon (,.,. ) ( mn ma,.,.) (., jmn,. ) (., jma,. ) (.,., kmn ) (.,., k ) ma an (,.,. ) ( mn ma,.,.) (., jmn,. ) (., jma,. ) (.,., kmn ) (.,., k ) ma an (,.,. ) ( mn ma,.,.) (., jmn,. ) (., jma,. ) (.,., kmn ) (.,., k ) ma mn, ma, j mn, j ma, k mn an k ma beng he lm of he regular oman an beng he number of cell n he (C)PML laer SOURCE FONCTION The ource gen b a eloc ecor n he plane k k ource, he ame a Komach an Marn (7) for alaon. A he pon ource, we ju a: ( ) ρ ( α ) a ( ) ep a ( ) n f * Equaon o ( ource,j ource,k ource ) an: ( α ) a ( ) ep a ( ) ( ) ρ co f * Equaon o ( ource,j ource,k ource ), beng he me, α, f, a an beng ource parameer gen n he coe an ρ beng he en STRUCTURE OF THE CODE The rucure of he compuaonal coe he followng: Inalaon BRGM/RP-559-FR Progre repor
23 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Calculaon of he alue of (), α (), κ (), (j), α (j), κ (j), (k), α (k) an κ (k) Begnnng of he me loop Tme (n/) (*) for XMIN-,,XMAX { for j YMIN-,,YMAX { for k ZMIN-,,ZMAX { f (,j,k, n he CPML) { Calculaon of he j n he CPML wh n orer FD } ele { Calculaon of he j n he regular oman wh n orer FD } } } } Tme n for XMIN-,,XMAX { for j YMIN-,,YMAX { for k ZMIN-,,ZMAX { f (,j,k n he CPML) { Calculaon of he n he CPML wh n orer FD } ele { Calculaon of he n he regular oman wh n orer FD } } BRGM/RP-559-FR Progre repor
24 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon } } Ang a eloc ecor a he ource Impong he bounar conon Sang he alue n fle for he emogram an he napho Incremen of n an reurn o (*) BRGM/RP-559-FR Progre repor
25 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 4. Numercal e 4.. PML To e our compuaonal coe, we coner he ame moel a Komach an Marn (7). We coner wo moel: he e of he fr moel *64*64 meer, repreenng a oman much longer han we (.e., a hn lce). In h moel, he wae wll propagae a grang ncence along he ege, whch ma caue purou reflecon. To hghlgh hee reflecon, we wll coner a econ, larger moel (6*64*64 meer), whch bounare are far enough from he ource o ao an purou reflecon, an we wll compare he reul of he wo moel. The hn moel cree ung a gr comprng pon * 64 pon * 64 pon. The large moel ecree ung a gr comprng 6 * 64 * 64 pon. The e of a gr cell m. The compreonal wae pee Vp 33 m. -, he hear wae pee V Vp/ m. - an he en ρ 8 kg.m -3. The me ep.6 an we perform he mulaon for 5 me ep,.e. four econ. The pon ource a eloc ecor locae n ource 79 m, ource 47 m an ource 39 m. The me funcon gen b Equaon an. The parameer for he ource are: α 35, f 7, a π*f,. / f an f 7 H Source funcon 5 ource (m/) 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 () Fgure - Source me funcon mpoe a pon ource The aon where we compue emogram are locae a he followng pon: BRGM/RP-559-FR Progre repor 3
26 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Saon number X coornae (m) Y coornae (m) Z coornae (m) Saon locae a he mle of he moel where here no purou reflecon. Saon an 3 are locae along he ege a grang an er grang ncence, where purou reflecon coul or he emogram. For he ampng coeffcen n he PML, we ue he followng profle: ( ) L N where he ance beween he pon where compue an he begnnng of he PML laer, an L he hckne of he PML: L *, N an: ( N ) log( R ) p c wh Rc.. L Dampng profle () 6 4 () (m) Fgure 3 Dampng profle (). In he regular oman ( < ), he ampng profle null an he wae equaon no mofe. In he aborbng bounar ( > ), he ampng profle ncreae wh an he wae equaon mofe n orer o hae eponenall-ecang plane wae oluon. Aborbng laer of hckne gr pon are mplemene on he e of he moel. 4 BRGM/RP-559-FR Progre repor
27 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon The reul are hown hereafer n Fgure 4 o 7. Ob - V - FD Ob - V - FD,, Thn lce Large lce,,,8 Thn lce Large lce,6,4 V -, V, -, -, -,4 -,3 -,6,,4,63,84,5,6,47,68,89,,3,5,7,93 3,4 3,35 3,56 3,77 3,98 me,,4,63,84,5,6,47,68,89,,3,5,7,93 3,4 3,35 3,56 3,77 3,98 me Ob - V - FD,,8,6,4 V, -, -,4 -,6 Thn lce Large lce,,44,67,89,,33,56,78,,45,67,89 3, 3,34 3,56 3,78 me Fgure 4 Tme eoluon of he (lef), (rgh) an (boom) componen of he hreemenonal eloc ecor a he fr aon of he numercal oluon wh clacal PML bounare for he hn lce (blue) an he large lce (pnk). A he fr aon, relael far from he begnnng of he PML laer an wh non-grang ncence, he agreemen almo perfec for an. For, here numercal peron an he eloc no eacl null. BRGM/RP-559-FR Progre repor 5
28 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Ob - V - FD Ob - V - FD,, Thn lce Large lce,,8,6 Thn lce Large lce,4 V V, -, -, -, -,4 -,3 -,6,,4,63,84,5,6,47,68,89,,3,5,7,93 3,4 3,35 3,56 3,77 3,98 me,,43,64,85,6,8,49,7,9,3,34,55,77,98 3,9 3,4 3,6 3,83 me Ob - V - FD,, V -, -, -,3 -,4 Thn lce Large lce,,44,66,88,,3,53,75,97,9,4,63,85 3,7 3,9 3,5 3,73 3,95 me Fgure 5 Tme eoluon of he (lef), (rgh) an (boom) componen of he hreemenonal eloc ecor a he econ aon of he numercal oluon wh clacal PML bounare for he hn lce (blue) an he large lce (pnk). A he econ aon, purou ocllaon, n parcular for he S-wae of he componen, ar o appear for an. For, here numercal peron an he eloc no eacl null. 6 BRGM/RP-559-FR Progre repor
29 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Ob 3 - V - FD Ob 3 - V - FD,5 Thn lce,5 Large lce, Large lce,4 Large lce,3,5, V V, -,5 -, -, -, -,5 -,3,,43,64,85,6,8,49,7,9,3,34,55,77,98 3,9 3,4 3,6 3,83 me,,43,64,85,6,8,49,7,9,3,34,55,77,98 3,9 3,4 3,6 3,83 me Ob 3 - V - FD,3,, V -, -, -,3 -,4 Thn lce Large lce,,44,66,88,,3,53,75,97,9,4,63,85 3,7 3,9 3,5 3,73 3,95 me Fgure 6 Tme eoluon of he (lef), (rgh) an (boom) componen of he hreemenonal eloc ecor a he hr aon of he numercal oluon wh clacal PML bounare for he hn lce (blue) an he large lce (pnk). A he hr aon, he ocllaon become large, he P-wae no correcl calculae an he hape of he S-wae compleel ore for an. For, here numercal peron an he eloc no eacl null. BRGM/RP-559-FR Progre repor 7
30 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Fgure 7 Snapho of he (lef), (mle) an (rgh) componen of he hreemenonal eloc ecor of he numercal oluon wh clacal PML bounare for he hn lce, a me ep 375, 65, 875, 5, 375 an 65 (op o boom), a 39 m. The hree o ncae he poon of he aon a whch he emogram repreene n Fgure 4 o 6 are recore. Spurou wae appear a grang ncence along he ege of he moel an en purou energ back no he man oman. We fn ha ome large reflecon wae reman epecall a aon an 3. In fac, he napho how ome rong reflecon a he bounar cloe o he ource n he begnnng. We alo obere he ocllaon ocke n he PML laer, whch propagae parallel o he bounar. Th he problem of he PML preoul pone ou (e.g. Komach an Marn, 7). 8 BRGM/RP-559-FR Progre repor
31 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 4.. CPML To e he CPML, we ue he ame approach a for he PML. We chooe o make α ar n a lnear fahon beween mamal alue α ma a he begnnng of he CPML an ero a op. We ake α ma π*f. We chooe κ conan, equal o. The reul are hown hereafer n Fgure 8 o. Ob - V - FD Ob - V - FD,, Thn lce Large lce,,,8 Thn lce Large lce,6,4 V -, V, -, -, -,4 -,3 -,6,,4,63,84,5,6,47,68,89,,3,5,7,93 3,4 3,35 3,56 3,77 3,98 me,,4,63,84,5,6,47,68,89,,3,5,7,93 3,4 3,35 3,56 3,77 3,98 me Ob - V - FD,,8,6,4 V, -, -,4 -,6 Thn lce Large lce,,44,67,89,,33,56,78,,45,67,89 3, 3,34 3,56 3,78 me Fgure 8 Tme eoluon of he (lef), (rgh) an (boom) componen of he hreemenonal eloc ecor a he fr aon of he numercal oluon wh CPML bounare for he hn lce (blue) an he large lce (pnk). A he fr aon, relael far from he begnnng of he CPML laer an wh non-grang ncence, he agreemen almo perfec for an. For, here numercal peron an he eloc no eacl null. BRGM/RP-559-FR Progre repor 9
32 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Ob - V - FD Ob - V - FD,, Thn lce Large lce,,8,6 Thn lce Large lce,4 V V, -, -, -, -,4 -,3 -,6,,4,63,84,5,6,47,68,89,,3,5,7,93 3,4 3,35 3,56 3,77 3,98 me,,43,64,85,6,8,49,7,9,3,34,55,77,98 3,9 3,4 3,6 3,83 me Ob - V - FD,, V -, -, -,3 -,4 Thn lce Large lce,,44,66,88,,3,53,75,97,9,4,63,85 3,7 3,9 3,5 3,73 3,95 me Fgure 9 Tme eoluon of he (lef), (rgh) an (boom) componen of he hreemenonal eloc ecor a he econ aon of he numercal oluon wh CPML bounare for he hn lce (blue) an he large lce (pnk). A he econ aon, wh grang ncence an raher cloe o he begnnng of he CPML, he agreemen reman ecellen for an. For, here numercal peron an he eloc no eacl null. 3 BRGM/RP-559-FR Progre repor
33 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Ob 3 - V - FD Ob 3 - V - FD,5, Thn lce Large lce,5,4 Thn lce Large lce,3,5, V V, -,5 -, -, -, -,5 -,3,,43,64,85,6,8,49,7,9,3,34,55,77,98 3,9 3,4 3,6 3,83 me,,43,64,85,6,8,49,7,9,3,34,55,77,98 3,9 3,4 3,6 3,83 me Ob 3 - V - FD,3,, V -, -, -,3 -,4 Thn lce Large lce,,44,66,88,,3,53,75,97,9,4,63,85 3,7 3,9 3,5 3,73 3,95 me Fgure Tme eoluon of he (lef), (rgh) an (boom) componen of he hreemenonal eloc ecor a he hr aon of he numercal oluon wh CPML bounare for he hn lce (blue) an he large lce (pnk). A he hr aon, n he ffcul cae of er grang ncence, of a long ance of propagaon an of a aon locae cloe o he begnnng of he CPML laer, he agreemen reman er afacor for an. For, here numercal peron an he eloc no eacl null. BRGM/RP-559-FR Progre repor 3
34 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon Fgure Snapho of he (lef), (mle) an (rgh) componen of he hreemenonal eloc ecor of he numercal oluon wh CPML bounare for he hn lce, a me ep 375, 65, 875, 5, 375 an 65 (op o boom) a 39 m. The hree o ncae he poon of he aon a whch he emogram repreene n Fgure 8 o are recore. No purou wae of gnfcan amplue ble, een a grang ncence. In he CPML cae, we are able o ge he encal reul a Komach an Marn (7) well mproe n erm of elmnaon of purou reflecon wh repec o he preou one (PML). We o no hae an more he arfcal wae reflecon een n he napho. We confrm here ha he CPML aborbng conon work much beer. 3 BRGM/RP-559-FR Progre repor
35 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 5. Concluon The behaour of he Perfecl mache Laer a grang ncence well mproe wh he new Conoluonnal PML approach. I can be ueful when ung wae propagaon n hn lce or when he ource locae near he ege of he moel. The co of CPML n memor orage mlar o he clacal PML. BRGM/RP-559-FR Progre repor 33
36
37 Implemenng a CPML n a 3D fne-fference coe for he mulaon of emc wae propagaon 6. Reference Collno F., Togka C. () - Applcaon of he PML aborbng laer moel o he lnear elaonamc problem n anoropc heerogeneou mea. Geophc, 66(), Grae R.W. (996) - Smulang Semc Wae propagaon n 3D Elac Mea Ung Saggere-Gr Fne Dfference. Bull. Sem. Soc. Am., 86, 9-6. Komach D., Marn R. (7) - An unpl conoluonal Perfecl Mache Laer mproe a grang ncence for he emc wae equaon. Geophc, n pre, -. Leaner A.R. (988) Fourh-orer fne-fference P-SV emogram. Geophc, 53, Maaraga R. (976) - Dnamc of an epanng crcular faul. Bull. Sem. Soc. Am., 65, Vreu J. (986) - P-SV wae propagaon n heerogeneou mea: eloc-re fne-fference meho. Geophc, 5, BRGM/RP-559-FR Progre repor 35
38
39
40 Scenfc an Techncal Cenre Deelopmen Plannng an Naural Rk Don 3, aenue Claue-Gullemn - BP Orléan Cee France Tel.: 33 ()
Implementing a Convolutional Perfectly Matched Layer in a finite-difference code for the simulation of seismic wave propagation in a 2D elastic medium
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 Implemenng a Conoluonal Perfecl Mached
More informationSound decay in a rectangular room with specular and diffuse reflecting surfaces
Soun ecay n a recangular room wh pecular an ffue reflecng urface Nkolay Kanev Anreyev Acouc Inue, Mocow, Rua Summary Soun ecay n room uually efne by aborpon an caerng propere of urface In ome cae oun ecay
More informationELEC 201 Electric Circuit Analysis I Lecture 9(a) RLC Circuits: Introduction
//6 All le courey of Dr. Gregory J. Mazzaro EE Elecrc rcu Analy I ecure 9(a) rcu: Inroucon THE ITADE, THE MIITAY OEGE OF SOUTH AOINA 7 Moulre Sree, harleon, S 949 V Sere rcu: Analog Dcoery _ 5 Ω pf eq
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 informationA. Inventory model. Why are we interested in it? What do we really study in such cases.
Some general yem model.. Inenory model. Why are we nereed n? Wha do we really udy n uch cae. General raegy of machng wo dmlar procee, ay, machng a fa proce wh a low one. We need an nenory or a buffer or
More informationRandomized Perfect Bipartite Matching
Inenive Algorihm Lecure 24 Randomized Perfec Biparie Maching Lecurer: Daniel A. Spielman April 9, 208 24. Inroducion We explain a randomized algorihm by Ahih Goel, Michael Kapralov and Sanjeev Khanna for
More informationFX-IR Hybrids Modeling
FX-IR Hybr Moeln Yauum Oajma Mubh UFJ Secure Dervave Reearch Dep. Reearch & Developmen Dvon Senor Manaer oajma-yauum@c.mu.jp Oaka Unvery Workhop December 5 h preenaon repreen he vew o he auhor an oe no
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 information(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function
MACROECONOMIC THEORY T J KEHOE ECON 87 SPRING 5 PROBLEM SET # Conder an overlappng generaon economy le ha n queon 5 on problem e n whch conumer lve for perod The uly funcon of he conumer born n perod,
More informationESS 265 Spring Quarter 2005 Kinetic Simulations
SS 65 Spng Quae 5 Knec Sulaon Lecue une 9 5 An aple of an lecoagnec Pacle Code A an eaple of a knec ulaon we wll ue a one denonal elecoagnec ulaon code called KMPO deeloped b Yohhau Oua and Hoh Mauoo.
More informationAlgorithms and Data Structures 2011/12 Week 9 Solutions (Tues 15th - Fri 18th Nov)
Algorihm and Daa Srucure 2011/ Week Soluion (Tue 15h - Fri 18h No) 1. Queion: e are gien 11/16 / 15/20 8/13 0/ 1/ / 11/1 / / To queion: (a) Find a pair of ube X, Y V uch ha f(x, Y) = f(v X, Y). (b) Find
More informationELEC 6041 LECTURE NOTES WEEK 3 Dr. Amir G. Aghdam Concordia University
ecre Noe Prepared b r G. ghda EE 64 ETUE NTE WEE r. r G. ghda ocorda Uer eceraled orol e - Whe corol heor appled o a e ha co of geographcall eparaed copoe or a e cog of a large ber of p-op ao ofe dered
More informationNONLOCAL BOUNDARY VALUE PROBLEM FOR SECOND ORDER ANTI-PERIODIC NONLINEAR IMPULSIVE q k INTEGRODIFFERENCE EQUATION
Euroean Journal of ahemac an Comuer Scence Vol No 7 ISSN 59-995 NONLOCAL BOUNDARY VALUE PROBLE FOR SECOND ORDER ANTI-PERIODIC NONLINEAR IPULSIVE - INTEGRODIFFERENCE EQUATION Hao Wang Yuhang Zhang ngyang
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 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 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 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 informationAlgorithmic Discrete Mathematics 6. Exercise Sheet
Algorihmic Dicree Mahemaic. Exercie Shee Deparmen of Mahemaic SS 0 PD Dr. Ulf Lorenz 7. and 8. Juni 0 Dipl.-Mah. David Meffer Verion of June, 0 Groupwork Exercie G (Heap-Sor) Ue Heap-Sor wih a min-heap
More informationCHAPTER 5: MULTIVARIATE METHODS
CHAPER 5: MULIVARIAE MEHODS Mulvarae Daa 3 Mulple measuremens (sensors) npus/feaures/arbues: -varae N nsances/observaons/eamples Each row s an eample Each column represens a feaure X a b correspons o he
More information6.302 Feedback Systems Recitation : Phase-locked Loops Prof. Joel L. Dawson
6.32 Feedback Syem Phae-locked loop are a foundaional building block for analog circui deign, paricularly for communicaion circui. They provide a good example yem for hi cla becaue hey are an excellen
More informationGraphene nanoplatelets induced heterogeneous bimodal structural magnesium matrix composites with enhanced mechanical properties
raphene nanoplaele nce heerogeneo bmoal rcral magnem marx compoe wh enhance mechancal propere Shln Xang a, b, Xaojn Wang a, *, anoj pa b, Kn W a, Xaoh H a, ngy Zheng a a School of aeral Scence an ngneerng,
More informationChapter 9 - The Laplace Transform
Chaper 9 - The Laplace Tranform Selece Soluion. Skech he pole-zero plo an region of convergence (if i exi) for hee ignal. ω [] () 8 (a) x e u = 8 ROC σ ( ) 3 (b) x e co π u ω [] ( ) () (c) x e u e u ROC
More informationu(t) Figure 1. Open loop control system
Open loop conrol v cloed loop feedbac conrol The nex wo figure preen he rucure of open loop and feedbac conrol yem Figure how an open loop conrol yem whoe funcion i o caue he oupu y o follow he reference
More informationNetwork Flow. Data Structures and Algorithms Andrei Bulatov
Nework Flow Daa Srucure and Algorihm Andrei Bulao Algorihm Nework Flow 24-2 Flow Nework Think of a graph a yem of pipe We ue hi yem o pump waer from he ource o ink Eery pipe/edge ha limied capaciy Flow
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 informationLesson 2 Transmission Lines Fundamentals
Lesson Transmsson Lnes Funamenals 楊尚達 Shang-Da Yang Insue of Phooncs Technologes Deparmen of Elecrcal Engneerng Naonal Tsng Hua Unersy Tawan Sec. -1 Inroucon 1. Why o scuss TX lnes srbue crcus?. Crera
More informationToday s focus. Bayes rule explained INF Multivariate classification Anne Solberg
Toay focu INF 4300 006 Mulvarae clafcaon Anne Solberg anne@fuono From a -menonal feaure vecor =[, ] T Gven K fferen clae k=, K Compue he probably ha belong o cla k Pk = p kpk/con How houl he mulvarae eny
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 informationA Demand System for Input Factors when there are Technological Changes in Production
A Demand Syem for Inpu Facor when here are Technologcal Change n Producon Movaon Due o (e.g.) echnologcal change here mgh no be a aonary relaonhp for he co hare of each npu facor. When emang demand yem
More informationMultivariable Dynamic Model and Robust Control of a Voltage-Source Converter for Power System Applications
Mularable Dynamc Moel an Robu Conrol of a Volage-Source Conerer for Power Syem Applcaon Ahmareza Tabeh an Reza Iraan Affne Conroller Parameerzaon for Decenralze Conrol Oer Banach Space Mchael Rokowz an
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 informationOptimal Filtering for Linear Discrete-Time Systems with Single Delayed Measurement
378 Hong-Guo Inernaonal Zhao, Journal Huan-Shu of Conrol, Zhang, Auomaon, Cheng-Hu an Zhang, Syem, an vol. Xn-Mn 6, no. Song 3, pp. 378-385, June 28 Opmal Flerng for Lnear Dcree-me Syem h Sngle Delaye
More informationProblem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow.
CSE 202: Deign and Analyi of Algorihm Winer 2013 Problem Se 3 Inrucor: Kamalika Chaudhuri Due on: Tue. Feb 26, 2013 Inrucion For your proof, you may ue any lower bound, algorihm or daa rucure from he ex
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 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 informationSSRG International Journal of Thermal Engineering (SSRG-IJTE) Volume 4 Issue 1 January to April 2018
SSRG Inernaonal Journal of Thermal Engneerng (SSRG-IJTE) Volume 4 Iue 1 January o Aprl 18 Opmal Conrol for a Drbued Parameer Syem wh Tme-Delay, Non-Lnear Ung he Numercal Mehod. Applcaon o One- Sded Hea
More informationIntroduction to Congestion Games
Algorihmic Game Theory, Summer 2017 Inroducion o Congeion Game Lecure 1 (5 page) Inrucor: Thoma Keelheim In hi lecure, we ge o know congeion game, which will be our running example for many concep in game
More informationPHYSICS 151 Notes for Online Lecture #4
PHYSICS 5 Noe for Online Lecure #4 Acceleraion The ga pedal in a car i alo called an acceleraor becaue preing i allow you o change your elociy. Acceleraion i how fa he elociy change. So if you ar fro re
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 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 informationLaplace Transformation of Linear Time-Varying Systems
Laplace Tranformaon of Lnear Tme-Varyng Syem Shervn Erfan Reearch Cenre for Inegraed Mcroelecronc Elecrcal and Compuer Engneerng Deparmen Unvery of Wndor Wndor, Onaro N9B 3P4, Canada Aug. 4, 9 Oulne of
More informationControl Systems. Mathematical Modeling of Control Systems.
Conrol Syem Mahemacal Modelng of Conrol Syem chbum@eoulech.ac.kr Oulne Mahemacal model and model ype. Tranfer funcon model Syem pole and zero Chbum Lee -Seoulech Conrol Syem Mahemacal Model Model are key
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 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 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 informationH = d d q 1 d d q N d d p 1 d d p N exp
8333: Sacal Mechanc I roblem Se # 7 Soluon Fall 3 Canoncal Enemble Non-harmonc Ga: The Hamlonan for a ga of N non neracng parcle n a d dmenonal box ha he form H A p a The paron funcon gven by ZN T d d
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 informationLaplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)
Laplace Tranform Maoud Malek The Laplace ranform i an inegral ranform named in honor of mahemaician and aronomer Pierre-Simon Laplace, who ued he ranform in hi work on probabiliy heory. I i a powerful
More informationNumerical Study of Large-area Anti-Resonant Reflecting Optical Waveguide (ARROW) Vertical-Cavity Semiconductor Optical Amplifiers (VCSOAs)
USOD 005 uecal Sudy of Lage-aea An-Reonan Reflecng Opcal Wavegude (ARROW Vecal-Cavy Seconduco Opcal Aplfe (VCSOA anhu Chen Su Fung Yu School of Eleccal and Eleconc Engneeng Conen Inoducon Vecal Cavy Seconduco
More informationDesign of Controller for Robot Position Control
eign of Conroller for Robo oiion Conrol Two imporan goal of conrol: 1. Reference inpu racking: The oupu mu follow he reference inpu rajecory a quickly a poible. Se-poin racking: Tracking when he reference
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 informationCooling of a hot metal forging. , dt dt
Tranen Conducon Uneady Analy - Lumped Thermal Capacy Model Performed when; Hea ranfer whn a yem produced a unform emperaure drbuon n he yem (mall emperaure graden). The emperaure change whn he yem condered
More informationEfficient Estimators for Population Variance using Auxiliary Information
Global Joural of Mahemacal cece: Theor ad Praccal. IN 97-3 Volume 3, Number (), pp. 39-37 Ieraoal Reearch Publcao Houe hp://www.rphoue.com Effce Emaor for Populao Varace ug Aular Iformao ubhah Kumar Yadav
More informationNovel Technique for PID Tuning by Particle Swarm Optimization
vel Technque for PID Tunng by Parcle Swarm Opmzaon *S. Eaer Selvan, Sehu Subramanan, S. Theban Solomon Abrac An aemp ha been mae by ncorporang ome pecal feaure n he convenonal parcle warm opmzaon (PSO)
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 informationHigher-order Graph Cuts
Example: Segmenaon Hgher-orer Graph Hroh Ihkawa 石川博 Deparmen of omper Scence & Engneerng Waea Unery 早稲田大学 Boyko&Jolly IV 3 Example: Segmenaon Local moel ex.: Moel of pxel ale for each kn of e Pror moel
More informationEngineering Circuit Analysis 8th Edition Chapter Nine Exercise Solutions
Engnrng rcu naly 8h Eon hapr Nn Exrc Soluon. = KΩ, = µf, an uch ha h crcu rpon oramp. a For Sourc-fr paralll crcu: For oramp or b H 9V, V / hoo = H.7.8 ra / 5..7..9 9V 9..9..9 5.75,.5 5.75.5..9 . = nh,
More information, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max
ecure 8 7. Sabiliy Analyi For an n dimenional vecor R n, he and he vecor norm are defined a: = T = i n i (7.) I i eay o how ha hee wo norm aify he following relaion: n (7.) If a vecor i ime-dependen, hen
More informationSample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems
Sample Final Exam Covering Chaper 9 (final04) Sample Final Exam (final03) Covering Chaper 9 of Fundamenal of Signal & Syem Problem (0 mar) Conider he caual opamp circui iniially a re depiced below. I LI
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 informationANALYSIS OF SOME SAFETY ASSESSMENT STANDARD ON GROUNDING SYSTEMS
ANAYSIS OF SOME SAFETY ASSESSMENT STANDARD ON GROUNDING SYSTEMS Shang iqun, Zhang Yan, Cheng Gang School of Elecrical and Conrol Engineering, Xi an Univeriy of Science & Technology, 710054, Xi an, China,
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 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 informationMATHEMATICAL MODEL OF THYRISTOR INVERTER INCLUDING A SERIES-PARALLEL RESONANT CIRCUIT
78 Avance n Elecrcal an Elecronc Engneerng MATHEMATIA MODE OF THYRISTOR INVERTER INUDING A SERIESPARAE RESONANT IRUIT M. uf, E. Szycha Faculy of Tranpor, Techncal Unvery of Raom, Polan ul. Malczewkego
More informationRevision: June 12, E Main Suite D Pullman, WA (509) Voice and Fax
.: apacors Reson: June, 5 E Man Sue D Pullman, WA 9963 59 334 636 Voce an Fax Oerew We begn our suy of energy sorage elemens wh a scusson of capacors. apacors, lke ressors, are passe wo-ermnal crcu elemens.
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 information2 shear strain / L for small angle
Sac quaons F F M al Sess omal sess foce coss-seconal aea eage Shea Sess shea sess shea foce coss-seconal aea llowable Sess Faco of Safe F. S San falue Shea San falue san change n lengh ognal lengh Hooke
More informationPhysics 3 (PHYF144) Chap 3: The Kinetic Theory of Gases - 1
Physcs (PYF44) ha : he nec heory of Gases -. Molecular Moel of an Ieal Gas he goal of he olecular oel of an eal gas s o unersan he acroscoc roeres (such as ressure an eeraure ) of gas n e of s croscoc
More informationPhysics 240: Worksheet 16 Name
Phyic 4: Workhee 16 Nae Non-unifor circular oion Each of hee proble involve non-unifor circular oion wih a conan α. (1) Obain each of he equaion of oion for non-unifor circular oion under a conan acceleraion,
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 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 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 information1 Motivation and Basic Definitions
CSCE : Deign and Analyi of Algorihm Noe on Max Flow Fall 20 (Baed on he preenaion in Chaper 26 of Inroducion o Algorihm, 3rd Ed. by Cormen, Leieron, Rive and Sein.) Moivaion and Baic Definiion Conider
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 informationCS344: Introduction to Artificial Intelligence
C344: Iroduco o Arfcal Iellgece Puhpa Bhaacharyya CE Dep. IIT Bombay Lecure 3 3 32 33: Forward ad bacward; Baum elch 9 h ad 2 March ad 2 d Aprl 203 Lecure 27 28 29 were o EM; dae 2 h March o 8 h March
More informationA L A BA M A L A W R E V IE W
A L A BA M A L A W R E V IE W Volume 52 Fall 2000 Number 1 B E F O R E D I S A B I L I T Y C I V I L R I G HT S : C I V I L W A R P E N S I O N S A N D TH E P O L I T I C S O F D I S A B I L I T Y I N
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 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 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 informationCHAPTER 7: SECOND-ORDER CIRCUITS
EEE5: CI RCUI T THEORY CHAPTER 7: SECOND-ORDER CIRCUITS 7. Inroducion Thi chaper conider circui wih wo orage elemen. Known a econd-order circui becaue heir repone are decribed by differenial equaion ha
More informationOn Line Supplement to Strategic Customers in a Transportation Station When is it Optimal to Wait? A. Manou, A. Economou, and F.
On Line Spplemen o Sraegic Comer in a Tranporaion Saion When i i Opimal o Wai? A. Mano, A. Economo, and F. Karaemen 11. Appendix In hi Appendix, we provide ome echnical analic proof for he main rel of
More informationSTUDY PROGRAM: UNIT 1 AND UNIT
IUIT ANAYSIS I MODUE ODE: EIAM4 STUDY POGAM: UNIT AND UNIT UT aal Unery of Technology EIAM4 haper : Fr Order rcu Page -. FIST ODE IUITS. Summary of Bac rcu oncep and onenon eor, capacor and nducor olage
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 informationP a g e 5 1 of R e p o r t P B 4 / 0 9
P a g e 5 1 of R e p o r t P B 4 / 0 9 J A R T a l s o c o n c l u d e d t h a t a l t h o u g h t h e i n t e n t o f N e l s o n s r e h a b i l i t a t i o n p l a n i s t o e n h a n c e c o n n e
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 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 informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 31 Signal & Syem Prof. Mark Fowler Noe Se #27 C-T Syem: Laplace Tranform Power Tool for yem analyi Reading Aignmen: Secion 6.1 6.3 of Kamen and Heck 1/18 Coure Flow Diagram The arrow here how concepual
More informationOP = OO' + Ut + Vn + Wb. Material We Will Cover Today. Computer Vision Lecture 3. Multi-view Geometry I. Amnon Shashua
Comuer Vson 27 Lecure 3 Mul-vew Geomer I Amnon Shashua Maeral We Wll Cover oa he srucure of 3D->2D rojecon mar omograh Marces A rmer on rojecve geomer of he lane Eolar Geomer an Funamenal Mar ebrew Unvers
More informationEE202 Circuit Theory II
EE202 Circui Theory II 2017-2018, Spring Dr. Yılmaz KALKAN I. Inroducion & eview of Fir Order Circui (Chaper 7 of Nilon - 3 Hr. Inroducion, C and L Circui, Naural and Sep epone of Serie and Parallel L/C
More informationThe Stochastic Jacobian in Affine Term Structure Models and its applications
he Sochac Jacoban n Affne erm Srucure Moel an applcaon Marno Graell Unverà egl Su Verona,anESILV Clauo ebal Unverà egl Su Verona Sepember 13, 2004 Abrac In Affne erm Srucure Moel (ASM) he ochac Jacoban
More informationMatrix reconstruction with the local max norm
Marx reconrucon wh he local max norm Rna oygel Deparmen of Sac Sanford Unvery rnafb@anfordedu Nahan Srebro Toyoa Technologcal Inue a Chcago na@cedu Rulan Salakhudnov Dep of Sac and Dep of Compuer Scence
More informationExperience with Optical Flow in Colour Video Image Sequences
Compuer Scence Deparmen of The Unvers of Aucklan CITR a Tamak Campus (hp://www.cr.aucklan.ac.nz) CITR-TR-102 Ocober 2001 perence wh Opcal Flow n Colour Veo Image Sequences John Barron 1 an Renhar Klee
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 informationLinear Algebra Primer
Linear Algebra rimer And a video dicuion of linear algebra from EE263 i here (lecure 3 and 4): hp://ee.anford.edu/coure/ee263 lide from Sanford CS3 Ouline Vecor and marice Baic Mari Operaion Deerminan,
More informationJohn Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany
Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy
More informationTP B.2 Rolling resistance, spin resistance, and "ball turn"
echnical proof TP B. olling reiance, pin reiance, and "ball urn" upporing: The Illuraed Principle of Pool and Billiard hp://billiard.coloae.edu by Daid G. Alciaore, PhD, PE ("Dr. Dae") echnical proof originally
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 information1.B Appendix to Chapter 1
Secon.B.B Append o Chper.B. The Ordnr Clcl Here re led ome mporn concep rom he ordnr clcl. The Dervve Conder ncon o one ndependen vrble. The dervve o dened b d d lm lm.b. where he ncremen n de o n ncremen
More informationDiscussion Session 2 Constant Acceleration/Relative Motion Week 03
PHYS 100 Dicuion Seion Conan Acceleraion/Relaive Moion Week 03 The Plan Today you will work wih your group explore he idea of reference frame (i.e. relaive moion) and moion wih conan acceleraion. You ll
More informationNON-HOMOGENEOUS SEMI-MARKOV REWARD PROCESS FOR THE MANAGEMENT OF HEALTH INSURANCE MODELS.
NON-HOOGENEOU EI-AKO EWA POCE FO THE ANAGEENT OF HEATH INUANCE OE. Jacque Janen CEIAF ld Paul Janon 84 e 9 6 Charlero EGIU Fax: 32735877 E-mal: ceaf@elgacom.ne and amondo anca Unverà a apenza parmeno d
More information