# Stability of the Lagrange Points, L 4 and L 5

Save this PDF as:

Size: px
Start display at page:

## Transcription

1 Stability of the Lagrange Points, L 4 and L 5 Thomas Greenspan January 7, 014 Abstract A proof of the stability of the non collinear Lagrange Points, L 4 and L 5. We will start by covering the basics of stability, stating a theorem without proof) with a few corollaries, and then turn to the Lagrange points, proving first the stability of all Lagrange points in the z-direction and then restricting our attention to the the points, L 4 and L 5. Introduction The three body problem is one that has been studied for many centuries. It consists of considering bodies, subject to only to mutually attracting forces determined by the inverse square force, gravitation) and to solve for the motions of the three bodies. This task can be extremely complicated and, following Poincaré s example, most mathematicians and physicists have interested themselves in the finding periodic solutions to the problem, that is, where xt) = xt + T ) where T is the period. Indeed not only does this offer a useful framework from within which to work, but also describes states that are most applicable and useful to real world situations. Finding these periodic solutions is by no means trivial. It is a problem that has been worked on for centuries and still is in modern times, one of the latest publications dealing with it as recent as March of Among the first to interest themselves in and find solutions to the body problem were Leonhard Euler in 1765 and Joseph Lagrange in 177. In his publication, Essai sur le problème des trois corps Essay on the -Body Problem), Lagrange proposed a method that had never been used until then, that of considering only the distances between the three bodies rather than their absolute positions. Through this method, he found that there are exactly 5 different configurations the three bodies can be arranged in so that their movement is both circular and periodic. Given the initial position of two of the masses usually the largest masses), the 5 different locations that the third body can be in such that the solution is circular and periodic are called the 5 Lagrange or Libration points. Lagrange himself did not actually believe that instances of the 5 Libration points existed in the real world, as he states at the beginning of his derivation. Indeed, it was not until 1906 that the first example, an asteroid sharing Jupiter s orbit but ahead of it by close to 60 with respect to the sun, was discovered by professor Max Wolf at Heidelberg 8. Because the asteroid and those that were subsequently discovered at the L 4 and L 5 points of the Jupter-sun system were named after characters in the Illiad, asteroids in L 4 and L 5 of any system are generally referred to as Trojan asteroids. 1

2 Figure 1: Lagrange Points The three first points, L 1, L, L are all aligned with the earth-sun vector. L 4 and L 5 form an equilateral triangle with the sun and earth. The arrows indicate potentials and show the potential peaks, wells and saddles. To find these points, Lagrange needed to make some approximations indeed, the system is known today to be chaotic otherwise). The results nevertheless are often quite good at accurately describing dynamics of our solar system. In our case we consider the circular restricted -body problem as described below. The circular restricted -body problem The assumptions are as follows: Two of the masses, m 1 and m, are much heavier than the third one which we thus consider to be negligible. The center of mass is thus on a line between m 1 and m. The masses follow a circular orbit around the center of mass. We will describe our system by putting it in the form: d Xt) = AXt) dt 1 Stability and its consequences In any dynamical system with one or more equilibrium points, it is important to know whether the equilibrium points are stable, that is, whether a point at equilibrium can be subjected to small perturbations and still stay close to the equilibrium point. We use the following definition: Definition An equilibrium solution, x 0, to a dynamical system is considered: stable if for every small ε > 0, there exists a δ > 0 such that every solution, xt), with initial conditions xt 0 ) x 0 < δ is such that xt) x 0 < ε t t 0. asymptotically stable if, in addition to being stable, there exists δ 0 > 0 such that every solution, xt), with initial conditions xt 0 ) x 0 < δ 0 is such that xt) x 0 as t. In words, an equilibrium point is stable if, given a small distance ε, there always exists a distance, δ, such that any solution with initial conditions within δ of the equilibrium point will always stay within ε of the equilibrium for any time t. As we will see, Lagrange Points do not fulfill the much stronger condition of asymptotically stability but they are stable in some cases.

