Abstract A theoretical model is developed for steady magnetohydrodynamic viscous flow resulting from a moving semi-infinite flat plate in an electrically conducting nanofluid. Thermal radiation and magnetic induction effects are included in addition to thermal convective boundary conditions. Buongiorno’s two-component nanoscale model is deployed, which features Brownian motion and thermophoresis effects. The governing nonlinear boundary layer equations are converted to nonlinear ordinary differential equations by using suitable similarity transformations. The transformed system of differential equations is solved numerically, employing the spectral relaxation method (SRM) via the MATLAB R2018a software. SRM is a simple iteration scheme that does not require any evaluation of derivatives, perturbation, and linearization for solving a nonlinear system of equations. Effects of embedded parameters such as sheet velocity parameter$\lambda$, magnetic field parameter$\beta$, Prandtl number$Pr$, magnetic Prandtl number$Prm$, thermal radiation parameter$Rd$, Lewis number$Le$, Brownian motion parameter$Nb$, and thermophoresis parameter$Nt$ on velocity, induced magnetic field, temperature, and nanoparticle concentration profiles are investigated. The skin-friction results, local Nusselt number, and Sherwood number are also discussed for various values of governing physical parameters. To show the convergence rate against iteration, residual error analysis has also been performed. The flow is strongly decelerated, and magnetic induction is suppressed with greater magnetic body force parameter, whereas temperature is elevated due to extra work expended as heat in dragging the magnetic nanofluid. Temperatures are also boosted with increment in nanoscale thermophoresis parameter and radiative parameter, whereas they are reduced with higher wall velocity, Brownian motion, and Prandtl numbers. Both hydrodynamic and magnetic boundary layer thicknesses are reduced with greater reciprocal values of the magnetic Prandtl number Prm. Nanoparticle (concentration) boundary layer thickness is boosted with higher values of thermophoresis and Prandtl number, whereas it is diminished with increasing wall velocity, nanoscale Brownian motion parameter, radiative parameter, and Lewis number. The simulations are relevant to electroconductive nanomaterial processing.