Hostname: page-component-cd9895bd7-7cvxr Total loading time: 0 Render date: 2024-12-18T10:06:08.810Z Has data issue: false hasContentIssue false
  • English
  • Français

Efficient Energy Stable Schemes with Spectral Discretization in Space for Anisotropic Cahn-Hilliard Systems

Published online by Cambridge University Press:  03 June 2015

Feng Chen  and
Jie Shen
Feng Chen*
Affiliation:
Department of Mathematics, Purdue University, West Lafayette, IN 47907-1957, USA
Jie Shen*
Affiliation:
Department of Mathematics, Purdue University, West Lafayette, IN 47907-1957, USA
*
Email:[email protected]
Corresponding author.Email:[email protected]
Get access

Abstract

We develop in this paper efficient and robust numerical methods for solving anisotropic Cahn-Hilliard systems. We construct energy stable schemes for the time discretization of the highly nonlinear anisotropic Cahn-Hilliard systems by using a stabilization technique. At each time step, these schemes lead to a sequence of linear coupled elliptic equations with constant coefficients that can be efficiently solved by using a spectral-Galerkin method. We present numerical results that are consistent with earlier work on this topic, and also carry out various simulations, such as the linear bi-Laplacian regularization and the nonlinear Willmore regularization, to demonstrate the efficiency and robustness of the new schemes.

Keywords

65N3565M1235K55Spectral-Galerkinphase-fieldanisotropicCahn-Hilliardstabilizationcoupled elliptic equations
Type
Research Article
Information
Communications in Computational Physics , Volume 13 , Issue 5 , May 2013 , pp. 1189 - 1208
Copyright
Copyright © Global Science Press Limited 2013

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Allen, S. M. and Cahn, J. W.A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening. Acta Metall. Mater., 27:1085–1095,1979.Google Scholar
[2]Caffarelli, L. A. and Muler, N. E.An 1∞ bound for solutions of the Cahn-Hilliard equation. Archive for Rational Mechanics and Analysis, 133(2): 129–144, 1995.CrossRefGoogle Scholar
[3]Cahn, J. W. and Hilliard, J. E.Free energy of a nonuniform system, I: Interfacial free energy. J. Chem. Phys., 28:258,1958.Google Scholar
[4]Canuto, C., Hussaini, M. Y., Quarteroni, A., and Zang, T. A.Spectral methods. Scientific Computation. Springer-Verlag, Berlin, 2006. Fundamentals in single domains.Google Scholar
[5]Chen, F. and Shen, J.Efficient spectral-Galerkin methods for systems of coupled second-order equations and their applications. J. Comput. Phys., 231(15):5016–5028, June 2012.Google Scholar
[6]Giorgi, E. DeSome remarks on T-convergence and least squares method. In Composite media and homogenization theory (Trieste, 1990), volume 5 of Progr. Nonlinear Differential Equations Appl., pages 135–142. Birkhauser Boston, Boston, MA, 1991.Google Scholar
[7]Du, Q., Liu, C., Ryham, R., and Wang, X.A phase field formulation of the Willmore problem. Nonlinearity, 18(3):1249–1267,2005.CrossRefGoogle Scholar
[8]Eggleston, J. J., McFadden, G. B., and Voorhees, P. W.A phase-field model for highly anisotropic interfacial energy. Phys. D, 150(1-2):91–103,2001.CrossRefGoogle Scholar
[9]Eyre, D. J.Unconditionally gradient stable time marching the Cahn-Hilliard equation. In Computational and mathematical models of microstructural evolution (San Francisco, CA, 1998), volume 529 of Mater. Res. Soc. Sympos. Proc., pages 39–46. MRS, Warrendale, PA, 1998.Google Scholar
[10]Gurtin, M. E.Thermomechanics of Evolving Phase Boundaries in the Plane. Oxford University Press, USA, June 1993.CrossRefGoogle Scholar
[11]Ratz, A. and Voigt, A.Higher order regularization of anisotropic geometric evolution equations in three dimensions. Journal of Computational and Theoretical Nanoscience, 3:560–564, August 2006.Google Scholar
[12]Shen, J.Efficient spectral-Galerkin method I: direct solvers of second- and fourth-order equations using Legendre polynomials. SIAM J. Sci. Comput., 15(6):1489–1505,1994.CrossRefGoogle Scholar
[13]Shen, J., Tang, T., and Wang, L.-L.Spectral Methods: Algorithms, Analysis and Applications, volume 41 of Springer Series in Computational Mathematics. Springer, 2011.Google Scholar
[14]Shen, J., Wang, C., Wang, X., and Wise, S. M.Second-order convex splitting schemes for gradient flows with Ehrlich-Schwoebel type energy: Application to thin film epitaxy. SIAM Journal on Numerical Analysis, 50(1):105–125,2012.CrossRefGoogle Scholar
[15]Shen, J. and Yang, X.Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Discrete and Continuous Dynamical Systems, 28(4):1669–1691, June 2010.Google Scholar
[16]Spencer, B. J.Asymptotic solutions for the equilibrium crystal shape with small corner energy regularization. Physical Review E, 69(1):011603,2004.CrossRefGoogle ScholarPubMed
[17]Taylor, J. E. and Cahn, J. W.Diffuse interfaces with sharp corners and facets: phase field models with strongly anisotropic surfaces. Phys. D, 112(3-4):381–411,1998.CrossRefGoogle Scholar
[18]Torabi, S., Lowengrub, J., Voigt, A., and Wise, S.A new phase-field model for strongly anisotropic systems. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Science, 465(2105):1337–1359, May 2009.CrossRefGoogle Scholar
[19]Wise, S. M., Wang, C., and Lowengrub, J. S.An energy-stable and convergent finite-difference scheme for the phase field crystal equation. SIAM J. Numer. Anal., 47(3):2269–2288,2009.CrossRefGoogle Scholar
[20]Wise, S., Kim, J., and Lowengrub, J.Solving the regularized, strongly anisotropic Cahn-Hilliard equation by an adaptive nonlinear multigrid method. Journal of Computational Physics, 226(1):414–446, September 2007.CrossRefGoogle Scholar
[21]Yang, X., Feng, J. J., Liu, C., and Shen, J.Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method. J. Comput. Phys., 218:417–428, 2006.Google Scholar
[22]Zhu, J., Chen, L.-Q., Shen, J., and Tikare, V.Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: Application of a semi-implicit fourier spectral method. Physical Review E, 60(4):3564,1999.Google Scholar