Research Article Matched Asymptotic Expansions to the Circular Sitnikov Problem with Long Period

Transcription

1 Journal of Applied Mathematics Volume 0, Article ID 79093, 5 pages doi:0.55/0/79093 esearch Article Matched Asymptotic Expansions to the Circular Sitnikov Problem with Long Period Marco osales-vera, High Temperatures Department, Instituto de Innovación en Minería y Metalurgia, Ahumada 3, 8308 Santiago, Chile Department of Physics, Universidad Gabriela Mistral, Avenida icardo Lyon 77, Santiago, Chile Correspondence should be addressed to Marco osales-vera, mrosa00@im.codelco.cl eceived 3 August 0; evised November 0; Accepted 3 November 0 Academic Editor: oberto Barrio Copyright q 0 Marco osales-vera. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The circular Sitnikov problem is revisited, using matched asymptotic expansions. In the case of large oscillation periods, approximate analytical expressions for the period and the orbit of the third body are found. The results are compared with those described in the literature and show that the movement of the third body can be well described by two analytical solutions, the inner and outer solutions.. Introduction The Sitnikov problem is a particular case of the restricted three-body problem. The configuration of the Sitnikov problem is defined by: two point-like bodies of equal masses called primaries orbiting around their common centre of mass due to their mutual gravitational forces, and a third body of negligible mass moving along a line, perpendicular to the orbital plane of the primaries, going through their barycentre and performs oscillations along of the straight line Figure. The circular Sitnikov problem is obtained when the two primary bodies describe circular orbits. This situation was originally discussed by MacMillan, demonstrating that this problem can be reduced to elliptic integrals. Later, several authors have analyzed the general and circular Sitnikov problem. The case when the two primaries bodies describe elliptical orbits was first investigated by Sitnikov, who proved the existence of oscillatory solutions for this problem. Subsequently, Moser 3 proved the existence the chaotic orbits. Perdios and Markellos examined the stability and bifurcations of straight line motions of the third body. Dvorak 5 studied, by numerical means, motion of the planetoid limited

2 Journal of Applied Mathematics m M M Figure : Sitnikov problem geometry. to a small region around the barycentre of the primaries and found that invariant curves exist for very small oscillations centering the barycentre. For bounded orbits with moderate amplitudes, Hagel 6 showed that the equations of motion can be expressed in a polynomial form. Liu and Sun 7 derived a mapping model to investigate the problem. Wodnar 8 introduced a new formulation for the equation of motion by using the true anomaly of the primaries as independent variable. Belbruno et al. 9 derived analytical expressions for the circular case MacMillan problem using elliptic functions. Faruque 0 found approximate analytical solutions for small oscillations around the barycenter and moderate eccentricities. ecently, Dvorak studied the complete phase space numerically. Kovács and Érdi reported the extended phase space of the Sitnikov problem by using a stroboscopic map and computing escape times. The instance when the mass of the third body is nonnegligible extended Sitnikov problem has been studied by Dvorak and Sun 3. Soulis et al. also studied the situation in which the third body could move off the z-axis generalized extended Sitnikov problem. In 009, Hagel 5 derived an analytical expression for the perihelion motion of the primaries, when the third body has a finite mass. Bountis and Papadakis 6 analyzed the problem for an arbitrary number of primary masses, Sidorenko 7 complemented the previous work on the alternation of stability and instability in the family of vertical motions. ecently, uzza and Lhotka 8 performed the implementation of a numerical construction of high order normal form near elliptic orbits. In this work, we assume that the energy of the third body is close to the escape energy. In this situation, the period of the orbit is large, because the body goes to great distances from the origin and then return. The numerical calculation of the period becomes complex, mainly because computational times are very long and it becomes necessary to resort to asymptotic analysis to support the numerical calculation.. Equations of Motion The equation of motion of the third body is z z GM r t z, 3/.

