Abstract: In order to reveal the heat transfer mechanism of heated tobacco products (HTPs) in the process of smoking, a porous media-based gas flow model and a gas-solid two-phase local heat balance model were established to simulate the temperature distribution of a leading HTP in the process of smoking, and the accuracy of the heat transfer models was verified by measuring the surface temperature of the tobacco stick with infrared thermal imaging method and the exit airflow or aerosol temperature of each functional section of the stick with a thermocouple. The results showed that: 1) The error between the measured and simulated values of the maximum stick surface temperature was less than 5 ℃, indicating that the models were practically accurate. 2) The overall trend of the exit airflow temperatures at the center of each functional segment of the stick were consistent with the experimental and simulated values, but there were differences in the absolute values due to the mass transfer occurred in the flow process of the aerosol with high heat capacity generated during puffing. 3) The simulation results of the temperature fields and airflow fields of the HTP showed that the temperature distribution of the heating device was not uniform, the local temperature in the upper part of the heating device was higher, especially near the stick where the temperature reached a maximum value of 85 ℃ after puffing. The maximum flow velocity was in the cavity section and at the inlet of the airflow during puffing, the airflow zone was not close to the shell of the heating device, which reduced the heat dissipation effect of the inlet cold airflow on the heating device surface. 4) The calculation results of the heat dissipation rate of the HTP showed that during the whole smoking process, the heat dissipation rate of the surface of the heating device increased relatively smoothly with time, and the heat dissipation rate curve of the exposed surface of the stick had a zigzag upward trend with the maximum total heat loss of about 0.95 W.