The simulation is based on the unsteady-state nonlinear diffusion of moisture and heat with respect to the orthotropic nature of wood. The moisture gradient and the temperature gradient are set as driving forces for mass and heat transfer. The model respects the dependence of material coefficients on temperature and moisture, the Soret effect, the Duffour effect and the anisotropic nature of wood. A matrix form of the multiphysics model for FEM solver was derived and numerical simulations were performed.