3 Beyond characterizing the dynamical system in question, considering stability is extremely important for any real-world applications. Indeed it is virtually impossible for any object to be precisely at an actual equilibrium point. Thus it is crucial to know what an object will do when close to an equilibrium for these are the actual dynamics that the object will necessarily follow. Before starting, we need a theorem and a few corollaries. Theorem 1.1 Given a system X t) = AXt) where A has distinct paired complex eigenvalues, α 1 +iβ 1, α 1 iβ 1..., α k +iβ k, α k iβ k. Let T be the matrix such that T 1 AT is in canonical form: where T 1 AT = B i = B 1... ) αi β i β i α i Then the general solution of X t) = AXt) is T Y t) where a 1 e α1t cos β 1 t + b 1 e α1t sin β 1 t a 1 e α1t sin β 1 t + b 1 e α1t cos β 1 t Y t) =. a k e αkt cos β k t + b k e αkt sin β k t a k e αkt sin β k t + b k e αkt cos β k t The proof of this theorem is beyond the scope of this paper, but can be found in chp. 6 of Differential Equations, Dynamical Systems, and an Introduction to Chaos 4. Corollary 1. A point in a dynamical system whose matrix of equations of motion, A, has an eigenvalue with positive real part is unstable. Proof Let m be the index s.t. α m > 0. To show the point is stable, it suffices to show that the point is unstable in only one direction. To do so we let a i = b i = 0 i m and assume that either a m 0 or b m 0 or both). This is effectively considering small perturbations of the equilibrium point in only one direction. We denote c i,j = T ) i,j where T is the matrix described in theorem 1.1). We begin by noting that T is invertible so every column must have at least one non-zero entry and no column is the multiple of another column. This means that, for every 1 j k, there exists an i such that c i,j 1 0 or c i,j 0 and c i,j 1 ±c i,j if this second condition were not the case then we would have T ) j 1 = ±T ) j which is impossible). Suppose the point is stable. Then we have T Y 0) T Y t) < ε for some ε and all t > 0. For any row, n, we get: ε > a mc n,m 1 c n,m ) e αmt c n,m 1 a m cos β m t + b m sin β m t) + c n,m a m cos β m t + b m sin β m t) 1.) There are two cases to consider: B k 1.1)

4 Case 1: b m 0. Let n be the index s.t. c n,m 1 0 or c n,m 0 and c n,m 1 c n,m. For simplicity we consider only the values t = kπ + π/)/β i for k Z. Our equation then becomes: ε > e αmt b m cn,m 1 + c n,m + am c n,m 1 c n,m ) As t, e αmt b m cn,m 1 + c n,m dominates and the value goes to ± depending on the sign of b m cn,m 1 + c n,m. Thus the solution is clearly not stable. Case : b m = 0. Similarly to the first case, let n be the index s.t. c n,m 1 0 or c n,m 0 and c n,m 1 c n,m. Again for simplicity, consider only the values t = kπ)/β i for k Z. Our equation then becomes: ε > e αmt a m cn,m 1 c n,m + am c n,m 1 c n,m ) As t, e αmt a m cn,m 1 c n,m dominates and the value goes to ± depending on the sign of a m cn,m 1 c n,m. Thus the solution is clearly not stable. Thus in either case, the point is not stable.. Corollary 1. A point in a dynamical system whose matrix of equations of motion, A, has purely imaginary non-zero) eigenvalues is stable. Proof Suppose that α i = 0 i in 1.1). As seen in 1.) the components of the solution becomes: k c n,i 1 a i cos β i t + b i sin β i t) + c n,i a i cos β i t + b i sin β i t) i=1 Since time, t, only appears in sines and cosines, it is clear that these values are completely bounded. Thus by choosing a 1,..., a k and b 1,..., b k appropriately i.e. small enough), we can make the solutions stay within ε of the initial state for any ε > 0 and we thus conclude that the points are stable.. The Lagrange Points Recall from lecture that the first step we used to finding the Lagrange Points is to consider the bodies in a rotating reference frame such that the two heavier masses m 1 and m do not move. Furthermore, since the movement of the three planets is planar, we can define the plane they move in to be the x-y plane with both heavier masses on the x-axis with the origin at the center of mass). Let R be the distance between m 1 and m, we then have the positions of the two heavier masses as: r 1 = m ) ) R m1 R, 0, 0 r =, 0, 0.1) m 1 + m m 1 + m Proposition.1 The angular frequency of the rotating reference frame, Ω, is given by: Ω R = GM 1 + M ).) This only holds because we are considering the circular restricted -body problem as discussed above and follows directly from Kepler s third law for more information see Landau & Livshits p. 1). 4

