A mathematical model for the sintering of ceramics with bimodal pore size distributions at intermediate and final stages is developed. It considers the simultaneous effects of coarsening by surface diffusion, and densification by grain boundary diffusion and lattice diffusion. This model involves population balances for the pores in different zones determined by each porosimetry peak, and is able to predict the evolution of pore size distribution function, surface area, and porosity over time. The model is experimentally validated for the sintering of lime and it is reliable in predicting the so called “initial induction period” in sintering, which is due to a decrease in intra‐aggregate porosity offset by an increase inter‐aggregate porosity. In addition, a novel methodology for determination of mechanisms based on the analysis of the pore size distribution function is proposed, and with this, it was demonstrated that lattice diffusion is the controlling mechanism in the CaO sintering. © 2016 American Institute of Chemical Engineers AIChE J , 63: 893–902, 2017