MAGNUM - A Fortran Library for the Calculation of Magnetic Configurations

Size: px
Start display at page:

Download "MAGNUM - A Fortran Library for the Calculation of Magnetic Configurations"

Transcription

1 CRYO/6/34 September, 3, 6 MAGNUM - A Fortran Lbrary for the Calculaton of Magnetc Confguratons L. Bottura Dstrbuton: Keywords: P. Bruzzone, M. Calv, J. Lster, C. Marnucc (EPFL/CRPP), A. Portone (EFDA- CSU) Magnetc feld Summary Ths note reports general formulae for the calculaton of magnetc feld generated by conductors n plane -D, axsymmetrc -D and 3-D confguratons, n the presence of current and magnetzaton sources. All formulae have been programmed n numercally stable routnes, collected n a lbrary named MAGNUM. Introducton In ths note we revew the state of the art of ntegral formulae for the calculaton of the magnetc feld n lnear and magnetsed meda. The formulae have been collected n a lbrary of FORTRAN routnes that compute the magnetc feld and the vector potental generated by dstrbutons of current and magnetzaton n several confguratons of nterest for magnet desgn and analyss. The routnes are organsed accordng to the feld source (current or magnetzaton), and spatal symmetry, for the followng confguratons: plane -D axsymmetrc -D general 3-D Most of the calculatons are based on analytcal formulae for the feld, and numercal ntegraton or functon approxmaton s lmted to the mnmum strctly necessary. In partcular, all plane -D and 3-D confguratons are computed analytcally. Axsymmetrc -D confguratons requre a numercal approxmaton to ellptc ntegrals and numercal ntegraton n one dmenson. 1

2 We descrbe below the confguratons, as well as the relevant analytcal formulae programmed n the lbrary. The am of ths note s to descrbe the workng prncple of the routnes, and the bascs of magnetostatcs are ntentonally left to the numerous textbooks and references quoted on the subjects. Plane, -D confguratons Long straght current lnes normal to a plane (x,y), or magnetc moments whch are unform n one space dmenson and have components only n the (x,y) plane generate magnetc feld confguratons that are -D and plane. The magnetc feld only has the two components n the plane (x,y), and the vector potental only has a component normal to the plane, along z. Under ths assumpton we can make an extensve use of the complex formalsm ntroduced by Beth (Beth, 1966). The poston of the current and magnetzaton sources s defned n the complex plane ζ = x + y. 1 Fgure 1. Plane -D reference frame and confguratons consdered for the calculaton of magnetc feld and vector potental of current and magnetzaton sources. As dscussed n (Beth, 1966) a complex, analytcal functon B = B y + B x can be defned from the x- and y-components of the magnetc feld vector. We ndcate below the z-component of the current densty n a current carryng conductor wth J, and the total current wth I. By analogy wth the defnton of the complex feld functon, we defne a complex magnetzaton functon M = M y + M x, where the two components M y and M x are the magnetzaton denstes, unform n z, n a magnetsed materal. Fnally, when computng the feld of small magnetc moments, we make use of the complex magnetc moment per unt length p = p y + p x, where the two components p y and p x are the magnetc moments per unt length, unform n z, n the magnetsed materal. 1 Note that we ndcate wth ζ the complex varable, to dstngush t from the axs coordnate z n the general 3-D frame (x,y,z).

3 The advantage of ths defnton s that, as B s analytcal, t can be expanded n a seres that converges n a crcle that does not contan current sources: n)1 # " B = + C n R ( (1) n=1 ref where C n s called complex harmonc, or multpole coeffcent, of order n, and R ref s a reference radus that has a pure normalzaton functon. The complex multpole coeffcents have a real and an magnary part: C n = B n + A n (). The real part B n s the so-called normal multpole, whle the magnary part, A n s the skew multpole. Magnetc feld calculaton Flamentary current The magnetc feld generated n a pont ζ P = x P + y P by a flamentary current I located at ζ C = x C + y C s gven by: B = µ I "# (3) where the dstance ρ s defned as: " = # C # P Unform current densty wth polygonal boundary (4) We consder the case of a polygonal conductor, delmted by a seres of N pecewse straght segments defned n the complex plane through the poston of the N vertces ζ. A unform current densty J flows nsde the polygonal conductor normal to the complex plane, see Fg.. 3

