Energy balance invariance for interacting particle systems

Size: px
Start display at page:

Download "Energy balance invariance for interacting particle systems"

Transcription

1 Z. ngew. Mh. Phys. 60 (009) /09/ DOI 0.007/s c 008 Birkhäuser Verlg, Bsel Zeischrif für ngewne Mhemik un Physik ZAMP Energy blnce invrince for inercing pricle sysems Arsh Yvri n Jerrol E. Mrsen Absrc. This pper suies he principle of invrince of blnce of energy n is consequences for sysem of inercing pricles uner groups of rnsformions. Blnce of energy n is invrince is firs exmine in Euclien spce. Unlike he cse of coninuous mei, i is shown h conservion n blnce lws o no follow from he ssumpion of invrince of blnce of energy uner ime-epenen isomeries of he mbien spce. However, he posule of invrince of blnce of energy uner rbirry iffeomorphisms of he mbien (Euclien) spce, oes yiel he correc conservion n blnce lws. These ies re hen exene o he cse when he mbien spce is Riemnnin mnifol. Pirwise inercions in he cse of geoesiclly complee Riemnnin mbien mnifols re efine by ssuming h he inercion poenil explicily epens on he pirwise isnces of pricles. Posuling blnce of energy n is invrince uner rbirry ime-epenen spil iffeomorphisms yiels blnce of liner momenum. I is seen h pirwise forces re irece long ngens o geoesics heir en poins. One lso obins iscree version of he Doyle Ericksen formul, which reles he mgniue of inernl forces o he re of chnge of he ineromic energy wih respec o iscree meric h is rele o he bckgroun meric. Keywors. Coninuum mechnics, pricle mechnics, energy blnce, covrince.. Inroucion This pper is concerne wih blnce of energy for sysem of inercing pricles n fins connecion beween iscree blnce lws n invrince, somehing h hs been suie horoughly in he seing of coninuum mechnics (see [4, 9, 6] n references herein). Pricle mechnics is normlly formule in R 3 vi blnce lws such s blnce of momenum; however, he liner srucure of R 3 cn someimes obscure imporn geomeric informion; for exmple, blnce of liner momenum is no covrin noion in h i looks rher ifferen in curviliner coorines; see [9, 6]. Of course one wy o overcome his is o mke use of generlize Deice o he memory of Shhrm Kvinpour ( ). Reserch suppore by he Georgi Insiue of Technology. Reserch prilly suppore by he Cliforni Insiue of Technology n NSF-ITR Grn ACI

2 74 A. Yvri n J. E. Mrsen ZAMP coorines n Lgrngin mechnics. However, noher ineresing lernive pproch, h is lso use in geomeric coninuum mechnics, is o mke use of invrince properies of blnce of energy. A mjor gol of his pper is o revisi his issue o bring he heory more ino line wih wh one oes in geomeric coninuum mechnics. In geomeric coninuum mechnics, one usully works wih wo configurions reference configurion n curren configurion. The configurion spce is he mnifol of mps beween hese wo configurions [9]. The curren configurion hs cler physicl inerpreion; i is wh one cn see in he lborory. Equilibri of he curren configurion in he seing of coninuum mechnics correspon o locl minim of he corresponing pricle sysem in he cse of qusisic eformions. In he corresponing pricle sysem he resuln force on ech pricle is blnce by is inernl forces. In geomeric coninuum mechnics, he curren configurion evolves in Riemnninn mnifol (S, g) [9, 6]. The reference configurion on he oher hn hs less cler physicl inerpreion. In he riionl remens of elsiciy i is usully ssume h here is well-efine sress-free reference configurion. This is no lwys rue s ws noice by Eckr few eces go [, ]. For n elsic boy, in generl, one cn hve se of nurl configurions. So, sress-free reference configurion coul be ny of hese nurl configurions. In he corresponing omic sysem nurl configurion woul be locl minimum of he energy (or free energy) in he bsence of exernl forces. I shoul be noe h reference configurion is in some sense rbirry n one cn choose ifferen reference configurions for he sme problem (see [7] for some iscussions on his). We shoul emphsize h we re no criicizing he exising remens of pricle mechnics in Euclien spce. As mer of fc, his woul be he nurl wy of formuling, for exmple, moleculr sysems h re embee in Euclien spce. However, one shoul noe h even for he simple exmple for he clssicl rigi boy in roion bou is cener of mss, he configurion spce is mnifol, i.e. SO(3) n he kineic energy meric is non-euclien. The oher hing o noe is h suying mechnicl sysems geomericlly cn give nonrivil insigh ino mechnics of pricles in Euclien spce. A goo exmple of his is he reucion proceure for suying he ynmics of he shpe spce of molecule; h is, he spce obine when one elimines rnslions n roions. When his is one, one obins ynmics on non-euclien spce n geomeric mehos re criicl when unerking such suy. There re mny exmples of his in he lierure in he works of, for exmple, Lilejohn n Iwi. For concree exmple pplie o he nlysis of conformion chnges in Argon-6, see [5] n references herein. This pper is srucure s follows. In we briefly review he clssicl heorem of Green, Nghi n Rivlin n he covrince ies in elsiciy. In 3 we sr wih blnce of energy for sysem of inercing pricles embee in Euclien spce n suy he consequences of posuling is invrince uner ifferen

3 Vol. 60 (009) Energy blnce invrince for inercing pricle sysems 75 groups of rnsformions of he mbien spce. 4 suies hese issues when he mbien spce is Riemnnin mnifol. Consequences of covrince of energy blnce re invesige in eil. Conclusions re given in 5.. Energy blnce, he Green Nghi Rivlin heorem n covrince in elsiciy In every mechnicl sysem here re blnce lws n some ssocie conserve quniies. For exmple, in elsiciy, one hs blnce of liner n ngulr momen n conservion of mss. One cn buil coninuum heory by posuling he relevn blnce lws. However, one cn lwys quesion he significnce of ech blnce lw n wheher hey hve ny inrinsic menings. Mos engineering heories, incluing he heory of elsiciy, re riionlly buil wih he implici ssumpion h he mbien spces re Euclien. I urns ou h he srucure of Euclien spce cn obscure he covrince of blnce lws since hese blnce lws look rher ifferen in curviliner coorines. A more nurl wy of builing fiel heories is o ssume h he mbien spces re mnifols. This wy, for exmple, one nurlly obins covrin heory n, in iion, one cn more esily sepre he meric epenen n inepenen relions. I hs long been known h here is eep connecion beween blnce lws n symmeries (see, for insnce, [0]). Noeher s heorem, for exmple, ienifies conserve quniy for ech locl symmery of he unerlying Lgrngin ensiy. In he cse of globl symmeries in coninuum mechnics Green n Rivlin [4] showe h posuling blnce of energy n is invrince uner spil isomeries of he Euclien mbien spce, one cn obin conservion of mss n blnces of liner n ngulr momen. Moive by his observion n he fc h blnce of energy cn lwys be inrinsiclly efine Mrsen n Hughes [9] evelope covrin heory of elsiciy by posuling blnce of energy n is spil covrince. This ssumpion resuls in conservion of mss, blnce of liner n ngulr momen n he Doyle Ericksen formul. Simo n Mrsen [4] erive meril version of Doyle Ericksen formul in erms of he roe sress ensor. Recenly Yvri e l. [6] suie he covrince conceps in elsiciy in some eil, possibiliy of meril covrince of energy blnce n he connecion beween covrince n Noeher s heorem. In he cse of linerize elsiciy, recenly i ws shown h covrince of linerize energy blnce cn give ll he fiel equions of linerize elsiciy [9]. See lso [8] for connecion beween energy blnce of iscreize soli n is blnce lws. To he bes of our knowlege, here is no sysemic suy of possible connecion beween invrince n blnce lws for pricle sysems. In his pper, we firs consier he cse in which he mbien spce is Euclien n show h here re some suble ifferences beween coninuum elsiciy n sysem of

