A flux‐limited numerical method for solving the MHD equations to simulate propulsive plasma flows

For numerical simulations to be effective tools in plasma propulsion research, a high-order accurate solver that captures MHD shocks monotonically and works reliably for strong magnetic fields is needed. For this purpose, a characteristics-based scheme for the MHD equations, with flux limiters to improve spatial accuracy, has been developed. In this method, the symmetric form of the MHD equations, accounting for waves propagating in all directions, are solved. The required eigensystem of axisymmetric MHD equations, with appropriate normalization, is presented. This scheme was validated with unsteady (Riemann problem) and force-free equilibrium (Taylor state) test cases, as well as with measured current density patterns in a magnetoplasmadynamic thruster. Copyright © 2001 John Wiley & Sons, Ltd.

[1]  M. Lapointe,et al.  Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters , 1992 .

[2]  P. Roe,et al.  A Solution-Adaptive Upwind Scheme for Ideal Magnetohydrodynamics , 1999 .

[3]  A. Jameson.  ANALYSIS AND DESIGN OF NUMERICAL SCHEMES FOR GAS DYNAMICS, 1: ARTIFICIAL DIFFUSION, UPWIND BIASING, LIMITERS AND THEIR EFFECT ON ACCURACY AND MULTIGRID CONVERGENCE , 1995 .

[4]  M. Brio,et al.  An upwind differencing scheme for the equations of ideal magnetohydrodynamics , 1988 .

[5]  N. Aslan.  TWO‐DIMENSIONAL SOLUTIONS OF MHD EQUATIONS WITH AN ADAPTED ROE METHOD , 1996 .

[6]  Edgar Y. Choueiri,et al.  Scaling of Thrust in Self-Field Magnetoplasmadynamic Thrusters , 1998 .

[7]  J. B. Taylor,et al.  Relaxation of toroidal plasma and generation of reverse magnetic fields , 1974 .

[8]  Philip L. Roe,et al.  An upwind scheme for magnetohydrodynamics , 1995 .

[9]  Luigi Martinelli,et al.  Calculations of viscous flows with a multigrid method , 1987 .

[10]  Uri Shumlak,et al.  An Implicit Scheme for Nonideal Magnetohydrodynamics , 1997 .

[11]  T. Miyasaka,et al.  AMA 99-2432 Numerical Prediction of Onset Phenomenon in a 2-Dimensional Axisymmetric MPD Thruster , 1999 .

[12]  Gérard Gallice,et al.  Roe Matrices for Ideal MHD and Systematic Construction of Roe Matrices for Systems of Conservation Laws , 1997 .

[13]  J. Brackbill,et al.  The Effect of Nonzero ∇ · B on the numerical solution of the magnetohydrodynamic equations☆ , 1980 .

[14]  Kenneth G. Powell,et al.  AN APPROXIMATE RIEMANN SOLVER FOR MAGNETOHYDRODYNAMICS (That Works in More than One Dimension) , 1994 .

[15]  Bram van Leer,et al.  Design of Optimally Smoothing Multi-Stage Schemes for the Euler Equations , 1989 .

[16]  E. Choueiri.  Anomalous resistivity and heating in current-driven plasma thrusters * , 1999 .

[17]  Eliahu Haym Niewood,et al.  An explanation for anode voltage drops in an MPD thruster , 1993 .

[18]  R. Jahn,et al.  Physics of Electric Propulsion , 1968 .

[19]  P. Roe.  CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .

[20]  M. Martinez-Sanchez,et al.  Two-dimensional numerical simulation of MPD flows , 1987 .

[21]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[22]  H. Seifert,et al.  Rocket Propulsion Elements , 1963 .

[23]  D Lileikis,et al.  Effects of azimuthal propellant injection asymmetry on MPD thruster performance using the Mach3 code , 1995 .

[24]  Monika Auweter-Kurtz,et al.  Numerical Modeling of the Flow Discharge in MPD Thrusters , 1987 .

[25]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[26]  Culbert B. Laney,et al.  Computational Gasdynamics: Waves , 1998 .

[27]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .