A Second Order Energy Stable Linear Scheme for a Thin Film Model Without Slope Selection

Weijia Li,Wenbin Chen,Cheng Wang,Yue Yan,Ruijian He

DOI: https://doi.org/10.1007/s10915-018-0693-y

2018-03-10

Journal of Scientific Computing

Abstract:In this paper we present a second order accurate, energy stable numerical scheme for the epitaxial thin film model without slope selection, with a mixed finite element approximation in space. In particular, an explicit treatment of the nonlinear term, ∇u1+|∇u|2\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\frac{\nabla u}{1+|\nabla u|^2}$$\end{document}, greatly simplifies the computational effort; only one linear equation with constant coefficients needs to be solved at each time step. Meanwhile, a second order Douglas–Dupont regularization term, AτΔ2(un+1-un)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$A\tau \varDelta ^2 ( u^{n+1} - u^n)$$\end{document}, is added in the numerical scheme, so that an unconditional long time energy stability is assured. In turn, we perform an ℓ∞(0,T;L2)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\ell ^\infty (0,T; L^2)$$\end{document} convergence analysis for the proposed scheme, with an O(τ2+hq)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$O (\tau ^2 + h^q)$$\end{document} error estimate derived. In addition, an optimal convergence analysis is provided for the nonlinear term using Qq\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$Q_q$$\end{document} finite elements, which shows that the spatial convergence order can be improved to q+1\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$q+1$$\end{document} on regular rectangular mesh. A few numerical experiments are presented, which confirms the efficiency and accuracy of the proposed second order numerical scheme.

mathematics, applied