4 Fgure. Polygonal conductor wth unfom current densty n the drecton normal to the plane. The feld s computed n a pont z p that can be n an arbtrary poston n the plane. The feld generated at any pont ζ P n the complex plane s gven by the surface ntegral: B = µ J " S 1 ( ) ds # # P (5) The above ontegral can be transformed n the followng lne ntegral (Beth, 1966), (Beth, 1967), (Halbach, 197): B = 4" µ J # # P d# (6) # # P where the varable of ntegraton ζ descrbes the contour of the polygon n postve (ant-clockwse) drecton. Introducng the dstance vector ρ, as defned n Eq. (4), and recallng that the contour of the polygon between two vertces -1 and s a straght lne, t s possble to smplfy the above ntegral: r B = N 4" µ J # # d# (7). =1 r 1 where we mplctly assume that ρ = ρ N. Introducng the auxlary quanttes: "# = # # 1 "# = # # 1 (8) (9) we can solve the contour ntegral: B = N 4" µ J # + 1# 1 # ),( ( ln ln 1 ) + (1). # =1 4

5 The expresson above provdes an exact and closed form soluton for the feld generated by a polygonal conductor wth unform current densty. Its evaluaton s straghtforward for ponts ζ P outsde the conductor, but requres care n the monotonc treatment of the argument of the complex natural logarthm for ponts ζ P placed nsde and on the contour of the conductor. Localsed magnetc moment The magnetc feld generated n a pont ζ P = x P + y P by a localsed magnetc moment per unt length p wth components (p x,p y ) and located at z M = x M + y M s gven by: B = "µ p # (11) where the dstance ρ s defned smlarly to Eq. (4). " = # M # P Unform magnetzaton wth polygonal boundary (1) As for the magnetc feld, we consder the polygonal conductor of Fg.. The materal n the polygon has a unform magnetzaton M wth components (M x, M y ). The feld generated at any pont ζ P n the complex plane s gven by the surface ntegral: B = " µ M 1 # ( " P ) ds (13) S For ponts external to the polygon, the above ntegral can be transformed n the followng lne ntegral (Beth, 1966; Beth, 1967; Halbach, 197): B = " 1 4# µ 1 M d (14) " P wth the same conventons as Eq. (6). For ponts nternal to the polygon the result s the same, but the magnetc feld n ths case s gven by the result of Eq. (14) plus the contrbuton of the magnetzaton, µ M. We ntroduce the dstance vector r, as defned n Eq. (4), and we wrte the above ntegral decomposng over the sdes of the polygon: In general the logarthm of a complex nuber ζ can be wrtten as: log(ζ) = log( ζ ) + arg(ζ) whle the log of the module of the complex number s well behaved and always sngle-valued, the argument can flp by π dependng on the evaluaton method. The care that s needed s to make sure that under all crcumstances the argument of the terms n the lne ntegral s a monotoncal ncreasng or decreasng functon wth no jumps. 5

6 B = " 1 N 4# µ 1 M d (15). =1 "1 where we mplctly assume that ρ = ρ N. Introducng the auxlary quanttes: "# = # # 1 "# = # # 1 (16) (17) we can solve the contour ntegral: B = " 1 N 4# µ M ),( ( ln " ln "1 ) + (18). =1 The expresson above but requres the same care n the monotonc treatment of the argument of the complex natural logarthm for ponts ζ P placed nsde and on the contour of the conductor as Eq. (1). Vector potental calculaton Flamentary current The vector potental generated n a pont ζ P = x P + y P by a flamentary current I located at ζ C = x C + y C s gven by: + I A z = "Re µ # ln + 1 (., / - ) (19) where the dstance ρ s defned as n Eq. (4). Unform current densty wth polygonal boundary The vector potental generated by a polygonal conductor as n Fg.., carryng a unform current dstrbuton s gven by the followng surface ntegral: A z = " 1 4# µ J ln ( " P ) + ln [ ( " P )]ds (). S that can be transformed (Halbach, 197) n the lne ntegral: A z = 8" µ J # ( # P ) ln (# # P ) + ln # [ ( # P )]d# (1). As prevously, we ntroduce the dstance vector ρ as n Eq. (4), and we break the ntegral n a summaton along the straght sdes of the polygon: 6

7 r A z = N 8" µ J ( # ln # + ln# )d# (). =1 r 1 The ntegral can be solved and leads to the followng closed form for the vector potental of a polygonal conductor wth unform current densty: A z = N 8" µ J, # 1 # # 1 # ) /-( + # ln # 1. # =1 [ ( ) # 1 ( ln # 1 1) ] # # # ln# 1 ) ( + # 1 ln # 1 ) 3 ( # # # ln # 1 ) ( + # 1 ( 1 ln# ) (3) that can be evaluated wth the same precautons taken for the magnetc feld (see prevous sectons). Localsed magnetc moment The vector potental generated n a pont ζ P = x P + y P by a localsed magnetc moment p wth components (p x,p y ) and located at ζ m = x m + y m s gven by: p ( A z = "Re µ ) # (4) where the dstance ρ s defned as n Eq. (1). Unform magnetzaton wth polygonal boundary The vector potental generated by a polygonal dstrbuton of unform magnetzaton, wth the polygon defned as n Fg.., s gven by the followng surface ntegral: p A z = "Re µ ( # S 1 ( ) ds " P ) + (5) Apart for a constant, the term n parentheses, and n partcular the ntegral, s the same as that of the magnetc feld of a unform current densty n a polygon, Eq. (1). Hence we can use the same procedure, and the fnal result s: - N A z = "Re 4# µ p + "1 " "1 ).,( ( ln " ln "1 ) + 1 / =1 (6). 7

8 Harmoncs calculaton Flamentary current A lne of current I located at a poston ζ C = x C + y C n the complex plane generates the followng harmoncs: I R C n = "µ ref ( #R ref ) C n (7). Localsed magnetc moment A localsed magnetc moment p wth components (p x,p y ) and located at a poston ζ m = x m + y m n the complex plane generates the followng harmoncs: C n = "nµ p #R ref R ref ( ) n +1 (8). Axsymmetrc, -D confguratons In the case of axsymmetrc crcular flamentary or massve cols, the current lnes are normal to the plane (R,z) n clndrcal coordnates, and magnetc moments have components only n the (R,z). Both currents and magnetc moments are constant n the thrd space dmenson θ. The resultng magnetc feld has only the two components n the plane (R,z), and the vector potental has only one component normal to the (R,z) plane, along θ. We lmt ourselves to flamentary currents or magnetc moments, or to polygons n -D plane wth constant current densty and volume magnetzaton. Fgure 3. Cylndrcal reference frame for the axsymmetrc -D confguraton, and geometres consdered for the calculaton of magnetc feld and vector potental of current and magnetzaton sources. 8

9 The calculaton of felds n ths confguraton nvarably results n expresson nvolvng the ellptc ntegrals, and n partcular the complete ellptc ntegral of frst knd K and E, defned as follows: K( k) = E k d" (9) 1# k sn " / / ( ) 1/ ( ) 1/ d# (3) ( ) = 1" k sn # Magnetc feld calculaton Flamentary loop current The magnetc feld generated n a pont (R P,z P ) by a flamentary loop current I wth radus R C and placed at an heght z C s gven by: B R = "µ I B z = µ " I z 1 # R p R C + R P [( ) + z ] K k 1/ 1 [( R C + R P ) + #z ] K k 1/ ( ) " R C + R P + z ( ) + z E ( k ) R C " R P ( ) + R C R P #z ( ) + #z E ( k ) R C R P ( ) ( ) (31) (3) where: k = 4R C R P ( ) + "z R C + R P (33) and the dstance Δz s defned as: "z = z P # z C (34) Unform current densty loop wth rectangular boundary NOTE: Ths feld prmtve s presently not documented Unform current densty loop wth polygonal boundary NOTE: Ths feld prmtve s presently not documented Localsed magnetc moment loop NOTE: Ths feld prmtve s presently not documented 9

10 Unform magnetzaton loop wth rectangular boundary NOTE: Ths feld prmtve s presently not documented Unform magnetzaton loop wth polygonal boundary NOTE: Ths feld prmtve s presently not documented Vector potental calculaton Flamentary current loop The vector potental generated n a pont (R P,z P ) by a flamentary loop current I wth radus R C and placed at an heght z C s gven by: I A " = µ k# R c R p ) ( 1/ + 1 k - ) K k, ( ( ) E( k) where k s defned as n Eq. (33).. (35) / Unform current densty loop wth rectangular boundary NOTE: Ths feld prmtve s presently not documented Unform current densty loop wth polygonal boundary NOTE: Ths feld prmtve s presently not documented Localsed magnetc moment loop NOTE: Ths feld prmtve s presently not documented Unform magnetzaton loop wth rectangular boundary NOTE: Ths feld prmtve s presently not documented Unform magnetzaton loop wth polygonal boundary NOTE: Ths feld prmtve s presently not documented 3-D confguratons In the case of general, 3-D confguratons, the current and magnetzaton vectors can have arbtrary orentaton n the (x, y, z) space. The resultng magnetc feld and vector potental also have three non-zero components. 1

11 Fgure 4. Reference frame for the 3-D confguraton, and geometres consdered for the calculaton of magnetc feld and vector potental of current and magnetzaton sources. For 3-D confguratons we restrct the lbrary of modellng elements to the case of localsed sources (flament current I [A] or pont-lke magnetc moment P [Am ]), and to a general volume element wth plane faces and constant current densty J [A/m] or constant volume magnetzaton M [A/m]. Furthermore, although the volume element could have an arbtrary number of nodes (the equatons below reflect ths case), the mplementaton has been done only for a 8-node hexahedron, whch s the most useful modellng element for col wndng packs. To smplfy the notaton, t s useful to defne local reference frames. In the case of a current flament or a localsed magnetc moment we defne a local reference frame (x,y,z ) orented such that the z axs has the drecton of the current, or of the magnetc moment, and centered n the center of the current flament, or n the locaton of the magnetc moment. The orentaton of the other two axes, x and y, s nessental, as the results are nvarant for a rotaton around z. Ths reference frame s schematcally shown n Fg. 5. Fgure 5. Local reference frames on faces and sdes of a source volume. 11