4 76 A. Yvri n J. E. Mrsen ZAMP inercing pricles. In priculr, we will show h he Green Nghi Rivlin (GNR) Theorem fils for pricle sysems. We hen efine pirwise inercions when he mbien spce is Riemnnin mnifol n suy he consequences of spil covrince of energy blnce. We will lso obin iscree version of he Doyle Ericksen formul. 3. Blnce of energy for pricle sysems in Euclien spce We firs consier sysem of inercing pricles in Euclien mbien spce n show h, unlike clssicl elsiciy, blnce lws o no follow from invrince of energy blnce uner ime-epenen isomeries of he Euclien mbien spce. In oher wors, he GNR Theorem fils. We hen show h moifie GNR rgumen using n enlrge group of rnsformions works n oes give he blnce lws for he pricle sysem. The filure of he GNR heorem Here, we look blnce of energy for n rbirry collecion of oms L n mke connecion beween invrince of blnce of energy for sysem of pricles n he clssicl Green Nghi Rivlin Theorem for coninuous mei. For he ske of simpliciy, we resric ourselves o pirwise inercions. Suppose we re given collecion L of pricles which hs he configurion {x i ()} R n ime. Blnce of energy for L cn be wrien s φ ij ( x i x j ) + m iẋ i ẋ i = F i ẋ i, (3.) where is he snr inner prouc of R n, m i is he mss of pricle i n F i is he exernl force on pricle i. Blnce of energy cn be simplifie o re f ij (v i v j ) + m i v i i = F i v i, (3.) where v i = ẋ i, i = ẍ i, f ij = φ ( ij x i x j ) x i x j x i x j x i x j, (3.3) n we hve use he snr frmework of clssicl mechnics in which he pricle msses re ssume o be ime inepenen. Noe h for pirwise inercions f ji = f ij. (3.4) Le us now posule h blnce of energy for L is invrin uner imeepenen rigi rnslion of he mbien spce S = R 3. Consier ξ : S S,

5 Vol. 60 (009) Energy blnce invrince for inercing pricle sysems 77 where x i = ξ (x i ) = x i + ( 0 )w, (3.5) for some rbirry w R 3. This implies h v i = v i + w for ll i L. (3.6) Also ime = 0 i is ssume h he boy forces F i rnsform s follows (see [9]): m i i F i = m i i F i for ll i L. (3.7) Invrince of blnce of energy mens h φ ( x i x j ) + m iv i v i = F i x i, (3.8) where he prime quniies re rele o he unprime ones hrough Crn s clssicl spceime heory. Simplifying (3.8) n evluing i ime = 0, we obin f ij (v i v j ) + m i (v i + w) i = F i (v i + w). (3.9) Subrcing (3.) from (3.9) yiels m i w i = F i w. (3.0) Becuse w is rbirry one cn conclue h m i i. (3.) F i = Eq. (3.) is nohing bu Newon s secon lw for he collecion of pricles. I is seen h he bove posule oes no give he known governing equions for ech pricle. Inse, i gives blnce of ol liner momenum for he whole collecion of pricles. Le us now look (3.) n ry o rewrie i for n rbirry subse M L. This is of course no lwys possible becuse in he collecion L, ech pricle inercs wih ll he oher pricles n blnce of energy for subcollecion cnno be wrien unmbiguously, in generl. Noe lso h energy my no even be conserve loclly []. In oher wors, for nonlocl sysems loclize blnce lw, in generl, involves so-clle resiul erm. In he cse of pirwise inercions, one my be empe o wrie he blnce

6 78 A. Yvri n J. E. Mrsen ZAMP of energy for subcollecion M L s follows. { φ ij ( x i x j ) + i,j M j L\M + φ ij( x i x j ) } m iẋ i ẋ i = F i ẋ i, (3.) where we hve ssume h he energy of he bon beween he pricles i n j is eqully shre by hem in he cse of pirwise inercions. Noe lso h unlike clssicl coninuum mechnics energy ensiy is no one-poin funcion, i.e., for given pricle energy hs conribuions from ll he oher pricles in he collecion, in generl. Blnce of energy (3.) cn be simplifie o re i,j M f ij (v i v j ) + j L\M f ij (v i v j ) + m i v i i = F i ẋ i. (3.3) Now le us posule h his blnce lw is invrin uner n rbirry imeepenen rigi rnslion of he mbien spce, i.e. i,j M f ij (v i v j )+ j L\M f ij (v i v j )+ m iv i i = F i x i. (3.4) Subrcing (3.3) from (3.4) evlue = 0 yiels f ij (w i w j )+ f ij (w i w j )+ m i w i i = F i w i, i,j M j L\M (3.5) where w i = ξ (x i ). Now if ξ is rigi rnslion, hen w i = c, for ll i L n hence m i c i = F i c. (3.6) Becuse c is rbirry, his hen implies h F i = m i i. (3.7) M is rbirry so we cn choose M = {i} n hence F i = m i i. (3.8) Noe h φij ( x i x j ) is he energy of he pir (i, j). There is fcor in he firs erm of he lef-hn sie of (3.) becuse in he sum boh (i, j) n (j, i) re coune. On he oher hn, he fcor ppers in he secon erm becuse i is ssume h he energy φ ij( x i x j ) is shre eqully beween he oms i n j.

7 Vol. 60 (009) Energy blnce invrince for inercing pricle sysems 79 This mens h he ineromic forces re self-equilibre, i.e. f ij = 0, (3.9) j L which is no lwys rue! This shows h he version of energy blnce given by (3.) is no invrin uner ime-epenen rigi rnslions, in generl. This goes bck o he nonlocl nure of inercions in pricle sysem n lso he presence of self inercions in given subse. In oher wors, blnce of energy for he whole sysem shoul be consiere. See [5] n [6] for similr iscussions. In summry, he GNR Theorem fils for pricle sysem. Remrk. In nonlocl sysem, one cn unmbiguously wrie globl blnce of energy, e.g. (3.) for our pricle sysem, bu pssge from globl energy blnce o locl energy blnce or o blnce of energy for n rbirry meril subse my no be unique [, 3]. In wh follows from here on we only wrie globl blnce of energy n suy is invrince (covrince) properies. The success of moifie GNR heorem Le us now consier lrger group of rnsformions for energy blnce invrince. Inse of consiering blnce of energy for M L n is invrince uner rigi rnslions, le us look blnce of energy for L n consier n rbirry C iffeomorphism ξ : S S. Thus, inse of (3.0) we hve where Noe h Hence f ij (w i w j ) + m i w i i = F i w i, (3.0) w i = ξ (x i ). (3.) f ij (w i w j ) = f ij w i. (3.) j L f ij w i = (F i m i i ) w i. (3.3) j L This cn be rewrien s [ F i j L f ij + m i i ] w i = 0. (3.4)

8 730 A. Yvri n J. E. Mrsen ZAMP Now ssuming h ξ is such h he w j = 0, j i, one obins F i + j L f ij = m i i, (3.5) which is wh one expecs. Le us now consier ime-epenen rigi roion in he mbien spce, i.e. ξ (x i ) = e ( 0)Ω x i (3.6) for some skew-symmeric mrix Ω. This mens h w i = Ωx i for ll i L. Obviously, posuling invrince uner rigi roions of he mbien spce oes no give ny new blnce lws s rigi roion is jus specil cse of iffeomorphisms consiere bove in (3.0). In oher wors, blnce of ngulr momenum is rivilly sisfie for collecion of inercing pricles. In summry, we hve prove he following proposiion. Proposiion 3.. For collecion of pricles wih pirwise inercions, posuling blnce of energy n is invrince uner isomeries of he mbien Euclien spce is no enough o fin ll he blnce lws (equions of moion). If, however, one posules blnce of energy n is invrince uner rbirry iffeomorphisms of Euclien spce, hen hese equions re ll obine. Remrk. In his pper, we o no consier inepenen roions, n in generl microsrucure, for pricles. However, he rgumens presene here cn be exene o more complice pricle sysems, e.g. pricle sysems wih nisoropic inercions, ec. 4. Energy blnce for pricle sysems on Riemnnin mnifols To pu he previous resuls in more generl frmework n lso for he ske of clriy, compleeness, n inellecul sisfcion, le us ssume h he pricles move on Riemnnin mnifol. Alhough his my seem oo bsrc firs glnce, i is nurl o sk how he blnce lws look when he mbien spce is Riemnnin. This is lso moive by he previous covrin formulions of coninuum elsiciy. We shoul lso menion h here hve been recen effors in formuling mechnics on non-euclien spces, e.g. [, 3]. Suppose he inercing pricles lie in Riemnnnin mnifol (S, g), which we ssume is geoesiclly complee. For he ske of simpliciy, le us ssume h only wo-boy inercions re presen. Consier wo pricles i, j L, which in he curren configurion lie in S, i.e., x i, x j S. We ssume h he poenil energy of hese wo pricles explicily epens on heir relive isnce in he Riemnin mnifol. Of course, here my be oher possibiliies bu his ssumpion woul be goo sring poin. As physicl sysem represene by his moel, le us consier finie number of pricles lying on wo-mnifol M R 3 connece o ech oher by some

9 Vol. 60 (009) Energy blnce invrince for inercing pricle sysems 73 springs h re consrine o lie in M. Energy of his sysem is obviously n explici funcion of lenghs of he springs. Given ny wo pricles, energy of he corresponing spring is minimize when he spring coincies wih he geoesic connecing he wo poins on M. Wh we fin covrinly in he following is blnce of forces projece on ngen spces of M pricle posiions. As noher moivion for suying his moel, we shoul menion he recen work by Koni n Sun [7]. Koni n Sun [7] consier weighe finie grph X = (V, E), where V n E re he verex n ege ses, respecively, n consier piecewise smooh mp Φ from X o Riemnnin mnifol (Y, g). They enoe he resricion of his mp o e E by Φ e (), [0, ]. Then, hey efine energy of his mp s E(Φ) = e E m E (e) 0 Φe, Φ e, (4.) where m E is weigh funcion efine on E. They show h mp Φ is criicl mp for his energy if n only if Φ e is geoesic for every e E. For such mp, energy is n explici funcion of geoesic lenghs. Remrk. Noe h, for generic Riemnnin mnifol n wo rbirry poins on he mnifol, here my be more hn one n even infiniely mny isnce minimizing geoesics connecing he wo poins. An exmple woul be he norh n souh pole of sphere. Here, we ssume h in given configurion of he pricle sysem here re no pirs of pricles lying on such pir of poins. Le us enoe he geoesic connecing x i n x j by l ij g, where i is cler from his noion h he geoesic explicily epens on he meric g. This curve hs prmerizion l : [, b] S n is lengh is efine s L g (l ij g ) = s lij g (s), s lij g (s) g s, (4.) where.,. g is he inner prouc inuce from he Riemnnin meric g. Therefore, for x i, x j S ( φ ij = φ ij L(l ij g ) ). (4.3) The ol inercion energy is efine s e ( {x i }, g ) = j L φ ij ( Lg (l ij g ) ). (4.4) Noe h we cn hink of g ij := L g (l ij g ) s iscree meric for he collecion {x i } S. Of course, his iscree meric is n explici funcion of he bckgroun Riemnnin meric g. Now blnce of energy cn be wrien s φ ij ( Lg (l ij g ) ) + m i ẋi, ẋ i g = F i, ẋ i g. (4.5)

10 73 A. Yvri n J. E. Mrsen ZAMP Noe h F i T x is is he exernl force on om i. Figure. A spil chnge of frme for sysem of inercing pricles. Le us ssume h uner spil chnge of frme ξ : S S blnce of energy is invrin, i.e. φ ij ( L g (l g ij ) ) + ẋ i m i, ẋ i g = F i, ẋ i. (4.6) g Noe h g = ξ g (see Fig. ). We ssume h he pirwise poenil rnsforms onsorilly, i.e. ( x i, x j, g ) ( = φ ij x i, x j, ξ g ). (4.7) φ ij Equivlenly φ ij ( Lg (l ij g ) ) = φ ij ( L ξ g(l ij ξ g) ). (4.8)

11 Vol. 60 (009) Energy blnce invrince for inercing pricle sysems 733 Blnce of energy for he new frming = 0 res m ẋi i + w i, i + ( ) φ g ij L ξ g(l ij ξ g) = F i, ẋ i + w i, =0 g (4.9) where we ssume h F i m i i = ξ (F i m i i ) [9]. Now he nonrivil sk is o simplify he secon erm on he lef hn sie of Eq. (4.9). Noe h given wo poins in Riemnnin mnifol he geoesic n is lengh boh explicily epen on he meric g. Noe lso h ( ) φ ij L ξ g(l ij ξ g) =0 = φ ij g ij ( ) L ξ g(l ij ξ g). (4.0) =0 Before proceeing ny furher, le us firs simplify he blnce of energy for he originl frme, i.e., he lef-hn sie of Eq. (4.5) φ ij ( Lg (l ij g ) ) = φ ij g ij L g(l ij g ). (4.) Consier he geoesic joining he poins x i (), x j () S. As hese poins re ime epenen, he geoesic joining hem woul be ime epenen s well, i.e., l = l(, s) := l ij g (, s), where s is he curve prmeer n is ime. Noe h L g(l ij g ) = l s, l s s, (4.) where l s = sl(, s) is he velociy of he prmeerize geoesic. We ssume h he curve l is prmeerize by rc lengh, i.e., i hs uni spee everywhere. Thus l s, l s s = l s, l s s = l s D l s, l s s, (4.3) where D is he covrin erivive long he geoesic n for vecor fiel long he curve l i is efine s D V = l Ṽ, (4.4) where Ṽ is n exension of V o S, is he Riemnnin connecion corresponing o he meric g (see [8] for more eils), n l = l(, s). I is esy o show h [8] D l s = D s l (4.5) s in locl coorine chr {x } hey boh hve he following represenion ( x c D l s = D s l = s + x x b ) s γc b c, (4.6)

12 734 A. Yvri n J. E. Mrsen ZAMP where γb c re Chrisoffel coefficiens of he Riemnnin connecion. Using his propery n he fc h l s = we obin b ( ) L g(l ij g ) = D s l, l s s = s l, l s l, D s l s s. (4.7) Bu becuse l is geoesic we hve D s l s = 0 n hence b L g(l ij g ) = s l, l s s = l, l s b (4.8). Noe h l (, ) = ẋ i, l (b, ) = ẋ j (4.9) n le us enoe he velociy vecors of he geoesic poins x i n x j by l s (, ) = ˆ ij T x is, l s (b, ) = ˆ ji T x j S. (4.0) Therefore L g(l ij g ) = ẋ j,ˆ ji ẋ i,ˆ ij. (4.) Now blnce of energy cn be wrien s φ ij ( ẋj,ˆ ji ẋ i,ˆ ij ) + m i i, ẋ i = F i, ẋ i. (4.) g ij Noe h he pricle ccelerion i is he covrin ime erivive of ẋ i. Le us now look blnce of energy for chnge of frme ξ : S S such h ξ =0 = I. Consier fmily of geoesics l(, s) joining he poins x i () = ξ (x i ()), x j () = ξ (x j ()) S. Moive by (4.8), for fixe, l(, s) is he geoesic wih respec o he meric g = ξ g joining hese wo poins. Assume h he geoesic l(s) = l( 0, s) is prmeerize by rc lengh, i.e., i hs uni spee everywhere. Noe h D s ls (, s) = 0 is sisfie for ech geoesic, where D s is he covrin erivive long he geoesic l(, s) wih respec o he meric g. The ime re of chnge of he lengh of his fmily is simplifie o re L g ( l ij g ) = ls, l s s = l s (, s) D ls, l s s, (4.3) where D is he covrin ime erivive long he curve l(, s) wih respec o he meric g. Noe h sill we hve he relion D ls = D s l s in he locl represenion (4.6) he only ifference woul be he -epennce of he Chrisoffel symbols. Thus L g ( l ij g ) = = l s (, s) l s (, s) D s l, l s ( s l, l s s ) l, D s ls s. (4.4)

13 Vol. 60 (009) Energy blnce invrince for inercing pricle sysems 735 Using he geoesic equion n evluing he bove relion = 0, we obin L g ( l ij g ) = l, l b s. (4.5) =0 Noe h l ( 0, ) = ẋ i + w i, l ( 0, b) = ẋ j + w j, (4.6) where w = ξ. Also l s ( 0, ) = ˆ ij T x is, ls ( 0, b) = ˆ ji T x j S. (4.7) Thus L g ( l ij g ) = ẋ j + w j,ˆ ji ẋ i + w i,ˆ ij. (4.8) =0 Blnce of energy for he new frming ime = 0 res φ ij ( ẋj + w j,ˆ ji ẋ i + w i,ˆ ij ) + m i i, ẋ i + w i g ij = F i, ẋ i + w i. (4.9) Subrcing (4.) from he bove blnce equion n noing h he vecor fiel w is rbirry, we obin he following blnce lw j L φ ij g ij ˆ ij + F i = m i i for ll i L. (4.30) I is seen h ssuming h pricles re embee in Riemnnin mnifol n posuling energy blnce n is invrince uner rbirry spil iffeomorphisms resuls in blnce of liner momenum. In erms of he number of blnce lws, hving Riemnnin mbien spce oes no give us ny new relions. However, efining more generl pirwise inercion we see h inercion force i L ue o he pricle j L is irece long he ngen o he geoesic joining x i, x j S he poin x i n is mgniue is equl o he re of chnge of he poenil energy of i, j L wih respec o he iscree meric g ij. Thus, we cn hink of f ij = φ ij = e (4.3) g ij g ij s iscree Doyle Ericksen formul. Noe h becuse ˆ ij n ˆ ji lie in wo ifferen ngen spces, he relion f ji = f ij is meningless, in generl. However, f ij = f ji. In summry, we hve prove he following proposiion. Proposiion 4.. Assuming h blnce of energy for sysem of pirwisely inercing pricles in n mbien Riemnnin mnifol is spilly covrin is equivlen o blnce of liner momenum.

