Abstract During selective laser sintering (SLS), the microstructure evolution and local temperature variation interact mutually. Application of conventional isothermal sintering model is thereby insufficient to describe SLS. In this work, we construct our model from entropy level, and derive the non-isothermal kinetics for order parameters along with the heat transfer equation coupled with microstructure evolution. Influences from partial melting and laser-powder interaction are also addressed. We then perform 3D finite element non-isothermal phase-field simulations of the SLS single scan. To confront the high computation cost, we propose a novel algorithm analogy to minimum coloring problem and manage to simulate a system of 200 grains with grain tracking algorithm using as low as 8 non-conserved order parameters. Specifically, applying the model to SLS of the stainless steel 316L powder, we identify the influences of laser power and scan speed on microstructural features, including the porosity, surface morphology, temperature profile, grain geometry, and densification. We further validate the first-order kinetics of the transient porosity during densification, and demonstrate the applicability of the developed model in predicting the linkage of densification factor to the specific energy input during SLS.