Abstract

The objective of the present work was to simulate three-dimensional drying of wood, based on the concept of the effective diffusion coefficient. For this, we used a conventional drying process of radiata pine in a controlled environment (dry and wet bulb temperatures of 44 and 36 oC, respectively). This experiment allowed us to gather data on: a) the spatial distribution of moisture content at different drying times, b) the drying curve, and c) the surface emission coefficient. A differential, partial, non- linear, second-order mathematical model was used, and the exponential functions of the effective diffu- sion coefficient characterized the three-dimensional orthotropic (radial, tangential, and longitudinal) transport of moisture content in the wood. This mathematical model was integrated numerically through the control volume finite element method, which contemplates: a) tetrahedral elements for dis- cretization, b) implicit Euler method for the time dif- ferential, and c) Gauss-Seidel with successive over- relaxation method to resolve the linear equation sys- tem. We compared these results for the three- dimensional spatial distributions of moisture content after 22, 44, 66, and 88 (h) of drying, and the result- ing drying curves were in good agreement with the experiments.