EXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS FOR EACH OF MULTIPLE BODIES

Size: px
Start display at page:

Download "EXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS FOR EACH OF MULTIPLE BODIES"

Transcription

1 The Astrophysical Journal, 687:1294Y1302, 2008 November 10 # The American Astronomical Society. All rights reserved. Printed in U.S.A. EXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS FOR EACH OF MULTIPLE BODIES Da-Zhu Ma, Xin Wu, and Shuang-Ying Zhong Department of Physics, Nanchang University, Nanchang , China; xwu@ncu.edu.cn Received 2008 May 26; accepted 2008 July 9 ABSTRACT For each object of an n-body problem in planetary dynamics, orbital elements except the mean anomaly are directly determined by five independently slow-varying quantities or quasi-integrals, which include the Keplerian energy, the three components of the angular momentum vector, and the z-component of the Laplace vector. The mean anomaly depends on the mean motion specified by the Keplerian energy. Decreasing integration errors of these quasiintegrals at every integration step means improving the accuracy of all the elements to a great extent. Because of this, we take reference values of these quantities in terms of the integral invariant relations as control sources of the errors and then give an extension of Nacozy s idea of manifold correction. The technique is almost the same as the linear transformation method of Fukushima in its explicit validity of correcting all elements, if the adopted basic integrators can give a necessary precision to the stabilizing sources considered. Especially it plays a more important role in significantly suppressing the growth of numerical errors in high eccentricities. Subject headinggs: celestial mechanics methods: n-body simulations methods: numerical 1. INTRODUCTION A pure Keplerian two-body problem is the simplest case in celestial mechanics and dynamical astronomy. In the relative coordinates, there are seven conserved quantities, involving the Keplerian energy K, the three components (P x, P y, and P z )ofthe Laplace vector P, and the three components (L x, L y, and L z )of the angular momentum vector L. They are closely associated with orbital elements. The Keplerian energy determines directly the semimajor axis as well as the mean motion, and the eccentricity is given by the magnitude of the Laplace vector, namely, the Laplace integral P (=jpj). In addition, the inclination and the longitude of the ascending node depend on the magnitude L (=jlj) of the angular momentum vector and the three components L x, L y, and L z. Finally, the argument of perihelion can be obtained from the z-component P z of the Laplace vector and the inclination. It should be worth noting that there exist two identical relations among the seven quantities. One is the orthogonality of two vectors, P = L ¼ 0. The other is the relation P 2 2KL 2 ¼ 2, where is the gravitational constant of the two-body problem ( Fukushima 2003c). So, there are only five independent integrals, K, L x, L y, L z, and P z. These facts imply that the consistency of the five independent integrals over the whole course of the numerical integration leads to the high accuracy of the orbital elements. The manifold correction of Nacozy (1971), as a pioneering work, is considered to be a preferable and convenient tool for arriving at this aim. First, Nacozy considered the case of a known integral. The motion remains on the hypersurface of the integral in phase space, but jumps to a different hypersurface for the presence of a numerical error. In order to compensate for the deviation, Nacozy used a Lagrange multiplier to find the minimum value of the square of the magnitude of a correction vector subjected to the constraint. Adding the least-squares correction to the numerical solution is a basic concept of Nacozy s approach. From the geometrical point of view, the method pulls the integrated orbit back to the original integral hypersurface along the shortest path or the perpendicular to the hypersurface. It is usually called the steepest descent method As emphasized, the corrected orbit does not lie exactly on the hypersurface. In fact, the constraint through the correction is only accurate to the second order of the uncorrected counterpart. Similar to this idea, a scheme on how to keep all invariants (like those in the two-body problem above) during a process of numerical integration was also given by Nacozy. In principle, the manifold correction is able to treat an n-body gravitational system with 10 constants of motion, made of the total energy, the total angular momentum, and integrals of the center of mass, in an inertial coordinate system. Nacozy (1971) found that an application of the correction method for maintaining the 10 integrals of a 25-body problem to numerical integrations of the system is successful at obtaining a significant gain in precision. However, Hairer et al. (1999) pointed out that the approach with the constancy of both the total energy and the total angular momentum fails to work well in a five-body integration of the Sun and four outer planets. Numerical explorations and analytical interpretations of the results like this have appeared in a series of references (e.g., Wu et al. 2006, 2007; Wu & He 2006). The reason stated in these articles is mainly that corrections of integrals, such as the total energy or the total angular momentum, in a full system do not bring great corrections to independent integrals for each subsystem or to individual varying quantities resembling the individual energy or the individual angular momentum for each object. This tells one that it is more important to correct individual nonconstant quantities of each body of an n-body problem, corresponding to the integrals of the above pure two-body problem, than to stabilize the 10 integrals of the n-body system. Unfortunately, Nacozy s original method of manifold correction has difficulty in treating the case of dissipative systems, because it is merely limited to the use of conservative systems with many constraints. If the individual varying quantities that need to be corrected have more accurate reference values at any time, the quantities integrated can be adjusted to the reference values by corrections of the positions and velocities. Namely, the manifold correction method is still valid in this case. For an illustration, the related reference values are obtained from integral invariant relations of the varying quantities

2 EXTENDING NACOZY S APPROACH 1295 (Szebehely & Bettis 1971; Huang & Innanen 1983; Mikkola & Innanen 2002), in which the equations of motion and the time derivatives of the quantities are integrated together. Based on this point, several kinds of extensions to Nacozy s approach have been developed recently. Next, let us introduce and remark on some of them according to different correction aims. Case 1: Corrections of individual Keplerian energy. Clearly, stabilization or correction of the Keplerian energy of each planet or asteroid in an n-body problem is rather significant, since it is the best help in the struggle against Lyapunov s instability of Keplerian motion so as to monitor the accumulation of the in-track error of the orbit ( Baumgarte 1972; Avdyushev 2003). Meanwhile, the accuracy of the semimajor axis and the mean anomaly can be raised drastically. Although the Keplerian energy is no longer invariant, it varies slowly. In this sense, it is still viewed as a quasi-integral. As mentioned above, its reference value at every step is given by the integral invariant relation. This just provides a chance to correct the Keplerian energy. From the theoretical point of view, there are various paths of manifold correction of the Keplerian energy. Two types of corrections are worth noting. One is the rigorous methods for exactly satisfying the relation of the Keplerian energy. The scaling method of Liu & Liao (1988, 1994), the single scaling method of Fukushima (2003a), and the velocity scaling method and the position scaling method of Ma et al. (2008b) are some typical examples. In light of the constraints of a two-body problem, Liu & Liao used two distinct scale factors, and, to adjust the position and velocity of the body, respectively. The two factors obey the relation 2 ¼ 1. Unlike this correction, Fukushima s single scaling method adopts the same scale factor for the integrated position and velocity. The factor is given by solving a certain cubic equation related to the Keplerian energy. As a refined version of the single scaling method, the velocity scaling method (or the position scaling method) uses only a scale factor for the integrated velocity vector (or the integrated position vector). Here, the determination of the scale factor becomes simpler and more explicit. The other is the approximate methods with the least-squares correction of the Keplerian energy. As an example, a direct extension to Nacozy s approach, which Wu et al. (2007) gave, is to only correct three components of the velocity vector. As Ma et al. (2008b) concluded, in spite of the five corrections from different directions, they all are almost effective in the sense of drastically improving the semimajor axis as well as the mean anomaly when the uncorrected integrator is full of sufficient precision. In addition, they are nearly the same at a negligible increase of the computational cost. In particular, both the velocity scaling method and the method of Wu et al. (2007) are the most convenient to apply. Case 2: Corrections of individual Keplerian energy and Laplace integral. A small object is closer to the centered body at the pericenter when orbital eccentricity becomes larger. This gives rise to the fast accumulation of numerical errors. To diminish it, Fukushima (2003b) gave a dual spatial scale transformation to the integrated positions and velocities of each object, where both the Keplerian energy relation and another functional relation associated with the Laplace integral are exactly satisfied. The dual scaling method is also a rigorous method. By using the Laplace integral instead of the functional relation, an approximate method with velocity corrections to the two integrals or quasiintegrals is presented ( Ma et al. 2008a). Still, the approximate scheme is nearly as valid as the rigorous method in raising the accuracy of the semimajor axis and the eccentricity when an uncorrected integrator has sufficient precision. For emphasis, the former is superior to the latter in the correction of eccentricity in general. This is because the former turns out to give the Laplace integral (equivalently, the eccentricity) a direct correction, while the latter only gives an indirect correction to the Laplace integral. Case 3: Corrections of individual Keplerian energy, Laplace integral, and angular momentum vector. In order to suppress the growth of integration errors in the inclination and the longitude of the ascending node, Fukushima (2003c) designed a rotation for consistency of the individual orbital angular momentum vector. Thus, the combination of the dual scaling and the rotation is completely in keeping with the individual Keplerian energy, Laplace integral, and angular momentum vector during the numerical integration, so that the errors in all the orbital elements of each body can be decreased. As the author claimed, this rotational operation is independent of the application of the dual scaling. It means that one can use a dual spatial scale transformation to the integrated positions and velocities with correction of both the Keplerian energy and the Laplace integral, and then carry out a space rotational transformation to the corrected ones for satisfying the angular momentum. From the theoretical point of view, the method with such independent operations of the two transformations, called a two-step correction method, does not keep the simultaneous consistency of these quasi-integrals. For the completeness of the correction theory, Fukushima (2004) constructed such a transformation as the combination of a singleaxis rotation and a linear transformation of the type (x; v)! ½s X x; s V (v x)š so that the three parameters s X, s V,and satisfy exactly the relations associated with the three quantities K, P, andl in a one-step correction. This is a linear transformation method. The main purpose of the present paper is to give an approximate method with the least-squares correction of the five independent quasi-integrals of K, L x, L y, L z, and P z like cases 1 and 2. Our paper is organized as follows. Section 2 gives a new extension to Nacozy s approach with integral invariant relations. Then by some numerical simulations we evaluate and compare it with the linear transformation method of Fukushima (2004) and the method of Wu et al. (2007) in x 3. Finally, x 4 concludes our results. 2. A NEW EXTENSION TO NACOZY S CORRECTION SCHEME Following the manifold correction algorithm of Nacozy (1971) with the integral invariant relations (Szebehely & Bettis1971), we generalize this idea to maintain the consistency of the Keplerian energy, the angular momentum vector, and the z-component of the Laplace vector for the case of a perturbed Keplerian problem. Several details are described in the following The Manifold Correction Algorithm of Nacozy Suppose an m-dimensional dynamical system has s integrals, i (x) ¼ c i (i ¼ 1; :::;s). That is to say, i (x) ¼ i (x) c i ¼ 0 for a true state vector x in the phase space. But, i (h) 6¼ 0due to a numerical solution h yielding errors in the computation. Hence, one gets a nonzero functional vector of the form e ¼ ½ 1 (h); :::; s (h) Š T, where the T superscript indicates transpose. ð1þ

3 1296 MA, WU, & ZHONG Vol. 687 Let E be an s ; m matrix, with the ith row and jth column i (h)/@ j ( j ¼ 1; :::;m). Nacozy (1971) applied Lagrangian multipliers to find a correction vector #h to the solution h, such that the following corrected solution with x ¼ h þ #h ð2þ #h ¼ W 1 E T EW 1 E T 1e ð3þ becomes closer to the true solution x than h because i (x) 2 i (h): Here, W is a weighting matrix, and the matrix E is required to have rank s. In a word, x, adjusting each integral (eq. [1]) with equation (4), is just what one wants. Next, we shall consider the application of the correction scheme The Algorithm Applied to a Pure Keplerian Problem A pure Keplerian problem in the relative coordinates can be simplified to a one-body problem with the Keplerian energy K ¼ 1 2 v 2 r : ð4þ ð5þ Besides this, there are integrals of the Laplace vector and angular momentum vector in the forms P ¼ v < L r r; L ¼ r < v: Here, we specify that r ¼ (x; y; z), v ¼ (ẋ; ẏ; ż), G(M þ m), and r ¼jrj denote the position, velocity, gravitational constant, and radius, respectively. It should be emphasized that these integrals determine directly five orbital elements, namely, the semimajor axis a, the eccentricity e, the inclination I, the longitude of ascending node, and the argument of perihelion!. Relations between the integrals and the elements are described as ð6þ ð7þ a ¼ 2K ; e ¼ P ; ð8þ L z L x I ¼ arccos ; ¼ arctan ; ð9þ L L y P z! ¼ arcsin : ð10þ e sin I For an illustration, the mean anomaly M, which includes information on the location in the orbit, is associated with the mean motion specified by the Keplerian energy. The orientation of on the orbital plane should be given by the signs of L x and L y. Similarly, the orientation of! is based on the sign of sin! and that of P x cos þ P y sin. Therefore, the conservation of K, L x, L y, L z, and P z means the improvement of the five elements (note that if K and L are controlled, then P is also controlled automatically). Naturally, the mean anomaly can be improved. Matching with the quantities in x 2.1, the state vector x ¼ (x; y; z; ẋ; ẏ; ż) and five constants, c 1 ;:::;c 5, are respectively K 0, L 0x, L 0y, L 0z, and P 0z,as initial values of K, L x, L y, L z, and P z. Now, the application of Nacozy s approach to the case becomes easy by y E z 0 1 K K 0 L x L 0x e ¼ L y L 0y ; ð11þ B L z L 0z A P z P y : z z A The concrete expressions of all the components of E are of the following forms, 0 x y z 1 r 3 r 3 r 3 ẋ ẏ ż 0 ż ẏ 0 z y E ¼ ż 0 ẋ z 0 x ; ð13þ B ẏ ẋ 0 y x 0 A Q 1 Q 2 Q 3 Q 4 Q 5 Q 6 with Q 1 ¼ ẋżþxz/r 3, Q 2 ¼ ẏż þ yz/r 3, Q 3 ¼ ẋ 2 þ ẏ 2 (x 2 þy 2 )/r 3, Q 4 ¼ 2zẋ xż, Q 5 ¼ yżþ2zẏ,andq 6 ¼ xẋ yẏ. Obviously, the matrix E is of rank 5. It should be mentioned that the components of the E matrix are of different dimensions, so that the resulting minimization process like the steepest descent method becomes disastrous. It is important to adopt the choice of nondimensionalization in the actual application of Nacozy s approach. For nondimensionalization, we carry out scale transformations, such as r! r/r 0, v! v/v 0, t! t/t, etc., where constants r 0, v 0, and T are the initial radius, the magnitude of the initial velocity, and the average orbital period, respectively. The weight matrix W is a unit matrix through the operation of nondimensionalization. Note that the original quantities are still used before and after the implementation of the manifold correction. Of course, the correction method can further be extended to treat the perturbed case The Algorithm Extended to a Perturbed Keplerian Problem When a perturbing acceleration a exists, the perturbed one body evolves according to the equation of motion dv dt ¼ r 3 r þ a: ð14þ

