Engineers and researchers in the field of heat transfer are in search of new techniques to optimize the performance of energy devices through heat transfer enhancement. For that reason, the flow analysis involving nanoparticles is assumed to be one of the most important techniques for enhancing heat transfer systems. Being a two-component system, different numerical approaches are available to model the thermo-fluids behavior of nanofluids. The current study describes the investigation of Buongiorno model for evaluation of transient flow and heat transfer of Williamson nanofluids under the impacts of velocity slip and magnetic field. In addition, the effects of heat generation/absorption and convective heat transfer have been employed in this analysis. The Boussinesq-approximations are implemented to obtain the governing conservation equations for nanofluids transport phenomenon. The leading equations of the modeled physical problem have been solved using a MATLAB code, based on Nachtsheim–Swigert shooting iteration scheme. The novelty of the current investigation is the existence of multiple solutions for the nanofluids flow past a radially shrinking surface by considering the effects of uniform suction at the wall. In this review, it is seen that the investigated physical parameters affect remarkably on the nanofluid stream function, temperature and nanoparticles concentration. The skin friction, local Nusselt number and the local Sherwood number and some samples of velocity, temperature and nanoparticle concentration profiles are discussed and presented with the help of graphs. From this study, it is observed that the existence range of dual solutions is significantly raised by the higher Weissenberg number.