The city scale buoyancy-driven flow is crucial for the wind and thermal environment in urban areas, exerting a significant impact on the city ventilation, pollutants dispersion and heat removal. The intra-urban variations of hot spots and heat flux are drawing more attention on beating the urban heat. This study aims at analyzing the influences of non-uniform heat flux and background wind on the city-scale heat removal ability and underlying mechanisms. The large eddy simulation model was used and verified with water tank experiments. Results show that positive heat flux and negative heat flux in the rural areas will enhance and depress the strength of the urban heat dome flow respectively compared to the adiabatic condition. The high-speed regions are wider and stronger when the rural area has positive heat flux condition. Under positive rural heat flux condition, the time-average mixed layer height in urban areas across the city center plane is nearly twice that under negative rural heat flux condition. Besides, the city-scale heat transfer coefficients in positive rural heat flux condition are higher than those under negative rural heat flux and adiabatic conditions. For cases of non-uniform heat flux in urban areas, the differences of time-average mixed layer height and heat transfer coefficients are not obvious. The heat transfer coefficients decrease initially and then increase with increasing background wind speed indicating the non-linear interaction between approaching wind and urban heat dome flow. Moreover, the growth rate of the heat transfer coefficient decreases when the non-dimensional wind speed exceeds 1.19.