Abstract:The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (<a class="link-external link-http" href="http://mfem.org" rel="external noopener nofollow">this http URL</a>). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in [Chacon, et al, JCP 2002] to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime ($\ge 10^7$) with the successful resolution of plasmoid instabilities and thin current sheets.

Adaptive mesh strategies for the spectral element method

Mesh Smoothing for the Spectral Element Method

A New High-Order Spectral Difference Method for Simulating Viscous Flows on Unstructured Grids with Mixed-Element Meshes

A spectral element method for meshes with skinny elements

Adaptive Hermite Spectral Methods in Unbounded Domains

Adaptive residual refinement in an RBF finite difference scheme for 2D time-dependent problems

A multi-mesh finite element method for Lagrange elements of arbitrary degree

Chebyshev Finite Spectral Method with Extended Moving Grids

High-resolution solutions of nonlinear, advection-dominated problems using a generalized finite element method

Spectral Analysis of High Order Continuous FEM for Hyperbolic PDEs on Triangular Meshes: Influence of Approximation, Stabilization, and Time-Stepping

Anisotropic space-time adaptation for reaction-diffusion problems

Mesh update techniques for free-surface flow solvers using spectral element method

A spectral element method for fluid dynamics: Laminar flow in a channel expansion

A posteriori error analysis and mesh adaptivity for a virtual element method solving the Stokes equations

Efficient adaptivity for simulating cardiac electrophysiology with spectral deferred correction methods

Spectral Element Simulations of Interactive Particles in a Fluid

Finite element analysis of a spectral problem on curved meshes occurring in diffusion with high order boundary conditions

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

A Three-Dimensional Spectral Boundary Element Algorithm for Interfacial Dynamics in Stokes Flow