METHODS: The governing equations of fluid flow have been transformed in the form of ordinary differential equations. These equations have been solved by two methods namely, shooting method and three-stage Lobatto IIIa formula.
RESULTS: The effects of different parameters on temperature, velocity, concentration profiles, skin friction coefficient, Sherwood number, and reduced Nusselt number were obtained and presented graphically. It was noticed that four solutions existed at definite ranges of the parameters for high suction over both surfaces for the first time. The results of the stability analysis revealed that only the first solution is more stable and possess physical reliability compared to the remaining solutions.
CONCLUSION: The graphs also indicated that the fluid velocity decreases as the thermophoresis parameter increases but the opposite behavior observed for both temperature and concentration profiles in the first solution. Furthermore, it was detected that the concentration profile declined at the higher values of the Brownian motion parameter.