Magnetic and Langmuir Probe Measurements on the Plasmoid Thruster Experiment (PTX)

Summary

This paper presents magnetic field and plasma diagnostic measurements performed on the Plasmoid Thruster Experiment (PTX), a pulsed inductive accelerator utilizing a conical theta-pinch coil. An array of B-dot magnetic probes and a quadruple Langmuir probe were constructed and utilized to determine magnetic field structures, electron number density, electron temperature, and bulk plasma flow velocity ratios for both argon and hydrogen propellants. Experimental results confirmed the formation and acceleration of plasmoids with characteristic magnetic field reversals, indicating the viability of inductive plasma acceleration without electrode erosion.

Title and Abstract

Magnetic and Langmuir Probe Measurements on the Plasmoid Thruster Experiment (PTX)

Syri J. Koelfgen* University of Alabama in Huntsville (UAH), Propulsion Research Center, S225 Technology Hall, Huntsville, AL 35899

Richard Eskridge,† Michael H. Lee,‡ and Adam Martin§ NASA Marshall Space Flight Center (MSFC), Propulsion Research Center TD40, Huntsville, AL 35812

Clark W. Hawk,¶ and Peter Fimognari# University of Alabama in Huntsville (UAH), Propulsion Research Center, S225 Technology Hall, Huntsville, AL 35899

The Plasmoid Thruster Experiment (PTX) operates by inductively producing plasmoids in a conical theta-pinch coil and ejecting them at high velocity. A plasmoid is a plasma with an imbedded closed magnetic field structure. The shape and magnetic field structure of the translating plasmoids have been measured with an array of magnetic field probes. Six sets of two B-dot probes were constructed for measuring Bz and Bθ, the axial and azimuthal components of the magnetic field. The probes are wound on a square G-10 form, and have an average (calibrated) NA of 9.37 x 10^-5 m^2, where N is the number of turns and A is the cross-sectional area. The probes were calibrated with a Helmholtz coil, driven by a high-voltage pulser to measure NA, and by a signal generator to determine the probe’s frequency response. The plasmoid electron number density ne, electron temperature Te, and velocity ratio v/cm (where v is the bulk plasma flow velocity and cm is the ion thermal speed) have also been measured with a quadruple Langmuir probe. The Langmuir probe tips are 10 mm long, 20-mil diameter stainless steel wire, housed in a 6-inch long 4-bore alumina rod. Measurements on PTX with argon and hydrogen from the magnetic field probes and quadruple Langmuir probe will be presented in this paper.

Nomenclature

Nomenclature A = cross-sectional area of B-dot probes A1, A2, A3, A4 = collection area of quadruple Langmuir probe electrodes 1, 2, 3, and 4 B, B = magnetic field Bz = axial component of magnetic field Bθ = azimuthal component of magnetic field cm = ion thermal speed ec = electron charge eeλmfp = electron-electron collision mean free path iiλmfp = ion-ion collision mean free path I1, I2, I3, I4 = current collected by quadruple Langmuir probe electrodes 1, 2, 3, and 4 IB = Bohm current It = current collected at Langmuir probe tip J = current density Je = electron saturation current density Ji = ion saturation current density L = Langmuir probe electrode length λD = Debye length me = electron mass Mi = ion mass ne = electron number density N = number of turns of B-dot probe rp = Langmuir probe electrode radius s = clearance between Langmuir probe electrodes t = time t0 = initial time Te = electron temperature in eV v = bulk plasma velocity vBohm = Bohm velocity V = voltage induced in B-dot probe coil V1, V2, V3, V4 = potential of Langmuir probe electrodes 1, 2, 3, and 4 Vd2 = potential difference between probes 2 and 1 Vd3 = applied potential difference between probes 3 and 1 Vd4 = applied potential difference between probes 4 and 1 Vout = B-dot probe passive integrator output signal Vp = plasma potential

  • NASA/MSFC GSRP Fellow and UAH Graduate Research Assistant, AIAA Student Member. † NASA/MSFC Aerospace Systems Engineer, AIAA Member. ‡ NASA/MSFC Electrical Engineer. § NASA/MSFC Physicist. ¶ UAH Propulsion Research Center Director and Mechanical & Aerospace Engineering Professor, AIAA Fellow.

UAH Graduate Research Assistant.

I. Introduction

