Optimal error estimates for a fully discrete, highly efficient decoupled scheme for the 2D/3D diffuse interface two-phase MHD flows
This work provides incremental improvements in numerical methods for simulating complex fluid dynamics in engineering or physics applications.
The authors tackled the problem of simulating two-phase magnetohydrodynamics flows by deriving optimal error estimates for a fully discrete, decoupled finite element method, achieving unconditional energy stability and validating results with numerical examples.
In this paper, we derive optimal L2- and H1-norm error estimates for a fully discrete convex-splitting decoupled finite element method (FEM) for the two-phase diffuse interface magnetohydrodynamics (MHD) system. We use the semi-implicit backward Euler scheme in time and employ the standard inf-sup stable Taylor--Hood or Mini elements to discretize the velocity and pressure. Furthermore, we apply a pressure-correction scheme to decouple the velocity from the pressure. The optimal error estimates are obtained via novel Ritz and Stokes quasi-projection techniques. In addition, the unconditional energy stability of the proposed scheme is ensured. Numerical examples are presented to validate the theoretical analysis.