A constrained transport divergence-free finite element method for incompressible MHD equations
发布时间:2021年10月11日
浏览次数:5113
发布者: Yu Feng
主讲人: 张东航
活动时间: 从 2021-10-15 10:00 到 11:30
场地: 北京国际数学研究中心,全斋全29教室
In this talk, we study finite element method for three-dimensional incompressible resistive magnetohydrodynamic equations, in which the velocity, the current density, and the magnetic induction are divergence-free. It is desirable that the discrete solutions should also satisfy divergence-free conditions exactly especially for the momentum equations. Inspired by constrained transport method, we devise a new stable mixed finite element method that can achieve the goal. We also prove the well-posedness of the discrete solutions. To solve the resulting linear algebraic equations, we propose a GMRES solver with an augmented Lagrangian block preconditioner. By numerical experiments, we verify the theoretical results and demonstrate the quasi-optimality of the discrete solver with respect to the number of degrees of freedom. A comparison with other discretization using lid driven cavity is also given.