The slope-selection epitaxial thin-film growth model was formulated as the $ L^2 $-gradient flow of an energy functional with a nonlinear potential of the surface slope. We developed a second-order, linear, and energy-stable scheme for this model based on a linear convex splitting in which the nonlinear potential term was treated explicitly together with an auxiliary term that ensured the convexity of the explicit part. The scheme was constructed using a second-order strong-stability-preserving implicit–explicit Runge–Kutta method. We explicitly identified the range of Runge–Kutta coefficients for which the original discrete energy decay property held, and proved that the scheme was unconditionally energy-stable with respect to the original discrete energy functional. Numerical results were presented to verify the accuracy, computational efficiency, and long-time energy stability of the scheme.
Citation: Hyun Geun Lee. A second-order linear energy-stable scheme for slope-selection epitaxial thin-film growth[J]. AIMS Mathematics, 2026, 11(4): 9989-10003. doi: 10.3934/math.2026412
The slope-selection epitaxial thin-film growth model was formulated as the $ L^2 $-gradient flow of an energy functional with a nonlinear potential of the surface slope. We developed a second-order, linear, and energy-stable scheme for this model based on a linear convex splitting in which the nonlinear potential term was treated explicitly together with an auxiliary term that ensured the convexity of the explicit part. The scheme was constructed using a second-order strong-stability-preserving implicit–explicit Runge–Kutta method. We explicitly identified the range of Runge–Kutta coefficients for which the original discrete energy decay property held, and proved that the scheme was unconditionally energy-stable with respect to the original discrete energy functional. Numerical results were presented to verify the accuracy, computational efficiency, and long-time energy stability of the scheme.
| [1] |
M. Ortiz, E.A. Repetto, H. Si, A continuum model of kinetic roughening and coarsening in thin films, J. Mech. Phys. Solids, 47 (1999), 697–730. https://doi.org/10.1016/S0022-5096(98)00102-1 doi: 10.1016/S0022-5096(98)00102-1
|
| [2] |
D. Moldovan, L. Golubović, Interfacial coarsening dynamics in epitaxial growth with slope selection, Phys. Rev. E, 61 (2000), 6190–6214. https://doi.org/10.1103/PhysRevE.61.6190 doi: 10.1103/PhysRevE.61.6190
|
| [3] |
B. Li, J.-G. Liu, Thin film epitaxy with or without slope selection, Eur. J. Appl. Math., 14 (2003), 713–743. https://doi.org/10.1017/S0956792503003764 doi: 10.1017/S0956792503003764
|
| [4] |
J. Shen, C. Wang, X. Wang, S. M. Wise, 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. https://doi.org/10.1137/110837841 doi: 10.1137/110837841
|
| [5] |
D.J. Eyre, Unconditionally gradient stable time marching the Cahn–Hilliard equation, MRS Proc., 529 (1998), 39–46. https://doi.org/10.1557/PROC-529-39 doi: 10.1557/PROC-529-39
|
| [6] |
Q. Du, R. A. Nicolaides, Numerical analysis of a continuum model of phase transition, SIAM J. Numer. Anal., 28 (1991), 1310–1322. https://doi.org/10.1137/0728091 doi: 10.1137/0728091
|
| [7] |
X. Yang, J. Zhao, Q. Wang, Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method, J. Comput. Phys., 333 (2017), 104–127. https://doi.org/10.1016/j.jcp.2016.12.040 doi: 10.1016/j.jcp.2016.12.040
|
| [8] |
X. Yang, Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends, J. Comput. Phys., 327 (2016), 294–316. https://doi.org/10.1016/j.jcp.2016.08.011 doi: 10.1016/j.jcp.2016.08.011
|
| [9] |
W. Feng, C. Wang, S. M. Wise, Z. Zhang, A second-order energy stable backward differentiation formula method for the epitaxial thin film equation with slope selection, Numer. Methods Partial Differential Eq., 34 (2018), 1975–2007. https://doi.org/10.1002/num.22378 doi: 10.1002/num.22378
|
| [10] |
Q. Cheng, J. Shen, X. Yang, Highly efficient and accurate numerical schemes for the epitaxial thin film growth models by using the SAV approach, J. Sci. Comput., 78 (2019), 1467–1487. https://doi.org/10.1007/s10915-018-0853-6 doi: 10.1007/s10915-018-0853-6
|
| [11] |
J. Shen, J. Xu, J. Yang, The scalar auxiliary variable (SAV) approach for gradient flows, J. Comput. Phys., 353 (2018), 407–416. https://doi.org/10.1016/j.jcp.2017.10.066 doi: 10.1016/j.jcp.2017.10.066
|
| [12] |
K. Glasner, S. Orizaga, Improving the accuracy of convexity splitting methods for gradient flow equations, J. Comput. Phys., 315 (2016), 52–64. https://doi.org/10.1016/j.jcp.2016.03.042 doi: 10.1016/j.jcp.2016.03.042
|
| [13] |
S. Orizaga, T. Witelski, IMEX methods for thin-film equations and Cahn–Hilliard equations with variable mobility, Comput. Mater. Sci. 243 (2024), 113145. https://doi.org/10.1016/j.commatsci.2024.113145 doi: 10.1016/j.commatsci.2024.113145
|
| [14] |
H. Song, Energy SSP-IMEX Runge–Kutta methods for the Cahn–Hilliard equation, J. Comput. Appl. Math., 292 (2016), 576–590. https://doi.org/10.1016/j.cam.2015.07.039 doi: 10.1016/j.cam.2015.07.039
|
| [15] |
H.G. Lee, Stability condition of the second-order SSP-IMEX-RK method for the Cahn–Hilliard equation, Mathematics, 8 (2020), 11. https://doi.org/10.3390/math8010011 doi: 10.3390/math8010011
|
Figures(8)
Hyun Geun Lee. A second-order linear energy-stable scheme for slope-selection epitaxial thin-film growth[J]. AIMS Mathematics, 2026, 11(4): 9989-10003. doi: 10.3934/math.2026412
DownLoad: