An analytical approach for treating problems involving oscillatory heat source is presented. The transient temperature profile involving circular, rectangular, and parabolic heat sources undergoing oscillatory motion on a semi-infinite body is determined by integrating the instantaneous solution for a point heat source throughout the area where the heat source acts with an assumption that the body takes all the heat. An efficient algorithm for solving the governing equations is developed. The results of a series simulations are presented, covering a wide range of operating parameters including a new dimensionless frequency and the dimensionless oscillation amplitude , whose product can be interpreted as the Peclet number involving oscillatory heat source, . Application of the present method to fretting contact is presented. The predicted temperature is in good agreement with published literature. Furthermore, analytical expressions for predicting the maximum surface temperature for different heat sources are provided by a surface-fitting method based on an extensive number of simulations.