alexa Mathematical Modelling of Thin Layer Drying Kinetics of Biodegradable Pellets | Open Access Journals
ISSN: 2157-7110
Journal of Food Processing & Technology
Like us on:
Make the best use of Scientific Research and information from our 700+ peer reviewed, Open Access Journals that operates with the help of 50,000+ Editorial Board Members and esteemed reviewers and 1000+ Scientific associations in Medical, Clinical, Pharmaceutical, Engineering, Technology and Management Fields.
Meet Inspiring Speakers and Experts at our 3000+ Global Conferenceseries Events with over 600+ Conferences, 1200+ Symposiums and 1200+ Workshops on
Medical, Pharma, Engineering, Science, Technology and Business

Mathematical Modelling of Thin Layer Drying Kinetics of Biodegradable Pellets

Kulsum Jan, Riar CS and Saxena DC*

Department of Food Engineering and Technology, Sant Longowal Institute of Engineering & Technology, Longowal – 148 106, Punjab, India

*Corresponding Author:
Saxena DC
Professor and Head, Department of Food Engineering and Technology
Sant Longowal Institute of Engineering & Technology
Longowal – 148 106, Punjab, India
Tel: +91-98-156-088-59
E-mail: [email protected]

Received date: July 16, 2014; Accepted date: August 30, 2014; Published date: September 03, 2014

Citation: Jan K, Riar CS, Saxena DC (2014) Mathematical Modelling of Thin Layer Drying Kinetics of Biodegradable Pellets. J Food Process Technol 5:370. doi: 10.4172/2157-7110.1000370

Copyright: © 2014 Jan K, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Visit for more related articles at Journal of Food Processing & Technology


The deoiled rice bran, industrial waste (paddy husk, potato peels and banana peels) and plasticizers were added to develop pellets after extrusion for their further use in preparing molds (pots) for horticultural use. Pellets were prepared by extrusion technology using glycerol and cashew nut shell liquid as plasticizers and were dried in a chamber with constant air circulation (1.5 m/sec). Drying kinetics of pellets made from different formulations was studied at different temperatures (60, 70 and 80°C). In declining rate period, moisture transfer from pellets was described by applying 12 different mathematical models, and the corresponding rate constants (k) were calculated. The models in each case were finally selected to explain the drying behaviour based on Chi-square, RMSE and R2. Results indicated that removal of moisture from these pellets was rapid in initial hours and increased with increase in temperature as well. The effect of temperature on rate constant (k) was interpreted according to Arrhenius law. Drying rate and hence, k values were found to be significantly affected by temperature. The energy value increased upon increasing the plasticizer concentration (CNSL and Glycerol) in the formulation. Also, the energy required in the drying of pellets with CNSL as a plasticizer was higher as compared to those with the glycerol.


Deoiled rice bran; Paddy husk; Potato peels; Banana peels; Cashew nut shell liquid


Nursery owners or growers use pots and cell trays of different materials, sizes, shapes, and colours to suit crop species, growing methods, and marketing strategies [1]. At nurseries and greenhouses, seeds, bulbs, and young plants are planted in cell trays or pots containing growing substrate to start the growth under uniform cultivation conditions. However, one of the limitations of impermeable and rigid containers is that roots tend to encircle the outer perimeter of the pot which results in reduced plant growth, health, and survival once transplanted [2]. Most of the containers used for transplanting are made of non-renewable oil-based raw materials such as polystyrene, polyethylene and polypropylene with suitable mechanical properties, chemical and microbial degradation resistance as well as durability. But at the same time disposal and recycling process of post-use plastic pots determine their uncontrolled combustion with the subsequent emission of toxic substances both into the atmosphere and into the soil. Hence, a valid alternative to the employment of petroleum based thermoplastic pots may be represented by the use of biodegradable pots [1-5]. Once buried, the biodegradable pots are subjected to biodegradation process, being transformed in biomass and inorganic products. Now, the attention is being focussed on the development of novel biodegradable and cost-competitive pots or multiple nurseries. These pots are biocomposites whose continuous phase is characterized by biopolymers coming from renewable and available origin, such as polysaccharides [6,7] and whose solid phase, dispersed within the polymeric matrix, is represented by natural fillers and fibres coming from wastes of agro-food and textile processing industries [8]. There is a substantial interest in replacing some or all of the synthetic plastics by biodegradable polymers in several applications. Industrial production of pellets by extrusion process for manufacture of biodegradable products is relatively new and therefore the effectiveness of the availability of pellets rely upon understanding the storage stability of these pellets. A novel kind of low cost, biodegradable and medium water-resistant composite from corn protein was developed by Wu et al. [4]. wherein glycerol-plasticized pots showed good biodegradability and promising water resistance [8].

