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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:351

Article

Li, Lingxiao^1,2  Zheng, Weiying^1,2 

[1] Chinese Acad Sci, Acad Math & Syst Sci, Inst Computat Math & Sci Engn Comp, NCMIS,LSEC, Beijing 100190, Peoples R China

[2] Univ Chinese Acad Sci, Sch Math Sci, Beijing 100049, Peoples R China

关键词: Incompressible magnetohydrodynamic equations;    Mixed finite element method;    Block preconditioner;    grad-div stabilization;    Newton's method;   

DOI  :  10.1016/j.jcp.2017.09.025

来源: Elsevier

【 摘 要 】

In this paper, we propose a Newton-Krylov solver and a Picard-Krylov solver for finite element discrete problem of stationary incompressible magnetohydrodynamic equations in three dimensions. Using a mixed finite element method, we discretize the velocity and the pressure by H-1(Omega)-conforming finite elements and discretize the magnetic field by H(curl, Omega)-conforming edge elements. An efficient preconditioner is proposed to accelerate the convergence of GMRES method for solving linearized discrete problems. By extensive numerical experiments, we demonstrate the robustness of the Newton-Krylov solver for relatively large physical parameters and the optimality with respect to the number of degrees of freedom. Moreover, the numerical experiments show that the Newton-Krylov solver is more robust than the Picard-Krylov solver for large Reynolds number. (C) 2017 Elsevier Inc. All rights reserved.

【 授权许可】

Free