Mathematical Applications in Modern Science

Size: px
Start display at page:

Download "Mathematical Applications in Modern Science"

Transcription

1 Effect of Perturbatons n the Corols Centrfugal Forces on the Locaton Stablty of the Equlbrum Solutons n Robe s Restrcted Problem of + Bodes BHAVNEET KAUR Lady Shr Ram College for Women Department of Mathematcs Unversty of Delh INDIA bhavneet.lsr@gmal.com RAJIV AGGARWAL Sr Aurobndo College Department of Mathematcs Unversty of Delh INDIA rajv agg97@yahoo.com Abstract: In ths problem, one of the prmares of mass m s a Roche Ellpsod flled wth a homogeneous ncompressble flud of densty ρ. The smaller prmary of mass m s a pont mass outsde the Ellpsod. The thrd the fourth bodes of mass m m 4 respectvely are small sold spheres of densty ρ ρ 4 respectvely nsde the Ellpsod, wth the assumpton that the mass the radus of the thrd the fourth body are nfntesmal. We assume that m s descrbng a crcle around m. The masses m m 4 mutually attract each other, do not nfluence the motons of m m but are nfluenced by them. In ths paper, we have studed the effect of perturbatons n the corols centrfugal forces on the locaton stablty of the equlbrum ponts n Robe s restrcted problem of + bodes wth the bgger prmary a Roche ellpsod Chrashekhar,987. We have taken nto consderaton all the three components of the pressure feld n dervng the expresson for the buoyancy force vz due to the own gravtatonal feld of the flud that orgnatng n the attracton of m that arsng from the centrfugal force. In ths paper, equlbrum solutons of m m 4 ther lnear stablty are analyzed. We have proved that there exst only sx equlbrum solutons of the system, provded they le wthn the Roche Ellpsod. In a system the prmares are consdered as earth-moon m, m 4 as cebergs, the equlbrum solutons of m m 4 respectvely when the dsplacement s gven n the drecton of x axs or x axs are unstable. Robe s Restrcted Problem; Roche Ellpsod; Equlbrum Soluton; Stablty; Buoyancy Force. Introducton Robe 977 has nvestgated a new knd of restrcted three-body problem n whch one of the prmares of mass m s a rgd sphercal shell flled wth a homogeneous ncompressble flud of densty ρ. The smaller prmary s a mass pont m outsde the shell. The thrd body of mass m, supposed movng nsde the shell, s a small sold sphere of densty ρ, wth the assumpton that the mass the radus of the thrd body are nfntesmal. He further assumed that the mass m descrbes a Kepleran orbt around the mass m. He has proved that the centre of the frst prmary, s the only equlbrum soluton for all values of the densty parameter K, mass parameter µ, eccentrcty parameter e. Further, he has dscussed the lnear stablty of ths equlbrum soluton. He has explctly dscussed two cases. In the frst case, the orbt of m around m s crcular n the second case, the orbt s ellptc, but the shell s empty.e. no flud nsde t or the denstes of m m are equal. In each case, the doman of stablty has been nvestgated for the whole range of parameters occurrng n the problem. In the above problem, the exstence of only one equlbrum soluton namely, centre of the frst prmary s dscussed. Hallan & Rana 00 studed the exstence the lnear stablty of all the equlbrum solutons n the Robe s restrcted three-body problem. They have proved that besdes the centre, there are other equlbrum solutons whch exst only when K 0 the second prmary moves around the frst n a crcular orbt. In addton to the exstng collnear equlbrum soluton, they have shown the exstence of one more whch s stable only under certan restrctons on K µ. There are two trangular equlbrum solutons whch are always unstable. In addton, there are an nfnte number of equlbrum solutons lyng on a crcle centre at the second prmary, provded the solutons are nsde the sphercal shell. The crcular equlbrum solutons are always unstable. ISBN:

