Global Optimization of Impulsive Interplanetary Transfers
|
|
- Jocelin Oliver
- 6 years ago
- Views:
Transcription
1 Global Optimization of Impulsive Interplanetary Transfers R. Armellin, Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano Taylor Methods and Computer Assisted Proofs Barcelona, June, 3 7, 2008
2 Motivation Space activities are expensive: Ariane 5 launch cost: 200 M$ Allowed Spacecraft Mass: kg = Cost per kilogram: $/kg Propellant represents the main contribution to s/c mass: Propellant is on average 40% of spacecraft mass we want to reduce the required propellant The goal of the trajectory design is to find the best solution in terms of propellant consumption while still achieving the mission goals
3 Outline Dynamical Model Patched-Conics Approximation Two-Impulse Transfers Ephemerides Evaluation Lambert s Problem Solution Differential Algebra Based Global Optimization Rigorous Global Optimization with COSY-GO Multiple Gravity Assist (MGA) space pruning
4 Dynamical Model: 2-Body Problem The 2-Body Problem considers two point masses in mutual orbit about each other m 1 The relative motion of the two masses is governed by: m 2 r = k r 3 r E.g. m1 m2 Sun Spacecraft Analytical solutions exist for the 2-Body Problem: Conic Arcs r = r(θ) explicit t = t(θ) implicit (Kepler s equation)
5 Patched-Conics Approximation The whole interplanetary transfer is divided in several arcs Each arc is the solution of a 2-Body Problem considering the spacecraft and only one other planet at a time E.g.: 2-impulse Earth-Mars transfer 3 conic arcs Earth escape Heliocentric phase Mars capture
6 2-Impulse Planet-to-Planet Transfer 2-impulse Earth-Mars transfer has been selected as first benchmark problem Applied for preliminary design of Earth-Mars (any planet-to-planet transfer) interplanetary transfers Objective function characterized by several comparable local minima Future benchmark problems Multiple Gravity Assist interplanetary transfers E.g.: Cassini-Huygens ] U A [ Y Saturn 01/12/2005 Cassini-Huygens trajectory (Alternative solution) Earth 25/10/1997 Jupiter 17/02/2001 Earth (GA) 0 18/08/1999 Venus (GA2) 24/06/1999 Venus (GA1) 19/05/ X [AU]
7 Optimization Problem The optimization variables are the time of departure t 0 and the time of flight ttof The positions of the starting and arrival planets are computed through the ephemerides evaluation: (re, ve) = eph(t0, Earth) and (rm, vm) = eph(t0 +ttof, Mars) The starting velocity v 1 and the final one v2 are computed by solving the Lambert s problem
8 Optimization Problem The parking velocity and desired final velocities are v c E = µ E /r c E v c M = µ M /r c M The pericenter velocities of the escape and arrival hyperbola v p 1 = 2µ E /r c E + v 1 2 v p 2 = 2µ M /r c M + v 2 2 Objective function V = V 1 + V 2 Constraint: V 1 < V 1,max
9 Ephemerides Evaluation Polynomial interpolations of accurate planetary ephemerides (JPL-Horizon) are used for the preliminary phase of the space trajectory design Given an epoch and a celestial body, its orbital parameters (a, e, i, Ω, ω, M) The nonlinear equation can be analytically evaluated M = E e sin E is solved for the eccentric anomaly E The relation delivers tan θ 2 = 1 + e 1 e tan E 2 θ (Kepler s Eq) The position and the velocity (r, v) of the celestial body in inertial frame reference frame are computed We have to solve an implicit equation: Kepler s equation
10 Lambert s Problem (1/2) Given: initial position r 1 final position r 2 time of flight t tof Find the initial velocity, v1, the spacecraft must have to reach r2 in ttof The solution of the BVP exploits the analytical solution of the 2-body problem Given r1, r2 and ttof there exists only one conic arc connecting the two points in the given time
11 Lambert s Problem (1/2) Several algorithms have been developed for the identification and characterization of the resulting conic arc We used an algorithm developed by Battin (1960) A nonlinear equation must be solved (Lagrange s equation for the time of flight): f(x) = log(a(x)) log(t tof ) = 0 in which A(x) = a(x) 3/2 ((α(x) sin(α(x))) (β(x))), ( ) s c s β(x) = 2 arcsin, and a(x) = 2a(x) 2(1 x 2 ) The value of s and c depend on r 1 and r2, so the nonlinear equation depends both on t0 and ttof
12 DA Solution of Parametric Implicit eqs Search the solution of p [p l, p u ] f(x, p) = 0 for p belonging to Use classical methods (e.g., Newton) to compute x0 solution of f(x, p 0 ) = 0 Initialize [x] = x 0 + x and [p] = p 0 + p as DA variables and expand f = M( x, p) Build the following map and invert it: ( ) f p = ( ) ( ) ( ) [M] x x [I p ] p p Force f = 0 so obtaining the Taylor expansion of of the solution w.r.t. the parameter: = x = x( p) ( ) 1 ( ) [M] f [I p ] p
13 Example: Mars Ephemerides Epoch interval: 40 days (a) Errors on position, (a), and velocity, (b), between the DA and the point-wise evaluation of Mars ephemerides Errors drastically decrease when the order of the Taylor series increases (b)
14 Example: Objective Function The DA evaluation of the planetary ephemerides and the Lambert s problem solution enables the Taylor expansion of the objective function Taylor representation of the objective function Taylor representation error w.r.t. point-wise evaluation Box width: 40 days
15 Earth-Mars Direct Transfer Search space: [1000, 6000] [100, 600] Maximum departure impulse: V1 < 5 km/s Platform: Pentium IV 3.06 GHz laptop Objective function overview 16
16 DA Based Global Optimizer (1/2) DA based global optimization algorithm: Subdivide the search space in subintervals Suitably initialize the value of V opt t tof X For each subinterval : Initialize t0 and ttof as DA variables and compute a Taylor expansion of the objective function V and the constraint V1 on X Bound the value of V 1 on IF X min V 1 > V 1,max Bound the value of V on IF min V > V opt X X discard discard X X t 0
17 DA Based Global Optimizer (2/2) Build and invert the map of the objective function gradient: ( t0 V ttof V Localize the zero-gradient point IF ) x / X = M ( t0 t tof ) discard ( t0 t tof ) x = (t 0, t tof ) X = M 1 ( t0 V ttof V ) Evaluate V = V ( x ) V < V opt V opt x X IF update, and store and If necessary, a more accurate identification of the actual optimum can be finally achieved using a higher order DA computation on the last stored subinterval X x
18 Earth-Mars Direct Transfer Search space: [1000, 6000] [100, 600] Maximum departure impulse: V1 < 5 km/s Platform: Pentium IV 3.06 GHz laptop Solution 1: 10-day boxes + 5th order Pruning + Global Opt: s V opt = km/s x* = [ , ]
19 Earth-Mars Direct Transfer Search space: [1000, 6000] [100, 600] Maximum departure impulse: V1 < 5 km/s Platform: Pentium IV 3.06 GHz laptop Solution 1: 10-day boxes + 5th order Pruning + Global Opt: s V opt = km/s x* = [ , ] Solution 2: 100-day boxes + 5th order Pruning + Global Opt: 0.55 s V opt = km/s x* = [ , ]
20 Verified GO of Earth-Mars Transfer Implicit equations can be solved in a verified way enabling the Taylor Model evaluation of the objective function COSY-GO is applied for the global optimization of an impulsive Earth-Mars transfer Number of steps: Computation time: s Enclosure of the minimum: [ , ] km/s Enclosure of the solution: t 0 [ , ] t tof [ , ]
21 Planet-to-Planet Transfer Number of steps: Computation time: s Enclosure of the minimum: [ , ] km/s Enclosure of the solution: t 0 [ , ] t tof [ , ] Most of the computational time is spent in splitting box containing discontinuities CutOff [km/s] Boxes containing discontinuities with size lower than a given threshold are rejected The solution is mathematically not rigorous Earth-Venus Iteration # x 10 4
22 MGA Space Pruning MGA transfers problems are characterized by High number of variables (+1 for each GA) Constraints on gravity assist radii and VGA Increased multimodality and more discontinuities Jupiter, T3, V3 Earth, T1, V1 Mars, T2, VGA Validated solution is almost impractical (we have obtained solution for 1 GA) The objective function is simply given by V = V1 + VGAi + Vn IDEA: use DA expansion and polynomial bounders just to prune the search space
23 DA Based Pruning of MGA For the first arc: i = 1 Subdivide T 1 and T 2 in subintervals For each subinterval X : X - Compute the expansion of V 1 on X - Bound the polynomial expansion of V 1 on X - IF min V 1 > V 1,max discard X and all its subsequent combinations For each gravity assist: i = 2,..., n 1 Subdivide T i+1 in subintervals For each subinterval remaining from the previous arc and for each combination with the subinterval on T i+1, X - Compute the expansion of V GAi and r pi on X - Bound the polynomial expansion of and on X V GAi r pi
24 DA Based Pruning of MGA - IF min V GAi > V GAi,max or discard X and all its subsequent combinations max r pi < r pi,min For the last arc: i = n For each subinterval remaining from the previous arc V n - Compute the expansion of on - Bound the polynomial expansion of on min V n > V n,max - IF discard X V n X X
25 Cassini-like transfer Problem Definition: Planet E Search Space [day] Int. size [day] V Cut-off [km/s] Dep. Epoch [MJD2000] V tof V tof E tof J tof S tof
26 Cassini-like transfer Tot. boxes Remaining 1085 (0.003%) Obj Fcn km/s This work has been done under ESA Ariadna contract Final report available at ariadna/completed.htm
27 Global Optimization of Impulsive Interplanetary Transfers R. Armellin, Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano Taylor Methods and Computer Assisted Proofs Barcelona, June, 3 7, 2008
28 Conclusions and Future Work DA and TM global optimizers are effective tools for the global optimization of planet-to-planet transfers Efficient management of regions with singularities is needed for TM global optimization with COSY-GO Validated management of nonlinear constraints will be required to apply COSY-GO to MGA transfers DA is a promising technique for search space pruning of high dimensional problems such Multiple Gravity Assist (MGA) interplanetary transfers
29 Verified Implicit Eq Solution - 1D Suppose to have the (n +1) differentiable function f over the domain D = [-1, 1] and its n-th order Taylor model P(x) + I so that f(x) P (x) + I for all x D Consider the enclosure R of P(x) + I over D and suppose P (x) >d >0 on D with P(0) = 0 Find the Taylor Model C(y) + J on R so that any solution of the problem f(x) = y lies in C(y) + J Algorithm: First compute C(y), the n-th order polynomial inversion of P(x), so that P (C(y)) = n y Using Taylor model computation, obtain where P (C(y)) y + J includes the terms of order exceeding n in P(C(y)), and thus scales with at least order n +1 J
30 Verified Implicit Eq Solution - 1D Use the consequences of small correction x to C(y) to find the rigorous reminder J for C(y) so that all the solutions of f (x) = y lie in C(y) +J. According to the mean value theorem: f(c(y) + x) y P (C(y) + x) y + I for suitable ξ [C(y), C(y) + x] Since P is bounded below by d, the set x P (ξ) + I + J will never contain the zero except for the interval J = I + J which is the desired interval d = P (C(y)) + x P (ξ) y + I y + J + x P (ξ) y + I = x P (ξ) + I + J
31 Verified Implicit Eq Solution - vd Let P(x) + I be a n-th order Taylor model of the (n +1) times differentiable function f over the domain f (x) P(x) + I for all x D Indicate with L(x) the linear part of P(x) D = [ 1, 1] v so that: Instead of the original problem and in analogy with the 1D case, consider the problem of finding a verified enclosure of the inverse of L 1 f where L is analytically inverted P(x) + J The Taylor model enclosure of L 1 f over D is: P + J = L 1 (P + I)
32 Verified Implicit Eq Solution - vd It is worth observing that: P 1 = x 1 + h.o.t. we can bound P 1 from below x 1 P 2 = x 2 + h.o.t. we can bound P 2 from below x 2 etc. Consequently we can proceed as in the 1D case on When the solution has been obtained for right-compose with L 1 Application to the solution of f (x, p) = 0: { y = f (x, p) p = p L 1 f L 1 f Once a validated inversion of the system is achieved, just set y = 0
33 Orbital parameters The orbital parameters are: (a, e, i, Ω, ω, θ) The position and the velocity (r, v) in cartesian coordinates are obtained from the orbital parameters by simple algebraic relations
Rigorous Global Optimization of Impulsive Space Trajectories
Rigorous Global Optimization of Impulsive Space Trajectories P. Di Lizia, R. Armellin, M. Lavagna K. Makino, M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16 19, 2006 Motivation
More informationGlobal Trajectory Optimisation: Can We Prune the Solution Space when Considering Deep Space Maneuvers?
Global Trajectory Optimisation: Can We Prune the Solution Space when Considering Deep Space Maneuvers? Final Report Authors: F. Bernelli-Zazzera, M. Lavagna, R. Armellin, P. Di Lizia, F. Topputo Affiliation:
More informationVerified High-Order Optimal Control in Space Flight Dynamics
Verified High-Order Optimal Control in Space Flight Dynamics R. Armellin, P. Di Lizia, F. Bernelli-Zazzera K. Makino and M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) L06: Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 Problem Statement? Hint #1: design the Earth-Mars transfer using known concepts
More informationTwo-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra
Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Politecnico di Milano Department of Aerospace Engineering Milan, Italy Taylor Methods and Computer Assisted Proofs
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 6. Interplanetary Trajectories 6.1 Patched conic method 6.2 Lambert s problem
More informationFlight and Orbital Mechanics
Flight and Orbital Mechanics Lecture slides Challenge the future 1 Flight and Orbital Mechanics AE-104, lecture hours 1-4: Interplanetary flight Ron Noomen October 5, 01 AE104 Flight and Orbital Mechanics
More informationRelevant Applications of Differential Algebra in Astrodynamics
AstroNet-II International Final Conference Tossa de Mar 15-19 June 2015 Relevant Applications of Differential Algebra in Astrodynamics Pierluigi Di Lizia in collaboration with: Roberto Armellin, Alessandro
More informationThe B-Plane Interplanetary Mission Design
The B-Plane Interplanetary Mission Design Collin Bezrouk 2/11/2015 2/11/2015 1 Contents 1. Motivation for B-Plane Targeting 2. Deriving the B-Plane 3. Deriving Targetable B-Plane Elements 4. How to Target
More informationApophis Encounter 2029: Differential Algebra and Taylor Model Approaches. 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids
Apophis Encounter 9: Differential Algebra and Taylor Model Approaches 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids 7- April 9 Granada, Spain R. Armellin (1), P. Di Lizia (1),
More informationInterplanetary Mission Opportunities
Interplanetary Mission Opportunities Introduction The quest for unravelling the mysteries of the universe is as old as human history. With the advent of new space technologies, exploration of space became
More informationINTERPLANETARY AND LUNAR TRANSFERS USING LIBRATION POINTS
INTERPLANETARY AND LUNAR TRANSFERS USING LIBRATION POINTS Francesco Topputo (), Massimiliano Vasile () and Franco Bernelli-Zazzera () () Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano,
More informationEarth-Mars Halo to Halo Low Thrust
Earth-Mars Halo to Halo Low Thrust Manifold Transfers P. Pergola, C. Casaregola, K. Geurts, M. Andrenucci New Trends in Astrodynamics and Applications V 3 June / -2 July, 28 Milan, Italy Outline o Introduction
More informationValidated integration of Solar System Dynamics. 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids
Validated integration of Solar System Dynamics 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids 27-3 April 29 Granada, Spain P. Di Lizia (1), R. Armellin (1), F. Bernelli Zazzera
More informationASEN 5050 SPACEFLIGHT DYNAMICS Interplanetary
ASEN 5050 SPACEFLIGHT DYNAMICS Interplanetary Prof. Jeffrey S. Parker University of Colorado Boulder Lecture 29: Interplanetary 1 HW 8 is out Due Wednesday, Nov 12. J2 effect Using VOPs Announcements Reading:
More informationMassimiliano Vasile, Stefano Campagnola, Paolo Depascale, Stefano Pessina, Francesco Topputo
A Toolbox for Preliminary Massimiliano Vasile, Stefano Campagnola, Paolo Depascale, Stefano Pessina, Francesco Topputo Mission Analysis and Design PAMSIT IMAGO ATOM-C EPIC Massimiliano Vasile, Stefano
More informationInterplanetary Travel
Interplanetary Travel Interplanetary Travel Concept Patched Conic Hypothesis Departure & Arrival Manoeuvres Interplanetary Travel Concept Interplanetary travel is concerned with motion of manmade objects
More informationBenchmarking different global optimisation techniques for preliminary space trajectory design
58th International Astronautical Congress, Hyderabad, 24-28 September 27. Copyright 27 by Advanced Concepts Team. Published by the IAF, with permission and Released to IAF or IAA to publish in all forms.
More informationInterplanetary Preliminary Mission Design: GA and AGA sequences optimisation
Interplanetary Preliminary Mission Design: GA and AGA sequences optimisation Michelle LAVAGNA Amalia ERCOLI FINZI Angelo POVOLERI 1 Tool purpose Space mission design Mission analysis System design (concurrent
More information10 Orbit and Constellation Design Selecting the Right Orbit
Orbit and Constellation Design Selecting the Right Orbit.7 Design of Interplanetary Orbits Faster Trajectories Ron Noomen, Delft University of Technology Using the recipe given in Table -9, one can compute
More informationL eaving Earth and arriving at another planet or asteroid requires
Designing Interplanetary Transfers L eaving Earth and arriving at another planet or asteroid requires a spacecraft to implement a sequence of manoeuvres. These include changes of velocity needed to escape
More informationInterplanetary Mission Analysis
Interplanetary Mission Analysis Stephen Kemble Senior Expert EADS Astrium stephen.kemble@astrium.eads.net Page 1 Contents 1. Conventional mission design. Advanced mission design options Page 1. Conventional
More informationCelestial Mechanics I. Introduction Kepler s Laws
Celestial Mechanics I Introduction Kepler s Laws Goals of the Course The student will be able to provide a detailed account of fundamental celestial mechanics The student will learn to perform detailed
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationPatched Conic Interplanetary Trajectory Design Tool
Copyright by Martin James Brennan 2011 The Thesis committee for Martin James Brennan Certifies that this is the approved version of the following thesis: Patched Conic Interplanetary Trajectory Design
More informationOverview of Astronautics and Space Missions
Overview of Astronautics and Space Missions Prof. Richard Wirz Slide 1 Astronautics Definition: The science and technology of space flight Includes: Orbital Mechanics Often considered a subset of Celestial
More informationTrajectory Optimisation under Uncertainties
Trajectory Optimisation under Uncertainties Final Report Authors: F. Bernelli-Zazzera, M. Lavagna, R. Armellin, P. Di Lizia, A. Morselli Affiliation: Department of Aerospace Engineering, Politecnico di
More informationGPU Accelerated Genetic Algorithm for Multiple Gravity-Assist and Impulsive V Maneuvers
GPU Accelerated Genetic Algorithm for Multiple Gravity-Assist and Impulsive V Maneuvers Sam Wagner, Brian Kaplinger, and Bong Wie Asteroid Deflection Research Center, Iowa State University, Ames, IA, USA
More information1 st ACT Competition on Global Trajectory Optimisation Solution of Team #7
1 st ACT Competition on Global Trajectory Optimisation Solution of Team #7 Régis Bertrand Thierry Ceolin Richard Epenoy CNES-DCT/SB/MO CS-ESPACE/FDS CNES-DCT/SB/MO Contents Presentation of team #7 Available
More informationResults found by the CNES team (team #4)
3 rd Global Trajectory Optimisation Competition (GTOC3) organized by the Aerospace Propulsion Group of the Dipartimento di Energetica at Politecnico di Torino Results found by the CNES team (team #4) Presented
More informationPatch Conics. Basic Approach
Patch Conics Basic Approach Inside the sphere of influence: Planet is the perturbing body Outside the sphere of influence: Sun is the perturbing body (no extra-solar system trajectories in this class...)
More information5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References...
1 The Two-Body Problem... 1 1.1 Position of the Problem... 1 1.2 The Conic Sections and Their Geometrical Properties... 12 1.3 The Elliptic Orbits... 20 1.4 The Hyperbolic and Parabolic Trajectories...
More informationProblem Description for the 6 th Global Trajectory Optimisation Competition. Anastassios E. Petropoulos 1
Problem Description for the 6 th Global Trajectory Optimisation Competition Anastassios E. Petropoulos 1 Outer Planet Mission Analysis Group Jet Propulsion Laboratory California Institute of Technology
More informationASTRIUM. Interplanetary Path Early Design Tools at ASTRIUM Space Transportation. Nathalie DELATTRE ASTRIUM Space Transportation.
Interplanetary Path Early Design Tools at Space Transportation Nathalie DELATTRE Space Transportation Page 1 Interplanetary missions Prime approach: -ST has developed tools for all phases Launch from Earth
More informationFrom Particle Accelerators to Celestial Dynamics: An Introduction to Differential Algebra
AstroNet II Summer School 2014 Friday, 27 June 2014 Zielona Góra (Poland) From Particle Accelerators to Celestial Dynamics: An Introduction to Differential Algebra 1, Pierluigi di Lizia 1, Roberto Armellin
More informationLecture D30 - Orbit Transfers
J. Peraire 16.07 Dynamics Fall 004 Version 1.1 Lecture D30 - Orbit Transfers In this lecture, we will consider how to transfer from one orbit, or trajectory, to another. One of the assumptions that we
More informationThe optimization of electric propulsion (EP) trajectories for interplanetary missions is a quite difficult
Indirect Optimization Method for Low-Thrust Interplanetary Trajectories IEPC-27-356 Presented at the 3 th International Electric Propulsion Conference, Florence, Italy Lorenzo Casalino, Guido Colasurdo
More informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationABOUT COMBINING TISSERAND GRAPH GRAVITY-ASSIST SEQUENCING WITH LOW-THRUST TRAJECTORY OPTIMIZATION
ABOUT COMBINING TISSERAND GRAPH GRAVITY-ASSIST SEQUENCING WITH LOW-THRUST TRAJECTORY OPTIMIZATION Volker Maiwald German Aerospace Center (DLR) Institute of Space Systems Department of System Analysis Space
More informationFrom the Earth to the Moon: the weak stability boundary and invariant manifolds -
From the Earth to the Moon: the weak stability boundary and invariant manifolds - Priscilla A. Sousa Silva MAiA-UB - - - Seminari Informal de Matemàtiques de Barcelona 05-06-2012 P.A. Sousa Silva (MAiA-UB)
More informationTrajectories to the Moons (incl. a trajectory for an Enceladus orbiter)
Trajectories to the Moons (incl. a trajectory for an Enceladus orbiter) Stefano Campagnola 1,Ryan Russell 2, and Nathan Strange 3 JPL-Caltech/CDS Meetings, August 3rd 2009 1 Aerospace and Mechanical Engineering,
More informationORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)
ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) PROBLEM 1 A) Summarize the content of the three Kepler s Laws. B) Derive any two of the Kepler
More informationResearch Article Analysis of Electric Propulsion System for Exploration of Saturn
Mathematical Problems in Engineering Volume 29, Article ID 75637, 4 pages doi:.55/29/75637 Research Article Analysis of Electric Propulsion System for Exploration of Saturn Carlos Renato Huaura Solórzano,
More informationFinal Rankings and Brief Descriptions of the Returned Solutions and Methods Used for the 2 nd Global Trajectory Optimisation Competition
Final Rankings and Brief Descriptions of the Returned Solutions and Methods Used for the nd Global Trajectory Optimisation Competition Anastassios E. Petropoulos Outer Planets Mission Analysis Group Jet
More informationORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM
ORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM JORGE K. S. FORMIGA 1,2 and ANTONIO F B A PRADO 2 National Institute for Space Research -INPE 1 Technology Faculty-FATEC-SJC
More informationLOTNAV. A Low-Thrust Interplanetary Navigation Tool: The Trajectory Reconstruction Module
LOTNAV A Low-Thrust Interplanetary Navigation Tool: The Trajectory Reconstruction Module Juan Luis Cano González Mission Analysis Division Deimos Space S.L. -1- The LOTNAV Tool The Low-Thrust Interplanetary
More informationAA 528 Spacecraft Dynamics and Control. Mehran Mesbahi Aeronautics & Astronautics Winter 2017 University of Washington
AA 528 Spacecraft Dynamics and Control Mehran Mesbahi Aeronautics & Astronautics Winter 2017 University of Washington Spacecraft dynamics and control What is this class all about? what is in the name?
More informationDEFLECTING HAZARDOUS ASTEROIDS FROM COLLISION WITH THE EARTH BY USING SMALL ASTEROIDS
DEFLECTING HAZARDOUS ASTEROIDS FROM COLLISION WITH THE EARTH BY USING SMALL ASTEROIDS N. Eismont (1), M. Boyarsky (1), A. Ledkov (1), B.Shustov (2), R. Nazirov (1), D. Dunham (3) and K. Fedyaev (1) (1)
More informationThird Body Perturbation
Third Body Perturbation p. 1/30 Third Body Perturbation Modeling the Space Environment Manuel Ruiz Delgado European Masters in Aeronautics and Space E.T.S.I. Aeronáuticos Universidad Politécnica de Madrid
More informationISAS MERCURY ORBITER MISSION TRAJECTORY DESIGN STRATEGY. Hiroshi Yamakawa
ISAS MERCURY ORBITER MISSION TRAJECTORY DESIGN STRATEGY Hiroshi Yamakawa Institute of Space and Astronautical Science (ISAS) 3-1-1 Yoshinodai, Sagamihara, Kanagawa, 229-851 Japan E-mail:yamakawa@pub.isas.ac.jp
More informationInterplanetary Trajectory design for Rosetta and Solar Orbiter
Jose Rodriguez-Canabal Memorial Interplanetary Trajectory design for Rosetta and Solar Orbiter Johannes Schoenmaekers H / Mission Analysis Section (HSO-GFA) Senior Adviser ESA / ESOC Email: johannes.schoenmaekers@esa.int
More informationState Vector Reference Body Transformations Throughout The Solar System
Introduction One of the more mind-boggling programming challenges encountered during interplanetary trajectory design software development is how to transform a specified inertial position and velocity
More informationChapter 4. Integrated Algorithm for Impact Lunar Transfer Trajectories. using Pseudo state Technique
Chapter 4 Integrated Algorithm for Impact Lunar Transfer Trajectories using Pseudo state Technique Abstract A new integrated algorithm to generate the design of one-way direct transfer trajectory for moon
More informationMASTER'S THESIS. Multi-gravity assist design tool for interplanetary trajectory optimisation
MASTER'S THESIS 2009:062 Multi-gravity assist design tool for interplanetary trajectory optimisation Luleå University of Technology Master Thesis, Continuation Courses Advanced material Science and Engineering
More informationDEVELOPMENT OF A MULTIPURPOSE LOW THRUST INTERPLANETARY TRAJECTORY CALCULATION CODE
AAS 03-667 DEVELOPMENT OF A MULTIPURPOSE LOW THRUST INTERPLANETARY TRAJECTORY CALCULATION CODE Tadashi Sakai 1 and John R. Olds 2 A multipurpose low thrust interplanetary trajectory calculation code has
More informationSolar Orbiter Ballistic Transfer Mission Analysis Synthesis
European Space Agency Agence Spatiale Européenne directorate of operations and infrastructure ground systems engineering department mission analysis office MAO Working Paper No. 483 Issue 1, Rev. 0 Solar
More informationDesign of a Multi-Moon Orbiter
C C Dynamical A L T E C S H Design of a Multi-Moon Orbiter Shane D. Ross Control and Dynamical Systems and JPL, Caltech W.S. Koon, M.W. Lo, J.E. Marsden AAS/AIAA Space Flight Mechanics Meeting Ponce, Puerto
More informationarxiv: v1 [physics.pop-ph] 29 May 2009
Gravity assist in 3D like in Ulysses mission Herbert Morales Escuela de Física, Universidad de Costa Rica, San José, Costa Rica arxiv:0905.4788v1 [physics.pop-ph] 29 May 2009 (Dated: October 22, 2018)
More informationEffect of Coordinate Switching on Translunar Trajectory Simulation Accuracy
Effect of Coordinate Switching on Translunar Trajectory Simulation Accuracy Mana P. Vautier Auburn University, Auburn, AL, 36849, USA This paper focuses on the affect of round-off error in the accurate
More informationWhat is the InterPlanetary Superhighway?
What is the InterPlanetary Superhighway? Kathleen Howell Purdue University Lo and Ross Trajectory Key Space Technology Mission-Enabling Technology Not All Technology is hardware! The InterPlanetary Superhighway
More informationEscape Trajectories from Sun Earth Distant Retrograde Orbits
Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists30, pp. Pd_67-Pd_75, 06 Escape Trajectories from Sun Earth Distant Retrograde Orbits By Yusue OKI ) and Junichiro KAWAGUCHI ) ) Department of Aeronautics
More informationEnd-Of-Life Disposal Concepts for Lagrange-Point and Highly Elliptical Orbit Missions
End-Of-Life Disposal Concepts for Lagrange-Point and Highly Elliptical Orbit Missions Executive summary of the main study and the study extension Version 1.0 12 February 2015 ESA/ESOC contract No. 4000107624/13/F/MOS
More informationSpacecraft Dynamics and Control
Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 5: Hyperbolic Orbits Introduction In this Lecture, you will learn: Hyperbolic orbits Hyperbolic Anomaly Kepler s Equation,
More informationAstrodynamics of Interplanetary Cubesats
Astrodynamics of Interplanetary Cubesats F. Topputo Politecnico di Milano, Italy 06 F. Topputo, Politecnico di Milano. All rights reserved. icubesat 06, 5th Interplanetary CubeSat Workshop 4-5 May 06,
More informationSUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS. Sandro da Silva Fernandes. Cleverson Maranhão Porto Marinho
SUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS Sandro da Silva Fernandes Instituto Tecnológico de Aeronáutica, São José dos Campos - 12228-900 - SP-Brazil, (+55) (12) 3947-5953 sandro@ita.br Cleverson
More informationProblem Description for the 8 th Global Trajectory Optimisation Competition. Anastassios E. Petropoulos 1
1 Background Problem Description for the 8 th Global Trajectory Optimisation Competition Anastassios E. Petropoulos 1 Outer Planet Mission Analysis Group Mission Design and Navigation Section Jet Propulsion
More informationINDIRECT PLANETARY CAPTURE VIA PERIODIC ORBIT ABOUT LIBRATION POINTS
6 th International Conference on Astrodynamics Tools and Technique (ICATT) INDIRECT PLANETARY CAPTURE VIA PERIODIC ORBIT LI Xiangyu, Qiao Dong, Cui Pingyuan Beijing Institute of Technology Institute of
More informationOrbital Dynamics and Impact Probability Analysis
Orbital Dynamics and Impact Probability Analysis (ISAS/JAXA) 1 Overview This presentation mainly focuses on a following point regarding planetary protection. - How to prove that a mission satisfies the
More informationSome Methods for Global Trajectory Optimisation
Some Methods for Global Trajectory Optimisation used in the First ACT Competition on Global Trajectory Optimisation European Space Agency Team 11: Jet Propulsion Laboratory California Institute of Technology
More informationD R A F T. P. G. Hänninen 1 M. Lavagna 1 Dipartimento di Ingegneria Aerospaziale
Multi-Disciplinary Optimisation for space vehicles during Aero-assisted manoeuvres: an Evolutionary Algorithm approach P. G. Hänninen 1 M. Lavagna 1 p.g.hanninen@email.it, lavagna@aero.polimi.it 1 POLITECNICO
More informationA Study of the Close Approach Between a Planet and a Cloud of Particles
A Study of the Close Approach Between a Planet a Cloud of Particles IIAN MARTINS GOMES, ANTONIO F. B. A. PRADO National Institute for Space Research INPE - DMC Av. Dos Astronautas 1758 São José dos Campos
More informationA Simple Semi-Analytic Model for Optimum Specific Impulse Interplanetary Low Thrust Trajectories
A Simple Semi-Analytic Model for Optimum Specific Impulse Interplanetary Low Thrust Trajectories IEPC-2011-010 * Presented at the 32nd International Electric Propulsion Conference, Wiesbaden Germany David
More informationList of Tables. Table 3.1 Determination efficiency for circular orbits - Sample problem 1 41
List of Tables Table 3.1 Determination efficiency for circular orbits - Sample problem 1 41 Table 3.2 Determination efficiency for elliptical orbits Sample problem 2 42 Table 3.3 Determination efficiency
More informationA space probe to Jupiter
Problem 3 Page 1 Problem 3 A space probe to Jupiter We consider in this problem a method frequently used to accelerate space probes in the desired direction. The space probe flies by a planet, and can
More informationConic Sections in Polar Coordinates
Conic Sections in Polar Coordinates MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Introduction We have develop the familiar formulas for the parabola, ellipse, and hyperbola
More informationCHAPTER 8 PLANETARY MOTIONS
1 CHAPTER 8 PLANETARY MOTIONS 8.1 Introduction The word planet means wanderer (πλάνητες αστέρες wandering stars); in contrast to the fixed stars, the planets wander around on the celestial sphere, sometimes
More informationOrbit Characteristics
Orbit Characteristics We have shown that the in the two body problem, the orbit of the satellite about the primary (or vice-versa) is a conic section, with the primary located at the focus of the conic
More informationChapter 8. Precise Lunar Gravity Assist Trajectories. to Geo-stationary Orbits
Chapter 8 Precise Lunar Gravity Assist Trajectories to Geo-stationary Orbits Abstract A numerical search technique for designing a trajectory that transfers a spacecraft from a high inclination Earth orbit
More informationMAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30.
MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30. Guidelines: Please turn in a neat and clean homework that gives all the formulae that you have used as well as details that
More informationAnalysis of V, Leveraging for Interplanetary Missions
Analysis of V, Leveraging for Interplanetary Missions Jon A. ~irns* and James M. Longuskit Purdue University, West Lufayette, Indiana 47907-1282 V, leveraging can significantly reduce the launch energy
More informationOptimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission
Optimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission Deepak Gaur 1, M. S. Prasad 2 1 M. Tech. (Avionics), Amity Institute of Space Science and Technology, Amity University, Noida, U.P.,
More informationFeedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field
Vol. 4, No. 4, 23 Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Ashraf H. Owis Department of Astronomy, Space and Meteorology, Faculty of Science, Cairo University Department
More information1 Summary of Chapter 2
General Astronomy (9:61) Fall 01 Lecture 7 Notes, September 10, 01 1 Summary of Chapter There are a number of items from Chapter that you should be sure to understand. 1.1 Terminology A number of technical
More informationLesson 11: Orbital Transfers II. 10/6/2016 Robin Wordsworth ES 160: Space Science and Engineering: Theory and ApplicaCons
Lesson 11: Orbital Transfers II 10/6/2016 Robin Wordsworth ES 160: Space Science and Engineering: Theory and ApplicaCons ObjecCves Introduce concept of sphere of influence Study the patched conics approach
More informationInterplanetary Trajectory Optimization using a Genetic Algorithm
Interplanetary Trajectory Optimization using a Genetic Algorithm Abby Weeks Aerospace Engineering Dept Pennsylvania State University State College, PA 16801 Abstract Minimizing the cost of a space mission
More informationSpiral Trajectories in Global Optimisation of Interplanetary and Orbital Transfers
Spiral Trajectories in Global Optimisation of Interplanetary and Orbital Transfers Final Report Authors: M. Vasile, O. Schütze, O. Junge, G. Radice, M. Dellnitz Affiliation: University of Glasgow, Department
More informationEarth-Mars transfers with ballistic escape and low-thrust capture
Earth-Mars transfers with ballistic escape and low-thrust capture G. Mingotti, F. Topputo, F. Bernelli-Zazzera To cite this version: G. Mingotti, F. Topputo, F. Bernelli-Zazzera. Earth-Mars transfers with
More informationMassimiliano Vasile. The Netherlands. Abstract
A GLOBAL OPTIMIZATION ALGORITHM FOR SPACE TRAJECTORY DESIGN Massimiliano Vasile European Space & Technology Centre ESA/ESTEC, SERA-A, Keplerlaan 1 2201 AZ Noordwijk, The Netherlands Massimiliano.Vasile@esa.int
More informationGravitation and the Waltz of the Planets
Gravitation and the Waltz of the Planets Chapter Four Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets
More informationGravitation and the Waltz of the Planets. Chapter Four
Gravitation and the Waltz of the Planets Chapter Four Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets
More informationMethod for Rapid Interplanetary Trajectory Analysis using V Maps with Flyby Options
Method for Rapid Interplanetary Trajectory Analysis using V Maps with Flyby Options The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationSURVEY OF GLOBAL OPTIMIZATION METHODS FOR LOW- THRUST, MULTIPLE ASTEROID TOUR MISSIONS
AAS 07-211 SURVEY OF GLOBAL OPTIMIZATION METHODS FOR LOW- THRUST, MULTIPLE ASTEROID TOUR MISSIONS INTRODUCTION Kristina Alemany *, Robert D. Braun Electric propulsion has recently become a viable option
More informationUlrich Walter. Astronautics. The Physics of Space Flight. 2nd, Enlarged and Improved Edition
Ulrich Walter Astronautics The Physics of Space Flight 2nd, Enlarged and Improved Edition Preface to Second Edition Preface XVII Acknowledgments XIX List of Symbols XXI XV 1 Rocket Fundamentals 1 1.1 Rocket
More informationGravitation and the Motion of the Planets
Gravitation and the Motion of the Planets 1 Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets go around
More informationThe Heliocentric Model of Copernicus
Celestial Mechanics The Heliocentric Model of Copernicus Sun at the center and planets (including Earth) orbiting along circles. inferior planets - planets closer to Sun than Earth - Mercury, Venus superior
More informationResearch Publications at Politecnico di Milano
RE.PUBLIC@POLIMI Research Publications at Politecnico di Milano This is the published version of: D. Menzio, C. Colombo The Combined Lambert-Tisserand Method Applied to the Single Flyby Problem Paper presented
More informationSwing-By Maneuvers for a Cloud of Particles with Planets of the Solar System
Swing-By Maneuvers for a Cloud of Particles with Planets of the Solar System VIVIAN MARTINS GOMES, ANTONIO F. B. A. PRADO National Institute for Space Research INPE - DMC Av. Dos Astronautas 1758 São José
More informationMission Trajectory Design to a Nearby Asteroid
Mission Trajectory Design to a Nearby Asteroid A project present to The Faculty of the Department of Aerospace Engineering San Jose State University in partial fulfillment of the requirements for the degree
More informationAbstract A COUPLED INTERPLANETARY ENTRY, DESCENT AND LANDING TARGETING PROCEDURE. Jeremy David Shidner, M.S., National Institute of Aerospace
Abstract Title of Thesis: A COUPLED INTERPLANETARY ENTRY, DESCENT AND LANDING TARGETING PROCEDURE Jeremy David Shidner, M.S., 006 Thesis Directed by: Dr. Robert Tolson National Institute of Aerospace In
More informationMission Overview. EAGLE: Study Goals. EAGLE: Science Goals. Mission Architecture Overview
Mission Overview OPAG Meeting November 8 th, 2006 Ryan Anderson & Daniel Calvo EAGLE: Study Goals Develop a set of science goals for a flagship mission to Enceladus Investigate mission architectures that
More information