14 736 A. Yvri n J. E. Mrsen ZAMP The proof of converse of his proposiion is similr o h of he coninuum version [6]. Remrks. () From Eq. (4.30) i is seen h he bckgroun meric eners he blnce of liner momenum only hrough he iscree meric g ij. (b) Noe h Eq. (4.30) is he ynmic version of Eq. () in [7]. Exmple. We show in his exmple h wh we jus erive for generl Riemnnin mnifol is reuce o he clssicl resuls when he mbien spce is Euclien. In his cse he geoesic joining x i, x j S = R n hs he following prmerizion l(s) = bxi x j b + xj x i s s [, b]. (4.3) b Becuse he curve is prmeerize by he rc lengh we mus hve b = x j x i. Thus ˆ ij = xj x i x j x i. (4.33) Also in his cse g ij = r ij := x j x i. Hence he following clssicl blnce of liner momenum is recovere j L φ r ij x i x j x j x i + Fi = m i i for ll i L. (4.34) Noe h in his cse he relion f ji = f ij hols s ngen spce o R n ny poin cn be ienifie wih R n. Remrk. Blnce of liner momenum for pricle sysem on Riemnnin mnifol cn of course lso be obine using Hmilon s principle of les cion which is noher covrin pproch vi Lgrngin mechnics. 5. Concluing remrks This pper suie he connecion beween blnce lws n energy blnce invrince for sysem of inercing pricles. I ws shown h, unlike clssicl elsiciy, posuling invrince of energy blnce uner isomeries of he (Euclien) mbien spce is no enough o obin he blnce lws. Inse, if one posules invrince of energy blnce uner rbirry iffeomorphisms, hen one recovers ll he blnce lws. This shows funmenl ifference beween coninuum n sysem of inercing pricles n cn be ssocie wih he nonlocl nure of inercions in sysem of pricles. Blnce of energy for sysem of pricles embee in Riemnnin mnifol ws lso invesige vi generlize form of pirwise inercions by ssuming

15 Vol. 60 (009) Energy blnce invrince for inercing pricle sysems 737 h he pirwise poenil energy of pricles explicily epens on heir pirwise isnces he lenghs of he geoesics joining hem. This efiniion nurlly reuces o he clssicl noion of pirwise inercions in Euclien spce. Posuling blnce of energy n is spil covrince, shows h one cn obin geomeric version of blnce of liner momenum. For pricle i L, blnce of liner momenum is wrien in he ngen spce of S x i S. I ws observe h in his generl seing, relion like f ji = f ij woul be meningless bu inse one hs he meningful relion f ji = f ij. Defining iscree meric s he pirwise isnces, we showe h f ij = e g ij, (5.) which cn be hough of s iscree version of he Doyle Ericksen formul. References [] Eckr, C., Remrks concerning forces on line efecs, Physicl Review 73 (948), [] Eelen, D. G. B. n Lws, N., Thermoynmics of sysems wih nonlocliy, Archive for Rionl Mechnics n Anlysis 43() (97), [3] Eringen, A. C. n Eelen, D. G. B., Nonlocl elsiciy, Inernionl Journl of Engineering Science 0 (97), [4] Green, A. E. n Rivlin, R. S., On Cuchy s equions of moion, ZAMP 5 (964), [5] Gurin, M. E. n Willims, W. O., On he firs lw of hermoynmics, Archive for Rionl Mechnics n Anlysis 4 (97), [6] Gurin, M. E. n Willims, W. O., On coninuum hermoynmics wih muul boy forces n inernl relion, ZAMP (97), [7] Koni, M. n Sun, T., Snr relizions of crysl lices vi hrmonic mps, Trnscions of he Americn Mhemicl Sociey 353 (000), 0. [8] Lee, J. M., Riemnnin Mnifol An Inroucion o Curvure, Springer-Verlg, New York 997. [9] Mrsen, J. E. n T. J. R. Hughes, Mhemicl Founions of Elsiciy, Dover, New York 983. [0] Mrsen, J. E. n T. Riu, Inroucion o Mechnics n Symmery, Springer, New York 003. [] Rjgopl, K. R. n A. R. Srinivs, On he hermomechnics of merils h hve muliple nurl configurions Pr II: Twinning n soli o soli phse rnsformion, ZAMP 55(6) (984), [] Shchepeilov, A. V., Reucion of he wo-boy problem wih cenrl inercion on simply connece spces of consn secionl curvure. Journl of Physics A-Mhemicl n Generl 3 (998), [3] Shchepeilov, A. V., Clculus n Mechnics on Two-Poin Homogenous Riemnnin Spces (Lecure Noes in Physics), Springer, New York 983. [4] Simo, J. C. n J. E. Mrsen, On he roe sress ensor n he meril version of he Doyle Ericksen formul, Archive for Rionl Mechnics n Anlysis 86 (984), 3 3. [5] Yno, T., W.-S. Koon, J.E. Mrsen, n I. Kevrekiis, Gyrion-rius ynmics in srucurl rnsiions of omic clusers, Journl of Chemicl Physics 6 (007), 7. [6] Yvri, A., J. E. Mrsen n M. Oriz, On he spil n meril covrin blnce lws in elsiciy, Journl of Mhemicl Physics 47 (006), 85.

16 738 A. Yvri n J. E. Mrsen ZAMP [7] Yvri, A., M. Oriz n K. Bhchry, A heory of nhrmonic lice sics for nlysis of efecive crysls, Journl of Elsiciy 86 (007), [8] Yvri, A., On geomeric iscreizion of elsiciy, Journl of Mhemicl Physics 49 (008), 36. [9] Yvri, A. n A. Ozkin, Covrince in linerize elsiciy, Zeischrif für Angewne Mhemik un Physik (ZAMP) (008), DOI 0.007/s Arsh Yvri School of Civil n Environmenl Engineering Georgi Insiue of Technology Aln, GA 3033 USA e-mil: rsh.yvri@ce.gech.eu Jerrol E. Mrsen Conrol n Dynmicl Sysems Cliforni Insiue of Technology Psen, CA 95 USA (Receive: My 9, 008) Publishe Online Firs: Ocober 0, 008 To ccess his journl online:

e t dt e t dt = lim e t dt T (1 e T ) = 1

e t dt e t dt = lim e t dt T (1 e T ) = 1 Improper Inegrls There re wo ypes of improper inegrls - hose wih infinie limis of inegrion, nd hose wih inegrnds h pproch some poin wihin he limis of inegrion. Firs we will consider inegrls wih infinie

More information

September 20 Homework Solutions

September 20 Homework Solutions College of Engineering nd Compuer Science Mechnicl Engineering Deprmen Mechnicl Engineering A Seminr in Engineering Anlysis Fll 7 Number 66 Insrucor: Lrry Creo Sepember Homework Soluions Find he specrum

More information

Contraction Mapping Principle Approach to Differential Equations

Contraction Mapping Principle Approach to Differential Equations epl Journl of Science echnology 0 (009) 49-53 Conrcion pping Principle pproch o Differenil Equions Bishnu P. Dhungn Deprmen of hemics, hendr Rn Cmpus ribhuvn Universiy, Khmu epl bsrc Using n eension of

More information

An object moving with speed v around a point at distance r, has an angular velocity. m/s m

An object moving with speed v around a point at distance r, has an angular velocity. m/s m Roion The mosphere roes wih he erh n moions wihin he mosphere clerly follow cure phs (cyclones, nicyclones, hurricnes, ornoes ec.) We nee o epress roion quniiely. For soli objec or ny mss h oes no isor

More information

14. The fundamental theorem of the calculus

