American Institute of Mathematical Sciences

Advanced
November  2013, 18(9): 2211-2238. doi: 10.3934/dcdsb.2013.18.2211

A mixed discontinuous Galerkin, convex splitting scheme for a modified Cahn-Hilliard equation and an efficient nonlinear multigrid solver

Andreas C. Aristotelous 1, , Ohannes Karakashian 2, and  Steven M. Wise 3,

1. 

Statistical and Applied Mathematical Sciences Institute (SAMSI), Research Triangle Park, NC 27709, United States
2. 

Department of Mathematics, University of Tennessee, Knoxville, TN 37996, United States
3. 

Department of Mathematics, The University of Tennessee, Knoxville, TN 37996-0612, United States

Received  June 2013 Revised  August 2013 Published  September 2013

In this paper we devise and analyze a mixed discontinuous Galerkin finite element method for a modified Cahn-Hilliard equation that models phase separation in diblock copolymer melts. The time discretization is based on a convex splitting of the energy of the equation. We prove that our scheme is unconditionally energy stable with respect to a spatially discrete analogue of the continuous free energy of the system, unconditionally uniquely solvable, and convergent in the natural energy norm with optimal rates. We describe an efficient nonlinear multigrid solver for advancing our semi-implicit scheme in time and conclude the paper with numerical tests confirming the predicted convergence rates and suggesting the near-optimal complexity of the solver.
Keywords: convergence, discontinuous Galerkin, diblock copolymer, symmetric interior penalty, nonlinear multigrid., convex splitting, Cahn-Hilliard equation, energy stability.
Mathematics Subject Classification: Primary: 65M60, 65M12; Secondary: 65M5.
Citation: Andreas C. Aristotelous, Ohannes Karakashian, Steven M. Wise. A mixed discontinuous Galerkin, convex splitting scheme for a modified Cahn-Hilliard equation and an efficient nonlinear multigrid solver. Discrete & Continuous Dynamical Systems - B, 2013, 18 (9) : 2211-2238. doi: 10.3934/dcdsb.2013.18.2211
References:
[1]

A. Aristotelous, "Adaptive Discontinuous Galerkin Finite Element Methods for a Diffuse Interface Model of Biological Growth,", PhD thesis, (2011).   Google Scholar

[2]

D. N. Arnold, An interior penalty finite element method with discontinuous elements,, SIAM Journal on Numerical Analysis, 19 (1982), 742.  doi: 10.1137/0719052.  Google Scholar

[3]

D. N. Arnold, F. Brezzi, B. Cockburn and L. D. Marini, Unified analysis of discontinuous Galerkin methods for elliptic problems,, SIAM Journal on Numerical Analysis, 39 (2002), 1749.  doi: 10.1137/S0036142901384162.  Google Scholar

[4]

I. Babuška and M. Zlámal, Nonconforming elements in the finite element method with penalty,, SIAM Journal on Numerical Analysis, 10 (1973), 863.  doi: 10.1137/0710071.  Google Scholar

[5]

G. A. Baker, Finite element methods for elliptic equations using nonconforming elements,, Math. Comp., (1977), 45.  doi: 10.1090/S0025-5718-1977-0431742-5.  Google Scholar

[6]

A. Baskaran, Z. Hu, J. S. Lowengrub, C. Wang, S. M. Wise and P. Zhou, Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation,, Journal of Computational Physics, 250 (2013), 270.  doi: 10.1016/j.jcp.2013.04.024.  Google Scholar

[7]

J. H. Bramble, "Multigrid Methods,", Research Notes in Mathematics Series. Chapman and Hall/CRC, (1993).   Google Scholar

[8]

J. W. Cahn, On spinodal decomposition,, Acta Metallurgica, 9 (1961), 795.  doi: 10.1016/0001-6160(61)90182-1.  Google Scholar

[9]

R. Choksi, M. Maras and J. F. Williams, 2D phase diagram for minimizers of a Cahn-Hilliard functional with long-range interactions,, SIAM Journal on Applied Dynamical Systems, 10 (2011), 1344.  doi: 10.1137/100784497.  Google Scholar

[10]

R. Choksi and X. Ren, On a derivation of a density functional theory for microphase separation of diblock copolymers,, Journal of Statistical Physics, 113 (2003), 151.  doi: 10.1023/A:1025722804873.  Google Scholar

[11]

P. G. Ciarlet, "Introduction to Numerical Linear Algebra and Optimisation,", Cambridge University Press, (1989).   Google Scholar

