This study presents a mathematical model that delineates the radially expanding axisymmetric discharge of an electrically conductive fluid over a surface, taking into account the effects of the Soret number. The dynamics of the flow are examined as the surface experiences exponential radial expansion. To transform the governing nonlinear partial differential equations into standard derivative forms, similarity transformations are applied. The flow dynamics are further investigated using the Successive Linearization Method. To achieve accurate solutions that converge effectively to the complete numerical solution, the Chebyshev spectral method is employed to solve the resulting linear system. Previous research is cited to support the findings related to the distribution of velocity, temperature, and concentration, emphasizing the convergence and accuracy of the solution while considering the influence of various fluid parameters.