4 No. 2, 2008 EXTENDING NACOZY S APPROACH 1297 Fig. 1. Variations of errors in all elements of a pure Keplerian orbit with time. In this case, K, P, and L in equations (5)Y(7) remain no longer invariant, but vary slowly with time in the following forms (Fukushima 2003c), dp dt dk dt ¼ v = a; ¼ 2(a = v)r (r = a)v (r = v)a; ð16þ dl dt ¼ r < a: ð15þ ð17þ Equations (15) Y(17) are so-called integral invariant relations about these quasi-integrals of the Keplerian energy, the Laplace vector, and the angular momentum vector. Note that the righthand sides of equations (15)Y(17) are small quantities, so the left-hand sides would use K ¼ K K 0, P ¼ P P 0,and L ¼ L L 0 rather than K, P, and L in order to greatly reduce the round-off errors. The values of K, P, and L from direct integration of equations (14)Y(17) should be more accurate than the corresponding values of them by equations (5)Y(7) with a numerical solution of equation (14). We mark them as K, P, and L in sequence. They are regarded as standard or reference values for correcting the quasi-integrals. Let the constants K 0, L 0x, L 0y, L 0z, and P 0z in equation (11) give way to K, L x, L y, L z, and Pz, respectively. Then the manifold correction algorithm of Nacozy does still work for the perturbed Keplerian problem. Hereafter, we call this extension method 1 (M1). It is worth emphasizing that the generalization, M1, is also suitable to the case of multiple bodies. As stated above, M1 gives corrections to the integrated velocities and positions, which approximately and synchronously satisfy the relations of the Keplerian energy, the three components of the angular momentum vector, and the z-component of the Laplace vector like equation (4). Before this method, a

5 1298 MA, WU, & ZHONG rigorous method with correction of the Keplerian energy, the magnitude of the angular momentum vector, and the Laplace integral was the linear transformation method ( M2) of Fukushima (2004). See this article for more information. In sum, a typical differencebetweenm1andm2liesinthatm1isapproximate to maintain the consistency of the relations of the five quantities, but M2 is rigorous to satisfy the three quantities. Thus, it is strongly desired to evaluate the validity of M1 and to especially know whether the difference has an effect on enhancing the quality of orbit integrations. In-depth numerical explorations will be necessary. For an overall comparison of various manifold corrections in an n-body problem, we consider the approximately individual Keplerian energycorrectionmethod(m3) of Wu et al. (2007). 3. EVALUATION OF THE METHOD In our numerical tests, we mainly take pure Keplerian twobody problems with various eccentricities as test models of numerical simulations. In addition, we pay attention to the case of a three-body problem made of the Sun, Jupiter, and Saturn Pure Keplerian Orbits Let us begin with the simplest case, numerical integration of a pure Keplerian problem (eq. [5]) with a massless particle. We select an orbit with initial conditions of x ¼ 1, y ¼ 0:005, z ¼ 0:005, ẋ ¼ 0:005, ẏ ¼ 1, and ż ¼ 0:005, and a fifth-order Runge-Kutta (RK5) algorithm with an invariant time step being 1/100 of the orbital period T, asabasicintegrator. Errors in all the Keplerian elements of RK5 and its corrections, M1, M2, and M3, are plotted in Figure 1. It is clear in Figure 1a that RK5 gains the linear growth of the error in the semimajor axis a with respect to time, but any of the three corrections achieves the error on the order of the machine epsilon. This displays that the three corrections have almost the same performance in controlling the Keplerian energy. Similarly, they should not produce a great difference in the error of the mean anomaly M (Fig. 1f ). However, the three corrections give distinct answers to the errors in the eccentricity e and the argument of pericenter!,as shown in Figures 1b and 1d. Both M1 and M2 can greatly improve the accuracy of the eccentricity. In addition, M2 is slightly better than M1, but M3 (like RK5) is the worst. This shows that the adjustments of K and L lead to the automatic adjustment of P for M1, and it is successful at correcting! by controlling the z-component of the Laplace vector. As to M3, it is only used to correct the Keplerian energy. On the other hand, it can be seen from Figures 1c and 1e that both M1 and M2 play the same role in the adjustment of the orbital inclination I and the longitude of the ascending node. This shows that M1 and M2 are equally valid in keeping with the angular momentum vector. Besides the accuracies of the orbital elements, we observe the relative position error of the orbit for each of the above methods. At the end of integration, both M1 and M3 raise the orbital accuracy by 4 orders of magnitude, as shown in Figure 2. As a slight difference between M1 and M3, M1 becomes slightly better. In particular, M2 is the best. Now, let us concentrate on the eccentricity dependence of element error for a Keplerian orbit when one of the corrections is applied. We fix initial orbital elements a ¼ 2 AU, I ¼ 20, ¼ 50,! ¼ 30, and M ¼ 40, but vary e from 0 to 0.7 at intervals of The integrator used is still the same as that in Figure 1, and each orbit stops computing after 10 4 periods. Plotted are the errors in the orbital elements as functions of eccentricity Fig. 2. Relative position errors for the same orbit from Fig. 1. for the numerical schemes above. Like Figure 1, Figure 3 illustrates that M1, M2, and M3 have the same effectiveness in improving a and M. Compared with M3, both M1 and M2 drastically decrease the errors in I and, down to the limit of the machine epsilon. There are not even any differences in the errors of e and! between M1 and M2. On the other hand, the magnitude of the error in a does not depend on the eccentricity for any of the three corrections; neither do the magnitudes of the errors in I and for M1 and M2. Especially for M1 and M2, the accuracy of e and! seems to become better as e gets larger. But the error in M increases with e for each correction. In order to evaluate the case of the high eccentricities near 1, we use a fifth- to sixth-order Runge- Kutta-Fehlberg algorithm of variable time step as a basic test tool and obtain results (not plotted) similar to those of Figure 3. But there is only one difference, namely, that the error of M does not depend on the eccentricity at all for any correction. For the computations of various orbits in Figure 3, let us provide a comparison of computational cost, say, CPU times. At the end of numerical integration for all 70 orbits, the CPU times of RK5, M1, M2, and M3 are 94, 193, 103, and 96 s, respectively. The costs of the correction methods M2 and M3 increase slightly. As expected, the computational efficiency of M3 is the best, and that of M1 is the worst because M1 needs some additional time to solve the inverse matrix (EE T ) 1. A series of numerical experiments about manifold corrections are carried out in the case that the uncorrected integrators, e.g., RK5, give the accuracy of a certain integral on the order of not larger than 10 8 in a double-precision environment. These tests turn out to show that M1 and M2 are almost the same in improving the orbital elements. Without loss of generality, we shall compare M1 with M2 using real physical models Three-Body Problem We apply RK5 and its corrections M1, M2, and M3 to a general three-body problem of the Sun, Jupiter, and Saturn as an experiment model. The related physical parameters and initial conditions are quoted from those at J in JPL s planetary ephemeris, DE405. The fixed time step adopted is days, about 1/120 of Jupiter s orbital period. A 12th-order Cowell method is used to get a higher precision reference orbit for the construction of the errors in the orbital elements and the positions of each planet. In the heliocentric coordinate system, individual quantities K, P, andl givenbyequations(5)y(7)

