`Journal of Applied MathematicsVolume 2012 (2012), Article ID 479093, 15 pageshttp://dx.doi.org/10.1155/2012/479093`
Research Article

## Matched Asymptotic Expansions to the Circular Sitnikov Problem with Long Period

1High Temperatures Department, Instituto de Innovación en Minería y Metalurgia, Ahumada 341, 8320181 Santiago, Chile
2Department of Physics, Universidad Gabriela Mistral, Avenida Ricardo Lyon 1177, Santiago, Chile

Received 3 August 2012; Revised 22 November 2012; Accepted 23 November 2012

Copyright © 2012 Marco Rosales-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.

#### Abstract

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.

#### 1. Introduction

The Sitnikov problem [1] 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 1).

Figure 1: Sitnikov problem geometry.

The circular Sitnikov problem is obtained when the two primary bodies describe circular orbits. This situation was originally discussed by MacMillan [2], 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 [1], who proved the existence of oscillatory solutions for this problem. Subsequently, Moser [3] proved the existence the chaotic orbits. Perdios and Markellos [4] examined the stability and bifurcations of straight line motions of the third body. Dvorak [5] studied, by numerical means, motion of the planetoid limited 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 [10] found approximate analytical solutions for small oscillations around the barycenter and moderate eccentricities. Recently, Dvorak [11] studied the complete phase space numerically. Kovács and Érdi [12] 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 [13]. Soulis et al. [14] also studied the situation in which the third body could move off the -axis (generalized extended Sitnikov problem). In 2009, Hagel [15] derived an analytical expression for the perihelion motion of the primaries, when the third body has a finite mass.

Bountis and Papadakis [16] analyzed the problem for an arbitrary number of primary masses, Sidorenko [17] complemented the previous work on the alternation of stability and instability in the family of vertical motions. Recently, Ruzza and Lhotka [18] 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.

#### 2. Equations of Motion

The equation of motion of the third body is where is the distance between the primary bodies and the distance of the third body from the center of mass and representing time. The distance between the primary bodies is given by where is the eccentric anomaly of the primaries, is the eccentricity of the orbits of the primaries. When the primaries move in circular orbits, , (2.1) is reduced to Equation (2.3) represents the circular Sitnikov problem, also known as the problem of MacMillan. Integrating (2.3) with the initial conditions: and we obtain Defining the parameter as Equation (2.4) is transformed into Making the change of variables , and then, allows us to find the following expression to describe the position of the third body [2]: The period is given by This integral can be performed using elliptic integrals, but it is impossible to reverse (2.7) to find an explicit expression for . However, it is possible to find an analytical approximation for the period of the orbit and 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 (2.4), the following value for the escape velocity is found: . In this case, we have .

#### 3. Asymptotic Expansions

When the third body has a speed close to the escape velocity, that is, , (2.6) can be written as

The mechanical total energy () of the third body, can be written in terms of as: , bounded orbit are obtained when . In the case , an unbounded orbit is obtained, 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.1. Near Field

In the regions close to the origin, we have: , thus, the solution of (3.1) can be approximated to the following inner asymptotic expansion: The validity of the expansion (3.2) will be verified in the matching process (see Section 3.3). Introducing the asymptotic expansion (3.2) in (3.1), we obtain the following equation for the leading order and the first order correction: The initial conditions are: and .

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 This integral can be calculated analytically where is the hypergeometric function. From the initial condition , we have , (3.4) can be transformed to more useful as follows: Introducing (3.3) into (3.7), we obtain the following equation for first order perturbations: The solution of (3.8) is given by the expression This integral can be expressed in terms of the hypergeometric function as

##### 3.2. Far Field

In the region far from the origin, the parameter has an important role in (3.1), because defines the maximum distance reached by the third body from the origin: . Moreover, in this region we have and (3.1) can be approximated by This differential equation represents the outer solution. Integrating this equation, the following expression is obtained: This integral can also be calculated analytically, that is: where the value of the constant must be calculated by matching between the inner and outer solutions.

##### 3.3. Matching Procedure

In the region , the inner and outer solutions are valid, then (3.6) and (3.13) must make matching in it.

