A two-dimensional cellular automata model was developed to describe the static recrystallization (SRX) arising from the subgrain growth, the driving force of which is dependent on boundary energy and local curvature. At the same time, the subgrain boundary energy and mobility rely on the boundary misorientation angle. On the basis, a deterministic switch rule was adopted to simulate the subgrain growth and kinetics of recrystallization quantitatively to provide an insight into the grain boundary bulging nucleation mechanism. Microstructure evolutions during SRX in different cases were simulated by the developed model. At the beginning of the simulation, the initial polycrystalline microstructure which contains large number of uniformly distributed subgrains in every pre-existing grain was prepared using simple assumption based on experimental observations. Then, both homogeneous and inhomogeneous subgrain growth phenomena were captured by the simulation with different inter-subgrain misorientation, which showed continuous and discontinuous recrystallization, respectively. The effects of initial mean subgrain radius, distribution of initial subgrains, distribution of inter-subgrain misorientations, and annealing temperature on the recrystallization kinetics were also investigated.