On the Multidomain Bivariate Spectral Local Linearisation Method for Solving Systems of Nonsimilar Boundary Layer Partial Differential Equations

In this work, a novel approach for solving systems of nonsimilar boundary layer equations over a large time domain is presented. The method is a multidomain bivariate spectral local linearisation method (MD-BSLLM), Legendre-Gauss-Lobatto grid points, a local linearisation technique, and the spectral collocation method to approximate functions defined as bivariate Lagrange interpolation. The method is developed for a general system of nonlinear partial differential equations. The use of the MD-BSLLM is demonstrated by solving a system of nonlinear partial differential equations that describe a class of nonsimilar boundary layer equations. Numerical experiments are conducted to show applicability and accuracy of the method. Grid independence tests establish the accuracy, convergence, and validity of the method. The solution for the limiting case is used to validate the accuracy of the MD-BSLLM. The proposed numerical method performs better than some existing numerical methods for solving a class of nonsimilar boundary layer equations over large time domains since it converges faster and uses few grid points to achieve accurate results. The proposed method uses minimal computation time and its accuracy does not deteriorate with an increase in time.