A three-dimensional model of transient bioheat transfer in the lower extremity during cryotherapy

The purpose of this research was to create a computational model of the human thigh undergoing cryotherapy. The tissue temperatures were measured for five cold pack temperatures of -8°C, -4°C, 0°C, 4°C, and 8°C in addition to six different time intervals of cold application and ice removal. The depth of cold penetration and duration of local tissue cooling were investigated at 10 points during 30 min of application and 7 h of post-application. The model was created in CATIA, using a mid-axial cut of the human thigh MRI without pathology. After validation by the available clinical data, this research applied the finite-volume discretization method to solve bioheat transfer equations. A 16°C decrease in the cold pack temperature reduced the tissue temperatures located 1 and 2 cm below the fat by almost 3.34°C and 1.4°C, respectively, after 30 min of cold application. It took the tissues 10-15 min to start cooling down, and the temperature reached its plateau after 100 min. Thirty minutes of cold application declined the superficial tissue and deep tissue temperatures near the bone by 22.59°C and 0.48°C, respectively. Intense cryotherapy led to an insignificant change in the deep tissue temperature at 2 cm and deeper below the fat tissue. After ice removal, tissues continued cooling down for about 8 min until 40 min, depending on the tissue depth. This study proposed a 100-min cold therapy with 10 min of ice removal to optimize tissue cooling.