The thermal energy generated by viscous dissipation during extrusion combines with the shearing effect. Several studies have reported chemical changes during extrusion cooking and related them to product functional qualities e.g. expansion volume, water solubility and color [4,9].

In order to study the consequences of the drying on the pellets it's necessary to possess a deep data of the mass transfer parameters and also the drying kinetics. The mathematical models have been established to be helpful in the new design and/or improvement of drying systems and analysis of mass transfer phenomena concerned throughout drying [10]. The optimal control of a drying process is thus critical and needs complete information regarding the drying behaviour of materials. Drying kinetics may be described using the transport properties such as thermal conductivity, thermal diffusivity, moisture diffusivity as well as interface heat and mass transfer coefficients [11].

Previous research is available on the dilatometry studies of phosphotungstic acid pellets during hydration and dehydration processes [12]. Bassene et al. [13] performed thin layer drying kinetics of pellets prepared from millet flour at different drying air temperatures and found that modified Henderson and Pabis was the best model to describe the behaviour of pellets from millet flour.

Since literature is scanty on the drying behavior/ characteristics of pellets used for development of biodegradable products, the present study was undertaken to describe the drying mechanism of pellets developed from different industrial byproducts, wastes and plasticizers at different drying temperatures.

The objectives of the present study are:

• To study the applicability of several mathematical models selected from literature to the drying of biodegradable pellets of different formulations

• To fit the drying data into the most suitable models by appropriate regression analysis in order to investigate the best model for describing the drying characteristics of the pellets.

Materials and Methods

Procurement of raw material

Deoiled rice bran used for the present study was kindly donated by M/s. AP Solvex Ltd., Dhuri (Punjab, India) which contain high level of protein helpful for binding of material to other ingredients. Glycerol was used as a reference plasticizer in this study was of analytical grade (M/s. Merck Specialities Pvt. Ltd., Mumbai, India). Other plasticizer used in this study was cashew nut shell liquid kindly donated by M/s. Allen Petrochemicals Pvt Lt., Meerut (India). Paddy husk was also provided from M/s. AP Solvex Ltd., Dhuri. Banana peels and potato peels were procured locally from local juice vendors and Hostel Mess as well as Canteen, respectively.

Sample preparation

Preparation of banana and potato peel powder: Paddy husk was ground in a Laboratory grinding mill (M/s. Philips India Limited, Kolkata, India) and passed through 60 mesh screen to obtain fine powder. Banana peels and potato peels were washed with water impregnated with citric acid (1%) and KMnSO4 (0.5%) to check the enzymatic browning. The peels were cut into transverse slices of about 2 mm thickness and dried in tray drier at 65°C overnight. The dried peels were ground in a Laboratory grinding mill as mentioned above and passed through 60 mesh screen to obtain powders. All powders were stored in airtight HDPE plastic packs for further use.

Extrusion process for the preparation of pellets

