This research uses numerical modeling to study the jet cooling properties of hot‐rolled seamless steel tubes under axial movement and circumferential rotation. The orthogonal test is used to examine the effects of axial velocity V, rotation speed R, jet distance S, and jet Reynolds number Re on the surface temperature distribution of steel tubes. The findings demonstrate that the temperature of the steel tube exhibits a “similar sinusoidal function” distribution along the circumference, with various angular frequencies and starting phases as a result of various cooling sites. It is determined that the axial velocity has the greatest influence on the temperature drop during the cooling process of steel tube, followed by the jet Reynolds number, the jet distance, and the rotation speed. When the axial velocity is 0.3–0.7 m s−1, the rotation speed is 30–45 rpm, the jet distance is 75–100 mm, and the jet Reynolds number is 38 000–51 000, the steel tube cooling efficiency is the best. The mathematical models of temperature drop ΔT and axial velocity V, rotation speed R, jet distance S, and Reynolds number Re are established. The maximum deviation of temperature drop between numerical simulation and mathematical model is only 1.02 °C with high prediction accuracy.