5 Solving the results found in lecture, We find our Lagrange points, L 4 and L 5 : ) ) ) ) R m1 m L 4 =, m 1 + m R, 0 R m1 m, L 5 =, m 1 + m R, 0.) To find out the dynamical stability of motion near the equilibrium points, we need to look at the generalized potential about the mass m. Because we are in a rotating frame we must add both the coriolis acceleration and the centrifugal acceleration. Let r be the position vector of m and Ω the angular velocity such that Ω = 0, 0, Ω). Furthermore let d 1 and d be the distances between m and m 1 and m respectively: d 1 = x + m R m 1 + m ) + y + z, d = x m 1R m 1 + m ) + y + z.4) Proposition. The total acceleration, r, and the generalized potential, U are given by r = Gm 1r r 1 ) d 1 Gm r r ) d Ω ṙ Ω Ω r.5) U = Gm 1 Gm Ωxẏ yẋ) Ω d 1 d x + y ).6) The second to last and last elements correspond to the coriolis and centrifugal accelerations respectively for further information on these forces see Landau & Livshits p. 18 1). Note that, in addition to position, the potential is dependent on velocity and, although it has no effect on the position of the Lagrange points, it must be taken into account when looking at their stability. To do so we separate the components dependent on velocity: ẍ = U = U + Ωxẏ yẋ) = Gm 1 Gm Ω d 1 d x + y ).7) Reducing.5) by components we get: Gm 1 x + ÿ = Gm 1y d 1 z = Gm 1z d 1 d 1 ) mr m 1+m Gm y d Gm z d Gm x ) m1r m 1+m d + Ωẏ + Ω x = U + Ωẏ.8) x Ωẋ + Ω y = U Ωẋ.9) y = U z.10) Finally, it can be useful to rewrite the generalized potential as a sum of partial derivatives. Using the Taylor series expansion of U around the Lagrange point x 0, y 0, z 0 ) we get to second order): U = U 0+U xx x 0 ) + U yy y 0 ) + U zz z 0 ) + 1 U xx x x 0 ) + U yyy y 0 ) + U zzz z 0 ) +U xyx x 0 )y y 0 ) + U xzx x 0 )z z 0 ) + U yzy y 0 )z z 0 ) 5

6 where U 0 = U and U X x0,y = U X for any variable X. Note, however, that U xz = 0,z 0) x0,y 0,z 0) U yz = 0 and furthermore we have that by definition, U x = U y = U z = 0 for Lagrange points. We thus get: U = U U xx x x 0 ) + U yyy y 0 ) + U zzz z 0 ) + U xyx x 0 )y y 0 ).11) Linearization and stability in the z-direction To analyze the stability about the equilibrium points in a system such as the -body problem, we linearize the equations of motion and look at the small perturbations. Note that because of how we have set up the reference frame, the Lagrange points are fixed points. We thus have the following equations: Plugging in these values into.11) then gives us: U = U Using.8)-.10) and.4) we get our equations: x = x 0 + δx ẋ = δẋ.1) y = y 0 + δy ẏ = δẏ.) z = z 0 + δz ż = δż.) U xx δx) + U yyδy) + U zzδz) + U xyδxδy.4) δẍ = U xxδx U xyδy + Ωδẏ.5) δÿ = U yyδy U xyδx Ωδẋ.6) δ z = U zzδz.7) We thus get the following: x δx δx y d z dt ẋ = d δy δy δz dt δẋ = δz U xx U xy 0 0 Ω 0 δẋ ẏ δẏ U xy U yy 0 Ω 0 0 δẏ ż δż 0 0 U zz δż.8) We start by considering only the z-direction. Theorem.1 All Lagrange points are stable in the z-direction. Proof Let δx = δy = 0. By.5)-.7) we can see that δz and δż are independent of δx, δẋ, δy, and δẏ and vice-versa so we can reduce our matrix to ) ) ) d δz 0 1 δz = dt δż U zz 0 δż 6

7 By.10) we have that U zz = Gm 1 d 1 + Gm d Furthermore, since d 1, d are distances we have d 1, d > 0 and thus U zz > 0. The eigenvalues of the above matrix are:.9) ±i U zz.10) Since these are always imaginary, we have by corollary 1. that the point is stable and thus conclude that all Lagrange points are stable in the z-direction.. In the same way and for the same reason) that we were able to consider only the z-direction, we can do the opposite and consider only the x and y directions. Our equation becomes: δx δx d δy dt δẋ = δy U xx U xy 0 Ω δẋ.11) δẏ U xy U yy Ω 0 δẏ It is this equation that we consider for the rest of the paper. 4 Stability of L 4 and L 5 The discovery of asteroids or other astrological bodies at the L 4 and L 5 points of almost every planet-sun system seems to indicate quite strongly that these points are indeed very stable. However, this is not always the case and requires that the heaviest mass, m 1 be significantly heavier than the second heaviest mass, m. We have the following theorem: Theorem 4.1 The Lagrange points, L 4 and L 5 are stable in all directions if and only if m ) m Proof We start by using.7) to compute the partial derivatives found in.11). Note that for L 4 and L 5, d 1 = d = R. Evaluating our partial double derivatives at L 4 and L 5 gives: ) ) U xx = Gm 1 d + Gm Gm 1 x + mr m 1+m Gm x m1r m 1+m 1 d d 5 1 d 5 Ω = Gm 1 + m ) R = 1 Gm 1 + m ) 4 R Ω Gm 1 m1+m )R m 1+m ) ) ) + m1+m )R m m 1+m ) R 5 Ω x= R ) m1 m m 1 +m 7