In the blend preparation, ingredients like deoiled rice bran, paddy husk powder, potato peel powder, banana peel powder, cashew nut shell liquid and glycerol were used at levels shown in Table 1. The moisture was adjusted by sprinkling the distilled water on all the dry ingredients. All the ingredients were weighed and then mixed in a Hobart mixer (M/s. Continental Equipment India Pvt. Ltd, Delhi, India) for 20 min at 328 rpm using sigmoid shape blade. This mixture was then passed through an 80 mesh sieve to eliminate the lumps formed due to addition of moisture. After mixing, samples were stored in polyethylene bags at room temperature for 24h. Extrusion was performed using a co-rotating twin-screw extruder (M/s. Basic Technology Pvt. Ltd., Kolkata, India). The twin screw extruder was kept running for suitable period of time to stabilize the set temperatures and samples were then fed in to feed hopper at the feed rate of 4 kg/h for easy and non-choking operation. Screw speed was set at 125 rpm and the suitable quantity of water was added to makeup the moisture content to 10%. The horizontal split design of the extruder allowed the collection of samples inside the extruder channel with a ‘dead-stop’ run. Samples were taken along the extruder channel and made into pellets using suitable die and high speed cutter. The pellets had on average about 7% moisture content. The preliminary experiments carried out to study the moisture isotherms (unpublished) showed that the pellets being hygroscopic absorbed moisture to about 30% in different storage conditions. This required drying of pellets to suitable moisture content before using injection moulding to prepare pots.

S. No. Model Name Model Equation References
1 Newton MR = exp(−kt) O’Callaghan et al. [18]
2 Page MR = exp(−ktn) Doymaz,[19,20]
3 Thompson t = aln (MR)+ b [ln(MR)]2 Jompob et al. [15]
4 Henderson and Pabis MR = a exp(−kt) Ibrahim et al. [21,22]
5 Modified Henderson and Pabis MR = a exp(−kt) + b exp(−gt) + c exp(−ht) Akpinar [23]
6 Logarithmic MR = a exp(−kt) + c Balasubramanian et al. [24]
7 Two-term exponential MR = a exp(−kt) + (1-a) exp(−kat) Akpinar [23]
8 Two-term MR = a exp(−kt) + b exp(−gt) Ertekin et al. [25]
9 Verma MR = a exp(−kt) + (1-a) exp(−gt) Yaldiz et al. [26]
10 Midilli MR = a exp(−ktn) + bt Midilli et al. [27]
11 Diffusion approach MR = a exp(−kt) + (1-a) exp(−kbt) Cihan et al. [28]
12 Wang and Singh MR = at2 + bt + 1 Mohapatra et al. [29]

Table 1: Mathematical model of drying curves.

Drying kinetics of pellets

The drying experiments carried out in a tray dryer (M/s. Balaji Enterprises, Saharanpur, India) at selected three air drying temperatures (60, 70 and 80°C) to study the drying characteristics of the pellets. The drying cabinet was equipped with an electrical heater, a fan, and temperature indicators. Air velocity was set at 1.5 m/s. When the dryer reached steady state conditions (i.e. when desired temperature was reached) the pellets were spread onto the trays. Moisture loss was recorded at every 15 min intervals during drying. Drying was continued until the moisture content of the sample reached to the equilibrium moisture content. The drying was continued till weight became constant. The experiments were conducted in duplicates.

Drying kinetics model

Various types of mathematical models have been used to describe drying of pellets, ranging from theoretical models based on classical diffusion theory and simplified forms of these to purely empirical models as shown in Table 1.

Statistical evaluation

The results obtained were statistically evaluated using Excel worksheet of MS office version 2010 (MS Office) and Statistica-log software package version 7 (StatSoft Inc., OK, USA). The regression analysis was performed using the STATISTICA computer program. Nonlinear regression, which used to evaluate goodness of fit of the mathematical models to the experimental data are coefficient of determination (R2), the reduced chi-square (χ2) and % Root Mean Square Error (RMSE) was used for data analysis and calculated as below:



Results and Discussion

Drying rate

The variation of moisture content with drying time is obtained for each sample made from different formulations. From these results, drying rates were calculated and then plotted against time. Plots of MR against drying time showed a linear relationship as shown in representative graphs in Figure 1, which indicates that the entire drying curve can be represented approximately by a diffusion model. In all the cases, the drying rate decreased continuously throughout the drying period. The drying took place in the falling rate period. Similar type of observation was reported by Singh et al. [14] wherein the drying rate of sweet potato slices decreased continuously throughout the drying period.


Figure 1: Drying rate curves for pellets made from deoiled rice bran (75 g), potato peel powder (75 g), CNSL (9 ml) at different temperatures.

