We present a numerical solution for the mathematical modeling of the hot pressing
process applied to medium density fiberboard. The model is based on the work of
Humphrey (1982), Humphrey and Bolton (1989), and Carvalho and Costa (1998) with some
modifications and extensions in order to take into account mainly the convective effects
on the phase change term and also a conservative numerical treatment of the resulting
system of partial differential equations.