12 Fgure 6. Local reference frames on faces and sdes of a source volume. In the case of a volume element we defne two frames, one on the plane faces and one on the sdes of the element. For the -th plane face delmtng the volume we defne the local reference frame (x,y,z ) wth z axs orented towards the external of the element. Ths reference frame s shown schematcally n Fg. 6. Note that wth ths choce the z component of the dstance vector between a feld calculaton pont and a source pont on the face s a constant. In turn, each face s delmted by straght sdes. For the j-th sde we defne a local reference frame (x,y,z ) wth z axs parallel to z defned above, and y parallel to the sde. Ths reference frame s also shown schematcally n Fg. 6. Note that wth ths choce the y component of the dstance vector between a feld calculaton pont and a source pont on the face s a constant. Furthermore, the equatons descrbng the sde becomes y = and z =. In general, the versors n the drecton of the axes of the local reference frames ((x,y,z ) or (x,y,z )) are ndcated as t, s, and n respectvely. The transformaton of a vector g from the cartesan frame (x, y, z) to the cartesan frame (t, s, n) s obtaned by the followng matrx relatons: " g x g ( x,y,z) = g y # g z " g t g ( t,s,n) = g s # g n " t x s x n x " = t y s y n y # t z s z n z # g t g s g n " t x t y t z " = s x s y s z # n x n y n z # = Tg ( t,s,n) (36) g x g y g z = T (1 g ( x,y,z) (37) 1

13 Magnetc feld calculaton Flamentary current The magnetc feld generated n a pont (x P,y P,z P ) by a straght flament of current I wth extremes (x C1, y C1, z C1 ) and (x C, y C, z C ) s gven n the local coordnate frame (x, y, z ) orented along the flament and shown n Fg. 5: B x = "µ I 4#, y p, x p +y, p, +, L " z p x p +y p + L " z p) ( + - L + z / p / x p +y p + L / + z p) / (./ (38) B y = µ I 4", x p, x p +y, p, +, L # z p x p +y p + L # z p) ( + - L + z / p / x p +y p + L / + z p ) / (./ (39) B z = (4) where L s the total length of the current flament, and x P, y P, z P are the coordnates of the feld pont n the local reference frame (x, y, z ). The components of the feld n the global reference frame (x, y, z) are obtaned by rotaton of the above values usng Eq. (36). Unform current densty n a sold volume The magnetc feld generated by a current densty dstrbuton n an arbtrary volume V s gven by: B = µ 4" J # r dv (41) V r 3 where r s the vector from the source pont (x Q, y Q, z Q ) to the feld pont (x P, y P, z P ): r = [ x P " x Q y P " y Q z P " z Q ] (4) We make the hypothess that the current densty s constant n the volume, and that the volume s delmted by plane faces. In ths case the calculaton of (41) can be performed analytcally usng the procedure devsed n (Colle, 1976) to frst reduce volume ntegrals to surface ntegrals, and then reduce surface ntegrals to lne ntegrals. We sketch here the process, wthout enterng n the detals of the method descrbed n the above reference. 13

14 As the current densty s constant, the ntegral n Eq. (4.6) can be wrtten: B = µ 4" J # r r dv (43). 3 V The vector n the ntegral can be wrtten as follows 3 : r r = "# 1 ) 3 r( (44). If we ndcate wth S the surface boundng the volume V, and wth n ts normal pontng towards the outsde of the element, we can use the followng property of vector functons: # "g dv = # gn ds (45) V S whch relates the volume ntegral of the dvergence of g to the flux of g over the surface S. We can transform the ntegral n Eq. (4.8) as follows: B = " µ 4# J ) ( n ) S 1 + r ds ), -) (46). Above we have ndcated wth S the -th plane surface composng the boundary of the volume element, and wth n ts normal (constant on the surface S ). As the ntegral s nvarant to a change of reference frame, we can choose that each surface ntegral s performed on the local reference frame (x, y, z ) shown n Fg. 6. In ths frame a pont on the surface S has coordnates (x Q, y Q, ), and the feld pont (x P, y P, z P ) s hence at a constant dstance along z. We can reduce further the dmenson of the ntegral by usng the followng vector relaton: " # g ds = g # s dl (47) S l where l s the curve enclosng the surface S, s s the normal to the curve, n the plane of S, pontng towards the outsde of the surface, and the dvergence s ntended as taken n the plane (.e. no dervatve n the drecton n normal to S). To apply Eq. (47) we use the dentty: 1 r = " # r r + z P ) ( (48). 3 the gradent s taken wth respect to the source pont (.e. the runnng varable n the ntegraton). 14

