A numerical investigation of the 3D unsteady flow and heat transfer performance induced by multi-piezoelectric fans on cooling asymmetrical concave surfaces is performed using dynamic meshing scheme according to its vibration fitting function. The effects of the relative curvature （K_r） and vibration phase difference （φ） on heat removal performance are described. The results show that the asymmetrical nature of the time-averaged convective heat transfer coefficient distributions is mainly due to a combination of the interaction of neighboring fans in the spanwise direction and the asymmetry of the surfaces in the chordwise direction. In particular， when K_r is decreased to 2， the asymmetry of distribution almost entirely disappears， and the heat transfer coefficient in the vibration envelopes is almost twice that of K_r=6. Since the characteristic scale of the shedding vortex excited by multi-fan vibrating in-phase （φ=0°， the phase difference between adjacent fans is 0°） is larger than that of out-of-phase （φ=180°， the phase difference between adjacent fans is 180°）， that the heat transfer capacity in vibration envelope under in-phase （φ=0°） is significantly better than that under out-of-phase （φ=180°）. In addition， the gap vortices generated by the in-phase vibration （φ=0°） also increase the heat transfer capactity in the adjacent fan gap.