3 Journal of Applied Mathematics 3 where rt is the distance between the primary bodies and zt the distance of the third body from the center of mass and t representing time. The distance between the primary bodies is given by rt e cos Et,. where Et is the eccentric anomaly of the primaries, e is the eccentricity of the orbits of the primaries. When the primaries move in circular orbits, e 0,. is reduced to z z GM z. 3/.3 Equation.3 represents the circular Sitnikov problem, also known as the problem of MacMillan. Integrating.3 with the initial conditions: żt 0 v 0 and zt 0 0we obtain ż z/ ) / v 0.. Defining the parameter β as v 0 β. GM/.5 Equation. is transformed into ż β z/ )..6 / Making the change of variables z/ /u, and then η u/β, allows us to find the following expression to describe the position of the third body : GM t 3 ηη η0 k η ) η ) k η )dη; k β..7 The period is given by GM 3 T 0 k η ) η ) k η )dη..8

4 Journal of Applied Mathematics This integral can be performed using elliptic integrals, but it is impossible to reverse.7 to find an explicit expression for zt. However, it is possible to find an analytical approximation for the period of the orbit and zt by matching asymptotic expansions when the energy of the third body is close to the energy of escape. The concept of escape velocity used in this paper is the minimum velocity at the origin to obtain an unbounded orbit. From., the following value for the escape velocity is found: v 0esc GM/. In this case, we have β. 3. Asymptotic Expansions When the third body has a speed close to the escape velocity, that is, β ε,.6 can be written as ż z/ ) ε. 3. / The mechanical total energy E of the third body, can be written in terms of ε as: E GMm/ε, bounded orbit are obtained when 0 <ε<. In the case ε 0, an unbounded orbit is obtained, E 0 is the minimum energy of the third body to escape of the gravitational field of the primary bodies. If ε is a positive small parameter, the orbit of the third body is bounded but with a large amplitude, then the equation can be solved by matching between two asymptotic solutions. 3.. Near Field In the regions close to the origin, we have: / z/ / ε, thus, the solution of 3. can be approximated to the following inner asymptotic expansion: z in t z 0 t εz t. 3. The validity of the expansion 3. will be verified in the matching process see Section 3.3. Introducing the asymptotic expansion 3. in 3., we obtain the following equation for the leading order and the first order correction: ż 0 z 0 / ), 3.3 / ż GM ż 0 z z 0 z 0 / ). 3. 3/ The initial conditions are: z 0 t 0 0andz t 0 0.

5 Journal of Applied Mathematics 5 Equation 3.3 represents the movement of a particle passing through the origin with a speed equal to the escape velocity. Integrating 3.3 gives z0 ) ) / dz 0 t c. 3.5 This integral can be calculated analytically z 0 z0 ) ) / z F, 3 ; 3 ; z0 ) ) t c 3, 3.6 where F a, b; c; z is the hypergeometric function. From the initial condition z 0 t 0 0, we have c 0, 3. can be transformed to more useful as follows: dz dz 0 ż 0 GM z z 0 z 0 / ). 3.7 / Introducing 3.3 into 3.7, we obtain the following equation for first order perturbations: dz dz 0 z z 0 z 0 / ) z0 ) ) /. 3.8 The solution of 3.8 is given by the expression z z 0 z 0 / ) / z0 ) ) 3/ dz This integral can be expressed in terms of the hypergeometric function as z z 0 z 0 / ) / 5 z 0 z0 ) ) 3/ 3 5 z 0 F, ; 3 ; z0 ) )) Far Field In the region far from the origin, the parameter ε has an important role in 3., because ε defines the maximum distance reached by the third body from the origin: z max /ε ε /ε. Moreover, in this region we have z/ and3. can be approximated by ż ) z ε. 3.

6 6 Journal of Applied Mathematics This differential equation represents the outer solution. Integrating this equation, the following expression is obtained: z dz εz t c. 3. This integral can also be calculated analytically, that is: εz )) εz εz tan εz ε 3/ εz εz t c, 3.3 where the value of the constant c must be calculated by matching between the inner and outer solutions Matching Procedure In the region z/ /ε, the inner and outer solutions are valid, then 3.6 and 3.3 must make matching in it. In the region z/, 3.6 has the following asymptotic representation: z0 ) 3/ 3 z 0 ) ) Γ/Γ3/ t, 3. 3 Γ3/Γ 3 where we have used the following asymptotic expansion for the hypergeometric function: F a, b; c; z Γb aγc ΓbΓc a z a Γa bγc ΓaΓc b z b ) a a c a b z b b c a bz ). 3.5 In the region z/ /ε, the asymptotic solution 3.3 can be expanded in power series of εz/, then we obtain the equation of first order in εz/, namely z ) 3/ Oε t 3 c In the leading order, the matching between 3. and 3.7 delivers the following value for the constant c : ) c π Γ. 3.7

