Consistency,stability and convergence of the finite-difference equations for flow about a rotating sphere in an axial stream

The finite-difference equations which have previously been developed to solve the problem of laminar boundary layer flow about a rotating sphere in an axial stream are analysed according to the available numerical stability theories. This analysis is necessary to determine the restrictions on velocities and mesh sizes required to obtain a convergent numerical solution. Convergence can be achieved if both consistency and stability of the finite-difference equations are fulfilled. The analysis reported in the present paper shows that the developed finite-difference equations are consistent with their original partial differential equations. Also, the analysis proves that the developed finite-difference procedure is numerically stable for all mesh sizes as long as the downstream meridional velocity is non-negative, i.e.as long as no flow reversals occur within the domain of solution.