[12]

C. Collins, J. Shen and S. M. Wise, Unconditionally stable finite difference multigrid schemes for the Cahn-Hilliard-Brinkman equation,, Commun. Comput. Phys., 13 (2013), 929.  doi: 10.4208/cicp.171211.130412a.  Google Scholar

[13]

J. Douglas and T. Dupont, Interior penalty procedures for elliptic and parabolic Galerkin methods,, In, (1976), 207.   Google Scholar

[14]

C. M. Elliott and S. Zheng, On the Cahn-Hilliard equation,, Archive for Rational Mechanics and Analysis, 96 (1986), 339.  doi: 10.1007/BF00251803.  Google Scholar

[15]

C. M. Elliott, The Cahn-Hilliard model for the kinetics of phase separation,, In J.F. Rodrigues, (1988), 1.   Google Scholar

[16]

D. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation,, In J. W. Bullard, (1998), 1686.   Google Scholar

[17]

X. Feng and O. A. Karakashian, Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn-Hilliard equation of phase transition,, Math. Comput., 76 (2007), 1093.  doi: 10.1090/S0025-5718-07-01985-0.  Google Scholar

[18]

X. Feng and S. M. Wise, Analysis of a Darcy-Cahn-Hilliard diffuse interface model for the Hele-Shaw flow and its fully discrete finite element approximation,, SIAM Journal on Numerical Analysis, 50 (2012), 1320.  doi: 10.1137/110827119.  Google Scholar

[19]

J. Gopalakrishnan and G. Kanschat, A multilevel discontinuous Galerkin method,, Numerische Mathematik, 95 (2003), 527.  doi: 10.1007/s002110200392.  Google Scholar

[20]

W. Hackbusch, "Multi-Grid Methods and Applications,", Springer Series in Computational Mathematics. Springer, (2010).   Google Scholar

[21]

M. R. Hanisch, Multigrid preconditioning for the biharmonic Dirichlet problem,, SIAM Journal on Numerical Analysis, 30 (1993), 184.  doi: 10.1137/0730009.  Google Scholar

[22]

Z. Hu, S. M. Wise, C. Wang and J. S. Lowengrub, Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation,, Journal of Computational Physics, 228 (2009), 5323.  doi: 10.1016/j.jcp.2009.04.020.  Google Scholar

[23]

O. A. Karakashian and W. N. Jureidini, A nonconforming finite element method for the stationary Navier-Stokes equations,, SIAM Journal on Numerical Analysis, 35 (1998), 93.  doi: 10.1137/S0036142996297199.  Google Scholar

[24]

D. Kay, V. Styles and E. Süli, Discontinuous Galerkin finite element approximation of the Cahn-Hilliard equation with convection,, SIAM Journal on Numerical Analysis, 47 (2009), 2660.  doi: 10.1137/080726768.  Google Scholar

[25]

T. Ohta and K. Kawasaki, Equilibrium morphology of block copolymer melts,, Macromolecules, 19 (1986), 2621.  doi: 10.1021/ma00164a028.  Google Scholar

[26]

P. Percell and M. F. Wheeler, A local residual finite element procedure for elliptic equations,, SIAM Journal on Numerical Analysis, 15 (1978), 705.  doi: 10.1137/0715047.  Google Scholar

[27]

U. Trottenberg, C. W. Oosterlee and A. Schüller, "Multigrid,", Academic Press, (2005).   Google Scholar

[28]

C. Wang, X. Wang and S. M. Wise, Unconditionally stable schemes for equations of thin film epitaxy,, Discrete and Continuous Dynamical Systems - Series A (DCDS-A), 28 (2010), 405.  doi: 10.3934/dcds.2010.28.405.  Google Scholar

[29]

C. Wang and S. M. Wise, An energy stable and convergent finite-difference scheme for the modified phase field crystal equation,, SIAM Journal on Numerical Analysis, 49 (2011), 945.  doi: 10.1137/090752675.  Google Scholar

[30]

G. N. Wells, E. Kuhl and K. Garikipati, A discontinuous Galerkin method for the Cahn-Hilliard equation,, Journal of Computational Physics, 218 (2006), 860.  doi: 10.1016/j.jcp.2006.03.010.  Google Scholar

[31]

M. F. Wheeler, An elliptic collocation-finite element method with interior penalties,, SIAM Journal on Numerical Analysis, 15 (1978), 152.  doi: 10.1137/0715010.  Google Scholar

