Freereturn trajectory around the moon /
Flying to the moon with Apollo
This article is part of a serie: Freereturn trajectory around the moon
 Flying to the moon with Apollo
About 50 years ago, humans landed on the Moon for the first time. I initially wanted to write on this subject because I like rockets and space, but this anniversary is a nice coincidence.
Going to the Moon is not an easy endeavor. It requires an outstanding number of subsystems and components to work flawlessly, and many maneuvers to execute right as planned. Each time a rocket takes off, an engineering miracle happens.
I have studied one particular aspect of this engineering problem, the trajectory, and I will share with you my findings. The goal of this article is to answer this simple question: How did Apollo go to the Moon?
Fly me to the moon
Let me play among the stars
Let me see what spring is like
On a Jupiter and Mars.– Bart Howard, 1954
What is a freereturn trajectory
Requirements
The objective is pretty simple, and can be traced back from one very important requirement: try not to kill astronauts^{1}.
The trajectory that takes the astronauts to the Moon should allow a safe return to the Earth in case of any problem, like an engine failure. Actually, the trajectory should be free, from a propulsive viewpoint, hence its name: Freereturn trajectory.
How does a freereturn trajectory works?
In theory, there is no friction in space, therefore no energy loss. That means that when a spacecraft is orbiting a body (for instance, the Earth), at some point it will come back at the exact same point. In real life, there are many forces (atmospheric drag, gravitational pull due to the imperfect shape of the Earth, solar pressure, etc) that disturb an orbit, and their magnitude depends on the location and physical properties of the spacecraft.
The spacecraft traveling away from a primary body (in our case, the Earth), follows a path that will be altered by the gravitational pull of the secondary body (in our case, the Moon). A freereturn trajectory works by predicting the trajectory, taking into account the gravitational influence of the Moon.
Of course, in practice no trajectory can be truly free and without propulsion, there will be some small midcourse corrections. Although Apollo missions had many midcourse maneuvers planned, the ones that were small enough were skipped.
The following graphs show the trajectory without (standard elliptic orbit) and with (free return trajectory) the gravitational pull of the Moon.
The total journey is much shorter with the Moon (7.5 vs 13 days), because the spacecraft accelerate as it falls towards the Moon. During this time elapsed, the Moon travels around 25% vs 50% of its orbit, which is shown in grey. The outbound part last 46 days, during which the Moon travels around 60 degrees (a full revolution is performed in 27.3 days). It can be seen that without the Moon, the orbit is not perturbed: the spacecraft comes back at the same position and velocity. With the Moon, the orbit is perturbed: the second orbit will be different from the first one.
This type of trajectory was known much earlier than Apollo missions. It was used for the first time in 1959 (while Apollo 11 happened in 69), when the Soviet probe Luna 3 had to send back to Earth the photographs it had taken^{2}.
Apollo 8, 10 and 11 used a freereturn trajectory. Due to the landingsite restrictions, subsequent Apollo missions used a hybrid trajectory that was effectively a freereturn. They then performed a midcourse maneuver to change to a transLunar trajectory that was not a free return.
Calculating the trajectory
Introduction
While twobody orbital mechanics (one spacecraft an another body) is easy, threebody orbital mechanics (spacecraft, Earth and Moon) is hard: there is no simple equations to model the trajectories. There are two main solutions:
 Splitting the trajectory into several two body approximations. This method is called patched conics^{3}.
 Using an iterative method, also called restricted circular three body (RC3B) problem^{4}.
