A model, translated into an algorithm, to simulate nonisothermal sintering of polydispersed glass undergoing concurrent crystallisation is proposed in this article. The model is based on the two classical sintering stages proposed by Fenkel and Mackenzie-Shuttleworth. The input parameters are green density, particle size distribution, surface tension, viscosity, number of surface nucleation sites and crystal growth velocity. The models have been tested using physical parameter and densification rates for two distinct glasses: a devitrifying cordierite glass having a narrow distribution of jagged particles and an alumino-bososilicate glass, ABS, having widely polydispersed jagged particles. The experimental densities of the compacts saturate at slightly different values than predicted.