Modeling Some Drying Characteristics of Garlic She
Modeling Some Drying Characteristics of Garlic She
Modeling Some Drying Characteristics of Garlic She
net/publication/266508199
CITATIONS READS
17 2,151
3 authors:
Behdad Shadidi
Buali Sina University
13 PUBLICATIONS 202 CITATIONS
SEE PROFILE
All content following this page was uploaded by Kamran Salari on 01 September 2016.
Abstract
Amiri Chayjan R., Salari K., Shadidi B., 2012. Modeling some drying characteristics of garlic sheets under
semi fluidized and fluidized bed conditions. Res. Agr. Eng., 58: 73–82.
Thin layer drying properties of high moisture garlic sheets under semi fluidized and fluidized bed conditions with high
initial moisture content (about 154.26% d.b.) were studied. Air temperatures of 50, 60, 70 and 80°C were applied to
garlic samples. Among the applied models, Page model was the best to predict the thin layer drying behavior of garlic
sheets. Using this model, correlation coefficient (R2) was high for all drying cases. The computed values of Deff were
between 3.38 × 10–10 and 2.54 × 10–9 m2/s during the falling rate drying. Values of Deff for garlic sheets were also in-
creased with increasing in input air temperature. Activation energy values were varied between 51.32 and 60.58 kJ/mol
for 50 to 80°C, respectively. The specific energy consumption (SEC) for garlic specimens was placed in the range of
0.316 × 106 and 0.979 × 106 kJ/kg from 50 to 80°C, respectively. An increase in air temperature caused decrease in
SEC value. Application of semi fluidized bed convective drying with temperature between 50 and 60°C was suitable to
produce dried garlic.
Garlic (Allium sativum L.) is an important crop in of the garlic will quickly appeared after dormancy.
the world. Due to its therapeutic properties it has When sprouting is started, the nutritional storage
been cultivated in many countries. Garlic is also in the garlic will be decreased and due to wither-
usually used as a flavoring agent; it may be used in ing the stem and trunk, garlic quality is obviously
the shape of powder or granule as a valuable condi- reduced. Drying is one of the best methods for pre-
ment for foods. serving the garlic quality.
When the garlic bulb is cut or split, pungency It is proved that the quality of dried food products
of flavor is diffused, while fresh garlic bulb has no are strongly affected by applied drying methods and
distinct pungency. Process of drying garlic is usu- various physical, chemical and biological changes
ally evaluated based on many indices such as: dry- may be created in the food material (Krokida et
ing time, effective moisture diffusivity, activation al. 2000). In other words, some properties of food-
energy, specific consumption of energy and color. stuffs such as color, structure, aroma compounds
Applying of proper drying method is a key to pro- and nutritional substance were changed or deterio-
ducing high quality dried garlic with a higher dif- rated. These changes may tend to reduce the prod-
fusivity, lower drying time and energy and better uct quality (Pezzutti, Crapiste 1997).
appearance as similar to fresh garlic as possible. Hot-air drying of garlic slices in a common fixed
If the environmental conditions are suitable and bed method is unfortunately not suitable due to a
the rest period passed about 60–80 days, the sprout significant decrease in the quality of dried product
73
Vol. 58, 2012, No. 2: 73–82 Res. Agr. Eng.
related to the fresh one. Applying high temperature high moisture garlic sheets and to compute the effec-
(about 60°C) in a fixed bed drying causes an in- tive moisture diffusivity, activation energy and specif-
crease in drying period, energy consumption, color ic energy consumption of high moisture garlic sheets
degradation and mass transfer. during falling rate of semi fluidized and fluidized bed
In recent years, fluidized bed drying was inves- thin layer drying method and their relation to input
tigated and utilized as a new method for obtaining parameters such as air temperature and velocity.
dried foodstuffs with high quality (Poomsa-Ad et
al. 2002; Cubillos, Reyes 2003; Amiri Chayjan
et al. 2009; Gazor 2009). Fluidized bed drying is MATERIAL AND METHODS
rapid and more uniform compared with fixed bed
drying (Soponronnarit et al. 1997). Fluidized Fluidized bed drying condition
bed drying was employed for drying of some ag-
ricultural grain products such as: broad beans In order to apply fluidized bed condition in dry-
(Hashemi et al. 2009), milky mushroom (Arumu- ing of garlic sheets, maximum pressure drop against
ganathan et al. 2009), rough rice (Amiri Chay- air flow velocity in fluidization curve was deter-
jan et al. 2009), green beans (Souraki, Mowla mined (Kunii, Levenspiel 1991). A digital appa-
2007) and corn (Soponronnarit et al. 1997). ratus (Standard ST-8897, Standard Instruments
Suspending of grain particles in air flow is known Co., Ltd., Hong Kong, China) utilized for recording
as fluidization. At the beginning of air passing both static pressure drop (kPa) and outlet air ve-
through grain bed, a fixed bed will be created. locity (m/s). About 25 g garlic sample was used in
With gradually increasing in air flow rate, a bed fluidization experiments. First experimental point
of fluffy material is obtained, namely minimum was determined at air velocity about 2.8 m/s as a
fluidized bed (semi fluidized bed); afterward with semi fluidized bed condition. Two other points
another increase in air flow rate, bubbling fluid- were in fluidized bed condition with air velocities
ized bed and transportation phenomenon would of 4.26 and 5.4 m/s, respectively. These points and
be observed. At the semi fluidized bed, maximum related fluidization curve were shown in Fig. 1. The
value of pressure drop is obtained and particles net static pressure of garlic bed for points of A, B
weight is equal to frictional force between bed and C were obtained 0.048, 0.037 and 0.026 kPa.
particles (Kunii, Levenspiel 1991; Soponron- As can be observed, net static pressure for point A
narit et al. 1997). (Fig. 1) is maximized and is defined as semi fluid-
Mathematical modeling of garlic drying using ized condition.
defined models can precisely predict the drying
kinetics in a drying system (Kaleta, Górnicki
2010). Effective moisture diffusivity is affected by Experimental setup
air temperature and velocity (Senadeera et al.
2000, 2003). All effective parameters on the mass Fresh garlic was purchased from a local farm
transfer phenomenon in drying process are repre- in September 2010. Garlic samples were kept in a
sented by effective moisture diffusivity (Hashemi
et al. 2009). Activation energy is important for es- 0.06
timation of minimum energy requirement for car- 0.05
Pressure drop (kPa)
74
Res. Agr. Eng. Vol. 58, 2012, No. 2: 73–82
refrigerator at 3 ± 1°C. Ambient air temperature The slope (B) is determined by plotting time
and relative air humidity changed from 26 to 32°C against ln(MR) as follows:
and from 22 to 31%, respectively. During the dry- π 2 Deff
ing experiments, inlet and outlet temperatures of B= (4)
4 L2
the drying chamber and the ambient air relative
humidity and temperature were measured. An ex- Activation energy was obtained using an Arrhe-
perimental fluidized bed dryer utilized to perform nius type equation (Babalis, Belessiotis 2004):
the experiments. Applied air temperatures in the E
experiments were 50, 60, 70, and 80°C. In total, Deff = D0 exp − a (5)
RT
12 experiments were conducted. Initial and final
moisture contents determination of garlic samples For calculating Ea, Eq. (5) can be linearized as fol-
was conducted using gravimetric method at 103°C lows:
(Vázquez et al. 1999). Drying process was started
E 1
at initial moisture content of about 154.26% (d.b.) ln( Deff ) = ln( D0 ) − a (6)
and finished at final moisture content of about 5% R T
(d.b.). where:
Ea – activation energy (kJ/mol)
R – universal gas constant (8.3143 kJ/mol.K)
Modeling T – absolute air temperature (K)
D0 – pre-exponential factor of the equation (m2/s)
Fick’s second law for moisture diffusion in infi-
nite slab was employed for modeling. Due to minor After plotting of 1/T against ln(Deff) using Eq. (6),
diameter of garlic sheets, that were much smaller six fitted models were obtained as straight lines
than their major diameter, we therefore assumed with the slope of B1.
that the garlic sheets are similar to infinite slab. E
After solution of Fickian equation for an infinite B1 = a
R (7)
slab, the following equation was obtained and can
describe effective moisture diffusivity changes in Simplified form of Eq. (3) for all models can be
garlic sheets: written as exponential equations as follows:
M − Me 8 ∞
1 − Deff (2n − 1) 2 π 2t M − Me
MR = =
M 0 − M e π2
∑ (2n − 1) 2
exp
4 L2
(1) MR =
M0 − Me
= a exp(−kt ) (8)
n =1
where: They are mathematical models for prediction of
MR – moisture ratio (decimal) garlic drying kinetics in thin layer mode. Some of
M – moisture content at any time (kgwater/kgdry mater) these models were employed in this study (Table 1).
M0 – initial moisture content (kgwater/kgdry mater) The value of equilibrium moisture content (Me) is
n = 1, 2, 3, ... – number of terms taken into consideration very small compared to M and M0 (Hassan-Beygi
t – time of drying (s) et al. 2009). Thus (M – Me)/(M0 – Me) was changed
Deff – effective moisture diffusivity (m2/s) to M/M0. Therefore all models in Table 1 and the
L – thickness of garlic sheet (m) basic Eq. (2) can be written as follows:
M
For a long drying period, only the first term of MR = (9)
M0
Eq. (1) can be considered. This simplification has
no negative effect on the final answer (Ramesh et Specific energy consumption (SEC) was calculat-
al. 2001): ed using the following model (Zhang et al. 2002):
M − Me 8 D π 2t
ln (MR) = ln = ln 2 − eff 2 (2) SEC =
(C Pa + C Pv ha )Qt (Tin − Tam )
(10)
M0 − Me π 4L mvVh
The simplified form of MR is as follows: where:
SEC – specific energy consumption (kJ/kg)
8 π 2 Deff t CPv,CPa – specific heat capacity of vapor and air, respec-
MR = 2 exp 2
(3)
π 4L tively, (1,004.16 and 1,828.8 J/kg°C)
75
Vol. 58, 2012, No. 2: 73–82 Res. Agr. Eng.
Table 1. Thin layer drying models used in thin layer modeling of high moisture garlic
1
a, b, c, k, k0, k1, n – drying constants
76
Table 2. Comparative indices of statistical model parameters for thin layer drying of high moisture garlic sheets
R2 χ2 RMSE
Air temperature
Model
(°C)
2.8 m/s 4.26 m/s 5.4 m/s 2.8 m/s 4.26 m/s 5.4 m/s 2.8 m/s 4.26 m/s 5.4 m/s
Res. Agr. Eng.
77
Vol. 58, 2012, No. 2: 73–82
Vol. 58, 2012, No. 2: 73–82 Res. Agr. Eng.
Table 3. Coefficients of Page model for prediction of Effect of bed condition on Deff
kinetic drying of high moisture garlic sheets
Variations of Deff against air temperature at dif-
Bed condition Temperature (°C) k n
ferent bed conditions are plotted in Fig. 3. Three
50 3.8534 1.1019 exponential models were fitted to the obtained val-
Semi fluidized 60 4.6402 1.0508 ues of Deff. These models as well as related R2 values
bed (2.8 m/s) 70 9.7840 1.2351 are represented in Tables 5 and 6. Results showed
80 28.748 1.4898
that the maximum value of Deff was obtained at the
highest air temperature level; they also indicated
50 3.2959 1.0781 that applying fluidized bed condition caused a de-
Fluidized bed 60 4.6924 1.0972 crease in Deff at upper air temperatures. Drying air
(4.26 m/s) 70 9.7901 1.2573 contact with garlic sheets at semi fluidized bed was
80 24.009 1.4791 the most effective due to its highest values of Deff.
Also at lower air temperatures, no significant dif-
50 5.3186 1.2938
ference was observed between Deff values of bed
Fluidized bed 60 5.5939 1.0706 conditions (air velocity levels), as applying mini-
(5.4 m/s) 70 10.997 1.2182 mum fluidized bed condition with lower air veloc-
80 11.681 1.1118 ity was even more suitable. Four quadratic models
were applied to fit on the obtained moisture diffu-
for calculation of Deff values. These values for all bed sivity. Applied models and related R2 values for four
conditions and temperatures are presented in Table 4. air temperatures are presented in Table 6.
Maximum value of Deff (2.54 × 10–9 m2/s) among all
drying cases was obtained for semi fluidized condi-
tion with air velocity of 2.8 m/s and air temperature Computation of activation energy
of 80°C. This is due to the most effective contact be-
tween garlic sheets and air flow. Also minimum value Fig. 4 shows the curve of ln(Deff) against 1/T for
of Deff (3.38 × 10–10 m2/s) belonged to fluidized condi- drying of garlic sheets in falling rate period. Ac-
tion with air velocity of 5.4 m/s and air temperature of tivation energy (Ea) was calculated using Eq. (6).
50°C. This is due to the low energy rate transferred to The obtained values of Ea for all bed conditions as
the garlic sheet and also low effective contact in tur- well as R2 values are represented in Table 7. Val-
bulent flow of fluidized bed. ues of Ea for food and agricultural crops generally
Values of Deff for high moisture garlic sheets are varied between 12.7 and 130 kJ/mol (Aghbashlo
greatly affected by drying air temperature. As can et al. 2008). Minimum and maximum values of
be seen in Table 4, an increase of 10°C in input Ea for high moisture garlic sheets varied between
drying temperature, cause about twice increase in 51.32 and 60.58 kJ/mol for all bed conditions. Two
Deff values. This pattern was observed in previous modes of moisture in agricultural materials are
studies, such as: peaches (Kingsly et al. 2007) and surface and chemical. Because most of the water in
plums (Goyal et al. 2007). high moisture garlic in falling rate period is in the
Time (h)
Time (hour)
0.0
0 0.5 1 1.5 2 2.5 3 3.5
-0.5 5.4 m/s
-1.0
ln(MR)
ln(MR)
-1.5 80 °C
70 °C
60 °C
-2.0 50 °C
78
Res. Agr. Eng. Vol. 58, 2012, No. 2: 73–82
Table 4. Effective moisture diffusivity and correlation coefficient for three bed conditions at different temperatures
form of chemical absorption, relatively more en- was increased as air temperature was decreased.
ergy is consumed to exhaust water and undesirable Air velocity has a direct relation with SEC, as an in-
change in chemical properties is noticeable in this crease in air velocity caused an intensive increase in
period (Amiri Chayjan et al. 2009). SEC. The minimum value of SEC (0.316 × 106 kJ/kg)
Maximum value of Ea belonged to the minimum was obtained for minimum fluidized bed condition
fluidized bed condition with air velocity 2.8 m/s. and drying air temperatures 90°C. The maximum
With increasing air velocity, the activation energy value of SEC (0.979 × 106 kJ/kg) was obtained for
was decreased. A linear equation is fitted to the cal- fluidized bed condition with air velocity of 5.4 m/s
culated data of Ea against air velocity as follows: and drying air temperature 50°C. Results indicated
that an increase in drying temperature affect SEC
Ea = –3.5451v – 70.673 R2 = 0.9950 (14) inversely. With increasing air velocity, effective con-
tact between air and garlic sheets was contrarily in-
creased but an output energy loss was increased and
Specific energy consumption SEC was therefore increased. Three linear models
were applied to fit the SEC data in semi fluidized bed
Computed SEC values for drying of samples are and fluidized bed conditions as follows:
presented in Fig. 5. It can be observed that the SEC SEC = 5,298.1T + 106 R2 = 0.9941 (2.8 m/s) (15)
Table 5. Fitted exponential models to Deff values for different bed conditions
Table 6. Fitted power models to Deff value for different air temperatures
Table 7. Activation energy and related correlation coefficient for different bed conditions
Bed condition Semi fluidized bed (V = 2.8 m/s) Fluidized bed (V = 4.26 m/s) Fluidized bed (V = 5.4 m/s)
7,286.3 6,795.4 6,172.1
Equation ln(Deff) = 0.76 – ln(Deff) = 0.766 – ln(Deff) = –2.7285 –
T T T
Ea (kJ/mol) 60.58 56.50 51.32
R2 0.9846 0.9888 0.9976
79
50 3.72E-10 50 3.63E-10 0.8661 3.9E-10
60 6.12E-10 60 5.97E-10 1.3433 6.05E-10
70 1.17E-09 70 1.07E-09 2.1641 9.75E-10
80 2.54E-09 80 2.19E-09 3.8081 1.71E-09
1.50E-09
2.8 1.00E-09
2.8 4.26 4.26 5.4 5.4
FFP FFP FFP FFP FFP FFP
0.003096 -21.7112 0.003096 -21.7377 0.003096 -21.808
0.003003 5.00E-10
-21.2135 0.003003 -21.2398 0.003003 -21.296
0.002915 -20.5686 0.002915 -20.6523 0.002915 -20.749
0.002833 -19.7894 0.002833 -19.9417 0.002833 -20.1839
0.00E+00
50 60 70 80
Air temperature (°C)
1/T
1/T
-19.5
0.0028 0.00285 0.0029 0.00295 0.003 0.00305 0.0031 0.00315
-20
4.26 m/s
ln(Deff
ln(D
5.4 m/s
-21
-21.5
Fig. 4. ln(Deff) versus 1/T at different levels
of air velocities for thin-layer drying of high
-22 moisture garlic
SEC = –6,085.7T + 106 R2 = 0.9850 (4.26 m/s) (16) as the best method applied in this study. Colour
analysis was accomplished on dried garlic samples
SEC = –6,596.6T + 106 R2 = 0.9899 (5.4 m/s) (17) from semi fluidized method. Results showed that
the garlic sheets dried at 50 and 60°C have accept-
able colour. But its color at 70 and 80°C changed to
Comparison between dried garlic sheets brown or black and transformation also occurred.
If a fluidized bed industrial dryer is designed for
With regard to calculated indices such as Deff, Ea garlic drying, then quality of final dry product
and SEC, semi fluidized bed condition was selected should be considered.
7.0E+05
6.0E+05
5.0E+05
4.0E+05
3.0E+05 Fig. 5. Specific energy consumption for
2.0E+05 thin layer drying of high moisture garlic
40 50 60 70 80 90 sheets at different air temperatures of
Temperature (°C) semi fluidized and fluidized bed drying
80
Res. Agr. Eng. Vol. 58, 2012, No. 2: 73–82
81
Vol. 58, 2012, No. 2: 73–82 Res. Agr. Eng.
Ramesh M.N., Wolf W., Tevini D., Jung G., 2001. Influence bed dryer with energy carrier: modeling with and without
of processing parameters on the drying of spice paprika. shrinkage. Journal of Food Engineering, 88: 9–19.
Journal of Food Engineering, 49: 63–72. Vázquez G., Chenlo F., Moreira R., Costoyas A., 1999.
Senadeera W., Bhandari B.G., Wijesinghe B., 2000. Physi- The dehydration of garlic: desorption isotherms and model-
cal properties and fluidization behavior of fresh green bean ling of drying kinetics. Drying Technology, 17: 1095–1108.
particulates during fluidized bed drying. Transactions of Wang C.Y., Singh R.P., 1978. Use of variable equilibrium
the IChemE, 78(C): 43–47. moisture content in modelling rice drying. Transactions
Senadeera W., Bhandari B.G., Wijesinghe B., 2003. of the ASAE, 11: 668–672
Influence of shapes of selected vegetables materials on Yaldiz O., Ertekin C., Uzun H.I., 2001. Mathematical
drying kinetics during fluidize bed drying. Journal of Food modelling of thin layer solar drying of sultana grapes.
Engineering, 58: 277–283. Energy, 26: 457–465.
Sharaf-Elden Y.I., Blaisdell J.L., Hamdy M.Y., 1980. A Zhang Q., Yang S.X., Mittal G.S., YI S., 2002. Prediction
model for ear corn drying. Transactions of the ASAE, 23: of performance indices and optimal parameters of rough
1261–1265. rice drying using neural network. Biosystems Engineering,
Soponronnarit S., Pongtornkulpanich A., Prachaya- 83: 281–290.
warakorn S., 1997. Drying characteristics of corn in
fluidized bed dryer. Drying Technology, 15: 1603–1615. Received for publication June 7, 2011
Souraki B.A., Mowla D., 2007. Axial and radial moisture Accepted after corrections August 23, 2011
diffusivity in cylindrical fresh green beans in a fluidized
Corresponding author:
82