Research Article Variation of the Equator due to a Highly Inclined and Eccentric Disturber

Size: px
Start display at page:

Download "Research Article Variation of the Equator due to a Highly Inclined and Eccentric Disturber"

Transcription

1 Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 009, Article ID , 10 pages doi: /009/ Research Article Variation of the Equator due to a Highly Inclined and Eccentric Disturber Clair do Nascimento and adashi Yokoyama Universidade Estadual Paulista, IGCE-DEMAC Caixa Postal 178, CEP Rio Claro, SP, Brazil Correspondence should be addressed to adashi Yokoyama, tadashi@ms.rc.unesp.br Received 18 June 009; Accepted 19 October 009 Recommended by Dane Quinn he New Horizons project is currently in progress, but until recently years ago, a mission to Neptune-riton System predicted in the mid term was one of the priorities of NASA s Solar System Exploration theme. Any way, it is important to increase our knowledge on the dynamics of the inner satellites of Neptune, since according to some authors, this was a key question in the mission to Neptune-riton system. In a previous work, we presented the expansion of the disturbing function for the dynamics of this system. Here we derive the averaged classical equations of the precession of the equator for this problem. he highly inclined and retrograde orbit of riton makes this problem very unusual. herefore, the classical truncations in the inclinations are not acceptable, so that the precession equations must be obtained in closed form for the inclination. With a significant mass and due to its distance from Neptune, which is continuously decreasing, riton should exert in the future, important precession on Neptune s equator. he effects of this precession on the inner satellites are shown, including some resonant cases predicted in the future. Although riton s orbit is almost circular, no expansion in the eccentricity is needed. Copyright q 009 C. do Nascimento and. Yokoyama. his is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 1. Introduction Some years ago, the exploration of the Neptune system in the next ten years was posed as an important priority of NASA s Solar System Exploration theme 1. Although the New Horizons project is now the top priority for the far reaches of the solar system, most probably, to design a mission to explore the Neptune-riton system will be just a mater of time. he authors 1 raised some important questions about the dynamics stability of the arcs and inner satellites of Neptune. In order to give support to the mission, we must increase our knowledge on the dynamics of this system. In this work we derive the classical equations of the precession for the equator of Neptune. Basically, this equator is disturbed by the presence of riton since it has a remarkable mass: m /M N m, M N are

2 Mathematical Problems in Engineering z Z Ĩ Equator y Y X Ω χ x Reference plane Figure 1: XYZ:inertialsystem,xyz: moving system fixed on the equator. the masses of riton and Neptune, resp.. Moreover, riton s inclination with respect to Neptune s equator is quite unusual In a short time span, this precession should not be so important on the dynamics of the inner satellites. However, for a longer time more than 600 years, thiseffect must be investigated. his work is organized in the following way: in Section the averaged precessional equations for Neptune s equator considering arbitrary inclination obtained. Some authors hypothesize that in the past, the eccentricity of riton e, might have been very high. herefore, the precessional equations are also obtained for arbitrary e.insection 3, we discuss the definition of the invariable plane and the existence of a corresponding equilibrium solution when we assume circular orbit for riton. Considering the current dynamics of Neptune-riton system, several authors have improved the ephemerides and data of the small satellites e.g., Owen et al. 3. Here, some numerical simulations showing the effects of the precession, particularly for the future : 1 and 3 : 1 resonances involving riton and Proteus, are presented. Section 4 is devoted to the conclusions.. he Ĩ and Ω equations Let XYZ be an arbitrary reference system inertial and a moving system xyz fixed on the equator of Neptune see Figure 1. We define Ĩ inclination of the equator and Ω node of the equator 4. Due to the presence of a massive satellite riton, the equator of Neptune will be disturbed. Let be the associated torque due to riton and let W be the corresponding oblateness potential of Neptune: W GM Nm J R N r 3 1 3sin β,.1 where G is the gravitation constant, M N and m are the masses of Neptune and riton, respectively. he oblateness constant is J, β is the riton s latitude and r is the radial distance

3 Mathematical Problems in Engineering 3 Neptune-riton in the xyz system. he equatorial radius of Neptune is indicated by R N. hen: L, M, N r W,. where is the usual gradient operator with respect to riton coordinates and L, M, andn are the components of the torque, that is: W L y W M x W N x z z y W z W z W y y, x, x,.3 where x,y and z are the cartesian coordinates of riton with respect to the moving system xyz. On account of the rotational symmetry, we can show that N 0. o this end, since β is the latitude we have: sin β z /r. hen a direct calculation shows that W depends on riton coordinates only through r and z /r, therefore: W W r, z W x r W r, r x y z 1/, x r, W y W r y r..4 Substituting these last two relations in N given in.3, wegetn 0. Now, let w p N,q N,w N be the rotation vector of Neptune. According to Euler s dynamical equations: ṗ N C A A w Nq N L A, q N C A A w Np N M A,.5 ẇ N N C 0, where A B and C are the moments of inertia of the planet. From the geometric equations, w can be expressed in terms of the inertial XYZ system: p N sinĩ sin χ Ω cos χ Ĩ, q N sinĩ cos χ Ω sin χ Ĩ,.6 w N cosĩ Ω χ.