15 where we have defned the vector r wth components: r= [ x P "x Q y P "y Q ] (49). Usng Eqs. (47) and (48), we have that Eq. (45) yelds: B = " µ - 4# J /. n ) / ( j l j r r + z P 1 / s j dl, + 3 / (5). Above we have ndcated wth l j the j-th straght lne of the boundary of the surface S, and wth s ts normal (constant on the lne l ). We choose now to ntegrate n the local reference frame (x, y, z ) shown n Fg. 6. In ths frame a pont on the lne l has coordnates (x Q,, ), and the feld pont (x P, y P, z P ) s hence at a constant dstance both along y and z. Furthermore, the product r s j s the dstance of the feld pont along y, and we can hence wrte: B = µ, 4" J #. n 1 - ( y P /. j l j r + z P ). dl + 1. (51). The last step s to solve the lne ntegral above. The general soluton s: 1 " r + z dx = ln x + r ( ) + z x z y tg#1 where we have ndcated wth x ) # tg #1 )) = I 1 x, y,z yr ( y (( ( ) (5) r = x + y + z (53). We can use the above result to wrte the followng expresson for the magnetc feld: B = " µ, 4# J. - n ( y P I 1 x Q "x P, y P,z P /. j [ ( ) " I 1 ( x Q1 "x P, y P,z P )] ) (54) where we have ndcated wth x Q1 and x Q the coordnates of the begnnng and end of the lne l j, and we have made use of the fact that n the reference frames (x, y, z ) and (x, y, z ) the z-coordnate of the fld pont s the same,.e. we have z P = z P. 15

16 Localsed magnetc moment The magnetc feld generated n a pont (x P,y P,z P ) by a localsed magnetc moment P located at (x M, y M, z M ) s gven n the local coordnate frame (x, y, z ) orented along the magnetc moment and shown n Fg. 5: B " = 3µ P 4# B " = 3µ P 4# B " = µ P 4# x P z P x P +y ( P +z P ) (55) 5 / y P z P x P +y ( P +z P ) (56) 5 / z P x P y P x P +y ( P +z P ) (57) 5 / where P s the module of the magnetc moment and x P, y P, z P are the coordnates of the feld pont n the local reference frame (x, y, z ). The components of the feld n the global reference frame (x, y, z) are obtaned by rotaton of the above values usng Eq. (36). Unform magnetzaton n a sold volume The magnetc feld generated by a unform magnetzaton n a volume can be obtaned n dfferent ways. We choose here the followng expresson: µ M " µ ) B = ( " µ ) 4# M r r 3 dv V 4# M r r 3 dv V nsde the volume V outsde the volume V (58) where we used common conventons n vector calculus, and n partcular the defnton of the operator: M " # = M x x o +M y y o +M z z o (59). Examnng Eq. (58) we see that the volume ntegral s the same as already solved for the magnetc feld,.e. Eq. (43). It s hence possble to use the result obtaned there, and n partcular Eq. (54), to wrte: 16

17 " µ 4# M r r dv = 3 V = " µ. 4# M ( / n y P 1 ) j I 1 x Q "x P, y P,z P [ ( ) " I 1 ( x Q1 "x P, y P,z P )] + - 3, 4 (6) The gradent s ntended to be taken wth respect to the feld pont. The dervaton of the expresson under the summatons can be smplfed by recallng that: " "x = " "x "x "x + " "y "y "x + " "z "z "x = t x " "x + s x " "y + n x " "z (61). where we used the components of the versors t, s, n n the frame (x, y, z) as defned n Eqs. (36) and (37). It s hence possble to derve the expresson under the summaton n the local frame (x, y, z ) and transform the dervatves n the frame (x, y, z) usng Eq. (61) and analogous for the other drectons to solve Eq. (6). Vector potental calculaton Flamentary current The vector potental generated n a pont (x P,y P,z P ) by a straght flament of current I wth extremes (x C1, y C1, z C1 ) and (x C, y C, z C ) s gven n the local coordnate frame (x, y, z ) orented along the flament and shown n Fg. 5: A x = A y = ) # I A z = µ 4" ln z p + L + x # p +y p + z p + L + ( + ( #,ln z p, L + x # p +y p + z p, L - ( / ( /. (6) (63) (64) where L s the total length of the current flament, and x P, y P, z P are the coordnates of the feld pont n the local reference frame (x, y, z ). The components of the vector potental n the global reference frame (x, y, z) are obtaned by rotaton of the above values usng Eq. (36). 17