7 Journal of Applied Mathematics 7 Introducing 3.8 into 3.3, the following expression is obtained to zt in the far field: ε 3/ εz )) εz εz tan εz t π Γ ). 3.8 We may proceed further by combining the inner solution 3.6 and the outer solution 3.8 into a single uniform approximation z unif t: z unif t z in t z out t z match t, 3.9 where z match t is obtained from 3., given: zt match 3 ) /3 t 3 π Γ ) / The quarter-period of the orbit is obtained by evaluating 3.8 in z /ε, thus we get the following expression for the period: 3 T π ε3/ 3 π Γ ) Matching in Higher Order To carry out the asymptotic matching for a first order in ε is more convenient to write 3. in the inner region as follows: z/ ) / dz ε z/ t c. 3. Expanding 3. in powers of ε,weobtain: n0 z ) ) /n a n ε n dz This integral can also be calculated analytically, that is: a n ε n z F, n ; 3 z ) ) ; n0 From the initial condition z 0 t 0 0, we have c 0. t c. 3.3 t c. 3.

8 8 Journal of Applied Mathematics In the region z/, 3. has this asymptotic representation ) 3 n Γ /n 3 Γ a n ε Γ /n n0 Γ3/ Γ/ n0 n Γ/n 3 a n ε Γ/n 7 z ) /n3 t. 3.5 To perform a more accurate estimation of the period, a more accurate estimate of the behavior of z for large z is necessary, a better approximation of 3. for large values of z leads to ż ε. 3.6 z/ /z Separating variables in 3.5, and after algebraic manipulation, leads to the following approximate integral: z z εz / dz t c. 3.7 Finally, the following expression is found for improved outer solution: ε 3/ ε z ) ε ) ε z /ε ε εz tan /ε εz) t c. 3.8 In the region z/ /ε, the asymptotic solution 3.8 can be expanded in power series of εz/, then we obtain the asymptotic expansion: z 3/ z 5/ 3 5 ) 3 ε t c. 3.9 The matching between 3.5 and 3.9 delivers the following value for the constant c: ) 3 Γ 3/ c Γ Γ / Γ 3 ) Γ 5/ Γ 3/ ε O ε ). 3.30

9 Journal of Applied Mathematics 9 zt) ε = GM 3 t Figure : The function zt with ε 0., the solid line corresponds to numerical integration of 3., the dashed line represents the outer solution 3.3. Introducing 3.30 into 3.8, the following expression is obtained to zt in the far field ε 3/ ε z ) ε ) εz /ε ε εz tan /ε εz) 3 t Γ ) Γ 3/ Γ / Γ 3 ) Γ 5/ Γ 3/ ε. 3.3 The quarter-period of the orbit is obtained by evaluating 3.3 in z /ε /ε,thus obtaining the following expression for the period including corrections to order ε : 3 T π ε3/ 3 π Γ ) π 3 5 Γ / ε 5 5 π Γ ) ε esults In this section, the results are expressed using dimensionless variables, the dimensionless time τ GM/ 3 t, and the dimensionless distance η z/. In this representation, the results are general and it is not necessary to give numerical values for the constants G, M and. Figure presents the function zt with ε 0.. The solid line corresponds to the solution obtained by numerical integration of 3., the dashed line represents the outer solution 3.3. It can be seen that, although ε is not very small, the analytical solution is quite accurate. Figure 3 shows zt with ε 0., in the first quarter period of the orbit, where the outer solution 3.3 and the inner solution 3.6 can be observed dashed line, the numerical

10 0 Journal of Applied Mathematics Inner solution ε = 0. Numerical solution zt) Outer solution GM 3 t Figure 3: The function zt in the first quarter period of the orbit, with ε 0.. Solid line corresponds to numerical integration of 3., the dashed line represents the outer and the inner solution ε = 0. zt) GM 3 t Figure : The function zt with ε 0., the solid line corresponds to numerical integration of 3., the dashed line represents the outer solution 3.3. solution is represented with a solid line. Thus, in the region close to the origin of coordinates, the outer solution should be replaced by the inner solution. In Figure 3 you can see that the matching between the two solution is in a region centered around z/.5. Figure shows the function zt con ε 0.. The solid line corresponds to the solution obtained by numerical integration of 3., the dashed line represents the outer