8 U yy = ) Gm 1 d + Gm 1 d Gm 1y d 5 Gm y 1 d 5 Ω y=± = Gm 1 + m ) R Gm 1 4 R + Gm 4 R R 5 Ω = 5 Gm 1 + m ) 4 R Ω U xy = Gm 1 x + d 5 1 mr m 1+m ) y ) R + Gm = Gm R 1 R 5 = Gm 1 m ) 4R = 4 κ ± Gm x ) R R d 5 Gm 1 + m ) R R ) m1r m 1+m y x= R ) m1 m m 1 +m,y=± R where κ ± = ±m 1 m )/m 1 + m ). Using.) we get: Our matrix from.11) thus becomes U xx = 4 Ω U yy = 9 4 Ω U xy = 4 κ ±Ω 4.) Ω 4 κ ±Ω 0 Ω 4.) 4 κ ±Ω 9 4 Ω Ω 0 The matrix has the four following eigenvalues: λ ± = ±i Ω 7κ ± σ ± = ±i Ω + 7κ ± 4.4) By corollary 1. and corollary 1., we have that all the eigenvalues must all be imaginary otherwise there will necessarily be at least one eigenvalue with Reλ) > 0, making the point unstable) so 7κ ± must be completely real. Note that since κ ± 1, 7κ ± so the only remaining condition is: After a little algebraic manipulation this gives us m m 1 7κ ± 0 4.5) 7 7 = which gives us the expected condition for the stability of the Lagrange points, L 4 and L ) 8

9 5 Conclusion This result would in fact be completely unexpected and surprising were it not for the discovery of actual instances of astrological bodies at these points in our Solar system. Indeed, by looking at the second partial derivatives in 4.) we find that U xx = U xx = 4 Ω < 0 and U yy = U yy = 9 4 Ω < 0. This would indicate that L 4 and L 5 are at peaks local maxima) of the potential in the x-y plane and would thus imply that that these points are extremely unstable. What gives L 4 and L 5 their stability is simply the coriolis force discussed previously. Initially, a body at L 4 and L 5 start moving away from the equilibrium point but, as the body pick up speed, the coriolis force takes effect sending the body into an effective orbit around the Lagrange point. Because of this effect, the areas around L 4 and L 5 that are effectively stable are in fact quite large as is illustrated by Figure 1. Indeed, this is the reason that so many trojan asteroids exist, some more than 5 off of 60, where the L 4 and L 5 are located. No satellites have been placed at these locations unlike L 1 and L hosting the SOHO satellite 7 and the WMAP satellite 6 respectively despite their inherent instability), however, they were visited in 009 by the STEREO satellites. As the only stable Lagrange points, L 4 and L 5 are unique phenomena in the solar system and the areas around them have been and are of great interest for example as possible places of origin of the moon or locations from which to better observe solar storms 5) to the astrophysics community. References 1 E. M. Lifshitz L. D. Landau. Course of Theoretical Physics, volume 1. Elsevier Inc., rd edition, 011. Joseph Lagrange. Le Problème des trois corps Richard Montgomery. A new solution to the three-body problem. Technical report, Notices Amer. Math Soc, Robert L. Devaney Morris W. Hirsch, Stephen Smale. Differential Equations, Dynamical Systems, and an Introduction to Chaos. Number Elsevier Inc., rd edition, NASA. Stereo. html. 6 NASA. Wilkinson microwave anisotropy probe. 7 NASA. Solar and heliospheric observatory. January Seth B. Nicholson. The trojan asteroids. Leaflet of the Astronomical Society of the Pacific, 8:9, Milovan & Dmitra sinović. V. Suvakov. Three classes of newtonian three-body planar periodic orbits. Phys. Rev. Lett., 110:11401, Mar 01. 9

### I ve Got a Three-Body Problem

I ve Got a Three-Body Problem Gareth E. Roberts Department of Mathematics and Computer Science College of the Holy Cross Mathematics Colloquium Fitchburg State College November 13, 2008 Roberts (Holy Cross)

### The Three Body Problem

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

### PURE Math Residents Program Gröbner Bases and Applications Week 3 Lectures

PURE Math Residents Program Gröbner Bases and Applications Week 3 Lectures John B. Little Department of Mathematics and Computer Science College of the Holy Cross June 2012 Overview of this week The research

### Abstract. M 1 r 1 +M 2 r 2 = 0.

Lagrangian Points Dennis Westra July 5, 07 Abstract In this text we shortly address some questions about the Lagrangian points. We show that there are five such points and then consider the stability of

### Lagrange Points and You. Robert Bond

Lagrange Points and You Robert Bond November 6, 2013 No one can make the argument that mankind has been merely a passive spectator on our now nearly 200,000 year old orbital journey around our sun and

### Three objects; 2+1 problem

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

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

### EE222 - Spring 16 - Lecture 2 Notes 1

EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued

### B5.6 Nonlinear Systems

B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre

### Lecture Notes for PHY 405 Classical Mechanics

Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter

### RELATIVE EQUILIBRIA IN THE

CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 10, Number 1, Spring 2003 RELATIVE EQUILIBRIA IN THE CHARGED n-body PROBLEM Based on an invited presentation at the annual meeting of the Canadian Applied

### On Central Configurations of the n Body Problem

of the n Body Problem Universidade Federal de Itajubá UNIFEI E-mail: lfmelo@unifei.edu.br Texas Christian University, Fort Worth GaGA Seminar August 23, 2017 Plan of the talk: 1 Introduction: Newtonian

### AB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT

FLÁIO SILESTRE DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT LECTURE NOTES LAGRANGIAN MECHANICS APPLIED TO RIGID-BODY DYNAMICS IMAGE CREDITS: BOEING FLÁIO SILESTRE Introduction Lagrangian Mechanics shall be

### Lecture 11: Binary stars

Lecture 11: Binary stars Senior Astrophysics 2018-04-11 Senior Astrophysics Lecture 11: Binary stars 2018-04-11 1 / 31 Outline 1 Observed characteristics 2 Interacting binaries 3 Gravity in a rotating

### PHYS The Three-Body Problem Special Lecture: Wednesday October 11, 2006, Juhan Frank, LSU. 1 The Three-Body Problem in Astronomy

PHYS 7 - The Three-Body Problem Special Lecture: Wednesday October,, Juhan Frank, LSU The Three-Body Problem in Astronomy The classical Newtonian three-body gravitational problem occurs in Nature exclusively

### Chaotic transport through the solar system

The Interplanetary Superhighway Chaotic transport through the solar system Richard Taylor rtaylor@tru.ca TRU Math Seminar, April 12, 2006 p. 1 The N -Body Problem N masses interact via mutual gravitational

### Connecting orbits and invariant manifolds in the spatial three-body problem

C C Dynamical A L T E C S H Connecting orbits and invariant manifolds in the spatial three-body problem Shane D. Ross Control and Dynamical Systems, Caltech Work with G. Gómez, W. Koon, M. Lo, J. Marsden,

### INTERPLANETARY TRANSFER TRAJECTORIES USING THE INVARIANT MANIFOLDS OF HALO ORBITS. A Thesis. presented to

INTERPLANETARY TRANSFER TRAJECTORIES USING THE INVARIANT MANIFOLDS OF HALO ORBITS A Thesis presented to the Faculty of California Polytechnic State University, San Luis Obispo In Partial Fulfillment of

### PHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010

Name: PHYS2330 Intermediate Mechanics Fall 2010 Final Exam Tuesday, 21 Dec 2010 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively

### Invariant Manifolds of Dynamical Systems and an application to Space Exploration

Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated

### Today in Astronomy 111: multiple-body systems

Today in Astronomy 111: multiple-body systems Two simple topics in systems of three bodies: Transfer orbits, and how asteroids get knocked from orbit to orbit I The two-body potential without Coriolis

### M2A2 Problem Sheet 3 - Hamiltonian Mechanics

MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,

### Periodic orbits high above the ecliptic plane in the solar sail 3-body problem

Periodic orbits high above the ecliptic plane in the solar sail 3-body problem Thomas Waters Department of Mechanical Engineering, University of Strathclyde, Glasgow. 4th October 2006 In conjunction with

### Uniform Circular Motion

Circular Motion Uniform Circular Motion Uniform Circular Motion Traveling with a constant speed in a circular path Even though the speed is constant, the acceleration is non-zero The acceleration responsible

### A plane autonomous system is a pair of simultaneous first-order differential equations,

Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

### Interplanetary Trajectory Design using Dynamical Systems Theory