6 Fig. 3. Eccentricity dependence of element error for a Keplerian orbit with the uncorrected integrator, RK5, of a fixed time step.

7 Fig. 4. Errors of all elements for Jupiter in the three-body problem of the Sun, Jupiter, and Saturn.

8 EXTENDING NACOZY S APPROACH 1301 Fig. 5. Same as Fig. 4, but for Saturn. vary slowly with time and satisfy equations (15)Y(17). As stated above, M1, M2, and M3 are still able to work in this sense. Unlike Figure 1, Figures 4 and 5 show that the three corrections arrive at the same effect of improving all the orbital elements for each planet by a great measure. The main agreement of M1, M2, and M3 lies in the errors in a, M, I, and. But M1 as well as M2 is better than M3 in correcting the errors of the eccentricity and the argument of pericenter (see Figs. 4b and 4d or Figs. 5b and 5d for more details). In a word, Figures 4 and 5 show that M1 is basically the same as M2 in the improvement of each orbital element. Here, three main points are worth noting. 1. The correction of a (besides M ) is the most valid of all the orbital elements. The reason can be found from equations (15)Y (17). Generally speaking, a is a small quantity, and jvj is rather smaller than jrj. Thus, the right hand-side of equation (16) is larger than that of equation (15), but smaller than that of equation (17). This seems to show that the accuracy of the reference value K, the accuracy of the reference value P, and the accuracy of the reference value L change from high to low. 2. M1 with M2 does not have an advantage over M3 in controlling the errors of I and, because L is nearly equal to the numerical value L given by equation (7). 3. The validity of orbital corrections for Jupiter is more explicit than that for Saturn. This is because the orbital accuracies for Jupiter are worse than those for Saturn in the uncorrected case with the same step size. The relative position error of each planet in Figure 6 is unlike that of the one-body in Figure 2. As an emphasis, here the three

9 1302 MA, WU, & ZHONG Fig. 6. Relative position errors of two planets in the three-body problem. corrections, M1, M2, and M3, are nearly the same after a short time of integration. 4. CONCLUSIONS Based on the integral invariant relations and the manifold correction of Nacozy, some correction schemes, with the adjustments of the varying Keplerian energy, Laplace integral, and/or angular momentum vector for a perturbed Keplerian one-body problem or each of multiple bodies, were recently developed along two different paths. One is the rigorous methods. One of the typical works is the linear transformation method (M2) of Fukushima. The other is the approximate methods with the leastsquares corrections. Along this line, the present paper has given a one-step correction method ( M1), where five independently slowvarying quantities or quasi-integrals, the Keplerian energy, the three components of the angular momentum vector, and the z-component of the Laplace vector, are approximately and simultaneously satisfied. The consistency of these quasi-integrals means the improvement of all the orbital elements of each object. By comparing between M1 and M2, we found that they have almost the same performance in improving all the orbital elements when the uncorrected integrators give the chosen integral a necessary precision on the order of not larger than 10 8 in a double-precision environment. The authors thank the referee for valuable suggestions. The research is supported by the Natural Science Foundation of China under contract It is also supported by the Science Foundation of Jiangxi Province ( ) and the Program for Innovative Research Team of Nanchang University. REFERENCES Avdyushev, V. A. 2003, Celest. Mech. Dyn. Astron., 87, 383 Liu, L., & Liao, X.-H , Celest. Mech. Dyn. Astron., 59, 221 Baumgarte, J. 1972, Celest. Mech., 5, 490 Ma, D.-Z., Wu, X., & Liu, F. Y. 2008a, Int. J. Mod. Phys. C, in press Fukushima, T. 2003a, AJ, 126, 1097 Ma, D.-Z., Wu, X., & Zhu. J.-F. 2008b, NewA, 13, b, AJ, 126, 2567 Mikkola, S., & Innanen, K. 2002, AJ, 124, c, AJ, 126, 3138 Nacozy, P. E. 1971, Ap&SS, 14, , AJ, 127, 3638 Szebehely, V., & Bettis, D. G. 1971, Ap&SS, 13, 365 Hairer, E., Lubich, C., & Wanner, G. 1999, Geometric Numerical Integration Wu, X., & He, J.-Z. 2006, Int. J. Mod. Phys. C, 17, 1613 ( Berlin: Springer) Wu, X., Huang, T.-Y., Wan, X.-S., & Zhang, H. 2007, AJ, 133, 2643 Huang, T.-Y., & Innanen, K. 1983, AJ, 88, 870 Wu, X., Zhu, J. F., He, J. Z., & Zhang, H. 2006, Comput. Phys. Commun., 175, Liu, L., & Liao, X.-H. 1988, Chinese Astron. Astrophys., 12, 26 15

EFFICIENT INTEGRATION OF HIGHLY ECCENTRIC ORBITS BY QUADRUPLE SCALING FOR KUSTAANHEIMO-STIEFEL REGULARIZATION

EFFICIENT INTEGRATION OF HIGHLY ECCENTRIC ORBITS BY QUADRUPLE SCALING FOR KUSTAANHEIMO-STIEFEL REGULARIZATION The Astronomical Journal, 128:3108 3113, 2004 December # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFICIENT INTEGRATION OF HIGHLY ECCENTRIC ORBITS BY QUADRUPLE SCALING

More information

EFFICIENT ORBIT INTEGRATION BY SCALING FOR KEPLER ENERGY CONSISTENCY Toshio Fukushima

EFFICIENT ORBIT INTEGRATION BY SCALING FOR KEPLER ENERGY CONSISTENCY Toshio Fukushima The Astronomical Journal, 126:1097 1111, 2003 August # 2003. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFICIENT ORBIT INTEGRATION BY SCALING FOR KEPLER ENERGY CONSISTENCY

More information

Orbit Determination Using Satellite-to-Satellite Tracking Data

