Abstract
A parallel, fully coupled, nonlinearly implicit Newton-Krylov-Schwarz algorithm is proposed for the numerical simulation of a magnetic reconnection problem described by a system of resistive Hall magnetohydrodynamics equations in slab geometry. A key component of the algorithm is a restricted additive Schwarz preconditioner defined for problems with doubly periodic boundary conditions. We show numerically that with such a preconditioned nonlinearly implicit method the time step size is no longer constrained by the CFL number or the convergence of the Newton solver. We report the parallel performance of the algorithm and software on machines with thousands of processors.
Original language | English (US) |
---|---|
Pages (from-to) | 1919-1936 |
Number of pages | 18 |
Journal | Journal of Computational Physics |
Volume | 225 |
Issue number | 2 |
DOIs | |
State | Published - Aug 10 2007 |
Externally published | Yes |
Keywords
- Fully implicit
- Magnetohydrodynamics
- Newton-Krylov
- Parallel computing
- Restricted additive Schwarz preconditioner
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics