Search

Equipping a spacecraft with multiple solar-powered electric engines (of the same or different types) compounds the task of optimal trajectory design due to presence of both real-valued inputs (power input to each engine in addition to the direction of thrust vector) and discrete variables (number of active engines). Each engine can be switched on/off independently and "optimal" operating power of each engine depends on the available solar power, which depends on the distance from the Sun. Application of the Composite Smooth Control (CSC) framework to a heliocentric fuel-optimal trajectory optimization from the Earth to the comet 67P/Churyumov-Gerasimenko is demonstrated, which presents a new approach to deal with multiple-engine problems. Operation of engine clusters with 4, 6, 10 and even 20 engines of the same type can be optimized. Moreover, engine clusters with different/mixed electric engines are considered with either 2, 3 or 4 different types of engines. Remarkably, the CSC framework allows us 1) to reduce the original multi-point boundary-value problem to a two-point boundary-value problem (TPBVP), and 2) to solve the resulting TPBVPs using a single-shooting solution scheme and with a random initialization of the missing costates. While the approach we present is a continuous neighbor of the discontinuous extremals, we show that the discontinuous necessary conditions are satisfied in the asymptotic limit. We believe this is the first indirect method to accommodate a multi-mode control of this level of complexity with realistic engine performance curves. The results are interesting and promising for dealing with a large family of such challenging multi-mode optimal control problems., Comment: 45 pages, 11 figures; Preprint submitted to Acta Astronautica

Efficient performance of a number of engineering systems is achieved through different modes of operation - yielding systems described as "hybrid", containing both real-valued and discrete decision variables. Prominent examples of such systems, in space applications, could be spacecraft equipped with 1) a variable-$I_{\text{sp}}$, variable-thrust engine or 2) multiple engines each capable of switching on/off independently. To alleviate the challenges that arise when an indirect optimization method is used, a new framework --- Composite Smooth Control (CSC) --- is proposed that seeks smoothness over the entire spectrum of distinct control inputs. A salient aftermath of the application of the CSC framework is that the original multi-point boundary-value problem can be treated as a two-point boundary-value problem with smooth, differentiable control inputs; the latter is notably easier to solve, yet can be made to accurately approximate the former hybrid problem. The utility of the CSC framework is demonstrated through a multi-year, multi-revolution heliocentric fuel-optimal trajectory for a spacecraft equipped with a variable-$I_{\text{sp}}$, variable-thrust engine., Comment: 51 pages, 9 figures; pre-print submitted to Acta Astronautica

A central problem in orbit transfer optimization is to determine the number, time, direction and magnitude of velocity impulses that minimize the total impulse. This problem was posed in 1967 by T. N. Edelbaum, and while notable advances have been made, a rigorous means to answer Edelbaum's question for multiple-revolution maneuvers has remained elusive for over five decades. We revisit Edelbaum's question by taking a bottom-up approach to generate a minimum-fuel switching surface. Sweeping through time profiles of the minimum-fuel switching function for increasing admissible thrust magnitude, and in the high-thrust limit, we find that the continuous thrust switching surface reveals the $N$-impulse solution. It is also shown that a \textit{fundamental} minimum-thrust solution plays a pivotal role in our process to determine the optimal minimum-fuel maneuver for all thrust levels. We find the answer to Edelbaum's question is not generally unique, but is frequently a set of equal-$\Delta v$ extremals. We further find, when Edelbaum's question is refined to seek the number of finite duration thrust arcs for a specific rocket engine, that a unique extremal is usually found. Numerical results demonstrate the ideas and their utility for several interplanetary and Earth-bound optimal transfers that consist of up to eleven impulses or, for finite thrust, short thrust arcs. Another significant contribution of the paper can be viewed as a unification in astrodynamics where the connection between impulsive and continuous-thrust trajectories are demonstrated through the notion of \textit{optimal switching surfaces}, Comment: 69 pages, 79 figures, 13 tables, Preprint submitted to the Journal of the Astronautical Sciences, also the Keynote paper at the AAS John L. Junkins Astrodynamics Symposium, Proceedings in press

We have developed a new method for solving low-thrust fuel-optimal orbit transfer problems in the vicinity of a large body (planet or asteroid), considering a high-fidelity spherical harmonic gravity model. The algorithm is formulated via the indirect optimization method, leading to a two-point boundary value problem (TPBVP). We make use of a hyperbolic tangent smoothing law for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and thus promote convergence. The TPBVP is solved using the method of particular solutions (MPS) shooting method and Picard-Chebyshev numerical integration. Application of Picard-Chebyshev integration affords an avenue for increased efficiency that is not available with step-by-step integrators. We demonstrate that computing the particular solutions with only a low-fidelity force model greatly increases the efficiency of the algorithm while ultimately achieving near machine precision accuracy. A salient feature of the MPS is that it is parallelizable, and thus further speedups are available. It is also shown that, for near Earth orbits and over a small number of en-route revolutions around the Earth, only the zonal perturbation terms are required in the costate equations to obtain a solution that is accurate to machine precision and optimal to engineering precision. The proposed framework can be used for trajectory design around small asteroids and also for orbit debris rendezvous and removal tasks.