18 Unform current densty n a sold volume We use here the same technque as used for the magnetc feld calculaton. The the vector potental A generated by an arbtrary dstrbuton of current J n a volume V s: A = µ 4" # J r dv (65) V where r s the vector from the source pont (x Q, y Q, z Q ) to the feld pont (x P, y P, z P ) defned as n Eq. (4). Under the hypothess of constant current densty, we have: A = µ 4" J # 1 r dv (66). V The scalar n the ntegral can be wrtten as follows: 1 r = 1 " # r ) r ( (67). and n accordance wth the relaton Eq. (45) we can transform the volume ntegral of Eq. (67) n the followng surface ntegral: A = µ J ( 4" )( S r # n r ( ds+,( (68) Above we have ndcated wth S the -th plane surface composng the boundary of the volume element, and wth n ts normal (constant on the surface S ). We choose that each surface ntegral s performed on the local reference frame (x, y, z ) shown n Fg. 6. In ths frame the scalar product r n s constant, and s equal to z P, the z component of the dstance of the feld pont from the surface. Based on ths, the surface ntegral becomes: A = µ 4" J 1 ( 1 ) z P # r ds + S (69) We note now that the surface ntegral n Eq. (69) s dentcal to that already solved for the magnetc feld generated by a unform current densty n a volume element, n Eq. (46). We can use the procedure already detaled there to obtan the followng summaton of lne ntegrals n the local reference frame (x, y, z ): A = µ + 4" J -, z P y P # j l j r + z P (/ - dl ) 1 - (7). 18

19 Fnally, usng the defnton of the prmtve n Eq. (5), we have: A = µ + 4" J -, 1.- z P y P I 1 x Q #x P, y P,z P j [ ( ) # I 1 ( x Q1 #x P,y P,z P )] (/- ) 1 - (71). where we have ndcated wth x Q1 and x Q the coordnates of the begnnng and end of the lne l j, and we have made use of the fact that n the reference frames (x, y, z ) and (x, y, z ) the z-coordnate of the fld pont s the same,.e. we have z P = z P. Localsed magnetc moment The vector potental generated n a pont (x P,y P,z P ) by a localsed magnetc moment P located at (x M, y M, z M ) s gven n the local coordnate frame (x, y, z ) orented along the magnetc moment and shown n Fg. 5: A x = " A y = A z = µ P y P 4# x P +y P +z P ( ) µ P x P 4" x P +y P +z P ( ) (7) (73) (74) The vector potental generated n a pont (x P,y P,z P ) by a localsed magnetc moment P located at (x M, y M, z M ) s gven n the local coordnate frame (ξ, η, ζ) orented along the magnetc moment: where P s the module of the magnetc moment and x P, y P, z P are the coordnates of the feld pont n the local reference frame (x, y, z ). The components of the vector potental n the global reference frame (x, y, z) are obtaned by rotaton of the above values usng Eq. (36). Unform magnetzaton n a sold volume The the vector potental A generated by an arbtrary dstrbuton of magnetzaton M n a volume V s: A = µ 4" M # r dv (75) V r 3 where r s the vector from the source pont (x Q, y Q, z Q ) to the feld pont (x P, y P, z P ) defned as n Eq. (4). Under the hypothess of constant magnetzaton densty, we have: 19

20 A = µ 4" M # r r dv (76) 3 V The above ntegral s dentcal to the expresson for the magnetc feld generated by a constant current densty n a volume (see Eq. (43)), where the vector potental takes the place of the magnetc feld and the magnetzaton takes the place of the current densty. Under the hypothess of a volume delmted by plane faces, and usng the local coordnate frame (x, y, z ) as shown n Fg. 6, we have: A = " µ, 4# M. - n ( y P I 1 x Q "x P, y P,z P /. j [ ( ) " I 1 ( x Q1 "x P, y P,z P )] ) (77) where the prmtve I 1 s defned n Eq. (5), and we have ndcated wth x Q1 and x Q the coordnates of the begnnng and end of the lne l j. References (Beth, 1966) R.A. Beth, Complex Representaton and Computaton of Two- Dmensonal Magnetc Felds, J. Appl. Phys, 37 (7), , (Beth, 1967) R.A. Beth, An Integral Formula for Two-Dmensonal Felds, J. Appl. Phys, 38 (1), , (Colle, 1976) (Halbach, 197) C.J. Colle, Magnetc Felds and Potentals of Lnearly Varyng Current or Magnetsaton n a Plane Bounded Regon, Proceedngs of Compumag Conference, Oxford, 86-95, 1976 K. Halbach, Felds and Frst Order Perturbaton Effects n Two- Dmensonal Conductor Domnated Magnets, Nucl. Inst. and Meth., 78, , 197.

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Three-dimensional eddy current analysis by the boundary element method using vector potential

Three-dimensional eddy current analysis by the boundary element method using vector potential Physcs Electrcty & Magnetsm felds Okayama Unversty Year 1990 Three-dmensonal eddy current analyss by the boundary element method usng vector potental H. Tsubo M. Tanaka Okayama Unversty Okayama Unversty

More information

Supplemental document

Supplemental document Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor Terje Haukaas University of British Columbia, Vancouver  The Q4 Element Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

PES 1120 Spring 2014, Spendier Lecture 6/Page 1

PES 1120 Spring 2014, Spendier Lecture 6/Page 1 PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

More information

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as: 1 Problem set #1 1.1. A one-band model on a square lattce Fg. 1 Consder a square lattce wth only nearest-neghbor hoppngs (as shown n the fgure above): H t, j a a j (1.1) where,j stands for nearest neghbors

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

From Biot-Savart Law to Divergence of B (1)

From Biot-Savart Law to Divergence of B (1) From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

For all questions, answer choice E) NOTA" means none of the above answers is correct.

For all questions, answer choice E) NOTA means none of the above answers is correct. 0 MA Natonal Conventon For all questons, answer choce " means none of the above answers s correct.. In calculus, one learns of functon representatons that are nfnte seres called power 3 4 5 seres. For

More information

Electricity and Magnetism - Physics 121 Lecture 10 - Sources of Magnetic Fields (Currents) Y&F Chapter 28, Sec. 1-7

Electricity and Magnetism - Physics 121 Lecture 10 - Sources of Magnetic Fields (Currents) Y&F Chapter 28, Sec. 1-7 Electrcty and Magnetsm - Physcs 11 Lecture 10 - Sources of Magnetc Felds (Currents) Y&F Chapter 8, Sec. 1-7 Magnetc felds are due to currents The Bot-Savart Law Calculatng feld at the centers of current

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

CIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56

CIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56 CIVL 8/77 Chapter 0 - Isoparametrc Formulaton 4/56 Newton-Cotes Example Usng the Newton-Cotes method wth = ntervals (n = 3 samplng ponts), evaluate the ntegrals: x x cos dx 3 x x dx 3 x x dx 4.3333333

More information

Physics 114 Exam 3 Spring Name:

Physics 114 Exam 3 Spring Name: Physcs 114 Exam 3 Sprng 015 Name: For gradng purposes (do not wrte here): Queston 1. 1... 3. 3. Problem 4. Answer each of the followng questons. Ponts for each queston are ndcated n red. Unless otherwse

More information

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom ------ Castglano s frst theorem ---- Examples 8.4 Prncpal of statonary

More information

Army Ants Tunneling for Classical Simulations

Army Ants Tunneling for Classical Simulations Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Numerical Transient Heat Conduction Experiment

Numerical Transient Heat Conduction Experiment Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use

More information

Math1110 (Spring 2009) Prelim 3 - Solutions

Math1110 (Spring 2009) Prelim 3 - Solutions Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Celestial Mechanics. Basic Orbits. Why circles? Tycho Brahe. PHY celestial-mechanics - J. Hedberg

Celestial Mechanics. Basic Orbits. Why circles? Tycho Brahe. PHY celestial-mechanics - J. Hedberg PHY 454 - celestal-mechancs - J. Hedberg - 207 Celestal Mechancs. Basc Orbts. Why crcles? 2. Tycho Brahe 3. Kepler 4. 3 laws of orbtng bodes 2. Newtonan Mechancs 3. Newton's Laws. Law of Gravtaton 2. The

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

A linear imaging system with white additive Gaussian noise on the observed data is modeled as follows:

A linear imaging system with white additive Gaussian noise on the observed data is modeled as follows: Supplementary Note Mathematcal bacground A lnear magng system wth whte addtve Gaussan nose on the observed data s modeled as follows: X = R ϕ V + G, () where X R are the expermental, two-dmensonal proecton

More information

Physics 4B. Question and 3 tie (clockwise), then 2 and 5 tie (zero), then 4 and 6 tie (counterclockwise) B i. ( T / s) = 1.74 V.

Physics 4B. Question and 3 tie (clockwise), then 2 and 5 tie (zero), then 4 and 6 tie (counterclockwise) B i. ( T / s) = 1.74 V. Physcs 4 Solutons to Chapter 3 HW Chapter 3: Questons:, 4, 1 Problems:, 15, 19, 7, 33, 41, 45, 54, 65 Queston 3-1 and 3 te (clockwse), then and 5 te (zero), then 4 and 6 te (counterclockwse) Queston 3-4

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2 P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons

More information

Lagrangian Field Theory

Lagrangian Field Theory Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Chapter 4 The Wave Equation

Chapter 4 The Wave Equation Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as

More information

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review Chapter Four States of Stress Part Three When makng your choce n lfe, do not neglect to lve. Samuel Johnson Revew When we use matrx notaton to show the stresses on an element The rows represent the axs

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

