Entropy Stable Numerical Schemes for Two-Fluid Plasma Equations

Two-fluid ideal plasma equations are a generalized form of the ideal MHD equations in which electrons and ions are considered as separate species. The design of efficient numerical schemes for the these equations is complicated on account of their non-linear nature and the presence of stiff source terms, especially for high charge to mass ratios and for low Larmor radii. In this article, we design entropy stable finite difference schemes for the two-fluid equations by combining entropy conservative fluxes and suitable numerical diffusion operators. Furthermore, to overcome the time step restrictions imposed by the stiff source terms, we devise time-stepping routines based on implicit-explicit (IMEX)-Runge Kutta (RK) schemes. The special structure of the two-fluid plasma equations is exploited by us to design IMEX schemes in which only local (in each cell) linear equations need to be solved at each time step. Benchmark numerical experiments are presented to illustrate the robustness and accuracy of these schemes.

[1]  Philip L. Roe,et al.  Affordable, entropy-consistent Euler flux functions II: Entropy production at shocks , 2009, J. Comput. Phys..

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

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

[4]  E. Tadmor.  Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems , 2003, Acta Numerica.

[5]  Timothy J. Barth,et al.  Numerical Methods for Gasdynamic Systems on Unstructured Meshes , 1997, Theory and Numerics for Conservation Laws.

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

[7]  Xiangxiong Zhang,et al.  On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes , 2010, J. Comput. Phys..

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

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

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

[11]  Eitan Tadmor,et al.  Arbitrarily High-order Accurate Entropy Stable Essentially Nonoscillatory Schemes for Systems of Conservation Laws , 2012, SIAM J. Numer. Anal..

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

[13]  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..

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

[15]  Eitan Tadmor,et al.  The numerical viscosity of entropy stable schemes for systems of conservation laws. I , 1987 .

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

[17]  C. M. Dafermos,et al.  Hyberbolic [i.e. Hyperbolic] conservation laws in continuum physics , 2005 .

[18]  C. Dafermos.  Hyberbolic Conservation Laws in Continuum Physics , 2000 .

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