14. The fundamental theorem of the calculus 4. The funmenl heorem of he clculus V 20 00 80 60 40 20 0 0 0.2 0.4 0.6 0.8 v 400 200 0 0 0.2 0.5 0.8 200 400 Figure : () Venriculr volume for subjecs wih cpciies C = 24 ml, C = 20 ml, C = 2 ml n (b) he

More information

4.8 Improper Integrals

4.8 Improper Integrals 4.8 Improper Inegrls Well you ve mde i hrough ll he inegrion echniques. Congrs! Unforunely for us, we sill need o cover one more inegrl. They re clled Improper Inegrls. A his poin, we ve only del wih inegrls

More information

Physics 100: Lecture 1

Physics 100: Lecture 1 Physics : Lecure Agen for Toy Aice Scope of his course Mesuremen n Unis Funmenl unis Sysems of unis Conering beween sysems of unis Dimensionl Anlysis -D Kinemics (reiew) Aerge & insnneous elociy n ccelerion

More information

( ) 2 a b ab. To do this, we are to use the Ricci identity (which we use to evaluate the RHS) and the properties of the Lie derivative.

( ) 2 a b ab. To do this, we are to use the Ricci identity (which we use to evaluate the RHS) and the properties of the Lie derivative. Exercise [9.6] This exercise sks s o show h he ccelerion of n (infiniesiml volme mesre V long he worlline he volme s cener e o he effecs of spceime crvre is given by: D V = R V ( b b To o his, we re o

More information

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6. [~ o o :- o o ill] i 1. Mrices, Vecors, nd Guss-Jordn Eliminion 1 x y = = - z= The soluion is ofen represened s vecor: n his exmple, he process of eliminion works very smoohly. We cn elimine ll enries

More information

0 for t < 0 1 for t > 0

0 for t < 0 1 for t > 0 8.0 Sep nd del funcions Auhor: Jeremy Orloff The uni Sep Funcion We define he uni sep funcion by u() = 0 for < 0 for > 0 I is clled he uni sep funcion becuse i kes uni sep = 0. I is someimes clled he Heviside

More information

Average & instantaneous velocity and acceleration Motion with constant acceleration

Average & instantaneous velocity and acceleration Motion with constant acceleration Physics 7: Lecure Reminders Discussion nd Lb secions sr meeing ne week Fill ou Pink dd/drop form if you need o swich o differen secion h is FULL. Do i TODAY. Homework Ch. : 5, 7,, 3,, nd 6 Ch.: 6,, 3 Submission

More information

f t f a f x dx By Lin McMullin f x dx= f b f a. 2

f t f a f x dx By Lin McMullin f x dx= f b f a. 2 Accumulion: Thoughs On () By Lin McMullin f f f d = + The gols of he AP* Clculus progrm include he semen, Sudens should undersnd he definie inegrl s he ne ccumulion of chnge. 1 The Topicl Ouline includes

More information

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration.

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration. Moion Accelerion Pr : Consn Accelerion Accelerion Accelerion Accelerion is he re of chnge of velociy. = v - vo = Δv Δ ccelerion = = v - vo chnge of velociy elpsed ime Accelerion is vecor, lhough in one-dimensionl

More information

More on Magnetically C Coupled Coils and Ideal Transformers

More on Magnetically C Coupled Coils and Ideal Transformers Appenix ore on gneiclly C Couple Coils Iel Trnsformers C. Equivlen Circuis for gneiclly Couple Coils A imes, i is convenien o moel mgneiclly couple coils wih n equivlen circui h oes no involve mgneic coupling.

More information

A Kalman filtering simulation

A Kalman filtering simulation A Klmn filering simulion The performnce of Klmn filering hs been esed on he bsis of wo differen dynmicl models, ssuming eiher moion wih consn elociy or wih consn ccelerion. The former is epeced o beer

More information

Physics Worksheet Lesson 4: Linear Motion Section: Name:

Physics Worksheet Lesson 4: Linear Motion Section: Name: Physics Workshee Lesson 4: Liner Moion Secion: Nme: 1. Relie Moion:. All moion is. b. is n rbirry coorine sysem wih reference o which he posiion or moion of somehing is escribe or physicl lws re formule.

More information

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment Mgneosics Br Mgne As fr bck s 4500 yers go, he Chinese discovered h cerin ypes of iron ore could rc ech oher nd cerin mels. Iron filings "mp" of br mgne s field Crefully suspended slivers of his mel were

More information

Mathematical Modeling

Mathematical Modeling ME pplie Engineering nlsis Chper Mhemicl Moeling Professor Ti-Rn Hsu, Ph.D. Deprmen of Mechnicl n erospce Engineering Sn Jose Se Universi Sn Jose, Cliforni, US Jnur Chper Lerning Ojecives Mhemicl moeling

More information

Characteristic Function for the Truncated Triangular Distribution., Myron Katzoff and Rahul A. Parsa

Characteristic Function for the Truncated Triangular Distribution., Myron Katzoff and Rahul A. Parsa Secion on Survey Reserch Mehos JSM 009 Chrcerisic Funcion for he Trunce Tringulr Disriuion Jy J. Kim 1 1, Myron Kzoff n Rhul A. Prs 1 Nionl Cener for Helh Sisics, 11Toleo Ro, Hysville, MD. 078 College

More information

3D Transformations. Computer Graphics COMP 770 (236) Spring Instructor: Brandon Lloyd 1/26/07 1

3D Transformations. Computer Graphics COMP 770 (236) Spring Instructor: Brandon Lloyd 1/26/07 1 D Trnsformions Compuer Grphics COMP 770 (6) Spring 007 Insrucor: Brndon Lloyd /6/07 Geomery Geomeric eniies, such s poins in spce, exis wihou numers. Coordines re nming scheme. The sme poin cn e descried

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > for ll smples y i solve sysem of liner inequliies MSE procedure y i = i for ll smples

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > 0 for ll smples y i solve sysem of liner inequliies MSE procedure y i i for ll smples

More information

Think of the Relationship Between Time and Space Again

Think of the Relationship Between Time and Space Again Repor nd Opinion, 1(3),009 hp://wwwsciencepubne sciencepub@gmilcom Think of he Relionship Beween Time nd Spce Agin Yng F-cheng Compny of Ruid Cenre in Xinjing 15 Hongxing Sree, Klmyi, Xingjing 834000,

More information

PHYSICS 1210 Exam 1 University of Wyoming 14 February points

PHYSICS 1210 Exam 1 University of Wyoming 14 February points PHYSICS 1210 Em 1 Uniersiy of Wyoming 14 Februry 2013 150 poins This es is open-noe nd closed-book. Clculors re permied bu compuers re no. No collborion, consulion, or communicion wih oher people (oher

More information

Chapter Three Systems of Linear Differential Equations

Chapter Three Systems of Linear Differential Equations Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n

More information

Probability, Estimators, and Stationarity

Probability, Estimators, and Stationarity Chper Probbiliy, Esimors, nd Sionriy Consider signl genered by dynmicl process, R, R. Considering s funcion of ime, we re opering in he ime domin. A fundmenl wy o chrcerize he dynmics using he ime domin

More information

Chapter 2. Motion along a straight line. 9/9/2015 Physics 218

Chapter 2. Motion along a straight line. 9/9/2015 Physics 218 Chper Moion long srigh line 9/9/05 Physics 8 Gols for Chper How o describe srigh line moion in erms of displcemen nd erge elociy. The mening of insnneous elociy nd speed. Aerge elociy/insnneous elociy

More information

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008)

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008) MATH 14 AND 15 FINAL EXAM REVIEW PACKET (Revised spring 8) The following quesions cn be used s review for Mh 14/ 15 These quesions re no cul smples of quesions h will pper on he finl em, bu hey will provide

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Swiss Federl Insiue of Pge 1 The Finie Elemen Mehod for he Anlysis of Non-Liner nd Dynmic Sysems Prof. Dr. Michel Hvbro Fber Dr. Nebojs Mojsilovic Swiss Federl Insiue of ETH Zurich, Swizerlnd Mehod of

More information

P441 Analytical Mechanics - I. Coupled Oscillators. c Alex R. Dzierba

P441 Analytical Mechanics - I. Coupled Oscillators. c Alex R. Dzierba Lecure 3 Mondy - Deceber 5, 005 Wrien or ls upded: Deceber 3, 005 P44 Anlyicl Mechnics - I oupled Oscillors c Alex R. Dzierb oupled oscillors - rix echnique In Figure we show n exple of wo coupled oscillors,

More information

Some Inequalities variations on a common theme Lecture I, UL 2007

Some Inequalities variations on a common theme Lecture I, UL 2007 Some Inequliies vriions on common heme Lecure I, UL 2007 Finbrr Hollnd, Deprmen of Mhemics, Universiy College Cork, fhollnd@uccie; July 2, 2007 Three Problems Problem Assume i, b i, c i, i =, 2, 3 re rel

More information

3. Renewal Limit Theorems

3. Renewal Limit Theorems Virul Lborories > 14. Renewl Processes > 1 2 3 3. Renewl Limi Theorems In he inroducion o renewl processes, we noed h he rrivl ime process nd he couning process re inverses, in sens The rrivl ime process

More information

( ) ( ) ( ) ( ) ( ) ( y )

( ) ( ) ( ) ( ) ( ) ( y ) 8. Lengh of Plne Curve The mos fmous heorem in ll of mhemics is he Pyhgoren Theorem. I s formulion s he disnce formul is used o find he lenghs of line segmens in he coordine plne. In his secion you ll

More information

Mathematics 805 Final Examination Answers

Mathematics 805 Final Examination Answers . 5 poins Se he Weiersrss M-es. Mhemics 85 Finl Eminion Answers Answer: Suppose h A R, nd f n : A R. Suppose furher h f n M n for ll A, nd h Mn converges. Then f n converges uniformly on A.. 5 poins Se

More information

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q).

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q). INTEGRALS JOHN QUIGG Eercise. Le f : [, b] R be bounded, nd le P nd Q be priions of [, b]. Prove h if P Q hen U(P ) U(Q) nd L(P ) L(Q). Soluion: Le P = {,..., n }. Since Q is obined from P by dding finiely

