3

I have to build an optimal trajectory between two orbits, which are in one plane, but with different semi-major axis, argument of perigee and eccentricity.

I've read different articles regarding the low-thrust optimization techniques, and now I'm looking for an example solution.

I would appreciate for a link/paper describing an example solution of the above problem.

UPDATE

I have an impulsive solution of this problem. Is it possible to convert the impulsive solution to the low-thrust?

Leeloo
  • 833
  • 11
  • 29
  • 1
    I have not seen a conclusive proof even for the instantaneous variant: https://space.stackexchange.com/questions/16931/is-the-ideal-transfer-between-two-arbitrary-planar-orbits-always-a-bi-tangential, despite spending a lot of rep trying to attract a proof. – SE - stop firing the good guys Sep 01 '20 at 16:02
  • 1
    Also, the geometry of such orbits: https://math.stackexchange.com/questions/2168522/ellipse-tangential-to-two-arbitrary-ellipses-sharing-a-focus – SE - stop firing the good guys Sep 01 '20 at 16:04
  • 2
    @SE-stopfiringthegoodguys - A bi-elliptical transfer to which you alluded involves three instantaneous burns, while a bi-tangential transfer involves only two. (A Hohmann transfer is a special case of a a bi-tangential transfer.) There are cases where a bi-elliptical transfer requires less delta V than a bi-tangential transfer, at the expense of taking more time, sometimes significantly more time. – David Hammen Sep 01 '20 at 16:22
  • 1
    @DavidHammen bi-elliptical transfers (with infinite apoapsis always minimizing delta-v), is in a sense also a bi-tangential transfer. It's just the "easy" solution. – SE - stop firing the good guys Sep 01 '20 at 18:40
  • @SE-stopfiringthegoodguys - A bi-elliptical transfer with infinite apoapsis is for all practical purposes a non-solution as it takes an infinite amount of time. – David Hammen Sep 02 '20 at 03:47
  • @uhoh Removing the check mark on that one is probably a good idea in any case. Not sure it's worth spending more on at the moment. I plan doing a more systematic series of bounties myself later to clear up a couple of related orb-mech problem. – SE - stop firing the good guys Sep 02 '20 at 08:01
  • Have you already addressed what it is that you wish to optimise? e.g. lowest propellant mass, or lowest propellant mass for a given transfer time, or lowest transfer time for a given propellant mass. This will fundamentally change the problem. – Puffin Sep 03 '20 at 13:33

1 Answers1

1

There are too many variables, some of them infinitely variable, to arrive at the optimal burn plan. And then you need to incorporate that burn plan into something that can be carried out with a reasonably priced computer (US $200K is "reasonably priced) than can handle high level radiation. This is not going to be a state of the art computer.

What you can do is use techniques that come close to the theoretical optimum but with constraints that keep the profile manageable by a not quite up to date computer. Mischa Kim's PhD thesis suggests using Adaptive Simulated Annealing. Markov Chain Monte Carlo is very similar to simulated annealing, and is perhaps easier to set up. Another related approach is a particle filter. All use a variant of the Metroplis algorithm, aka the Metropolis-Hastings algorithm.

Suppose you know that it will take multiple orbits from the initial orbit to the final desired orbit. In this case, you do not want the spacecraft firing continuously. You instead want the spacecraft to take advantage of the Oberth effect. In the case of very low thrust, this means firing for roughly 2/3 of the orbit, roughly centered around periapsis, plus perhaps a final finite thrust when the spacecraft nears the desired orbit.

I would do this in stages. First assume impulsive burns at periapsis, with the impulse equivalent to what could be achieved by the finite thrusters over the course of 2/3 of an orbit, with a penalty for gravity losses. Use a search technique to find a tentative multi-burn solution. This becomes the initial guess for the multiple finite burn solution. The transformation from impulsive burns to finite burns will miss the target.

So now you use similar techniques to refine the transformation. The nice thing about the Metropolis algorithm, and its variants, is that the randomness applied to tentative solutions kicks things away from local minima that in fact are lousy solutions.

David Hammen
  • 74,662
  • 5
  • 185
  • 283
  • I have an impulsive solution. Is there a way and would it be correct to convert it to a low-thrust solution? – Leeloo Sep 02 '20 at 12:49