Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Towards a Scalable Fully-Implicit Fully-coupled Resistive MHD Formulation with Stabilized FE Methods

Journal Article · · Journal of Computational Physics
This paper explores the development of a scalable, nonlinear, fully-implicit stabilized unstructured finite element (FE) capability for 2D incompressible (reduced) resistive MHD. The discussion considers the implementation of a stabilized FE formulation in context of a fully-implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton-Krylov methods, which are preconditioned using fully-coupled algebraic multilevel preconditioners. These preconditioners are shown to enable a robust, scalable and efficient solution approach for the large-scale sparse linear systems generated by the Newton linearization. Verification results demonstrate the expected order-of-accuracy for the stabilized FE discretization. The approach is tested on a variety of prototype problems, including both low-Lundquist number (e.g., an MHD Faraday conduction pump and a hydromagnetic Rayleigh-Bernard linear stability calculation) and moderately-high Lundquist number (magnetic island coalescence problem) examples. Initial results that explore the scaling of the solution methods are presented on up to 4096 processors for problems with up to 64M unknowns on a CrayXT3/4. Additionally, a large-scale proof-of-capability calculation for 1 billion unknowns for the MHD Faraday pump problem on 24,000 cores is presented.
Research Organization:
Oak Ridge National Laboratory (ORNL)
Sponsoring Organization:
SC USDOE - Office of Science (SC)
DOE Contract Number:
AC05-00OR22725
OSTI ID:
985756
Journal Information:
Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: 229 Vol. 229; ISSN 1090-2716; ISSN 0021-9991
Country of Publication:
United States
Language:
English

Similar Records

Towards a Scalable Fully-Implicit Fully-coupled Resistive MHD Formulation with Stabilized FE Methods
Journal Article · Wed Jun 03 00:00:00 EDT 2009 · Journal of Computational Physics · OSTI ID:985484

Initial performance of fully-coupled AMG and approximate block factorization preconditioners for solution of implicit FE resistive MHD.
Conference · Tue Jun 01 00:00:00 EDT 2010 · OSTI ID:1020538

Scalable implicit incompressible resistive MHD with stabilized FE and fully-coupled Newton–Krylov-AMG
Journal Article · Tue Feb 09 19:00:00 EST 2016 · Computer Methods in Applied Mechanics and Engineering · OSTI ID:1310309