Selection of drying models

The convective drying of biological materials in the falling rate period is diffusion controlled process and may be represented by Fick’s second law of diffusion. Different types of mathematical models have been used to describe drying of foodstuffs, ranging from theoretical models based on classical diffusion theory and simplified forms of these to purely empirical models [15]. The drying data were converted to moisture ratios which were then fitted to twelve drying models at temperatures 60, 70 and 80°C and the parameters for validation of different equations were shown as in Table 2. The higher value for R2 and the lower values for χ2 and RMSE indicated the better fitness of model Togrul [16], and the best model based on lowest RMSE and χ2 values as well as highest R2 were shown in Table 2 for each formulation.

S.No. Raw material (g) Plasticizer(ml) Code Models Parameters χ2 RMSE R2
1 Deoiled Rice Bran (75) Paddy
Husk (75)
Gly (9) A160 Midilli n: 1.436 ; a: 0.997 ;b: -0.005 ;k: 0.903 0.004 0.001 0.999
2 Gly (15) A260 Logarithmic a: 1.449 ; c:-0.528 ; k: 0.309 0.028 0.036 0.983
3 Gly (21) A360 Verma a:252.609 ; g: 0.616 ; k: 0.617 0.016 0.030 0.990
4 Potato
CNSL (9) B160 Midilli n: 0.735 ; a: 1.002 ; b: -0.108 ; k: 0.615         0.004 0.011 0.998
5 CNSL (15) B260 Midilli n: 1.513 ; a: 0.989 ; b:-0.061 ; k: 0.296 0.003 0.020 0.997
6 CNSL (21) B360 Midilli n: 0.210 ; a: 18.764 ; b: -0.077 ; k: 3.354 0.010 0.014 0.997
7 Banana
Peel (75)
CNSL (9) C160 Diffusion a: 1.587 ; b: 0.999 ; k: 1.435                            0.709 0.060 0.969
8 CNSL (15) C260 Midilli n: 0.220 ; a: 1.000 ; b: -0.185 ; k: 0.519 0.002 0.007 0.999
9 CNSL (21) C360 Midilli n: 0.321 ; a: 0.999 ; b: -0.152 ; k: 0.657 0.011 0.013 0.998
10 Paddy
Husk (75)
CNSL (9) D160 Logarithmic a: 0.778 ; c: .035 ; k: 0.53 0.267 0.082 0.868
11 CNSL (15) D260 Page n: 1.535 ; k: 0.337 0.111 0.044 0.979
12 CNSL (21) D360 Midilli n: 0.639 ; a: 1.005 ; b: -0.030 ; k: 0.117 0.008 0.014 0.997
1 Deoiled Rice Bran (75) Paddy
Husk (75)
Gly (9) A170 Midilli n: 1.114 ; a: 1.000 ; b: -0.090 ; k: 0.989 0.008 0.014 0.997
2 Gly (15) A270 Logarithmic a: 2.013 ; c: -1.041 ; k: 0.406 0.004 0.017 0.995
3 Gly (21) A370 Modified Henderson g:0.660; c: 0.443; h: 0.660; a: 0.100;b:0.444; k:0.7305 0.100 0.048 0.968
4 Potato
Peel (75)
CNSL (9) B170 Diffusion a: 1.000 ; b: -3.879 ; k: 0.898 0.013 0.024 0.993
5 CNSL (15) B270 Logarithmic a: 2.010; c: -1.092 ; k: 0.363 0.038 0.041 0.975
6 CNSL (21) B370 Logarithmic a: 8.806 ; c: -7.820 ; k: 1.989 0.019 0.028 0.990
7 Banana
Peel (75)
CNSL (9) C170 Modified Henderson g: 1.929; c:-186.91; h:2.097; a:89.824;b:98.09; k:1.251 0.032 0.028 0.988
8 CNSL (15) C270 Logarithmic a: 1.327 ; c: -0.377 ; k: 0.665 0.019 0.033 0.985
9 CNSL (21) C370 Logarithmic a: 1.690 ; c: -0.771 ; k: 0.415 0.041 0.053 0.959
10 Paddy
Husk (75)
CNSL (9) D170 Logarithmic a: 10.338 ; c: -9.332 ; k: 2.272 0.006 0.003 0.999
11 CNSL (15) D270 Diffusion a: 0.562 ; b: 25.597 ; k: 0.807 2.880 0.378 0.861
12 CNSL (21) D370 Midilli n: 0.711 ; a : 0.999 ; b: -0.065 ; k: 0.516 0.016 0.020 0.989
1 Deoiled RiceBran (75) Paddy
Husk (75)
Gly (9) A180 Diffusion a: 1.000 ; b: -5.849 ; k: 1.225 0.028 0.023 0.993
2 Gly (15) A280 Midilli n: 0.450 ; a: 0.999 ; b: -0.188 ; k: 0.496 0.019 0.028 0.989
3 Gly (21) A380 Midilli n: -0.079 ; a: -0.079 ; b: -0.079 ; k: 1.210 0.009 0.019 0.996
4 Potato
Peel (75)
CNSL (9) B180 Midilli n: 1.692 ; a: 0.980 ; b: 0.016 ; k: 1.163 0.026 0.0299 0.991
5 CNSL (15) B280 ModifiedHenderson g:1.160; c:0.491; h:1.160; a: 0.091; b:0.491; k:0.689 0.017 0.028 0.990
6 CNSL (21) B380 Page n: 1.281 ; k: 1.312 0.038 0.043 0.983
7 Banana
Peel (75)
CNSL (9) C180 Midilli n: 0.373 ; a: 0.999 ; b: -0.17 ; k: 0.987 0.020 0.026 0.991
8 CNSL (15) C280 ModifiedHenderson g:2.768; c:-152.73; h:2.519; a:81.280; b:72.458;k:0.894 0.010 0.014 0.998
9 CNSL (21) C380 Verma a: 4.076 ; g: 0.220 ; k: 0.345 0.010 0.023 0.995
10 Paddy
Husk (75)
CNSL (9) D180 Midilli g:2.485; c:-382.10; h:2.672; a:183.85; b:199.25;k:2.272 0.063 0.048 0.975
11 CNSL (15) D280 ModifiedHenderson g:2.10; c:0.222; h:35.59; a: 2.665; b:-1.888 ; k:1.665 0.010 0.013 0.998
12 CNSL (21) D380 Midilli n: 0.481 ; a: 0.999 ; b: -0.081 ; k: 0.870 0.009 0.018 0.994

