Abstract

In this paper, we develop a novel second order in time, decoupled, energy stable finite element scheme for simulation of Cahn-Hilliard-Hele-Shaw system. The idea of scalar auxiliary variable approach is introduced to handle the nonlinear bulk. An operator-splitting strategy is utilized to fully decouple the coupled Cahn-Hilliard equation and Darcy equation. A full discretization is built in the framework of Galerkin finite element method. The unique solvability of numerical solution and preservation of energy law are rigorously established. Numerical experiences are recorded to illustrate the features of the designed numerical method, verify the theoretical results and conduct realistic applications.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call