A comprehensive heat and mass transfer computational model is used to analyse the intricate two-way coupling arising from the activated chemical reactions involved in the heat treatment of wood. The 2D version of a drying code known as TransPore is used to simulate the coupled heat and mass transfer phenomena. This code accounts for the internal pressure in the porous medium. The pyrolysis model describing the chemical reactions occurring in the main constituents within the cell walls of wood (cellulose, hemicelluloses and lignins) is derived using data taken from the literature. Refined computational strategies were required to address the two-way coupling between the heat and mass transfer and chemical mechanisms, including the thermal activation of the chemical reactions, together with the treatment of heat sources (or sinks) and the production of volatiles. The experimental set-up allows the overall weight loss, and the internal temperature and pressure at specific locations within the board to be determined during processing. The reported simulations highlight that the model is able to capture two particular phenomena observed during the heat treatment of wood: the double pressure peak due to water evaporation and volatiles production; and the temperature overshoot during the heat treatment phase.