The present paper performed a numerical study on two-phase nanofluid (NFs) flow in a circular heatsink for cooling several LEDs. The heatsink is symmetrically designed and has two inlets and four outlets. Six heat sources or LEDs are placed on the circumference of a circle and a heat source is also mounted in the center of the heatsink. By varying the diameter of the circle, the side length of the heat sources, and the input velocity of the NFs, one may estimate the values of thermal resistance (THR), temperature uniformity (TUY) on the heatsink, heat transfer coefficient (HTC), and pressure drop in the heatsink. The finite element and two-phase mixture method are utilized for NFs simulations. It demonstrate that the heat source placed in the middle has a lower temperature than other heat sources. The results are most significantly affected by changing the NFs' velocity. The value of dimensionless temperature increases and subsequently decreases as the sides of the heat sources get longer. The dimensionless temperature first decreases and then increases as the distance between the heat sources and the heatsink's center increases. The amount of THR is high when the heat sources' side length or velocity values are large.