Mathematics > Numerical Analysis
[Submitted on 5 May 2011 (v1), last revised 25 Apr 2013 (this version, v2)]
Title:Analysis of a Darcy-Cahn-Hilliard Diffuse Interface Model for the Hele-Shaw Flow and its Fully Discrete Finite Element Approximation
View PDFAbstract:In this paper we present PDE and finite element analyses for a system of partial differential equations (PDEs) consisting of the Darcy equation and the Cahn-Hilliard equation, which arises as a diffuse interface model for the two phase Hele-Shaw flow. We propose a fully discrete implicit finite element method for approximating the PDE system, which consists of the implicit Euler method combined with a convex splitting energy strategy for the temporal discretization, the standard finite element discretization for the pressure and a split (or mixed) finite element discretization for the fourth order Cahn-Hilliard equation. It is shown that the proposed numerical method satisfies a mass conservation law in addition to a discrete energy law that mimics the basic energy law for the Darcy-Cahn-Hilliard phase field model and holds uniformly in the phase field parameter $\epsilon$. With help of the discrete energy law, we first prove that the fully discrete finite method is unconditionally energy stable and uniquely solvable at each time step. We then show that, using the compactness method, the finite element solution has an accumulation point that is a weak solution of the PDE system. As a result, the convergence result also provides a constructive proof of the existence of global-in-time weak solutions to the Darcy-Cahn-Hilliard phase field model in both two and three dimensions. Finally, we propose a nonlinear multigrid iterative algorithm to solve the finite element equations at each time step. Numerical experiments based on the overall solution method of combining the proposed finite element discretization and the nonlinear multigrid solver are presented to validate the theoretical results and to show the effectiveness of the proposed fully discrete finite element method for approximating the Darcy-Cahn-Hilliard phase field model.
Submission history
From: Xiaobing Feng Dr. [view email][v1] Thu, 5 May 2011 19:10:36 UTC (288 KB)
[v2] Thu, 25 Apr 2013 16:13:24 UTC (545 KB)
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.