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.

A stabilized second order exponential time differencing multistep method for thin film growth model without slope selection

Energy stability and error estimates of exponential time differencing schemes for the epitaxial growth model without slope selection

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

A Second-Order Energy Stable Backward Differentiation Formula Method for the Epitaxial Thin Film Equation with Slope Selection

A Second-Order Exponential Time Differencing Multi-step Energy Stable Scheme for Swift–Hohenberg Equation with Quadratic–Cubic Nonlinear Term

Finite Speed of Propagation and Algebraic Time Decay of Solutions to A Generalized Thin Film Equation

Global-in-time energy stability analysis for the exponential time differencing Runge-Kutta scheme for the phase field crystal equation

Structure-preserving, energy stable numerical schemes for a liquid thin film coarsening model

New Epitaxial Thin‐film Models and Numerical Approximation

An Exponential Time Differencing Runge–Kutta Method ETDRK32 for Phase Field Models

An Accurate and Efficient Algorithm for The Time-fractional Molecular Beam Epitaxy Model with Slope Selection

Analysis of the second order BDF scheme with variable steps for the molecular beam epitaxial model without slope selection

A fourth-order exponential time differencing scheme with dimensional splitting for non-linear reaction-diffusion systems

Thin film deposition with time varying temperature

A class of fourth-order nonlinear parabolic equations modeling the epitaxial growth of thin films

An efficient second order stabilized scheme for the two dimensional time fractional Allen-Cahn equation

A Second-Order, Global-in-Time Energy Stable Implicit-Explicit Runge–Kutta Scheme for the Phase Field Crystal Equation

Convergence analysis of exponential time differencing scheme for the nonlocal Cahn-Hilliard equation

The operator splitting method with semi-implicit spectral deferred correction for molecular beam epitaxial growth models

A stabilized semi-implicit Fourier spectral method for nonlinear space-fractional reaction-diffusion equations

Asymptotic Behaviour of Time Stepping Methods for Phase Field Models