In the region , (3.6) has the following asymptotic representation: where we have used the following asymptotic expansion for the hypergeometric function: In the region , the asymptotic solution (3.13) can be expanded in power series of , then we obtain the equation of first order in , namely In the leading order, the matching between (3.14) and (3.17) delivers the following value for the constant : Introducing (3.18) into (3.13), the following expression is obtained to in the far field: We may proceed further by combining the inner solution (3.6) and the outer solution (3.18) into a single uniform approximation : where is obtained from (3.14), given: The quarter-period of the orbit is obtained by evaluating (3.18) in , thus we get the following expression for the period:

##### 3.4. Matching in Higher Order

To carry out the asymptotic matching for a first order in is more convenient to write (3.1) in the inner region as follows: Expanding (3.22) in powers of , we obtain: This integral can also be calculated analytically, that is: From the initial condition , we have .

In the region , (3.24) has this asymptotic representation To perform a more accurate estimation of the period, a more accurate estimate of the behavior of for large is necessary, a better approximation of (3.11) for large values of leads to Separating variables in (3.25), and after algebraic manipulation, leads to the following approximate integral: Finally, the following expression is found for improved outer solution: In the region , the asymptotic solution (3.28) can be expanded in power series of , then we obtain the asymptotic expansion: The matching between (3.25) and (3.29) delivers the following value for the constant : Introducing (3.30) into (3.28), the following expression is obtained to in the far field The quarter-period of the orbit is obtained by evaluating (3.31) in , thus obtaining the following expression for the period (including corrections to order ):

#### 4. Results

In this section, the results are expressed using dimensionless variables, the dimensionless time , and the dimensionless distance . In this representation, the results are general and it is not necessary to give numerical values for the constants , and .

Figure 2 presents the function with . The solid line corresponds to the solution obtained by numerical integration of (3.1), the dashed line represents the outer solution (3.31). It can be seen that, although is not very small, the analytical solution is quite accurate.

Figure 2: The function with , the solid line corresponds to numerical integration of (3.1), the dashed line represents the outer solution (3.31).

Figure 3 shows with , in the first quarter period of the orbit, where the outer solution (3.31) and the inner solution (3.6) can be observed (dashed line), the numerical 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 .

Figure 3: The function in the first quarter period of the orbit, with . Solid line corresponds to numerical integration of (3.1), the dashed line represents the outer and the inner solution.

Figure 4 shows the function con . The solid line corresponds to the solution obtained by numerical integration of (3.1), the dashed line represents the outer solution (3.31). 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 4: The function with , the solid line corresponds to numerical integration of (3.1), the dashed line represents the outer solution (3.23).

Figure 5 shows with in the region of matching between the inner solution (3.6) and the outer solution (3.31), where it can be noticed that the matching is in a region centered around .

Figure 5: The function with . The solid line corresponds to numerical integration of (3.1), the dashed line represents the outer solution and the inner solution.

As the parameter , the position of the region of matching between the two solutions will be centered approximately in a position , so, in the case , the orbit is unbounded and the inner solution is the exact solution to the orbit.

Figure 6 shows with , in the first quarter period of the orbit, the dashed line represents the uniform approximation (3.19) and the solid line denotes the numerical solution. It can be seen that the numerical and the uniform approximation are indistinguishable on this scale.

Figure 6: The function in the first quarter period of the orbit, with , solid line corresponds to numerical integration of (3.1), the dashed line represents the uniform approximation.

Defining the relative error between the numerical solution and the uniform approximation as , and taking the temporal average of a quarter period, the following values for the mean relative error are obtained (Table 1).

Table 1: Relative error between the numerical solution and the uniform approximation.

Figure 7 shows the period of the orbit as a function of the parameter defined as . The solid line corresponds to the period obtained by numerical integration of Macmillan equation (2.8), the dashed line represents the asymptotic solution (3.32).

Figure 7: The period of the orbit as a function of the parameter .

Figure 7 indicates that the asymptotic expansion (3.32) for the period of the third body matched very well with the numerical calculation of (2.8), the for the period for values of k close to .

The relative error, defined as , between the asymptotic solution and the numerical solution of the period of oscillations as a function of the parameter can be seen in Figure 8. The results show that the asymptotic solution for the period of oscillation (3.32) gives good results in a wide range of values of the parameter . The asymptotic solution is valid in the region , that is, when is close to . In this region the relative error is minor than 0.005. The region of small amplitudes in the movement of the third body, corresponds to , or equivalently . In this region the relative error of the asymptotic solution is close to 0.04.

