The rheological properties of the non-Newtonian fluids make them highly beneficial for various industrial and engineering applications because their unique characteristics offer advantages in a wide range of fields. A study on the behavior of Jeffrey fluid with viscous dissipation on a bi-directional stretching surface is conducted. Homogeneous and heterogeneous reactions in electrically conducting Jeffrey fluid are considered. Governing equations and boundary conditions are derived and transformed into nondimensional forms using similarity transformations. MATLAB's bvp4c solver is employed to solve and analyze these equations. Graphical representations and tables investigate the impact of various flow parameters, including temperature, concentration, velocity, skin friction, heat, and mass transfer rates. Results indicate that with an increase in the Prandtl number, temperature field decline. The velocity, thermal, and solutal profiles drop when the Deborah number increases. Moreover, comparing homogeneous and heterogeneous parameters reveals a decrease in the surface's heat and mass transfer rate in the presence of thermal transport and concentration profile variations at different values of dimensionless variables.