A fully divergence-free finite element method for magnetohydrodynamic equations
-
Hiptmair, Ralf
Seminar of Applied Mathematics, Swiss Federal Institute of Technology, Zurich, CH-8092, Switzerland
-
Li, Lingxiao
NCMIS, LSEC, ICMSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, School of Mathematical Science, Beijing 100190, P. R. China
-
Mao, Shipeng
NCMIS, LSEC, ICMSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, School of Mathematical Science, Beijing 100190, P. R. China
-
Zheng, Weiying
NCMIS, LSEC, ICMSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, School of Mathematical Science, Beijing 100190, P. R. China
Published in:
- Mathematical Models and Methods in Applied Sciences. - World Scientific Pub Co Pte Lt. - 2018, vol. 28, no. 04, p. 659-695
English
We propose a finite element method for the three-dimensional transient incompressible magnetohydrodynamic equations that ensures exactly divergence-free approximations of the velocity and the magnetic induction. We employ second-order semi-implicit timestepping, for which we rigorously establish an energy law and, as a consequence, unconditional stability. We prove unique solvability of the linear systems of equations to be solved in every timestep. For those we design an efficient preconditioner so that the number of preconditioned GMRES iterations is uniformly bounded with respect to the number of degrees of freedom. As both meshwidth and timestep size tend to zero, we prove that the discrete solutions converge to a weak solution of the continuous problem. Finally, by several numerical experiments, we confirm the predictions of the theory and demonstrate the efficiency of the preconditioner.
-
Language
-
-
Open access status
-
green
-
Identifiers
-
-
Persistent URL
-
https://sonar.ch/global/documents/229370
Statistics
Document views: 61
File downloads: