We developed a mathematical model to estimate the increase in firefighters' core body temperature from energy expenditure (EE) measured by accelerometry to prevent heat illness during firefighting. Wearing firefighter personal protective equipment, seven male subjects aged 23-42 years underwent a graded walking test on a treadmill while esophageal temperature (T es ) and skin temperature were measured with thermocouples and EE was measured with a tri-axial accelerometer. To estimate the increase in T es from EE, we proposed a mathematical model composed of the heat capacity of active muscles (C 1 , kcal•°C −1 ), the heat capacity of the sum of resting muscles and skin (C 2 ), the resistance to heat flux from C 1 to C 2 (R 1 ,°C•min•kcal −1 ), and the resistance from C 2 to the skin surface (R 2 ). We determined the parameters while minimizing the differences between the estimated and measured changes in T es profiles during graded walking. We found that C 1 and C 2 in individuals were highly correlated with their body weight (kg) and body surface area (m 2 ), respectively, whereas R 1 and R 2 were similar across subjects. When the profiles of measured T es (y) and estimated T es (x) were pooled in all subjects, they were almost identical and were described by a regression equation without an intercept, y = 0.96x (r = 0.96, P < 0.0001), with a mean difference of − 0.01 ± 0.12°C (mean ± SD) ranging from − 0.18 to 1.56°C of the increase in T es by Bland-Altman analysis. Thus, the model can be used for firefighters to prevent heat illness during firefighting.