In this work, steady, magnetized 2D incompressible Jeffrey fluid flow on a stretched curved surface is modeled, considering the ferric oxide \(\left(F{e}_{2}{O}_{3}\right)\) , graphene oxide \(\left(GO\right)\) , and zirconium dioxide \(\left(Zr{O}_{2}\right)\) , immersed in the base fluid ethylene glycol (EG). The model includes, among other things, viscous dissipation, mixed convection, magnetohydrodynamics, and Joule heating. This work operates the motion properties of a non-Newtonian fluid by applying a magnetic force perpendicular to the flow. This mathematical formulation is based on boundary layer approximations and the assumption of a low magnetic Reynolds number. Nonlinear partial differential equations are created from the governing system by applying the proper non-similar transformations. The local non-similarity approach is utilized for simulation up to the second truncation level using the MATLAB built-in solver bvp4c. The temperature and velocity profiles of nanofluids are shown in a variety of graphs and numerical tables. In addition, the skin friction coefficient and Nusselt number are calculated numerically, taking into consideration the effects of the Deborah number, slip parameter, Eckert number, and Prandtl number. A greater Eckert number causes the Nusselt number to rise, whereas a higher thermophoretic parameter causes it to fall. Improving heat transfer efficiency in industrial heat exchangers, cooling systems, chemical processing, biomedical engineering, aerospace engineering, renewable energy systems, the automotive sector, and manufacturing processes are just a few of the practical uses for the present study.