Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 120736, 8 pages

http://dx.doi.org/10.1155/2013/120736

Research Article

Cortex Effect on Vacuum Drying Process of Porous Medium

Zhijun Zhang,¹ Shiwei Zhang,¹ Tianyi Su,¹ and Shuangshuang Zhao²

¹School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110004, China

²Shenyang Aircraft Design and Research Institute, Shenyang 110035, China

Received 25 February 2013; Accepted 5 April 2013

Academic Editor: Jun Liu

Copyright © 2013 Zhijun Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Corns, fruits, and vegetables are usually used as porous medium in drying process. But in fact, it must be considered as the cortex effect on mass transfer because the mass transfer of cortex is very difficult than inner medium. Based on the theory of heat and mass transfer, a coupled model for the porous medium vacuum drying process with cortex effect is constructed. The model is implemented and solved using COMSOL software. The water evaporation rate is determined using a nonequilibrium method with the rate constant parameter that has been studied. The effects of different vapor pressures (1000, 5000, and 9000 Pa), initial moisture contents (0.3, 0.4, and 0.5 water saturation), drying temperatures (323, 333, and 343 K), and intrinsic permeability for cortex part (10⁻¹³, 10⁻¹⁴, 10⁻¹⁵ m²) on vacuum drying process were studied. The results facilitate a better understanding of the porous medium vacuum drying process that nearer to the reality.

1. Introduction

Scientists and engineers in China are currently studying vacuum drying equipment that could be used in corn drying [1–3]. However, the corn vacuum drying theory remains unclear. Hypothesized that corn is a porous medium, the vacuum drying of corn is a complicated heat and mass transfer process that has been the subject of intensive research [4–7]. All vacuum drying models have to address the water phase change during numerical solving. In one method, the vapor pressure is equal to its equilibrium value [8–11]. Another method is nonequilibrium method [12–16]. As the porous medium, the heat and mass transfer in vacuum drying process has been studied by nonequilibrium method [17]. The water evaporation rate is determined using a nonequilibrium method with the rate constant parameter . values of 1, 10, 1000, and 10000 are simulated. The effects of different vapor pressures (1000, 5000, and 9000 Pa), initial moisture contents (0.4, 0.5, and 0.6 water saturation), drying temperatures (323, 333, and 343 K), and intrinsic permeability (10⁻¹³, 10⁻¹⁴, 10⁻¹⁵ m²) are studied. It was observed that the temperature increased quickly at the start of drying and then lowered gradually. As the drying process continued, the temperature increased slowly. In the absence of free water, temperature increased rapidly. As the drying process concluded, the temperature remained unchanged. The water evaporation rate could not be obtained during the porous medium vacuum drying process. The rate constant parameter is essential to the nonequilibrium method. When , the simulation of the drying process was not evidently affected. Vapor pressure and heat transfer affected the transfer of mass. A similar effect was found in the initial moisture and the heat temperature. Intrinsic permeability had a greater effect on the drying process. In the study above, it is as the uniform porous medium. But in fact, it must be considered as the cortex effect on mass transfer because the mass transfer of cortex is more difficult than inner medium. The study of about cortex effect is very few.

In this paper, heat and mass transfer of porous medium with cortex in the vacuum drying process is implemented by using a nonequilibrium method. The effects of vapor pressure, initial moisture, heat temperature, and cortex intrinsic permeability on the drying process were then examined.

2. Physical Model

A physical one-dimensional (1D) model that explains the drying process is shown in Figure 1. The heat and mass transfer is considered only in the direction. Because consider the cortex effect. The heat and mass transfer is through the surface of medium that is the cortex position. The total height of the porous medium is 1.05 cm, and the cortex height is 0.05 cm.

Figure 1: 1D model of porous medium with cortex.

3. Mathematical Model

