Highly Efficient and Energy Stable Schemes for the 2D/3D Diffuse Interface Model of Two-Phase Magnetohydrodynamics
Su Haiyan,Zhang Guo-Dong
DOI: https://doi.org/10.1007/s10915-021-01741-3
2021-01-01
Journal of Scientific Computing
Abstract:In this paper, we consider a two-phase MHD (magnetohydrodynamics) model based on the phase-field method. The governing model consists of the phase field equations, the Navier–Stokes equations and the Maxwell equations via the convection term, surface stress force, the Lorentz force and Ohm’s law. The model couples the phase field function, velocity, magnetic field, pressure and the introduced Lagrange multiplier, which results in a strong coupled, nonlinear and double saddle points type system. To solve the model efficiently, we propose two totally decoupled, linear and unconditionally energy stable schemes. The given schemes combine semi-implicit stabilization/IEQ methods for the phase-field equations, projection method for the double saddle points MHD equations (Schötzau in Numer Math 96:771-800, 2004) with some subtle implicit-explicit treatments for nonlinear coupled terms. For space discretization, the magnetic field is approximated by the $$H(\mathrm {curl})$$ -conforming edge element. Then, the coupled, nonlinear and saddle point type model is separated into a series of small elliptic type problems. We also rigorously prove the unconditional energy stabilities of the proposed schemes in both temporal discretization and full discretization cases. To our knowledge, it is the first totally decoupled and unconditionally stable scheme for the two-phase MHD problem. Finally, we implement some numerical examples, including accuracy tests, energy stability tests, 2D and 3D two-phase electromagnetic driven cavity benchmark problems, to verify the stability and convergence of our schemes.