4 4 Mathematical Problems in Engineering aking the time derivative of p N, ṗ N cosĩ Ĩ sin χ Ω sinĩ χ cos χ Ω sin χ Ω sin χ χ Ĩ cos χ Ĩ..7 Assuming that the polar component of the rotation of the planet is dominant, the following quantities: Ĩ Ω, Ĩ e Ω e Ω can be neglected when compared to χ 5, 6. hen, ṗ N sinĩ χ cos χ Ω sin χ χ Ĩ..8 Substituting ṗ N e q N in the first Euler s equation, L A Ω χ cos χ Ĩ sin Ĩ χ sin χ C A w N Ω sinĩ cos χ sin χ Ĩ w N Ω sinĩ cos χ sin χ Ĩ..9 Substituting w N cos Ĩ Ω χ and neglecting terms factored by Ω and Ĩ Ω, we have C A w N sinĩ cos χ Ω sin χ Ĩ L A..10 Similarly, putting q N and p N into the second Euler dynamical equation, we obtain: C A w N sinĩ sin χ Ω cos χ Ĩ M A..11 Now we can eliminate L and M in the above equations. Defining W 0 Gm M N J R N /, we have W W 0 3z r 5 1, r 3 [ W 3z W 0 y r r 5 1 r 3 ] y r,.1 [ W 3z W 0 z r r 5 1 r 3 ] z r W z,

