Abstract
The primary focus of this paper is on studying different numerical models for drying of wet wood particles. More specifically, the advantages and disadvantages of the models, with respect to numerical efficiency, stability, and accuracy, are investigated. The two basic models that are studied in detail are the thermal drying model and the kinetic rate drying model. The drying models have been implemented in an in-house simulation tool that solves for drying and devolatilization of a one-dimensional cylindrical wood log. It is found that the choice of drying model can significantly influence the computational time associated with the thermal conversion. Furthermore, the occurrence of numerical pressure oscillations in the thermal drying model has been found and investigated. The numerical oscillations are reduced by introducing an evaporation fraction, fevap. When the thermal drying model is applied, the drying zone is very thin, commonly only including one grid point, which can result in numerical instabilities. The evaporation fraction allows the smearing of the drying zone by reducing the heat flux used for evaporation of liquid water and using the residual heat flux for heating the grid points. Reducing the evaporation fraction also resulted in reduced CPU times. It was found that model accuracy was not significantly influenced by the choice of drying model.