A new one-dimensional heat conduction model for predicting stem heating during fires is presented. The model makes use of moisture and temperature dependent thermal properties for bark and wood. Also, the thermal aspects of the processes of bark swelling, desiccation, and devolatilization are treated in an approximate fashion. Simulation with a surface flux boundary condition requires that these phenomena be accounted for. Previous models have used temperature-time boundary conditions, which prevents them from being directly coupled to fire behavior models. This model uses a flux-time profile for its boundary condition, making it possible to eventually couple it to fire behavior models. Cambial mortality predictions are made through the incorporation of a cell mortality model. The model was developed and validated with laboratory experiments on four species.