[32]

S. M. Wise, Unconditionally stable finite difference, nonlinear multigrid simulation of the Cahn-Hilliard-Hele-Shaw system of equations,, Journal of Scientific Computing, 44 (2010), 38.  doi: 10.1007/s10915-010-9363-4.  Google Scholar

[33]

S. M. Wise, C. Wang and J. S. Lowengrub, An energy-stable and convergent finite-difference scheme for the phase field crystal equation,, SIAM Journal on Numerical Analysis, 47 (2009), 2269.  doi: 10.1137/080738143.  Google Scholar

[34]

Y. Xia, Y. Xu and C. W. Shu, Local discontinuous Galerkin methods for the Cahn-Hilliard type equations,, Journal of Computational Physics, 227 (2007), 472.  doi: 10.1016/j.jcp.2007.08.001.  Google Scholar

show all references

References:
[1]

A. Aristotelous, "Adaptive Discontinuous Galerkin Finite Element Methods for a Diffuse Interface Model of Biological Growth,", PhD thesis, (2011).   Google Scholar

[2]

D. N. Arnold, An interior penalty finite element method with discontinuous elements,, SIAM Journal on Numerical Analysis, 19 (1982), 742.  doi: 10.1137/0719052.  Google Scholar

[3]

D. N. Arnold, F. Brezzi, B. Cockburn and L. D. Marini, Unified analysis of discontinuous Galerkin methods for elliptic problems,, SIAM Journal on Numerical Analysis, 39 (2002), 1749.  doi: 10.1137/S0036142901384162.  Google Scholar

[4]

I. Babuška and M. Zlámal, Nonconforming elements in the finite element method with penalty,, SIAM Journal on Numerical Analysis, 10 (1973), 863.  doi: 10.1137/0710071.  Google Scholar

[5]

G. A. Baker, Finite element methods for elliptic equations using nonconforming elements,, Math. Comp., (1977), 45.  doi: 10.1090/S0025-5718-1977-0431742-5.  Google Scholar

[6]

A. Baskaran, Z. Hu, J. S. Lowengrub, C. Wang, S. M. Wise and P. Zhou, Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation,, Journal of Computational Physics, 250 (2013), 270.  doi: 10.1016/j.jcp.2013.04.024.  Google Scholar

[7]

J. H. Bramble, "Multigrid Methods,", Research Notes in Mathematics Series. Chapman and Hall/CRC, (1993).   Google Scholar

[8]

J. W. Cahn, On spinodal decomposition,, Acta Metallurgica, 9 (1961), 795.  doi: 10.1016/0001-6160(61)90182-1.  Google Scholar

[9]

R. Choksi, M. Maras and J. F. Williams, 2D phase diagram for minimizers of a Cahn-Hilliard functional with long-range interactions,, SIAM Journal on Applied Dynamical Systems, 10 (2011), 1344.  doi: 10.1137/100784497.  Google Scholar

[10]

R. Choksi and X. Ren, On a derivation of a density functional theory for microphase separation of diblock copolymers,, Journal of Statistical Physics, 113 (2003), 151.  doi: 10.1023/A:1025722804873.  Google Scholar

[11]

P. G. Ciarlet, "Introduction to Numerical Linear Algebra and Optimisation,", Cambridge University Press, (1989).   Google Scholar

[12]

C. Collins, J. Shen and S. M. Wise, Unconditionally stable finite difference multigrid schemes for the Cahn-Hilliard-Brinkman equation,, Commun. Comput. Phys., 13 (2013), 929.  doi: 10.4208/cicp.171211.130412a.  Google Scholar

[13]

J. Douglas and T. Dupont, Interior penalty procedures for elliptic and parabolic Galerkin methods,, In, (1976), 207.   Google Scholar

[14]

C. M. Elliott and S. Zheng, On the Cahn-Hilliard equation,, Archive for Rational Mechanics and Analysis, 96 (1986), 339.  doi: 10.1007/BF00251803.  Google Scholar

[15]

C. M. Elliott, The Cahn-Hilliard model for the kinetics of phase separation,, In J.F. Rodrigues, (1988), 1.   Google Scholar

[16]

D. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation,, In J. W. Bullard, (1998), 1686.   Google Scholar

[17]

X. Feng and O. A. Karakashian, Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn-Hilliard equation of phase transition,, Math. Comput., 76 (2007), 1093.  doi: 10.1090/S0025-5718-07-01985-0.  Google Scholar

