This paper analytically solves the Newtonian fluid movement over a semi-infinite porous expanding/contracting boundary with multiple slips in the company of carbon nanotubes (CNTs), thermal radiation, and a magnetic field. There are two distinct kinds of CNT models used, namely single-walled CNTs (SWCNT) and multiple-walled CNTs (MWCNT), with the possibility of fluid sliding in accordance with Navier’s and second-order slips. The dimensionless and similarity transformations with the radiative heat flux approximation are applied to the transformation of governing equations into ordinary differential equations (ODEs). Since the flow field solution contains a shift exponential function, the balanced heat equation between temperature and thermal radiation leads to a linear ODE with exponential-type coefficients. The result of the energy equation is then expressed as a function of incomplete gamma functions using the Laplace transform. Significant performances of SWCNTs and MWCNTs to the fluid flow and heat transfer are evaluated by varying physical parameters such as Prandtl number, Brinkman ratio, viscosity ratio, thermal radiation, and multiple slip parameters. These findings have industrial applications, such as polymer manufacture and oil extraction, as well as medicinal applications.