Design of low-energy transfers in cislunar space using sequences of lobe dynamics

This paper proposes a graph-based framework that leverages lobe dynamics within the circular restricted three-body problem to systematically design and optimize low-energy cislunar transfer trajectories, which are subsequently refined using the bicircular restricted four-body problem.

The Gist

Imagine you are trying to drive a car from one city to another, but the roads don't exist. Instead, you are floating in space, and the "roads" are invisible, winding currents created by the gravity of the Earth and the Moon. Sometimes these currents flow smoothly, but often they are chaotic, like a river with whirlpools and rapids.

This paper is about a new, clever way to navigate these invisible currents to get from Earth to the Moon using the least amount of fuel possible.

Here is the breakdown of their method using simple analogies:

1. The Problem: The "Whirlpool" of Space

In the space between Earth and the Moon (cislunar space), gravity is a tug-of-war. If you try to fly straight there, you burn a lot of fuel. If you try to "surf" the gravity, you might get stuck in a loop or thrown off course.

Scientists have known for a while that there are "tubes" of gravity (like invisible tunnels) that can carry spacecraft. However, these tubes are hard to find and connect. It's like trying to find a specific, hidden tunnel in a massive, dark cave system without a map.

2. The New Idea: The "Lobe" Map

The authors propose looking at the space not as a few big tunnels, but as a collection of small, swirling pockets called lobes.

• The Analogy: Imagine a river with many small eddies (swirls). A boat can get caught in a swirl, spin around, and then get spit out into a different part of the river.
• The Science: In space, these "swirls" are regions where a spacecraft can get temporarily trapped by the gravity of the Moon, spin around, and then be flung toward a new destination. The authors call these lobes.

3. The Method: Building a "Train Network"

The genius of this paper is how they connect these swirling pockets. Instead of trying to calculate one giant, perfect path from start to finish (which is a nightmare for computers), they built a graph (a network map).

• The Nodes (Stations): Each "station" on their map is the center of one of these swirling lobes.
• The Edges (Tracks): The lines connecting the stations represent the natural paths a spacecraft can take from one swirl to the next.
• The Strategy: They created a system where the spacecraft hops from one lobe to another.
  • Step 1: Jump into a swirl (Lobe A).
  • Step 2: Let the chaos of the swirl carry you naturally to the edge of the next swirl (Lobe B).
  • Step 3: Give a tiny nudge (a small engine burn) to catch the next swirl.
  • Repeat until you reach the Moon.

Think of it like a game of Connect the Dots, but the dots are swirling whirlpools, and you are trying to find the path that requires the least amount of "pushing."

4. Why "Lobe Dynamics" is Better

Previous methods tried to find a single, long, smooth tube. But tubes are rigid; if you miss the entrance by a tiny bit, you miss the whole ride.

Lobe dynamics is more flexible. It's like hopping from stone to stone across a river. Even if you miss a stone slightly, you can adjust your jump to land on the next one. The authors' method uses a computer to check thousands of these "stone-to-stone" combinations to find the absolute cheapest route.

5. The Result: A Fuel-Efficient Trip

They tested this on a trip from Low Earth Orbit (LEO) to Low Lunar Orbit (LLO).

• The Old Way: Usually requires a lot of fuel to break free and aim directly.
• The New Way: By hopping through these chaotic swirls, they found a path that takes about 192 days (a bit slower, but much cheaper) and saves a massive amount of fuel.

They even tested this in a more realistic model that includes the Sun's gravity (which acts like a giant wind blowing on your boat). Even with the Sun's interference, their "whirlpool hopping" method worked, proving it's a robust way to travel.

Summary

Imagine you are a leaf floating down a river.

• Old Method: You try to swim straight to the destination, fighting the current.
• This Paper's Method: You let the river's natural swirls carry you. You just need to know exactly which swirl to jump into and when to make a tiny adjustment to catch the next one.

The authors created a "map" of these swirls and a computer program to find the best sequence of jumps. This allows spacecraft to travel to the Moon like a surfer riding the waves of gravity, rather than a rocket fighting through them.

Technical Summary

1. Problem Statement

Designing fuel-efficient low-energy transfers in cislunar space is challenging due to the highly chaotic nature of the multi-body gravitational environment (specifically the Earth-Moon system).

• Current Limitations:
  • Numerical Approaches: While effective at finding Pareto-optimal solutions, they require accurate initial guesses and suffer from high computational costs.
  • Geometrical/Dynamical Approaches: Traditional methods like Tube Dynamics (utilizing invariant manifolds of libration points $L_1/L_2$) are excellent for global transport between bodies but struggle with departures from the primary body (Earth) and inter-system connections.
  • Lobe Dynamics: Previous applications of lobe dynamics (based on resonant orbits) were limited to specific resonant transitions or moon-to-moon transfers, lacking a systematic framework to combine multiple lobe sequences for complex trajectory design.
• Objective: To propose a systematic method for designing low-energy, fuel-optimal Earth-to-Moon transfers by combining sequences of lobe dynamics associated with multiple resonant orbits, overcoming the limitations of single-resonance or tube-only approaches.

2. Methodology

The authors propose a graph-based framework that integrates dynamical systems theory with combinatorial optimization.

A. Dynamical Models

• Circular Restricted Three-Body Problem (CR3BP): Used as the primary model for trajectory design and lobe identification. The Earth-Moon system is modeled with a mass ratio $\mu \approx 0.01215$.
• Bicircular Restricted Four-Body Problem (BCR4BP): Used in the final validation stage to account for the Sun's periodic perturbation on the Earth-Moon system.