More information

GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS

GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS - TAMKANG JOURNAL OF MATHEMATICS Volume 5, Number, 7-5, June doi:5556/jkjm555 Avilble online hp://journlsmhkueduw/ - - - GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS MARCELA

More information

MTH 146 Class 11 Notes

MTH 146 Class 11 Notes 8.- Are of Surfce of Revoluion MTH 6 Clss Noes Suppose we wish o revolve curve C round n is nd find he surfce re of he resuling solid. Suppose f( ) is nonnegive funcion wih coninuous firs derivive on he

More information

REAL ANALYSIS I HOMEWORK 3. Chapter 1

REAL ANALYSIS I HOMEWORK 3. Chapter 1 REAL ANALYSIS I HOMEWORK 3 CİHAN BAHRAN The quesions re from Sein nd Shkrchi s e. Chper 1 18. Prove he following sserion: Every mesurble funcion is he limi.e. of sequence of coninuous funcions. We firs

More information

Location is relative. Coordinate Systems. Which of the following can be described with vectors??

Location is relative. Coordinate Systems. Which of the following can be described with vectors?? Locion is relive Coordine Sysems The posiion o hing is sed relive o noher hing (rel or virul) review o he physicl sis h governs mhemicl represenions Reerence oec mus e deined Disnce mus e nown Direcion

More information

Solutions to Problems from Chapter 2

Solutions to Problems from Chapter 2 Soluions o Problems rom Chper Problem. The signls u() :5sgn(), u () :5sgn(), nd u h () :5sgn() re ploed respecively in Figures.,b,c. Noe h u h () :5sgn() :5; 8 including, bu u () :5sgn() is undeined..5

More information

Chapter Direct Method of Interpolation

Chapter Direct Method of Interpolation Chper 5. Direc Mehod of Inerpolion Afer reding his chper, you should be ble o:. pply he direc mehod of inerpolion,. sole problems using he direc mehod of inerpolion, nd. use he direc mehod inerpolns o

More information

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples.

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples. Improper Inegrls To his poin we hve only considered inegrls f(x) wih he is of inegrion nd b finie nd he inegrnd f(x) bounded (nd in fc coninuous excep possibly for finiely mny jump disconinuiies) An inegrl

More information

can be viewed as a generalized product, and one for which the product of f and g. That is, does

can be viewed as a generalized product, and one for which the product of f and g. That is, does Boyce/DiPrim 9 h e, Ch 6.6: The Convoluion Inegrl Elemenry Differenil Equion n Bounry Vlue Problem, 9 h eiion, by Willim E. Boyce n Richr C. DiPrim, 9 by John Wiley & Son, Inc. Someime i i poible o wrie

More information

3 Motion with constant acceleration: Linear and projectile motion

3 Motion with constant acceleration: Linear and projectile motion 3 Moion wih consn ccelerion: Liner nd projecile moion cons, In he precedin Lecure we he considered moion wih consn ccelerion lon he is: Noe h,, cn be posiie nd neie h leds o rie of behiors. Clerl similr

More information

HW #1 Solutions. Lewis Structures: Using the above rules, determine the molecular structure for Cl2CO. Hint: C is at the center.

HW #1 Solutions. Lewis Structures: Using the above rules, determine the molecular structure for Cl2CO. Hint: C is at the center. HW # Soluions Cron Mss Prolem: ssuming n erge surfce pressure of m, n erge ropospheric emperure of 55 K, n glol CO mixing rio of 385 ppm, wh is he curren mospheric Cron reseroir (in unis of g m -? Compre

More information

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang jordnmcd Eigenvlue-eigenvecor pproch o solving firs order ODEs -- ordn norml (cnonicl) form Insrucor: Nm Sun Wng Consider he following se of coupled firs order ODEs d d x x 5 x x d d x d d x x x 5 x x

More information

RESPONSE UNDER A GENERAL PERIODIC FORCE. When the external force F(t) is periodic with periodτ = 2π

RESPONSE UNDER A GENERAL PERIODIC FORCE. When the external force F(t) is periodic with periodτ = 2π RESPONSE UNDER A GENERAL PERIODIC FORCE When he exernl force F() is periodic wih periodτ / ω,i cn be expnded in Fourier series F( ) o α ω α b ω () where τ F( ) ω d, τ,,,... () nd b τ F( ) ω d, τ,,... (3)

More information

D.I. Survival models and copulas

D.I. Survival models and copulas D- D. SURVIVAL COPULA D.I. Survival moels an copulas Definiions, relaionships wih mulivariae survival isribuion funcions an relaionships beween copulas an survival copulas. D.II. Fraily moels Use of a

More information

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function ENGR 1990 Engineering Mhemics The Inegrl of Funcion s Funcion Previously, we lerned how o esime he inegrl of funcion f( ) over some inervl y dding he res of finie se of rpezoids h represen he re under

More information

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points)

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points) Delf Universiy of Technology Fculy of Civil Engineering nd Geosciences Srucurl echnics Secion Wrie your nme nd sudy numer he op righ-hnd of your work. Exm CT5 Theory of Elsiciy Fridy Ocoer 00, 9:00 :00

More information

Version 001 test-1 swinney (57010) 1. is constant at m/s.

Version 001 test-1 swinney (57010) 1. is constant at m/s. Version 001 es-1 swinne (57010) 1 This prin-ou should hve 20 quesions. Muliple-choice quesions m coninue on he nex column or pge find ll choices before nswering. CubeUniVec1x76 001 10.0 poins Acubeis1.4fee

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method IOSR Journl of Mhemics (IOSR-JM) e-issn: 78-578. Volume 5, Issue 3 (Jn. - Feb. 13), PP 6-11 Soluions for Nonliner Pril Differenil Equions By Tn-Co Mehod Mhmood Jwd Abdul Rsool Abu Al-Sheer Al -Rfidin Universiy

More information

Procedia Computer Science

Procedia Computer Science Procedi Compuer Science 00 (0) 000 000 Procedi Compuer Science www.elsevier.com/loce/procedi The Third Informion Sysems Inernionl Conference The Exisence of Polynomil Soluion of he Nonliner Dynmicl Sysems

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

NMR Spectroscopy: Principles and Applications. Nagarajan Murali Advanced Tools Lecture 4

NMR Spectroscopy: Principles and Applications. Nagarajan Murali Advanced Tools Lecture 4 NMR Specroscop: Principles nd Applicions Ngrjn Murli Advnced Tools Lecure 4 Advnced Tools Qunum Approch We know now h NMR is rnch of Specroscop nd he MNR specrum is n oucome of nucler spin inercion wih

More information

Conservation Law. Chapter Goal. 6.2 Theory

Conservation Law. Chapter Goal. 6.2 Theory Chpter 6 Conservtion Lw 6.1 Gol Our long term gol is to unerstn how mthemticl moels re erive. Here, we will stuy how certin quntity chnges with time in given region (sptil omin). We then first erive the

More information

Convergence of Singular Integral Operators in Weighted Lebesgue Spaces

Convergence of Singular Integral Operators in Weighted Lebesgue Spaces EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 10, No. 2, 2017, 335-347 ISSN 1307-5543 www.ejpm.com Published by New York Business Globl Convergence of Singulr Inegrl Operors in Weighed Lebesgue

More information

Introduction to continuum mechanics-ii

Introduction to continuum mechanics-ii Ioan R. IONESCU LSPM, Universiy Paris 3 (Sorbonne-Paris-Cié) July 3, 27 Coninuum kinemaics B uneforme/reference configuraion, B eforme/acual configuraion Moion: x(, X ) = X + u(, X ) wih x(, B ) = B. Maerial/Lagrangian

More information

LAGRANGIAN AND HAMILTONIAN MECHANICS WITH FRACTIONAL DERIVATIVES

LAGRANGIAN AND HAMILTONIAN MECHANICS WITH FRACTIONAL DERIVATIVES LAGRANGIAN AND HAMILTONIAN MEHANIS WITH FRATIONAL DERIVATIVES EMIL POPESU 2,1 1 Asronomicl Insiue of Romnin Acdemy Sr uiul de Argin 5, 40557 Buchres, Romni 2 Technicl Universiy of ivil Engineering, Bd

More information

SOME USEFUL MATHEMATICS

SOME USEFUL MATHEMATICS SOME USEFU MAHEMAICS SOME USEFU MAHEMAICS I is esy o mesure n preic he behvior of n elecricl circui h conins only c volges n currens. However, mos useful elecricl signls h crry informion vry wih ime. Since

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak .65, MHD Theory of Fusion Sysems Prof. Freidberg Lecure 9: The High e Tokmk Summry of he Properies of n Ohmic Tokmk. Advnges:. good euilibrium (smll shif) b. good sbiliy ( ) c. good confinemen ( τ nr )

More information

Physic 231 Lecture 4. Mi it ftd l t. Main points of today s lecture: Example: addition of velocities Trajectories of objects in 2 = =

Physic 231 Lecture 4. Mi it ftd l t. Main points of today s lecture: Example: addition of velocities Trajectories of objects in 2 = = Mi i fd l Phsic 3 Lecure 4 Min poins of od s lecure: Emple: ddiion of elociies Trjecories of objecs in dimensions: dimensions: g 9.8m/s downwrds ( ) g o g g Emple: A foobll pler runs he pern gien in he

More information

Neural assembly binding in linguistic representation

Neural assembly binding in linguistic representation Neurl ssembly binding in linguisic represenion Frnk vn der Velde & Mrc de Kmps Cogniive Psychology Uni, Universiy of Leiden, Wssenrseweg 52, 2333 AK Leiden, The Neherlnds, vdvelde@fsw.leidenuniv.nl Absrc.

More information

Integral representations and new generating functions of Chebyshev polynomials

Integral representations and new generating functions of Chebyshev polynomials Inegral represenaions an new generaing funcions of Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 186 Roma, Ialy email:

More information

ME 141. Engineering Mechanics

ME 141. Engineering Mechanics ME 141 Engineeing Mechnics Lecue 13: Kinemics of igid bodies hmd Shhedi Shkil Lecue, ep. of Mechnicl Engg, UET E-mil: sshkil@me.bue.c.bd, shkil6791@gmil.com Websie: eche.bue.c.bd/sshkil Couesy: Veco Mechnics

More information

FM Applications of Integration 1.Centroid of Area

FM Applications of Integration 1.Centroid of Area FM Applicions of Inegrion.Cenroid of Are The cenroid of ody is is geomeric cenre. For n ojec mde of uniform meril, he cenroid coincides wih he poin which he ody cn e suppored in perfecly lnced se ie, is

More information

5.1-The Initial-Value Problems For Ordinary Differential Equations

5.1-The Initial-Value Problems For Ordinary Differential Equations 5.-The Iniil-Vlue Problems For Ordinry Differenil Equions Consider solving iniil-vlue problems for ordinry differenil equions: (*) y f, y, b, y. If we know he generl soluion y of he ordinry differenil

More information

Transforms II - Wavelets Preliminary version please report errors, typos, and suggestions for improvements

Transforms II - Wavelets Preliminary version please report errors, typos, and suggestions for improvements EECS 3 Digil Signl Processing Universiy of Cliforni, Berkeley: Fll 007 Gspr November 4, 007 Trnsforms II - Wveles Preliminry version plese repor errors, ypos, nd suggesions for improvemens We follow n

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

1.0 Electrical Systems

1.0 Electrical Systems . Elecricl Sysems The ypes of dynmicl sysems we will e sudying cn e modeled in erms of lgeric equions, differenil equions, or inegrl equions. We will egin y looking fmilir mhemicl models of idel resisors,

More information

5.4, 6.1, 6.2 Handout. As we ve discussed, the integral is in some way the opposite of taking a derivative. The exact relationship

5.4, 6.1, 6.2 Handout. As we ve discussed, the integral is in some way the opposite of taking a derivative. The exact relationship 5.4, 6.1, 6.2 Hnout As we ve iscusse, the integrl is in some wy the opposite of tking erivtive. The exct reltionship is given by the Funmentl Theorem of Clculus: The Funmentl Theorem of Clculus: If f is

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Motion in a Straight Line

Motion in a Straight Line Moion in Srigh Line. Preei reched he mero sion nd found h he esclor ws no working. She wlked up he sionry esclor in ime. On oher dys, if she remins sionry on he moing esclor, hen he esclor kes her up in

More information

A Simple Method to Solve Quartic Equations. Key words: Polynomials, Quartics, Equations of the Fourth Degree INTRODUCTION

A Simple Method to Solve Quartic Equations. Key words: Polynomials, Quartics, Equations of the Fourth Degree INTRODUCTION Ausrlin Journl of Bsic nd Applied Sciences, 6(6): -6, 0 ISSN 99-878 A Simple Mehod o Solve Quric Equions Amir Fhi, Poo Mobdersn, Rhim Fhi Deprmen of Elecricl Engineering, Urmi brnch, Islmic Ad Universi,

More information

Hermite-Hadamard-Fejér type inequalities for convex functions via fractional integrals