2 Szebehely 967 consdered the effect of small perturbatons n the corols force on the stablty of equlbrum ponts n the classcal restrcted problem, keepng the centrfugal force constant proved that the collnear ponts reman unstable. The range of stablty of the trangular ponts ncreases or decreases dependng upon whether the change ɛ n the corols force s postve or negatve thereby establshng that the corols force s a stablzng force. In the same problem, Bhatnagar & Hallan 978 studed the effect of the changes ɛ ɛ n the corols centrfugal forces respectvely on the stablty of equlbrum ponts concluded that the collnear ponts reman unstable for the trangular ponts, the range of stablty ncreases or decreases dependng upon whether the pont ɛ, ɛ le n one or the other of the two parts n whch the lne 6ɛ 9ɛ = 0 dvdes the ɛ, ɛ plane. Shrvastava Garan 99 studed the effect of small perturbatons n the corols centrfugal forces on the locaton of equlbrum ponts n the Robe s crcular restrcted three body problem when the denstes ρ ρ are equal evaluated the concomtant shft n the locaton of the equlbrum pont. In dervng the expresson for the buoyancy force Robe assumed that the pressure feld of the flud ρ has sphercal symmetry about the centre of the shell. He has taken nto account one of the three components of the pressure feld, that s, due to the own gravtatonal feld of the flud. Plastno & Plastno 995 nvestgated the effect of the remanng two components that orgnatng n the attracton of m, that arsng from the centrfugal force. They ncorporated these two components of pressure feld nto the dynamcs of the Robe model by consderng that the flud m adopts the shape of a Roche ellpsod. They have shown that the effect of the buoyancy force mght be thought as beng equvalent to a perturbaton of the corols force that the equlbrum pont s always stable,when the densty of the thrd body s greater than that of the surroundng medum. Gordano 997 studed the effect of drag force on the stablty of equlbrum ponts, both n the Robe s problem the problem studed by Plastno. Hallan & Rana 00 studed the effect of small perturbatons ɛ ɛ n the corols centrfugal forces respectvely on the locaton lnear stablty of the equlbrum ponts n the Robe s crcular problem wth denstes ρ ρ equal to zero vz, K equal to zero. They proved that there s only one equlbrum pont whch les on the lne jonng the centre of the shell the second prmary les to the rght or left of centre of the shell accordng as ɛ s postve or negatve. They have further studed the lnear stablty of the equlbrum pont. Hallan & Rana 00 later studed the effect of perturbatons n the corols centrfugal forces on the locaton ans stablty of the equlbrum ponts n the Robe s crcular restrcted problem wth densty parameter havng arbtrary value. Many authors have worked on problem of + bodes. Sgnfcant work s that of Whpple 984. He studed equlbrum solutons of the restrcted problem of + bodes. He further studed the lnear stablty of all the equlbrum solutons. Bhavneet & Rajv 0 extended the Robe s restrcted threebody problem to + bodes. Bhavneet & Rajv 0 studed the Robe s restrcted problem of + bodes when the bgger prmary s a Roche Ellpsod. They studed the equlbrum solutons of the nfntesmal masses analysed ther lnear stablty. In ths paper, we study the effect of perturbatons n the corols centrfugal forces on the locaton stablty of the equlbrum ponts n Robe s restrcted problem of + Bodes. Statement of the Problem Equatons of Moton In ths paper, one of the prmares of mass m s a Roche Ellpsod flled wth homogeneous ncompressble flud of densty ρ. The second prmary s a mass pont m m > m outsde the Ellpsod. The thrd the fourth body of mass m m 4 respectvely are small sold spheres of densty ρ ρ 4 respectvely nsde the Ellpsod, wth the assumpton that the mass radus of the thrd the fourth body are nfntesmal. Let R be the dstance between the centres of mass of m m. We assume that m descrbes a crcular orbt of radus R around m wth constant angular velocty ω. The masses m m 4 mutually attract each other, do not nfluence the motons of m m but are nfluenced by them. As n the case of classcal restrcted problem Szebehely 967,we adopt a unformly rotatng coordnate system Ox x x wth orgn of the coordnate system at the centre of the bgger prmary,ox pontng towards m Ox x beng the orbtal plane of m around m. The coordnate system Ox x x s as shown n the Fgure. Let the synodc system of coordnates ntally concdent wth the nertal system rotate wth angular velocty ω. Ths s the same as the angular velocty of m whch s descrbng a crcle around m. The equatons of moton of m m 4 n the Robe s crcular restrcted problem of + bodes n dmensonless cartesan coordnates Bhavneet & Rajv ISBN:

3 M x, x, x x O 0,0,0 M 4 x 4, x 4, x 4 R x Fgure : Geometry of the Robe s restrcted problem of + bodes when the bgger prmary m s a Roche Ellpsod the smaller prmary a pont mass 0are U = [ µj + x x + B + µ R j + µ { ρ [B + ρ { + µ = U, = U, = U } ] } + { m j µ j =, D = ρ m + m ρ B = πgρ I A A } ], 4, j =, 4; j. 5 A To acheve them we fx the unts such that m +m =, R =. We choose t n such a way that G =. The quantty µ becomes numercally equal to the rato m m + m. In the new unts ω =. The perturbatons n the corols centrfugal forces are expressed wth the help of parameters α, β, the unperturbed value of each beng unty. We take α β as α = + ɛ, ɛ, β = + ɛ, ɛ. Consequently, the equatons of moton of m m 4 can be wrtten as or α + α βx βx = U, 6 = U, 7 = U α = V + α, 8, 9 = V, 0 = V V =U + β x µ β + x = µ j [ + D B + { + R j { +µ ]. + + µ β + ɛ + ɛ = µ j } } The Equatons of moton of m m 4 can be rewrtten as xj = µ j R j + [ D + µ C + ɛ ], + + ɛ = µ j xj R j + [ D µ C + ɛ ], xj R j + D µ C 4, j =, 4; j. 5 C l = πgρ A l 6 µ = µ A l l =,,. 7 ISBN:

