Development and qualification of Hall thruster KM-60 and the flow control unit

Similar documents
OPERATION PECULIARITIES OF HALL THRUSTER WITH POWER kw AT HIGH DISCHARGE VOLTAGES

PlaS-40 Development Status: New Results

Development of stationary plasma thruster SPT-230 with discharge power of kw

High-impulse SPT-100D thruster with discharge power of kw

Experimental study of a high specific impulse plasma thruster PlaS-120

Development and Research of the Plasma Thruster with a hollow magnet Anode PlaS-40

Parametric family of the PlaS-type thrusters: development status and future activities

Multi-Mode Thruster with Anode Layer Development Status

Progress in Testing of QM and FM HEMP Thruster Modules

Thrust Balance Characterization of a 200W Quad Confinement Thruster for High Thrust Regimes

Investigation of SPT Performance and Particularities of it s Operation with Kr and Kr/Xe Mixtures *+

Operating Envelopes of Thrusters with Anode Layer

Research and Development of High-Power, High-Specific-Impulse Magnetic-Layer-Type Hall Thrusters for Manned Mars Exploration

Development of Low-Power Cylindrical type Hall Thrusters for Nano Satellite

INVESTIGATION OF THE POSSIBILITY TO REDUCE SPT PLUME DIVERGENCE BY OPTIMIZATION OF THE MAGNETIC FIELD TOPOLOGY IN THE ACCELERATING CHANNEL

Micro-Cathode Arc Thruster Development and Characterization

FLIGHT TESTING OF VOLUME-IONIZATION ION THRUSTERS ON METEOR SATELLITE

Development of Microwave Engine

Commissioning of the Aerospazio s vacuum facilities with Safran s Hall Effect Thruster

Characterization of an adjustable magnetic field, low-power Hall Effect Thruster

Research and Development of Very Low Power Cylindrical Hall Thrusters for Nano-Satellites

Multiple Thruster Propulsion Systems Integration Study. Rusakol, A.V..Kocherpin A.V..Semenkm A.V.. Tverdokhlebov S.O. Garkusha V.I.

Flight Demonstration of Electrostatic Thruster Under Micro-Gravity

VHITAL-160 thermal model building and 3D thermal modeling carrying-out

Development Status of 200mN class Xenon Hall Thruster of MELCO

SPT Operation in Machine-Gun Mode

Experimental Investigation of Magnetic Field Topology Influence on Structure of Accelerating Layer and Performance of Hall Thruster

Miniature Vacuum Arc Thruster with Controlled Cathode Feeding

Abstract. Introduction

Abstract. Objectives. Theory

OPERATIONAL CHARACTERISTICS OF CYLINDRICAL HALL THRUSTERS

New 2d Far Field Beam Scanning Device at DLR s Electric Propulsion Test Facility

Figure 1, Schematic Illustrating the Physics of Operation of a Single-Stage Hall 4

PPS 1350-G Performance assessment with permanent magnets

ALCATEL SPACE PLASMA PROPULSION SUBSYSTEM QUALIFICATION STATUS

Telemicroscopy Erosion Measurements of 5 kw-class Hall Effect Thruster Channel Walls

RESEARCH ON TWO-STAGE ENGINE SPT-MAG

In-situ temperature, grid curvature, erosion, beam and plasma characterization of a gridded ion thruster RIT-22

Pole-piece Interactions with the Plasma in a Magnetic-layertype Hall Thruster

Cold Gas Thruster Qualification for FORMOSAT 5

Electric Propulsion Activity in Russia *

PPS 1350-G Qualification status h

Study of Low Power TAL Characteristics

Advanced laboratory for testing plasma thrusters and Hall thruster measurement campaign

DESIGN AND PRELIMINARY CHARACTERIZATION OF A 5 KW HALL THRUSTER PROTOTYPE

BPT-4000 Hall Thruster Extended Power Throttling Range Characterization for NASA Science Missions

Operation Characteristics of Diverging Magnetic Field Electrostatic Thruster

Plasma Behaviours and Magnetic Field Distributions of a Short-Pulse Laser-Assisted Pulsed Plasma Thruster

The µnrit-4 Ion Engine: a first step towards a European mini-ion Engine System development.

Performance Measurements of a High Powered Quad Confinement Thruster.

Status of the THALES High Efficiency Multi Stage Plasma Thruster Development for HEMP-T 3050 and HEMP-T 30250

Applied-Field MPD Thruster with Magnetic-Contoured Anodes

Development of a Two-axis Dual Pendulum Thrust Stand for Thrust Vector Measurement of Hall Thrusters

Sitael Low Power Hall Effect Thrusters for Small Satellites

Evaluation of Quasi-Steady Operation of Applied Field 2D- MPD Thruster using Electric Double-Layer Capacitors

Geometry optimization and effect of gas propellant in an electron cyclotron resonance plasma thruster

Very High I sp Thruster with Anode Layer (VHITAL): An Overview

Effects of Background Pressure on the NASA 173M Hall Current Thruster Performance

Some results of the small power SPT models creation

Assessment of the Azimuthal Homogeneity of the Neutral Gas in a Hall Effect Thruster using Electron Beam Fluorescence

Electric Rocket Engine System R&D

THE PPS 1350 PROGRAM

Ten-Ampere-Level, Direct Current Operation of Applied-Field Magnetoplasmadynamics (MPD) Thruster using LaB 6 Hollow Cathode

Characterization of a Cylindrical Hall Thruster with Permanent Magnets

3 A NEW FRENCH FACILITY FOR ION PROPULSION RESEARCH

HEMP Thruster Assembly Performance with increased Gas Tubing Lengths of Flow Control Unit

A review of plasma thruster work at the Australian National University

Alternative Neutralization Technique for a 40 Watt Quad Confinement Thruster

- 581 IEPC the ion beam diagnostics in detail such as measurements of xenon with double charges, Introduction

Development and Test of XR-150, a New High-Thrust 100 W Resistojet

Investigation of a 5 kw class Hall-effect thruster operating with different xenon-krypton mixtures

Development of Micro-Vacuum Arc Thruster with Extended Lifetime

Endurance Testing of HEMPT-based Ion Propulsion Modules for SmallGEO

Sergey O. Tverdokhlebov

Development and Testing of a New Type of MPD Thruster #

Effect of a Plume Reduction in Segmented Electrode Hall Thruster. Y. Raitses, L.A. Dorf, A A. Litvak and N.J. Fisch

Plasma Properties Inside a Small Hall Effect Thruster

Langmuir Probe Measurements of a Magnetoplasmadynamic Thruster

Kinetic simulation of the stationary HEMP thruster including the near field plume region

Subscale Lifecycle Test of Thermal Arcjet Thruster TALOS for the Lunar Mission BW1

Hall Thruster Modifications for Reduced Power Operation

The division of energy sources and the working substance in electric propulsioncan determines the range of applicability of electro jet propulsion sys

Computer modeling and comparison with diagnostic results of a downscaled High Efficiency Multi Stage Plasma Thruster

Initial performance characterisation of a plasma thruster employing magnetic null regions

INTEGRAL AND SPECTRAL CHARACTERISTICS OF ATON STATIONARY PLASMA THRUSTER OPERATING ON KRYPTON AND XENON

The Quad Confinement Thruster - Preliminary Performance Characterization and Thrust Vector Control

ETS-Ⅷ Ion Engine and its Operation on Orbit

Characterization of a Colloid Thruster Performing in the micro-newton Thrust Range.

Development of IPPLM s Krypton HET

Realization of Low Frequency Oscillation Free Operation in a Hall Thruster

S thruster, which is widely used at present on satellites with 1 kw power on bord 3. The main advantages of the

Propulsion means for CubeSats

Downscaling a HEMPT to micro-newton Thrust levels: current status and latest results

(b) Analyzed magnetic lines Figure 1. Steady state water-cooled MPD thruster.

Experimental Study of a 1-MW-Class Quasi-Steady-State Self-Field Magnetoplasmadynamic Thruster

Characteristics of Side by Side Operation of Hall Thruster

Experimental Investigations of a Krypton Stationary Plasma Thruster

Influence of Electrode Configuration of a Liquid Propellant PPT on its Performance

The Effect of the Discharge Chamber Structure on the Performance of a 5 cm-diameter ECR Ion Thruster

Ultra-Low Power Stationary Plasma Thruster

Transcription:

Development and qualification of Hall thruster KM-60 and the flow control unit IEPC-2013-055 Presented at the 33rd International Electric Propulsion Conference, The George Washington University Washington, D.C. USA A.N.Kostin 1), A.S.Lovtsov 2), A.I.Vasin 3), V.V.Vorontsov 4) SSC FSUE Keldysh Research Centre, Onezhskaya, 8, Moscow, Russia, 125438 Abstract: Qqualification of correction unit consists of a thruster KM-60 and a flow control unit has been completed. Overview of the correction unit tests including lifetime tests is presented in this paper. The maximum correction unit operating time during the test is more than 4120 hours. At that the shown total trust impulse more than 6 10 5 N s at 8357 cycles. Average specific impulse shown during maximum lifetime test is more than 1860 s. Nomenclature = coefficient characterizing the reduction of thrust due to the angular divergence of the ion beam E = coefficient characterizing the reduction of thrust due to the non monoenergetic ion beam = gas efficiency = the proportion of the ion current to the discharge current CU = correction unit FCU = flow control unit PPCU=power processing and control unit I. Introduction Nowadays increasing of spacecrafts lifetime and on-board power systems capacity are among main directions in progress of spacecraft systems engineering. As a result, new requirements for propulsion system specific impulse and lifetime arise. Hall thruster KM-60 and flow control unit (FCU) developed at the Keldysh Research Centre are intended for use on board the satellite platform "Express-1000" as correction unit (CU). CU has increased specific impulse in relation to existing models, which exceeds 2000 seconds in the beginning of the operation. 1 Leading engineer, department of electrophysics 2 Head of department of electrophysics, lovtsov@kerc.msk.ru 3 Leading staff scientist, department of electrophysics 4 Leading engineer, department of electrophysics 1

II. Development and Qualifications of the Correction Unit Feature of the developed correction unit is long-term work with a high specific impulse at a relatively low power. Development, manufacturing and qualification of the thruster and FCU have been conducted at Keldysh Research Centre. А. The Main Parameters of the Thruster and the FCU Main parameters of the thruster and FCU are presented in Table 1. Table 1 Main parameters of the Thruster Main parameters of the FCU Thrust 42 mn Propellant Xenon Discharge voltage 500 V Working inlet pressure 1.75 0.15 kgf/cm 2 Discharge current 1.8 А Valves power supply voltage Specific Impulse (acceptance >2000 s At startup (0.1-1 s) 6.3 0.3 V tests) On hold 2.2 0.2 V Average specific impulse (at >1860 s Mass 0.325 kg the lifetime) Total thrust impulse >380 kn s Current in the coils of magnetic 1.5 2.5 А system Mass 3.1 kg In addition, there are a number of requirements to the CU, determined by its functional purpose, such as: large number of cycles (up to 8250 switching cycles); high level of vibration loads (random vibration with a standard deviation up to 25g); high level of single shocks reaching 2500 g. Appearance of correction unit is shown in Fig. 1. а) б) Figure 1. Appearance of correction unit based on thruster KM-60 (a - KM-60 thruster, b - flow control unit) KM-45 was used as prototype of KM-60 thruster having similar type of magnetic system, discharge chamber and cathode units. This thruster had passed lifetime test for more than 1000 hours before starting of KM-60 development. 2

B. Test facilities The fire tests of CU were carried out on KVU-35 and KVU-90 test facilities located at KeRC. Cryogenic vacuum test facility KVU-35 has the main chamber section diameter equal to 3 meters and total volume equal to 35 m 3. It is equipped with cryopumps with total capacity of 42 m 3 /s. The facility is equipped with a control and measurement system, which provides test operation in a fully automatic mode. The main purpose of the facility is carrying out of Hall thrusters lifetime tests. The facility is equipped with a system for remote insulators profile measurement and for video surveillance, which allows monitoring status of the thruster during lifetime tests. The KM-60 thruster, mounted in the vacuum chamber KVU-35 is shown in Fig. 2 while setting diagnostic systems. Figure 2. Thruster KM-60 in the KVU-35 facility while setting diagnostic systems. The vacuum chamber of KVU-90 facility has diameter of 3.8 m and volume of 90 m 3. It is equipped with cryogenic pumps with a total capacity for xenon equal to 140 m 3 /s. The camera is equipped with two prechambers, which allow operational testing and replacement of thrusters without letting-to-air of the main part of chamber. The facility is used for fire tests of Hall and ion thrusters. Appearance of KVU-90 is shown in Fig. 3. Figure 3. Appearance of KVU-90 vacuum chamber. 3

