The present study investigated the steady magnetohydrodynamics of the axisymmetric flow of a incompressible, viscous, electricity-conducting nanofluid with convective boundary conditions and thermo-diffusion over a radially stretched surface. The nanoparticles’ volume fraction was passively controlled on the boundary, rather than actively controlled. The governing non-linear partial differential equations were transformed into a system of nonlinear, ordinary differential equations with the aid of similarity transformations which were solved numerically, using the very efficient variational finite element method. The coefficient of skin friction and rate of heat transfer, and an exact solution of fluid flow velocity, were contrasted with the numerical solution gotten by FEM. Excellent agreement between the numerical and exact solutions was observed. The influences of various physical parameters on the velocity, temperature, and solutal and nanoparticle concentration profiles are discussed by the aid of graphs and tables. Additionally, authentication of the convergence of the numerical consequences acquired by the finite element method and the computations was acquired by decreasing the mesh level. This exploration is significant for the higher temperature of nanomaterial privileging technology.