The porous medium consists of a continuous rigid solid phase, an incompressible liquid phase (free water), and a continuous gas phase that is assumed to be a perfect mixture of vapor and dry air, considered as ideal gases. For a mathematical description of the transport phenomenon in a porous medium, we adopt a continuum approach, wherein macroscopic partial differential equations are achieved through the volume averaging of the microscopic conservation laws. The value of any physical quantity at a point in space is given by its average value on the averaging volume centered at this point.

The moisture movement of the inner porous medium is liquid water and vapor movement; that is, the liquid water could become vapor, and the vapor and liquid water are moved by the pressure gradient. The heat and mass transfer theory could be found in everywhere [8].

The compressibility effects of the liquid phase are negligible, and the phase is homogeneous:

The solid phase is rigid and homogeneous:

The gaseous phase is considered an ideal gas. This phase ensures that

The assumption of the local thermal equilibrium between the solid, gas, and liquid phases involves

Mass conservation equations are written for each component in each phase. Given that the solid phase is rigid, the following is given:

The averaged mass conservation of the dry air yields

For vapor,

For free water,

For water, the general equation of mass conservation is obtained from the sum of the conservation equations of vapor () and free water (). The general equation is written as follows:

For the Darcy flow of vapor,

For the Darcy flow of air,

where the gas and free water velocity is given by

The effective diffusion coefficient [8] is given by

The vapor fraction in mixed gas is given by

The pressure moving the free water is given by

For capillary pressure,

The saturation of free water and gas is

Free water relative permeability is given by

Gas relative permeability is given by

The water phase change rate is expressed as

Water saturation vapor pressure is given by

By considering the hypothesis of the local thermal equilibrium, the energy conservation is reduced to a unique equation:

4. Boundary Condition and Parameters

The air pressure on the external surface of the porous medium is fixed, and the boundary condition for air is given by

The boundary condition for vapor at the surface of the porous medium is given by

To simulate the vapor pressure of the vacuum drying chamber effect on the drying process, four different vapor pressure boundary values are used.

The boundary condition for free water at the top of the porous medium is

The boundary condition at the surface of the porous medium is

Three different values are used in the simulation.

The initial moisture of the porous medium is represented by the liquid water saturation; different initial water saturation values are used. To compare the effects, drying base moisture content (d.b.) was also used, as shown in (9). The water phase change rate is used as 1000 that has been studied before [17]. Intrinsic permeability of inner medium is 10⁻¹³ m². Intrinsic permeability of cortex is 10⁻¹⁴, 10⁻¹⁵ m². In order to compare the results, the intrinsic permeability of cortex 10⁻¹³, that is no cortex effect is simulated. In order to easily converge in, the parameter changer is used a smooth method. The modeling parameters are shown in Table 1.

Table 1: Parameters used in the simulation process.

Parameter
Symbol
Value
Unit

Rate constant parameter

1000
s⁻¹

Intrinsic permeability of inner medium

10⁻¹³
m²

Intrinsic permeability of cortex

10⁻¹³, 10⁻¹⁴, 10⁻¹⁵
m²

Initial water saturation

0.5, 0.4, 0.3

Vapor pressure of vacuum drying chamber

1000, 5000, 9000
Pa

Heat temperature

323, 333, 343
K

Porosity

0.615

Solid density

476
kg m⁻³

Air pressure of vacuum drying chamber

0.001
Pa

5. Numerical Solution

COMSOL Multiphysics 3.5a was used to solve the set of equations. COMSOL is an advanced software used for modeling and simulating any physical process described by partial derivative equations. The set of equations introduced above was solved using the relative initial and boundary conditions of each. COMSOL offers three possibilities for writing the equations: (1) using a template (the Fick law and the Fourier law), (2) using the coefficient form (for mildly nonlinear problems), and (3) using the general form (for most nonlinear problems). Differential equations in the coefficient form were written using an unsymmetric-pattern multifrontal method. We used a direct solver for sparse matrices (UMFPACK), which involves significantly more complicated algorithms than solvers used for dense matrices. The main complication is the need to handle the fill-in in factors and efficiently.

A two-dimensional (2D) grid was used to solve the equations using COMSOL Multiphysics 3.5a. Given the symmetry condition setting at the left and the right sides, The 2D is applied to the the 1D model shown in Figure 1. The mesh consists of elements (2D), and time stepping is 1 (from 0 s to 100 s of solution), 5 (from 100 s to 200 s of solution), 20 (from 200 s to 1000 s of solution), 30 (from 1000 s to 2000 s of solution), 40 (from 2000 s to 4000 s of solution), 50 (from 4000 s to 20000 s of solution), and 100 (from 20000 s to 50000 s of solution). Several grid sensitivity tests were conducted to determine the sufficiency of the mesh scheme and to ensure that the results are grid independent. The maximum element size was established as . A backward differentiation formula was used to solve time-dependent variables. Relative tolerance was set to , whereas absolute tolerance was set to . The simulations were performed using a Tongfang PC with Intel Core 2 Duo processor with 3.0 GHz processing speed, and 4096 MB of RAM running Windows 7.

6. Results and Discussion

6.1. Effect of Boundary Condition

In the study before [17], the heat and mass transfer direction is the same. That is, the heat boundary is on the bottom of the medium. The heat transfer is from the bottom to top. The mass transfer boundary is on the top of bottom. The mass transfer is from the bottom to top by the pressure driving. But in this study, the cortex is considered. And the heat transfer and mass transfer are on the reversed direction as Figure 1. The heat and mass transfer with dryer is on the surface of the medium. The both condition was simulated to be compared as Figure 2. The simulation parameters are , ,  K,  Pa,  m², and no cortex. The moisture curve is just little different. The drying rate is little larger at the drying initial stage for same direction heat and mass transfer, and then it is little lower. The other condition is just reversed. The drying rate is little lower at the initial stage for reversed direction heat and mass transfer, and then it is little larger. The drying time for both is 6 hours.

Figure 2: The effect of boundary condition, heat, and mass transfer is same and is reverse direction.

6.2. Effect of Cortex Resistance

Intrinsic permeability of porous medium is an inherent property and cannot be changed, and measuring it is difficult. Intrinsic permeability has a greater effect because the transfer of free water and vapor is affected by (13). Usually, the cortex intrinsic permeability is less than inner medium. Figure 3 is the moisture curves of no cortex ( and , , and  m², resp.), with cortex ( and and  m² resp.). The other simulation parameters are , ,  K, and  Pa. The drying time became evidently longer as the intrinsic permeability was reduced because the moisture movement velocity was lowered at the same pressure gradient. The different is not so obviously as the  m², that is lower 10-fold than inner medium. But when  m², the drying time is increased to 2-fold nearly. In order to compare, the no cortex,  m², and  m² was shown. Because the numerical calculation is not convergence, the cortex  m²,  m² result not gotten. The reason is that the material parameter changer is too big.

Figure 3: Moisture curves of no cortex (, , , and  m²), (, ,  m²).

6.3. Effect of Vapor Pressure in Vacuum Drying Chamber

The pressure of a vacuum drying chamber, especially vapor pressure, plays an important role in the vacuum drying process and is also linked to the drying cost. The moisture curves of , 5000, and 9000 Pa are shown in Figures 4, 5, and 6 with different cortex intrinsic permeability. The other simulation parameters are , , and  K. The vapor pressure has a greater effect on the drying process; a lower vapor pressure results in greater pressure degradation. The movements of free water and vapor as well as the free water evaporation rate are quicker, as given by (13) and (21), respectively. But the pressure effect is lower by the cortex intrinsic permeability reduced.

Figure 4: Moisture curves of different , 5000, 9000 Pa,  m²,  m².

Figure 5: Moisture curves of different , 5000, 9000 Pa,  m²,  m².

