
Research Article

Mathematical Modelling of Thin Layer Drying of Salak

S.P. Ong
and
C.L. Law


ABSTRACT

Thin layer drying kinetics of salak (Salacca edulis) in hot air drying was investigated for temperature range of 50 to 80°C. Salak fruits were prepared in two forms; slices with membrane layer and slices without membrane layer on surface. A mathematical model was determined to describe the moisture reduction process during isothermal drying. Theoretical model of Fick’s second law, six semitheoretical models and one empirical model that commonly employed in food dehydration were evaluated. Results showed that the MidilliKucuk equation gave the best prediction to the drying kinetics evidenced by coefficient of determination, R^{2} ranging from 0.99760.9993. A new approach was used for the estimation of drying rate and the commencement point of second fallingrate period by using a semitheoretical model. Diffusivity coefficients of moisture transfer were found ranging from 4.19x10^{11} to 2.58x10^{10} m^{2} sec^{1}. It appears that moisture diffusion is the controlling mechanism of the thin layer drying and it is an increasing function of temperature. Activation energy was determined at 45.27 and 39.78 kJ mol^{1} for sample with and without surface membrane, respectively. Higher activation energy indicates that the membrane layer is a strong barrier to moisture movement across the sample surface. 




INTRODUCTION
Salak (Salacca edulis) or also known as snake fruit belongs to Arecaceae
family. This palm is native in Indonesia and Malaysia. Nevertheless, cultivation
of salak palm can also be found in Thailand and Philippines recently (AbadGarcía
et al., 2007). Salak fruit is a good source for dietary fibres and
carbohydrate (Lestari et al., 2003). This peculiar
fruit tastes sweet to subacid like combination of pineapple and jackfruit at
its ripe stage but immature fruit tastes sour and astringent. Despite the unique
smell and taste, salak contains valuable bioactive antioxidants such as vitamin
A, vitamin C and phenolic compounds (Setiawan et al.,
2001; Leong and Shui, 2002;
Shui and Leong, 2005; Leontowicz et al., 2007).
However, salak has a short shelf life of less than a week due to rapid ripening
and degradation of the bioactive ingredients.
Drying is one of the common methods used to preserve agricultural products.
Removal of moisture from produces may deactivate enzymes or microorganisms that
often cause undesired chemical reaction and lead to quality deterioration (Jayaraman
and Das Gupta, 2006). Drying kinetics is important in the analysis of moisture
migration process in a solid material. Many thermophysical properties and transport
properties that are usually integrated in a drying model can be determined from
the analysis of drying kinetics (MarinosKouris and Maroulis,
2006). Moreover, movement of moisture within a food material during drying
is a complex process with various diffusion mechanisms. Variation in composition
and complicated biopolymer structure of the food material may result in unique
drying characteristic for different food product (Zogzas
et al., 1996). Study on the drying kinetics may allow understanding
of the controlling mechanism during drying and hence influence of drying parameters
on dryer design and dried product quality can be determined.
Nemerous studies have been conducted to the drying kinetics and product quality
of foodstuffs that undergoing isothermal convective drying, for instance garlic,
pumpkin, grapes, apple and orange skin (Madamba et al.,
1996; Sawhney et al., 1999; Doymaz,
2007; Kaya et al., 2007). However, to date
there is no literature report on the drying kinetics of salak fruit which is a
good source for dietary antioxidant but unexploited. On the other hand, several
mathematical models have been proposed for the evaluation of moisture transport
behavior by either applying theoretical (Fick’s second law of diffusion),
semitheoretical and empirical modelling (Crank, 1975;
Jayas et al., 1991; Kabganian et al., 2002;
Demir et al., 2007; Midilli
et al., 2002; Doymaz, 2005). Theoretical model
can be used to get an insight on the mechanism of moisture transfer during the
fallingrate period in drying. However, simplified assumptions such as constant
diffusivity and onedimension liquid diffusion sometimes resulted in inadequate
prediction to the moisture distribution. Meanwhile, a semitheoretical model may
not able to explain the exact mechanism of moisture transport but it often gives
good estimation by incorporating lump values of other effects into the model parameters.
The objective of this study is to investigate the kinetics of thin layer drying of salak and to determine mathematical model for the moisture distribution profile of the convective hot air drying. This study also aims to extend the use of semitheoretical model in the estimation of drying rate at different drying temperature and determination of commencement point of second falling rate period. MATERIALS AND METHODS Sample preparation: Salak fruits (Salacca edulis) were purchased from a local fruit supplier. Fruits with similar size (ca. 57x5 cm) and skin colour were selected for all experiments to ensure consistency in the samples used. Fresh fruits were kept in refrigerator at about 10°C with storage period not longer than seven days. Salak slices (3x2x2 mm) were prepared in two forms, with surface membrane and without surface membrane. Sample without surface membrane was prepared by removing the membrane layer manually from the fruit surface before slicing. Drying experiments: Convective hot air drying experiments were performed by using an oven with forced air circulation (Memmert, UFP500, Germany). Heating elements were integrated evenly around the four inner walls to create uniform temperature distribution inside the chamber. Fresh air was preheated in a prewarming chamber and continuously added to the air inside the chamber. Meanwhile, exhaust air was discharged through an opening situated in the middle of back wall. Four levels of hot air temperature (50, 60, 70 and 80°C) were used. Samples (about 100 g) were arranged on flat perforated trays (65x35 cm with 1 cm diameter openings) and placed in the middle section of the oven. The weight of the sample was recorded periodically during the drying with 15 min interval in the first 2 h, 30 min interval in the subsequent next 2 h and hourly thereafter by using a digital balance (Adventurer OHAUS, AR3130, USA). Experiment was terminated when constant weight was obtained in three consecutive measurements. The moisture content obtained at this stage was marked as the Equilibrium Moisture Content (EMC). Dry Matter (DW) weight of the material was determined by drying the sample in oven at 105°C for 24 h.
Mathematical modelling
Moisture ratio: Moisture content was determined as a dimensionless
parameter denoted as Moisture Ratio (MR). Equation 1 relates
the sample moisture content in real time (X_{t}) to the initial moisture
content (X_{o}) and equilibrium moisture content (X_{e}) (Akpinar
et al., 2003):
Theoretical model: Fick’s second law was solved analytically for
diffusion in a plane sheet by assuming onedimension diffusion and constant
diffusivity (Crank, 1975). Ten terms (n = 10) from the
Fick’s series in Eq. 2 was used to estimate the moisture
content by conducting nonlinear regression analysis using Microsoft Excel SOLVER
tool (Microsoft Office Professional 2003, USA). Effective diffusivity constant
(D_{eff}) and sample thickness (L) are the model parameters.
Semitheoretical and empirical model: Semitheoretical and empirical models that commonly applied for fruits and vegetables drying were adopted from literatures as shown in Table 1. Three statistical parameters, coefficient of determination (R^{2}), Chisquare (χ^{2}) and Root Mean Square Error (RMSE), were used to assess the goodness of fit of the models to experiment data.
Drying rate and fallingrate period: Drying rate at a given temperature
was estimated by using first derivation of MidilliKucuk model equation while
gradient of the drying rate was calculated by conducting second derivation to
the same model.
Table 1: 
Semitheoretical and empirical models 

Effective diffusivity and activation energy: Effective diffusivity
constant (D_{eff}) can be determined from the linearized first term
from the Fick’s infinite equation as shown in Eq. 3 by
assuming that constant diffusivity applied to the internal moisture movement
during first fallingrate period.
Natural logarithm of MR was plotted versus t and D_{eff} can be determined
from the slope. Relationship between D_{eff} and temperature can be
evaluated with Arrhenius equation as shown in Eq. 4.
RESULTS AND DISCUSSION
Drying kinetics: Moisture content of samples at four levels of temperature
(50, 60, 70 and 80°C) were monitored over the drying period. Initial moisture
content of the salak fruits were recorded at 4.24±0.39 kg water kg DW^{1}.
Figure 1a and b show the moisture ratio
versus time for salak slices with and without surface membrane, respectively.
Both graphs show exponential trend for the drying curves and it can be observed
that the samples reached EMC in a shorter time at higher temperature. Many literatures
showed similar results when conducting hot air drying to agricultural products
such as red chillies, pineapple and potatoes (Arora et
al., 2006; Nicoleti et al., 2001; Rossello
et al., 1992). Higher drying temperature would result in higher drying
kinetics. However, membrane layer on the fruit surface is a barrier to moisture
transfer. Higher moisture content was found in salak with surface membrane when
compared with those without membrane at a given time period. Nevertheless, the
effect of surface barrier was less pronounce when drying process was conducted
at higher temperature. Moisture contents for both samples (with and without
membrane) were found close to each other for the hot air drying at temperature
80°C.
Mathematical modelling: Mathematical modelling of the drying curves
under various temperatures was conducted by using nonlinear regression analysis
coupled with generalized reduced gradient algorithm (GRG2). All the model constants
were calculated based on the iterative method and the statistical parameters
for the estimation are shown in Table 2. Based on the obtained
values for the three statistical parameters, it can be seen that MidilliKucuk
equation is the best model to describe the drying kinetics of salak slices at
all tested temperatures with R^{2}, χ^{2} and RMSE values
ranging from 0.99760.9993, 0.0001 0.0003 and 0.00680.0128, respectively.
 Fig. 1: 
