Analysis of a mixed semi-implicit/implicit algorithm for low-frequency two-fluid plasma modeling

A temporally staggered algorithm for advancing solutions of the two-fluid plasma model is analyzed with von Neumann's method and with differential approximation. The implicit leapfrog algorithm [C.R. Sovinec et al., J. Phys. Conf. Series 16 (2005) 25-34] is found to be numerically stable at arbitrarily large time-step when the advective, Hall, and gyroviscous terms are temporally centered in their respective advances and the coefficient of the semi-implicit operator meets the criterion found for basic hyperbolic systems. Numerical instability with forward or backward differencing of advection is evident as an ill-posed equation in the differential approximation for a simplified system. At large time-step, the accuracy of the algorithm is comparable to that of the Crank-Nicolson method for all plane waves except the parallel mode that is sensitive to the ion cyclotron resonance. An implementation reproduces theoretical results on the transition from resistive magnetohydrodynamics to two-fluid reconnection in a sheared-slab linear tearing mode. A nonlinear three-dimensional computation in toroidal geometry shows an increasing exponentiation rate of kinetic energy as magnetic reconnection from an internal kink mode changes from current-sheet to 'X-point' geometry.

[1]  Lawrence F. Shampine.  Stability of the leapfrog/midpoint method , 2009, Appl. Math. Comput..

[2]  Luis Chacon,et al.  An optimal, parallel, fully implicit Newton–Krylov solver for three-dimensional viscoresistive magnetohydrodynamicsa) , 2008 .

[3]  James Drake,et al.  Two-fluid theory of collisionless magnetic reconnection , 1997 .

[4]  R. Denton,et al.  Role of dispersive waves in collisionless magnetic reconnection. , 2001, Physical review letters.

[5]  L. Margolin,et al.  On balanced approximations for time integration of multiple time scale systems , 2003 .

[6]  J. Linker,et al.  Stability of Algorithms for Waves with Large Flows , 1999 .

[7]  D. A. Knoll,et al.  A 2D high-ß Hall MHD implicit nonlinear solver , 2003 .

[8]  Sin-Chung Chang,et al.  A critical analysis of the modified equation technique of Warming and Hyett , 1990 .

[9]  Scott Kruger,et al.  Nonlinear extended magnetohydrodynamics simulation using high-order finite elements , 2005 .

[10]  Y. Shokin.  Analysis of conservative properties of the difference schemes by the method of differential approximation , 1981 .

[11]  R. F. Warming,et al.  The modified equation approach to the stability and accuracy analysis of finite-difference methods , 1974 .

[12]  Olivier Czarny,et al.  Bézier surfaces and finite elements for MHD simulations , 2008, J. Comput. Phys..

[13]  C. Sovinec,et al.  Global axisymmetric simulations of two-fluid reconnection in an experimentally relevant geometry , 2008 .

[14]  M. Shay,et al.  The Hall fields and fast magnetic reconnection , 2008 .

[15]  Stephen C. Jardin,et al.  A high-order implicit finite element method for integrating the two-fluid magnetohydrodynamic equations in two dimensions , 2007, J. Comput. Phys..

[16]  D. Schnack,et al.  Absence of complete finite-Larmor-radius stabilization in extended MHD. , 2008, Physical review letters.

[17]  Dalton D. Schnack,et al.  Semi-implicit magnetohydrodynamic calculations , 1987 .

[18]  Scott Kruger,et al.  Computational modeling of fully ionized magnetized plasmas using the fluid approximation , 2005 .

[19]  K Lerbinger,et al.  A new semi-implicit method for MHD computations , 1991 .

[20]  C. W. Hirt.  Heuristic stability theory for finite-difference equations☆ , 1968 .

[21]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[22]  P. Catto,et al.  Drift-ordered fluid equations for field-aligned modes in low-β collisional plasma with equilibrium pressure pedestals , 2003 .

[23]  Carol S. Woodward,et al.  A fully implicit numerical method for single-fluid resistive magnetohydrodynamics , 2006, J. Comput. Phys..

[24]  Douglas S. Harned,et al.  Accurate semi-implicit treatment of the hall effect in magnetohydrodynamic computations , 1989 .

[25]  Charlson C. Kim,et al.  Hall magnetohydrodynamics simulations of end-shorting induced rotation in field- reversed configurations , 2007 .

[26]  James Demmel,et al.  SuperLU_DIST: A scalable distributed-memory sparse direct solver for unsymmetric linear systems , 2003, TOMS.

[27]  P. Fischer,et al.  High-Order Methods for Incompressible Fluid Flow , 2002 .

[28]  Boris B. Kadomtsev,et al.  Reviews of Plasma Physics , 2012 .

[29]  Francesco Pegoraro,et al.  Magnetic reconnection in electron magnetohydrodynamics , 1992 .

[30]  Edward J Caramana.  Derivation of implicit difference schemes by the method of differential approximation , 1991 .

[31]  Steven J. Plimpton,et al.  Nonlinear magnetohydrodynamics simulation using high-order finite elements , 2004 .

[32]  J. M. Sanz-Serna,et al.  Studies in Numerical Nonlinear Instability I. Why do Leapfrog Schemes Go Unstable , 1985 .

[33]  C. Hegna,et al.  Two-fluid tearing instability in force-free magnetic configuration , 2004 .

[34]  Stephen C. Jardin,et al.  Calculations of two-fluid magnetohydrodynamic axisymmetric steady-states , 2009, J. Comput. Phys..

[35]  S. I. Braginskii.  Transport Processes in a Plasma , 1965 .

[36]  Nuno F. Loureiro,et al.  An iterative semi-implicit scheme with robust damping , 2007, J. Comput. Phys..

[37]  Hinrich Lütjens,et al.  The XTOR code for nonlinear 3D simulations of MHD instabilities in tokamak plasmas , 2008, J. Comput. Phys..

[38]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[39]  Parametric Analysis of the Two-Fluid Tearing Instability , 2009 .

[40]  A. Y. Aydemir,et al.  Nonlinear studies of m = 1 modes in high-temperature plasmas , 1992 .

[41]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .