Krylov Smoothing for Fully-Coupled AMG Preconditioners for VMS Resistive MHD

Lin, Paul; Shadid, John N.; Tsuji, Paul H.

doi:10.1007/978-3-030-30705-9_24

Cited by 5 publications

(1 citation statement)

References 20 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The next example we consider is the resistive tearing mode [3,25,50,51]. The tearing mode is one of the most fundamental resistive instabilities.…”

Section: Tearing Modementioning

confidence: 99%

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

Tang,

Chacon,

Kolev

et al. 2021

Preprint

Self Cite

View full text Add to dashboard Cite

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 (mfem.org). 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 [Chacón, 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 nonconforming 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 (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

show abstract

“…The next example we consider is the resistive tearing mode [3,25,50,51]. The tearing mode is one of the most fundamental resistive instabilities.…”

Section: Tearing Modementioning

confidence: 99%