# Cortex Effect on Vacuum Drying Process of Porous Medium

Authors: Zhijun Zhang; Shiwei Zhang; Tianyi Su; Shuangshuang Zhao
Journal: Mathematical Problems in Engineering
 (2013)

---

## 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 parameterKr 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−13, 10−14, 10−15 m2) on vacuum drying process were studied. The results facilitate a better understanding of the porous medium vacuum drying process that nearer to the reality.


---

## Body

## 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 Kr. Kr 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−13, 10−14, 10−15 m2) 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 Kr≥1000, 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 Figure1. The heat and mass transfer is considered only in the y 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:(1)ρ¯l=constant.
The solid phase is rigid and homogeneous:
(2)ρ¯s=constant.
The gaseous phase is considered an ideal gas. This phase ensures that
(3)ρ¯a=maPaRT¯,ρ¯v=mvP¯vRT¯,P¯g=P¯a+P¯v,ρ¯g=ρ¯a+ρ¯v.The assumption of the local thermal equilibrium between the solid, gas, and liquid phases involves(4)T¯s=T¯g=T¯l=T¯.Mass conservation equations are written for each component in each phase. Given that the solid phase is rigid, the following is given:(5)∂ρ¯s∂t=0.The averaged mass conservation of the dry air yields(6)∂(ε·Sgρ¯a)∂t+∇·(ρ¯aV¯a)=0.For vapor,(7)∂(ε·Sgρ¯v)∂t+∇·(ρ¯vV¯v)=I˙.For free water,(8)∂(ε·Swρ¯l)∂t+∇·(ρ¯lV¯l)=-I˙.For water, the general equation of mass conservation is obtained from the sum of the conservation equations of vapor (v) and free water (l). The general equation is written as follows:
(9)∂W∂t+∇·{1ρ¯s(ρ¯lV¯l+ρ¯vV¯v)}=0,(10)W=ε·Swρ¯l+ε·Sgρ¯v(1-ε)ρ¯s.For the Darcy flow of vapor,(11)ρ¯vV¯v=ρ¯vV¯g-ρ¯gDeff·∇ω¯.For the Darcy flow of air,(12)ρ¯aV¯a=ρ¯aV¯g-ρ¯gDeff·∇ω¯,
where the gas and free water velocity is given by
(13)V¯g=k·krgμg·(∇P¯g-ρgg→),V¯l=k·krlμl·(∇P¯l-ρlg→).
The effective diffusion coefficient [8] is given by
(14)Deff=DB=.The vapor fraction in mixed gas is given by(15)ω¯=ρvρg.
The pressure moving the free water is given by
(16)P¯l=P¯g-P¯c.
For capillary pressure,
(17)P¯c=56.75×103(1-Sl)exp(1.062Sl).
The saturation of free water and gas is
(18)Sg+Sl=1.
Free water relative permeability is given by
(19)krl={(Sl-Scr1-Scr)3Sw>Scr0Sw≤Scr.
Gas relative permeability is given by
(20)krg=Sg.
The water phase change rate is expressed as
(21)I˙=Krmv(aωPsat-Pv)SgεRT.
Water saturation vapor pressure is given by
(22)Psat=101325760×10(8.07131-(1730.63/(233.426+(T-273)))).By considering the hypothesis of the local thermal equilibrium, the energy conservation is reduced to a unique equation:(23)∂ρ¯h¯∂t+∇·{(ρ¯aV¯ah¯a+ρ¯vV¯vh¯v+ρ¯lV¯lh¯l-λe·∇T¯-ΔH·I˙)}=0,λe=(1-ε)λs+ε(Sl+Sg(ωλv+(1-ω)λa)),ρ¯h¯=ρ¯sh¯s+ε·Sgρ¯ah¯a+ε·Sgρ¯vh¯v+ε·Slρ¯lh¯l.


## 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(24)Pa=Pav.The boundary condition for vapor at the surface of the porous medium is given by(25)Pv=Pvb.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(26)n·(-D∇Sw)=0.The boundary condition at the surface of the porous medium is(27)T=Th.Three differentTh 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−13 m2. Intrinsic permeability of cortex is 10−14, 10−15 m2. In order to compare the results, the intrinsic permeability of cortex 10−13, 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
K
r
1000
s−1
Intrinsic permeability of inner medium
k
in
10−13
m2
Intrinsic permeability of cortex
k
cor
10−13, 10−14, 10−15
m2
Initial water saturation
S
l
0
0.5, 0.4, 0.3
Vapor pressure of vacuum drying chamber
P
v
b
1000, 5000, 9000
Pa
Heat temperature
T
h
323, 333, 343
K
Porosity
ε
0.615
Solid density
ρ
s
476
kg m−3
Air pressure of vacuum drying chamber
P
a
b
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 factorsL and U 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 Figure1. The mesh consists of 2×200 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 1e-4. A backward differentiation formula was used to solve time-dependent variables. Relative tolerance was set to 1e-3, whereas absolute tolerance was set to 1e-4. 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 Sl0=0.5, Kr=1000, Th=323 K, Pvb=1000 Pa, kin=10-13 m2, 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 (kin and kcor=10-13, 10-14, and 10-15 m2, resp.), with cortex (kin=10-13 and kcor=10-14 and 10-15 m2 resp.). The other simulation parameters are Sl0=0.5, Kr=1000, Th=323 K, and Pvb=1000 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 kcor=10-14 m2, that is lower 10-fold than inner medium. But when kcor=10-15 m2, the drying time is increased to 2-fold nearly. In order to compare, the no cortex, kcor=10-14 m2, and 10-15 m2 was shown. Because the numerical calculation is not convergence, the cortex kcor=10-16 m2, kin=10-13 m2 result not gotten. The reason is that the material parameter changer is too big.Figure 3
Moisture curves of no cortex (kin=10-13, kcor=10-13, 10-14, and 10-15 m2), (kcor=10-14, 10-15, kin=10-13 m2).


### 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 ofPvb=1000, 5000, and 9000 Pa are shown in Figures 4, 5, and 6 with different cortex intrinsic permeability. The other simulation parameters are Sl0=0.5, Kr=1000, and Th=323 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 differentPvb=1000, 5000, 9000 Pa, kcor=10-15 m2, kin=10-13 m2.Figure 5
Moisture curves of differentPvb=1000, 5000, 9000 Pa, kcor=10-14 m2, kin=10-13 m2.Figure 6
Moisture curves of differentPvb=1000, 5000, 9000 Pa, no cortex, kin=10-13 m2.


### 6.4. Effect of Initial Moisture Content

The effect of initial moisture content on the moisture curve is shown in Figures7, 8, and 9 for Sl0=0.5, 0.4, 0.3 with different cortex intrinsic permeability. To compare the results, moisture is represented by the moisture ratio M/M0. The other parameters are Kr=1000, Pvb=1000 Pa, and Th=323 K. The drying time is about 10 hours, 6.5 hours, and 6 hours for Sl0=0.5, 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 kcor=10-15 m2. The drying rate is little larger for Sl0=0.5, and the drying rate is almost the same for Sl0=0.4, 0.3.Figure 7
Moisture curves at differentSl0=0.5, 0.4, 0.3, kcor=10-15 m2, kin=10-13 m2.Figure 8
Moisture curves at differentSl0=0.5, 0.4, 0.3, kcor=10-14 m2, kin=10-13 m2.Figure 9
Moisture curves at differentSl0=0.5, 0.4, 0.3, no cortex, kin=10-13 m2.


### 6.5. Effect of Heat Temperature

The effect of heat temperature on moisture is shown in Figures10, 11, and 12 for the heat temperature Th=323, 333, and 343 K. The other parameters are Sl0=0.5, Kr=1000, and Pv=1000 Pa. The effect of heat temperature is obviously in all simulation. The drying time is more than 14 hours for Th=323 K, kcor=10-15 m2, and kin=10-13 m2. And it is about 10 hours and 8 hours for Th=333, 343 K, kcor=10-15 m2, and kin=10-13 m2. The drying time is more than 9 hours for Th=333 K, kcor=10-15 m2, and kin=10-13 m2. It is about 7 hours and 6 hours for Th=333, 343 K, kcor=10-15 m2, and kin=10-13 m2. The drying time is more than 7 hours for Th=333 K, kcor=10-15 m2, and kin=10-13 m2. It is about 6 hours and 5 hours for Th=333, 343 K, kcor=10-15 m2, and kin=10-13 m2.Figure 10
Moisture curves at differentTh=323, 333, 343 K, kcor=10-15 m2, kin=10-13 m2.Figure 11
Moisture curves at differentTh=323, 333, 343 K, kcor=10-14 m2, kin=10-13 m2.Figure 12
Moisture curves at differentTh=323, 333, 343 K, no cortex, kin=10-13 m2.


## 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 Sl0=0.5, Kr=1000, Th=323 K, Pvb=1000 Pa, kin=10-13 m2, 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 (kin and kcor=10-13, 10-14, and 10-15 m2, resp.), with cortex (kin=10-13 and kcor=10-14 and 10-15 m2 resp.). The other simulation parameters are Sl0=0.5, Kr=1000, Th=323 K, and Pvb=1000 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 kcor=10-14 m2, that is lower 10-fold than inner medium. But when kcor=10-15 m2, the drying time is increased to 2-fold nearly. In order to compare, the no cortex, kcor=10-14 m2, and 10-15 m2 was shown. Because the numerical calculation is not convergence, the cortex kcor=10-16 m2, kin=10-13 m2 result not gotten. The reason is that the material parameter changer is too big.Figure 3
Moisture curves of no cortex (kin=10-13, kcor=10-13, 10-14, and 10-15 m2), (kcor=10-14, 10-15, kin=10-13 m2).


## 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 ofPvb=1000, 5000, and 9000 Pa are shown in Figures 4, 5, and 6 with different cortex intrinsic permeability. The other simulation parameters are Sl0=0.5, Kr=1000, and Th=323 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 differentPvb=1000, 5000, 9000 Pa, kcor=10-15 m2, kin=10-13 m2.Figure 5
Moisture curves of differentPvb=1000, 5000, 9000 Pa, kcor=10-14 m2, kin=10-13 m2.Figure 6
Moisture curves of differentPvb=1000, 5000, 9000 Pa, no cortex, kin=10-13 m2.


## 6.4. Effect of Initial Moisture Content

The effect of initial moisture content on the moisture curve is shown in Figures7, 8, and 9 for Sl0=0.5, 0.4, 0.3 with different cortex intrinsic permeability. To compare the results, moisture is represented by the moisture ratio M/M0. The other parameters are Kr=1000, Pvb=1000 Pa, and Th=323 K. The drying time is about 10 hours, 6.5 hours, and 6 hours for Sl0=0.5, 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 kcor=10-15 m2. The drying rate is little larger for Sl0=0.5, and the drying rate is almost the same for Sl0=0.4, 0.3.Figure 7
Moisture curves at differentSl0=0.5, 0.4, 0.3, kcor=10-15 m2, kin=10-13 m2.Figure 8
Moisture curves at differentSl0=0.5, 0.4, 0.3, kcor=10-14 m2, kin=10-13 m2.Figure 9
Moisture curves at differentSl0=0.5, 0.4, 0.3, no cortex, kin=10-13 m2.


## 6.5. Effect of Heat Temperature

The effect of heat temperature on moisture is shown in Figures10, 11, and 12 for the heat temperature Th=323, 333, and 343 K. The other parameters are Sl0=0.5, Kr=1000, and Pv=1000 Pa. The effect of heat temperature is obviously in all simulation. The drying time is more than 14 hours for Th=323 K, kcor=10-15 m2, and kin=10-13 m2. And it is about 10 hours and 8 hours for Th=333, 343 K, kcor=10-15 m2, and kin=10-13 m2. The drying time is more than 9 hours for Th=333 K, kcor=10-15 m2, and kin=10-13 m2. It is about 7 hours and 6 hours for Th=333, 343 K, kcor=10-15 m2, and kin=10-13 m2. The drying time is more than 7 hours for Th=333 K, kcor=10-15 m2, and kin=10-13 m2. It is about 6 hours and 5 hours for Th=333, 343 K, kcor=10-15 m2, and kin=10-13 m2.Figure 10
Moisture curves at differentTh=323, 333, 343 K, kcor=10-15 m2, kin=10-13 m2.Figure 11
Moisture curves at differentTh=323, 333, 343 K, kcor=10-14 m2, kin=10-13 m2.Figure 12
Moisture curves at differentTh=323, 333, 343 K, no cortex, kin=10-13 m2.


## 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 thekcor=10-14 m2, that is lower 10-fold than inner medium. But when kcor=10-15 m2, 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 kcor=10-15 m2. The drying rate is little larger for Sl0=0.5, and the drying rate is almost the same for Sl0=0.4, 0.3. The effect of heat temperature is obviously in all simulation.



---

*Source: 120736-2013-04-27.xml*