A novel particle-to-particle contact thermal resistance model is proposed in this study, based on the concept of analogy between force and heat transfer. The distribution of heat flux on the contact surface is assumed to resemble the distribution of stress, eliminating the issue of temperature singularity at the contact edge. The model is validated against existing theories and experiments, showing good agreement. The combined use of the model and the thermal discrete element method is applied to analyze temperature and effective thermal conductivity distributions in a pebble bed within a high-temperature test unit. The average effective thermal conductivity obtained from thermal conduction is found to be 2.99 and 2.61 W/(m·K) for power inputs of 20 kW and 82 kW, respectively. An increase in the outer wall temperature from 200 °C to 1000 °C results in an approximate 20 % increase in the effective thermal diffusivity.