Hot-pressing process modeling for medium density fiberboard (MDF)

Noberto M. Nigro,Mario A. Storti

DOI: https://doi.org/10.48550/arXiv.math/0010173

2001-06-23

Abstract:In this paper we present a numerical solution for the mathematical modeling of the hot-pressing process applied to medium density fiberboard. The model is based in the work of Humphrey[82], Humphrey and Bolton[89] and Carvalho and Costa[98], 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.

Numerical Analysis,Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

This paper attempts to solve the heat and mass transfer problems involved in the hot - pressing process of medium - density fiberboard (MDF). Specifically, the main research objectives include:

1. **Improve heat conduction efficiency**: By introducing the convective effect of water vapor evaporation and condensation, improve the heat transfer model from the press plate to the center of the board.
2. **Optimize parameter adjustment**: Predict and optimize the process parameters (such as water content, press plate cycle time and its temperature history) that affect the final product to obtain the desired density distribution.
3. **Improve the numerical model**: Based on existing models (such as the work of Humphrey [1982] and Carvalho and Costa [1998]), propose a more accurate energy balance equation and use a conservative numerical treatment method to solve the system of partial differential equations.

### Main content overview

#### 1. Mathematical model of hot - pressing
The paper describes the multiphase model, energy balance equation, and mass balance equation of steam and gas mixtures in the hot - pressing process. These equations take into account the complex heat and mass transfer phenomena inside the material, especially the influence of the evaporation and condensation processes.

- **Energy balance equation**:
\[
\rho_s C_p\frac{\partial T}{\partial t}=\nabla\cdot\left(k\nabla T - \rho_v V_g(C_{pv}T+\lambda + Q_l)\right)-\dot{m}(\lambda + Q_l)
\]
where \(T\) is the temperature, \(t\) is the time, \(C_p\) is the specific heat capacity, \(k\) is the thermal conductivity, \(\rho_s\) is the density of the dry board, \(\rho_v\) is the density of water vapor, \(V_g\) is the volume - averaged gas velocity, \(C_{pv}\) is the specific heat capacity of water vapor, \(\dot{m}\) is the evaporation rate, \(\lambda\) is the latent heat of vaporization of free water, and \(Q_l\) is the heat of adsorption.

- **Moisture content balance equation**:
\[
\rho_s\frac{\partial H}{\partial t}=\nabla\cdot[\epsilon D_v\nabla\rho_v - V_g\rho_v]
\]

- **Air mass balance equation**:
\[
\epsilon\frac{\partial\rho_a}{\partial t}=\nabla\cdot[\epsilon D_a\nabla\rho_a - V_g\rho_a]
\]

#### 2. Numerical method
To solve the above system of partial differential equations, the author adopts the finite element method and uses the SUPG (Streamline Upwind - Petrov Galerkin) scheme to stabilize the numerical scheme and avoid spurious oscillations. In addition, an implicit time - integration method is also implemented to improve numerical stability.

#### 3. Physical and transport properties
The calculation methods of physical properties such as thermal conductivity, permeability, diffusion coefficient of water vapor in air, and saturated vapor pressure are discussed in detail to ensure the accuracy and reliability of the model.

### Conclusion
The numerical model proposed in this paper not only takes into account the convective effect in the evaporation and condensation processes but also uses a conservative numerical treatment method to solve the system of partial differential equations. The two - dimensional simulation results show that this improvement can better estimate the boundary effect and improve the heat conduction efficiency in the hot - pressing process.

Through these improvements, researchers can more accurately predict and optimize the hot - pressing process of medium - density fiberboard, thereby reducing production costs and improving product quality.