Moisture ratio of salak slices at different temperature (a)
with surface membrane and (b) without surface membrane 
In agreement with other literatures, MidilliKucuk model often gave better
prediction on the moisture distribution profile among other semitheoretical
models (Ait Mohamed et al., 2005;
HacIhafIzoglu et al., 2008; Karaaslan and Tunçer,
2008). On the other hand, Fick’s equation (theoretical) resulted in
lower value of R^{2} as compared with semitheoretical models, while
the Wang and Singh model (empirical) presented the lowest R^{2} value
among the theoretical and semitheoretical models.
Drying rate and fallingrate period: High value of R^{2} and
low value of χ^{2} and RMSE obtained from the regression analyses
indicate that the MidilliKucuk model fitted well with the experiments data.
Thus, it can be used to represent the drying curve for salak slices under isothermal
drying and the drying rate can be calculated by applying first derivative to
the MidilliKucuk equation. Drying rate curves for salak slices with and without
surface membrane are shown in Fig. 2a and b,
respectively. It was observed hat there was no constantrate period in all drying
rate curves.
Table 2: 
Results of statistical parameters estimated from regression
analyses for temperature of 5080°C 

 Fig. 2: 
Drying rate curves for salak slices (a) with surface membrane
and (b) without surface membrane 
All drying experiments comprised of only fallingrate period. The rate of drying
fell steadily with time. It appears that higher drying rate was obtained at
high temperature and sample without surface membrane attained higher drying
rate when compared with those with membrane layer at a given temperature.
Many literatures suggested that moisture diffusion in solid material during
fallingrate period is diffusioncontrolled and moisture diffusion in first
fallingrate period is more prominent and significant as compared with diffusion
in second fallingrate period. This is because second fallingrate period is
the period of unsaturated surface evaporation and usually the amount of moisture
removed is relatively small but the time required may be long. Hence, moisture
diffusivity can be estimated from the first fallingrate period and the calculated
effective diffusivities are lump values over the range of moisture content (Geankoplis,
1993). Nevertheless, in most cases, the commencement of second fallingrate
period is difficult to be estimated from the drying rate curve due to unobvious
curve change especially when data points are closed to each other at low moisture
content.
Thus, plots for drying rate gradient against drying time were constructed by
applying second derivation to the MidilliKucuk model. Figure
3a and b show the curves for the dying rate gradient with
and without membrane layer, respectively. Overall, drying rate gradient was
found increased with time before reaching a constant value and decreased with
time thereafter. The obtained constant value is an indication of first fallingrate
period due to the characteristic of the first fallingrate curve which is usually
a straight line with slope.
 Fig. 3: 