С. Sequence of correction unit development and qualification The initial geometry optimization of the acceleration channel of the discharge chamber, magnetic field topology, parametric tests were performed on a laboratory model of the KM-60. The experience obtained during development and investigation of Hall thrusters at Keldysh Research Centre was used 2. Parameters of the KM-60 laboratory model in the range of discharge voltage from 280 to 680 Volts are presented in Fig. 4 and Fig. 5. Figure 4. The dependence of KM-60 thrust on discharge voltage Figure 5. The dependence of KM-60 specific impulse on discharge voltage Qualification of KM-60 thruster and FCU has been conducted in accordance with a complex program of experimental development. Testing of the KM-60 and FCU has been performed separately with the exception of joint and lifetime tests. In accordance with program the following tests have been performed: autonomous tests of 2 cathode units on the number of cycles (20 000 cycles); the first stage of the CU autonomous tests (two engineering models of KM-60 and FCU); the second stage of the CU autonomous tests (two engineering models of KM-60 and FCU); joint test of CU (KM-60 + FCU) with the laboratory model of the PPCU; qualification tests of CU (three qualifying models of KM-60 and FCU); full-scale lifetime test of CU - 4100 hours. D. Results of qualification Climatic, mechanical, thermal cycling, firing tests and 500-hour lifetime test have been performed at the first stage of autonomous tests. Modifications of the CU construction have been performed at the stage of autonomous tests to eliminate the identified faults and improve specific impulse and lifetime. Operation of thruster after exposure to climatic, mechanical, thermal, transportation loads has been confirmed at the second stage of autonomous tests. Lifetime test of CU with duration equal to 513 hours has been performed. Variation of thrust and specific impulse during this test is shown in Fig. 6, 7. 4

Figure 6. Dependence of thrust on operating time at the second stage of autonomous tests. Figure 7. Dependence of specific impulse on operating time at the second stage of autonomous tests. Parameters of the thruster ion beam before and after lifetime test were investigated during analysis of thruster performance degradation. The structure of anode efficiency 3 defined on the basis of measurement data before and after the lifetime test is shown in Table 2. 5

КМ-60 EM 0 hours КМ-60 EM 513 hours 0.93 0.94 E 0.94 0.92 0.89 0.83 0.73 0.71 Table 2 The dependence of the efficiency components on the operating time with the assumption that the 2 factor varies linearly during lifetime is shown in Fig. 8. One can see from Fig. 8 that the gas efficiency starts to fall steadily after approximately 250 hours of thruster operation, which is the reason for reduction of specific impulse. Since the specific impulse depends only on the value of the product it can be concluded that the main problem is reduction of gas efficiency. The question of the causes of such reduction is analyzed in Ref. 4. Reduction of 2 during the test was not significant and was caused by the degradation of the energy component (lowering of the spectrum mean energy). Figure 8. Changes in the structure of the anode efficiency during lifetime test Modification of thruster design has been made with taking into account of additional tests results. It was based on the results of calculations of magnetic field topology and the position of the ionization zone in the accelerating channel 5,6. Parametric tests of the upgraded KM -60 have shown high initial values of specific impulse (about 2000 seconds). In this regard, it was decided to conduct short lifetime test with equivalent operating time of at least 1400 hours. The short lifetime test has been carried out in several stages according to the method described in Ref. 7. The parameters of the upgraded KM -60 have been measured in the equivalent period of operating time up to 1400 hours. Thrust was equal to (40 ± 1) mn, average specific impulse as the interpolation results for the tests was equal to (1860 ± 50) s. The predicted average specific impulse for the full lifetime was approximately equal to (1800 ± 50) s. The thruster KM -60 was admitted to the qualification tests based on the positive results of the shortened test. Three correction units have been manufactured for qualification tests and they are shown in Fig. 9. Figure 9. Qualification models of correction unit. Parametric, climatic, thermal cycling tests, random and sinusoidal vibration tests, shocks impact, transportation tests as well as measuring the thrust vector have been conducted during qualification of CU. Further each of the three CU has been subjected to the lifetime test with duration not less than 500 hours. The additional correction unit has been manufactured for carrying out of full-scale lifetime test and it also has passed through all qualification procedure just before starting of lifetime test. Thus, 4 KM-60 thrusters and 4 FCU have passed through the qualification test procedure. The results of lifetime tests for all CU were similar so only results of full-scale test are presented below. 6

