This article examines the natural convection Couette flow in a hydrodynamic viscous fluid that is electrically conductive due to the thermal radiation effect. The governing flow in this study is modelled in the form of partial differential equations (PDEs) in dimensional form with initial and boundary conditions and the Couette fluid model is also be used to characterize the fluid behavior. Then, using suitable non-dimensional quantities, the governing non-linear PDEs are transformed. Since the flow governing equations of the problem under study are extremely complex and complicated, techniques that complement experimental and theoretical fluid dynamics by providing alternative potentially cheaper means of testing fluid flow systems is employed. Therefore, the Finite Element Method (FEM) is used after discretization of the PDEs. With Graphs and tables, the effects of embedded thermo physical parameters of engineering interests associated with the flow quantities viz. velocity, temperature, concentration of the fluid were examined through series of numerical experiments and discussed. This research also analyzes and compares the results obtained by Omokhuale and Jabaka (2022b). It is interesting to report that an excellent agreement was established, thereby authenticating and validating the accuracy of FEM as a strong tool. According to the results of this study, the actions of thermal radiation on the thermal and momentum boundary layers for increasing values are significant, also, increasing the magnetic field parameter impedes the fluid movement due to the Lorentz force action.