Figure 6: Moisture curves of different , 5000, 9000 Pa, no cortex,  m².

6.4. Effect of Initial Moisture Content

The effect of initial moisture content on the moisture curve is shown in Figures 7, 8, and 9 for , 0.4, 0.3 with different cortex intrinsic permeability. To compare the results, moisture is represented by the moisture ratio . The other parameters are ,  Pa, and  K. The drying time is about 10 hours, 6.5 hours, and 6 hours for , 0.4, 0.3, respectively. In the same intrinsic permeability, the drying time is almost the same. The initial moisture has a less effect when initial moisture is 0.4 and 0.3, especially cortex intrinsic permeability  m². The drying rate is little larger for , and the drying rate is almost the same for , 0.3.

Figure 7: Moisture curves at different , 0.4, 0.3,  m²,  m².

Figure 8: Moisture curves at different , 0.4, 0.3,  m²,  m².

Figure 9: Moisture curves at different , 0.4, 0.3, no cortex,  m².

6.5. Effect of Heat Temperature

The effect of heat temperature on moisture is shown in Figures 10, 11, and 12 for the heat temperature , 333, and 343 K. The other parameters are , , and  Pa. The effect of heat temperature is obviously in all simulation. The drying time is more than 14 hours for  K,  m², and  m². And it is about 10 hours and 8 hours for , 343 K,  m², and  m². The drying time is more than 9 hours for  K,  m², and  m². It is about 7 hours and 6 hours for , 343 K,  m², and  m². The drying time is more than 7 hours for  K,  m², and  m². It is about 6 hours and 5 hours for , 343 K,  m², and  m².

Figure 10: Moisture curves at different , 333, 343 K,  m²,  m².

Figure 11: Moisture curves at different , 333, 343 K,  m²,  m².

Figure 12: Moisture curves at different , 333, 343 K, no cortex,  m².

7. Conclusion

A coupled model of porous medium vacuum drying with cortex effect based on the theory of heat and mass transfer was implemented in this paper. The drying rate is little higher at the drying initial stage for the same direction of heat and mass transfer, and then it is little lower. The drying time became evidently longer as the intrinsic permeability was reduced because the moisture movement velocity was lowered at the same pressure gradient. The difference is not so obvious as the  m², that is lower 10-fold than inner medium. But when  m², the drying time is increased to 200% nearly. The vapor pressure has a greater effect on the drying process; a lower vapor pressure results in greater pressure degradation. But the pressure effect is lower by the reduction of the cortex intrinsic permeability. For initial moisture effect, in the same intrinsic permeability, the drying time is almost the same. The initial moisture is a less effect when initial moisture is 0.4 and 0.3, especially cortex intrinsic permeability  m². The drying rate is little larger for , and the drying rate is almost the same for , 0.3. The effect of heat temperature is obviously in all simulation.

Nomenclature

: | Diagonal tensor
: | Diffusivity (m² s⁻¹)
: | Diffusion tensor (m² s⁻¹)
: | Gravity vector (m s⁻²)
: | Intrinsic averaged enthalpy (Jkg⁻¹)
: | Water phase rate (kg s⁻¹ m⁻³)
: | Intrinsic permeability (m²)
: | Relative permeability
: | Mass (kg)
: | Outer unit normal to the product
: | Pressure (Pa)
: | Capillary pressure (Pa)
: | Universal Gas constant (J kmol⁻¹ K⁻¹)
: | Saturation
: | Time (s)
: | Temperature (K)
: | Moisture content (in dry basis).

Greek Letters

: | Latent of phase change (Jkg⁻¹)
: | Effective thermal conductivity tensor (W)
: | Viscosity (kg )
: | Density (kg)
: | Vapor fraction.

Subscripts

: | Dry air
: | Gas
: | Liquid
: | Solid
: | Vapor
sat: | Vapor saturation
in: | Inner medium
cor: | Cortex.

Mathematical Operators

