A Diffuse Interface Model and Semi-Implicit Energy Stable Finite Element Method for Two-Phase Magnetohydrodynamic Flows
In this paper, we propose a diffuse interface model and finite element approximation for two-phase magnetohydrodynamic (MHD) flows with different viscosities and electric conductivities. An energy stable scheme, which is based on the finite element method for the spatial discretization and first order semi-implicit scheme combined with convex splitting method for the temporal discretization, is proposed to solve this new model. The numerical scheme is proved to be mass-conservative and energy law preserving. By Leray-Schauder fixed point theorem, the existence of solutions to the numerical scheme is shown. The uniqueness of the numerical solutions is obtained. Utilizing the stability of the numerical scheme and the compactness method, the existence of the weak solutions to the two-phase MHD model is established as well. Furthermore, given more regularity on the weak solution, the convergence of the numerical scheme is derived. Finally, numerical experiments are provided to verify the theoretical results and validate the proposed model.
J. Yang et al., "A Diffuse Interface Model and Semi-Implicit Energy Stable Finite Element Method for Two-Phase Magnetohydrodynamic Flows," Computer Methods in Applied Mechanics and Engineering, vol. 356, pp. 435-464, Elsevier B.V., Nov 2019.
The definitive version is available at https://doi.org/10.1016/j.cma.2019.07.022
Mathematics and Statistics
Center for High Performance Computing Research
Keywords and Phrases
Diffuse interface models; Energy stability; Finite element method; Semi-implicit scheme; Two-phase MHD flows
International Standard Serial Number (ISSN)
Article - Journal
© 2019 Elsevier B.V., All rights reserved.
01 Nov 2019