Design and Performance Testing of a 1-cm Miniature Radio- Frequency Ion Thruster

Similar documents
Micro RF Ion Engine for Small Satellite Applications

Performance and Thermal Characteristics of High-Power Hydrogen Arcjet Thrusters with Radiation-Cooled Anodes for In-Space Propulsion

ARIANEGROUP ORBITAL PROPULSION ROBERT-KOCH-STRASSE TAUFKIRCHEN GERMANY

DLR s Electric Propulsion Test Facility the First Three Years of Thruster Operation

An Overview of Electric Propulsion Activities in China

Variable Specific Impulse High Power Ion Thruster

L-3 Communications ETI Electric Propulsion Overview

For permission to copy or to republish, contact the copyright owner named on the first page. For AIAA-held copyright, write to AIAA Permissions

QinetiQ Electric Propulsion

SMALLSAT PROPULSION. Pete Smith, Roland McLellan Marotta UK Ltd, Cheltenham, and Dave Gibbon SSTL, Guildford, UK.

Resistojet Thrusters for Auxiliary Propulsion of Full Electric Platforms

ESA Propulsion Laboratory at ESTEC

The Electric Propulsion Development in LIP

RDT&E BUDGET ITEM JUSTIFICATION SHEET (R-2 Exhibit) June 2001

Performance Characteristics of Low-Power Arcjet Thruster Systems with Gas Generators for Water

ESA Propulsion Laboratory

1,000-hour Running of a 20-mN Ion Thruster with Pyrolytic Graphite Grids IEPC

Evaluation of 25-cm XIPS Thruster Life for Deep Space Mission Applications

Experimental Testing of a Rotating Detonation Engine Coupled to Nozzles at Conditions Approaching Flight

High-Voltage Terminal Test of Test Stand. for 1-MV Electrostatic Accelerator

25 B43 B43.1 THE MEASUREMENT OF e/m BY THE BAINBRIDGE METHOD

Performance and Flatness of a Multiple-Cathode, Rectangular Ion Thruster Discharge Chamber

This chapter gives details of the design, development, and characterization of the

Experimental Investigation of Hot Surface Ignition of Hydrocarbon-Air Mixtures

Oil Free Vacuum Pumps for Mass Spectrometry

Development of Low-thrust Thruster with World's Highest Performance Contributing to Life Extension of Artificial Satellites

AUTOMOTIVE EMC TEST HARNESSES: STANDARD LENGTHS AND THEIR EFFECT ON RADIATED EMISSIONS

Annular-Geometry Ion Engine: Concept, Development Status, and Preliminary Performance

CFD Investigation of Influence of Tube Bundle Cross-Section over Pressure Drop and Heat Transfer Rate

Electric Propulsion Electronics Activities in Astrium Germany

EPIC Gap analysis and results

Development of a 50,000-s, Lithium-fueled, Gridded Ion Thruster

EXTENDED GAS GENERATOR CYCLE

PNEUMATIC HIGH SPEED SPINDLE WITH AIR BEARINGS

Life and Operating Range Extension of the BPT 4000 Qualification Model Hall Thruster

NEXT Single String Integration Test Results

APPLICATION OF VARIABLE FREQUENCY TRANSFORMER (VFT) FOR INTEGRATION OF WIND ENERGY SYSTEM

AMBR* Engine for Science Missions

Aerospace Engineering Aerospace Vehicle System. Introduction of Propulsion Engineering

Results of a 2000-Hour Wear Test of the NEXIS Ion Engine

Multipulse Detonation Initiation by Spark Plugs and Flame Jets

Theoretical and Experimental Investigation of Compression Loads in Twin Screw Compressor

Propulsion Solutions for CubeSats and Applications

Cassini-Huygens Power Conversion Technology

Formation Flying Experiments on the Orion-Emerald Mission. Introduction

Supersonic Combustion Experimental Investigation at T2 Hypersonic Shock Tunnel

Results of a 35-cm Xenon Ion Thruster Endurance Test *

Wireless Energy Transfer Through Magnetic Reluctance Coupling

Development of a Nitrous Oxide Monopropellant Thruster

NASA s Electric Propulsion Program

Design Rules and Issues with Respect to Rocket Based Combined Cycles

719. Diagnostic research of rotor systems with variable inertia moment

MIRI Cooler System Design Update