[18]

X. Feng and S. M. Wise, Analysis of a Darcy-Cahn-Hilliard diffuse interface model for the Hele-Shaw flow and its fully discrete finite element approximation,, SIAM Journal on Numerical Analysis, 50 (2012), 1320.  doi: 10.1137/110827119.  Google Scholar

[19]

J. Gopalakrishnan and G. Kanschat, A multilevel discontinuous Galerkin method,, Numerische Mathematik, 95 (2003), 527.  doi: 10.1007/s002110200392.  Google Scholar

[20]

W. Hackbusch, "Multi-Grid Methods and Applications,", Springer Series in Computational Mathematics. Springer, (2010).   Google Scholar

[21]

M. R. Hanisch, Multigrid preconditioning for the biharmonic Dirichlet problem,, SIAM Journal on Numerical Analysis, 30 (1993), 184.  doi: 10.1137/0730009.  Google Scholar

[22]

Z. Hu, S. M. Wise, C. Wang and J. S. Lowengrub, Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation,, Journal of Computational Physics, 228 (2009), 5323.  doi: 10.1016/j.jcp.2009.04.020.  Google Scholar

[23]

O. A. Karakashian and W. N. Jureidini, A nonconforming finite element method for the stationary Navier-Stokes equations,, SIAM Journal on Numerical Analysis, 35 (1998), 93.  doi: 10.1137/S0036142996297199.  Google Scholar

[24]

D. Kay, V. Styles and E. Süli, Discontinuous Galerkin finite element approximation of the Cahn-Hilliard equation with convection,, SIAM Journal on Numerical Analysis, 47 (2009), 2660.  doi: 10.1137/080726768.  Google Scholar

[25]

T. Ohta and K. Kawasaki, Equilibrium morphology of block copolymer melts,, Macromolecules, 19 (1986), 2621.  doi: 10.1021/ma00164a028.  Google Scholar

[26]

P. Percell and M. F. Wheeler, A local residual finite element procedure for elliptic equations,, SIAM Journal on Numerical Analysis, 15 (1978), 705.  doi: 10.1137/0715047.  Google Scholar

[27]

U. Trottenberg, C. W. Oosterlee and A. Schüller, "Multigrid,", Academic Press, (2005).   Google Scholar

[28]

C. Wang, X. Wang and S. M. Wise, Unconditionally stable schemes for equations of thin film epitaxy,, Discrete and Continuous Dynamical Systems - Series A (DCDS-A), 28 (2010), 405.  doi: 10.3934/dcds.2010.28.405.  Google Scholar

[29]

C. Wang and S. M. Wise, An energy stable and convergent finite-difference scheme for the modified phase field crystal equation,, SIAM Journal on Numerical Analysis, 49 (2011), 945.  doi: 10.1137/090752675.  Google Scholar

[30]

G. N. Wells, E. Kuhl and K. Garikipati, A discontinuous Galerkin method for the Cahn-Hilliard equation,, Journal of Computational Physics, 218 (2006), 860.  doi: 10.1016/j.jcp.2006.03.010.  Google Scholar

[31]

M. F. Wheeler, An elliptic collocation-finite element method with interior penalties,, SIAM Journal on Numerical Analysis, 15 (1978), 152.  doi: 10.1137/0715010.  Google Scholar

[32]

S. M. Wise, Unconditionally stable finite difference, nonlinear multigrid simulation of the Cahn-Hilliard-Hele-Shaw system of equations,, Journal of Scientific Computing, 44 (2010), 38.  doi: 10.1007/s10915-010-9363-4.  Google Scholar

[33]

S. M. Wise, C. Wang and J. S. Lowengrub, An energy-stable and convergent finite-difference scheme for the phase field crystal equation,, SIAM Journal on Numerical Analysis, 47 (2009), 2269.  doi: 10.1137/080738143.  Google Scholar

[34]

Y. Xia, Y. Xu and C. W. Shu, Local discontinuous Galerkin methods for the Cahn-Hilliard type equations,, Journal of Computational Physics, 227 (2007), 472.  doi: 10.1016/j.jcp.2007.08.001.  Google Scholar

[1]

Tomáš Roubíček. Cahn-Hilliard equation with capillarity in actual deforming configurations. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 41-55. doi: 10.3934/dcdss.2020303
[2]

Leilei Wei, Yinnian He. A fully discrete local discontinuous Galerkin method with the generalized numerical flux to solve the tempered fractional reaction-diffusion equation. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020319
[3]

Zexuan Liu, Zhiyuan Sun, Jerry Zhijian Yang. A numerical study of superconvergence of the discontinuous Galerkin method by patch reconstruction. Electronic Research Archive, 2020, 28 (4) : 1487-1501. doi: 10.3934/era.2020078
[4]

Yue Feng, Yujie Liu, Ruishu Wang, Shangyou Zhang. A conforming discontinuous Galerkin finite element method on rectangular partitions. Electronic Research Archive, , () : -. doi: 10.3934/era.2020120
[5]

Gang Luo, Qingzhi Yang. The point-wise convergence of shifted symmetric higher order power method. Journal of Industrial & Management Optimization, 2021, 17 (1) : 357-368. doi: 10.3934/jimo.2019115
[6]

Pierluigi Colli, Gianni Gilardi, Jürgen Sprekels. Deep quench approximation and optimal control of general Cahn–Hilliard systems with fractional operators and double obstacle potentials. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 243-271. doi: 10.3934/dcdss.2020213
[7]

Hirokazu Ninomiya. Entire solutions of the Allen–Cahn–Nagumo equation in a multi-dimensional space. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 395-412. doi: 10.3934/dcds.2020364
[8]

Teresa D'Aprile. Bubbling solutions for the Liouville equation around a quantized singularity in symmetric domains. Communications on Pure & Applied Analysis, 2021, 20 (1) : 159-191. doi: 10.3934/cpaa.2020262
[9]

Oussama Landoulsi. Construction of a solitary wave solution of the nonlinear focusing schrödinger equation outside a strictly convex obstacle in the $ L^2 $-supercritical case. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 701-746. doi: 10.3934/dcds.2020298
[10]

Eduard Feireisl, Elisabetta Rocca, Giulio Schimperna, Arghir Zarnescu. Weak sequential stability for a nonlinear model of nematic electrolytes. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 219-241. doi: 10.3934/dcdss.2020366
[11]

Gervy Marie Angeles, Gilbert Peralta. Energy method for exponential stability of coupled one-dimensional hyperbolic PDE-ODE systems. Evolution Equations & Control Theory, 2020  doi: 10.3934/eect.2020108
[12]

Scipio Cuccagna, Masaya Maeda. A survey on asymptotic stability of ground states of nonlinear Schrödinger equations II. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020450
[13]

Lei Liu, Li Wu. Multiplicity of closed characteristics on $ P $-symmetric compact convex hypersurfaces in $ \mathbb{R}^{2n} $. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020378
[14]

Reza Chaharpashlou, Abdon Atangana, Reza Saadati. On the fuzzy stability results for fractional stochastic Volterra integral equation. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020432
[15]

Shengxin Zhu, Tongxiang Gu, Xingping Liu. AIMS: Average information matrix splitting. Mathematical Foundations of Computing, 2020, 3 (4) : 301-308. doi: 10.3934/mfc.2020012
[16]

Hassan Mohammad. A diagonal PRP-type projection method for convex constrained nonlinear monotone equations. Journal of Industrial & Management Optimization, 2021, 17 (1) : 101-116. doi: 10.3934/jimo.2019101
[17]

Xin-Guang Yang, Lu Li, Xingjie Yan, Ling Ding. The structure and stability of pullback attractors for 3D Brinkman-Forchheimer equation with delay. Electronic Research Archive, 2020, 28 (4) : 1395-1418. doi: 10.3934/era.2020074
[18]

Xuefei He, Kun Wang, Liwei Xu. Efficient finite difference methods for the nonlinear Helmholtz equation in Kerr medium. Electronic Research Archive, 2020, 28 (4) : 1503-1528. doi: 10.3934/era.2020079
[19]

Thierry Cazenave, Ivan Naumkin. Local smooth solutions of the nonlinear Klein-gordon equation. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020448
[20]

Zuliang Lu, Fei Huang, Xiankui Wu, Lin Li, Shang Liu. Convergence and quasi-optimality of $ L^2- $norms based an adaptive finite element method for nonlinear optimal control problems. Electronic Research Archive, 2020, 28 (4) : 1459-1486. doi: 10.3934/era.2020077

2019 Impact Factor: 1.27

Tools

Metrics

  • PDF downloads (154)
  • HTML views (0)
  • Cited by (27)

Other articles
by authors

[Back to Top]

Copyright © 2020 American Institute of Mathematical Sciences