Wednesday, June 10, 2015

Mass parameter and ITN

It seems like every other post I'm singing the praises of EML2. I'm also enthusiastic about L1 and L2 necks for big moons orbiting gas giants.

So why do I diss the Sun Earth L2 or the Sun Mars L1? 

Robert Walker put it fairly well:

I've no idea why you think there's some essential difference between e.g. transfers between moons in the Jupiter system and transfer between planets around the sun. Mathematically it's the same situation, multiple masses around a central planet or sun. Obviously the moons of Jupiter are larger compared with Jupiter than planets are compared to the sun, and the orbits are far shorter. But they still have Hohmann transfer orbits, and hill spheres, and lagrange points, and these tubes, which lead out from the lagrange points. 

It's The μ



The big reason is mass parameter. This quantity is often denoted μ in discussion of 3 body mechanics.


It's common to choose units so that mass of central and orbiting body sum to one.

For example if central mass were 90% percent of the system's mass, central mass would be .9 and orbiting body would be .1. In this case μ would be 1/10.

The small the μ, the closer L1 and L2 get to the orbiting body.

What paths do payloads follow when nudged away from the orbiting body at L1 or L2? Well, we can notice a few things about L1 and L2:

The L1 and L2 have the same ω (angular velocity) as the orbiting body.

L1 and L2 are collinear with central and orbiting bodies.

It just so happens I have a diagram of collinear points all having the same ω. It's what I use to model vertical tethers. By scaling this diagram it could also be used to model space elevators. (Space elevators are a special case of vertical tether where the tether foot coincides with planet surface and circular orbit of the balancing point coincides with planet synchronous orbit):

Eccentricity Vertical Tether Conics = |1 - r3|

Release a payload from any point on the tether and the path will be conic section having eccentricity
|1 - r3| where distance from center to balancing point is 1 and a point's distance from center is r.

But this diagram was derived using 2 body mechanics. When nudged away from the orbiting body's Hill Sphere, the payload will quickly enter a regions where the central body gravity dominates and conic sections are fairly accurate.

But while in the neighborhood of the Hill Sphere, there's a short interval when the path should be modeled using both the accelerations of central and orbiting body:



While falling away from the moon near L1,  a payload surges ahead while the moon tugs it backwards and a little up. This has the effect of lowering the apo and periapsis as well as rotating line of apsides in a prograde direction. A payload nudged from L2 away from the moon will lag behind the moon. While in the lavender region, the moon pulls the L2 payload forward boosting peri and apo-apsis as well as rotating line of apsides forward. L2 necks throw higher and L2 drops lower than corresponding points from a vertical tether.

Here are orbital sims for various mass parameters where colored pellets are nudged with slightly different velocities from L1 and L2:

μ = .1

μ = .001

μ = .00001

μ = .000001

Notice as μ shrinks, the orbits get closer and closer to what the tether model would suggest. For μ = .000001, apogee is very close to point of release and eccentricity of ellipses is approaching |1 - r3|. As the Hill Sphere shrinks the lavender two body zone gets thinner and the 2 body model becomes increasingly accurate.

μ for sun-earth is .00000304 and μ for sun-Mars is .000000323. The paths from the planets' L1 and L2 necks don't go far and there's not much variation.

What About Gravity Assists? Just Look At Rosetta

"Well sure, nudging a payload from SEL2 doesn't get us much past a 1.07 A.U. aphelion" an ITN defender replies. "But earth gravity assists can boost that aphelion. Look at Rosetta's March 2005 gravity assist -- it boosted an earth like orbit to an aphelion past Mars."

So let's look at the Rosetta gravity assist.



We can see the March 2005 gravity assist gets Rosetta past Mars orbit. And the orbit from launch to gravity assist looks pretty earth like, right?

No. 

The orbit from launch to gravity assist is an ~.9 x 1.1 A.U. ellipse with an eccentricity of around .11. When r = 1 A.U., flight path is about 5 degrees. As it approaches earth, Vinfinity is about 2.6 km/s:



In contrast, Vinf of an orbit departing from SEL1 or 2 will be about .3 km/s. Rosetta's initial orbit couldn't be accomplished via a WSB from an L1 or L2 neck.

Further, a payload from SEL2 remains outside earth's orbit, it does not cross. The closest it comes it .01 A.U. during which time it's flight path is zero. The same is true of an orbit nudged from SEL1:


Synodic Period

Another thing to consider is synodic period. A way to think of synodic period is how often one runner laps another as they race about a circular track. If both runners are going nearly the same speed, it will take a long time.

Synodic period of orbiting bodies is |(T1 * T2)/(T1 - T2)| where T1 and T2 are bodies' orbital periods. Orbital period of a 1.01 x 1.06 ellipse is 1.053 years. Synodic period is 19.88 years. So the payload wouldn't even come close to the earth until almost two decades later!

When the payload finally does lap the earth 19.88 years later, it will be 43º from perihelion.


Instead of being .01 AU from earth, the payload will be more like .02 or .03 A.U. from earth. It won't get close to earth until nearly 8 synodic periods later. That's about 160 years.

The tinier the μ, the bigger the synodic periods of payloads released from L1 and L2 necks.

Synodic Period with a big μ

In closing I'll take a look at synodic period of something dropped with from Earth Moon L1. When it comes to μ's the earth moon's .012 is the 900 pound gorilla of the solar system. It's the biggest I know of except for Pluto Charon's .104.

Dropping from EML1, payloads fall into an approximately 100,000 x 300,000 km orbit:


Period is about 11 days. Synodic period is 20 days. So within a month's time these pellets will fly by the moon when they're near apogee:




EML1 and EML2 can do lots of stuff within a fairly short time. I have seen some crazy stuff running earth moon sims.

But zoom out and it gets boring. I've let sun earth sims run for centuries without seeing any drama. The L1 and L2 necks for the sun/rocky planets are a bunch of duds.

A few mass parameters

Here are a few mass parameters for central and orbiting bodies:


Pluto/Charon 1.043E-01
Earth/Moon 1.216E-02
Sun/Jupiter 9.545E-04
Sun/Saturn 2.856E-04
Saturn/Titan 2.374E-04
Jupiter/Ganymede 7.789E-05
Jupiter/Callisto 5.684E-05
Sun/Neptune 5.153E-05
Jupiter/Io 4.700E-05
Sun/Uranus 4.366E-05
Jupiter/Europa 2.526E-05
Saturn/Rhea 4.046E-06
Sun/Earth 3.039E-06
Sun/Venus 2.448E-06
Saturn/Dione 1.935E-06
Saturn/Tethys 1.091E-06
Sun/Mars 3.229E-07
Saturn/Enceladus 1.935E-07
Sun/Mercury 1.659E-07
Saturn/Mimas 7.037E-08
Mars/Phobos 1.682E-08
Sun/Pluto& Charon 7.149E-09
Mars/Deimos 2.803E-09
Sun/Ceres 4.741E-10

I hope I'll some time to play with the Pluto/Charon 3 body system. I believe there are some wonderful possibilities in this setting.

Does the ITN include Gravity Assists?

Given big enough  μ's, some WSBs can snake by another body which can lend a gravity assist.

Does that mean gravity assists are part of the ITN? No. We've been using gravity assists for many decades. They were in common use before Ross' or Belbruno's techniques came on the scene.

Some claim astrogators are now using Belbruno's or Ross' techniques to find opportunities for gravity assists. I haven't seen any evidence of this. So far as I can tell mission planners are still finding such opportunities the old fashioned way: looking for needles in a hay stack. In other words with persistence and hard work.



No comments: