A Monte Carlo model for simulating final-stage sintering has been developed. This model incorporates realistic microstructural features (grains and pores), variable surface difusivity, grain-boundary diffusivity, and grain-boundary mobility. A preliminary study of a periodic array of pores has shown that the simulation procedure accurately reproduces theoretically predicted sintering kinetics under the restricted set of assumptions. Studies on more realistic final-stage sintering microstructure show that the evolution observed in the simulation closely resembles microstructures of real sintered materials over a wide range of diffusivity, initial porosity, and initial pore sizes. Pore shrinkage, grain growth, pore breakaway, and reattachment have all been observed. The porosity decreases monotonically with sintering time and scales with the initial porosity and diffusivity along the grain boundary. Deviations from equilibrium pore shapes under slow surface diffusion or fast grain-boundary diffusion conditions yield slower than expected sintering rates.