A decoupled finite element method for a modified Cahn-Hilliard-Hele-Shaw system

In this paper, a decoupled finite element method for a modified Cahn-Hilliard-Hele-Shaw system with double well potential is constructed. The proposed scheme is based on the convex splitting of the energy functional in time and uses the mixed finite element discretzation in space. The computation of the velocity $ \mathbf{u} $ is separated from the computation of the pressure $ p $ by using an operator-splitting strategy. A Possion equation is solved to update the pressure at each time step. Unconditional stability and error estimates are analyzed in detail theoretically. Furthermore, the theoretical part is verified by several numerical examples, whose results show that the numerical examples are consistent with the results of the theoretical part.