Figure 8: Relative error between the asymptotic solution and the numerical solution of the period as a function of the parameter .

In the region of small amplitudes in the movement of the third body, MacMillan [2] finds an expression for the period of oscillations of the third body, namely: Figure 9 allows us to see the relative error between the solution of MacMillan (4.1) and the numerical solution of the period of oscillations as a function of the parameter (dash line). What is also shown is the relative error between the asymptotic solution (3.32) and the numerical solution of the period of oscillations as a function of the parameter .

Figure 9: Relative error of the period, MacMillan ((4.1) dashed line), asymptotic solution ((3.32) solid line).

Figure 9 indicates that the MacMillan solution (4.1) works well for moderate values of , while the asymptotic solution (3.32) works well in the region where . Then, one can construct the following formula valid for the entire domain to the period of the third body: This solution has a relative error less than 0.012 for the entire domain .

#### 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.

#### References

1. K. Sitnikov, “Existence of oscillating motions for the three-body problem,” Doklady Akademii Nauk, vol. 133, no. 2, pp. 303–306, 1960.
2. W. D. MacMillan, “An integrable case in the restricted problem of three bodies,” Astronomical Journal, vol. 27, no. 625-626, pp. 11–13, 1911.
3. J. Moser, Stable and Random Motions in Dynamical Systems, Princeton University Press, Princeton, NJ, USA, 1973.
4. E. Perdios and V. V. Markellos, “Stability and bifurcations of Sitnikov motions,” Celestial Mechanics and Dynamical Astronomy, vol. 42, no. 1–4, pp. 187–200, 1987.
5. R. Dvorak, “Numerical results to the Sitnikov-problem,” Celestial Mechanics and Dynamical Astronomy, vol. 56, pp. 71–80, 1993.
6. J. Hagel, “A new analytic approach to the Sitnikov problem,” Celestial Mechanics and Dynamical Astronomy, vol. 53, no. 3, pp. 267–292, 1992.
7. J. Liu and Y. S. Sun, “On the Sitnikov problem,” Celestial Mechanics and Dynamical Astronomy, vol. 49, no. 3, pp. 285–302, 1990.
8. K. Wodnar, Predictability, Stability, and Chaos in N-Body Dynamical Systems, Plenum Press, New York, NY, USA, 1991, Edited by: A. E. Roy.
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. 1, pp. 99–129, 1994.
10. S. B. Faruque, “Solution of the Sitnikov problem,” Celestial Mechanics and Dynamical Astronomy, vol. 87, no. 4, pp. 353–369, 2003.
11. R. Dvorak, “The Sitnikov problem—a complete picture of phase space,” Publications of the Astronomy Department of the Eötvös University, vol. 19, pp. 129–140, 2007.
12. T. Kovács and B. Érdi, “The structure of the extended phase space of the Sitnikov problem,” Astronomische Nachrichten, vol. 328, pp. 801–804, 2007.
13. R. Dvorak and Y. S. Sun, “The phase space structure of the extended Sitnikov-problem,” Celestial Mechanics and Dynamical Astronomy, vol. 67, no. 1, pp. 87–106, 1997.
14. P. Soulis, T. Bountis, and R. Dvorak, “Stability of motion in the Sitnikov 3-body problem,” Celestial Mechanics and Dynamical Astronomy, vol. 99, no. 2, pp. 129–148, 2007.
15. J. Hagel, “An analytical approach to small amplitude solutions of the extended nearly circular Sitnikov problem,” Celestial Mechanics and Dynamical Astronomy, vol. 103, no. 3, pp. 251–266, 2009.
16. T. Bountis and K. E. Papadakis, “The stability of vertical motion in the $N$-body circular Sitnikov problem,” Celestial Mechanics and Dynamical Astronomy, vol. 104, no. 1-2, pp. 205–225, 2009.
17. 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. 109, no. 4, pp. 367–384, 2011.
18. S. Di Ruzza and C. Lhotka, “High order normal form construction near the elliptic orbit of the Sitnikov problem,” Celestial Mechanics and Dynamical Astronomy, vol. 111, no. 4, pp. 449–464, 2011.