4 Equlbrum Solutons µ = µ πgρ. 8 The equlbrum solutons of m m 4 are gven by.e., µ j µ j µ j wth µ 4 V = V xj R j xj R j xj R j = V = 0 =, 4. + { D + µ C + ɛ } = 0, 9 + { D µ C + ɛ } = 0, 0 + D µ C = 0., j =, 4; j. Equlbrum Solutons lyng on x axs The Equatons 0, are satsfed wth, x ; =, 4 equal to zero. Hence we determne ; =, 4 from the Equaton 9 by puttng, x ; =, 4 equal to zero, we get x x 4 µ x x 4 + x 4 x x 4 x + { D + µ C + ɛ } x = 0. { D4 + µ C + ɛ } x 4 = 0. Multplyng the Equaton by µ by µ 4 addng, we get x 4 = λx 4 λ = D µ 5 D 4 µ 4 Substtutng n the Equaton,we get x = ± [ ] µ 4, {D + µ C + ɛ } + λ Mathematcal Applcatons n Modern Scence 6 [ x 4 = λ µ 4 {D + µ C + ɛ } + λ ]. 7 Hence x, 0, 0 x 4, 0, 0 are the equlbrum solutons for m m 4 respectvely provded they le wthn the Roche s Ellpsod. Equlbrum Solutons lyng on x axs If x = x = 0, then [ x = ± x 4 = D [ µ D 4 µ 4 µ 4 {D µ C + ɛ } + λ ], 8 µ 4 {D µ C + ɛ } + λ ]. 9 Hence 0, x, 0 0, x 4, 0 are the equlbrum solutons for m m 4, respectvely, provded they le wthn the Roche s Ellpsod. Equlbrum Solutons lyng on x axs If x = x = 0, then [ x = ± µ 4 {D µ C + ɛ } + λ ], 0 x 4 = D [ ] µ µ 4. D 4 µ 4 {D µ C + ɛ } + λ Hence 0, 0, x 0, 0, x 4 are the equlbrum solutons for m m 4 respectvely provded they le wthn the Roche s Ellpsod. There exst two equlbrum solutons of the system lyng each on x axs or x axs or x axs. Hence, there exst sx equlbrum solutons of the system, provded they le wthn Roche Ellpsod. We observe that there are no other equlbrum solutons except these sx. The equlbrum solutons of m are shown n Fgure.We observe that the equlbrum solutons donot exst when ρ = ρ ; ρ = ρ 4 or D = 0, =, 4 or when ρ = ρ ; ρ ρ 4. 4 Stablty of Equlbrum Solutons Let m m 4 be dsplaced to x j,, ; j =, 4 α j = X j +α j, = are very small X j, ISBN:

5 x x R Fgure : Locaton of Equlbrum Solutons n the Robe s Restrcted Problem of + Bodes when the Bgger Prmary s a Roche Ellpsod the Smaller Prmary a pont mass when the perturbatons n the corols centrfugal forces are taken nto consderaton. Crcles denote the postons of m trangles denote the postons of m 4. X j, Xj ; j =, 4 are the equlbrum soluton of m m 4 respectvely. The varatonal equatons of m m 4 are α α α + α α =α α =α V x + α V + α α =α + α V x x V V x x V x o + α x V x o o, o + α V x o o, o + α V o, =, 4 x o 4 V = µ j x Rj + µ V j xj = x µ j xj R 5 j µ V j xj = x V = µ j x Rj + R 5 j µ V j xj = x V = µ j x Rj + R 5 j xj xj µ j xj R 5 j R 5 j µ j xj, j =, 4; j. R 5 j + D + µ C + ɛ. xj + D µ C + ɛ. + D µ C Here the upper suffx o denotes the evaluaton of the dervatves at the equlbrum solutons under consderaton. Stablty of Equlbrum Solutons lyng on x axs Let the equlbrum soluton x 4, 0, 0 x, 0, 0 of m m 4 be ds- x + α, α, α placed to x 4 + α 4, α4, α4. The varatonal equatons of m m 4 are α α α = α α + α α = α p α = α p + λ + p λ = D µ D 4 µ 4, + λ + λ p + λ + p 5 6 =, 4 7 p = D + µ C + ɛ = p + ɛ p = D µ C + ɛ = p + ɛ p = D µ C. ISBN:

6 From the Equaton 7, we ascertan that the moton of m m 4 parallel to x axs s stable when D + µ C + ɛ D µ C D + µ C + ɛ D µ C > + D µ D 4 µ 4 =, 4 n case ρ > ρ, ρ > ρ 4 8 < + D µ D 4 µ 4 =, 4 n case ρ < ρ, ρ < ρ 4 9 The remanng Equatons 5 6 admt solutons of the form α j = A j e Ljt, =,, ; j =, 4. The characterstc equatons of m m 4 are gven by L 4 + L 4 + ɛ Q + Q + Q Q = 0 Q = p Q = p + λ + λ p + λ. =, 4 40 Let Λ = L, we obtan Λ + Λ 4 + ɛ Q + Q + Q Q = 0 =, 4. 4 we must have S + S < 0, S S > 0, > 0. Case I:ρ < ρ, ρ < ρ 4. In ths case Q < 0. Now, S S > 0 f Q < 0.e. f, D + µ C + ɛ D µ C + ɛ < + D µ D 4 µ 4 =, 4 4 Ths mples that S + S < 0, > 0. Hence, the equlbrum solutons of m m 4 are stable f the Equaton 4 holds, else unstable. Case II:ρ > ρ, ρ > ρ 4. Thus Q > 0. Now, S S > 0 f Q > 0. The equlbrum solutons of m m 4 are stable f D µ C D + µ C λ + ɛ > 0, + λ + λ 4 + λ + λ D + µ C + D µ C + ɛ + λ + λ < 4 + ɛ, 6 + 4ɛ + p 4 + λ + λ 6ɛ p 4 + λ p + λ + λ 8 + ɛ p + λ + λ p + p + λ 8ɛ > 0 + λ The shaded regon n Fgure shows the regon n whch the equlbrum solutons of m are stable for λ = 0. ɛ 0 8, ɛ Let the roots of the Equaton 4 be S, S. Then, S + S = Q + Q 4 + ɛ, S S = Q Q. Dscrmnant of the Equaton 4 s = 6 + 4ɛ + Q Q 8 + ɛ Q + Q p The equlbrum solutons x, 0, 0 x 4, 0, 0 of m m 4 respectvely when the dsplacement s gven n the drecton of x axs or x axs are stable f S S are real negatve or p Fgure : The shaded regon corresponds to the regon of stablty for the equlbrum solutons of m lyng on x axs ISBN:

7 Stablty of Equlbrum Solutons lyng on x axs Let the equlbrum solutons 0, x, 0 0, x 4, 0 of m m 4 be dsplaced to α, x + α, α α 4, x4 + α 4, α4. The varatonal equatons of m m 4 are α α α = α α + α α = α α = α p p + λ + p + λ + λ p + λ + p =, The moton of m m 4 parallel to x axs s stable when p p > + λ =, 4 n case ρ > ρ, ρ > ρ 4 47 p p < + λ =, 4 n case ρ < ρ, ρ < ρ 4 48 The characterstc equatons of m m 4 are gven by L 4 + L 4 + ɛ Q + Q + Q Q = 0 Q = p Q = p + λ + λ p + λ. =, 4 49 Case I:ρ < ρ, ρ < ρ 4. In ths case D > 0, λ > 0, p < 0, p < 0 =, 4. Thus Q < 0. The equlbrum solutons x 4, 0, 0 x, 0, 0 of m m 4 respectvely when the dsplacement s gven n the drecton of x axs or x axs are stable f Q < 0,.e. f D µ C + ɛ D + µ C + ɛ < + λ =, 4 50 as then S + S < 0, > 0. Case II:ρ > ρ, ρ > ρ 4. In ths case D < 0, λ > 0, p > 0, p > 0 =, 4. Thus Q > 0. Now, S S > 0, f Q > 0. The regon that s common to Q > 0, Q + Q < 4 + ɛ > 0 s the stable regon for the equlbrum solutons x, 0, 0 x 4, 0, 0 of m m 4 respectvely when the dsplacement s gven n the drecton of x axs or x axs. The shaded regon n Fgure 4 shows the regon n whch equlbrum solutons of m are stable for λ = 0. ɛ 0 8, ɛ 0 9 Q + Q = p =6 + 4ɛ + + 6ɛ + λ 8 + ɛ p + λ + p + λ + ɛ + λ + λ 4 + λ + λ 4 + λ p p + λ + λ + λ p p p + λ 8ɛ + λ 5 Stablty of Equlbrum Solutons lyng on x axs Let the equlbrum soluton 0, 0, x 0, 0, x 4 of m m 4 be dsplaced to α, α, x + α α 4, α4, x4 + α 4. The varatonal equatons of m m 4 are α α α α + α = α α = α α = α p p + λ + p p + λ + p + λ + λ 5 5 =, 4 54 From the Equaton 54, we ascertan that the moton of m m 4 parallel to x axs s stable f p < 0 ISBN:

8 Mathematcal Applcatons n Modern Scence when the dsplacement s gven n the drecton of x axs or x axs are stable. If Q > 0 Q > 0, that s common the regon 40 0 p X\ to Q > 0, Q > 0, Q + Q < 4 + > 0 s the stable regon for the equlbrum solutons of m m4 when the dsplacement s gven n the drecton of x axs or x axs. Case II: Smlar analyss follows when ρ > ρ, ρ > ρ4. The shaded regon n Fgure 5 shows the regon n whch equlbrum solutons of m are stable for λ = , λ Q +Q = p +p p + 0 +λ +λ 40 p X\ Fgure 4: The shaded regon corresponds to the regon of stablty for the equlbrum solutons of m lyng on x axs. =6 + 4 p p p + p + or equvalently D > 0 =, 4.e, ρ < ρ ρ < ρ4. The characterstc equatons of m m4 from the Equaton 5 5 are gven by L4 + L 4 + Q + Q + Q Q = 0 0 p +λ!! p X\ =, 4 55 p X\ Q = p 0 p +λ Q = p p. +λ p X\ or equvalently Λ + Λ 4 + Q + Q + Q Q = 0 5 Fgure 5: The shaded regon corresponds to the regon of stablty for the equlbrum solutons of m lyng on x axs. =, wth Λ = L Case I:ρ < ρ, ρ < ρ4. In ths case D > 0, λ > 0, p < 0, p < 0, p < 0 =, 4. As p < 0, the moton of m m4 parallel to x axs s always stable. Now, S S > 0 f Q > 0 Q > 0 or f Q < 0 Q < 0. If Q < 0 Q < 0, then S +S < 0, > 0. In ths case, the equlbrum solutons of m m4 ISBN: Applcatons Applcaton Moton of ceberg n the Earth-Moon system. We take ρ = 07 kg/m, ρ = 804 kg/m, ρ4 = 80 kg/m We have assumed a = Equatoral radus of earth = km. a = km. 86 λ +λ

9 a = km. Roche Lmt 0 Equatoral radus of moon = 78. km. Polar radus of moon = 76.0 km. Unts of Mass: In dmensonless unts, m + m =. Mass of Earth m = kg. Mass of Moon m = kg. Thus, kg = new unt Therefore, 4 µ = m m + m = 0.0. m µ = = 0 8. m + m m 4 µ 4 = = 0 7. m + m Unts of Dstance: Dstance between Earth Moon = km = new unt Thus, metre = new unt. Hence, ρ = 9645 new unts. ρ = 7550 new unts. ρ 4 = 7607 new unts. D = ρ = 0.77 ρ D 4 = ρ = 0.67 ρ 4 a = a = a = C = C = C = We take ɛ 0 8, ɛ 0 9. Table gves the equlbrum solutons for m m 4 usng the above data. The numercal solutons are correct n all dgts prnted. Fgure 6 plots the equlbrum solutons of m lyng on x axs vs µ. The graph s plotted for the postve value of x. For the earth moon system, we have µ = 0.0. The hghlghted pont represents the equlbrum soluton of m lyng on x axs for ths value of µ. Fgure 7 plots the correspondng equlbrum soluton of m 4 lyng on x axs vs µ when the equlbrum soluton of m m 4 are related by the Equaton 4. Stablty of Equlbrum Solutons lyng on x axs x axs The moton of m m 4 parallel to x axs s stable. The equlbrum solutons of m m 4 respectvely when the dsplacement s gven n the drecton of x axs or x axs are unstable as they donot satsfy all the condtons of stablty. Stablty of Equlbrum Solutons lyng on x axs The moton of m m 4 parallel to x axs s unstable. The equlbrum solutons of m m 4 respectvely when the dsplacement s gven n the drecton of x axs or x axs are unstable. 5. Applcaton Moton of submarnes n the Earth-Moon System We take ρ = 07 kg/m, ρ = 00 kg/m, ρ 4 = 00 kg/m In new unts ρ = 9645 new unts. ρ = 00 new unts. ρ 4 = 69 new unts. D = ρ = ρ D 4 = ρ = 0.44 ρ 4 a = a = a = C = C = C = We observe that n ths case the equlbrum solutons of m m 4 donot exst. 6 Dscusson In ths paper, we have studed the effect of perturbatons n the corols centrfugal forces on the locaton stablty of the equlbrum ponts n Robe s restrcted problem of + bodes. Plastno & Plastno 995 revsed the Robe s restrcted three-body problem under the assumpton that the flud body assumes ISBN:

10 the shape of a Roche Ellpsod. They took nto consderaton all the three components of the pressure feld n dervng the expresson for the buoyancy force vz, due to the own gravtatonal feld of the flud, that orgnatng n the attracton of m that arsng from the centrfugal force. Bhavneet & Rajv 0 studed the equlbrum solutons lnear stablty of the mnor bodes m m 4 n the Robe s restrcted problem of + bodes when the bgger prmary s a Roche ellpsod.? studed the equlbrum solutons lnear stablty of the mnor bodes m m 4 n the Robe s restrcted problem of + bodes when the bgger prmary s a Roche ellpsod smaller prmary an oblate body. In the paper Bhavneet & Rajv 0, we had neglected the tdal deformaton of the flud body m due to m. Also the components of buoyancy force arsng from the attracton of m centrfugal forces were gnored. We had studed the moton of the system n two dmenson as here by extendng to three dmenson we have determned the equatons of moton of m m 4 n the x drecton as well. All these have been ncorporated n the present study. The equatons of moton, equlbrum solutons of m m 4 ther stablty are analyzed takng nto account the above consderatons, perturbatons n the corols centrfugal forces takng m a Roche Ellpsod. Our results are n contrast to those of Plastno & Plastno 995. In ther study, when ρ = ρ, every pont wthn the Roche Ellpsod s an equlbrum soluton n general 0, 0, 0 s the only equlbrum soluton. We have taken two nfntesmal masses wthn the Roche Ellpsod as Plastno & Plastno 995 has taken only one nfntesmal mass wthn the Roche Ellpsod. When there are small perturbatons ɛ ɛ n the corols centrfugal forces respectvely, n the Robe s restrcted problem of + bodes, the number of equlbrum solutons are the same as n the problem wth no perturbaton Bhavneet & Rajv 0, but postons of the equlbrum solutons have changed. The change n the corols force does not affect the postons of the equlbrum solutons. There exst two equlbrum solutons of the system lyng each on x axs or x axs or x axs. Hence, there exst sx equlbrum solutons of the system, provded they le wthn Roche Ellpsod. We observe that there are no other equlbrum solutons except these sx It may also be observed that the problem of + bodes can be extended to + n bodes takng n nfntesmal masses wthn the ellpsod. Ths problem has applcaton n the moton of submarnes n the earth moon setup. The problem can also be appled to the moton of ceberg. Earth s the only celestal body we know of whch has flud the test partcles wthn the flud are consdered as the submarnes. In the solar or the extra solar system, we donot know of any other celestal body havng lqud nsde t. Our entre work s applcable n solar or extra solar system when f such a system s dscovered. 7 Concluson In ths paper, we have studed the moton of two nfntesmal masses m m 4 supposed movng nsde m n three dmensons, takng m as a pont mass m a Roche Ellpsod. We have studed the effect of perturbatons n the corols centrfugal forces on the locaton stablty of the equlbrum ponts n Robe s restrcted problem of + bodes. We have taken nto consderaton all the three components of the pressure feld n dervng the expresson for the buoyancy force vz, due to the own gravtatonal feld of the flud, that orgnatng n the attracton of m that arsng from the centrfugal force.we have proved that there exst only sx equlbrum solutons of the system, provded they le wthn the Roche Ellpsod. In a system the prmares are consdered as earth-moon m, m 4 as ceberg, the equlbrum solutons of m m 4 respectvely when the dsplacement s gven n the drecton of x axs or x axs are unstable. Table : Equlbrum Solutons x, x 4 =,, of m m 4 lyng on x, x, x axs respectvely. x x 4 x axs ± x axs ± x axs ± ISBN:

