A conservative numerical method for the Cahn-Hilliard equation in complex domains

We propose an efficient finite difference scheme for solving the Cahn-Hilliard equation with a variable mobility in complex domains. Our method employs a type of unconditionally gradient stable splitting discretization. We also extend the scheme to compute the Cahn-Hilliard equation in arbitrarily shaped domains. We prove the mass conservation property of the proposed discrete scheme for complex domains. The resulting discretized equations are solved using a multigrid method. Numerical simulations are presented to demonstrate that the proposed scheme can deal with complex geometries robustly. Furthermore, the multigrid efficiency is retained even if the embedded domain is present.

[1]  Wei Shyy,et al.  COMPUTATIONS OF DROP DYNAMICS WITH THE IMMERSED BOUNDARY METHOD, PART 2: DROP IMPACT AND HEAT TRANSFER , 2003 .

[2]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[3]  J. E. Hilliard,et al.  Spinodal decomposition: A reprise , 1971 .

[4]  Harald Garcke,et al.  Finite Element Approximation of the Cahn-Hilliard Equation with Degenerate Mobility , 1999, SIAM J. Numer. Anal..

[5]  Richard Welford,et al.  A multigrid finite element solver for the Cahn-Hilliard equation , 2006, J. Comput. Phys..

[6]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[7]  C. M. Elliott,et al.  Numerical Studies of the Cahn-Hilliard Equation for Phase Separation , 1987 .

[8]  Junseok Kim,et al.  A numerical method for the Cahn–Hilliard equation with a variable mobility , 2007 .

[9]  D. Calhoun.  A Cartesian Grid Method for Solving the Two-Dimensional Streamfunction-Vorticity Equations in Irregular Regions , 2002 .

[10]  Yan Xu,et al.  Local discontinuous Galerkin methods for the Cahn-Hilliard type equations , 2007, J. Comput. Phys..

[11]  J. Lowengrub,et al.  Conservative multigrid methods for Cahn-Hilliard fluids , 2004 .

[12]  Daisuke Furihata,et al.  A stable and conservative finite difference scheme for the Cahn-Hilliard equation , 2001, Numerische Mathematik.

[13]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy and Free Energy of a Nonuniform System. III. Nucleation in a Two‐Component Incompressible Fluid , 2013 .

[14]  Peter W Voorhees,et al.  A phase-field model for highly anisotropic interfacial energy , 2001 .

[15]  Yinnian He,et al.  On large time-stepping methods for the Cahn--Hilliard equation , 2007 .

[16]  Junseok Kim,et al.  A numerical method for the ternary Cahn--Hilliard system with a degenerate mobility , 2009 .

[17]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[18]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[19]  Qiang Du,et al.  Coarsening Kinetics of a Two Phase Mixture with Highly Disparate Diffusion Mobility , 2010 .

[20]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[21]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[22]  Zhi-zhong Sun,et al.  A second-order accurate linearized difference scheme for the two-dimensional Cahn-Hilliard equation , 1995 .

[23]  D. Weitz,et al.  Geometrically mediated breakup of drops in microfluidic devices. , 2003, Physical review letters.

[24]  E. Mello,et al.  Numerical study of the Cahn–Hilliard equation in one, two and three dimensions , 2004, cond-mat/0410772.

[25]  J. Taylor,et al.  Overview no. 113 surface motion by surface diffusion , 1994 .

[26]  Wei Shyy,et al.  COMPUTATIONS OF DROP DYNAMICS WITH THE IMMERSED BOUNDARY METHOD, PART 1: NUMERICAL ALGORITHM AND BUOYANCY-INDUCED EFFECT , 2003 .

[27]  Charles M. Elliott,et al.  Kinetics of phase decomposition processes: numerical solutions to Cahn–Hilliard equation , 1990 .

[28]  Fengshan Bai,et al.  Numerical computations of coarsening in the one-dimensional Cahn-Hilliard model of phase separation , 1994 .

[29]  John W. Barrett,et al.  Finite element approximation of the Cahn-Hilliard equation with concentration dependent mobility , 1999, Math. Comput..

[30]  Glenn Ricart,et al.  The MOSAIC Internet browser , 1994 .

[31]  Krishna Garikipati,et al.  A discontinuous Galerkin method for the Cahn-Hilliard equation , 2006, J. Comput. Phys..

[32]  Marcus S. Day,et al.  Embedded Boundary Algorithms for Solving the Poisson Equation on Complex Domains , 1998 .

[33]  Charles M. Elliott,et al.  The Cahn–Hilliard equation with a concentration dependent mobility: motion by minus the Laplacian of the mean curvature , 1996, European Journal of Applied Mathematics.

[34]  D. J. Eyre,et al.  An Unconditionally Stable One-Step Scheme for Gradient Systems , 1997 .

[35]  Xingde Ye,et al.  The Legendre collocation method for the Cahn-Hilliard equation , 2003 .

[36]  Daisuke Furihata,et al.  A stable, convergent, conservative and linear finite difference scheme for the Cahn-Hilliard equation , 2003 .