I. Introduction THE plasmoid thruster is a pulsed inductive rocket that utilizes the J x B (Lorentz) force to accelerate plasmoids and produce thrust (Fig. 1). A plasmoid is a plasma with an internal magnetic field structure, also known as a compact toroid.1 Plasmoid accelerators have been used for fueling fusion devices2-4 and their application to space propulsion has been proposed and studied theoretically.5-7 The Plasmoid Thruster Experiment (PTX) is an experiment to evaluate the use of plasmoids in propulsion.8 The plasmoids are formed inductively in a conical theta-pinch coil, which eliminates the problems associated with discharging a large current across electrodes to produce and accelerate the plasma. Electrode erosion is a limiting factor in the functional lifetime of many pulsed electromagnetic accelerators, such as the magnetoplasmadynamic (MPD) thruster.5,9 The plasmoid thruster is an option for removing this limitation. Additionally, since the plasmoid has a self-contained magnetic field and is therefore not connected to any external magnetic field lines, the problem of detachment is reduced. Magnetic insulation of the plasma should also lead to reduced thermal losses to the walls. An actual plasmoid thruster would operate repetitively at a frequency of 10-100 Hz and would use solid-state switching. In the present experiment, only a single plasmoid is formed and accelerated; PTX is not yet optimized.

[Figure 1. Plasmoid in Conical θ-Pinch Coil]

II. PTX

II. PTX In PTX, the plasmoids are produced in a Pyrex tube, which serves as the vacuum chamber, situated inside of a single-turn conical theta-pinch coil with a half-angle of 17.5° (Fig. 2). The Pyrex tube is connected to a rectangular vacuum chamber for diagnosis of the exhaust plume. The theta-pinch coil is driven by a 560 nF, 40 kV capacitor bank that is switched by a Perkin-Elmer GP-32B spark-gap switch. The coil and capacitor bank constitute a tank (LC) circuit that rings sinusoidally with a peak current of 53 kA at 35 kV charge voltage (Fig. 3). The bank is fired by sending an optical trigger signal to a high-voltage pulser that then triggers the spark-gap switch. When triggered, the spark-gap closes, allowing the capacitor bank to discharge through the theta-pinch coil. A solenoid valve injects propellant into the chamber, before the bank is discharged at a time t0. (The time before the bank discharge at which propellant is injected is referred to as the puff valve delay.) The rapidly changing axial magnetic field creates an azimuthal electric field that preionizes the plasma. Simultaneously, the plasma is seeded with the magnetic field produced during the first half-cycle of the sinusoidal discharge. When the current swings through zero, the field is reversed and the chamber contains a cold, partially ionized plasma with an imbedded magnetic field. As the external magnetic field increases in the opposite direction, it reconnects with the bias field and the plasma is compressed and fully ionized.1 During the formation process, the interaction of the magnetic field with the plasma induces a large azimuthal current in the plasma. The resultant J x B force accelerates the plasmoid away from the coil and generates thrust. Argon and hydrogen gases were chosen to model heavyweight and lightweight propellant in PTX.

[Figure 2. Plasmoid Thruster Experiment] [Figure 3. PTX Capacitor Bank Discharge]

III. Diagnostics & A. Internal B-dot Probe Array

III. Diagnostics A variety of diagnostics are used on PTX, including: a. a high-speed Cordin 220B framing camera (with eight independently triggered CCD arrays; inter-frame time as low as 10 ns), for imaging and estimating the plasmoid velocity; b. an excluded flux array (an array of magnetic field probes and flux loops), for determining the shape of the plasmoid in the coil; c. a light pipe velocimeter (fiber optic cables fed into fast PIN photo diodes), for measuring velocity; d. a heterodyne, quadrature HeNe laser interferometer, for measuring the line-averaged electron density, ne; e. an array of magnetic field sensors (B-dot probes), for measuring axial and azimuthal magnetic field; and f. a quadruple Langmuir probe, for measuring the plasmoid electron temperature Te, number density ne, and velocity ratio v/cm, where v is the bulk plasma flow velocity and cm is the ion thermal speed. This paper will focus on the B-dot probe array and quadruple Langmuir probe.

A. Internal B-dot Probe Array An array of internal B-dot probes10,11 has been constructed and installed in the exhaust chamber, perpendicular to the z-axis (Fig. 4). Six sets of two B-dot probes measure the axial and azimuthal magnetic fields (Bz and Bθ) of the plasmoid as it translates past the probe array. The probes were wound with 36 gauge magnet wire on a 3.15 mm square G-10 form (Fig. 4). The probe signals are passively integrated (RC ≈ 20 μs) for improved dynamic range. The output signals from the passive integrator, Vout, are recorded with Acqiris digitizers (8-bit, 250 Msamples/s). The magnetic field measured by each probe is given by:

