A high-order implicit finite element method for integrating the two-fluid magnetohydrodynamic equations in two dimensions

We describe a new method for solving the time-dependent two-fluid magnetohydrodynamic (2F-MHD) equations in two dimensions that has significant advantages over other methods. The stream-function/potential representation of the velocity and magnetic field vectors, while fully general, allows accurate description of nearly incompressible fluid motions and manifestly satisfies the divergence condition on the magnetic field. Through analytic manipulation, the split semi-implicit method breaks the full matrix time advance into four sequential time advances, each involving smaller matrices. The use of a high-order triangular element with continuous first derivatives (C1 continuity) allows the Galerkin method to be applied without introduction of new auxiliary variables (such as the vorticity or the current density). These features, along with the manifestly compact nature of the fully node-based C1 finite elements, lead to minimum size matrices for an unconditionally stable method with order of accuracy h4. The resulting matrices are compatible with direct factorization using SuperLU_dist. We demonstrate the accuracy of the method by presenting examples of two-fluid linear wave propagation, two-fluid linear eigenmodes of a tilting cylinder, and of a challenging nonlinear problem in two-fluid magnetic reconnection.

[1]  A. H. Glasser,et al.  The SEL macroscopic modeling code , 2004, Comput. Phys. Commun..

[2]  J. Z. Zhu,et al.  The finite element method , 1977 .

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

[4]  Bruce D. Scott,et al.  Free-energy conservation in local gyrofluid models , 2005 .

[5]  M. Zlámal,et al.  SOME RECENT ADVANCES IN THE MATHEMATICS OF FINITE ELEMENTS , 1973 .

[6]  Michael Hesse,et al.  Geospace Environmental Modeling (GEM) magnetic reconnection challenge , 2001 .

[7]  G. Strang,et al.  An Analysis of the Finite Element Method , 1974 .

[8]  Richard D. Sydora,et al.  Magnetic reconnection driven by current repulsion , 1990 .

[9]  J. Breslau,et al.  Global extended magnetohydrodynamic studies of fast magnetic reconnection , 2003 .

[10]  I. N. Molchanov,et al.  On an Approach to Integrating Boundary Problems With a Non-Unique Solution , 1972, Inf. Process. Lett..

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

[12]  W. Park,et al.  Plasma simulation studies using multilevel physics models , 1999 .

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

[14]  Dana A. Knoll,et al.  An Implicit, Nonlinear Reduced Resistive MHD Solver , 2002 .

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

[16]  Stephen C. Jardin,et al.  A triangular finite element with first-derivative continuity applied to fusion MHD applications , 2004 .

[17]  N. M. Ferraro,et al.  Finite element implementation of Braginskii's gyroviscous stress with application to the gravitational instability , 2006 .

[18]  Stephen C. Jardin,et al.  Implicit solution of the four-field extended-magnetohydrodynamic equations using high-order high-continuity finite elementsa) , 2004 .

[19]  J. Whiteman.  The Mathematics of Finite Elements and Applications. , 1983 .