B. Lobe Dynamics & Effective Lobes

• Resonant Orbits: The method utilizes unstable resonant orbits (specifically 7:2 and 3:1 resonances with the Moon) to generate complex phase space structures.
• Poincaré Maps: A Periapsis Poincaré map is constructed in Delaunay elements to visualize the phase space. This map reveals regions of quasi-periodic motion (resonances) and chaotic zones.
• Lobe Identification: Stable and unstable manifolds of resonant orbits intersect to form lobes (enclosed regions).
• Effective Lobes: To ensure robustness against control and navigation errors, the authors define "Effective Lobes" as those with a radius ($r_L$) exceeding a specific threshold ($r^*_L = 0.002$). This ensures that a spacecraft entering a lobe remains within it despite small perturbations.
• Lobe Sequences: A sequence of lobes mapped by the system dynamics ($F$) represents a potential transport path.

C. Graph-Based Optimization

The trajectory design is formulated as a shortest-path problem on a weighted, directed graph:

• Nodes: Represent the departure orbit (LEO), arrival orbit (LLO), and the centroids of the identified effective lobes.
• Edges: Represent transfer paths between nodes.
  • Natural Arcs: Transfers within the same lobe sequence (cost $\approx 0$).
  • Controlled Arcs: Transfers between different lobe sequences or from/to circular orbits, calculated using targeting methods (impulsive $\Delta V$ to match states).
• Constraints:
  • Lobe Constraint: To ensure the exploitation of lobe dynamics, a spacecraft must traverse at least two adjacent lobes within a sequence (preventing "jumping" over the dynamics).
  • Weight Threshold: Edges with $\Delta V$ exceeding a threshold ($w^*$) are discarded to maintain low-energy characteristics.
• Optimization: An exhaustive search is performed to find the path with the minimum total $\Delta V$ (sum of edge weights) from LEO to LLO.

D. Refinement to BCR4BP

The optimal CR3BP trajectory serves as an initial guess for a multiple-shooting optimization in the BCR4BP. The trajectory is segmented, and the Sun's initial phase angle ($\theta^*_s$) is varied to find the global minimum fuel solution in the more realistic four-body environment.

3. Key Contributions

1. Systematic Combination of Lobe Dynamics: Unlike previous studies that focused on single resonant orbits, this paper introduces a method to systematically combine multiple lobe sequences (from different resonances) into a single transfer trajectory.
2. Graph-Theoretic Framework: The construction of a weighted directed graph allows for the efficient exploration of the combinatorial space of lobe sequences, transforming a complex chaotic trajectory design problem into a solvable optimization problem.
3. Robustness via Effective Lobes: The introduction of "Effective Lobes" with a radius threshold provides a practical mechanism to handle control and navigation errors, ensuring the theoretical chaotic paths are physically realizable.
4. Integration of Tube and Lobe Dynamics: The method demonstrates how lobe dynamics can bridge the gap between Earth departure and the lunar transfer tubes (often associated with $L_1$), a connection difficult to achieve with tube dynamics alone.

4. Results

The method was validated through a test problem and a full Earth-Moon transfer scenario.

• Test Problem (CR3BP):

  • Designed a transfer from a 7:2 stable resonant orbit to the stable manifold of the $L_1$ Lyapunov orbit.
  • Case 2 (With Targeting) proved superior to Case 1 (Natural propagation only), correcting propagation errors with small maneuvers ($|\Delta V| < 9$ m/s) and achieving a total cost of 153.25 m/s.
  • The constraint requiring the use of adjacent lobes reduced the search space without significantly increasing the optimal cost.

• Earth-Moon Transfer (CR3BP):

  • Trajectory: From Low Earth Orbit (LEO, 167 km) to Low Lunar Orbit (LLO, 100 km).
  • Performance: Total $\Delta V \approx$ 4274.7 m/s with a flight time of ~192 days.
  • Mechanism: The trajectory utilized a sequence of lobes (7:2 and 3:1 resonances) to gradually increase the semi-major axis and transition resonance ratios, effectively navigating the chaotic zone before entering the lunar sphere of influence.

• Earth-Moon Transfer (BCR4BP):

  • The CR3BP solution was refined in the Sun-Earth-Moon system.
  • Best Solution: With a Sun phase angle $\theta^*_s = 0$, the total $\Delta V$ was reduced to 3832.6 m/s (flight time ~193 days).
  • Comparison: This result is competitive with existing literature solutions for interior transfers, demonstrating the method's effectiveness in a realistic dynamical environment.

5. Significance

• New Paradigm for Trajectory Design: The paper moves beyond the traditional "tube dynamics" (libration point manifolds) approach, offering a structured way to utilize resonant lobe dynamics for complex transfers.
• Scalability: The graph-based approach is scalable and can be extended to other multi-body systems (e.g., Jupiter's moons) or different mission constraints (time vs. fuel).
• Practical Application: By validating the method in the BCR4BP and comparing it to known solutions, the authors prove that lobe dynamics is not just a theoretical tool but a viable strategy for designing real-world, fuel-efficient cislunar missions (relevant for Artemis and future lunar logistics).
• Insight on Weak Stability Boundary (WSB): The study clarifies the relationship between lobe dynamics and the WSB, showing that while WSBs indicate weak capture, lobe dynamics provides a more precise mechanism for predicting and controlling resonance transitions.

In conclusion, this work provides a robust, systematic, and computationally efficient framework for designing low-energy transfers by harnessing the chaotic transport properties of multiple resonant orbits, significantly advancing the state-of-the-art in cislunar trajectory design.