5 Mathematical Problems in Engineering 5 and so, L W y z W z y W z y,.13 where W z 6W 0 z /r 5. Similarly, M W zx.from.10 and.11, Ĩ Ω sin W z y cos χ x sin χ, Cw N Ĩ W z Cw N x cos χ y sin χ..14 If a sequence of three rotations, Ω, Ĩ, isχ are applied on the system XYZ, it will coincide with xyz. herefore we can write x, y,andz in terms of the well-known elliptic orbital elements of riton. x r cos f w cos Ω r sin f w sin Ω cos I, y r cos f w sin Ω r sin f w cos Ω cos I,.15 z r sin f w sin I, where r, f, w, Ω, w,andi refer to classical orbital elements of riton referred to X, Y, Z inertial system. Finally, averaging the resulting equations with respect to the mean anomaly of riton, the precessional equations, without any restriction on the inclination I or on the eccentricity e,are Ĩ Ω F 0 4a 3 1 e F 0 4a 3 1 e ] 3/ [sinĩ sin I sin Ω Ω cosĩ sin I sin Ω Ω, 3/ 1 3 cos I sin I cos I cos Ĩ cos Ω Ω cosĩ sinĩ sin I cosĩ cos Ω Ω,.16 where F 0 3Gm M N R N J Cw N..17 Notice that, as expected in this approximation, these two equations are decoupled from χ equation.

6 6 Mathematical Problems in Engineering R n 1 R n 1 Invariable plane t 1 S Ĩ Equator a riton s orbit t 1 π I b Figure : t 1 is perpendicular to the orbital plane of riton, n 1 is perpendicular to Neptune s equator, and R is the resulting vector which defines the invariable plane. 3. he Invariable Plane In this section, we define the reference plane. Let t 1 and n 1 be the orbital angular momentum of riton and rotation angular momentum of Netpune, respectively. hen R t 1 n 1 is a constant vector. herefore, the plane perpendicular to the direction of R is invariable see Figure, and it will be adopted as our reference plane. he dynamics of the system involving Ĩ, Ω, Ω,andI shows an interesting solution in this invariable plane. o see that, let us assume riton in a circular Keplerian orbit. In this case, we have t 1 m n a, where n is the mean motion, and m is the riton s mass. On the other hand, n 1 Cw N, where C is the polar moment of inertia of Neptune and w N is its angular velocity of rotation spin rate. Now,fromFigure, we have s sin π I, sin Ĩ scw m n a N 3.1 that is, the inclination Ĩ of the invariable plane is given by: sin Ĩ m n a sin I Cw N. 3. In order to have the complete set of equations that govern the motion of the equator, we have to consider.16 and also the corresponding equations for I and Ω. hey are di dt dω dt 1 R J n a sin I, Ω 1 n a sin I R J I, 3.3

7 Mathematical Problems in Engineering 7 where R J is the disturbing function due to the oblateness of Neptune which is written with respect to the following invariable plane: R J S 0 [ sin Ĩ 1 3 sin I 3 Ĩ 8 sin sin I cos Ω Ω 3 8 sin Ĩsin I cos Ω Ω ], 3.4 where GM N R N S 0 J [ a 3 1 e 3/ ]. 3.5 For Ω Ω π, we easily see that Ω S İ Ĩ 0, sinĩ sin Ĩ I, 3.6 Ω 3 GM N R N J 4 n a 5 sin I sin Ĩ I, where S 3Gm M N R N J Cw N a 3 herefore, Ω Ω 3GM NR N J 4a 3 sin 1 Ĩ I n a sin I m, 3.8 Cw N sinĩ that is, the term inside the brackets is automatically zero provided the reference plane is the invariable plane see 3.. In this case; Ω Ω 0, I and Ĩ are constants; Ω and Ω are synchronous. his means that the precession rate of the Neptune s equator is the same of the riton s node on the invariable plane 688 years. Itisalsoeasytoshowthat Ω Ω π defines a stable periodic solution of the system simple calculations show that such solution does not exist for Ω Ω 0. his means that it is very convenient to choose the invariable plane as the reference for our purposes.

8 8 Mathematical Problems in Engineering 1. 1 Inclination 0.6 Eccentricity Years Years a b Inclination 9 Inclination Years Years c d Figure 3: a, c : the lighter points were obtained taking into account the Ĩ and Ω equationsinthe integration, while in the thick points these contributions were neglected. b, d : time variation of the eccentricity b and inclination d for two orbits. Initially their paths are similar, but when their semimajor axes enter the : 1 resonance region, a jump appears in the eccentricity b, first in the orbit which contains the effect of the precession. he same occurs for the inclination d, that is, the precession accelerates the onset of the jump of the inclination. 4. Some Effects on the Inner Satellites Initially we tested our.16 for the case of Earth and Mars. As expected, the results agree quite well with very known values 6000 and years for Earth and Mars, resp.. For the Neptune s equator, considering a R N, our numerical experiments show 688 years for the period of precession, referred to the invariable plane Neptune-riton. his period is rather sensitive for the adopted value of the semimajor axis of riton. In this section we test the effect of the precession of Neptune s equator on an inner satellite of this planet. In Figure 3, we always integrate two orbits. he lighter points in the left panels, correspond to the orbits when we include the precessional effect, while the thick points come from the orbit

9 Mathematical Problems in Engineering 9 when we neglect.16. he inclinations are referred to the invariable plane. We consider a test particle Proteus disturbed by riton. In the integration of the equations of the motion, we used our averaged disturbing function R as described Yokoyama et al. 7 which is very efficient for this case. he initial conditions for Proteus are a 4.751, e , w 0 0, Ω 90 0, I 1 0 Figure 3 a and I 10 0 Figure 3 c,andλ 90 0, while for riton we consider a 14.35, e , w 0 0, Ω 70 0, I 157 0,andλ 70 0.heeffect of the precession is very clear: Figure 3 a, shows the variation of Proteus inclination which remains confined in a very narrow interval when precession equations are considered. In Figure 3 c, in opposition to the above panel, we take Proteus at high inclination hen, although the periods are different, the inclinations of both orbits undergo a large variation. Moreover, there is quite clear a periodic oscillation of about 670 years, which is related to the precession of the node Ω. Panels on the right, refer to the variation of the eccentricity Figure 3 b and inclination Figure 3 d of Proteus in the case of future : 1 resonance. As described by Yokoyama 8 and Yokoyama et al. 7, due to the tides, riton is falling down to the planet, so that, its orbital period will be in a 3 : 1 or : 1 orbital resonance with Proteus in the future. For these two panels, instead of the averaged equations, we integrated the exact equations and we applied a slow dissipation on the semimajor axis of riton in order to simulate the tidal effect. As pointed in Yokoyama et al. 7 the eccentricity suffers a strong variation and Proteus can be ejected. he effect of the precession is to anticipate the ejection. However, this is not a general rule, since some tests performed in the 3 : 1 resonance showed different results, that is, the precession can delay the ejection. More tests must be carried on in order to clarify these details. 5. Conclusion In this work, we derived the classical averaged equations for the precession of the Neptune s equator, extending for the case of arbitrary inclination and eccentricity of the disturber. Due to the high retrograde inclination of riton, these equations must be obtained in closed form for this inclination. If the invariable plane is chosen as the reference plane, the averaged equations of the precession in the case of e 0, satisfy an equilibrium solution where Ω and Ω remain synchronous while Ĩ and I are constants. he numerical experiments show that, when the precessional equations are considered, the inclination of an inner satellite of Neptune Proteus remains confined in a very narrow range, that is, the variation of the inclination with respect to the invariable plane is very small. However for high initial inclination 10 0 of the inner satellite, this confinement is not preserved anymore, that is, even considering the precessional effect, the amplitude of the variation of the inclination is significant. Additional simulations show that the precessional effects are also important in the future : 1 and 3 : 1 resonances, since the jump in inclination or eccentricity can be anticipated or delayed. Acknowledgments he authors thank FAPESP, CNPQ, and FUNDUNESP.

10 10 Mathematical Problems in Engineering References 1 H. B. Hammel, K. H. Baines, J. D. Cuzzi, et al., Exploration of the Netpune system, in he Future of Solar System Exploration, : Community Contributions to the NRC Solar System Exploration Decadal Survey, M. V. Sykes, Ed., vol. 7 of ASP Conference Series, pp , 00. P. Goldreich, N. Murray, P. Y. Longaretti, and D. Banfield, Neptune s story, Science, vol. 45, pp , W. M. Owen Jr., R. M. Vaughan, and S. P. Synnott, Orbits of the six new satellites of Neptune, Astronomical Journal, vol. 101, no. 4, pp , E. K. Woolard, Astronomical Papers of American Ephemeris, vol. 15, part 1, R. Kurth, Introduction to the Mechanics of the Solar System, Pergamon Press, London, UK, D. Lazzaro, Semi-Analytical heory for the Motion of Neptune s Satellites, Doctor thesis, Universidade de São Paulo, São Paulo, Brazil, Yokoyama, C. Nascimento, and M.. Santos, Inner satellites of Neptune I: the mutual disturbing function, Advances in Space Research, vol. 36, pp , Yokoyama, Possible effects of secular resonances in Phobos and riton, Planetary and Space Science, vol. 50, no. 1, pp , 00.

Searching for less perturbed elliptical orbits around Europa

Searching for less perturbed elliptical orbits around Europa Journal of Physics: Conference Series PAPER OPEN ACCESS Searching for less perturbed elliptical orbits around Europa To cite this article: J Cardoso dos Santos et al 2015 J. Phys.: Conf. Ser. 641 012011

More information

Research Article Dynamics of Artificial Satellites around Europa

Research Article Dynamics of Artificial Satellites around Europa Mathematical Problems in Engineering Volume 2013, Article ID 182079, 7 pages http://dx.doi.org/10.1155/2013/182079 Research Article Dynamics of Artificial Satellites around Europa Jean Paulo dos Santos

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

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

Astronomy 6570 Physics of the Planets. Precession: Free and Forced

Astronomy 6570 Physics of the Planets. Precession: Free and Forced Astronomy 6570 Physics of the Planets Precession: Free and Forced Planetary Precession We have seen above how information concerning the distribution of density within a planet (in particular, the polar

More information

Research Article Geodesic Effect Near an Elliptical Orbit

Research Article Geodesic Effect Near an Elliptical Orbit Applied Mathematics Volume 2012, Article ID 240459, 8 pages doi:10.1155/2012/240459 Research Article Geodesic Effect Near an Elliptical Orbit Alina-Daniela Vîlcu Department of Information Technology, Mathematics

More information

Third Body Perturbation

Third Body Perturbation Third Body Perturbation p. 1/30 Third Body Perturbation Modeling the Space Environment Manuel Ruiz Delgado European Masters in Aeronautics and Space E.T.S.I. Aeronáuticos Universidad Politécnica de Madrid

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

Chapter 5 - Part 1. Orbit Perturbations. D.Mortari - AERO-423

Chapter 5 - Part 1. Orbit Perturbations. D.Mortari - AERO-423 Chapter 5 - Part 1 Orbit Perturbations D.Mortari - AERO-43 Orbital Elements Orbit normal i North Orbit plane Equatorial plane ϕ P O ω Ω i Vernal equinox Ascending node D. Mortari - AERO-43 Introduction

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

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

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

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

Short-period planetary systems and their mysteries

Short-period planetary systems and their mysteries Short-period planetary systems and their mysteries Rosemary Mardling Monash Geneva 3 December 2014 Some open questions: gas giants How do hot jupiters arrive at their orbits? Are systems multiple systems

More information

Astronomy 6570 Physics of the Planets

Astronomy 6570 Physics of the Planets Astronomy 6570 Physics of the Planets Planetary Rotation, Figures, and Gravity Fields Topics to be covered: 1. Rotational distortion & oblateness 2. Gravity field of an oblate planet 3. Free & forced planetary

More information

CHANGING INCLINATION OF EARTH SATELLITES USING THE GRAVITY OF THE MOON

CHANGING INCLINATION OF EARTH SATELLITES USING THE GRAVITY OF THE MOON CHANGING INCLINATION OF EARTH SATELLITES USING THE GRAVITY OF THE MOON KARLA DE SOUZA TORRES AND A. F. B. A. PRADO Received 3 August 005; Revised 14 April 006; Accepted 18 April 006 We analyze the problem

More information

APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING

APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types

More information

13. Rigid Body Dynamics II

13. Rigid Body Dynamics II University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 13. Rigid Body Dynamics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License

More information

EART162: PLANETARY INTERIORS

EART162: PLANETARY INTERIORS EART162: PLANETARY INTERIORS Francis Nimmo Last Week Applications of fluid dynamics to geophysical problems Navier-Stokes equation describes fluid flow: Convection requires solving the coupled equations

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

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

Supplementary Materials for

Supplementary Materials for www.sciencemag.org/content/345/622/1317/suppl/dc1 Supplementary Materials for Chaotic dynamics of stellar spin in binaries and the production of misaligned hot Jupiters Natalia I. Storch, Kassandra R.

More information

Ph1a: Solution to the Final Exam Alejandro Jenkins, Fall 2004

Ph1a: Solution to the Final Exam Alejandro Jenkins, Fall 2004 Ph1a: Solution to the Final Exam Alejandro Jenkins, Fall 2004 Problem 1 (10 points) - The Delivery A crate of mass M, which contains an expensive piece of scientific equipment, is being delivered to Caltech.

More information

arxiv: v2 [astro-ph.ep] 28 Sep 2017

arxiv: v2 [astro-ph.ep] 28 Sep 2017 Noname manuscript No. (will be inserted by the editor) Tilting Styx and Nix but not Uranus with a Spin-Precession-Mean-motion resonance Alice C. Quillen 1, Yuan-Yuan Chen 1,2, Benoît Noyelles 3, & Santiago

More information

THIRD-BODY PERTURBATION USING A SINGLE AVERAGED MODEL

THIRD-BODY PERTURBATION USING A SINGLE AVERAGED MODEL INPE-1183-PRE/67 THIRD-BODY PERTURBATION USING A SINGLE AVERAGED MODEL Carlos Renato Huaura Solórzano Antonio Fernando Bertachini de Almeida Prado ADVANCES IN SPACE DYNAMICS : CELESTIAL MECHANICS AND ASTRONAUTICS,

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

HM NAUTICAL ALMANAC OFFICE NAO TECHNICAL N OTE

HM NAUTICAL ALMANAC OFFICE NAO TECHNICAL N OTE Council for the Central Laboratory of the Research Councils HM NAUTICAL ALMANAC OFFICE NAO TECHNICAL N OTE No. 73 2003 June Approximate ephemerides and phenomena of natural satellites by D.B. Taylor and

More information

Evolution of Mercury s obliquity

Evolution of Mercury s obliquity Icarus 181 2006) 327 337 www.elsevier.com/locate/icarus Evolution of Mercury s obliquity Marie Yseboodt, Jean-Luc Margot Department of Astronomy, Cornell University, Ithaca, NY 14853, USA Received 24 June

More information

Physics 106b/196b Problem Set 9 Due Jan 19, 2007

Physics 106b/196b Problem Set 9 Due Jan 19, 2007 Physics 06b/96b Problem Set 9 Due Jan 9, 2007 Version 3: January 8, 2007 This problem set focuses on dynamics in rotating coordinate systems (Section 5.2), with some additional early material on dynamics

More information

PRECESSIONS IN RELATIVITY

PRECESSIONS IN RELATIVITY PRECESSIONS IN RELATIVITY COSTANTINO SIGISMONDI University of Rome La Sapienza Physics dept. and ICRA, Piazzale A. Moro 5 00185 Rome, Italy. e-mail: sigismondi@icra.it From Mercury s perihelion precession

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

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

Chapter 11. Angular Momentum

Chapter 11. Angular Momentum Chapter 11 Angular Momentum Angular Momentum Angular momentum plays a key role in rotational dynamics. There is a principle of conservation of angular momentum. In analogy to the principle of conservation

More information

r( θ) = cos2 θ ω rotation rate θ g geographic latitude - - θ geocentric latitude - - Reference Earth Model - WGS84 (Copyright 2002, David T.

r( θ) = cos2 θ ω rotation rate θ g geographic latitude - - θ geocentric latitude - - Reference Earth Model - WGS84 (Copyright 2002, David T. 1 Reference Earth Model - WGS84 (Copyright 22, David T. Sandwell) ω spheroid c θ θ g a parameter description formula value/unit GM e (WGS84) 3.9864418 x 1 14 m 3 s 2 M e mass of earth - 5.98 x 1 24 kg

More information

Revision of the physical backgrounds of Earth s rotation

Revision of the physical backgrounds of Earth s rotation Revision of the physical backgrounds of Earth s rotation L. Völgyesi Budapest University of Technology and Economics, H-151 Budapest, Müegyetem rkp. 3. Hungary e-mail: lvolgyesi@epito.bme.hu bstract. Rotation

More information

Physics for Scientists and Engineers 4th Edition, 2017

Physics for Scientists and Engineers 4th Edition, 2017 A Correlation of Physics for Scientists and Engineers 4th Edition, 2017 To the AP Physics C: Mechanics Course Descriptions AP is a trademark registered and/or owned by the College Board, which was not

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

A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3:1 RESONANCE. Érica C. Nogueira, Othon C. Winter

A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3:1 RESONANCE. Érica C. Nogueira, Othon C. Winter A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3: RESONANCE Érica C. Nogueira, Othon C. Winter Grupo de Dinâmica Orbital e Planetologia UNESP -- Guaratinguetá -- Brazil Antonio F.B. de

More information

7. The gyroscope. 7.1 Introduction. 7.2 Theory. a) The gyroscope

7. The gyroscope. 7.1 Introduction. 7.2 Theory. a) The gyroscope K 7. The gyroscope 7.1 Introduction This experiment concerns a special type of motion of a gyroscope, called precession. From the angular frequency of the precession, the moment of inertia of the spinning

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

Dynamics of the Earth

Dynamics of the Earth Time Dynamics of the Earth Historically, a day is a time interval between successive upper transits of a given celestial reference point. upper transit the passage of a body across the celestial meridian

More information

OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES

OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES Gislaine de Felipe and Antonio Fernando Bertachini de Almeida Prado Instituto Nacional de Pesquisas Espaciais - São José dos Campos - SP - 12227-010

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

On Sun-Synchronous Orbits and Associated Constellations

On Sun-Synchronous Orbits and Associated Constellations On Sun-Synchronous Orbits and Associated Constellations Daniele Mortari, Matthew P. Wilkins, and Christian Bruccoleri Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843,

More information

arxiv: v1 [astro-ph.ep] 11 Aug 2016

arxiv: v1 [astro-ph.ep] 11 Aug 2016 Noname manuscript No. (will be inserted by the editor) Secular and tidal evolution of circumbinary systems. Alexandre C. M. Correia Gwenaël Boué Jacques Laskar arxiv:168.3484v1 [astro-ph.ep] 11 Aug 216

More information

Orbital Stability Regions for Hypothetical Natural Satellites

Orbital Stability Regions for Hypothetical Natural Satellites Orbital Stability Regions for Hypothetical Natural Satellites By Samantha RIEGER, 1) Daniel SCHEERES, 1) 1) Ann and H.J. Smead Department of Aerospace Engineering Sciences, University of Colorado Boulder,

