Accelerated convergence in numerical simulations of surface supersaturation for crystal growth in solution under steady-state conditions

This is an investigative paper which reports the results of comparisons of two numerical techniques for the solution of the Burton Cabrera and Frank (BCF) equation for the growth on crystal surfaces under steady state conditions. A successive over-relaxation (SOR) scheme for the equivalent finite difference equation gives rapid convergence to the static solution. It is known that a suitable choice of scattering parameters in a transmission line matrix (TLM) network analogue of the Laplace equation yields ultra-fast convergence. The results of numerical experiments which are reported here suggests that a similar situation also applies to the solution of the Poisson equation with shunt losses (the BCF equation), although the choice of optimum conditions appears to be different for different spatial positions within the solution space.