11 x Μ Fgure 6: Equlbrum soluton of m postve root lyng on x axs vs µ. The hghlghted pont represents the equlbrum soluton of m for µ = 0.0 x Μ Fgure 7: Equlbrum soluton of m 4 lyng on x axs vs µ. The hghlghted pont represents the equlbrum soluton of m 4 for µ = 0.0 ISBN:

12 References: Bhatnagar, K. B. Hallan, P.P.:978, Effect of Perturbatons n Corols Centrfugal Forces on the Stablty of Lbraton Ponts n the Restrcted Problem, Celestal Mechancs 8,05-. Bhavneet Kaur, Rajv Aggarwal, Robe s restrcted problem of + bodes when the bgger prmary s a Roche ellpsod, Acta Astronautca, Volume 89, AugustSeptember 0, Pages -7, ISSN , Szebehely, V.S.:967 Theory of orbts, Academc Press, New York. Szebehely V, Peters CF 967 Complete soluton of a general problem of three bodes, Astronomcal Journal 7:876-88, DOI 0.086/055. Wkpeda:0, Roche Lmt, Wkpeda. Whpple, A.L.:984, Equlbrum Solutons of the Restrcted Problem of + Bodes, Celestal Mechancs,7-94. Chrashekhar, S.:987 Ellpsodal Fgures of Equlbrum Chapter 8, Dover Publcaton Inc, New York. Gordano, C.M., Plastno, A.R. Plastno, A.:997, Robe s Restrcted Three-Body Problem wth Drag, Celestal Mechancs Dynamcal Astronomy 6, Hallan, P.P. Rana, N.:00, The Exstence Stablty of Equlbrum Ponts n the Robe s Restrcted Three-Body Problem, Celestal Mechancs Dynamcal Astronomy 79, Hallan, P.P. Rana, N.:00, Effect of Perturbatons n Corols Centrfugal Forces on the Locaton Stablty of Equlbrum Pont n the Robe s Crcular Restrcted Three-Body Problem, Planetary Space Scence 99, Hallan, P.P. Rana, N.:00, Effect of Perturbatons n the Corols Centrfugal Forces on the Locatons Stablty of the Equlbrum Pont n the Robe s Crcular Problem wth Densty Parameter Havng Arbtrary Value, Indan J. pure appl. Math., 47, Kaur, Bhavneet Aggarwal, R.:0 Robe s Problem: Its Extenson to + Bodes, Astrophyscs Space Scence 9, Plastno, A.R. & Plastno, A. :995, Robe s Restrcted Three-Body Problem Revsted, Celestal Mechancs Dynamcal Astronomy 6, Robe, H.A.G.:977, A New Knd of Three-Body Problem, Celestal Mechancs 6, 4-5. Shrvastava, A. K. Garan, D.: 99, Effect of perturbaton on the locaton of lbraton pont n the Robe restrcted problem of three bodes, Celest. Mech. & Dyn. Astr. 5, 677. ISBN:

Elshaboury SM et al.; Sch. J. Phys. Math. Stat., 2015; Vol-2; Issue-2B (Mar-May); pp

Elshaboury SM et al.; Sch. J. Phys. Math. Stat., 2015; Vol-2; Issue-2B (Mar-May); pp Elshabour SM et al.; Sch. J. Phs. Math. Stat. 5; Vol-; Issue-B (Mar-Ma); pp-69-75 Scholars Journal of Phscs Mathematcs Statstcs Sch. J. Phs. Math. Stat. 5; (B):69-75 Scholars Academc Scentfc Publshers

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

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

The Photogravitational Circular Restricted Four-body Problem with Variable Masses

The Photogravitational Circular Restricted Four-body Problem with Variable Masses The Photogravtatonal Crcular Restrcted Four-body Problem wth Varable Masses Abdullah A. Ansar Department of Mathematcs College of Scence Majmaah Unversty Al-Zulf Kngdom of Saud Araba Abstract a.ansar@mu.edu.sa

More information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

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

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

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

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

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

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

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

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

ˆ (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

χ 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

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding. Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays

More information

So far: simple (planar) geometries

So far: simple (planar) geometries Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector

More information

The classical spin-rotation coupling

The classical spin-rotation coupling LOUAI H. ELZEIN 2018 All Rghts Reserved The classcal spn-rotaton couplng Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 louaelzen@gmal.com Abstract Ths paper s prepared to show that a rgd

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

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

One Dimensional Axial Deformations

One Dimensional Axial Deformations One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the

More information

Errors in Nobel Prize for Physics (7) Improper Schrodinger Equation and Dirac Equation

Errors in Nobel Prize for Physics (7) Improper Schrodinger Equation and Dirac Equation Errors n Nobel Prze for Physcs (7) Improper Schrodnger Equaton and Drac Equaton u Yuhua (CNOOC Research Insttute, E-mal:fuyh945@sna.com) Abstract: One of the reasons for 933 Nobel Prze for physcs s for

More information

Solutions to Exercises in Astrophysical Gas Dynamics

Solutions to Exercises in Astrophysical Gas Dynamics 1 Solutons to Exercses n Astrophyscal Gas Dynamcs 1. (a). Snce u 1, v are vectors then, under an orthogonal transformaton, u = a j u j v = a k u k Therefore, u v = a j a k u j v k = δ jk u j v k = u j

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

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2 Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on

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

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

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

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

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

Study Guide For Exam Two

Study Guide For Exam Two Study Gude For Exam Two Physcs 2210 Albretsen Updated: 08/02/2018 All Other Prevous Study Gudes Modules 01-06 Module 07 Work Work done by a constant force F over a dstance s : Work done by varyng force

More information

Chapter 11 Angular Momentum

Chapter 11 Angular Momentum Chapter 11 Angular Momentum Analyss Model: Nonsolated System (Angular Momentum) Angular Momentum of a Rotatng Rgd Object Analyss Model: Isolated System (Angular Momentum) Angular Momentum of a Partcle

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

CHAPTER 10 ROTATIONAL MOTION

CHAPTER 10 ROTATIONAL MOTION CHAPTER 0 ROTATONAL MOTON 0. ANGULAR VELOCTY Consder argd body rotates about a fxed axs through pont O n x-y plane as shown. Any partcle at pont P n ths rgd body rotates n a crcle of radus r about O. The

More information

Chapter 20 Rigid Body: Translation and Rotational Motion Kinematics for Fixed Axis Rotation

Chapter 20 Rigid Body: Translation and Rotational Motion Kinematics for Fixed Axis Rotation Chapter 20 Rgd Body: Translaton and Rotatonal Moton Knematcs for Fxed Axs Rotaton 20.1 Introducton... 1 20.2 Constraned Moton: Translaton and Rotaton... 1 20.2.1 Rollng wthout slppng... 5 Example 20.1

More information

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017 17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran

More information

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics) CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

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

Physics 111: Mechanics Lecture 11

Physics 111: Mechanics Lecture 11 Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

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

Chapter 11: Angular Momentum

Chapter 11: Angular Momentum Chapter 11: ngular Momentum Statc Equlbrum In Chap. 4 we studed the equlbrum of pontobjects (mass m) wth the applcaton of Newton s aws F 0 F x y, 0 Therefore, no lnear (translatonal) acceleraton, a0 For

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

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

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

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

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11) Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the

More information

THE GRAVITATIONAL TWO-BODY PROBLEM

THE GRAVITATIONAL TWO-BODY PROBLEM THE GRAVITATIONAL TWO-BODY PROBLEM Danel J. Scheeres Department of Aerospace Engneerng Scences, The Unversty of Colorado, USA Keywords: Two-body problem, Kepleran soluton, dstrbuted bodes, rotatonal moton,

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 Smooth Length for SPH Modelling of High Speed Impact

Tensor Smooth Length for SPH Modelling of High Speed Impact Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru

More information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

Conservation of Angular Momentum = "Spin"

Conservation of Angular Momentum = Spin Page 1 of 6 Conservaton of Angular Momentum = "Spn" We can assgn a drecton to the angular velocty: drecton of = drecton of axs + rght hand rule (wth rght hand, curl fngers n drecton of rotaton, thumb ponts

More information

V.C The Niemeijer van Leeuwen Cumulant Approximation

V.C The Niemeijer van Leeuwen Cumulant Approximation V.C The Nemejer van Leeuwen Cumulant Approxmaton Unfortunately, the decmaton procedure cannot be performed exactly n hgher dmensons. For example, the square lattce can be dvded nto two sublattces. For

More information

1. Review of Mechanics Newton s Laws

1. Review of Mechanics Newton s Laws . Revew of Mechancs.. Newton s Laws Moton of partcles. Let the poston of the partcle be gven by r. We can always express ths n Cartesan coordnates: r = xˆx + yŷ + zẑ, () where we wll always use ˆ (crcumflex)

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Spring 2002 Lecture #13

Spring 2002 Lecture #13 44-50 Sprng 00 ecture # Dr. Jaehoon Yu. Rotatonal Energy. Computaton of oments of nerta. Parallel-as Theorem 4. Torque & Angular Acceleraton 5. Work, Power, & Energy of Rotatonal otons Remember the md-term

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

How Differential Equations Arise. Newton s Second Law of Motion

How Differential Equations Arise. Newton s Second Law of Motion page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons

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

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

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018 MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng

More information

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of

More information

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.

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

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors. SCALARS AND ECTORS All phscal uanttes n engneerng mechancs are measured usng ether scalars or vectors. Scalar. A scalar s an postve or negatve phscal uantt that can be completel specfed b ts magntude.

More information

ENGI 1313 Mechanics I

ENGI 1313 Mechanics I ENGI 11 Mechancs I Lecture 40: Center of Gravty, Center of Mass and Geometrc Centrod Shan Kenny, Ph.D., P.Eng. ssstant Professor Faculty of Engneerng and ppled Scence Memoral Unversty of Nefoundland spkenny@engr.mun.ca

More information

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis Part C Dynamcs and Statcs of Rgd Body Chapter 5 Rotaton of a Rgd Body About a Fxed Axs 5.. Rotatonal Varables 5.. Rotaton wth Constant Angular Acceleraton 5.3. Knetc Energy of Rotaton, Rotatonal Inerta

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

10/24/2013. PHY 113 C General Physics I 11 AM 12:15 PM TR Olin 101. Plan for Lecture 17: Review of Chapters 9-13, 15-16

10/24/2013. PHY 113 C General Physics I 11 AM 12:15 PM TR Olin 101. Plan for Lecture 17: Review of Chapters 9-13, 15-16 0/4/03 PHY 3 C General Physcs I AM :5 PM T Oln 0 Plan or Lecture 7: evew o Chapters 9-3, 5-6. Comment on exam and advce or preparaton. evew 3. Example problems 0/4/03 PHY 3 C Fall 03 -- Lecture 7 0/4/03

More information

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

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

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS 6th ICPT, Sapporo, Japan, July 008 APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS James MAINA Prncpal Researcher, Transport and Infrastructure Engneerng, CSIR Bult Envronment

More information

Infinitesimal Orbits Near Any Point on Fawzy Equilibrium Circle in the Three Dimensional Restricted Three Body Problem

Infinitesimal Orbits Near Any Point on Fawzy Equilibrium Circle in the Three Dimensional Restricted Three Body Problem Earth Scence Inda, eissn: 097 850 Vol. 6 (I, January, 0, pp. - http://www.earthscencenda.nfo/ Infntesmal Orbts Near ny Pont on Fawzy Equlbrum Crcle n the hree Dmensonal Restrcted hree Body Problem F..

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

coordinates. Then, the position vectors are described by

coordinates. Then, the position vectors are described by Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,

More information

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force. Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act

More information

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,

More information

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

More information

Work is the change in energy of a system (neglecting heat transfer). To examine what could

Work is the change in energy of a system (neglecting heat transfer). To examine what could Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes

More information

Significance of Dirichlet Series Solution for a Boundary Value Problem

Significance of Dirichlet Series Solution for a Boundary Value Problem IOSR Journal of Engneerng (IOSRJEN) ISSN (e): 5-3 ISSN (p): 78-879 Vol. 6 Issue 6(June. 6) V PP 8-6 www.osrjen.org Sgnfcance of Drchlet Seres Soluton for a Boundary Value Problem Achala L. Nargund* and

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

PY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg

PY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg PY2101 Classcal Mechancs Dr. Síle Nc Chormac, Room 215 D Kane Bldg s.ncchormac@ucc.e Lectures stll some ssues to resolve. Slots shared between PY2101 and PY2104. Hope to have t fnalsed by tomorrow. Mondays

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

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

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

CHAPTER 4. Vector Spaces

CHAPTER 4. Vector Spaces man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear

More information

Calculating the Quasi-static Pressures of Confined Explosions Considering Chemical Reactions under the Constant Entropy Assumption

Calculating the Quasi-static Pressures of Confined Explosions Considering Chemical Reactions under the Constant Entropy Assumption Appled Mechancs and Materals Onlne: 202-04-20 ISS: 662-7482, ol. 64, pp 396-400 do:0.4028/www.scentfc.net/amm.64.396 202 Trans Tech Publcatons, Swtzerland Calculatng the Quas-statc Pressures of Confned

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

11. Introduction to Liouville s Theorem Michael Fowler

11. Introduction to Liouville s Theorem Michael Fowler 11 Introducton to Louvlle s Theorem Mchael Fowler Paths n Smple Phase Spaces: the SHO and Fallng Bodes Let s frst thnk further about paths n phase space For example the smple harmonc oscllator wth 2 2

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information