Decoupled, linear, unconditionally energy stable and charge-conservative finite element method for a inductionless magnetohydrodynamic phase-field model

TL;DR

This paper introduces a novel linear, decoupled finite element method for a complex inductionless magnetohydrodynamics phase-field model, ensuring energy stability and charge conservation, with demonstrated accuracy and efficiency.

Contribution

The paper presents a new linear, decoupled finite element scheme that is unconditionally energy stable and charge-conservative for a multi-physics MHD phase-field model.

Findings

The scheme is unconditionally energy stable.

It conserves mass and charge.

Numerical experiments confirm accuracy and efficiency.

Abstract

In this paper, we consider the numerical approximation for a diffuse interface model of the two-phase incompressible inductionless magnetohydrodynamics problem. This model consists of Cahn-Hilliard equations, Navier-Stokes equations and Poisson equation. We propose a linear and decoupled finite element method to solve this highly nonlinear and multi-physics system. For the time variable, the discretization is a combination of first-order Euler semi-implicit scheme, several first-order stabilization terms and implicit-explicit treatments for coupling terms. For the space variables, we adopt the finite element discretization, especially, we approximate the current density and electric potential by inf-sup stable face-volume mixed finite element pairs. With these techniques, the scheme only involves a sequence of decoupled linear equations to solve at each time step. We show that the scheme is provably mass-conservative, charge-conservative and unconditionally energy stable. Numerical experiments are performed to illustrate the features, accuracy and efficiency of the proposed scheme.