The present study aims to provide a mathematical model of the Williamson fluid flow via a permeable stretching/shrinking sheet in the MHD boundary layer in the presence of a heat source, chemical reaction, and suction. This study is novel because it investigates the physical effects of thermal and solutal stratification on convective heat and mass transport using thermal radiation. The flow’s PDEs are numerically solved using the BVP4c approach and the pertinent similarity variables until a stable solution is found. Through visual analysis, the effects of dimensionless factors on temperature, velocity, and concentration profiles are examined. This encompasses the mass transfer rate, the heat transfer rate, and the coefficient of friction. The results of the present analysis are found to be consistent with those of previously published studies. The findings demonstrate that enhanced temperature and concentration profiles cause the Williamson, magnetic, and permeability parameters to rise in conjunction with a drop in the dimensionless velocity. In relation to temperature, the thermal stratification parameter exhibits the opposite tendency. Regarding the solutal stratification parameter, concentration profiles are seen to show the opposite trend. Lastly, the current work will have important implications for the removal of dust and viruses from viscoelastic fluid in bioengineering, the medical sciences, and medical equipment.