Table 2: Parameters of model equations for drying of pellets with formulations based on different raw material, plasticizers at different temperatures.

The results reveal that the best mathematical model for sample A160 is Midilli when R2 is 0.999 where the model equation is MR = a exp(−ktn) + bt and the parameters are described in Table 2. For different samples, values of R2 are obtained in the range of 0.969 to 1.000. Similar type of results was also reported by Augusta et al. [17] wherein they have determined the highest R2=0.983 and lowest χ2=0.002 for drying of cork planks. The lowest χ2 values for different samples also determine the best fit of the models as in case of B160, the value of χ2 is 0.004. For other samples, the range of the χ2 is in between 0.001 to 0.709. Bassene et al. [13] have reported the similar results of χ2 for drying of millet flour pellets. For all cases, lowest values of RMSE are in between 0.001 to 0.082%.

Effect of temperature variables on drying characteristics

The results suggested that the Arrhenius law might be applicable to relate the dependence of the rate constant on drying air temperature. The relationship is given by the following equation:

k = ko exp(-E/RT)

The computed values for activation energy for pellets are shown in Table 3. As evident that the energy required for drying of pellets containing CNSL is more than that of those pellets which contain glycerol as plasticizer. Moreover, it was also seen that the activation energy increased with the increase in the concentration of the CNSL. The dependence of the rate constant on drying air temperatures and types of plasticizers (glycerol and CNSL) are shown in representative graphs as in Figure 2.


Figure 2: Temperature dependence of the drying rate constant for pellets made from deoiled rice bran (75 g), paddy husk powder (75 g) and glycerol (9, 15, 21 ml) / CNSL (9, 15, 21 ml).

