This numerical investigation examines the behavior of a quiescent couple-stress fluid and the associated heat transfer phenomena on a stretching/shrinking cylinder subjected to specified surface temperature and heat flux boundary conditions. It is founded on the similarity solution that highly non-linear partial differential equations (PDEs) can be simplified to highly non-linear ordinary differential equations (ODEs). The problem involves several parameters, including couple stress, suction/injection, magnetic field inclination, stretching/shrinking, heat source/sink, radiation parameters, as well as the Hartmann, Darcy, Prandtl, and Eckert numbers. The reduced ODEs were numerically simulated utilizing Maple. The temperature profiles and local Nusselt numbers for various boundary conditions are compared using plots and tables. The results indicated that the impact of prescribed heat flux on reducing the local Nusselt number is significantly greater than that of prescribed surface temperature for a given Eckert number. Also, under PST, the local Nusselt number rises more with s than under PHF.



