Implicit solution of the four-field extended-magnetohydrodynamic equations using high-order high-continuity finite elements

Here we describe a technique for solving the four-field extended-magnetohydrodynamic (MHD) equations in two dimensions. The introduction of triangular high-order finite elements with continuous first derivatives (C1 continuity) leads to a compact representation compatible with direct inversion of the associated sparse matrices. The split semi-implicit method is introduced and used to integrate the equations in time, yielding unconditional stability for arbitrary time step. The method is applied to the cylindrical tilt mode problem with the result that a nonzero value of the collisionless ion skin depth will increase the growth rate of that mode. The effect of this parameter on the reconnection rate and geometry of a Harris equilibrium and on the Taylor reconnection problem is also demonstrated. This method forms the basis for a generalization to a full extended-MHD description of the plasma with six, eight, or more scalar fields.