Palamede, more than a microsatellite. Workshop on University Micro Satellites in Italy Rome, July 27, 2005

Demonstration of the XR-12 Hall Current Thruster

Flight Demonstration and Application of Electric Propulsion at CAST

Available online at ScienceDirect. Procedia CIRP 33 (2015 )

ARIANEGROUP ORBITAL PROPULSION ROBERT-KOCH-STRASSE TAUFKIRCHEN GERMANY

Development History and Current Status of DC-Type Ion Engines at JAXA

The Effects of Magnetic Circuit Geometry on Torque Generation of 8/14 Switched Reluctance Machine

Gujarat, India,

Conversion of Automotive Turbocharger to Gas Turbine

CRYOGENIC MOTORS FOR HERSCHEL/PACS AND JAMES WEBB/MIRI AND NIRSPEC

Dr. Jim Henry, P.E. Professor of Engineering University of Tennessee at Chattanooga 615 McCallie Avenue Chattanooga, TN Dr.

Proposal to establish a laboratory for combustion studies

CFD Analysis on a Different Advanced Rocket Nozzles

Additively Manufactured Propulsion System

Static Structural and Thermal Analysis of Aluminum Alloy Piston For Design Optimization Using FEA Kashyap Vyas 1 Milan Pandya 2

DESIGN AND EXPERIMENTATION OF TEST RIG TO CHARACTERIZE HYDROSTATIC DRIVEFOR LINEAR ACTUATOR

Cooling Enhancement of Electric Motors

Study on Flow Fields in Variable Area Nozzles for Radial Turbines

LPT6510 Pulse-tube Cooler for K applications

DEVELOPMENT STATUS OF NEXT: NASA S EVOLUTIONARY XENON THRUSTER

Study on Flow Characteristic of Gear Pumps by Gear Tooth Shapes

HYSYS System Components for Hybridized Fuel Cell Vehicles

Optimization of Three-stage Electromagnetic Coil Launcher

PERFORMANCE ESTIMATION AND ANALYSIS OF PULSE DETONATION ENGINE WITH DIFFERENT BLOCKAGE RATIOS FOR HYDROGEN-AIR MIXTURE

Integration Test of the High Voltage Hall Accelerator System Components

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Aeronautics and Astronautics

Effect of Compressor Inlet Temperature on Cycle Performance for a Supercritical Carbon Dioxide Brayton Cycle

100 kw Nested Hall Thruster System Development

Six Decades of Thrust - The Ariane Group Radiofrequency Ion Thrusters and Systems Family

D.A. Barnhart*, J.M. McCombet, D.L. Tilley$ Air Force Phillips Laboratory Edwards A.F.B., CA

PPU MK3 FOR 5 KW HALL EFFECT THRUSTERS 11TH EUROPEAN SPACE POWER CONFERENCE

All-Electric CubeSat Propulsion Technologies

ME3264: LAB 9 Gas Turbine Power System

An Advanced Compressor for Turbo-Brayton Cryocoolers

Aircraft Overview. Snowbird Human-Powered Ornithopter. Todd Reichert Cameron Robertson University of Toronto Institute for Aerospace Studies Canada

GEO Dual Mode PPU & LEO HEMPT PPU

INFLUENCE OF THE NUMBER OF NOZZLE HOLES ON THE UNBURNED FUEL IN DIESEL ENGINE

IMPACT OF SKIN EFFECT FOR THE DESIGN OF A SQUIRREL CAGE INDUCTION MOTOR ON ITS STARTING PERFORMANCES

Comparative Study and Analysis of Air Ejector Flow in Convergent and Convergent Divergent Nozzle of Aircraft

CONJUGATE HEAT TRANSFER ANALYSIS OF HELICAL COIL HEAT EXCHANGE USING CFD

IAC-15-C4.3.1 JET INDUCER FOR A TURBO PUMP OF A LIQUID ROCKET ENGINE

Cryocooler with Cold Compressor for Deep Space Applications

Flywheel energy storage retrofit system

Cooldown Measurements in a Standing Wave Thermoacoustic Refrigerator

Rotating Detonation Wave Stability. Piotr Wolański Warsaw University of Technology

Computational Fluid Dynamics in Torque Converters: Validation and Application

Development, Qualification and Delivery Status of the HEMPT based Ion Propulsion System for SmallGEO

Transcription:

Design and Performance Testing of a 1-cm Miniature Radio- Frequency Ion Thruster IEPC-009-167 Presented at the 31st International Electric Propulsion Conference University of Michigan Ann Arbor, Michigan USA Thomas A. Trudel, 1 Sven G. Bilén, and Michael M. Micci 3 The Pennsylvania State University, University Park, Pennsylvania 1680, USA This paper presents the design and performance testing of the 1-cm Miniature Radio- Frequency Ion Thruster (MRIT) currently under development at The Pennsylvania State University. The primary goal of the MRIT program is to produce progressively smaller, micronewton range, RF ion propulsion thrusters for use in high precision spacecraft attitude control and on miniaturized spacecraft. The MRIT is an inductively coupled RF ion thruster currently using argon gas as a propellant. The latest MRIT model is the smallest iteration to date with a conical plasma chamber that is 1 cm in both diameter and length. The total thruster diameter and length are each just over cm. The thruster has been tested in a cylindrical vacuum chamber at pressures between 10 4 and 10 6 Torr. The MRIT produces a maximum thrust and specific impulse of 59.0 µn and 5480 s, respectively. The overall performance of this MRIT iteration was characterized by a series of tests conducted at PSU from 008 009. Selected data and results are presented here. I. Introduction iniature ion propulsion devices show great promise as a means of providing finite attitude control for satellites Mand other spacecraft due to their high specific impulse and operational efficiencies as well as their exceedingly low fuel consumption. 1 A small RF ion thruster, like the Miniature Radio-Frequency Ion Thruster (MRIT) under development at The Pennsylvania State University, could produce low levels of thrust, on the order of 1 to 50 μn, at a precise thrust resolution providing consistent and accurate positioning capability over the lifetime of a spacecraft. These micronewton-scale thrusters could be used for applications such as routine satellite station keeping, precise attitude control in formation flying spacecraft, or potentially as a primary propulsion device on a miniature spacecraft. Specific uses such as drag compensation, orbit control and reconfiguration, disturbance cancellation, and controlled formation flying are commonly noted. The objective of the MRIT program is to produce progressively smaller, micronewton range, radio frequency (RF) ion propulsion thrusters. The previous iteration of the MRIT was a cylindrical thruster whose plasma chamber measured 1.5 cm in both diameter and length. It used a two-grid extraction system operating at a total voltage of 100 V to produce approximately 75 μn of thrust with an I sp of 400 s. 3 Concerns with the previous MRIT design included the accelerator grid directly intercepting ions accelerated by the screen grid and beam steering that was occurring due to grid misalignment and slight inaccuracies in grid construction. Solving these optics problems was one of the primary goals of this MRIT iteration. In addition, the change from a cylindrical to conical plasma chamber was made in an effort to increase the chamber volume to surface area ratio and thereby decrease discharge losses to the chamber walls and increase thruster efficiency. 4 The MRIT iteration presented in this paper continued the miniaturization of the previous MRIT and produced a smaller, more efficient, and more reliable thruster. The MRIT design presented in this paper is shown in Fig. 1. It has a conical plasma chamber 1.0 cm in both diameter and length. Total thruster length and diameter are each just over.0 cm. It is an inductively coupled RF ion 1 Undergraduate Student, Aerospace Engineering, tat5011@gmail.com. Associate Professor, Engineering Design, Electrical Engineering, and Aerospace Engineering, sbilen@psu.edu. 3 Professor, Aerospace Engineering, micci@psu.edu. 1

