In this study, the Spectral Relaxation Method (SRM) is used to solve the coupled highly nonlinear system of partial differential equations due to an unsteady flow over a stretching surface in an incompressible rotating viscous fluid in presence of binary chemical reaction and Arrhenius activation energy. The velocity, temperature and concentration distributions as well as the skin-friction, heat and mass transfer coefficients have been obtained and discussed for various physical parametric values. The numerical results obtained by (SRM) are then presented graphically and discussed to highlight the physical implications of the simulations.

Nonlinear partial differential equations (PDEs) modelling unsteady boundary-layer flows are solved by the spectral relaxation method (SRM) and the spectral quasilinearization method (SQLM). The SRM and SQLM are Chebyshev pseudospectral based methods that have been successfully used to solve nonlinear boundary layer flow problems described by systems of ordinary differential equations. In this paper application of these methods is extended, for the first time, to systems of nonlinear PDEs that model unsteady boundary layer flow. The new extension is tested on two problems: boundary layer flow caused by an impulsively stretching plate and a coupled four-equation system that models the problem of unsteady MHD flow and mass transfer in a porous space. Numerous simulation experiments are conducted to determine the accuracy and compare the computational performance of the proposed methods against the popular Keller-box finite difference scheme which is widely accepted as being one of the ideal tools for solving nonlinear PDEs that model boundary layer flow problems. The results indicate that the methods are more efficient in terms of computational accuracy and speed compared with the Keller-box.

There are currently only a few theoretical studies on convective heat transfer in polymer nanocomposites. In this paper, the unsteady incompressible flow of a polymer nanocomposite represented by an Oldroyd-B nanofluid along a stretching sheet is investigated. Recent studies have assumed that the nanoparticle fraction can be actively controlled on the boundary, similar to the temperature. However, in practice, such control presents significant challenges and in this study the nanoparticle flux at the boundary surface is assumed to be zero. We have used a relatively novel numerical scheme; the spectral relaxation method to solve the momentum, heat and mass transport equations. The accuracy of the solutions has been determined by benchmarking the results against the quasilinearisation method. We have conducted a parametric study to determine the influence of the fluid parameters on the heat and mass transfer coefficients.

A new multistage numerical method based on blending a Gauss‐Siedel relaxation method and Chebyshev pseudospectral method, for solving complex dynamical systems exhibiting hyperchaotic behavior, is presented. The proposed method, called the multistage spectral relaxation method (MSRM), is applied for the numerical solution of three hyperchaotic systems, namely, the Chua, Chen, and Rabinovich‐Fabrikant systems. To demonstrate the performance of the method, results are presented in tables and diagrams and compared to results obtained using a Runge‐Kutta‐(4,5)‐based MATLAB solver, ode45, and other previously published results.

Li et al. (SIAM J. Sci. Comput. 20:719–738, 1998) used the moving mesh partial differential equation (MMPDE) to solve a scaled Fisher’s equation and the initial condition consisting of an exponential function. The results obtained are not accurate because MMPDE is based on a familiar arc-length or curvature monitor function. Qiu and Sloan (J. Comput. Phys. 146:726–746, 1998) constructed a suitable monitor function called modified monitor function and used it with the moving mesh differential algebraic equation (MMDAE) method to solve the same problem of scaled Fisher’s equation and obtained better results. In this work, we use the forward in time central space (FTCS) scheme and the nonstandard finite difference (NSFD) scheme, and we find that the temporal step size must be very small to obtain accurate results. This causes the computational time to be long if the domain is large. We use two techniques to modify these two schemes either by introducing artificial viscosity or using the approach of Ruxun et al. (Int. J. Numer. Methods Fluids 31:523–533, 1999). These techniques are efficient and give accurate results with a larger temporal step size. We prove that these four methods are consistent for partial differential equations, and we also obtain the region of stability.