Convective drying analysis of three-dimensional porous solid by mass lumping finite element technique

• Published in: Heat and mass transfer 2008-02, Vol.44 (4), p.401-412
• Main Authors: Murugesan, K., Lo, D. C., Young, D. L., Chen, C. W., Fan, C. M.
• Format: Article
• Language: English
• Subjects: Applied sciences; Devices using thermal energy; Dryers; Energy; Energy. Thermal use of fuels; Engineering; Engineering Thermodynamics; Exact sciences and technology; Heat and Mass Transfer; Industrial Chemistry/Chemical Engineering; Original; Thermodynamics
• ISSN: 0947-7411; 1432-1181
• DOI: 10.1007/s00231-007-0260-9

A numerical analysis of convective drying of a 3D porous solid of brick material is carried out using the finite element method and mass lumping technique. The energy equation and moisture transport equations for the porous solid are derived based on continuum approach following Whitaker’s theory of drying. The governing equations are solved using the Galerkin’s weighted residual method, which convert the governing equations into discretized form of matrix equations. The resulting capacitance matrices are made diagonal matrices by following the classical row-sum mass lumping technique. Hence with the use of the Eulerian time marching scheme, the final equations are reduced to simple algebraic equations, which can be solved directly without using an equation solver. The proposed numerical scheme is initially validated with experimental results for 1D drying problem and then tested by application to convective drying of 3D porous solid of brick material for four different aspect ratios obtained by varying the cross section of the solid. The mass lumping technique could correctly predict the wet bulb temperature of the solid under evaporative drying conditions. A parametric study carried out for three different values of convective heat transfer coefficients, 15, 30 and 45 W/m 2  K shows an increased drying rate with increase in area of cross section and convective heat transfer coefficient. The proposed numerical scheme could correctly predict the drying behavior shown in the form of temperature and moisture evolutions.