The present investigation is devoted to analyze the buoyancy-driven flow behavior and associated thermal dissipation rate in a nanofluid-filled annular region with five different single source-sink and three different dual source-sink arrangements along the vertical surfaces. The remaining region on the vertical boundaries and horizontal surfaces are kept adiabatic. Numerical simulations have been performed by solving the dimensionless model equations by employing finite difference method. To analyze the impacts of the type of nanofluid, nanoparticle volume fraction, Rayleigh number, size and arrangement of sources and sinks, the results are graphically represented through streamline and isotherm contours, thermal profiles, average Nusselt number and cup-mixing temperature. Comparison of the results of current investigation reveal fairly good agreement with the existing benchmark results. The results showed that identifying an optimum location and length of source-sink with a proper selection of other control parameters can lead to enhanced thermal transport and thermal mixing in the enclosure. Also, the calculated enhancement ratio of the heat dissipation rate enhances with an increment in nanoparticle concentration.