The nanofluid and porous medium together are able to fulfill the requirement of high cooling rate in many engineering problems. So, here the impact of various shapes of nanoparticles on unsteady stagnation-point flow of Cu–H2O nanofluid on a flat surface in a porous medium is examined. Moreover, the thermal radiation and viscous dissipation effects are considered. The problem governing partial differential equations are converted into self-similar coupled ordinary differential equations and those are numerically solved by the shooting method. The computed results can reveal many vital findings of practical importance. Firstly, dual solutions exist for decelerating unsteady flow and for accelerating unsteady and steady flows, the solution is unique. The presence of nanoparticles affects the existence of dual solution in decelerating unsteady flow only when the medium of the flow is a porous medium. But different shapes of nanoparticles are not disturbing the dual solution existence range, though it has a considerable impact on thermal conductivity of the mixture. Different shapes of nanoparticles act differently to enhance the heat transfer characteristics of the base fluid, i.e., the water here. On the other hand, the existence range of dual solutions becomes wider for a larger permeability parameter related to the porous medium. Regarding the cooling rate of the heated surface, it rises with the permeability parameter, shape factor (related to various shapes of Cu-nanoparticles), and radiation parameter. The surface drag force becomes stronger with the permeability parameter. Also, with growing values of nanoparticle volume fraction, the boundary layer thickness (BLT) increases and the thermal BLT becomes thicker with larger values of shape factor. For decelerating unsteady flow, the nanofluid velocity rises with permeability parameter in the case of upper branch solution and an opposite trend for the lower branch is witnessed. The thermal BLT is thicker with radiation parameter. Due to the existence of dual solutions, a linear stability analysis is made and it is concluded that the upper branch and unique solutions are stable solutions.