Interplanetary Trajectory Design using Dynamical Systems Theory THESIS REPORT by Linda van der Ham 8 February 2012 The image on the front is an artist impression of the Interplanetary Superhighway [NASA,

### Stability of Lagrange Points: James Webb Space Telescope

Stability of Lagrange Points: James Webb Space Telescope Pearson, Kyle Alexander; Hanson, Jake; Cacolici, Gianna Nicole; Lejoly, Cassandra; Reynolds, Katherine ABSTRACT We present a detailed procedure

### 2x (x 2 + y 2 + 1) 2 2y. (x 2 + y 2 + 1) 4. 4xy. (1, 1)(x 1) + (1, 1)(y + 1) (1, 1)(x 1)(y + 1) 81 x y y + 7.

Homework 8 Solutions, November 007. (1 We calculate some derivatives: f x = f y = x (x + y + 1 y (x + y + 1 x = (x + y + 1 4x (x + y + 1 4 y = (x + y + 1 4y (x + y + 1 4 x y = 4xy (x + y + 1 4 Substituting

### Classical Mechanics. Luis Anchordoqui

1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body

### Welcome back to Physics 211. Physics 211 Spring 2014 Lecture Gravity

Welcome back to Physics 211 Today s agenda: Newtonian gravity Planetary orbits Gravitational Potential Energy Physics 211 Spring 2014 Lecture 14-1 1 Gravity Before 1687, large amount of data collected

### Welcome to CAMCOS Reports Day Fall 2010

Welcome to CAMCOS Reports Day Fall 2010 and Dynamics in Electrically Charged Binary Asteroid Systems Doug Mathews, Lara Mitchell, Jennifer Murguia, Tri Nguyen, Raquel Ortiz, Dave Richardson, Usha Watson,

### Physics 351 Wednesday, February 14, 2018

Physics 351 Wednesday, February 14, 2018 HW4 due Friday. For HW help, Bill is in DRL 3N6 Wed 4 7pm. Grace is in DRL 2C2 Thu 5:30 8:30pm. Respond at pollev.com/phys351 or text PHYS351 to 37607 once to join,

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

### 7 Planar systems of linear ODE

7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

### Regular Keplerian motions in classical many-body systems

Eur. J. Phys. 21 (2000) 1 18. Printed in the UK PII: S0143-0807(00)14581-7 Regular Keplerian motions in classical many-body systems Eugene I Butikov St Petersburg State University, St Petersburg, Russia

### Existence and stability of collinear equilibrium points in elliptic restricted three body problem with radiating primary and triaxial secondary

Modelling, Measurement and Control A Vol. 9, No., March, 08, pp. -8 Journal homepage: http://iieta.org/journals/mmc/mmc_a Existence and stability of collinear equilibrium points in elliptic restricted

### Periodic Orbits in Rotating Second Degree and Order Gravity Fields

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

### Welcome back to Physics 215

Welcome back to Physics 215 Today s agenda: More rolling without slipping Newtonian gravity Planetary orbits Gravitational Potential Energy Physics 215 Spring 2018 Lecture 13-1 1 Rolling without slipping

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

### The 3D restricted three-body problem under angular velocity variation. K. E. Papadakis

A&A 425, 11 1142 (2004) DOI: 10.1051/0004-661:20041216 c ESO 2004 Astronomy & Astrophysics The D restricted three-body problem under angular velocity variation K. E. Papadakis Department of Engineering

### The Planar, Circular, Restricted Four-Body Problem

The Planar, Circular, Restricted Four-Body Problem Gareth E. Roberts Julianne Kulevich Christopher J. Smith Department of Mathematics and Computer Science College of the Holy Cross NSF DMS-0708741 HC Faculty

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

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

### 5.1. Accelerated Coordinate Systems:

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

### Astronomy 111, Fall October 2011

Astronomy 111, Fall 011 4 October 011 Today in Astronomy 111: asteroids, perturbations and orbital resonances Leftovers: proofs of Kepler s second and third laws Elliptical orbits and center of mass More

### A NOTE ON THE DZIOBECK CENTRAL CONFIGURATIONS JAUME LLIBRE

This is a preprint of: A note on the Dziobek central configurations, Jaume Llibre, Proc. Amer. Math. Soc., vol. 43(8), 3587 359, 205. DOI: [http://dx.doi.org/0.090/s0002-9939-205-2502-6] A NOTE ON THE

### SPACECRAFT DYNAMICS NEAR A BINARY ASTEROID. F. Gabern, W.S. Koon and J.E. Marsden

PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 10 SPACECRAFT DYNAMICS NEAR A BINARY ASTEROID F. Gabern, W.S. Koon

### arxiv: v2 [physics.gen-ph] 16 Mar 2015

Dynamical description of Tychonian Universe arxiv:1304.7290v2 [physics.gen-ph] 16 Mar 2015 Luka Popov University of Zagreb, Department of Physics, Bijenička cesta 32, Zagreb, Croatia E-mail: lpopov@phy.hr

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

### Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 3, 2017

Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 3, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at positron.hep.upenn.edu/q351

### Bifurcations thresholds of halo orbits

10 th AstroNet-II Final Meeting, 15 th -19 th June 2015, Tossa Del Mar 1/23 spazio Bifurcations thresholds of halo orbits Dr. Ceccaroni Marta ceccaron@mat.uniroma2.it University of Roma Tor Vergata Work

### Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 6, 2015

Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 6, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at positron.hep.upenn.edu/q351

### Lecture 16. Gravitation

Lecture 16 Gravitation Today s Topics: The Gravitational Force Satellites in Circular Orbits Apparent Weightlessness lliptical Orbits and angular momentum Kepler s Laws of Orbital Motion Gravitational

### MSMS Basilio Bona DAUIN PoliTo

MSMS 214-215 Basilio Bona DAUIN PoliTo Problem 2 The planar system illustrated in Figure 1 consists of a bar B and a wheel W moving (no friction, no sliding) along the bar; the bar can rotate around an

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

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

### 2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

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

### Department of Mathematics IIT Guwahati

Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,

### Today. Laws of Motion. Conservation Laws. Gravity. tides

Today Laws of Motion Conservation Laws Gravity tides Newton s Laws of Motion Our goals for learning: Newton s three laws of motion Universal Gravity How did Newton change our view of the universe? He realized

### Part II. Dynamical Systems. Year

Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2

### A Thesis. Submitted to the Faculty. Purdue University. Kiarash Tajdaran. In Partial Fulfillment of the. Requirements for the Degree

i INCORPORATION OF MISSION DESIGN CONSTRAINTS IN FLOQUET MODE AND HAMILTONIAN STRUCTURE-PRESERVING ORBITAL MAINTENANCE STRATEGIES FOR LIBRATION POINT ORBITS A Thesis Submitted to the Faculty of Purdue

### Practice Problems for Final Exam

Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,

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

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

### Periodic orbits above the ecliptic in the solar sail. restricted 3-body problem

Periodic orbits above the ecliptic in the solar sail restricted 3-body problem Thomas J. Waters and Colin R. McInnes Department of Mechanical Engineering, University of Strathclyde, Glasgow, G1 1XJ, Scotland,

### Question 1: Spherical Pendulum

Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position

### Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325

Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0

### Physics 306, Honors Astronomy, Final Exam NAME

Physics 306, Honors Astronomy, Final Exam NAME You are graded on your work, with partial credit where it is deserved. 1. Clearly describe how you (as an astronomer) would measure each of the following:

### Nm kg. The magnitude of a gravitational field is known as the gravitational field strength, g. This is defined as the GM

Copyright FIST EDUCATION 011 0430 860 810 Nick Zhang Lecture 7 Gravity and satellites Newton's Law of Universal Gravitation Gravitation is a force of attraction that acts between any two masses. The gravitation

### Lecture Tutorial: Angular Momentum and Kepler s Second Law

2017 Eclipse: Research-Based Teaching Resources Lecture Tutorial: Angular Momentum and Kepler s Second Law Description: This guided inquiry paper-and-pencil activity helps students to describe angular

### Welcome back to Physics 215

Welcome back to Physics 215 Today s agenda: Gravity 15-2 1 Current assignments HW#15 due Monday, 12/12 Final Exam, Thursday, Dec. 15 th, 3-5pm in 104N. Two sheets of handwritten notes and a calculator

### Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.

Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics

### The Pennsylvania State University The Graduate School College of Engineering OPTIMAL EARTH RETURN TRANSFERS FROM LAGRANGE

The Pennsylvania State University The Graduate School College of Engineering OPTIMAL EARTH RETURN TRANSFERS FROM LAGRANGE POINT ORBITS USING PARTICLE SWARM OPTIMIZATION A Thesis in Aerospace Engineering

### Newton s Gravitational Law

1 Newton s Gravitational Law Gravity exists because bodies have masses. Newton s Gravitational Law states that the force of attraction between two point masses is directly proportional to the product of

### In this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:

Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this

### Potential/density pairs and Gauss s law

Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical

### The 350-Year Error in Solving the Newton - Kepler Equations

The 350-Year Error in Solving the Newton - Kepler Equations By Professor Joe Nahhas; Fall 1977 joenahhas1958@yahoo.com Abstract: Real time Physics solution of Newton's - Kepler's equations proves with

### NEW STACKED CENTRAL CONFIGURATIONS FOR THE PLANAR 6-BODY PROBLEM. Gokul Bhusal, Hamas Tahir, Dr. Xie. 8/03/2017 Summer Research 2017

NEW STACKED CENTRAL CONFIGURATIONS FOR THE PLANAR 6-BODY PROBLEM. Gokul Bhusal, Hamas Tahir, Dr. Xie. 8/03/2017 Summer Research 2017 Outline of the presentation: Celestial mechanics, Central configuration,

### Poincaré Map, Floquet Theory, and Stability of Periodic Orbits

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct

### Physics 141 Dynamics 1 Page 1. Dynamics 1

Physics 141 Dynamics 1 Page 1 Dynamics 1... from the same principles, I now demonstrate the frame of the System of the World.! Isaac Newton, Principia Reference frames When we say that a particle moves

### 3 Space curvilinear motion, motion in non-inertial frames

3 Space curvilinear motion, motion in non-inertial frames 3.1 In-class problem A rocket of initial mass m i is fired vertically up from earth and accelerates until its fuel is exhausted. The residual mass

### NOTES ON LINEAR ODES

NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary

### Gravity and the Orbits of Planets

Gravity and the Orbits of Planets 1. Gravity Galileo Newton Earth s Gravity Mass v. Weight Einstein and General Relativity Round and irregular shaped objects 2. Orbits and Kepler s Laws ESO Galileo, Gravity,

### A qualitative analysis of bifurcations to halo orbits

1/28 spazio A qualitative analysis of bifurcations to halo orbits Dr. Ceccaroni Marta ceccaron@mat.uniroma2.it University of Roma Tor Vergata Work in collaboration with S. Bucciarelli, A. Celletti, G.

### Hawking & the Universe

Hawking & the Universe This is a supplement to the lecture given on Sept.9 by Dr. Mounib El Eid, Physics department. It may motivate you to explore some of the issues presented here. There is no claim

### Lagrange Multipliers

Optimization with Constraints As long as algebra and geometry have been separated, their progress have been slow and their uses limited; but when these two sciences have been united, they have lent each

### SHAPE OF ORBITAL PATH

SHAPE OF ORBITAL PATH According to MATTER (Re-examined) Nainan K. Varghese, matterdoc@gmail.com http://www.matterdoc.info Abstract: In any system of macro bodies, relative considerations can provide only

### MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso

MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium

### Projectile Motion. Conceptual Physics 11 th Edition. Projectile Motion. Projectile Motion. Projectile Motion. This lecture will help you understand:

Conceptual Physics 11 th Edition Projectile motion is a combination of a horizontal component, and Chapter 10: PROJECTILE AND SATELLITE MOTION a vertical component. This lecture will help you understand:

### Theoretical physics. Deterministic chaos in classical physics. Martin Scholtz

Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton

### 11 Newton s Law of Universal Gravitation

Physics 1A, Fall 2003 E. Abers 11 Newton s Law of Universal Gravitation 11.1 The Inverse Square Law 11.1.1 The Moon and Kepler s Third Law Things fall down, not in some other direction, because that s

### Astrodynamics (AERO0024)

Astrodynamics (AERO0024) L04: Non-Keplerian Motion Gaëtan Kerschen Space Structures & Systems Lab (S3L) Non-Keplerian Motion 4 Dominant Perturbations Analytic Treatment Numerical Methods Concluding Remarks

### Gravitational radiation from a particle in bound orbit around a black hole; relativistic correction

Journal of Physics: Conference Series PAPER OPEN ACCESS Gravitational radiation from a particle in bound orbit around a black hole; relativistic correction To cite this article: Ashok Tiwari and Udayaraj

### Astro Lecture 12. Energy and Gravity (Cont d) 13/02/09 Habbal Astro Lecture 12 1

Astro 110-01 Lecture 12 Energy and Gravity (Cont d) 13/02/09 Habbal Astro110-01 Lecture 12 1 Energy due to movement of Kinetic Energy: object E k = ½ m v 2 13/02/09 Habbal Astro110-01 Lecture 12 2 Gravitational

### A space probe to Jupiter

Problem 3 Page 1 Problem 3 A space probe to Jupiter We consider in this problem a method frequently used to accelerate space probes in the desired direction. The space probe flies by a planet, and can

### BINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem

BINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem Michael Priolo Jerry Marsden, Ph.D., Mentor Shane Ross, Graduate Student, Co-Mentor Control and Dynamical Systems 107-81 Caltech,

### The Search for a Fundamental Theory of the Universe

The Search for a Fundamental Theory of the Universe Lecture 1- History & basic concepts, including Newton, Maxwell, Einstein & Quantum Mechanics Lecture 2 - Where are we now? General relativity & the Standard

### SOLUTIONS, PROBLEM SET 11

SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that

### Lecture PowerPoints. Chapter 6 Physics for Scientists and Engineers, with Modern Physics, 4 th edition Giancoli

Lecture PowerPoints Chapter 6 Physics for Scientists and Engineers, with Modern Physics, 4 th edition Giancoli 2009 Pearson Education, Inc. This work is protected by United States copyright laws and is