This article is concerned with the fluid mechanics of MHD steady 2D flow of Williamson fluid over a nonlinear stretching curved surface in conjunction with homogeneous-heterogeneous reactions with convective boundary conditions. An effective similarity transformation is considered that switches the nonlinear partial differential equations riveted to ordinary differential equations. The governing nonlinear coupled differential equations are solved by using MATLAB bvp4c code. The physical features of nondimensional Williamson fluid parameter
λ
, power-law stretching index
m
, curvature parameter
K
, Schmidt number
Sc
, magnetic field parameter
M
, Prandtl number
Pr
, homogeneous reaction strength
k
1
, heterogeneous reaction strength
k
2
, and Biot number
γ
are presented through the graphs. The tabulated form of results is obtained for the skin friction coefficient. It is noted that both the homogeneous and heterogeneous reaction strengths reduced the concentration profile.