An initial performance of solid oxide fuel cell deeply correlates with a composite sintering process during anode fabrication. In this study, we develop a three-dimensional kinetic Monte Carlo (KMC) code based on a q-state Potts model in order to numerically predict a microstructural evolution during NiO-YSZ composite sintering. Our simulations successfully predict the grain size-density trajectories and the tortuosity factor evolution during sintering that well match with our experimental measurements using a three-dimensional FIB-SEM tomography.