Contact thermal resistance between two individual silicon nanowires is investigated by nonequilibrium molecular dynamic simulations as a function of temperature, overlap, bonding strength and spacing between them. The results indicate that contact thermal resistance per unit area increases with temperature increasing. The increasing overlap leads to the increase of the contact areas, which enhances the per unit area contact thermal resistance. With a weakened interfacial van der Waals bonding strength, the contact thermal resistance per unit area increases significantly. Additionally, a method to verify the effect of the bonding strength is used by changing the interfacial spacing, and a reasonable result is observed.