Numerical simulation of MPD viscous flows

Magneto-plasma axisymmetric flows have been calculated by solving plasma flow equations and the equation of induction derived from the Maxwell's equation and the Ohm's law. In most of existing methods, the equation of induction was solved separately by a relaxation method as an elliptic equation, because the characteristic time of magnetic field is much shorter than that of flow field.
In the present study, the equation of induction is simultaneously calculated with the plasma flow equations. The diffusion term in the equation of induction is calculated by a point-implicit method coupled with the LU-SGS scheme to satisfy the convergence of the solution(Figurej.

References

Satoru YAMAMOTO and Fumito CHIBA, An Efficient Implicit Scheme of MPD Thruster Viscous Flow, Trans. JSME, 66(2000), 647-654(in Japanese).

Hironori Takeda and Satoru Yamamoto, Implicit Time-marching Solution of Partially Ionized Flows in Self-Field MPD Thruster, Transactions of the Japan Society for Aeronautical and Space Sciences, 44-146 (2002) 223-228.

Hironori TAKEDA and Satoru YAMAMOTO, Numerical Investigation of Supersonic MPD Viscous Flows with Ionization, JSME International Journal, 45(2002), 97-101.