More information

Today in Astronomy 111: rings, gaps and orbits

Today in Astronomy 111: rings, gaps and orbits Today in Astronomy 111: rings, gaps and orbits Gap sizes: the Hill radius Perturbations and resonances The variety of structures in planetary rings Spiral density waves Titan Bending waves Horseshoe and

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

An Introduction to Celestial Mechanics

An Introduction to Celestial Mechanics An Introduction to Celestial Mechanics This accessible text on classical celestial mechanics the principles governing the motions of bodies in the solar system provides a clear and concise treatment of

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

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

The in-plane Motion of a Geosynchronous Satellite under the Gravitational Attraction of the Sun, the Moon and the Oblate Earth

The in-plane Motion of a Geosynchronous Satellite under the Gravitational Attraction of the Sun, the Moon and the Oblate Earth J Astrophys Astr (1990) 11, 1 10 The in-plane Motion of a Geosynchronous Satellite under the Gravitational Attraction of the Sun, the Moon and the Oblate Earth Κ Β Bhatnagar Department of Mathematics,

More information

1/30. Rigid Body Rotations. Dave Frank

1/30. Rigid Body Rotations. Dave Frank . 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid

More information

Orbital resonances in the inner neptunian system I. The 2:1 Proteus Larissa mean-motion resonance

Orbital resonances in the inner neptunian system I. The 2:1 Proteus Larissa mean-motion resonance Icarus 188 (2007) 386 399 www.elsevier.com/locate/icarus Orbital resonances in the inner neptunian system I. The 2:1 Proteus arissa mean-motion resonance Ke Zhang, Douglas P. Hamilton Department of Astronomy,

More information

Beyond Janus & Epimetheus: Momentum Trading Among Co-Orbiting Satellite Groups

Beyond Janus & Epimetheus: Momentum Trading Among Co-Orbiting Satellite Groups Beyond Janus & Epimetheus: Momentum Trading Among Co-Orbiting Satellite Groups DOUG BALCOM U NIVERSITY OF WASHINGTON APPLIED MATHEMATICS Special Thanks to Sasha Malinsky Janus and Epimetheus: Momentum

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

2 so that the (time-dependent) torque, as measured in space-xed coordinates, is ~N = ma2! 2 (^x cos!t +^y sin!t) 10 = a p 5 2 ^z ~ F 2 ~F 1 = a p 5^z

2 so that the (time-dependent) torque, as measured in space-xed coordinates, is ~N = ma2! 2 (^x cos!t +^y sin!t) 10 = a p 5 2 ^z ~ F 2 ~F 1 = a p 5^z 1 PHYS 321 Homework Assignment #9 Due: Friday, 22 November 2002 (4 probs) 1. Problem 7-14, B&O p. 280. The situation is as shown above. We take the body-xed axes ^x 0 and ^y 0 to lie in the plane of the

