Calculating L2 Lagrangian Point: A Beginner's Guide

[Irrelativity]

Hi this is my first time posting on this forum. I have an question about Lagragian points.

I was trying to find L2 lagrangian point, a point that lies on the line defined by the two large masses, beyond the smaller of the two. Here, the gravitational forces of the two large masses balance the centrifugal force on the smaller mass. Let's say the larger object has mass M and the smaller object has mass m, and the distance between the two objects is R and the distance from the smaller object to L2 is r, then how do you guys derive a formula to calculate r? I was given a hint that I need to use approximation method meaning R is exceedingly greater than r such that (r/R)^2 is equal to zero. Again, Unlike L1 which is located between the two objects, L2 lies behind the smaller object.

At least if you can give me some hints, I would really appreciate. Thanks.

 

[D H]

Work in the frame with origin at the center of mass and rotating with the objects so that the location of the two massive objects in this frame is constant. The five Lagrangian points will also have a fixed location in this frame. The accelerations on a fixed point in a rotating frame are the gravitational accelerations toward the two bodies and the centrifugal force. The accelerations must sum to zero at some point if it is to remain fixed. The first goal is to find this polynomial. Once you have this polynomial, you can solve it using numerical methods, or you can find an approximate solution.

 

[rbj]

DH, i just wanted to be a little more anal about this, to include a couple of points that i needed to remember when i did this (for recreational purposes, i don't think i ever needed to do this for a class).

Work in the frame with origin at the center of mass

of the two massive objects. ignoring any masses of anything that might be living at the L2 point, we can easily model the two planets as point masses revolving around their common center of mass. and we don't have to think about them being affected by the test mass at L2.

and rotating with the objects so that the location of the two massive objects in this frame is constant. The five Lagrangian points will also have a fixed location in this frame. The accelerations on a fixed point in a rotating frame are the gravitational accelerations toward the two bodies and the centrifugal force. The accelerations must sum to zero at some point if it is to remain fixed.

the other thing i needed, but had to remember to something i only used in freshman physics or in my Statics and Dynamics class required for engineering was D'Alembert's principle. this centrifugal force is a fictitious force that we put in because of this.

The first goal is to find this polynomial. Once you have this polynomial, you can solve it using numerical methods, or you can find an approximate solution.

i think you can get an closed form solution, no? the order of the polynomial doesn't exceed 3, does it?

i think the L4 and L5 points are the cool ones. those are the ones that are "pits" or "depressions" at a local minimum of potential energy

 

[tony873004]

I made a javascript calculator for finding the L points. It's here:
http://orbitsimulator.com/formulas/LagrangePointFinder.html

 

[D H]

i think you can get an closed form solution, no? the order of the polynomial doesn't exceed 3, does it?

The equation for the L2 point is

\frac{G(M+m)}{R^3}\left(R\frac{M}{M+m}+r\right)\,-\,\frac{GM}{(R+r)^2}\,-\,\frac {Gm}{r^2}= 0

where M and m are the masses of the two massive objects (M > m), R is the distance between the two massive objects, and r, the unknown quantity, is the distance between the smaller object and the L2 point. Denoting α=m/M, x=r/R, this becomes

(1+\alpha)x+1-\frac{1}{(1+x)^2}-\frac {\alpha}{x^2}= 0

Clearing the fraction leads to

((1+\alpha)x+1)x^2(1+x)^2-x^2-\alpha(1+x)^2 = 0

which is a fifth order polynomial.

 

[Irrelativity]

The equation for the L2 point is

((1+\alpha)x+1)x^2(1+x)^2-x^2-\alpha(1+x)^2 = 0

which is a fifth order polynomial.

If this is the case, the above polynomial becomes r = -R/2 if you use (r/R)^2 = approximately zero. But the solution to this question is, I believe, r = R(m/3M)^(1/3). Am I missing something here?

 

[D H]

The equation for the L2 point is

((1+\alpha)x+1)x^2(1+x)^2-x^2-\alpha(1+x)^2 = 0

which is a fifth order polynomial.

If this is the case, the above polynomial becomes r = -R/2 if you use (r/R)^2 = approximately zero. But the solution to this question is, I believe, r = R(m/3M)^(1/3). Am I missing something here?

r=R\,(m/3M)^{1/3} is an approximate solution. How to get there from the fifth degree polynomial:

First rewrite the polynomial this way:

x^2((1+x)^3-1) -\alpha(1+x)^2(1-x^3)= 0

Expanding the first term,

x^3(3+3x+x^2) -\alpha(1+x)^2(1-x^3)= 0

Assuming x\ll 1, 3+3x+x^2\approx 3 and (1+x)^2(1-x^3)\approx 1[/tex], so<br /> <br /> 3x^3 \approx \alpha<br /> <br /> or<br /> <br /> x\approx\left(\frac{\alpha}3\right)^{1/3}

 

[Irrelativity]

r=R\,(m/3M)^{1/3} is an approximate solution. How to get there from the fifth degree polynomial:

First rewrite the polynomial this way:

x^2((1+x)^3-1) -\alpha(1+x)^2(1-x^3)= 0

Expanding the first term,

x^3(3+3x+x^2) -\alpha(1+x)^2(1-x^3)= 0

Assuming x\ll 1, 3+3x+x^2\approx 3 and (1+x)^2(1-x^3)\approx 1[/tex], so<br /> <br /> 3x^3 \approx \alpha<br /> <br /> or<br /> <br /> x\approx\left(\frac{\alpha}3\right)^{1/3}

<br /> <br /> Thanks for your help. But I still not comprehend the last approximation method part. How come x^3 still survives when the other terms vanish?

 

[D H]

Dropping the x² and higher power terms from the original polynomial leads to \alpha(1+2x) \approx 0, which means x \approx -1/2 or \alpha \approx 0 (or both). The former is nonsense and the latter is just a reconfirmation of the assumption \alpha \ll 1. In other words, the naive first-order approximation doesn't work. Refactoring the polynomial as a polynomial in alpha and then eliminating higher order terms in x in each term does work.