Patched conics
This method only takes into account the strongest gravity force to which the vessel is subject. The idea is to determine the sphere of influence of each body, and use standards orbital mechanics equations to compute the trajectory of the spacecraft, for the current sphere of influence. When the spacecraft transitions from one sphere of influence to another, the equations are recomputed with the new body.
Patched conics yields very approximate results, but it can be fine in most cases^{5}. Motion is deterministic and simple to calculate, which is very useful for rough mission design and back of the envelope studies.
In our case, when leaving the Earth, the spacecraft is assumed to travel only under classical two body dynamics. Its acceleration is dominated by the Earth until it reaches the Moon’s sphere of influence, at which point the Earth is ignored and only the Moon’s influence is taken into account.
The following image visually explains the concept of sphere of influence.
The formula to compute the size of the sphere of influence can be approximated with various amounts of precision. For more details, you can refer to the Wikipedia page.
Restricted three body problem
The other method, the restricted circular three body problem, can be fine tunned to reach the required level of precision (by changing the size of the iteration step). If correctly implemented, it can be very fast while retaining a very good level of precision.
The method is iterative: given a (relatively) small time increment (which can range from seconds to hours, depending on the time available for the computation and the precision required), the acceleration of the spacecraft is computed. In its simplest form, the acceleration is derived from the gravitational force of the bodies (Earth and Moon in our case), but other forces can also be taken into account, such as modeling the Moon’s true orbital motion, gravitation from other astronomical bodies, the nonuniformity of the Earth’s and Moon’s gravity, including solar radiation pressure, etc. Propagating the spacecraft motion in such a model is numerically intensive, but necessary for true mission accuracy.
This method does not have an analytic solution, and requires numerical calculation with methods such as Euler or RungeKutta integration. Even if you do not know what I am talking about, if you have a basic engineering background and were to implement this simulation method, you would very probably write an Euler integrator. I will not go into these various integrators (Euler, RungeKutta, etc) and optimizations details in this article, but I will keep them for a future article.
Implementing the RC3B method
In this section, a software implementation of the restricted circular three body (RC3B) problem will be presented. This method offers a good compromise, as it can be reasonably fast for initial studies, while keeping the possibility of increased accuracy (at the cost of an increased computation duration, of course). The initial implementation is also very simple, and the accuracy can latter be increased by taking into account additional perturbation forces.
In order to implement this and for the sake of simplicity, I will use a few assumptions and simplifications that are not exactly true in real life. They should not affect the results too much, and I will make sure to mention them along the way.
The entire code is around 250 lines in its final version, with comments, debug output, proper constants definition, helpers, etc. The interesting part is much shorter, and we will study it together, piece by piece. There are two small code snippets to set up things, and two others for the core of the algorithm.
Coordinate system
First, a coordinates system will be defined with an orthonormal basis. The simulation will be based on a 2D plane, with (x, y) coordinates expressed in meters. Even if that means working with huge numbers, you should always use standard units. Otherwise, sad things will happen.
Here is a first approximation: the real world is in 3 dimensions, and even if the orbit of the Moon is in a plane, the trajectory is not. The launchpad from where the rocket will take off will not be in the Moon orbit plane.
Python has awesome math/scientific libraries such as numpy
and scipy
, but I reimplement a Vector
class.
This will allow me to define custom methods such as from_polar_rad()
.
I will also implement basic operators thanks to Python’s magic methods, which will allow us to add/sub/mul/div our custom objects.


The choice of the origin of our basis is important, as it can make our computations easy or very hard.
Since the goal is to visualize the trajectory of the spacecraft with respect to the Earth and the Moon, these two bodies should be fixed in our coordinates system.
This leads us to use the barycenter EarthMoon as the origin (the spacecraft has a negligible mass).
The barycenter is easily calculated with: distance * mass2 / (mass1 + mass2)
.
To fix the Earth and the Moon in our coordinates system, the whole space will rotate at an angular speed of by 2*pi / (27.3*24*60*60)
radians per seconds (which corresponding to the angular speed of the Moon.
Initialization
Now that the foundation of our simulation are set, we can build from here.
First, I will define a class Body
for the Earth, the Moon and our spacecraft.
It will hold their position, velocity and acceleration (remember, an iterative method is used, where the acceleration is integrated twice to obtain the position).
This class will be presented just after this initialization code snippet.
The Earth’s and Moon’s properties are straightforward:
 The position (
earth_pos
andmoon_pos
) is on the xaxis, on the left of the barycenter for the Earth and on the right of the barycenter for the Moon.  The velocity (
earth_vel
andmoon_vel
) is a function of the distance to the barycenter and its sign is set for the body to have a counterclockwise direction (y for the Earth and +y for the Moon).
The spacecraft properties are a little more complex. For the ease of the calculations, its position will be based on the Earth. The injection define the position and speed applied to the spacecraft for it to transfer from a stable Low Earth Orbit to an orbit toward the Moon.


The spacecraft’s mass is not used, therefore it is set to 1
, for safety, to avoid any division by zero.
The spacecraft’s velocity direction and position are a little hard to understand with the numbers, but with the following diagram it will be easier.
Simulation’s main loop
The two previous code snippets have set up the simulation, now comes the first part of the core of the simulation.
This code snippet loops on the given simulation time step sim_dt
.
It works in two steps: first, the acceleration is computed, without updating any position or velocity.
Then, the velocity and position of each body is updated at the same time.
The loop would be a good place to save an history of the position of each body, in order to draw the trajectory once the simulation is finished.


A note about the while
keyword used instead of the perhaps more suitable for
keyword.
I used it because it felt more suited for this particular case.
The simulation and the export of the positions might not be synchronized, and I changed the logic many times.
This way of doing was easier to rework and experiment with than a more classical for
.
Euler integrator
This is the final code snippet, where all the magic happens. Do not skip this section too quickly, it is surprisingly easy to understand. Even if all the magic is here, it will not look like nonsensical incantations that are used to summon daemons.
As recalled from the previous code snippets, two functions are needed:
update_acc()
to update the accelerationupdate_state()
to update the velocity and position
The update_state()
function is very straightforward: with a basic engineering background, the Euler integrator should look familiar.
The update_acc()
function has a little more lines, but is also simple: it computes the acceleration by adding the gravitation from each body, using the equation for universal gravitation.


Checking the coherence of units is very helpful.
The equation for universal gravitation is expressed in Newton
which is m*kg/s**2
, but the acceleration is expressed in m/s**2
.
It is clearly needed to divide by kg
at some point, which makes sense to compute an acceleration.
There is another approximation here: the Earth is not a true sphere and the gravity is not constant everywhere. Real mission design will use a more precise model of the gravity of the Earth.
Full code
To keep the code short I have omitted a few non essential things: the definition of the constants (Earth and Moon properties, as well as some simulation parameters), the full Vector
class (which does what you would expect and nothing more) and the export of the state (to plot the trajectory, a rather boring piece of code).
There is definitely no hidden magic.
The code is a little long to show here (around 250 lines), but it is available here.
Output
When running this script, it outputs the following result.


Ignoring the debug prints, it says that it took a little more than 3 seconds^{6} to compute the final position of the spacecraft, which is (37 398 968 ; 320 058 357) meters.
The code is configured to simulate 10 days and output the state of the spacecraft every 50 seconds, which gives us 17280 points. these states are written to a file, which can then be used to plot the trajectory of the spacecraft. This allow to run a lengthy but precise simulation and plot it quickly without simulating it again and again.


Our example can be plotted with this export of the positions. The following two diagrams show the same simulation, one with the trajectory of the spacecraft in a fixed frame (with respect to the Earth), and the other in a rotating frame (in order to keep the Moon fixed).
The Moon looks surprisingly close, although it is still several days away …
Exploring the parameter space
Now that I have a simulation tool, I wanted to know which trajectories were viable or not for a journey to the Moon. Why not run a few thousands simulations and find out what is the best strategy for outerspace travel?
A spacecraft beginning its journey to the Moon can tune two parameters to adjust its trajectory:
 Translunar injection angle: where on the orbit (the angle MoonEarthspacecraft)
 Translunar injection speed (or deltav): how much energy is imparted to the spacecraft (the speed added to the initial orbital speed)
Here is again the schematic of the position and velocity of the spacecraft to really understand these two parameters.
Impact of the injection parameters
These two parameters are very chaotics: a small change will change the final position by a (much) large(r) amount. With just a few m/s or degrees, or even a tenth, the final position of the spacecraft will change by several hundreds of km. In real life, several midcourse correction maneuvers will take care of adjusting the trajectory to ensure a safe return and landing.
Here are a few examples showing the impact of varying these two parameters.
A crazy precision is required: an error as small as 0.1 m/s or 0.1 ° (0.1 deg corresponds to 1.5 sec in Low Earth Orbit) is the difference between a correct landing and big troubles (toasted and crashed, or simply frozen after a close Earth flyby: neither offer even the slightest survival chances). The diagrams shows a 200 km Low Earth Orbit. Do not try to find it on the first one, it is too close to the surface to be visible.
In order to perform a flyby of the Moon and come back to the Earth, the injection deltav must be between 3 050 and 3 125 m/s (corresponding to a velocity between 10 840 and 10 920 m/s with respect to the Earth, for trajectories passing between 2 000 and 20 000 km from the Moon^{7}). Empirically, the injection angle should be between 120 and 130 degrees.
Finding out how not to miss the Earth on the inbound journey
With this short explanation of the injection parameters, we can now try to answer the initial question: what is the outcome when these parameters are varied on a larger spectrum? What are the limits, and at which values does the spacecraft crashes on the Moon, injects into an heliocentric orbit or crashes on the Earth?
I ran a huge number of simulations (a few thousands) with different parameters: for each injection deltav and angle I saved the closest approach to the Moon and the Earth. At around one second per simulation, it took a few minutes to have a rough idea and less than an hour for detailed results.
The result is plotted in 3D below. The horizontal plane is for the two parameters, and the vertical axis is used for the distance to the Earth: red for crashes, green for skimming through the atmosphere, grey for a close flyby and finally black for very far.
The two plots show the closest approach to the Moon and to the Earth. Since we want to focus on the safe return of our astronauts, we will focus on the Earth closest approach. The closest approach to the Moon is far less interesting, and any trajectory that end in the Earth atmosphere will pass by the Moon at a great distance (~10 000 km). In real life, the Moon closest approach would be very critical, due to landing site precision and safety requirements.
The graph shows that the injection angle affects less the result than the injection deltav. The critical parameter is the injection deltav: the “safe” entry distance (in green) happens only in a verry narrow window around 3150 m/s. The distances at the edge of the 4 greyblack zones are the usual satellite altitudes: <5 000 km (Low Earth Orbit) for most observation satellites, 20 000 km (Medium Earth Orbit) for GPS and 36 000 km (Geosynchronous Equatorial Orbit) for telecom.
A final note on the “safe” entry distance aforementioned: the green section of the plot corresponds to an altitude between 0 and 100 km. This value was chosen as a tradeoff between fair approximation and being able to see something. Compared to real life, this is way too large, but otherwise we would not see anything. A real spacecraft would probably reenter with a projected closest approach between something like 30 and 50 km^{8}. Actually engineers do not focus on the closest approach (whose fancy name is perigee, when talking about the Earth), but rather on an entry interface. The entry interface specifications include the altitude, but also the speed and the angle of attack. A typical Apollo Earth entry interface was 6.5 degrees at 120 km altitude, with an acceptable range from 5.2 to 7.7 degrees^{9}.
Conclusion
In this article, I presented the trajectory used by the Apollo missions in 1969 and the following years. This freereturn trajectory allowed the astronauts to reach the Moon without compromising their safety in case of a serious issue. Apollo 13 took advantage of this feature to return safely to the Earth, without using the CSM (command and service module) main engine after the explosion of their oxygen tank, by fear that it could worsen the situation due to malfunction or even explosion. I also explained which methods can be used to simulate this trajectory, and how I implemented the one which uses the restricted circular three body problem. I showed how very small deviations of the input parameters (injection angle and deltav) can have a huge impact on the return closest approach to the Earth.
Source code and data
I am a firm believer of open source, therefore the main source code used for the simulation can be found here, as well as all the code used to generate the diagrams and graphs. Since the simulator is rather slow, I also included the raw data that I used in my plots, so you can quickly plot and play with the data. You are welcome to explore and improve my work and send me any result you may find worth sharing!
Additional material
When writing this article, I found a some interesting material that I tough were worth sharing. Here they are:
 “Fly me to the Moon” by Jack Crenshaw: an awesome explanation of freereturn trajectories, well worth reading even after this article  Original link and Webarchive link
 “Trajectories in the EarthMoon space with symmetrical free return properties” by Arthur J. Schwaniger: advanced analysis of freereturn trajectories, with some nice graphs  Orignal link
 “Introduction to Lunar FreeReturn Trajectories” by Max Fagin: if you like videos  Youtube link
 “A Simple Targeting Procedure For Lunar TransEarth Injection” by Shane Robinson and David Geller: hard maths and trajectory optimization  Original link
Upcoming article: optimization and benchmark
I initially wrote the simulator in Python, but was quickly put out by the slowness of the language. For calculation heavy applications like this one, Python is not the best choice. Nonetheless, it has allowed me to iterate quickly and promptly achieve a working prototype. It was also very easy to create plots of the results. I started trying to optimize the algorithm, but with Python it was hopeless. Since I am currently learning Rust, I started rewriting the simulator with this language. At this point I already had a few different versions and started comparing their speed. I went down a rabbit hole, and I plan to share my findings in an upcoming article that will probably be titled Freereturn trajectory around the Moon: optimization and benchmark.
 This was probably NASA’s second requirement, just after Get on the Moon before the Russians. ^{[return]}
 https://en.wikipedia.org/wiki/Freereturn_trajectory#Earth–Moon ^{[return]}
 https://en.wikipedia.org/wiki/Patched_conic_approximation ^{[return]}
 https://en.wikipedia.org/wiki/Translunar_injection#Restricted_circular_three_body_.28RC3B.29 ^{[return]}
 Actually, it is the method used in the well known game Kerbal Space Program ^{[return]}
 Do not focus to much on the 3 second figure, I am keeping the optimization and benchmark considerations for another article ;) ^{[return]}
 https://en.wikipedia.org/wiki/Freereturn_trajectory ^{[return]}
 To be honest I have no clue. 3050 km is the usual value when playing KSP. ^{[return]}
 space.stackexchange.com and history.nasa.gov/afj ^{[return]}