Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion
Journal Article
·
· SIAM Journal on Scientific Computing
- RIKEN Center for Computational Science, Kobe (Japan)
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States)
- Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States); University of New Mexico, Albuquerque, NM (United States)
- Sandia National Laboratories (SNL-CA), Livermore, CA (United States)
This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.
- Research Organization:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR). Scientific Discovery through Advanced Computing (SciDAC); USDOE Office of Science (SC), Fusion Energy Sciences (FES)
- Grant/Contract Number:
- 89233218CNA000001
- OSTI ID:
- 2478091
- Report Number(s):
- LA-UR--23-33406
- Journal Information:
- SIAM Journal on Scientific Computing, Journal Name: SIAM Journal on Scientific Computing Journal Issue: 5 Vol. 46; ISSN 1064-8275
- Publisher:
- Society for Industrial and Applied Mathematics (SIAM)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF
Initial performance of fully-coupled AMG and approximate block factorization preconditioners for solution of implicit FE resistive MHD.
Scalable implicit incompressible resistive MHD with stabilized FE and fully-coupled Newton–Krylov-AMG
Journal Article
·
Tue Aug 22 20:00:00 EDT 2023
· Computer Methods in Applied Mechanics and Engineering
·
OSTI ID:1998134
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