S.No. Raw material (g) Plasticizer (ml) Code E k0 R2
1 Deoiled rice bran (75) Paddy Husk (75) Gly (9) A1 14.828 187.2229 0.939
2 Gly (15) A2 23.176 1347.09 0.995
3 Gly (21) A3 32.732 79379.86 0.912
4 Potato Peel(75) CNSL (9) B1 31.21574 49266.52 0.992
5 CNSL (15) B2 41.052 756153.7 0.907
6 CNSL (21) B3 45.9207 4862192.3 0.998
7 Banana Peel (75) CNSL (9) C1 24.6917 4577.462 0.999
8 CNSL (15) C2 26.55741 7525.105 0.995
9 CNSL (21) C3 31.6472 145947.2 0.951
10 Paddy Husk (75) CNSL (9) D1 71.02318 7.02E+10 0.988
11 CNSL (15) D2 78.1798 4.46E+11 0.998
12 CNSL (21) D3 98.628 4.01E+14 0.939

Table 3: Effect of temperature on drying characteristicsa.


In this study, drying kinetics of pellets were studied at drying air temperatures of 60, 70 and 80°C using a constant air velocity of 1.5 m/s. Mathematical modelling allowed the estimation of drying parameters and different models were found fit for different samples of pellets made from different formulations. It was found that drying rate took place in the falling rate period. The drying rate constant was related to temperature using the Arrhenius relationship. It would be possible to attain faster drying with increased drying temperature up to 80°C with reduced time. It was also seen that the drying energy increased with increasing amount of CNSL in the pellets and the pellets made with glycerol took less energy than that of those made from CNSL.


Select your language of interest to view the total content in your interested language
Post your comment

Share This Article

Relevant Topics

Recommended Conferences

Article Usage

  • Total views: 11757
  • [From(publication date):
    September-2014 - Sep 25, 2017]
  • Breakdown by view type
  • HTML page views : 7966
  • PDF downloads :3791

Post your comment

captcha   Reload  Can't read the image? click here to refresh

Peer Reviewed Journals
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals
International Conferences 2017-18
Meet Inspiring Speakers and Experts at our 3000+ Global Annual Meetings

Contact Us

Agri, Food, Aqua and Veterinary Science Journals

Dr. Krish

[email protected]

1-702-714-7001 Extn: 9040

Clinical and Biochemistry Journals

Datta A

[email protected]

1-702-714-7001Extn: 9037

Business & Management Journals


[email protected]

1-702-714-7001Extn: 9042

Chemical Engineering and Chemistry Journals

Gabriel Shaw

[email protected]

1-702-714-7001 Extn: 9040

Earth & Environmental Sciences

Katie Wilson

[email protected]

1-702-714-7001Extn: 9042

Engineering Journals

James Franklin

[email protected]

1-702-714-7001Extn: 9042

General Science and Health care Journals

Andrea Jason

[email protected]

1-702-714-7001Extn: 9043

Genetics and Molecular Biology Journals

Anna Melissa

[email protected]

1-702-714-7001 Extn: 9006

Immunology & Microbiology Journals

David Gorantl

[email protected]

1-702-714-7001Extn: 9014

Informatics Journals

Stephanie Skinner

[email protected]

1-702-714-7001Extn: 9039

Material Sciences Journals

Rachle Green

[email protected]

1-702-714-7001Extn: 9039

Mathematics and Physics Journals

Jim Willison

[email protected]

1-702-714-7001 Extn: 9042

Medical Journals

Nimmi Anna

[email protected]

1-702-714-7001 Extn: 9038

Neuroscience & Psychology Journals

Nathan T

[email protected]

1-702-714-7001Extn: 9041

Pharmaceutical Sciences Journals

John Behannon

[email protected]

1-702-714-7001Extn: 9007

Social & Political Science Journals

Steve Harry

[email protected]

1-702-714-7001 Extn: 9042

© 2008-2017 OMICS International - Open Access Publisher. Best viewed in Mozilla Firefox | Google Chrome | Above IE 7.0 version