EXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS FOR EACH OF MULTIPLE BODIES
|
|
- Gyles Osborn Carter
- 6 years ago
- Views:
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
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 informationEFFICIENT 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 informationOrbit 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 informationintegral 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 informationPeriodic 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 informationDynamical 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 informationThe 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 informationExtending 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 informationLecture 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 informationA 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 informationA 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 informationRestricted 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 informationCalculation 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 informationarxiv: 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 informationAST111 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 information1 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 informationMOTIONS 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 informationSymplectic 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 informationAnalysis 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 informationThe 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 informationThe 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 informationThe 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 informationThe 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 informationReducing 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 informationThe 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 informationSTABILITY 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 informationTowards 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 informationMotion 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 informationAST111, 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 informationDynamics 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 informationLong 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 informationLecture 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 informationPlanetary 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 informationTides 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 informationObservational 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 informationResearch 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 informationorbits 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 informationACCURATE 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 informationASTRONOMY 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 informationInvestigation 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 informationarxiv: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 informationKEPLER 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 informationEccentricity 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 informationImportance 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 informationHill'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 informationTheory 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 informationOrbital 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 informationPhysics 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 informationThree 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 informationF. 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 informationOn 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 informationSatellite 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 informationGravitational 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 informationTHE 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 informationKeplerian 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 informationDIFFUSION 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 informationSatellite 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 informationPHYSICS. 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 informationAS3010: 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 informationTRANSFER 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 informationANNEX 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 informationResonance 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 informationAdvanced 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 informationACCURATE 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 informationCALCULATION 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 informationISIMA 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 informationUNIVERSITY 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 informationChapter 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 informationHAMILTONIAN 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 informationKozai-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 informationTargeting 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 informationAstrodynamics (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 informationKepler, 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 informationChapter 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 informationJoseph 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 informationA 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 informationA 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 informationUse 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 informationSimple 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 informationPartial 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 informationThe 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 =
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 informationTHE 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 informationFundamentals 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 informationAstronomy 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 informationLecture 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 informationLong-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 informationSatellite 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 informationDynamical 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 informationHOW 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 informationarxiv: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 informationas 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 information5.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 informationEFFECTS 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 informationAstrodynamics (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 informationCelestial 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 informationNUMERICAL 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 informationTests 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 informationA 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 informationA 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