Drying rate gradient for salak slices (a) with surface membrane
and (b) without surface membrane 
It appears that time period for the occurrence of the first fallingrate is
dependent on the temperature. For instance, the beginning points for the second
fallingrate period (also the end points for first fallingrate period) were
found at 120, 75, 60 and 45 min at temperature of 50, 60, 70 and 80°C, respectively,
for the drying of salak with membrane layer. It seems that at higher temperature
the first fallingrate period will be shorten.
Effective diffusivity and activation energy: Plots of lnMR versus t
for the four drying temperatures were constructed by using only the data fell
within the first fallingrate period and the results are shown in Fig.
4a and b. The constant diffusivity (D_{o}) and
activation energy (E_{a}) were determined from the yintercept and slope
of the linearized Arrhenius equation curve, respectively. Table
3 shows the summary of the D_{eff}, D_{o }and E_{a}.
The estimated D_{eff} ranging from 4.19x10^{11} to 2.58x10^{10}
m sec^{1} and the obtained values are fall within the range of values
reported in other literatures (Madamba et al., 1996;
Nicoleti et al., 2001; Biju Cletus and Carson,
2008; Garau et al., 2006).
 Fig. 4: 
Plot of lnMR versus t for salak slices (a) with surface membrane
and (b) without surface membrane 
Table 3: 
Effective diffusivity, constant diffusivity and activation
energy for salak with or without membrane 

Apparently, the effective diffusivity increased with drying temperature. On
the other hand, activation energy for sample with surface membrane was found
higher than those without surface membrane. This could be resulted by the existence
of the membrane on the sample surface that had become the resistance for the
moisture movement. Thus higher energy required to overcome the resistance.
CONCLUSION
Elevated air temperature resulted in higher drying kinetics. Effective diffusivity
constant was found increased with temperature and correlated well with the Arrhenius
equation. Higher activation energy showed that the surface membrane was a strong
barrier to moisture diffusion from internal to surrounding air. MidilliKucuk
model gave better prediction to the moisture distribution profile as compared
to other semitheoretical models.
ACKNOWLEDGMENTS The authors are grateful to the Ministry of Science, Technology and Innovation for the financial support through eScience funding and support from the University of Nottingham, Malaysia Campus.