11 Journal of Applied Mathematics 0 9 zt) Inner solution Outer solution ε = 0. Numerical solution GM 3 t Figure 5: The function zt with ε 0.. The solid line corresponds to numerical integration of 3., the dashed line represents the outer solution and the inner solution. solution 3.3. It can be seen that the numerical and analytical solution are indistinguishable on this scale. This result indicates, that for values of ε, the movement of the third body can be well described by two analytical solutions, the inner and outer solution. Figure 5 shows zt with ε 0. in the region of matching between the inner solution 3.6 and the outer solution 3.3, where it can be noticed that the matching is in a region centered around z/. As the parameter ε 0, the position of the region of matching between the two solutions will be centered approximately in a position /ε, so, in the case ε 0, the orbit is unbounded and the inner solution is the exact solution to the orbit. Figure 6 shows zt with ε 0., in the first quarter period of the orbit, the dashed line represents the uniform approximation 3.9 and the solid line denotes the numerical solution. It can be seen that the numerical and the uniform approximation are indistinguishable on this scale. Defining the relative error between the numerical solution and the uniform approximation as z num t z unif t/z num t, and taking the temporal average of a quarter period, the following values for the mean relative error are obtained Table. Figure 7 shows the period of the orbit as a function of the parameter k defined as k / ε. The solid line corresponds to the period obtained by numerical integration of Macmillan equation.8, the dashed line represents the asymptotic solution 3.3. Figure 7 indicates that the asymptotic expansion 3.3 for the period of the third body matched very well with the numerical calculation of.8, the for the period for values of k close to /. The relative error, defined as T asymp T num /T num, between the asymptotic solution and the numerical solution of the period of oscillations as a function of the parameter k can be seen in Figure 8. The results show that the asymptotic solution for the period of oscillation 3.3 gives good results in a wide range of values of the parameter ε. The asymptotic solution is valid in the region ε, that is, when k is close to /. In this region the relative error is minor than The region of small amplitudes in the movement of the third body,

12 Journal of Applied Mathematics ε = 0. 6 zt) GM 3 t Numerical solution Uniform approximation Figure 6: The function zt in the first quarter period of the orbit, with ε 0., solid line corresponds to numerical integration of 3., the dashed line represents the uniform approximation. Table : elative error between the numerical solution and the uniform approximation. ε Error corresponds to ε, or equivalently k. In this region the relative error of the asymptotic solution is close to 0.0. In the region of small amplitudes in the movement of the third body, MacMillan finds an expression for the period of oscillations of the third body, namely: GM T π 9 3 k 35 ) 6 k.. Figure 9 allows us to see the relative error between the solution of MacMillan. and the numerical solution of the period of oscillations as a function of the parameter k dash line. What is also shown is the relative error between the asymptotic solution 3.3 and the numerical solution of the period of oscillations as a function of the parameter k. Figure 9 indicates that the MacMillan solution. works well for moderate values of k, while the asymptotic solution 3.3 works well in the region where k > Then,

13 Journal of Applied Mathematics GM 3 t k = ε) Numerical solution Asymptotic solution Figure 7: The period of the orbit as a function of the parameter k Error k = ε) Figure 8: elative error between the asymptotic solution and the numerical solution of the period as a function of the parameter k. one can construct the following formula valid for the entire domain 0 <k</ tothe period of the third body: GM 3 π T π k 9 k 35 3/ ) 6 k ) 3 π Γ π 3 5 Γ / k 0.33 k )) 0.33 k<.. This solution has a relative error less than 0.0 for the entire domain 0 <k</.

