The current article aims to present a numerical analysis of MHD Williamson nanofluid flow maintained to flow through porous medium bounded by a non-linearly stretching flat surface. The second law of thermodynamics was applied to analyze the fluid flow, heat and mass transport as well as the aspects of entropy generation using Buongiorno model. Thermophoresis and Brownian diffusion is considered which appears due to the concentration and random motion of nanoparticles in base fluid, respectively. Uniform magnetic effect is induced but the assumption of tiny magnetic Reynolds number results in zero magnetic induction. The governing equations (PDEs) are transformed into ordinary differential equations (ODEs) using appropriately adjusted transformations. The numerical method is used for solving the so-formulated highly nonlinear problem. The graphical presentation of results highlights that the heat flux receives enhancement for augmented Brownian diffusion. The Bejan number is found to be increasing with a larger Weissenberg number. The tabulated results for skin-friction, Nusselt number and Sherwood number are given. A decent agreement is noted in the results when compared with previously published literature on Williamson nanofluids.