DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Monolithic Multigrid Methods for Magnetohydrodynamics

Abstract

The magnetohydrodynamics equations model a wide range of plasma physics applications and are characterized by a nonlinear system of partial differential equations that strongly couples a charged fluid with the evolution of electromagnetic fields. After discretization and linearization, the resulting system of equations is generally difficult to solve due to the coupling between variables and the heterogeneous coefficients induced by the linearization process. In this paper, we investigate multigrid preconditioners for this system based on specialized relaxation schemes that properly address the system structure and coupling. Here, three extensions of Vanka relaxation are proposed and applied to problems with up to 170 million degrees of freedom and fluid and magnetic Reynolds numbers up to 400 for stationary problems and up to 20,000 for time-dependent problems.

Authors:
ORCiD logo [1];  [2];  [3]; ORCiD logo [4]; ORCiD logo [5];  [3]
  1. Tufts Univ., Medford, MA (United States)
  2. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  3. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
  4. Univ. of Oxford (United Kingdom)
  5. Univ. of Newfoundland (Canada)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1841848
Report Number(s):
LLNL-JRNL-811901
Journal ID: ISSN 1064-8275; 1018684
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
SIAM Journal on Scientific Computing
Additional Journal Information:
Journal Volume: 43; Journal Issue: 5; Journal ID: ISSN 1064-8275
Publisher:
Society for Industrial and Applied Mathematics (SIAM)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; magnetohydrodynamics; monolithic multigrid; Vanka relaxation

Citation Formats

Adler, James H., Benson, Thomas R., Cyr, Eric C., Farrell, Patrick E., MacLachlan, Scott P., and Tuminaro, Ray S. Monolithic Multigrid Methods for Magnetohydrodynamics. United States: N. p., 2021. Web. doi:10.1137/20m1348364.
Adler, James H., Benson, Thomas R., Cyr, Eric C., Farrell, Patrick E., MacLachlan, Scott P., & Tuminaro, Ray S. Monolithic Multigrid Methods for Magnetohydrodynamics. United States. https://doi.org/10.1137/20m1348364
Adler, James H., Benson, Thomas R., Cyr, Eric C., Farrell, Patrick E., MacLachlan, Scott P., and Tuminaro, Ray S. Tue . "Monolithic Multigrid Methods for Magnetohydrodynamics". United States. https://doi.org/10.1137/20m1348364. https://www.osti.gov/servlets/purl/1841848.
@article{osti_1841848,
title = {Monolithic Multigrid Methods for Magnetohydrodynamics},
author = {Adler, James H. and Benson, Thomas R. and Cyr, Eric C. and Farrell, Patrick E. and MacLachlan, Scott P. and Tuminaro, Ray S.},
abstractNote = {The magnetohydrodynamics equations model a wide range of plasma physics applications and are characterized by a nonlinear system of partial differential equations that strongly couples a charged fluid with the evolution of electromagnetic fields. After discretization and linearization, the resulting system of equations is generally difficult to solve due to the coupling between variables and the heterogeneous coefficients induced by the linearization process. In this paper, we investigate multigrid preconditioners for this system based on specialized relaxation schemes that properly address the system structure and coupling. Here, three extensions of Vanka relaxation are proposed and applied to problems with up to 170 million degrees of freedom and fluid and magnetic Reynolds numbers up to 400 for stationary problems and up to 20,000 for time-dependent problems.},
doi = {10.1137/20m1348364},
journal = {SIAM Journal on Scientific Computing},
number = 5,
volume = 43,
place = {United States},
year = {Tue Feb 16 00:00:00 EST 2021},
month = {Tue Feb 16 00:00:00 EST 2021}
}

Works referenced in this record:

Robust preconditioners for incompressible MHD models
journal, July 2016


Coalescence of magnetic islands, sloshing, and the pressure problem
journal, March 2006

  • Knoll, D. A.; Chacón, L.
  • Physics of Plasmas, Vol. 13, Issue 3
  • DOI: 10.1063/1.2173515

Preconditioners for Mixed Finite Element Discretizations of Incompressible MHD Equations
journal, January 2017

  • Wathen, Michael; Greif, Chen; Schötzau, Dominik
  • SIAM Journal on Scientific Computing, Vol. 39, Issue 6
  • DOI: 10.1137/16M1098991

A New Approximate Block Factorization Preconditioner for Two-Dimensional Incompressible (Reduced) Resistive MHD
journal, January 2013

  • Cyr, Eric C.; Shadid, John N.; Tuminaro, Raymond S.
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 3
  • DOI: 10.1137/12088879X

Multigrid in H (div) and H (curl)
journal, April 2000

  • Arnold, Douglas N.; Falk, Richard S.; Winther, Ragnar
  • Numerische Mathematik, Vol. 85, Issue 2
  • DOI: 10.1007/PL00005386

Scalable implicit incompressible resistive MHD with stabilized FE and fully-coupled Newton–Krylov-AMG
journal, June 2016

  • Shadid, J. N.; Pawlowski, R. P.; Cyr, E. C.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 304
  • DOI: 10.1016/j.cma.2016.01.019

Numerical Analysis of Vanka‐Type Solvers for Steady Stokes and Navier–Stokes Flows
journal, January 2006

  • Manservisi, S.
  • SIAM Journal on Numerical Analysis, Vol. 44, Issue 5
  • DOI: 10.1137/060655407

On a local Fourier analysis for overlapping block smoothers on triangular grids
journal, July 2016


IMEX and exact sequence discretization of the multi-fluid plasma model
journal, November 2019


GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986

  • Saad, Youcef; Schultz, Martin H.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
  • DOI: 10.1137/0907058

Mixed finite element methods for stationary incompressible magneto?hydrodynamics
journal, February 2004


A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier–Stokes equations
journal, January 2008

  • Elman, Howard; Howle, V. E.; Shadid, John
  • Journal of Computational Physics, Vol. 227, Issue 3
  • DOI: 10.1016/j.jcp.2007.09.026

A Block Preconditioner for an Exact Penalty Formulation for Stationary MHD
journal, January 2014

  • Phillips, Edward G.; Elman, Howard C.; Cyr, Eric C.
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 6
  • DOI: 10.1137/140955082

Firedrake: Automating the Finite Element Method by Composing Abstractions
journal, December 2016

  • Rathgeber, Florian; Ham, David A.; Mitchell, Lawrence
  • ACM Transactions on Mathematical Software, Vol. 43, Issue 3
  • DOI: 10.1145/2998441

On Schwarz-type Smoothers for Saddle Point Problems
journal, August 2003


Choosing the Forcing Terms in an Inexact Newton Method
journal, January 1996

  • Eisenstat, Stanley C.; Walker, Homer F.
  • SIAM Journal on Scientific Computing, Vol. 17, Issue 1
  • DOI: 10.1137/0917003

Numerical performance of smoothers in coupled multigrid methods for the parallel solution of the incompressible Navier-Stokes equations
journal, January 2000


Efficient Mesh Management in Firedrake Using PETSc DMPlex
journal, January 2016

  • Lange, Michael; Mitchell, Lawrence; Knepley, Matthew G.
  • SIAM Journal on Scientific Computing, Vol. 38, Issue 5
  • DOI: 10.1137/15M1026092

Multigrid relaxation methods for systems of saddle point type
journal, December 2008


Vector-potential finite-element formulations for two-dimensional resistive magnetohydrodynamics
journal, January 2019

  • Adler, James H.; He, Yunhui; Hu, Xiaozhe
  • Computers & Mathematics with Applications, Vol. 77, Issue 2
  • DOI: 10.1016/j.camwa.2018.09.051

A Scalable Approximate Inverse Block Preconditioner for an Incompressible Magnetohydrodynamics Model Problem
journal, January 2020

  • Wathen, Michael; Greif, Chen
  • SIAM Journal on Scientific Computing, Vol. 42, Issue 1
  • DOI: 10.1137/19M1255409

Convergence of a B - E based finite element method for MHD models on Lipschitz domains
journal, April 2020


Monolithic Multigrid Methods for Two-Dimensional Resistive Magnetohydrodynamics
journal, January 2016

  • Adler, James H.; Benson, Thomas R.; Cyr, Eric C.
  • SIAM Journal on Scientific Computing, Vol. 38, Issue 1
  • DOI: 10.1137/151006135

Mixed finite elements in ?3
journal, September 1980


Island Coalescence Using Parallel First-Order System Least Squares on Incompressible Resistive Magnetohydrodynamics
journal, January 2013

  • Adler, J. H.; Brezina, M.; Manteuffel, T. A.
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 5
  • DOI: 10.1137/120880227

A mixed finite element method with exactly divergence-free velocities for incompressible magnetohydrodynamics
journal, November 2010

  • Greif, Chen; Li, Dan; Schötzau, Dominik
  • Computer Methods in Applied Mechanics and Engineering, Vol. 199, Issue 45-48
  • DOI: 10.1016/j.cma.2010.05.007

A combined conjugate gradient - multi-grid algorithm for the numerical solution of the Stokes problem
journal, October 1984


Solver Composition Across the PDE/Linear Algebra Barrier
journal, January 2018

  • Kirby, Robert C.; Mitchell, Lawrence
  • SIAM Journal on Scientific Computing, Vol. 40, Issue 1
  • DOI: 10.1137/17M1133208

Nonlinear, three-dimensional magnetohydrodynamics of noncircular tokamaks
journal, January 1976


Local Fourier analysis for multigrid with overlapping smoothers applied to systems of PDEs
journal, January 2011

  • MacLachlan, Scott P.; Oosterlee, C. W.
  • Numerical Linear Algebra with Applications, Vol. 18, Issue 4
  • DOI: 10.1002/nla.762

Implicit adaptive mesh refinement for 2D reduced resistive magnetohydrodynamics
journal, October 2008

  • Philip, Bobby; Chacón, Luis; Pernice, Michael
  • Journal of Computational Physics, Vol. 227, Issue 20
  • DOI: 10.1016/j.jcp.2008.06.029

A general nonlinear fluid model for reacting plasma-neutral mixtures
journal, July 2012

  • Meier, E. T.; Shumlak, U.
  • Physics of Plasmas, Vol. 19, Issue 7
  • DOI: 10.1063/1.4736975

A comparative study of efficient iterative solvers for generalized Stokes equations
journal, January 2008

  • Larin, Maxim; Reusken, Arnold
  • Numerical Linear Algebra with Applications, Vol. 15, Issue 1
  • DOI: 10.1002/nla.561

A Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations
journal, February 2011


Block Preconditioners for Stable Mixed Nodal and Edge finite element Representations of Incompressible Resistive MHD
journal, January 2016

  • Phillips, Edward G.; Shadid, John N.; Cyr, Eric C.
  • SIAM Journal on Scientific Computing, Vol. 38, Issue 6
  • DOI: 10.1137/16M1074084

An Augmented Lagrangian Preconditioner for the 3D Stationary Incompressible Navier--Stokes Equations at High Reynolds Number
journal, January 2019

  • Farrell, Patrick E.; Mitchell, Lawrence; Wechsung, Florian
  • SIAM Journal on Scientific Computing, Vol. 41, Issue 5
  • DOI: 10.1137/18M1219370

Numerical solution of saddle point problems
journal, April 2005


Block-implicit multigrid solution of Navier-Stokes equations in primitive variables
journal, July 1986


A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling
journal, January 2001

  • Amestoy, Patrick R.; Duff, Iain S.; L'Excellent, Jean-Yves
  • SIAM Journal on Matrix Analysis and Applications, Vol. 23, Issue 1
  • DOI: 10.1137/S0895479899358194

Preconditioners for the discretized time-harmonic Maxwell equations in mixed form
journal, January 2007

  • Greif, Chen; Schötzau, Dominik
  • Numerical Linear Algebra with Applications, Vol. 14, Issue 4
  • DOI: 10.1002/nla.515

Mixed finite elements for incompressible magneto-hydrodynamics
journal, July 2003


Nested Iteration and First-Order System Least Squares for Incompressible, Resistive Magnetohydrodynamics
journal, January 2010

  • Adler, J. H.; Manteuffel, T. A.; McCormick, S. F.
  • SIAM Journal on Scientific Computing, Vol. 32, Issue 3
  • DOI: 10.1137/090766905

Block-implicit multigrid calculation of two-dimensional recirculating flows
journal, November 1986


Higher-order finite element discretizations in a benchmark problem for incompressible flows
journal, January 2001

  • John, Volker; Matthies, Gunar
  • International Journal for Numerical Methods in Fluids, Vol. 37, Issue 8
  • DOI: 10.1002/fld.195

A robust solver for the finite element approximation of stationary incompressible MHD equations in 3D
journal, December 2017


A Simple and Efficient Segregated Smoother for the Discrete Stokes Equations
journal, January 2014

  • Gaspar, Francisco J.; Notay, Yvan; Oosterlee, Cornelis W.
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 3
  • DOI: 10.1137/130920630

On iterative methods for the incompressible Stokes problem
journal, March 2011

  • ur Rehman, M.; Geenen, T.; Vuik, C.
  • International Journal for Numerical Methods in Fluids, Vol. 65, Issue 10
  • DOI: 10.1002/fld.2235

Towards a scalable fully-implicit fully-coupled resistive MHD formulation with stabilized FE methods
journal, October 2010

  • Shadid, J. N.; Pawlowski, R. P.; Banks, J. W.
  • Journal of Computational Physics, Vol. 229, Issue 20
  • DOI: 10.1016/j.jcp.2010.06.018

Toward textbook multigrid efficiency for fully implicit resistive magnetohydrodynamics
journal, September 2010

  • Adams, Mark F.; Samtaney, Ravi; Brandt, Achi
  • Journal of Computational Physics, Vol. 229, Issue 18
  • DOI: 10.1016/j.jcp.2010.04.024

An Implicit, Nonlinear Reduced Resistive MHD Solver
journal, May 2002

  • Chacón, L.; Knoll, D. A.; Finn, J. M.
  • Journal of Computational Physics, Vol. 178, Issue 1
  • DOI: 10.1006/jcph.2002.7015

A Note on Preconditioning for Indefinite Linear Systems
journal, January 2000

  • Murphy, Malcolm F.; Golub, Gene H.; Wathen, Andrew J.
  • SIAM Journal on Scientific Computing, Vol. 21, Issue 6
  • DOI: 10.1137/S1064827599355153