14 Journal of Applied Mathematics Error k = ε) Figure 9: elative error of the period, MacMillan. dashed line, asymptotic solution 3.3 solid line. 5. Conclusions In this paper, the circular Sitnikov problem is revisited. By matching between asymptotic expansions, approximate analytical expressions for the period and the orbit of the third body in the case of large oscillation are found. The results show that for values ε, the movement of the third body can be well described by two analytical solutions, the inner and outer solution. The analytical asymptotic equation to the period of the third body matched very well with the numerical calculations. It also provides an improved expression for the period of the orbit valid for the entire domain of bounded orbits. The method of matched asymptotic expansions employed in this paper can be extended to more complex problems. For example, the elliptic Sitnikov problem, variable mass of the third body and the relativistic Sitnikov problem. In all these problems, when the speed of the third body is very close to the escape velocity, computational times are very long and it becomes necessary to resort to asymptotic analysis. eferences K. Sitnikov, Existence of oscillating motions for the three-body problem, Doklady Akademii Nauk, vol. 33, no., pp , 960. W. D. MacMillan, An integrable case in the restricted problem of three bodies, Astronomical Journal, vol. 7, no , pp. 3, 9. 3 J. Moser, Stable and andom Motions in Dynamical Systems, Princeton University Press, Princeton, NJ, USA, 973. E. Perdios and V. V. Markellos, Stability and bifurcations of Sitnikov motions, Celestial Mechanics and Dynamical Astronomy, vol., no., pp , Dvorak, Numerical results to the Sitnikov-problem, Celestial Mechanics and Dynamical Astronomy, vol. 56, pp. 7 80, J. Hagel, A new analytic approach to the Sitnikov problem, Celestial Mechanics and Dynamical Astronomy, vol. 53, no. 3, pp. 67 9, J. Liu and Y. S. Sun, On the Sitnikov problem, Celestial Mechanics and Dynamical Astronomy, vol. 9, no. 3, pp , 990.

15 Journal of Applied Mathematics 5 8 K. Wodnar, Predictability, Stability, and Chaos in N-Body Dynamical Systems, Plenum Press, New York, NY, USA, 99, Edited by: A. E. oy. 9 E. Belbruno, J. Llibre, and M. Ollé, On the families of periodic orbits which bifurcate from the circular Sitnikov motions, Celestial Mechanics and Dynamical Astronomy, vol. 60, no., pp. 99 9, S. B. Faruque, Solution of the Sitnikov problem, Celestial Mechanics and Dynamical Astronomy, vol. 87, no., pp , Dvorak, The Sitnikov problem a complete picture of phase space, Publications of the Astronomy Department of the Eötvös University, vol. 9, pp. 9 0, 007. T. Kovács and B. Érdi, The structure of the extended phase space of the Sitnikov problem, Astronomische Nachrichten, vol. 38, pp , Dvorak and Y. S. Sun, The phase space structure of the extended Sitnikov-problem, Celestial Mechanics and Dynamical Astronomy, vol. 67, no., pp , 997. P. Soulis, T. Bountis, and. Dvorak, Stability of motion in the Sitnikov 3-body problem, Celestial Mechanics and Dynamical Astronomy, vol. 99, no., pp. 9 8, J. Hagel, An analytical approach to small amplitude solutions of the extended nearly circular Sitnikov problem, Celestial Mechanics and Dynamical Astronomy, vol. 03, no. 3, pp. 5 66, T. Bountis and K. E. Papadakis, The stability of vertical motion in the N-body circular Sitnikov problem, Celestial Mechanics and Dynamical Astronomy, vol. 0, no. -, pp. 05 5, V. V. Sidorenko, On the circular Sitnikov problem: the alternation of stability and instability in the family of vertical motions, Celestial Mechanics and Dynamical Astronomy, vol. 09, no., pp , 0. 8 S. Di uzza and C. Lhotka, High order normal form construction near the elliptic orbit of the Sitnikov problem, Celestial Mechanics and Dynamical Astronomy, vol., no., pp. 9 6, 0.

16 Advances in Operations esearch Volume 0 Advances in Decision Sciences Volume 0 Mathematical Problems in Engineering Volume 0 Journal of Algebra Probability and Statistics Volume 0 The Scientific World Journal Volume 0 International Journal of Differential Equations Volume 0 Volume 0 Submit your manuscripts at International Journal of Advances in Combinatorics Mathematical Physics Volume 0 Journal of Complex Analysis Volume 0 International Journal of Mathematics and Mathematical Sciences Journal of Stochastic Analysis Abstract and Applied Analysis International Journal of Mathematics Volume 0 Volume 0 Discrete Dynamics in Nature and Society Volume 0 Volume 0 Journal of Journal of Discrete Mathematics Journal of Volume 0 Applied Mathematics Journal of Function Spaces Volume 0 Volume 0 Volume 0 Optimization Volume 0 Volume 0