Significant effect of the atmosphere influence on the specific parameters of the thruster had been observed owing to opening of the vacuum chamber for measurement of the erosion profiles. Therefore full-scale lifetime test has been performed without opening the vacuum chamber. Lifetime tests has been conducted in vacuum chamber KVU -35. Pressure in the chamber was maintained at (7.6 8.4) 10-3 Pa level by air calibrated vacuum gage scale in the course of the test. High purity Xe was used as the propellant. Xenon feed was provided by FCU. Flow measurement was performed using a commercial gas flow regulator in the measurement mode. Measurement error of xenon consumption was not exceeding 2.2%. Power supply of discharge was performed by standard facility power supply unit. During operation, the discharge voltage and discharge current were kept constant and equal to (500 ± 5) V and (1.8 ± 0.05) A correspondingly. Measurement error of the thrust was not exceeding 2.5%.The valves and regulators of FCU were supplied by the facility control system. The functioning of the CU including start procedure was performed in accordance with mission profile. Starts with spare cathode were performed not rare than once per every 100 hours. Initially it had been stated in the test program that the life test was to be divided into two phases. In the first phase it was planned to carry out 2500 hours test without opening the chamber. At least 5500 cycles were to be conducted on this stage. At the second stage, the remaining lifetime (up to 3750 hours) and number of cycles were supposed to be confirmed by shortened tests procedure. Taking into account influence of chamber opening on thruster parameters and difficulties associated with manufacturing with sufficient accuracy of insulators with erosion profiles for 2500 hours and more it was decided to continue the direct lifetime test without opening the vacuum chamber. Intense detachment of films formed on the anode and cylindrical surfaces of insulators was observed after 500 hours of operation. Detached fragments of films overlapped accelerating channel, causing a change in thruster operation, reducing its effectiveness. Reduction of the thrust by 10-20 % from the nominal value was observed at periods 800-1000 and 1200-1780 hours of operation. The most likely cause of trust reduction in these periods of operation is changing of the thruster operation mode due to films detachment and their accumulation in the cavity of the discharge chamber, as well as changes in geometry of the acceleration channel. Operation of thruster was optimized by varying the magnetic field in the acceleration channel during lifetime test. The main purpose of this optimizing was to obtain the maximum value of thrust and specific impulse, stable operation and development of on-board algorithm of magnetic field correction for the flight units. Dependence of current in the coils on operating time is shown in Fig. 10. Figure 10. Dependence of current in the coils of the magnetic system in course of 4100 hours lifetime test. Dependences of thrust and specific impulse on operating time are shown in Fig. 11 and 12 respectively. 7

Figure 11. Dependence of thrust on operating time during 4100 hours lifetime test. Figure 12. Dependence of specific impulse on operating time during 4100 hours lifetime test. Ceramics erosion profiles during lifetime tests were measured by contact method and by laser range profilometer placed directly in the vacuum chamber. Diagnostic system also allows remote monitoring of the discharge chamber elements during pauses of thruster operation, taking photos and videos without opening the vacuum chamber. The results of erosion profiles measurement are shown in Fig. 13. The results of erosion profiles prediction are also shown in this figure. Prediction had been made on the basis of first qualification model lifetime test with total operation time equal to 525 hours. 8

а) b) Figure 13. Profile changing of the internal (a) and external (b) insulators of the discharge chamber during the 4100 hours lifetime test and comparison with the prediction of erosion made on the basis of QM1 525 hours lifetime test. (* - Profiles measured with laser range profilometer) Profiles of the external insulator in different sections by azimuth not differ much after 4100 hours of operation. Four profiles for external insulator are shown in Fig. 14. The good agreement of the insulator erosion profiles indicates the azimuthal uniformity of propellant flow in the accelerating channel. 9