More information

arxiv: v1 [astro-ph.ep] 12 Dec 2016

arxiv: v1 [astro-ph.ep] 12 Dec 2016 Astronomy & Astrophysics manuscript no. manuscript_arxiv c ESO 2018 December 20, 2018 Frequency-dependent tidal dissipation in a viscoelastic Saturnian core and expansion of Mimas semi-major axis arxiv:1612.03664v1

More information

Lecture 1d: Satellite Orbits

Lecture 1d: Satellite Orbits Lecture 1d: Satellite Orbits Outline 1. Newton s Laws of Motion 2. Newton s Law of Universal Gravitation 3. Kepler s Laws 4. Putting Newton and Kepler s Laws together and applying them to the Earth-satellite

More information

Can we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth?

Can we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth? Can we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth? BARRETO L., FORTIN M.-A., IREDALE A. In this simple analysis, we compare the historical record

More information

imin...

imin... Pulsar Timing For a detailed look at pulsar timing and other pulsar observing techniques, see the Handbook of Pulsar Astronomy by Duncan Lorimer and Michael Kramer. Pulsars are intrinsically interesting

More information

THE LUNAR LIBRATION: COMPARISONS BETWEEN VARIOUS MODELS- A MODEL FITTED TO LLR OBSERVATIONS

THE LUNAR LIBRATION: COMPARISONS BETWEEN VARIOUS MODELS- A MODEL FITTED TO LLR OBSERVATIONS THE LUNAR LIBRATION: COMPARISONS BETWEEN VARIOUS MODELS- A MODEL FITTED TO LLR OBSERVATIONS J. CHAPRONT, G. FRANCOU SYRTE - Observatoire de Paris - UMR 8630/CNRS 61, avenue de l Observatoire 75014 Paris

More information

expression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects.

expression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects. RUSSIAN JOURNAL OF EARTH SCIENCES, English Translation, VOL, NO, DECEMBER 998 Russian Edition: JULY 998 On the eects of the inertia ellipsoid triaxiality in the theory of nutation S. M. Molodensky Joint

More information

Experiment Design and Performance. G. Catastini TAS-I (BUOOS)

