A decoupled, unconditionally energy-stable and structure-preserving finite element scheme for the incompressible MHD equations with magnetic-electric formulation
被引:3
作者:
Zhang, Xiaodi
论文数: 0引用数: 0
h-index: 0
机构:
Zhengzhou Univ, Henan Acad Big Data, Zhengzhou 450052, Peoples R China
Zhengzhou Univ, Sch Math & Stat, Zhengzhou 450001, Peoples R ChinaZhengzhou Univ, Henan Acad Big Data, Zhengzhou 450052, Peoples R China
Zhang, Xiaodi
[1
,2
]
Su, Haiyan
论文数: 0引用数: 0
h-index: 0
机构:
Xinjiang Univ, Coll Math & Syst Sci, Urumqi 830046, Peoples R ChinaZhengzhou Univ, Henan Acad Big Data, Zhengzhou 450052, Peoples R China
Su, Haiyan
[3
]
机构:
[1] Zhengzhou Univ, Henan Acad Big Data, Zhengzhou 450052, Peoples R China
[2] Zhengzhou Univ, Sch Math & Stat, Zhengzhou 450001, Peoples R China
[3] Xinjiang Univ, Coll Math & Syst Sci, Urumqi 830046, Peoples R China
MHD equations;
Decoupled scheme;
Finite element method;
Energy stable;
Structure-preserving;
DIFFUSE INTERFACE MODEL;
NUMERICAL-SOLUTION;
ERROR ANALYSIS;
APPROXIMATION;
STATIONARY;
D O I:
10.1016/j.camwa.2023.06.029
中图分类号:
O29 [应用数学];
学科分类号:
070104 ;
摘要:
In this paper, we propose a decoupled, unconditionally energy stable and structure-preserving finite element scheme for the incompressible magnetohydrodynamic (MHD) equations with magnetic-electric formulation. The spatial discretization is based on mixed finite element method and finite element exterior calculus, where the hydrodynamic unknowns are approximated by stable finite element pair, the magnetic field and electric field are discretized by H(div)-conforming face element and H(curl)-conforming edge element. The time marching is combining the first-order semi-implicit Euler scheme, a first-order stabilized term and some subtle implicit -explicit treatments for the coupling terms. As a result, the fully discrete scheme only needs to solve some linear sub-problems at each time step and yields an exactly divergence-free magnetic field on the discrete level. In particular, we decouple the computations of the magnetic field and electric field equivalently by utilizing the structure-preserving property. Furthermore, the unique solvability and unconditional stability of the scheme are proved. Numerical experiments are performed to verify the efficiency and accuracy of the proposed scheme.