The aim of this research is mainly concerned with the numerical examination of Darcy-Forchheimer relation in convective magnetohydrodynamic nanofluid flow bounded by non-linear stretching sheet. A visco-elastic and strictly incompressible liquid saturates the designated porous medium under the direct influence of the Darcy-Forchheimer model and convective boundary. The magnetic effect is taken uniformly normal to the flow direction. However, the model is bounded to a tiny magnetic Reynolds number for practical applications. Boundary layer formulations are taken into consideration. The so-formulated leading problems are converted into highly nonlinear ordinary problems using effectively modified transformations. The numerical scheme is applied to solve the governing problems. The outcomes stipulate that thermal layer receives significant modification in the incremental direction for augmented values of thermal radiation parameter Rd. Elevation in thermal Biot number γ1 apparently results a significant rise in thermal layer and associated boundary layer thickness. The solute Biot number is found to be an enhancing factor the concentration profile. Besides the three main profiles, the contour and density graphs are sketched for both the linear and non-linear cases. Furthermore, skin friction jumps for larger porosity and larger Forchheimer number. Both the heat and mass flux numbers receive a reduction for augmented values of the Forchheimer number. Heat flux enhances, while mass flux reduces, the strong effect of thermal Biot number. The considered problem could be helpful in any several industrial and engineering procedures, such as rolling, polymeric extrusion, continuously stretching done in plastic thin films, crystal growth, fiber production, and metallic extrusion, etc.