Figure 14. Profiles of the external insulator in different sections after 4100 hours of operation. After erosion of protective coating erosion of magnet pole pieces was started at approximately 1550 hour of operation. The appearance of magnet pole pieces after the lifetime is shown in Fig. 15. Figure 15. Appearance of KM-60 thruster after lifetime test. Figure 15 shows that erosion of the inner pole piece is symmetric but erosion of the outer one is significantly asymmetric with a maximum sputtering level near to the spare cathode. Traces of erosion on the surface of the cathode ignition electrodes, which were under floating potential during test, are virtually absent. The type of magnetic poles erosion shows that it is likely determined by the backward ion flow to the surface. Despite the erosion of the magnetic poles and ceramic insulators thruster was operable and its parameters were in concordance with technical requirements. Thrust vector measurements were conducted before and after 4100 hour lifetime test to determine the deflection angle of the thrust vector from the axis of the thruster, also divergence of the plasma plume and the energy structure of the ions were measured. Comparison of measured thrust vector deflection from the axis of the thruster is shown in Table 3. Measurement error of the thrust vector deflection angle does not exceed 12 angular minutes. 10

Table 3 Cathode Main Spare The deflection angle of the thrust vector before lifetime test, 37.0 38.0 angular minutes The deflection angle of the thrust after before lifetime test, angular minutes 35.7 42.2 After lifetime test the thrust vector deflection angle does not changes considerably. The maximum value of deflection is 42.2 angular minutes. Comparison of plasma plume divergence angle 90 for the main and spare cathodes at distances 700 and 1500 mm from the thruster are shown in table 4. Table 4 700 mm 1500 mm The value of the angular divergence of 90% of the ion current to the beginning of 4100 hours life test, degrees The value of the angular divergence of 90% of the ion current after 4100 hours life test, degrees Main cathode Spare cathode Main cathode Spare cathode 54 54 54 54 43 44 44 47 Thus, the divergence angle of the plasma plume after 4100 hours lifetime test has decreased by 7-11 degrees. Ion current density on the thruster axis has increased about 1.5 times. Analysis of ion energy distributions obtained at a distance of 1500 mm from the thruster have shown narrowing of the ion energy spectrum with increasing their average energy at small angles to the axis of the thruster and reduction of high-energy ions number with a drop of their average energy at the angles of 40-50 degrees from the thruster axis. III Conclusions Qualification of correction unit consists of Hall thruster KM-60 and flow control unit has been successfully completed. Lifetime tests of the correction unit were conducted during qualification. The maximum operation time of single correction unit is equal to 4120 hours at 8357 cycles, the total shown thrust impulse is 602 kn s with an average specific impulse -1862 seconds, thrust at the end of the test - 41.8 mn. At the end of lifetime test correction unit remained operable, and its parameters were in concordance with the specifications. References 1 O.Gorshkov, M.Belikov, V.Muravlev, A.Shagayda, K.Shanbhogue, S.Pre4mkumar, M.Jayaraman, V.Nandalan Development and qualification test of Electric Propulsion System based on the HET KM-45 for GSAT-4 Spacecraft, IAC-07-C4.4.01, 58 th IAC, Hyderabad, India, September 24-28, 2007. 2 O.A.Gorshkov, V.A. Muravlev, A.A. Shagayda Hall and Ion Electric Propulsion for Spacecraft, Mashinostroenie, Moscow, 2008 (in Russian). 3 O. A. Gorshkov and A. A. Shagayda. Determining the Efficiency of a Plasma Thruster with Closed Electron Drift // Technical Physics Letters, 2008, Vol. 34, No. 2, pp. 153 155. 4 A. A. Shagaida, O. A. Gorshkov and D. A. Tomilin. Influence of the erosion of the discharge channel wall on the efficiency of a stationary plasma thruster // Technical Physics, 2012, Vol. 57, No. 8, pp. 1083-1089. 5 N.Blinov, O.Gorshkov, A.Shagayda Experimental investigation of magnetic field topology influence on structure of acceleration layer and performance of hall thrusters, IEPC-2005-033, 29 International Electric Propulsion Conference, Princeton (USA), October 30 - November 04, 2005. 6 O. A. Gorshkov and A. A. Shagaida. The Criterion of Optimal Configuration of Magnetic Field in a Thruster with Closed Electron Drift // High Temperature, 2008, Vol. 46, No. 4, pp. 529 534. 7 Lovtsov A., Shagayda A. and Gorshkov O. Semi-Empirical Method of Hall Thrusters Lifetime Prediction, 42 nd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, AIAA Paper 2006-4661. 11