Abstract
The numerical simulation of the dynamics of the molecular beam epitaxy (MBE) growth is considered in this article. The governing equation is a nonlinear evolutionary equation that is of linear fourth order derivative term and nonlinear second order derivative term in space. The main purpose of this work is to construct and analyze two linearized finite difference schemes for solving the MBE model. The linearized backward Euler difference scheme and the linearized Crank-Nicolson difference scheme are derived. The unique solvability, unconditional stability and convergence are proved. The linearized Euler scheme is convergent with the convergence order of O(τ + h 2) and linearized Crank-Nicolson scheme is convergent with the convergence order of O(τ 2 + h 2) in discrete L 2-norm, respectively. Numerical stability with respect to the initial conditions is also obtained for both schemes. Numerical experiments are carried out to demonstrate the theoretical analysis.
Original language | English |
---|---|
Pages (from-to) | 1893-1915 |
Number of pages | 23 |
Journal | Numerical Methods for Partial Differential Equations |
Volume | 28 |
Issue number | 6 |
DOIs | |
Publication status | Published - 1 Nov 2012 |
Externally published | Yes |
Keywords
- convergence
- finite difference scheme
- linearized difference scheme
- molecular beam epitaxy
- stability
ASJC Scopus subject areas
- Analysis
- Numerical Analysis
- Computational Mathematics
- Applied Mathematics