Analysis of dual solution for MHD flow of Williamson fluid with slippage

This study investigates the numerical solutions of MHD boundary layer and heat transfer of the Williamson fluid flow on the exponentially vertical shrinking sheet, having variable thickness and thermal conductivity under effects of the velocity and thermal slip parameters. It is also assumed that shrinking/stretching velocity, as well as the wall temperature, has the exponential function form. In this study, the continuity, momentum and energy equations with buoyancy parameter and Hartmann number are incorporated especially in the Williamson fluid flow case. Similarity transformation variables have been employed to formulate the ordinary differential equations (ODEs) from partial differential equations (PDEs). The resultant ODEs are solved by shooting method with Runge Kutta of fourth order method in Maple software. The effects of the different applied non-dimensional physical parameters on the boundary layer and heat transfer flow problems are presented in graphs. The effects of Williamson parameter, Prandtl number, and slip parameters on velocity and temperature profiles have been thoroughly demonstrated and discussed. The numerical results show that the buoyancy force and the slip parameters contribute to the occurrence of the dual solutions on the boundary layer and heat transfer flow problems. Furthermore, the stability analysis suggests that the first solution is stable and physically possible.