Numerical analysis for the conjugate heat transfer of skin under various temperature conditions of contrast therapy

In this paper, the contrast therapy of skin was numerically investigated by solving the conjugate heat transfer problem. A finite volume method based on the SIMPLE algorithm was adopted to solve the axisymmetric incompressible Navier-Stokes equations, coupled with an energy equation. These 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 a stimulator for contrast therapy. From the numerical simulations, the regions with cold/hot threshold temperatures were found for five input temperature profiles. It was shown that the temperature varies mildly for different input profiles as the depth increases, owing to the Pennes effect. The input temperatures for effective hot/cold stimulation of dermis layer were found to be 47°C and 7°C, respectively. The present numerical results will be used for finding an optimal temperature profile of a stimulator for contrast therapy.