B = ∫ B-dot dt = (1 / NA) ∫ V dt = (1 / NA) (RC Vout + ∫ Vout dt) (1)

[Figure 4. Internal B-dot Probe Array: a) Internal B-dot Probe Array on PTX, b) Internal B-dot Probe Array]

The B-dot probes were calibrated with a custom-built Helmholtz coil (Fig. 5). The Helmholtz coil is a set of two coaxial current loops with equal radii, separated from each other by a distance equal to that radius. The radius and separation of the loops is approximately 2.5 inches. Equal amounts of current flowing in the same direction are applied to the loops, producing an extraordinarily uniform magnetic field inside the volume of the Helmholtz coil. The PTX Helmholtz coil was calibrated with a DC power supply, multimeter, and Hall probe to determine the magnetic field inside the Helmholtz coil as a function of applied current. To calibrate the individual B-dot probes in the array, the Helmholtz coil was driven with a 300 Volt pulser, and the resulting voltage responses of the Helmholtz coil and each B-dot probe were recorded with a digital oscilloscope. The calibrated quantity, NA, of the probes was found to be within 20% of the design value of NA = 10 x (3.15 mm)^2 = 9.92 x 10^-5 m^2. The frequency response of the B-dot probes was also assessed by applying a sine wave to the Helmholtz coil at frequencies ranging from 1 kHz to 12 MHz. The resulting Helmholtz coil signal (obtained with a Pearson current probe) and the B-dot probe responses were measured with the oscilloscope and combined to determine NA as a function of frequency.

[Figure 5. Helmholtz Coil]

III. Diagnostics - B. Quadruple Langmuir Probe

B. Quadruple Langmuir Probe A quadruple Langmuir probe12 has also been constructed (Fig. 6) to make local measurements of the plasmoid electron temperature Te, electron number density ne, and velocity ratio v/cm. Three of the four probes are aligned parallel to the plasma flow (the z-axis); the fourth probe is aligned perpendicular. The Langmuir probe tips are 10 mm long, 0.5 mm (20 mil) diameter stainless steel wires, housed in a 152 mm long, 4-bore alumina rod. The alumina rod is torr-sealed into a stainless steel tube and installed at the downstream end of the PTX exhaust chamber through a sliding O-ring seal, so that measurements may be taken along the z-axis.

[Figure 6. Quadruple Langmuir Probe]

An electrical schematic of the quadruple Langmuir probe is shown in Fig. 7. Probe 2 is unbiased, and so it floats to a potential such that it does not collect any current. Probes 3 and 4 are biased at approximately equal voltages (Vd3 and Vd4) relative to probe 1. Probe 1 is biased above the floating potential and below electron saturation, and emits current. Probes 3 and 4 are biased in the ion saturation region of the ideal probe characteristic, and collect current. For both probes 3 and 4, the bias voltage is produced by a variable DC power supply, with 20 μF of capacitance in parallel, in order to source the current on the fast time scale of the experiment.

[Figure 7. Quadruple Langmuir Probe Circuit]

Quadruple Langmuir probe theory12,13 is a combination of triple probe theory14,15 and a crossed electrostatic probe technique.16,17 The measured quantities for data analysis are the currents collected by probes 3 and 4 (I3 and I4), and the potential differences of probes 2, 3, and 4 relative to probe 1, called Vd2, Vd3, and Vd4 respectively, where

Vd2 = V2 - V1 (2)

and

Vd3 = V3 - V1 Vd4 = V4 - V1 (3)(a-c) Vd3 = Vd4

Simultaneous collection of I3, I4, Vd2, Vd3, and Vd4 yield measurements of Te, ne, and v/cm. Analysis of the probe response assumes that several conditions are met, including: a. the probe is operating in the collisionless sheath regime, i.e. the condition eeλmfp / λD >> 1 is met, where eeλmfp is the electron-electron collision mean free path and λD is the Debye length; b. the sheath around each probe electrode is thin, as prescribed by rp / λD >> 1, where rp is the Langmuir probe electrode radius; c. the probe is in the free molecular regime, so the effect of free-streaming particle collisions with the probe can be ignored, such that eeλmfp / rp > 1 and iiλmfp / rp > 1 are satisfied, where iiλmfp is the ion-ion collision mean free path; d. the end effect is minimized, which is achieved if the plasma velocity v (the heavy-particle velocity) is small compared to the Bohm velocity vBohm, so that there is not a large number of ions reaching the probe causing a peak in the measured ion current, such that the condition L > 50 λD (v / vBohm) is satisfied, where L is the length of the probe tip; e. the error introduced to Langmuir probe current measurements by probe tip currents is less than approximately 10%, i.e. It / IB = (1 / 2) (rp / L) (v / vBohm) ≤ 0.1; and f. the clearance s between probe electrodes is sufficient so that sheath interactions between adjacent electrodes is avoided, such that s / λD > 250. The PTX quadruple Langmuir probe was designed to ensure that these conditions were adequately met under expected operating conditions.

