Hostname: page-component-74656f9884-8r8nv Total loading time: 0 Render date: 2026-11-10T02:02:12.346Z Has data issue: false hasContentIssue false
  • English
  • Français

Parameter-Free Time Adaptivity Based on Energy Evolution for the Cahn-Hilliard Equation

Part of: Elliptic equations and systems Partial differential equations, boundary value problems Equations of mathematical physics and other areas of application Miscellaneous topics - Partial differential equations Partial differential equations, initial value and time-dependent initial-boundary value problems

Published online by Cambridge University Press:  17 May 2016

Fuesheng Luo ,
Tao Tang  and
Hehu Xie
Fuesheng Luo*
Affiliation:
The Third Institute of Oceanography, SOA, Xiamen 361005, China
Tao Tang*
Affiliation:
Department of Mathematics, Southern University of Science and Technology of China, Shenzhen, Guangdong 518055, China Department of Mathematics, Hong Kong Baptist University, Kowloon Tong, Hong Kong
Hehu Xie*
Affiliation:
LSEC, ICMSEC, NCMIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
*
*Corresponding author. Email addresses:luofusheng@tio.org.cn (F. Luo), tangt@sustc.edu.cn (T. Tang), hhxie@lsec.cc.ac.cn (H. Xie)
*Corresponding author. Email addresses:luofusheng@tio.org.cn (F. Luo), tangt@sustc.edu.cn (T. Tang), hhxie@lsec.cc.ac.cn (H. Xie)
*Corresponding author. Email addresses:luofusheng@tio.org.cn (F. Luo), tangt@sustc.edu.cn (T. Tang), hhxie@lsec.cc.ac.cn (H. Xie)

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the 'Save PDF' action button.

It is known that large time-stepping method are useful for simulating phase field models. In this work, an adaptive time-stepping strategy is proposed based on numerical energy stability and equi-distribution principle. The main idea is to use the energy variation as an indicator to update the time step, so that the resulting algorithm is free of user-defined parameters, which is different from several existing approaches. Some numerical experiments are presented to illustrate the effectiveness of the algorithms.

Keywords

Adaptive time-stepping methodCahn-Hilliard equationCrank-Nicolson schemeconvex splitting method

MSC classification

Secondary: 65N15: Error bounds 35J25: Boundary value problems for second-order elliptic equations 35Q99: None of the above, but in this section 35R99: None of the above, but in this section 65M12: Stability and convergence of numerical methods 65M70: Spectral, collocation and related methods

Information

Type
Research Article
Information
Communications in Computational Physics , Volume 19 , Issue 5 , May 2016 , pp. 1542 - 1563
Copyright
Copyright © Global-Science Press 2016 

References

[1]Cahn, J.W. and Hilliard, J. E., Free energy of a non-uniform system I: Interfacial free energy, J. Chem. Phys., 28 (1958), 258–267.Google Scholar
[2]Choo, S. M., Chung, S. K. and Kim, K. I., Conservative nonlinear difference scheme for the Cahn-Hilliard equation-II, Comput. Math. Appl., 39 (2000), 229–243.Google Scholar
[3]Eyre, D., An unconditionally stable one-step scheme for gradient systems, Preprint, 1998.Google Scholar
[4]Demmel, J. W., Applied Numerical Linear Algebra. Society for Industrial and Applied Mathematics, 1997.Google Scholar
[5]Doyle, J. C., Francis, B. A. and Tannenbaum, A. R., Feedback Control Theory, Macmillan Publishing Co., 1990.Google Scholar
[6]Du, Q. and Nicolaides, R. A., Numerical analysis of a continuum model of phase transition, SIAM J. Numer. Anal., 28 (1994), 1310–1322.Google Scholar
[7]Du, Q., Tian, L. and Ju, L., Finite element approximation of the Cahn-Hilliard equation on surfaces, Comp. Meth. Appl. Mech. Engr, 200 (2011), 2458–2470.Google Scholar
[8]Elliott, C. M. and French, D. A., Numerical studies of the Cahn-Hilliard equation for phase separation, IMA J. Appl. Math., 38 (1987), 97–128.Google Scholar
[9]Elliott, C. M. and Zheng, S. M., On the Cahn-Hilliard equation, Arch. Rat. Meth. Anal., 96 (1986), 339–357.Google Scholar
[10]Cueto-Felgueroso, L. and Peraire, J., A time-adaptive finite volume method for the Cahn-Hillard and Kuramoto-Sivashinsky equations, J. Comput. Phys., 227 (2008), 9985–10017.Google Scholar
[11]Feng, W., Yu, P., Hu, S., Liu, Z., Du, Q. and Chen, L., A Fourier spectral moving mesh method for the Cahn-Hilliard equation with elasticity, Commun. Comput. Phys., 5 (2009), 582–599.Google Scholar
[12]Feng, X., Tang, T. and Yang, J., Long time numerical simulations for phase-field problems using p-adaptive spectral deferred correction methods, SIAM J. Sci. Computing, 37 (2015), A271–A294.CrossRefGoogle Scholar
[13]Franklin, G., Powell, J. D. and Emami-Naeini, A., Feedback Control of Dynamic Systems, Prentice Hall, 2006.Google Scholar
[14]Furihata, D., A stable and conservative finite difference scheme for the Cahn-Hilliard equation, Numer. Math., 87 (2001), 675–699.CrossRefGoogle Scholar
[15]Gomez, H., Hughes, T.J.R., Provably unconditionally stable, second-order time accurate, mixed variational methods for phase-field models, J. Comput. Phys., 230 (2011), 5310–5327.Google Scholar
[16]He, L. and Liu, Y., A class of stable spectral methods for the Cahn-Hillard equation, J. Comput. Phys., 228 (2009), 5101–5110.CrossRefGoogle Scholar
[17]He, Y. N., Liu, Y. X. and Tang, T., On large time-stepping methods for the Cahn-Hilliard equation, Appl. Numer. Math., 57 (2007), 616–628.CrossRefGoogle Scholar
[18]Li, B. and Liu, J. G., Thin film epitaxy with or without slope selection, European J. Appl. Math., 14 (2003), 713–743.CrossRefGoogle Scholar
[19]Novick-Cohen, A. and Segel, L. A., Nonlinear aspects of the Cahn-Hilliard equation, Phys. D, 10 (1984), 277–298.Google Scholar
[20]Qiao, Z. H., Zhang, Z. R., and Tang, T., An adaptive time-stepping strategy for the molecular beam epitaxy models, SIAM. J. Sci. Comput., 33 (2011), 1395–1414.Google Scholar
[21]Shen, J., Tang, T. and Wang, L., Spectral Methods: Algorithms, Analysis and Applications, Springer-Verlag, Berlin Heidelberg, 2011.CrossRefGoogle Scholar
[22]Shen, J., Wang, C., Wang, X. and Wise, S., Second-order convex splitting schemes for gradient flows with Ehrlich-Schwoebel type energy: Application to thin film epitaxy, SIAM J. Numer. Anal., 50 (2012), 105–125.Google Scholar
[23]Shen, J. and Yang, X., Numerical approximations of Allen-Cahn and Cahn-Hilliard equations, DCDS, Series A, 28 (2010), 1669–1691.Google Scholar
[24]Söderlind, G., Automatic control and adaptive time-stepping, Numer. Algor., 31 (2001), 281–310.Google Scholar
[25]Sun, Z. Z., A second-order accurate linearized difference scheme for the two-dimensional Cahn-Hilliard equation, Math. Comput., 64 (1995), 1463–1471.Google Scholar
[26]Wells, G., Kuhl, E. and Garikipati, K., A discontinuous Galerkin method for the Cahn-Hillard equation, J. Comput. Phys., 218 (2006), 860–877.Google Scholar
[27]Wise, S. M., Unconditionally stable finite difference, nonlinear multigrid simulation of the Cahn-Hillard-Hele-Shaw system of equations, J. Sci. Comput., 44 (2010), 38–68.Google Scholar
[28]Xu, C. J. and Tang, T., Stability analysis of large time-stepping methods for epitaxial growth models, SIAM J. Numer. Anal., 44 (2006), 1759–1779.Google Scholar
[29]Zhang, Z. R. and Qiao, Z. H., An adaptive time-stepping strategy for the Cahn-Hilliard equation, Commun. Comput. Phys., 11 (2012), 1261–1278.Google Scholar
[30]Zhu, J., Chen, L., Shen, J. and Tikare, V., Coarsening kinetics from a variable mobility Cahn-Hilliard equation-application of semi-implicit Fourier spectral method, Phys. Review E, 60(4) (1999), 3564–3572.Google Scholar