Design of a Planetary Explorer
Third Year Project
Jack Andrews (ChCh) Aamir Aziz (New) Nicholas Baker (SJC) Alexander Chadwick (SJC) James Coates (BNC) Oliver Cohen (LMH) Liam Donovan (Hert) Hugo Grimmett (BNC) James Hawkes (Keb) Michael May (Worc) Joshua McFarlane (SJC)
2009-2010 Contents
1 Introduction 13
1.1 Mission to Titan ...... 13
1.2 Titan’s Characteristics ...... 13
1.3 Mission Components ...... 14
2 Delivery of the Orbiter 15
2.1 Introduction ...... 15
2.2 Orbital Mechanics ...... 15
2.2.1 Kepler’s Laws ...... 15
2.2.2 Energy Considerations ...... 17
2.2.3 Characterising Orbits ...... 18
2.3 Transfer Orbits ...... 20
2.4 The N-Body Problem ...... 20
2.4.1 A Computational Model ...... 20
2.4.2 Sourcing Planetary Data ...... 21
2.5 Simulating the Problem in Matlab ...... 22
2.6 Fly-by Manoeuvres ...... 22
2.6.1 The Origin of Gravity-Assist Velocity Increase ...... 22
1 CONTENTS
2.6.2 Calculating the Velocity Increase ...... 23
2.7 The Oberth Effect ...... 23
2.8 Time-to-Intercept & Mean Anomaly ...... 23
2.9 Targetting Manoeuvres ...... 24
2.10 Planning out the route to Titan ...... 26
2.11 Launch ...... 26
2.11.1 Calculating the Relative Launch Velocity ...... 28
2.11.2 Optimal Launch Window ...... 29
2.12 The Control System ...... 31
2.12.1 Implementing the Control System in Matlab ...... 31
2.13 The Simulation ...... 32
2.14 Conclusions ...... 35
3 Attitude Control System 36
3.1 Introduction ...... 36
3.1.1 Purpose ...... 36
3.1.2 Control Goal ...... 37
3.1.3 System Overview ...... 37
3.1.4 Aim ...... 38
3.2 Mathematical model of the orbiter dynamics ...... 38
3.2.1 Reference Frames ...... 38
3.2.2 Rigid body Dynamics ...... 39
3.2.3 Attitude Kinematics ...... 40
3.2.4 Actuator Dynamics ...... 40
3.2.5 Gravity Gradient Torque ...... 41
2 CONTENTS
3.2.6 Linearisation of the model ...... 42
3.3 Continuous-time Controller ...... 43
3.3.1 Plant Model ...... 43
3.3.2 Design Specifications ...... 44
3.3.3 PID Control ...... 45
3.3.4 Optimal Control ...... 47
3.4 Discrete-time Controller ...... 50
3.4.1 Discrete Plant Model ...... 50
3.4.2 Design by Emulation ...... 51
3.4.3 Estimator Design ...... 51
3.4.4 Regulator Design ...... 53
3.5 Conclusions ...... 54
4 Landing the Probe 56
4.1 Introduction ...... 56
4.2 Gathering Data ...... 56
4.2.1 Wind Speed ...... 57
4.2.2 Pressure and Density ...... 59
4.3 Simulation Engine ...... 60
4.3.1 Gravity ...... 60
4.3.2 Drag ...... 60
4.3.3 First Simulation ...... 61
4.3.4 Drift ...... 62
4.4 A Safe Landing ...... 62
4.4.1 Parachutes/Thrusters ...... 62
3 CONTENTS
4.4.2 Simulating Parachutes ...... 63
4.4.3 Test Data ...... 63
4.5 Analysis ...... 66
4.5.1 Test Plan ...... 66
4.5.2 The Need For Analysis ...... 67
4.5.3 Analysis Example - Drift ...... 68
4.5.4 Parameter Weighting ...... 70
4.5.5 Result Weighting ...... 72
4.5.6 Finding the Optimal Area ...... 72
4.6 Final Design ...... 75
4.6.1 Results ...... 76
4.7 Errors ...... 78
4.7.1 Systematic Errors ...... 78
4.7.2 Theoretical Wind Model ...... 80
4.8 Predicted Landing Zone ...... 81
5 Problems Faced by the Explorer During its Voyage Through Space 83
5.1 Introduction ...... 83
5.2 Problems ...... 83
5.2.1 Problems Encountered on Earth ...... 84
5.2.2 Problems Encountered whilst en Route to Titan ...... 84
5.2.3 Problems Encountered whilst in Orbit around Titan ...... 85
5.3 Effects ...... 87
5.3.1 Heating Effects ...... 87
5.3.2 Radiation Effects ...... 88
4 CONTENTS
5.3.3 Power Requirement Effects ...... 90
5.4 Solutions ...... 91
5.4.1 Solving the Heating Problem ...... 91
5.4.2 Solving the Radiation Problem ...... 94
5.4.3 Solving the Electrostatic Discharge Problem ...... 97
5.4.4 Solving the Power Demand Problem ...... 99
5.4.5 Solving the Environmental Hazard Problem ...... 101
5.5 Conclusion ...... 103
6 Lander Power Systems 104
6.1 Introduction ...... 104
6.1.1 Lander Power Specifications ...... 104
6.1.2 Design Parameters ...... 105
6.2 Power System Options ...... 105
6.2.1 Solar ...... 105
6.2.2 Radioisotope Thermoelectric Generators ...... 106
6.2.3 Wind ...... 107
6.2.4 Tidal ...... 108
6.2.5 Batteries ...... 108
6.2.6 Selection ...... 109
6.3 Wind Power ...... 110
6.3.1 Feasibility ...... 110
6.3.2 Implementation ...... 112
6.3.3 Deployment and Operation ...... 115
6.3.4 Voltage rectification ...... 116
5 CONTENTS
6.3.5 Conclusions ...... 118
6.4 Secondary Batteries ...... 119
6.4.1 Chemistry ...... 119
6.4.2 Transit to Titan ...... 122
6.4.3 Charging ...... 123
6.4.4 Conclusions ...... 125
6.5 Power Usage Routine ...... 127
6.6 Lander Heating ...... 127
7 UAV Design and Power 128
7.1 Introduction ...... 128
7.2 Operational Conditions on Titan ...... 129
7.2.1 Location ...... 129
7.2.2 Atmosphere and Terrain ...... 129
7.2.3 Meteorological activity ...... 130
7.3 Design Decisions ...... 131
7.3.1 Heavier or Lighter than Air flight? ...... 131
7.3.2 Balloon Type ...... 132
7.3.3 Propulsion ...... 134
7.3.4 Operating Height ...... 135
7.3.5 Materials ...... 136
7.4 Power Decisions ...... 136
7.4.1 Batteries ...... 137
7.4.2 Solar Power ...... 137
7.4.3 RTGs ...... 138
6 CONTENTS
7.4.4 Other options ...... 139
7.4.5 RHUs ...... 139
7.5 Balloon Thermodynamics ...... 140
7.6 Conclusion ...... 140
8 UAV Control 141
8.1 Introduction ...... 141
8.1.1 Nomenclature ...... 141
8.2 Hot Air Balloon Dynamics ...... 142
8.2.1 Atmospheric properties ...... 142
8.2.2 Balloon Equations of Motion ...... 143
8.2.3 Linearisation ...... 145
8.2.4 Numerical Parameters ...... 146
8.3 Modelling ...... 147
8.3.1 Linear Model ...... 147
8.3.2 Non-Linear Model ...... 148
8.3.3 Comparison of Models ...... 148
8.4 Thermodynamics ...... 151
8.4.1 Thermodynamic Model ...... 152
8.5 Titanic ballooning ...... 153
8.5.1 A More Refined Atmospheric Model ...... 154
8.5.2 Valving ...... 155
8.5.3 Building a Controller ...... 157
8.6 Testing the Controller ...... 158
8.6.1 Sinusoidal Demand ...... 158
7 CONTENTS
8.6.2 Entry Scenario ...... 158
8.6.3 Landing ...... 161
8.6.4 Disturbance Rejection ...... 163
8.7 Horizontal Control ...... 163
8.7.1 Wind Data ...... 163
8.7.2 Horizontal Dynamics ...... 165
8.7.3 The Horizontal Controller ...... 165
8.8 Conclusion ...... 169
9 Remote Sensing 171
9.1 Introduction ...... 171
9.2 Information to Sense ...... 171
9.2.1 Temperature ...... 172
9.2.2 Chemical ...... 173
9.2.3 Life ...... 175
9.2.4 Light ...... 176
9.2.5 Sound ...... 177
9.2.6 Acceleration ...... 177
9.2.7 Wind Speed ...... 179
9.3 Communications ...... 181
9.3.1 Radiation Hardening Chips ...... 181
9.3.2 Wired Communication ...... 182
9.3.3 Light Communication ...... 183
9.3.4 Wireless Protocols ...... 183
9.3.5 Localisation ...... 185
8 CONTENTS
9.4 Module Design ...... 186
9.4.1 Antenna Design ...... 186
9.4.2 Landing Protection ...... 188
9.4.3 Flotation ...... 190
9.4.4 Storage and Deployment ...... 191
9.5 Conclusion ...... 193
10 Deep space communications 194
10.1 Introduction ...... 194
10.1.1 Challenges and Requirements ...... 194
10.2 Link Budget and Signal to Noise Ratio ...... 195
10.2.1 Frequency selection ...... 198
10.2.2 Path loss ...... 200
10.2.3 Orbiter antenna ...... 200
10.2.4 Earth antenna ...... 201
10.2.5 Conclusion ...... 201
10.3 Modulation scheme ...... 201
10.3.1 Appropriate schemes ...... 202
10.3.2 Spectral efficiency ...... 204
10.3.3 Pulse shaping ...... 206
10.3.4 Bit error rate and noise rejection ...... 207
10.3.5 Microwave mixing ...... 209
10.4 Simulation of modulator/transmitter and receiver/demodulator ...... 209
10.4.1 Transmitter design ...... 209
10.4.2 Receiver design ...... 210
9 CONTENTS
10.4.3 Demodulation ...... 210
10.5 Hardware and software implementation ...... 214
10.5.1 Implementation of the digital stages ...... 214
10.5.2 Mixers ...... 216
10.5.3 Amplification ...... 216
10.5.4 Conclusion ...... 216
10.6 Internal communications ...... 217
10.6.1 Internal Protocol ...... 217
11 Navigation Systems 219
11.1 Introduction ...... 219
11.2 Coordinate Systems ...... 220
11.3 Orbiter Navigation ...... 221
11.3.1 Potential Methods and Analysis ...... 222
11.3.2 Software ...... 225
11.3.3 Hardware ...... 227
11.4 UAV and Lander Navigation ...... 228
11.4.1 Potential Methods and Analysis ...... 228
11.4.2 Conclusions ...... 233
12 Titanic Communication Systems 234
12.1 Overview of UAV, Lander and Orbiter Communication ...... 234
12.2 Antennas ...... 234
12.3 Amplifiers ...... 236
12.4 Noise ...... 236
10 CONTENTS
12.5 Range of Transmission ...... 237
12.6 Conclusions ...... 240
13 Imaging Systems 241
13.1 Introduction ...... 241
13.2 Systems Overview ...... 241
13.2.1 The Three Imaging Systems ...... 241
13.2.2 System Specialisations ...... 242
13.3 Base Imaging System ...... 242
13.3.1 CCD Sensor ...... 243
13.3.2 Computer Hardware ...... 243
13.3.3 Operating System ...... 244
13.3.4 Image Compression ...... 245
13.3.5 Image Storage ...... 248
13.3.6 Hierarchy of Transmission ...... 248
13.4 Lander Imaging System ...... 249
13.4.1 Camera Mounting ...... 249
13.4.2 Lamp ...... 249
13.4.3 Autofocusing Method ...... 250
13.4.4 Autofocus Noise Rejection ...... 251
13.4.5 Focus Window ...... 252
13.4.6 Climbing Search Algorithm ...... 253
13.5 Orbiter Imaging System ...... 253
13.5.1 Lens Filter ...... 254
13.5.2 Infinity Focus ...... 255
11 CONTENTS
13.5.3 Mapping Algorithm ...... 255
13.5.4 Modeling the Orbit ...... 256
13.5.5 Locating previous photos which have captured portions of the current view . . 257
13.5.6 Determining the useful region of a photo ...... 258
13.5.7 The complete Simulation ...... 261
13.5.8 Optimisation and Implementation ...... 263
13.5.9 Choosing parameters based on simulation ...... 263
13.6 UAV Imaging System ...... 266
13.6.1 Infinity Focus ...... 266
13.6.2 Lens Filter ...... 267
13.7 Conclusion ...... 267
Appendices 268
A Calculating Launch Windows In Matlab 269
B The PD Controller in Matlab 272
C Optimal Control: LQR Controller in Matlab 273
D Estimator Design in Matlab 275
E Orbiter Navigation Simulation in MATLAB 277
Bibliography 280
12 CHAPTER 1 - Hugo Grimmett
Introduction
1.1 Mission to Titan
This document details the technical specifications of a mission to explore the moon Saturn VI, more commonly known as Titan. The mission is to send various autonomous units to the distant moon in order to learn more about it. Having arrived at their destination, these units are to take detailed photographs of the surface, take measurements of radiation levels, temperature, wind speed and direction, and any seasonal change over the mission duration, which is estimated to be greater than
365 days.
1.2 Titan’s Characteristics
Titan is interesting because it is a very large body with a dense nitrogen-based atmosphere (with radius 0.404 times that of Earth’s [154]). This makes it difficult to take accurate readings (such as those mentioned above) from Earth using telescopes, and so a closer encounter is required. A previous mission which took readings on Titan, the Huygens probe from the Cassini orbiter, found methane lakes and methane clouds, but as of yet there has been no mission dedicated to learning more about this seemingly impenetrable moon. Very little is known about the composition of the core, although various hypotheses about layers of ice and ammonia-rich water exist [213]. However, it is certain that there is no metallic core, and hence no magnetic field.
13 Chapter 1. Introduction Hugo Grimmett
1.3 Mission Components
There are three units or parts to the mission apparatus: an orbiter, a lander, and a hot air balloon
(or Unmanned Aerial Vehicle, UAV). Once in a stable orbit, the orbiter will drop the lander near the Northern pole of Titan, which in turn will release the UAV during its descent. The lander will establish a permanent position, and the UAV will circle the moon at an altitude of approximately
10km. The orbiter will gather the data accumulated by its two counterparts and relay them back to
Earth.
There are eleven members working on this project, and they have been split into groups of three or four to work on specific domains. The domains are as follows:
• Deployment. This group is responsible for getting the orbiter to Titan, and making sure the
lander and orbiter are released correctly. Getting the lander to the correct landing velocity
requires carefully-timed parachute releases, which if performed incorrectly, could result in a
crash. This group is also responsible for analysing the wind characteristics on the moon and
modelling how the UAV will move over time.
• Power. The three units will require various power sources for electricity and heat, and all have
different requirements for duration, intensity and mass. The electrical systems used for sensing,
processing and communication will need to be kept above certain minimum temperatures in
order to function correctly.
• Navigation, communication and sensing. This involves the data gathered by the three units, and
how each datum is “titan-tagged”, or linked to a precise time and location on Titan corresponding
to when and where it was taken. These data need to be passed between the unit which gathered
them, combined with the time and location information which may reside on another unit, and
transmitted back to Earth via the orbiter.
This report has been compiled and formatted by James Coates.
14 CHAPTER 2 - Joshua McFarlane
Delivery of the Orbiter
2.1 Introduction
To achieve the successful delivery of the probe into an orbit around Titan the route from earth needs careful planning and simulation. In order to orbit around Titan the probe must gain energy to climb out of the Sun’s gravity well to the orbit radius of Saturn. Then it must shed energy to enter an orbit around Saturn at the same orbit radius as Titan; finally an orbit around Titan will be established by further shedding potential energy and becoming trapped in Titan’s gravity well. In order to devise a plan to achieve this the orbital mechanics of the solar system were carefully considered.
With the route planned, it will be implemented on the probe while in space using an on-board con- troller. To test the operation of this controller and the feasibility of the route a simulation was written in Matlab and the results analysed to assess possible improvements to the controller and route.
2.2 Orbital Mechanics
2.2.1 Kepler’s Laws
The trajectory of a body in free orbit around a much more massive body is determined by Keper’s laws of planetary motion: [95]
l r = (2.1) 1 + e cos ν
15 Chapter 2. Delivery of the Orbiter Joshua McFarlane
( ) d 1 r2 ν˙ = 0 (2.2) dt 2
4π2a3 P 2 = (2.3) GM
Where:
• r = the orbital radius
• l = the semi-latus rectum1 (see figure 2.1)
• e = the eccentricity of the orbit
• ν = the true anomaly - the angle through which the body has travelled since passing the
periapsis2 (see figure 2.1)
• P = the orbital period
• a = the semi-major axis3 (see figure 2.1)
• G = the universal gravitational constant (∼ 6.67 × 10−11 m3 kg−1 s−1)
• M = the mass of the central body
These equations are simplified from their original form with the assumption that one of the bodies has a much larger mass than the other and so the centre of mass of the system can be assumed to be coincident with that of the largest body.
The first equation states that orbits follow eliptical paths with the central, massive body at one of the focii. The second equation states that a line drawn between the body in orbit and the body which it is orbiting will sweep out equal areas in equal times (see figure 2.2.) The third equation states that the orbital period squared is proportional to the semi-major axis of the orbit cubed.
1 π The semi-latus rectum is the orbital radius of the body in consideration when ν = 2 2The periapsis is the point at which the orbiting body is closest to the central body 3The semi-major axis is the arithmetic mean of the apoapsis radius and the periapsis radius
16 Chapter 2. Delivery of the Orbiter Joshua McFarlane
Orbiting body Semi-latus rectum
Second focus Central body