III. Diagnostics - Langmuir Probe Equations (Electron Temperature, Density, Velocity)

  1. Electron Temperature This paper takes a slightly different approach to the derivations of the formulas for Te and ne, however the results are equivalent to those found in other works.12,13 To determine an expression for the electron temperature from the quadruple Langmuir probe measurements, we begin with equations describing the current collected by probes 1, 2, and 3:

I1 = Je A1 exp((V1 - Vp) / Te) - Ji A1 (4) I2 = Je A2 exp((V2 - Vp) / Te) - Ji A2 (5) -I3 = Je A3 exp((V3 - Vp) / Te) - Ji A3 (6)

where A1, A2, and A3 are the current collection areas of probes 1, 2, and 3, Vp is the plasma potential, Te is the electron temperature in units of electron volts (eV), and Je and Ji are the electron and ion saturation current densities, as given by:

Je = ne ec^(3/2) sqrt(Te / (2 π me)) (7) Ji = exp(-1/2) ne ec^(3/2) sqrt(Te / Mi) (8)

where ne is the electron number density, ec is the electron charge, me is the electron mass, and Mi is the ion mass. Since the current emitted by probe 1 is collected by both probes 3 and 4, and since I1 flows in the opposite direction of I3 and I4 (as shown in Fig. 7),

I1 = -(I3 + I4) = -I3 (1 + I4 / I3). (9)

Given that A1 = A2 = A3, and substituting Eq. (9) into Eq. (4), we have:

I1 = -Je A3 [ exp((V3 - Vp) / Te) - Ji / Je ] (1 + I4 / I3). (10)

Since probe 2 is floating and therefore I2 = 0, Eq. (5) reduces to

Ji / Je = exp((V2 - Vp) / Te) (11)

Equating Eq. (4) and Eq. (10), noting that A1 = A3, dividing by Je, and applying Eq. (11) yields:

exp((V1 - Vp) / Te) - exp((V2 - Vp) / Te) = -(1 + I4 / I3) [ exp((V3 - Vp) / Te) - exp((V2 - Vp) / Te) ]. (12)

Multiplying Eq. (12) by -exp(Vp / Te) and exp(-V1 / Te), and applying Eq. (2) and Eq. (3a), we are left with an expression for Te (in eV) as a function of the measured quantities I3, I4, and Vd2 as well as the applied voltage Vd3:

exp(Vd2 / Te) - 1 = (1 + I4 / I3) [ exp(Vd3 / Te) - exp(Vd2 / Te) ]. (13)

  1. Electron Density To derive an expression for the electron density, Eq. (6) is represented as:

-I3 = Je A3 exp(V3 / Te) exp(-Vp / Te) - Ji A3. (14)

Rearranging Eq. (11) into

exp(-Vp / Te) = (Ji / Je) exp(-V2 / Te), (15)

substituting this into Eq. (14), and using Eq. (2) and Eq. (3a), we are left with:

-I3 = Ji A3 [ exp((Vd3 - Vd2) / Te) - 1 ]. (16)

Substituting Eq. (8) into Eq. (16), and rearranging, yields the expression for ne:

ne = -I3 / [ exp(-1/2) ec^(3/2) A3 sqrt(Te / Mi) [ exp((Vd3 - Vd2) / Te) - 1 ] ]. (17)

  1. Velocity In addition to measurements of Te and ne, the quadruple Langmuir probe also provides a measurement of the velocity ratio v/cm from the collected currents I3 and I4. For Vd3 = Vd4, with probe 3 aligned with the flow vector and probe 4 perpendicular to it, the ratio of I3 to I4 has been shown to be17

I4 / I3 = (2 A4 / (sqrt(π) A3)) exp(-(v / cm)^2) [ sum_{n=0}^{inf} ( (v / cm)^(2n) / n! ) Γ(n + 3/2) ]^(1/2). (19)

The velocity ratio v/cm is obtained by iteration of Eq. (19).

IV. Experimental Results

IV. Experimental Results The B-dot probe array is mounted on the top of the PTX exhaust chamber, so the plasmoid magnetic field can be measured at various vertical positions in the chamber (Fig. 8). In Fig. 8(a), probe B_z3, the third axial B-dot probe in the array, is positioned on the centerline of the theta-pinch coil and chamber. Probe B_z1 is in the lowest vertical position (below the centerline) and B_z6 is in the highest vertical position. In Fig. 8(b), probe B_z5 is on the centerline.

[Figure 8. B-dot Probe Array Positions in Relation to the PTX Centerline (CL): (a) Probe B_z3 on CL, (b) Probe B_z5 on CL]

For the same operating conditions (argon at 6.5 psig manifold pressure injected at a 2600 μs puff valve delay), PTX was run with the B-dot probe array in the two different vertical positions on two different days. This was done so that the plasmoid magnetic fields could be mapped over the entire vertical cross-section of the chamber. On February 23, 2004 Shot 3, the B-dot probe array positioning shown in Fig. 8(a) was used, and on February 25, 2004 Shot 43, the configuration in Fig. 8(b) was used. The axial and azimuthal magnetic fields as a function of time were obtained for both shots. The magnetic field averaged over several data points about 18.65 μs (the time near the axial magnetic field peak) was then determined for each B-dot probe, so that the magnetic field could be graphed as a function of vertical position, as each B-dot probe corresponds to a different vertical position in the chamber. The axial magnetic field profile vs. vertical position for the two shots at 18.65 μs is shown in Fig. 9(a). A Cordin camera photo of the plasmoid translating toward the B-dot probe array is shown in Fig. 9(b). The light pipe velocimeter measured plasmoid velocities of 28 km/s on February 23, 2004 Shot 3 and 22 km/s on February 25, 2004 Shot 43. For the same operating conditions, the quadruple Langmuir probe was also used to measure Te, ne, and v/cm as a function of time (Fig. 10).

[Figure 9. Plasmoid Magnetic Field vs. Vertical Position and Cordin Camera Photograph, 6.5 psig Argon at 2600 μs delay: (a) Bz vs. Vertical Position, (b) Cordin Photo from 2-25-04 Shot 43] [Figure 10. Plasmoid Temperature and Density, May 19, 2004, Shot 7, 6.5 psig Argon at 2600 μs delay]

Figure 9(a) is consistent with the axial magnetic fields expected in a plasmoid. Figure 10 shows a maximum electron temperature of 7.6 eV and a maximum electron number density of 4.9 x 10^19 m^-3. At 18.65 μs, the temperature is 6.3 eV and the density is 4.0 x 10^19 m^-3. In determining the velocity ratio v/cm, the summation term in Eq. (19) reaches a constant value after 50 terms, so for n = 0..50, the peak velocity ratio for this test was approximately 3.2. For comparison to argon, hydrogen was used as a lightweight propellant. The axial and azimuthal magnetic fields in PTX for hydrogen at 38.6 psig manifold pressure, injected 2400 μs before the bank discharge, are shown in Fig. 11. Numbers corresponding to each B-dot probe are placed on their corresponding curve, where “1” denotes the magnetic field curve obtained from B-dot probe z1 or θ1. In Fig. 11 we see that the azimuthal magnetic field, Bθ, reverses direction, which is consistent with the field-reversal characteristic of plasmoids. The noise in the beginning of the magnetic field plots is pick-up from the discharge of the capacitor bank. Figure 12 shows the electron temperature and density profiles for similar operating conditions. The temperatures at the first and second peaks are 22 eV and 23 eV. The corresponding densities at the two peaks are 1.1 x 10^20 m^-3 and 1.2 x 10^20 m^-3. The velocity ratio v/cm on the first peak in Fig. 12 is approximately 0.4.

[Figure 11. Plasmoid Magnetic Fields, May 18, 2004, Shot 28, 38.6 psig Hydrogen at 2400 μs delay] [Figure 12. Plasmoid Temperature and Density, June 14, 2004, Shot 3, 38.7 psig Hydrogen at 2400 μs delay]

V. Discussion, Acknowledgments, and References

