We present a linear numerical scheme for a model of epitaxial thin film growth without slope selection. the PDE, which is a nonlinear, fourth-order parabolic equation, is the L2 gradient flow of the energy ∫Ω(-1/2 ln(1 + |ø|2) + ε2 2 |Ø(x)|2) dx. the idea of convex-concave decomposition of the energy functional is applied, which results in a numerical scheme that is unconditionally energy stable, i.e., energy dissipative. the particular decomposition used here places the nonlinear term in the concave part of the energy, in contrast to a previous convexity splitting scheme. as a result, the numerical scheme is fully linear at each time step and unconditionally solvable. Collocation Fourier spectral differentiation is used in the spatial discretization, and the unconditional energy stability is established in the fully discrete setting using a detailed energy estimate. We present numerical simulation results for a sequence of values ranging from 0.02 to 0.1. in particular, the long time simulations show the -log(t) decay law for the energy and the t 1/2 growth law for the surface roughness, in agreement with theoretical analysis and experimental/numerical observations in earlier works. © Springer Science+Business Media, LLC 2011.
W. Chen et al., "A Linear Energy Stable Scheme for a Thin Film Model Without Slope Selection," Journal of Scientific Computing, vol. 52, no. 3, pp. 546 - 562, Springer, Jan 2012.
The definitive version is available at https://doi.org/10.1007/s10915-011-9559-2
Mathematics and Statistics
Keywords and Phrases
Convexity Splitting; Energy Stability; Epitaxial Thin Film Growth; Fourier Collocation Spectral; Slope Selection
International Standard Serial Number (ISSN)
Article - Journal
© 2023 Springer, All rights reserved.
01 Jan 2012