Experiment Design and Performance. G. Catastini TAS-I (BUOOS) Experiment Design and Performance G. Catastini TAS-I (BUOOS) 10 EP and the GRT Einstein s General Relativity Theory Weak Equivalence Principle: all test particles at the same space-time point in a given

More information

Ceres Rotation Solution under the Gravitational Torque of the Sun

Ceres Rotation Solution under the Gravitational Torque of the Sun Ceres Rotation Solution under the Gravitational Torque of the Sun Martin Lara, Toshio Fukushima, Sebastián Ferrer (*) Real Observatorio de la Armada, San Fernando, Spain ( ) National Astronomical Observatory,

More information

Lesson 9. Luis Anchordoqui. Physics 168. Tuesday, October 24, 17

Lesson 9. Luis Anchordoqui. Physics 168. Tuesday, October 24, 17 Lesson 9 Physics 168 1 Static Equilibrium 2 Conditions for Equilibrium An object with forces acting on it but that is not moving is said to be in equilibrium 3 Conditions for Equilibrium (cont d) First

More information

Rotational Kinematics and Dynamics. UCVTS AIT Physics

Rotational Kinematics and Dynamics. UCVTS AIT Physics Rotational Kinematics and Dynamics UCVTS AIT Physics Angular Position Axis of rotation is the center of the disc Choose a fixed reference line Point P is at a fixed distance r from the origin Angular Position,

More information

previously Department of Solid Mechanics, Royal Institute of Technology, KTH, Sweden

previously Department of Solid Mechanics, Royal Institute of Technology, KTH, Sweden Hypergravity derived from the Le-density and applied to a model for explosion. Analysis of universal constants of gravity and light in conjunction with the J-space. Lena J-T Strömberg previously Department

More information

The Restricted 3-Body Problem

The Restricted 3-Body Problem The Restricted 3-Body Problem John Bremseth and John Grasel 12/10/2010 Abstract Though the 3-body problem is difficult to solve, it can be modeled if one mass is so small that its effect on the other two

More information

Long-term perturbations due to a disturbing body in elliptic

Long-term perturbations due to a disturbing body in elliptic Long-term perturbations due to a disturbing body in elliptic inclined orbit Xiaodong Liu 1, Hexi Baoyin, and Xingrui Ma 3 School of Aerospace, Tsinghua University, 100084 Beijing, China Email: liu-xd08@mails.tsinghua.edu.cn;

More information

Homework 2; due on Thursday, October 4 PY 502, Computational Physics, Fall 2018

Homework 2; due on Thursday, October 4 PY 502, Computational Physics, Fall 2018 Homework 2; due on Thursday, October 4 PY 502, Computational Physics, Fall 2018 Department of Physics, Boston University Instructor: Anders Sandvik PERTURBATION OF A GEO-STATIONARY ORBIT Here you will

More information

arxiv: v1 [astro-ph.ep] 23 Nov 2011

arxiv: v1 [astro-ph.ep] 23 Nov 2011 Draft version September 3, 18 Preprint typeset using L A TEX style emulateapj v. 11/1/9 PUMPING THE ECCENTRICITY OF EXOPLANETS BY TIDAL EFFECT Alexandre C.M. Correia Department of Physics, I3N, University

More information

Introduction to CLASSICAL MECHANICS

Introduction to CLASSICAL MECHANICS Introduction to CLASSICAL MECHANICS Introduction to CLASSICAL MECHANICS A.P. FRENCH Massachusetts Institute oj Technology M.G. EBISON The Institute oj Physics, London KLUWER ACADEMIC PUBLISHERS DORDRECHT

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

DYNAMICS OF THE LUNAR SPIN AXIS

DYNAMICS OF THE LUNAR SPIN AXIS The Astronomical Journal, 131:1864 1871, 2006 March # 2006. The American Astronomical Society. All rights reserved. Printed in U.S.A. DYNAMICS OF THE LUNAR SPIN AXIS Jack Wisdom Massachusetts Institute

More information

Resonant Cosmos. Resume. planet's rotation around a central body (star) in the plane of celestial equator and

Resonant Cosmos. Resume. planet's rotation around a central body (star) in the plane of celestial equator and V.D. Krasnov Email: apeyron7@yandex.ru Resonant Cosmos (Oscillatory motion and forming of the inclination of the planets' rotation plane) Resume A mechanism of originating and forming of the inclination

More information

Earth-Centered, Earth-Fixed Coordinate System

Earth-Centered, Earth-Fixed Coordinate System Fundamentals of Global Positioning System Receivers: A Software Approach James Bao-Yen Tsui Copyright 2000 John Wiley & Sons, Inc. Print ISBN 0-471-38154-3 Electronic ISBN 0-471-20054-9 CHAPTER FOUR Earth-Centered,

More information

Lecture 2c: Satellite Orbits

Lecture 2c: Satellite Orbits Lecture 2c: Satellite Orbits Outline 1. Newton s Laws of Mo3on 2. Newton s Law of Universal Gravita3on 3. Kepler s Laws 4. Pu>ng Newton and Kepler s Laws together and applying them to the Earth-satellite

More information

Chapter 2: Orbits and Launching Methods

