Robust Finite Volume Schemes for Two-Fluid Plasma Equations

Two-fluid plasma equations are derived by taking moments of Boltzmann equations. Ignoring collisions and viscous terms and assuming local thermodynamic equilibrium we get five moment equations for each species (electrons and ions), known as two-fluid plasma equations. These equations allow different temperatures and velocities for electrons and ions, unlike ideal magnetohydrodynamics equations. In this article, we present robust second order MUSCL schemes for two-fluid plasma equations based on Strang splitting of the flux and source terms. The source is treated both explicitly and implicitly. These schemes are shown to preserve positivity of the pressure and density. In the case of explicit treatment of source term, we derive explicit condition on the time step for it to be positivity preserving. The implicit treatment of the source term is shown to preserve positivity, unconditionally. Numerical experiments are presented to demonstrate the robustness and efficiency of these schemes.

[1]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[2]  Derek M. Causon,et al.  On the Choice of Wavespeeds for the HLLC Riemann Solver , 1997, SIAM J. Sci. Comput..

[3]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[4]  Claus-Dieter Munz,et al.  Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model , 2000 .

[5]  S. Baboolal.  Finite-difference modeling of solitons induced by a density hump in a plasma multi-fluid , 2001 .

[6]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[7]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[8]  Uri Shumlak,et al.  Approximate Riemann solver for the two-fluid plasma model , 2003 .

[9]  F. Bouchut.  Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws: and Well-Balanced Schemes for Sources , 2005 .

[10]  J. Goedbloed,et al.  Principles of Magnetohydrodynamics , 2004 .

[11]  C. Berthon,et al.  Stability of the MUSCL Schemes for the Euler Equations , 2005 .

[12]  Uri Shumlak,et al.  A high resolution wave propagation scheme for ideal Two-Fluid plasma equations , 2006, J. Comput. Phys..

[13]  Christophe Berthon,et al.  Why the MUSCL–Hancock Scheme is L1-stable , 2006, Numerische Mathematik.

[14]  S. Baboolal,et al.  Two-scale numerical solution of the electromagnetic two-fluid plasma-Maxwell equations: Shock and soliton simulation , 2007, Math. Comput. Simul..

[15]  E. A. Johnson,et al.  Collisionless Magnetic Reconnection in a Five-Moment Two-Fluid Electron-Positron Plasma , 2008 .

[16]  S. Baboolal,et al.  High-resolution numerical simulation of 2D nonlinear wave structures in electromagnetic fluids with absorbing boundary conditions , 2010, J. Comput. Appl. Math..

[17]  Uri Shumlak,et al.  A Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations , 2010, 1003.4542.

[18]  Uri Shumlak,et al.  Numerical Methods for Two-Fluid Dispersive Fast MHD Phenomena , 2011 .

[19]  Harish Kumar,et al.  Entropy Stable Numerical Schemes for Two-Fluid Plasma Equations , 2011, J. Sci. Comput..