V. Discussion The magnetic field plots obtained from the B-dot probe array are consistent with plasmoid magnetic field configurations,1 so PTX appears to be producing plasmoids. In the tests presented, the peak electron temperature, electron number density, and velocity ratio were 7.6 eV, 4.9 x 10^19 m^-3, and 3.2 for argon, and 23 eV, 1.2 x 10^20 m^-3, and 0.4 for hydrogen. Hydrogen produced higher temperatures and densities, and lower v/cm, than argon, as expected with a lighter gas. The two peaks in the hydrogen temperature and density plots may indicate that two plasmoids are being formed during a single shot. In addition, some of the azimuthal B-dot probes show magnetic fields changing direction twice, which may indicate that two plasmoids are being formed in successive cycles of the ringing discharge. A large volume of PTX data, which is sufficiently reproducible for similar operating conditions, has been collected and is being analyzed.

Acknowledgments Syri Koelfgen is funded by the NASA/MSFC Graduate Student Researchers Program (GSRP) fellowship. This research was conducted in the NASA/MSFC Propulsion Research Center (TD40) laboratories. The authors would like to thank Jeff Richeson, Tommy Reid, and Doug Galloway of Mainthia for their technical support of PTX.

References 1 Tuszewski, M., “Field Reversed Configurations – Review Paper,” Nuclear Fusion, Vol. 28, No. 11, 1988, pp. 2033-2092. 2 Hoffman, A. L., et al., “Inductive Field-Reversed Configuration Accelerator for Tokamak Fueling,” Fusion Technology, Vol. 36, Sept. 1999. 3 Miley, G. H., “Compact tori for alternate fuel fusion,” Nuclear Instruments and Methods in Physics Research, Vol. 207, No. 1-2, 15 March 1983, pp. 111-120. 4 Kawai, K., Pietrzyk, Z. A., and Hunter, H. T., “Generation of poloidally rotating spheromaks by the conical theta pinch,” Phys. Fluids, Vol. 30, No. 8, Aug. 1987, pp. 2561-2568. 5 Jahn, R. G., Physics of Electric Propulsion, McGraw-Hill, 1968. 6 Bourque, R. F., and Parks, P. B., “Time-Dependent Analysis of the Pulsed Plasmoid Electric Thruster,” AIAA/NASA/OAI Conference on Advanced SEI Technologies, AIAA Paper 91-3614, Cleveland, OH, 4-6 Sept. 1991. 7 Slough, J., “Propagating Magnetic Wave Plasma Accelerator (PMWAC) for Deep Space Exploration,” STTR Phase-1 Final Report, 2000. 8 Koelfgen, S. J., Eskridge, R., Hawk, C. W., Lee, M. H., Martin, A., and Smith, J. W., “A Plasmoid Thruster for Space Propulsion,” 39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, AIAA Paper 2003-4992, Huntsville, AL, 20-23 July 2003. 9 Myers, R. M., et al., “Cathode Phenomena in a Low-Power Magnetoplasmadynamic Thruster,” Journal of Propulsion and Power, Vol. 7, No. 5, Sept.-Oct. 1991, pp. 760-766. 10 Hutchinson, I. H., Principles of Plasma Diagnostics, Cambridge University Press, 1987. 11 Huddlestone, R. H., and Leonard, S. L., Plasma Diagnostic Techniques, Academic Press, 1965. 12 Burton, R. L., DelMedico, S. G., and Andrews, J. C., “Application of a Quadruple Probe Technique to MPD Thruster Plume Measurements,” Journal of Propulsion and Power, Vol. 9, No. 5, Sept.-Oct. 1995, pp. 771-777. 13 Bufton, S. A., “Exit Plane Plasma Measurements of a Low-Power Hydrazine Arcjet,” Ph.D. Dissertation, Univ. of Illinois, Urbana-Champaign, 1996. 14 Tilley, D. L., Kelly, A. J., and Jahn, R. G., “The Application of the Triple Probe Method to MPD Thruster Plumes,” 21st AIAA/DGLR/JSASS International Electric Propulsion Conference, AIAA Paper 90-2667, Orlando, FL, 18-20 July 1990. 15 Chen, S. L., and Sekiguchi, T., “Instantaneous Direct-Display System of Plasma Parameters by Means of Triple Probe,” Journal of Applied Physics, Vol. 36, No. 8, Aug. 1965, pp. 2363-2375. 16 Kanal, M., “Theory of Current Collection of Moving Cylindrical Probes,” Journal of Applied Physics, Vol. 35, No. 6, June 1964, pp. 1697-1703. 17 Johnson, B. H., and Murphree, D. L., “Plasma Velocity Determination by Electrostatic Probes,” AIAA Journal, Vol. 7, Oct. 1969, pp. 2028-2030.