The new visualization of applied-eld magnetoplasmadynamic thruster operation offered by numerical simulations with the MACH2 code is used to develop an analytic model. This model provides simple expressions for the thrust and voltage drop across the plasma that capture all of the trends and in most cases the magnitudes of a wide range of experimental data including tests at multiple laboratories on three propellants: argon, hydrogen, and lithium. It establishes that the thrust scales as the square root of the product of the applied-eld strength, discharge current, and mass-ow rate. Voltage scales linearly with applied-eld strength and is independent of discharge current and mass-ow rate. This model also provides additional insights that can be used to improve applied-eld magnetoplasmadynamicthruster performance.
NomenclatureA = atomic weight, amu a = cathode radius to electrode length ratio, r c / L B = magnetic induction, T I = discharge current, A j = current density, A/m 2 L = electrode length, m Çm = mass-ow rate, kg/s p = pressure, N/m 2 Q = cross section, m 2 R = electrode radius ratio, r a / r c R = speci c gas constant, m 2 /s 2 -K Rm = magnetic Reynolds number r a = anode radius, m r c = cathode radius, m T H = heavy particle temperature, K U = characteristic speed, IB/ l , m/s V = discharge voltage, V V emf = electromotive voltage, V V f = fall voltage, V V p = plasma voltage, V V RH = resistive + Hall voltage, V v = axial velocity, m/s w = azimuthal velocity, m/s Z = charge number a = degree of ionization g F = ow ef ciency k = mean free path, m l = viscosity coef cient, kg/m-s q = mass density, kg/m 3 ¾ = stress tensor, N/m 2 } = ionization factor, i [a i / j a j Z 2 i Z 2 j ] X = Hall parameter Subscripts r, h , z = physical dimensions: radial, azimuthal, axial