Chapter 2: Orbits and Launching Methods 9/20/ Chapter 2: Orbits and Launching Methods Prepared by Dr. Mohammed Taha El Astal EELE 6335 Telecom. System Part I: Satellite Communic ations Winter Content Kepler s First, Second, and Third Law Definitions

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

Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming

Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming IEPC-2015-299 /ISTS-2015-b-299 Presented at Joint Conference of 30th International Symposium on Space Technology

More information

THE ORBITAL MOTION IN LEO REGION: CBERS SATELLITES AND SPACE DEBRIS. Keywords: CBERS Satellites, Space Debris, Orbital Motion, Resonance.

THE ORBITAL MOTION IN LEO REGION: CBERS SATELLITES AND SPACE DEBRIS. Keywords: CBERS Satellites, Space Debris, Orbital Motion, Resonance. THE ORBITAL MOTION IN LEO REGION: CBERS SATELLITES AND SPACE DEBRIS Jarbas Cordeiro Sampaio (1a), Rodolpho Vilhena de Moraes (1b) and Sandro da Silva Fernandes (2c) (1) UNIFESP- Univ Federal de Sao Paulo,

More information

INTERNAL THALES ALENIA SPACE

INTERNAL THALES ALENIA SPACE Workshop Galileo Galilei (GG) and GGG lab prototype: state of the art and new possibilities Template reference : 100181670S-EN GG error budget from the end-to-end simulator developed at TAS-I Giuseppe

More information

Computational Problem: Keplerian Orbits

Computational Problem: Keplerian Orbits Computational Problem: Keplerian Orbits April 10, 2006 1 Part 1 1.1 Problem For the case of an infinite central mass and an orbiting test mass, integrate a circular orbit and an eccentric orbit. Carry

More information

is a revolution relative to a fixed celestial position. is the instant of transit of mean equinox relative to a fixed meridian position.

is a revolution relative to a fixed celestial position. is the instant of transit of mean equinox relative to a fixed meridian position. PERIODICITY FORMULAS: Sidereal Orbit Tropical Year Eclipse Year Anomalistic Year Sidereal Lunar Orbit Lunar Mean Daily Sidereal Motion Lunar Synodical Period Centenial General Precession Longitude (365.25636042

More information

Solar vs. Lunar Tides

Solar vs. Lunar Tides 1 2 3 4 Solar vs. Lunar Tides In the force equations M is the mass of the tide-causing object, r is the separation between the two objects. dr is the size of the object on which the tides are being raised.

More information

NAVIGATION & MISSION DESIGN BRANCH

NAVIGATION & MISSION DESIGN BRANCH c o d e 5 9 5 National Aeronautics and Space Administration Michael Mesarch Michael.A.Mesarch@nasa.gov NAVIGATION & MISSION DESIGN BRANCH www.nasa.gov Outline Orbital Elements Orbital Precession Differential

More information

Gravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the

Gravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the Gravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the Geocentric-Heliocentric debate (Brahe s data and Kepler s analysis)

More information

B.Sc. (Semester - 5) Subject: Physics Course: US05CPHY01 Classical Mechanics

B.Sc. (Semester - 5) Subject: Physics Course: US05CPHY01 Classical Mechanics 1 B.Sc. (Semester - 5) Subject: Physics Course: US05CPHY01 Classical Mechanics Question Bank UNIT: I Multiple choice questions: (1) The gravitational force between two masses is (a) Repulsive (b) Attractive

More information

Practice Test 3. Name: Date: ID: A. Multiple Choice Identify the choice that best completes the statement or answers the question.

Practice Test 3. Name: Date: ID: A. Multiple Choice Identify the choice that best completes the statement or answers the question. Name: Date: _ Practice Test 3 Multiple Choice Identify the choice that best completes the statement or answers the question. 1. A wheel rotates about a fixed axis with an initial angular velocity of 20

More information

Physical Dynamics (PHY-304)

Physical Dynamics (PHY-304) Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.

More information

Problem Solving Session 11 Three Dimensional Rotation and Gyroscopes Solutions

Problem Solving Session 11 Three Dimensional Rotation and Gyroscopes Solutions MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01 Problem Solving Session 11 Three Dimensional Rotation and Gyroscopes Solutions W14D3-1 Rotating Skew Rod Solution Consider a simple

More information

Precise modeling of the precession-nutation of Mars

Precise modeling of the precession-nutation of Mars Astron. Astrophys. 345, 8 97 (999) ASTRONOMY AND ASTROPHYSICS Precise modeling of the precession-nutation of Mars S. Bouquillon and J. Souchay Observatoire de Paris/DANOF, URA 5 du CNRS, 6 Avenue de l

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

Universal Gravitation

Universal Gravitation Universal Gravitation Newton s Law of Universal Gravitation Every particle in the Universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely

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

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

8.012 Physics I: Classical Mechanics Fall 2008

8.012 Physics I: Classical Mechanics Fall 2008 IT OpenCourseWare http://ocw.mit.edu 8.012 Physics I: Classical echanics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. ASSACHUSETTS INSTITUTE

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

Resonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009

Resonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009 Resonance Capture Alice Quillen University of Rochester Isaac Newton Institute Dec 2009 This Talk Resonance Capture extension to nonadiabatic regime Astrophysical settings: Dust spiraling inward via radiation

More information