Contrast therapy of skin was numerically investigated by solving the conjugate heat transfer problem of the skin under the stimulator that periodically imposes hot and cold stimulations on the skin. A finite volume method based on the SIMPLE algorithm was adopted to solve the axisymmetric incompressible Navier-Stokes equations coupled with an energy equation. Those equations are strongly coupled with the Pennes bio-heat equation in order to consider the effect of blood perfusion rate. We investigated the thermal response of skin at some selected depths for various input temperature profiles of the stimulator. The present numerical results will be used in designing the optimal temperature profile of the stimulator.