: | Gradient operator
: | Divergence operator.

Acknowledgment

This research was supported by the National Natural Science Foundation of China (Grant nos. 31000665 and nos. 51176027).

References

C. H. Xu, Z. J. Zhang, S. W. Zhang, and X. He, “Probe into the structure of tower continuous vacuum dryer,” in Proceedings of the 5th Asia-Pacific Drying Conference (ADC '07), pp. 1261–1267, August 2007.

Z. J. Zhang, C. H. Xu, S. W. Zhang, and X. He, “The study of corn low temperature continuous tower type vacuum dryer,” in Proceedings of the 5th Asia-Pacific Drying Conference (ADC '07), pp. 330–337, August 2007.

Z. Zhang, C. Xu, S. Zhang, and L. Zhao, “Computer simulation of flow field in tower continuous vacuum dryer,” in Proceedings of the International Conference on Computer Science and Information Technology (ICCSIT '08), pp. 534–538, September 2008.

I. Yasuaki and A. P. S. Selvadurai, Transport Phenomena in Porous Media, Aspects of Micro/Macro Behaviour, Springer, 2012.

A. K. Haghi, “Transport phenomena in porous media: a review,” Theoretical Foundations of Chemical Engineering, vol. 40, no. 1, pp. 14–26, 2006.

S. J. Kowalski, Drying of Porous Materials, Springer, 2007.

B. Jacob and Y. Bachmat, Introduction to Modeling of Transport Phenomena in Porous Media, Springer, 1990.

A. Erriguible, P. Bernada, F. Couture, and M. A. Roques, “Simulation of vacuum drying by coupling models,” Chemical Engineering and Processing, vol. 46, no. 12, pp. 1274–1285, 2007.

A. Erriguible, P. Bernada, F. Couture, and M. A. Roques, “Modeling of heat and mass transfer at the boundary between a porous medium and its surroundings,” Drying Technology, vol. 23, no. 3, pp. 455–472, 2005.

K. Murugesan, H. N. Suresh, K. N. Seetharamu, P. A. Aswatha Narayana, and T. Sundararajan, “A theoretical model of brick drying as a conjugate problem,” International Journal of Heat and Mass Transfer, vol. 44, no. 21, pp. 4075–4086, 2001.

P. Perré and I. W. Turner, “A dual-scale model for describing drier and porous medium interactions,” AIChE Journal, vol. 52, no. 9, pp. 3109–3117, 2006.

S. S. Torres, W. Jomaa, J.-R. Puiggali, and S. Avramidis, “Multiphysics modeling of vacuum drying of wood,” Applied Mathematical Modelling, vol. 35, no. 10, pp. 5006–5016, 2011.

S. S. Torres, J. R. Ramírez, and L. Méndez-Lagunas, “Modeling plain vacuum drying by considering a dynamic capillary pressure,” Chemical and Biochemical Engineering Quarterly, vol. 25, no. 3, pp. 327–334, 2011.

A. Warning, A. Dhall, D. Mitrea, and A. K. Datta, “Porous media based model for deep-fat vacuum frying potato chips,” Journal of Food Engineering, vol. 110, no. 3, pp. 428–440, 2012.

A. Halder, A. Dhall, and A. K. Datta, “An improved, easily implementable, porous media based model for deep-fat frying. Part I: model development and input parameters,” Food and Bioproducts Processing, vol. 85, no. 3 C, pp. 209–219, 2007.

A. Halder, A. Dhall, and A. K. Datta, “An improved, easily implementable, porous media based model for deep-fat frying. Part II: results, validation and sensitivity analysis,” Food and Bioproducts Processing, vol. 85, no. 3, pp. 220–230, 2007.

Z. Zhang and N. Kong, “Nonequilibrium thermal dynamic modeling of porous medium vacuum drying process,” Mathematical Problems in Engineering, vol. 2012, Article ID 347598, 22 pages, 2012.