thruster with a two-grid, 7-apertures-per-grid, extraction system. A schematic overview of the thruster is shown in Fig.. 5 During the previous MRIT iteration it was shown that the thruster operates at a higher efficiency using an RF frequency of 1.0 1.5 MHz rather than the previously used 13.56 MHz. Figure 1. MRIT size comparison. Figure. MRIT system diagram. In order to reduce costs, the laboratory model of the thruster was constructed from materials somewhat different from those typically used in ion thrusters. Most notably, the extraction grids were constructed from stainless steel and the assembly that houses them was constructed from several pieces of Teflon. The flight model of the MRIT will utilize molybdenum extraction grids and alumina ceramic extraction assembly components. The MRIT plasma chamber and Faraday cage surrounding the thruster were constructed from an alumina ceramic and 6061 aluminum alloy. It is likely these materials will remain the same in the flight model thruster. Testing with xenon propellant was also postponed to reduce costs. Argon is significantly cheaper and easier to obtain. The results presented in this paper were obtained using argon gas as a propellant, but testing using xenon gas will be carried out in the near future. Xenon gas is the intended thruster propellant. II. Testing Facility and Experimental Setup A. Test Facility The facility used to test the MRIT system is located at Penn State. The vacuum chamber used is shown in Fig. 3. The chamber is approximately 0.6 m in diameter and 1.0 m in depth. The vacuum system is a two-stage system that uses a BOC Edwards IPUP scroll pump and a CTI-Cryogenics Cryo-Torr 10 series cryopump to reach an ultimate pressure as low as 1.0 10 6 Torr. B. Experimental Equipment The pressure inside the vacuum chamber is measured using an MKS Series 999 Multi Sensor Pressure Transducer and an Inficon CC3 cold cathode vacuum gauge. The MKS Series 999 Sensor measures pressures from atmosphere to 10 10 Torr with an accuracy of 0.01 10 n Torr, where n is the order of magnitude at which work is being conducted. The Inficon CC3 measures pressures from 10 3 to 10 8 Torr with similar accuracy. Both sensors were used and their readings compared to ensure accurate pressure measurements. A Horiba Stec mass flow controller (MFC) was used in combination with an MKS 147B control box to regulate the flow of propellant into the thruster s ionization chamber. Cold flow tests were conducted frequently to ensure proper propellant flow rate within the plasma chamber. The MFC was set to various flow rates and the change in chamber pressure was recorded. Two Bertan 05B series high voltage sources were used to provide the positive and negative high voltage used on the screen and acceleration grids, respectively. The high voltage sources are internally grounded and, in addition, the cables carrying the high voltages are grounded to the metallic vacuum chamber via their Figure 3. Thermal vacuum chamber.

respective feed-throughs. To generate the RF field, a Hewlett Packard 3310A Arbitrary Waveform Generator was used to produce the necessary signal at a frequency of 1.5 MHz. It was then amplified by an RF Power Labs Model ML50 RF amplifier to provide the necessary range of RF power levels. The signal was tuned using an MFJ Deluxe Versa Tuner V and a 4:1 matching transformer to provide the lowest possible SWR. C. Experimental Setup The MRIT exhaust plume was analyzed using a 3.5-cm-long Faraday cup with a 0.1-cm aperture. The cup was placed 5 cm downstream of the MRIT exit grids and aligned axially with the center of the thruster. The Faraday cup was placed in this position to make all beam current density measurements except for those associated with beam divergence testing. To test for beam divergence the Faraday cup was placed 5 cm downstream of the thruster exit in a radial pattern from 0 degrees to 4 degrees on either side of its initial position. Beam current density measurements were taken every degrees throughout this range to produce a two dimensional beam current density profile for the MRIT exhaust plume. All other thruster parameters (e.g., total ion current, accelerator grid current, etc.) were recorded from the various laboratory power supplies, power meters, and mass flow control systems listed above. Thruster operational parameters were then calculated from this data and the beam current density profile attained from the Faraday cup testing. To test for proper beam neutralization a Langmuir probe was placed 5 cm and 10 cm downstream of the MRIT. The probe was excited by a potential of 50 to +50 V and the collected probe current was recorded. At each throttle point the required current through the neutralizer to produce a neutral beam (no current reading at LP location) was recorded. III. Experimental Results and Discussion A. Operational Parameters Before analysis of the MRIT exhaust plume and other performance characteristics could begin, the operational limits of the thruster and a reliable operating procedure had to be established. The primary characteristics that needed to be known were: functional propellant flow rates, minimum operational RF power level, and functional exit grid potential. A photograph of the operational thruster can be seen in Fig. 4. To determine the functional flow rate, first the necessary flow rate to achieve a chamber pressure between 10 and 10 3 Torr was calculated using Poiseuille s equation. 6 A flow rate of 0.03 sccm to 0.04 sccm was calculated for the current MRIT geometry. The MRIT successfully sustained operation with flow rates ranging to as low as 0.0 0.1 sccm. Maximum propellant efficiency occurred at the lowest functional flow rates. A 4:1 matching transformer was used between the RF amplifier and the thruster coil, with SWR values as low as 1.5 consistently obtained. It was found that, at an SWR of.0, the thruster could function with as little as 10 W of delivered power. Figure 4. MRIT firing. The two sets of tests above were both conducted with a total exit grid potential of 100 V. The screen grid was charged to +1000 V and the acceleration grid was set to 00 V. These were the standard potentials used on the previous iterations of the MRIT and produced consistent operation here. Arcing occurred between the grids at total voltages above 1700 V and the thruster would lose stability at total voltages below 500 V. In the results presented herein, the steady state operating conditions for the MRIT are: +1000 V and 00 V screen and acceleration grid potential, respectively, 0.035 sccm propellant flow rate, and a delivered RF power level of 15 W. B. Exhaust Plume Characteristics The exhaust plume of the MRIT was analyzed using the Faraday cup setup outlined in Section II.C. The average beam current density at steady state conditions was on the order of.0 ma/cm. At these conditions the MRIT produces a calculated thrust of.5 µn at a specific impulse of 096 s. A two dimensional beam current density profile for the MRIT is shown in Fig. 5. This profile was used to determine the beam divergence of the MRIT exhaust plume. The beam divergence is approximately 17 degrees. 3

Thrust values reported in this paper are calculated from the beam current measurements and the beam divergence value determined via experimentation..5 Beam Current Density (ma/cm^) 1.5 1 0.5 0-30 -0-10 0 10 0 30 Angular Position (degrees) Figure 5. Two-dimensional beam current density profile. C. Neutralizer A simple tungsten filament 1 mm wide by 1.5 cm in length was used to neutralize the ion beam expelled by the thruster. As is typical with electric propulsion systems, the neutralizer self-biases to automatically expel the required number of electrons to effectively neutralize the thruster exhaust beam. 7 A current of 3.75 A to 6.03 A at a potential of 3.11 V to 5.7 V was fed through the filament in order to neutralize the MRIT ion beam at all throttle points (on the order of 10 14 10 17 electrons/m 3 ). This is a simple, yet effective, method for neutralizing the exhausted ion beam. D. Electron Backstreaming Tests Testing was done to determine the electron backstreaming limit on the MRIT optics. The backstreaming limit was defined as the acceleration grid voltage at which a 1% increase in screen grid current was seen (due to backstreaming electrons). 8 Figure 6 below shows the results of these tests and the identified backstreaming limit acceleration grid potential of 60 V. Screen Grid Current (ma) Figure 6. Electron backstreaming limit. E. Throttling MRIT throttling is accomplished by varying the total voltage applied to the exit grids. Throttling by varying the delivered RF power was also investigated, but it was found that better thrust control and resolution was achieved when using the optics. Throttling using both systems will continue to be investigated as the MRIT electronics evolve. Results from the optics throttling tests are shown below in Fig. 7. The minimum and maximum thrust levels were 1.45 µn and 59.01 µn, respectively. Thrust intervals ranged from 4 µn to 10 µn..04.0 1.98 1.96 1.94 1.9 1.9 500 400 300 00 100 0 Acceleration Grid Potential (V) 4

70.00 60.00 Thrust (μn) 50.00 40.00 30.00 0.00 10.00 0.00 0 00 400 600 800 1000 Time (s) Figure 7. MRIT thrust vs. time (demonstrating throttling). IV. Performance Analysis The current MRIT iteration functioned more reliably than any previous iteration. The thruster produced a maximum calculated operational thrust of 59.0 µn with a specific impulse of 5480 s and a mass efficiency on the order of 60% to 80% depending on propellant flow rate. Figure 8 gives the propellant efficiency as a function of thrust at multiple propellant flow rates. 100.0 Mass Efficiency (%) 80.0 60.0 40.0 0.0 0.0 0.00 0.00 40.00 60.00 80.00 Thrust (μn) 0.03 0.04 0.05 0.06 0.07 0.08 0.09 Figure 8. MRIT mass efficiency vs. thrust at multiple propellant flow rates. The MRIT operated continuously for over 50 hours without any incident or fault in operation. During these 50 hours, the thruster was throttled up and down over the entire operational thrust range, shut down and re-ignited numerous times, and subjected to various changes in propellant flow rate, exit grid potential, and delivered RF power. In every case the MRIT operated without fail provided these parameters remained within the bounds of the operational parameters outlined in Section III.A. This is the longest and most successful series of testing any MRIT iteration has undergone. The power-to-thrust ratio was 300 W/mN with an electrical efficiency of 15.%. The power to thrust ratio is high and the electrical efficiency suffers accordingly. This is due in part to the exceedingly small size of the MRIT; however, as the MRIT flight electronics evolve and the transition from laboratory electronics is made, a more power efficient thruster is a primary goal. Total efficiency was 1.84%. The total efficiency suffered as a direct result of the low electrical efficiency. Once again, it is hoped the total efficiency will increase with advancements in the MRIT support systems. 5