Hermite-Hadamard-Fejér type inequalities for convex functions via fractional integrals Sud. Univ. Beş-Bolyi Mh. 6(5, No. 3, 355 366 Hermie-Hdmrd-Fejér ype inequliies for convex funcions vi frcionl inegrls İmd İşcn Asrc. In his pper, firsly we hve eslished Hermie Hdmrd-Fejér inequliy for

More information

EULER-LAGRANGE EQUATIONS. Contents. 2. Variational formulation 2 3. Constrained systems and d Alembert principle Legendre transform 6

EULER-LAGRANGE EQUATIONS. Contents. 2. Variational formulation 2 3. Constrained systems and d Alembert principle Legendre transform 6 EULER-LAGRANGE EQUATIONS EUGENE LERMAN Contents 1. Clssicl system of N prticles in R 3 1 2. Vritionl formultion 2 3. Constrine systems n Alembert principle. 4 4. Legenre trnsform 6 1. Clssicl system of

More information

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM Elecronic Journl of Differenil Equions, Vol. 208 (208), No. 50, pp. 6. ISSN: 072-669. URL: hp://ejde.mh.xse.edu or hp://ejde.mh.un.edu EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE

More information

How to prove the Riemann Hypothesis

How to prove the Riemann Hypothesis Scholrs Journl of Phsics, Mhemics nd Sisics Sch. J. Phs. Mh. S. 5; (B:5-6 Scholrs Acdemic nd Scienific Publishers (SAS Publishers (An Inernionl Publisher for Acdemic nd Scienific Resources *Corresonding

More information

Hamilton-Jacobi irrotational hydrodynamics for binary neutron star inspiral

Hamilton-Jacobi irrotational hydrodynamics for binary neutron star inspiral Hmilon-Jcoi irroionl hydrodynmics for inry neuron sr inspirl Chrlmpos Mrkkis Mhemicl Sciences, Universiy of Souhmpon work in progress in collorion wih: John Friedmn, Msru Shi, Nils Andersson, Nicls Moldenhuer,

More information

LECTURES ON RECONSTRUCTION ALGEBRAS I

LECTURES ON RECONSTRUCTION ALGEBRAS I LETURES ON REONSTRUTION ALGEBRAS I MIHAEL WEMYSS. Inroducion Noncommuive lgebr (=quivers) cn be used o solve boh explici nd non-explici problems in lgebric geomery, nd hese lecures will ry o explin some

More information

Exact Minimization of # of Joins

Exact Minimization of # of Joins A Quer Rewriing Algorihm: Ec Minimizion of # of Joins Emple (movie bse) selec.irecor from movie, movie, movie m3, scheule, scheule s2 where.irecor =.irecor n.cor = m3.cor n.ile =.ile n m3.ile = s2.ile

More information

A 1.3 m 2.5 m 2.8 m. x = m m = 8400 m. y = 4900 m 3200 m = 1700 m

A 1.3 m 2.5 m 2.8 m. x = m m = 8400 m. y = 4900 m 3200 m = 1700 m PHYS : Soluions o Chper 3 Home Work. SSM REASONING The displcemen is ecor drwn from he iniil posiion o he finl posiion. The mgniude of he displcemen is he shores disnce beween he posiions. Noe h i is onl

More information

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x)

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x) Properies of Logrihms Solving Eponenil nd Logrihmic Equions Properies of Logrihms Produc Rule ( ) log mn = log m + log n ( ) log = log + log Properies of Logrihms Quoien Rule log m = logm logn n log7 =

More information

MAT 266 Calculus for Engineers II Notes on Chapter 6 Professor: John Quigg Semester: spring 2017

MAT 266 Calculus for Engineers II Notes on Chapter 6 Professor: John Quigg Semester: spring 2017 MAT 66 Clculus for Engineers II Noes on Chper 6 Professor: John Quigg Semeser: spring 7 Secion 6.: Inegrion by prs The Produc Rule is d d f()g() = f()g () + f ()g() Tking indefinie inegrls gives [f()g

More information

PARABOLA. moves such that PM. = e (constant > 0) (eccentricity) then locus of P is called a conic. or conic section.

PARABOLA. moves such that PM. = e (constant > 0) (eccentricity) then locus of P is called a conic. or conic section. wwwskshieducioncom PARABOLA Le S be given fixed poin (focus) nd le l be given fixed line (Direcrix) Le SP nd PM be he disnce of vrible poin P o he focus nd direcrix respecively nd P SP moves such h PM

More information

A new model for limit order book dynamics

A new model for limit order book dynamics Anewmodelforlimiorderbookdynmics JeffreyR.Russell UniversiyofChicgo,GrdueSchoolofBusiness TejinKim UniversiyofChicgo,DeprmenofSisics Absrc:Thispperproposesnewmodelforlimiorderbookdynmics.Thelimiorderbookconsiss

More information

THREE IMPORTANT CONCEPTS IN TIME SERIES ANALYSIS: STATIONARITY, CROSSING RATES, AND THE WOLD REPRESENTATION THEOREM

THREE IMPORTANT CONCEPTS IN TIME SERIES ANALYSIS: STATIONARITY, CROSSING RATES, AND THE WOLD REPRESENTATION THEOREM THR IMPORTANT CONCPTS IN TIM SRIS ANALYSIS: STATIONARITY, CROSSING RATS, AND TH WOLD RPRSNTATION THORM Prof. Thoms B. Fomb Deprmen of conomics Souhern Mehodis Universi June 8 I. Definiion of Covrince Sionri

More information

Section 5: Chain Rule

Section 5: Chain Rule Chaper The Derivaive Applie Calculus 11 Secion 5: Chain Rule There is one more ype of complicae funcion ha we will wan o know how o iffereniae: composiion. The Chain Rule will le us fin he erivaive of

More information

1. Introduction. 1 b b

1. Introduction. 1 b b Journl of Mhemicl Inequliies Volume, Number 3 (007), 45 436 SOME IMPROVEMENTS OF GRÜSS TYPE INEQUALITY N. ELEZOVIĆ, LJ. MARANGUNIĆ AND J. PEČARIĆ (communiced b A. Čižmešij) Absrc. In his pper some inequliies

More information

1. Consider a PSA initially at rest in the beginning of the left-hand end of a long ISS corridor. Assume xo = 0 on the left end of the ISS corridor.

1. Consider a PSA initially at rest in the beginning of the left-hand end of a long ISS corridor. Assume xo = 0 on the left end of the ISS corridor. In Eercise 1, use sndrd recngulr Cresin coordine sysem. Le ime be represened long he horizonl is. Assume ll ccelerions nd decelerions re consn. 1. Consider PSA iniilly res in he beginning of he lef-hnd

More information

How to Prove the Riemann Hypothesis Author: Fayez Fok Al Adeh.

How to Prove the Riemann Hypothesis Author: Fayez Fok Al Adeh. How o Prove he Riemnn Hohesis Auhor: Fez Fok Al Adeh. Presiden of he Srin Cosmologicl Socie P.O.Bo,387,Dmscus,Sri Tels:963--77679,735 Emil:hf@scs-ne.org Commens: 3 ges Subj-Clss: Funcionl nlsis, comle

More information

PART V. Wavelets & Multiresolution Analysis

PART V. Wavelets & Multiresolution Analysis Wveles 65 PART V Wveles & Muliresoluion Anlysis ADDITIONAL REFERENCES: A. Cohen, Numericl Anlysis o Wvele Mehods, Norh-Hollnd, (003) S. Mll, A Wvele Tour o Signl Processing, Acdemic Press, (999) I. Dubechies,

More information

Application on Inner Product Space with. Fixed Point Theorem in Probabilistic

Application on Inner Product Space with. Fixed Point Theorem in Probabilistic Journl of Applied Mhemics & Bioinformics, vol.2, no.2, 2012, 1-10 ISSN: 1792-6602 prin, 1792-6939 online Scienpress Ld, 2012 Applicion on Inner Produc Spce wih Fixed Poin Theorem in Probbilisic Rjesh Shrivsv

More information

Physics 2A HW #3 Solutions

Physics 2A HW #3 Solutions Chper 3 Focus on Conceps: 3, 4, 6, 9 Problems: 9, 9, 3, 41, 66, 7, 75, 77 Phsics A HW #3 Soluions Focus On Conceps 3-3 (c) The ccelerion due o grvi is he sme for boh blls, despie he fc h he hve differen

More information

(b) 10 yr. (b) 13 m. 1.6 m s, m s m s (c) 13.1 s. 32. (a) 20.0 s (b) No, the minimum distance to stop = 1.00 km. 1.

(b) 10 yr. (b) 13 m. 1.6 m s, m s m s (c) 13.1 s. 32. (a) 20.0 s (b) No, the minimum distance to stop = 1.00 km. 1. Answers o Een Numbered Problems Chper. () 7 m s, 6 m s (b) 8 5 yr 4.. m ih 6. () 5. m s (b).5 m s (c).5 m s (d) 3.33 m s (e) 8. ().3 min (b) 64 mi..3 h. ().3 s (b) 3 m 4..8 mi wes of he flgpole 6. (b)

More information

Faraday s Law. To be able to find. motional emf transformer and motional emf. Motional emf

Faraday s Law. To be able to find. motional emf transformer and motional emf. Motional emf Objecie F s w Tnsfome Moionl To be ble o fin nsfome. moionl nsfome n moionl. 331 1 331 Mwell s quion: ic Fiel D: Guss lw :KV : Guss lw H: Ampee s w Poin Fom Inegl Fom D D Q sufce loop H sufce H I enclose

More information

A LIMIT-POINT CRITERION FOR A SECOND-ORDER LINEAR DIFFERENTIAL OPERATOR IAN KNOWLES

A LIMIT-POINT CRITERION FOR A SECOND-ORDER LINEAR DIFFERENTIAL OPERATOR IAN KNOWLES A LIMIT-POINT CRITERION FOR A SECOND-ORDER LINEAR DIFFERENTIAL OPERATOR j IAN KNOWLES 1. Inroducion Consider he forml differenil operor T defined by el, (1) where he funcion q{) is rel-vlued nd loclly

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

Chapter 2: Evaluative Feedback

Chapter 2: Evaluative Feedback Chper 2: Evluive Feedbck Evluing cions vs. insrucing by giving correc cions Pure evluive feedbck depends olly on he cion ken. Pure insrucive feedbck depends no ll on he cion ken. Supervised lerning is

More information