Orbit Determination Using Satellite-to-Satellite Tracking Data Chin. J. Astron. Astrophys. Vol. 1, No. 3, (2001 281 286 ( http: /www.chjaa.org or http: /chjaa.bao.ac.cn Chinese Journal of Astronomy and Astrophysics Orbit Determination Using Satellite-to-Satellite

More information

integral invariant relations is not limited to one or two such

integral invariant relations is not limited to one or two such The Astronomical Journal, 126:3138 3142, 2003 December # 2003. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFICIENT ORBIT INTEGRATION BY SCALING AND ROTATION FOR CONSISTENCY

More information

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic Orbits in Rotating Second Degree and Order Gravity Fields Chin. J. Astron. Astrophys. Vol. 8 (28), No. 1, 18 118 (http://www.chjaa.org) Chinese Journal of Astronomy and Astrophysics Periodic Orbits in Rotating Second Degree and Order Gravity Fields Wei-Duo Hu

More information

Dynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly

Dynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly First Kepler s Law The secondary body moves in an elliptical orbit, with the primary body at the focus Valid for bound orbits with E < 0 The conservation of the total energy E yields a constant semi-major

More information

The two-body Kepler problem

The two-body Kepler problem The two-body Kepler problem set center of mass at the origin (X = 0) ignore all multipole moments (spherical bodies or point masses) define r := r 1 r 2,r:= r,m:= m 1 + m 2,µ:= m 1 m 2 /m reduces to effective

More information

Extending the Patched-Conic Approximation to the Restricted Four-Body Problem

Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean

More information

Lecture 13. Gravity in the Solar System

Lecture 13. Gravity in the Solar System Lecture 13 Gravity in the Solar System Guiding Questions 1. How was the heliocentric model established? What are monumental steps in the history of the heliocentric model? 2. How do Kepler s three laws

More information

A Kozai-resonating Earth quasi-satellite

A Kozai-resonating Earth quasi-satellite MNRASL 437, L85 L89 (2014) Advance Access publication 2013 November 13 doi:10.1093/mnrasl/slt147 A Kozai-resonating Earth quasi-satellite M. Connors Athabasca University Observatories, Athabasca AB T9S

More information

A study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1

A study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1 Astronomy & Astrophysics manuscript no. Eris c ESO 2013 March 27, 2013 A study upon Eris I. Describing and characterizing the orbit of Eris around the Sun I. Breda 1 Faculty of Sciences (FCUP), University

More information

Restricted three body problems in the Solar System: simulations

Restricted three body problems in the Solar System: simulations Author:. Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain. Advisor: Antoni Benseny i Ardiaca. Facultat de Matemàtiques, Universitat de Barcelona, Gran Via de les Corts

More information

Calculation of Earth s Dynamic Ellipticity from GOCE Orbit Simulation Data

Calculation of Earth s Dynamic Ellipticity from GOCE Orbit Simulation Data Available online at www.sciencedirect.com Procedia Environmental Sciences 1 (1 ) 78 713 11 International Conference on Environmental Science and Engineering (ICESE 11) Calculation of Earth s Dynamic Ellipticity

More information

arxiv: v1 [astro-ph.ep] 3 Apr 2018

arxiv: v1 [astro-ph.ep] 3 Apr 2018 Astronomy& Astrophysics manuscript no. zanardi_ c ESO 28 September 2, 28 The role of the general relativity on icy body reservoirs under the effects of an inner eccentric Jupiter M. Zanardi, 2, G. C. de

More information

AST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000).

AST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000). AST111 PROBLEM SET 2 SOLUTIONS Home work problems 1. Angles on the sky and asteroid motion An asteroid is observed at two different times. The asteroid is located at RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch

More information

1 Summary of Chapter 2

1 Summary of Chapter 2 General Astronomy (9:61) Fall 01 Lecture 7 Notes, September 10, 01 1 Summary of Chapter There are a number of items from Chapter that you should be sure to understand. 1.1 Terminology A number of technical

More information

MOTIONS 1. INTRODUCTION 2. KOLMOGOROV-SMIRNOV REGULARIZED ORBITAL

MOTIONS 1. INTRODUCTION 2. KOLMOGOROV-SMIRNOV REGULARIZED ORBITAL THE ASTRONOMICAL JOURNAL, 120:3333È3339, 2000 December ( 2000. The American Astronomical Society. All rights reserved. Printed in U.S.A. LONG-TERM INTEGRATION ERROR OF KUSTAANHEIMO-STIEFEL REGULARIZED

More information

Symplectic Correctors for Canonical Heliocentric N-Body Maps

Symplectic Correctors for Canonical Heliocentric N-Body Maps Symplectic Correctors for Canonical Heliocentric N-Body Maps J. Wisdom Massachusetts Institute of Technology, Cambridge, MA 02139 wisdom@poincare.mit.edu Received ; accepted 2 ABSTRACT Symplectic correctors

More information

Analysis of frozen orbits for solar sails

Analysis of frozen orbits for solar sails Trabalho apresentado no XXXV CNMAC, Natal-RN, 2014. Analysis of frozen orbits for solar sails J. P. S. Carvalho, R. Vilhena de Moraes, Instituto de Ciência e Tecnologia, UNIFESP, São José dos Campos -

More information

The Heliocentric Model of Copernicus

The Heliocentric Model of Copernicus Celestial Mechanics The Heliocentric Model of Copernicus Sun at the center and planets (including Earth) orbiting along circles. inferior planets - planets closer to Sun than Earth - Mercury, Venus superior

More information

The orbit evolution of 32 plutinos over 100 million year

The orbit evolution of 32 plutinos over 100 million year A&A 368, 700 705 (2001) DOI: 10.1051/0004-6361:20010056 c ESO 2001 Astronomy & Astrophysics The orbit evolution of 32 plutinos over 100 million year X.-S. Wan 1 and T.-Y. Huang 1,2 1 Astronomy Department,

More information

The Kepler Problem and the Isotropic Harmonic Oscillator. M. K. Fung

The Kepler Problem and the Isotropic Harmonic Oscillator. M. K. Fung CHINESE JOURNAL OF PHYSICS VOL. 50, NO. 5 October 01 The Kepler Problem and the Isotropic Harmonic Oscillator M. K. Fung Department of Physics, National Taiwan Normal University, Taipei, Taiwan 116, R.O.C.

More information

The precession of the perihelion of Mercury explained by Celestial Mechanics of Laplace Valdir Monteiro dos Santos Godoi

The precession of the perihelion of Mercury explained by Celestial Mechanics of Laplace Valdir Monteiro dos Santos Godoi The precession of the perihelion of Mercury explained by Celestial Mechanics of Laplace Valdir Monteiro dos Santos Godoi valdir.msgodoi@gmail.com ABSTRACT We calculate in this article an exact theoretical

More information

Reducing round-off errors in symmetric multistep methods

Reducing round-off errors in symmetric multistep methods Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,

More information

The Three Body Problem

The Three Body Problem The Three Body Problem Joakim Hirvonen Grützelius Karlstad University December 26, 2004 Department of Engineeringsciences, Physics and Mathematics 5p Examinator: Prof Jürgen Füchs Abstract The main topic

More information

STABILITY OF HYPOTHETICAL TROJAN PLANETS IN EXOPLANETARY SYSTEMS

STABILITY OF HYPOTHETICAL TROJAN PLANETS IN EXOPLANETARY SYSTEMS STABILITY OF HYPOTHETICAL TROJAN PLANETS IN EXOPLANETARY SYSTEMS Bálint Érdi, Georgina Fröhlich, Imre Nagy and Zsolt Sándor Department of Astronomy Loránd Eötvös University Pázmány Péter sétány 1/A H-1117

More information

Towards stability results for planetary problems with more than three bodies

Towards stability results for planetary problems with more than three bodies Towards stability results for planetary problems with more than three bodies Ugo Locatelli [a] and Marco Sansottera [b] [a] Math. Dep. of Università degli Studi di Roma Tor Vergata [b] Math. Dep. of Università

More information

Motion under the Influence of a Central Force

Motion under the Influence of a Central Force Copyright 004 5 Motion under the Influence of a Central Force The fundamental forces of nature depend only on the distance from the source. All the complex interactions that occur in the real world arise

More information

AST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements

AST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends

More information

Dynamics of Charged Particles Around a Magnetized Black Hole with Pseudo Newtonian Potential

Dynamics of Charged Particles Around a Magnetized Black Hole with Pseudo Newtonian Potential Commun. Theor. Phys. 60 (2013) 433 438 Vol. 60, No. 4, October 15, 2013 Dynamics of Charged Particles Around a Magnetized Black Hole with Pseudo Newtonian Potential WANG Ying ( Ä), 1,2 WU Xin ( ), 1, and

More information

Long simulations of the Solar System: Brouwer s Law and chaos

Long simulations of the Solar System: Brouwer s Law and chaos ANZIAM J. 46 (E) pp.c1086 C1103, 2005 C1086 Long simulations of the Solar System: Brouwer s Law and chaos K. R. Grazier W. I. Newman James M. Hyman Philip. W. Sharp (received 28 October 2004, revised 18

More information

Lecture XIX: Particle motion exterior to a spherical star

Lecture XIX: Particle motion exterior to a spherical star Lecture XIX: Particle motion exterior to a spherical star Christopher M. Hirata Caltech M/C 350-7, Pasadena CA 95, USA Dated: January 8, 0 I. OVERVIEW Our next objective is to consider the motion of test

More information

Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune

Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune PASJ: Publ. Astron. Soc. Japan 54, 493 499, 2002 June 25 c 2002. Astronomical Society of Japan. Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune Tetsuharu FUSE Subaru Telescope,

More information

Tides and Lagrange Points

Tides and Lagrange Points Ast111, Lecture 3a Tides and Lagrange Points Arial view of Tidal surge in Alaska (S. Sharnoff) Tides Tidal disruption Lagrange Points Tadpole Orbits and Trojans Tidal Bulges Tides Tidal Force The force

More information

Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws

Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial

More information

Research Paper. Trojans in Habitable Zones ABSTRACT

Research Paper. Trojans in Habitable Zones ABSTRACT ASTROBIOLOGY Volume 5, Number 5, 2005 Mary Ann Liebert, Inc. Research Paper Trojans in Habitable Zones RICHARD SCHWARZ, 1 ELKE PILAT-LOHINGER, 1 RUDOLF DVORAK, 1 BALINT ÉRDI, 2 and ZSOLT SÁNDOR 2 ABSTRACT

More information

orbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG

orbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG A Science Enrichment Programme for Secondary 3-4 Students : Teaching and Learning Resources the and Spacecrafts orbits Moon, Planets Calculating the 171 of by Dr. Shiu-Sing TONG 172 Calculating the orbits

More information

ACCURATE KEPLER EQUATION SOLVER WITHOUT TRANSCENDENTAL FUNCTION EVALUATIONS

ACCURATE KEPLER EQUATION SOLVER WITHOUT TRANSCENDENTAL FUNCTION EVALUATIONS AAS ACCURATE KEPLER EQUATION SOLVER WITHOUT TRANSCENDENTAL FUNCTION EVALUATIONS Adonis Pimienta-Peñalver and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 14260-4400

More information

ASTRONOMY AND ASTROPHYSICS. Escape with the formation of a binary in two-dimensional three-body problem. I. Navin Chandra and K.B.

ASTRONOMY AND ASTROPHYSICS. Escape with the formation of a binary in two-dimensional three-body problem. I. Navin Chandra and K.B. Astron. Astrophys. 46 65 66 1999 ASTRONOMY AND ASTROPHYSICS Escape with the formation of a binary in two-dimensional three-body problem. I Navin Chandra and K.B. Bhatnagar Centre for Fundamental Research

More information

Investigation of evolution of orbits similar to that of (4179) Toutatis

Investigation of evolution of orbits similar to that of (4179) Toutatis Investigation of evolution of orbits similar to that of (4179) Toutatis Wm. Robert Johnston May 1994 Abstract This project compared the orbital evolution of asteroid (4179) Toutatis to that of several

More information

arxiv:physics/ v2 [physics.gen-ph] 2 Dec 2003

arxiv:physics/ v2 [physics.gen-ph] 2 Dec 2003 The effective inertial acceleration due to oscillations of the gravitational potential: footprints in the solar system arxiv:physics/0309099v2 [physics.gen-ph] 2 Dec 2003 D.L. Khokhlov Sumy State University,

More information

KEPLER S LAWS OF PLANETARY MOTION

KEPLER S LAWS OF PLANETARY MOTION KEPLER S LAWS OF PLANETARY MOTION In the early 1600s, Johannes Kepler culminated his analysis of the extensive data taken by Tycho Brahe and published his three laws of planetary motion, which we know

More information

Eccentricity pumping of a planet on an inclined orbit by a disc

Eccentricity pumping of a planet on an inclined orbit by a disc Mon. Not. R. Astron. Soc. 44, 49 414 21) doi:1.1111/j.1365-2966.21.16295.x Eccentricity pumping of a planet on an inclined orbit by a disc Caroline Terquem 1,2 and Aikel Ajmia 1 1 Institut d Astrophysique

More information

Importance of the study of extrasolar planets. Exoplanets Introduction. Importance of the study of extrasolar planets

Importance of the study of extrasolar planets. Exoplanets Introduction. Importance of the study of extrasolar planets Importance of the study of extrasolar planets Exoplanets Introduction Planets and Astrobiology (2017-2018) G. Vladilo Technological and scientific spin-offs Exoplanet observations are driving huge technological

More information

Hill's Approximation in the Three-Body Problem

Hill's Approximation in the Three-Body Problem Progress of Theoretical Physics Supplement No. 96, 1988 167 Chapter 15 Hill's Approximation in the Three-Body Problem Kiyoshi N AKAZA W A and Shigeru IDA* Department of Applied Physics, T.okyo Institute

More information

Theory of mean motion resonances.

Theory of mean motion resonances. Theory of mean motion resonances. Mean motion resonances are ubiquitous in space. They can be found between planets and asteroids, planets and rings in gaseous disks or satellites and planetary rings.

More information

Orbital Motion in Schwarzschild Geometry

Orbital Motion in Schwarzschild Geometry Physics 4 Lecture 29 Orbital Motion in Schwarzschild Geometry Lecture 29 Physics 4 Classical Mechanics II November 9th, 2007 We have seen, through the study of the weak field solutions of Einstein s equation

More information

Physics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top

Physics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem

More information

Three objects; 2+1 problem

Three objects; 2+1 problem Three objects; 2+1 problem Having conquered the two-body problem, we now set our sights on more objects. In principle, we can treat the gravitational interactions of any number of objects by simply adding

More information

F. Marzari, Dept. Physics, Padova Univ. Planetary migration

F. Marzari, Dept. Physics, Padova Univ. Planetary migration F. Marzari, Dept. Physics, Padova Univ. Planetary migration Standard model of planet formation based on Solar system exploration Small semimajor axes Large eccentricities The standard model Protostar +Disk

More information

On the Estimated Precession of Mercury s Orbit

On the Estimated Precession of Mercury s Orbit 1 On the Estimated Precession of Mercury s Orbit R. Wayte. 9 Audley Way, Ascot, Berkshire, SL5 8EE, England, UK e-mail: rwayte@googlemail.com Research Article, Submitted to PMC Physics A 4 Nov 009 Abstract.

More information

Satellite Communications

Satellite Communications Satellite Communications Lecture (3) Chapter 2.1 1 Gravitational Force Newton s 2nd Law: r r F = m a Newton s Law Of Universal Gravitation (assuming point masses or spheres): Putting these together: r

More information

Gravitational Effects of Rotating Bodies

Gravitational Effects of Rotating Bodies Mod. Phys. Lett. A16(2001) 789 794 Gravitational Effects of Rotating Bodies Yuan-Zhong Zhang a, Jun Luo b, and Yu-Xin Nie c a Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2735,

More information

THE exploration of planetary satellites is currently an active area

THE exploration of planetary satellites is currently an active area JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 9, No. 5, September October 6 Design of Science Orbits About Planetary Satellites: Application to Europa Marci E. Paskowitz and Daniel J. Scheeres University

More information

Keplerian Elements Tutorial

Keplerian Elements Tutorial Keplerian Elements Tutorial This tutorial is based on the documentation provided with InstantTrack, written by Franklin Antonio, N6NKF. Satellite Orbital Elements are numbers that tell us the orbit of

More information

DIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES

DIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES DIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES KLEOMENIS TSIGANIS and HARRY VARVOGLIS Section of Astrophysics, Astronomy and Mechanics, Department of Physics, University of Thessaloniki, 540 06 Thessaloniki,

More information

Satellite meteorology

Satellite meteorology GPHS 422 Satellite meteorology GPHS 422 Satellite meteorology Lecture 1 6 July 2012 Course outline 2012 2 Course outline 2012 - continued 10:00 to 12:00 3 Course outline 2012 - continued 4 Some reading

More information

PHYSICS. Chapter 13 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.

PHYSICS. Chapter 13 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc. PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 13 Lecture RANDALL D. KNIGHT Chapter 13 Newton s Theory of Gravity IN THIS CHAPTER, you will learn to understand the motion of satellites

More information

AS3010: Introduction to Space Technology

AS3010: Introduction to Space Technology AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital

More information

TRANSFER TO THE COLLINEAR LIBRATION POINT L 3 IN THE SUN-EARTH+MOON SYSTEM

TRANSFER TO THE COLLINEAR LIBRATION POINT L 3 IN THE SUN-EARTH+MOON SYSTEM TRANSFER TO THE COLLINEAR LIBRATION POINT L 3 IN THE SUN-EARTH+MOON SYSTEM HOU Xi-yun,2 TANG Jing-shi,2 LIU Lin,2. Astronomy Department, Nanjing University, Nanjing 20093, China 2. Institute of Space Environment

More information

ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS

ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.

More information

Resonance and chaos ASTRONOMY AND ASTROPHYSICS. II. Exterior resonances and asymmetric libration

Resonance and chaos ASTRONOMY AND ASTROPHYSICS. II. Exterior resonances and asymmetric libration Astron. Astrophys. 328, 399 48 (1997) ASTRONOMY AND ASTROPHYSICS Resonance and chaos II. Exterior resonances and asymmetric libration O.C. Winter 1,2 and C.D. Murray 2 1 Grupo de Dinâmica Orbital e Planetologia,

More information

Advanced Newtonian gravity

Advanced Newtonian gravity Foundations of Newtonian gravity Solutions Motion of extended bodies, University of Guelph h treatment of Newtonian gravity, the book develops approximation methods to obtain weak-field solutions es the

More information

ACCURATE KEPLER EQUATION SOLVER WITHOUT TRANSCENDENTAL FUNCTION EVALUATIONS

ACCURATE KEPLER EQUATION SOLVER WITHOUT TRANSCENDENTAL FUNCTION EVALUATIONS AAS 2-62 ACCURATE KEPLER EQUATION SOLVER WITHOUT TRANSCENDENTAL FUNCTION EVALUATIONS Adonis Pimienta-Peñalver and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 4260-4400

More information

CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION

CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION ARTIFICIAL SATELLITES, Vol. 50, No. 3 2015 DOI: 10.1515/arsa-2015-0008 CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION W. Góral, B. Skorupa AGH University

More information

ISIMA lectures on celestial mechanics. 1

ISIMA lectures on celestial mechanics. 1 ISIMA lectures on celestial mechanics. 1 Scott Tremaine, Institute for Advanced Study July 2014 The roots of solar system dynamics can be traced to two fundamental discoveries by Isaac Newton: first, that

More information

UNIVERSITY of CALIFORNIA SANTA CRUZ

UNIVERSITY of CALIFORNIA SANTA CRUZ UNIVERSITY of CALIFORNIA SANTA CRUZ NUMERICAL AND ANALYTIC RESEARCH INTO THE FORMATION OF THE HD 80606B PLANETARY SYSTEM A thesis submitted in partial satisfaction of the requirements for the degree of

More information

Chapter 2 Introduction to Binary Systems

Chapter 2 Introduction to Binary Systems Chapter 2 Introduction to Binary Systems In order to model stars, we must first have a knowledge of their physical properties. In this chapter, we describe how we know the stellar properties that stellar

More information

HAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS. 1. Introduction

HAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS. 1. Introduction HAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS ALESSANDRA CELLETTI 1 and LUIGI CHIERCHIA 2 1 Dipartimento di Matematica, Università di Roma Tor Vergata, Via della Ricerca Scientifica,

More information

Kozai-Lidov oscillations

Kozai-Lidov oscillations Kozai-Lidov oscillations Kozai (1962 - asteroids); Lidov (1962 - artificial satellites) arise most simply in restricted three-body problem (two massive bodies on a Kepler orbit + a test particle) e.g.,

More information

Targeting Mars. 19th June University of Ljubljana. Faculty of Mathematic and Physics. Seminar: 4. year-university program

Targeting Mars. 19th June University of Ljubljana. Faculty of Mathematic and Physics. Seminar: 4. year-university program University of Ljubljana Faculty of Mathematic and Physics Targeting Mars Seminar: 4. year-university program 19th June 016 Author: Jo²t Vrabi Koren Advisor: Prof. Dr. Tomaº Zwitter Ljubljana, June 016

More information

Astrodynamics (AERO0024)

Astrodynamics (AERO0024) Astrodynamics (AERO0024) L06: Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 Problem Statement? Hint #1: design the Earth-Mars transfer using known concepts

More information

Kepler, Newton, and laws of motion

Kepler, Newton, and laws of motion Kepler, Newton, and laws of motion First: A Little History Geocentric vs. heliocentric model for solar system (sec. 2.2-2.4)! The only history in this course is this progression: Aristotle (~350 BC) Ptolemy

More information

Chapter 9 Uniform Circular Motion

Chapter 9 Uniform Circular Motion 9.1 Introduction Chapter 9 Uniform Circular Motion Special cases often dominate our study of physics, and circular motion is certainly no exception. We see circular motion in many instances in the world;

More information

Joseph Castro Mentor: Nader Haghighipour

Joseph Castro Mentor: Nader Haghighipour ON THE POSSIBILITY OF ADDITIONAL PLANETS IN THE γ CEPHEI BINARY-PLANETARY SYSTEM Joseph Castro Mentor: Nader Haghighipour ABSTRACT Results of the simulations of the dynamical stability of additional hypothetical

More information

A REGION VOID OF IRREGULAR SATELLITES AROUND JUPITER

A REGION VOID OF IRREGULAR SATELLITES AROUND JUPITER The Astronomical Journal, 136:909 918, 2008 September c 2008. The American Astronomical Society. All rights reserved. Printed in the U.S.A. doi:10.1088/0004-6256/136/3/909 A REGION VOID OF IRREGULAR SATELLITES

More information

A SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION

A SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION 3 A SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION K.G. HADJIFOTINOU and JOHN D. HADJIDEMETRIOU Department of Physics, University of Thessaloniki, 540 06 Thessaloniki, Greece

More information

Use conserved quantities to reduce number of variables and the equation of motion (EOM)

Use conserved quantities to reduce number of variables and the equation of motion (EOM) Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.

More information

Simple algorithms for relative motion of satellites

Simple algorithms for relative motion of satellites Simple algorithms for relative motion of satellites Claudiu-Lucian Prioroc a, Seppo Mikkola a a Department of Physics and Astronomy, University of Turku, Tuorla Observatory, Väisäläntie 20, Piikkiö, FI

More information

Partial Derivatives of Observables with Respect to Two-Body Orbital Elements

Partial Derivatives of Observables with Respect to Two-Body Orbital Elements Partial Derivatives of Observables with Respect to Two-Body Orbital Elements Marc A. Murison U.S. Naval Observatory, Washington, DC murison@usno.vy.mil 6 November, 26 Abstract We calculate the partial

More information

The first term involves the cross product of two parallel vectors and so it vanishes. We then get

The first term involves the cross product of two parallel vectors and so it vanishes. We then get Physics 3550 Angular Momentum. Relevant Sections in Text: 3.4, 3.5 Angular Momentum You have certainly encountered angular momentum in a previous class. The importance of angular momentum lies principally

More information

(b) The period T and the angular frequency ω of uniform rotation are related to the cyclic frequency f as. , ω = 2πf =

(b) The period T and the angular frequency ω of uniform rotation are related to the cyclic frequency f as. , ω = 2πf = PHY 302 K. Solutions for problem set #9. Non-textbook problem #1: (a) Rotation frequency of 1 Hz means one revolution per second, or 60 revolutions per minute (RPM). The pre-lp vinyl disks rotated at 78

More information

THE PLANE OF THE KUIPER BELT

THE PLANE OF THE KUIPER BELT The Astronomical Journal, 127:2418 2423, 2004 April # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. THE PLANE OF THE KUIPER BELT Michael E. Brown Division of Geological

More information

Fundamentals of Astrodynamics and Applications

Fundamentals of Astrodynamics and Applications Fundamentals of Astrodynamics and Applications Third Edition David A. Vallado with technical contributions by Wayne D. McClain Space Technology Library Published Jointly by Microcosm Press Hawthorne, CA

More information

Astronomy 241: Review Questions #2 Distributed: November 7, 2013

Astronomy 241: Review Questions #2 Distributed: November 7, 2013 Astronomy 241: Review Questions #2 Distributed: November 7, 2013 Review the questions below, and be prepared to discuss them in class. For each question, list (a) the general topic, and (b) the key laws

More information

Lecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:

Lecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations

More information

Long-term integrations and stability of planetary orbits in our Solar system

Long-term integrations and stability of planetary orbits in our Solar system Mon. Not. R. Astron. Soc. 336, 483 500 (2002) Long-term integrations and stability of planetary orbits in our Solar system Takashi Ito and Kiyotaka Tanikawa Data Analysis Center and Theoretical Astrophysics

More information

Satellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN

Satellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,

More information

Dynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth.

Dynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth. Dynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth. P. Caspari ª ª Swinburne Astronomy Online, Student 6593658 Not yet received Not yet available online Abstract This paper announces

More information

HOW TO FIND SPATIAL PERIODIC ORBITS AROUND THE MOON IN THE TBP *

HOW TO FIND SPATIAL PERIODIC ORBITS AROUND THE MOON IN THE TBP * IJST, Transactions of Mechanical Engineering, Vol. 6, No. M1, pp 8-9 Printed in The Islamic Republic of Iran, 01 Shiraz University HOW TO FIND SPATIAL PERIODIC ORBITS AROUND THE MOON IN THE TBP * A. ARAM

More information

arxiv:astro-ph/ v2 31 Jan 2007

arxiv:astro-ph/ v2 31 Jan 2007 THE PRECESSION OF MERCURY AND THE DEFLECTION OF STARLIGHT FROM SPECIAL RELATIVITY ALONE Robert L. Kurucz Harvard-Smithsonian Center for Astrophysics arxiv:astro-ph/0608434v2 31 Jan 2007 60 Garden St, Cambridge,

More information

as the orbits of distant planetoids are expected to be randomized over billions of year by the gravity of the four giant planets.

as the orbits of distant planetoids are expected to be randomized over billions of year by the gravity of the four giant planets. Dynamics of Extreme Outer Solar System Objects Principle investigators/supervisors: Denis Gaidashev (Matematiska Institution); Erik Zackrisson (Institutionen fo r fysik och astronomi) 1. Purpose and aims

More information

5.1. Accelerated Coordinate Systems:

5.1. Accelerated Coordinate Systems: 5.1. Accelerated Coordinate Systems: Recall: Uniformly moving reference frames (e.g. those considered at 'rest' or moving with constant velocity in a straight line) are called inertial reference frames.

More information

EFFECTS OF DIRECT AND INDIRECT SOLAR RADIATION PRESSURE IN ORBITAL PARAMETERS OF GPS SATELITTES

EFFECTS OF DIRECT AND INDIRECT SOLAR RADIATION PRESSURE IN ORBITAL PARAMETERS OF GPS SATELITTES DOI: 10.2478/auom-2014-0039 An. Şt. Univ. Ovidius Constanţa Vol. 22(2),2014, 141 150 EFFECTS OF DIRECT AND INDIRECT SOLAR RADIATION PRESSURE IN ORBITAL PARAMETERS OF GPS SATELITTES Sergiu Lupu and Eugen

More information

Astrodynamics (AERO0024)

Astrodynamics (AERO0024) Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 6. Interplanetary Trajectories 6.1 Patched conic method 6.2 Lambert s problem

More information

Celestial Mechanics I. Introduction Kepler s Laws

Celestial Mechanics I. Introduction Kepler s Laws Celestial Mechanics I Introduction Kepler s Laws Goals of the Course The student will be able to provide a detailed account of fundamental celestial mechanics The student will learn to perform detailed

More information

NUMERICAL INTEGRATION OF THE RESTRICTED THREE BODY PROBLEM WITH LIE SERIES

NUMERICAL INTEGRATION OF THE RESTRICTED THREE BODY PROBLEM WITH LIE SERIES NUMERICAL INTEGRATION OF THE RESTRICTED THREE BODY PROBLEM WITH LIE SERIES ELBAZ. I. ABOUELMAGD 1, JUAN L.G. GUIRAO 2 AND A. MOSTAFA 3 Abstract. The aim of this work is to present some recurrence formulas

More information

Tests of General Relativity from observations of planets and spacecraft

Tests of General Relativity from observations of planets and spacecraft Tests of General Relativity from observations of planets and spacecraft Pitjeva E.V. Institute of Applied Astronomy, Russian Academy of Sciences Kutuzov Quay 10, St.Petersburg, 191187 Russia e-mail: evp@ipa.nw.ru

More information

A Classical Approach to the Stark-Effect. Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University

A Classical Approach to the Stark-Effect. Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University A Classical Approach to the Stark-Effect Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University Abstract The state of an atom in the presence of an external electric field is known

More information

A method for calculating probability of collision between space objects

A method for calculating probability of collision between space objects RAA 2014 Vol. 14 No. 5, 601 609 doi: 10.1088/1674 4527/14/5/009 http://www.raa-journal.org http://www.iop.org/journals/raa Research in Astronomy and Astrophysics A method for calculating probability of

More information