Numerical modelling of refill friction stir spot welding helps in the deep investigation of physics involved in the joining process. The current study investigates four variants of contact models and their influence on thermal prediction in the process. The contact models with constant shear and coulomb friction coefficients and temperature-dependent shear and coulomb friction coefficients were used in the numerical models to investigate the influence of friction coefficients in the prediction of thermal cycles in refill friction stir spot welding. The results from the numerical model are compared and validated with the experimental results from a previous study. The contact model with temperature-dependent coulomb friction coefficient was concluded to be more reliable when compared to the other three contact models. The temperature for this contact model at the centre of the weld is 505 °C, which differs from the temperature recorded in the experiment by 2 . The peak temperatures of numerical models with constant coulomb and shear friction coefficients differ from the experimental temperature at the weld centre by 5.17 and 16.8 respectively. © 2023 Wiley-VCH GmbH.