PHYS 1101 Practice problem set 12, Chapter 32: 21, 22, 24, 57, 61, 83 Chapter 33: 7, 12, 32, 38, 44, 49, 76

PHYS 1101 Practice problem set 12, Chapter 32: 21, 22, 24, 57, 61, 83 Chapter 33: 7, 12, 32, 38, 44, 49, 76 PHYS 1101 Practce problem set 1, Chapter 3: 1,, 4, 57, 61, 83 Chapter 33: 7, 1, 3, 38, 44, 49, 76 3.1. Vsualze: Please reer to Fgure Ex3.1. Solve: Because B s n the same drecton as the ntegraton path s

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

One Dimension Again. Chapter Fourteen

One Dimension Again. Chapter Fourteen hapter Fourteen One Dmenson Agan 4 Scalar Lne Integrals Now we agan consder the dea of the ntegral n one dmenson When we were ntroduced to the ntegral back n elementary school, we consdered only functons

More information

Pressure Measurements Laboratory

Pressure Measurements Laboratory Lab # Pressure Measurements Laboratory Objectves:. To get hands-on experences on how to make pressure (surface pressure, statc pressure and total pressure nsde flow) measurements usng conventonal pressuremeasurng

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m) 7.. = = 3 = 4 = 5. The electrc feld s constant everywhere between the plates. Ths s ndcated by the electrc feld vectors, whch are all the same length and n the same drecton. 7.5. Model: The dstances to

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order: 68A Solutons to Exercses March 05 (a) Usng a Taylor expanson, and notng that n 0 for all n >, ( + ) ( + ( ) + ) We can t nvert / because there s no Taylor expanson around 0 Lets try to calculate the nverse

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Physics for Scientists & Engineers 2

Physics for Scientists & Engineers 2 Equpotental Surfaces and Lnes Physcs for Scentsts & Engneers 2 Sprng Semester 2005 Lecture 9 January 25, 2005 Physcs for Scentsts&Engneers 2 1 When an electrc feld s present, the electrc potental has a

More information

Solution Thermodynamics

Solution Thermodynamics Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

PHY2049 Exam 2 solutions Fall 2016 Solution:

PHY2049 Exam 2 solutions Fall 2016 Solution: PHY2049 Exam 2 solutons Fall 2016 General strategy: Fnd two resstors, one par at a tme, that are connected ether n SERIES or n PARALLEL; replace these two resstors wth one of an equvalent resstance. Now

More information

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials Journal of Physcs: Conference Seres PAPER OPEN ACCESS A boundary element method wth analytcal ntegraton for deformaton of nhomogeneous elastc materals To cte ths artcle: Moh. Ivan Azs et al 2018 J. Phys.:

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

9 Characteristic classes

9 Characteristic classes THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14 APPROXIMAE PRICES OF BASKE AND ASIAN OPIONS DUPON OLIVIER Prema 14 Contents Introducton 1 1. Framewor 1 1.1. Baset optons 1.. Asan optons. Computng the prce 3. Lower bound 3.1. Closed formula for the prce

More information

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures IMS 2 Workshop Analytcal Gradent Evaluaton of Cost Functons n General Feld Solvers: A Novel Approach for Optmzaton of Mcrowave Structures P. Harscher, S. Amar* and R. Vahldeck and J. Bornemann* Swss Federal

More information

. The kinetic energy of this system is T = T i. m i. Now let s consider how the kinetic energy of the system changes in time. Assuming each.

. The kinetic energy of this system is T = T i. m i. Now let s consider how the kinetic energy of the system changes in time. Assuming each. Chapter 2 Systems of Partcles 2. Work-Energy Theorem Consder a system of many partcles, wth postons r and veloctes ṙ. The knetc energy of ths system s T = T = 2 mṙ2. 2. Now let s consder how the knetc

More information

kq r 2 2kQ 2kQ (A) (B) (C) (D)

kq r 2 2kQ 2kQ (A) (B) (C) (D) PHYS 1202W MULTIPL CHOIC QUSTIONS QUIZ #1 Answer the followng multple choce questons on the bubble sheet. Choose the best answer, 5 pts each. MC1 An uncharged metal sphere wll (A) be repelled by a charged

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0 Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the

More information

Magnetostatics. T i. Ursula van Rienen, Universität Rostock, FB Elektrotechnik und Informationstechnik, AG Computational Electrodynamics

Magnetostatics. T i. Ursula van Rienen, Universität Rostock, FB Elektrotechnik und Informationstechnik, AG Computational Electrodynamics Magnetostatcs grd G E dual grd G J Ch = j SM h = T wth h= h + S Φ curl = J Ch = j T dv grad ϕ = dv = q SM S Φ = SM h = q m m The basc equatons n magnetostatcs are smlar to the equatons n electrostatcs.

More information