We investigate the two-pion interferometry in ultrarelativistic heavy ion collisions in the granular source model of quark-gluon plasma droplets. The pion transverse momentum spectra and HBT radii of the granular sources agree well with the experimental data of the √ sNN = 200 GeV Au-Au and √ sNN = 2.76 TeV Pb-Pb most central collisions. In the granular source model the larger initial system breakup time may lead to the larger HBT radii Rout, R side , and R long . However, the large droplet transverse expansion and limited average relative emitting time of particles in the granular source lead to small ratios of the transverse HBT radii Rout/R side .PACS numbers: 25.75.Gz Hanbury-Brown-Twiss (HBT) interferometry is a useful tool to probe the space-time geometry of the particleemitting sources in high energy heavy ion collisions [1][2][3][4]. The experimental results of the HBT measurements for the Au-Au collisions at the high energies of the Relativistic Heavy Ion Collider (RHIC) indicate that it is hard to describe the source space-time by a simple evolution model [5][6][7][8]. HBT interferometry data provide strong constraints for the models of source space-time. Recently, the HBT measurement for the √ s N N = 2.76TeV Pb-Pb most central collisions at the Large Hadron Collider (LHC) is performed [9]. A consistent explanation to the HBT data of the LHC and RHIC experiments is required naturally for the source models, which will be helpful to understand the initial condition, source evolution, and particle freeze-out in ultrarelativistic heavy ion collisions. In Refs. [10,11], the granular source model of quarkgluon plasma (QGP) droplets [12] is developed to explain the RHIC HBT data [7,8]. In this work we investigate the two-pion HBT interferometry in ultrarelativistic heavy ion collisions in the granular source model of QGP droplets. Our results indicate that the granular source for the LHC Pb-Pb collisions may have the same initial droplet temperature and velocity formula as those for the RHIC Au-Au collisions, but a larger initial system breakup time. The consistent granular source model reproduces the pion transverse momentum spectra and HBT radii in the most central collisions of the RHIC [8,13,14] and LHC [9, 15] experiments.In ultrarelativistic heavy ion collisions, the system at central rapidity may reach a local equilibrium at a very short time τ 0 , then fast expand in the beam direction (z-axis). Because of the initial fluctuation, the local- * wnzhang@dlut.edu.cn equilibrium system is not uniform in space [16,17]. It may form many tubes along the beam direction during the fast longitudinal expansion, and finally fragment into many droplets (see Fig.1 of Ref.[10]) due to the "sausage" instability and surface tension [10]. On the other hand, the rapidly increased bulk viscosity in the QGP near the phase transition may also leads to the system breakup [18].We assume that the system fragments and forms a granular source of many QGP droplets at a time t 0 (> τ 0 ). On the basis of the Bjorken hypothe...