REFERENCES 
AbadGarcia, B., L.A. Berrueta, D.M. LopezMarquez, I. CrespoFerrer, B. Gallo and F. Vicente, 2007. Optimization and validation of a methodology based on solvent extraction and liquid chromatography for the simultaneous determination of several polyphenolic families in fruit juices. J. Chromatogr. A, 1154: 8796. CrossRef  Direct Link 
Ait Mohamed, L., M. Kouhila, A. Jamali, S. Lahsasni, N. Kechaou and M. Mahrouz, 2005. Single layer solar drying behaviour of Citrus aurantium leaves under forced convection. Energy Conver. Manage., 46: 14731483. CrossRef  Direct Link 
Akpinar, E.K., Y. Bicer and C. Yildiz, 2003. Thin layer drying of red pepper. J. Food Eng., 59: 99104. CrossRef  Direct Link 
Arora, S., S. Bharti and V.K. Sehgal, 2006. Convective drying kinetics of red chillies. Dry. Technol., 24: 189193. Direct Link 
Biju Cletus, A. and J.K. Carson, 2008. Drying curves and apparent diffusivity of New Zealand chestnut variety 1015. J. Food Eng., 85: 381386. Direct Link 
Crank, J., 1975. The Mathematics of Diffusion. 2nd Edn., Clarendon Press, Oxford, UK.
Demir, V., T. Gunhan and A.K. Yagcioglu, 2007. Mathematical modelling of convection drying of green table olives. Biosyst. Eng., 98: 4753. CrossRef  Direct Link 
Doymaz, I., 2005. Drying characteristics and kinetics of Okra. J. Food Eng., 69: 275279. CrossRef  Direct Link 
Doymaz, I., 2007. The kinetics of forced convective airdrying of pumpkin slices. J. Food Eng., 79: 243248. CrossRef  Direct Link 
Garau, M.C., S. Simal, A. Femenia and C. Rosselló, 2006. Drying of orange skin: Drying kinetics modelling and functional properties. J. Food Eng., 75: 288295. Direct Link 
Geankoplis, C.J., 1993. Transport Processes and Unit Operations. 3rd Edn., PrenticeHall, Inc., Englewood Cliffs, New Jersey.
HacIhafIzoglu, O., A. Cihan and K. Kahveci, 2008. Mathematical modelling of drying of thin layer rough rice. Food Bioproducts Process., 86: 268275. Direct Link 
Jayaraman, K.S. and D.K. Das Gupta, 2006. Drying of Fruits and Vegetables. In: Handbook of Industrial Drying, Mujumdar, A.S. (Ed.). CRC Press, Boca Raton, FL.
Jayas, D.S., S. Cenkowski, S. Pabis and W.E. Muir, 1991. Review of thinlayer drying and wetting equations. Dry. Technol., 9: 551588. CrossRef  Direct Link 
Kabganian, R., D.J. Carrier and S. Sokhansanj, 2002. Physical characteristics and drying rate of echinacea root. Dry. Technol., 20: 637649. Direct Link 
Karaaslan, S.N. and I.K. Tunçer, 2008. Development of a drying model for combined microwavefanassisted convection drying of spinach. Biosyst. Eng., 100: 4452. Direct Link 
Kaya, A., O. Aydin and C. Demirtas, 2007. Drying kinetics of red delicious apple. Biosyst. Eng., 96: 517524. CrossRef  Direct Link 
Leong, L.P. and G. Shui, 2002. An investigation of antioxidant capacity of fruits in Singapore markets. Food Chem., 76: 6975. Direct Link 
Leontowicz, M., H. Leontowicz, J. Drzewiecki, Z. Jastrzebski and R. Haruenkit et al., 2007. Two exotic fruits positively affect rat's plasma composition. Food Chem., 102: 192200. Direct Link 
Lestari, R., S.H. Keil and G. Ebert, 2003. Variation in fruit quality of different salak genotypes (Salacca Zalacca (Gaert.) Voss) from Indonesia. Deutscher TropentagTechnological and Institutional Innovations for Sustainable Rural Development, Göttingen, Germany.
Madamba, P.S., R.H. Driscoll and K.A. Buckle, 1996. The thinlayer drying characteristics of garlic slices. J. Food Eng., 29: 7597. Direct Link 
MarinosKouris, D. and Z.B. Maroulis, 2006. Transport Properties in the Drying of Solids. In: Handbook of Industrial Drying, Mujumdar, A.S. (Ed.). CRC Press, Boca Raton, FL.
Midilli, A., H. Kucuk and Z. Yapar, 2002. A new model for singlelayer drying. Dry. Technol., 20: 15031513. Direct Link 
Nicoleti, J.F., J. TelisRomero and V.R.N. Telis, 2001. Airdrying of fresh and osmotically pretreated pineapple slices: Fixed air temperature versus fixed slice temperature drying kinetics. Dry. Technol. Int. J., 19: 21752191. Direct Link 
Rossello, C., J. Canellas, S. Simal and A. Berna, 1992. Simple mathematical model to predict the drying rates of potatoes. J. Agric. Food Chem., 40: 23742378. Direct Link 
Sawhney, R.L., D.R. Pangavhane and P.N. Sarsavadia, 1999. Drying kinetics of single layer Thompson seedless grapes under heated ambient air conditions. Dry. Technol., 17: 215236. Direct Link 
Setiawan, B., A. Sulaeman, D.W. Giraud and J.A. Driskell, 2001. Carotenoid content of selected indonesian fruits. J. Food Comp. Anal., 14: 169176. Direct Link 
Shui, G. and L.P. Leong, 2005. Screening and identification of antioxidants in biological samples using highperformance liquid chromatographymass spectrometry and its application on salacca edulis reinw. J. Agric. Food Chem., 53: 880886. Direct Link 
Zogzas, N.P., Z.B. Maroulis and D. MarinosKouris, 1996. Moisture diffusivity data compilation in foodstuffs. Drying Technol. Int. J., 14: 22252253. Direct Link 



