Optimization of Low-Thrust Gravity Assist Trajectories Under Application of Tisserand’S Criterion
Total Page:16
File Type:pdf, Size:1020Kb
OPTIMIZATION OF LOW-THRUST GRAVITY-ASSIST TRAJECTORIES UNDER APPLICATION OF TISSERAND’S CRITERION Vom Fachbereich Produktionstechnik der UNIVERSITÄT BREMEN zur Erlangung des Grades Doktor-Ingenieur genehmigte Dissertation von Dipl.-Ing. Volker Maiwald Gutachter: Prof. Dr.-Ing. Andreas Rittweger Fachbereich Produktionstechnik, Universität Bremen Prof. Dr.-Ing. Bernd Dachwald Fachbereich Luft- und Raumfahrttechnik, Fachhochschule Aachen Tag der mündlichen Prüfung: 23. März 2018 I think we are going… because it’s in the nature of the human being to face challenges. - Neil Armstrong Acknowledgements My gratitude goes towards my colleagues and friends from the System Analysis Space Segment Department for their support and fellowship. I thank especially Dominik Quantius, Andy Braukhane, Daniel Schubert and Conrad Zeidler for their discussions, support and proof reading concerning my dissertation. I also thank Dr. Marco Sharringhausen, who always helped me whenever I had questions about mathematics and Dr. Wolfgang Seboldt for stepping up on a short notice during the final stages of this thesis and providing much appreciated insight and advice. I thank my parents, Dr. Ursula and Werner Maiwald for enabling me to follow my dreams – I did, thank you! – and my brothers Rüdiger and Tobias for their encouragement. And of course, I thank my incredible wife Inga who patiently endured my struggles with my work and supported me so much, despite the bumpy road far off an optimal path. Revisions for Publication This is a revised edition of the original dissertation. The changes have been made to allow a more self-enclosed reading of the dissertation. Next to editorial changes to address language and typing errors, the following changes have been made to the originally submitted version: 1) The paper by Chen, Kloster and Longuski (2008) has been included in the discussion of Chapters 2.4.2, 6 and 7.11.4. 2) The content about Simulated Annealing in Chapter 5 and its application in Chapter 7 have been expanded for clarification. 3) A thorough explanation about the shape-based trajectory model has been included in Chapter 5, to erase the need for lecture of the respective paper by Wall and Conway (2009). 4) A depiction of an example trajectory (as given in Fig. 5-4) in conjunction with Tisserand Graphs has been added, including a discussion, and can be found in Chapter 7.12 Abstract Exploring the solar system is becoming more and more challenging and the associated space missions are becoming similarly more demanding concerning Δv, targeted at distant bodies. Recent examples for such missions are New Horizons, which visited the dwarf planet Pluto and Rosetta’s rendezvous with the comet 67P/Churyumov-Gerasimenko. These missions would not have been possible without the application of gravity assists, a technique already used for Pioneer 10 and 11 and the Voyager missions. It is based on transferring energy between a planet or other flyby partner and a spacecraft, which results in trajectory changes without associated propellant consumption. Low-thrust propulsion, which is usually operating continuously for significant amounts of time during a mission, is very efficient and thus often acts as mission enabler. The effort - expressed in used propellant mass fraction – to achieve a certain mission goal is usually smaller for low-thrust propulsion than for a similar mission applying chemical propulsion. This is due to the large specific impulse (typically >3000s) associated with low-thrust propulsion. To further enhance our ability to explore the solar system it is therefore desirable to combine these two techniques. Currently, a number of methods exists to optimize low-thrust trajectories and these have been applied to gravity-assist scenarios as well. However, up to now a method for optimization of the actual sequence of gravity- assist partners is not publically available. Instead, the respective mission analyst has provided the sequence, based on assumptions, estimates and experience. The major goal of this dissertation is to close this gap and investigate if it is possible to reduce the amount of expert information needed and for optimization and thus increase the degree to which the search space is searched in completeness. Reducing the amount of involved experience and trajectory resp. sequence information, improves the chance of finding unforseen and worthwhile new mission sequences. This dissertation provides an analysis of how to incorporate the gravity- assist sequence into the optimization process and thus allow a complete, thorough and effective search of useful mission scenarios. First, it is analyzed how gravity-assist sequences are obtained for impulsive missions and whether the same technique can be similarly applied to low-thrust scenarios. A typical approach for gravity- assist planning of impulsive mission scenarios is the application of Tisserand’s Criterion. It is an energy relation, derived for observations of comets in the 19th century within the context of the circular restricted three-body system. This criterion can be represented visually in the form of so- called Tisserand graphs, which are used to map possible gravity-assist sequences. Furthermore, two sources of errors are analyzed: First, the error is analysed, which is caused by applying this relation to the actual solar system and thus diverting from the circular restricted three-body system. Second, the error is analysed, which is caused by adding thrust into the equations of motion. It is shown that the thrust force is negligible in the balance of forces for realistic mission designs. However, the thrust’s effect on the specific orbit energy cannot be neglected. Next, a correction term to apply Tisserand’s Criterion in a low-thrust situation is derived. This term however, does not allow an a priori evaluation of the mission sequence and therefore the modified Tisserand’s Criterion cannot be used to map out possible paths of gravity assists similarly as the unmodified Tisserad’s Criteron for impulsive missions. A method, based on a heuristic search, a shape-based trajectory model and the inclusion of possible gravity-assist partners as variable, is set up. The heuristic search circumvents the usage of Tisserand’s Criterion. As a next step to establish the method, the respective variable structure is analyzed. Some of the relevant trajectory variables are shown to be interdepent. For instance, the flight times of the individual trajectory legs between all encountered bodies have to add up to the total mission flight time. Consequently, two types of variables are defined: global and local. The first are independent and thus can be used by an evolutionary algorithm to evolve the mission candidates into better (based on solution fitness) solutions during the search. The latter are interdependent. A number of example calculations are conducted to assess the success and usefulness of applying this search method. It is shown that for a single gravity assist in an Earth to Jupiter mission the obtained results are reliable concerning the gravity-assist partner and date and successfully improve the non-gravity assist benchmark by more than 20%. Performance for multi-gravity-assist missions is worse, but still solution improvement and a dominance of the global variables can be observed. Further calculations introduce constraints based on the variable regions of maximum v change and show that the solution quality and reliability can thus be improved. Furthermore, the limitations of reproducing existing trajectories are shown and discussed. Δ Finally, the drawbacks of the applied shape-based trajectory model are expressed in a discussion concerning its relation to Tisserand graph diagrams. It is recommended to use a different trajectory model to address these drawbacks in future work. Zusammenfassung Die Herausforderung das Sonnensystem zu erforschen wird immer größer und die damit verbundene Planung und Durchführung von Raumfahrtmissionen entsprechend anspruchsvoller. Dies ist auch dadurch bedingt, dass die Ziele solcher Missionen immer schwerer zu erreichen sind. New Horizons, die einen Vorbeiflug am Zwergplaneten Pluto durchgeführt hat und Rosettas Rendezvous mit dem Kometen 67P/Tschurjumow-Gerasimenko sind zwei aktuelle Beispiele für solche anspruchsvollen Missionen. Ohne die Anwendung von sogenannten Schwungholmanövern wären sie nicht möglich gewesen. Diese Methode ist bereits bei den berühmten Missionen Pioneer 10 und 11 angewendet worden und ebenso während des Voyager-Programms. Sie beruht auf Energieaustauch zwischen Planet und Raumfahrzeug. Dieser Energieaustausch ermöglicht Bahnänderungen ohne Treibstoffaufwand und kann so dazu beitragen, dass eine Mission überhaupt durchführbar wird. Sogenannte Niedrigschubtriebwerke, die üblicherweise über lange Zeitabschnitte einer Mission kontinuierlich betrieben werden, stellen eine effiziente Möglichkeit zur Bahnänderung dar. Um eine bestimmte Mission durchführen zu können, ist der Treibstoffaufwand üblicherweise kleiner als für vergleichbare Mission, die mit chemischen, d.h. impulsiven Triebwerken geflogen werden. Dies ergibt sich aus dem hohen spezifischen Impuls (typischerweise > 3000s) von Niedrigschubtriebwerken. Um also die Erforschung des Sonnensystems voranzutreiben, ist es ratsam den Einsatz von Niedrigschubtriebwerken mit Schwungholmanövern zu kombinieren. Zur Optimierung von Niedrigschubbahnen existieren verschiedene Methoden und diese wurden auch bereits für Schwungholszenarien bei solchen Missionen eingesetzt. Allerdings