V. Conclusions and Future Work Initial testing of the current iteration of the MRIT showed that the MRIT thruster can operate at RF input power levels as low as 13 W and mass flow rates from as low as 0.0 to 0.1 sccm. Steady-state operational testing was conducted with an RF input power level of 15 W, a mass flow rate of 0.035 sccm, and a total exit grid potential of 100 V. Faraday cup testing conducted under these conditions produced results consistent with theoretical predictions. The MRIT produced 1.45 µn to 59.0 µn of thrust at throttling intervals ranging from 4 µn to 10 µn. At maximum thrust the thruster achieved an I sp of 5480 s. The maximum mass efficiency of the MRIT was just under 80%. Far lower electrical and total efficiencies were produced, and while reduction in these efficiencies is expected with reductions in scale, there is room for improvement. It will be a primary goal of subsequent MRIT iterations to produce more efficient flight electronics. Future work should focus on improving thruster efficiency, conversion to flight materials, and production of MRIT specific flight electronics. Flight model electronics will allow further RF and optic throttling testing to be done and will allow a better prediction of the practical performance of the thruster as these properties are primarily related to the design of the support electronics and not the thruster itself. Acknowledgments We would like to thank The Pennsylvania State University for their continued support of the MRIT program. In addition, we would like to thank students of the Penn State Electric Propulsion group for their input and ideas throughout this project. We would like to thank Prof. Mark Wharton for his assistance with the construction and troubleshooting of the MRIT RF system. Last, T.A.T. would like to thank Dr. Sven Bilén and Dr. Michael Micci for their help and support throughout this project. References 1 Wirz, R., Gale, M., Juergen, M., and Marrese, C., Miniature Ion Thrusters for Precision Formation Flying, AIAA Joint Propulsion Conference, AIAA 004-4115, AIAA, Fort Lauderdale, Florida, 004. Feili, D., Di Cara, D.M., Leiter, H.J., Del Amo, J.G., Loeb, H.W., Weis, St., Kirmse, D., Frigot, P.E., Orlandi, M., Müller, H., Meyer, B.K., The μnrit-4 Ion Engine: a first step towards a European mini-ion Engine System development, 30th International Electric Propulsion Conference, IEPC-007-18, Florence, Italy, September, 007. 3 Trudel, T.A., Design and Testing of a Miniature RF Ion Thruster, B.S. Thesis, The Pennsylvania State University, May 009. 4 Leiter, H.J., Loeb, H.W., and Schartner K.H., The RIT15 Ion Engines A Survey of the Present State of Radio Frequency Ion Thruster Technology and its Future Potentiality, 3rd International Conference on Spacecraft Propulsion, ESA SP-465, Cannes, France, October, 000. 5 Mistoco, V.F., Bilén, S.G., and Micci, M.M., Development and Chamber Testing of a Miniature Radio- Frequency Ion Thruster for Microspacecraft, 40th AIAA Joint Propulsion Conference, Ft. Lauderdale, FL, July, 004. 6 Mistoco, V.F., and Bilén, S.G., Numerical Modeling of a Miniature Radio Frequency Ion Thruster, 44th AIAA Joint Propulsion Conference, Hartford, CT, 008. 7 Wirz, R., Sullivan, R., Przybylowski, J., and Silva, M., Hollow Cathode and Low-Thrust Extraction Grid Analysis for a Miniature Ion Thruster, International Journal of Plasma Science and Engineering, Vol. 008, Article 69385, 008. 8 Goebel, D.M., and Katz, I., Fundamentals of Electric Propulsion: Ion and Hall Thrusters, 1 st edition, JPL California Institute of Technology, March 008, p. 197. 6