1-s2.0-s0260877405000415-main

Upload: le-quoc-tuan

Post on 05-Apr-2018

217 views

Category:

Documents


0 download

TRANSCRIPT

  • 8/2/2019 1-s2.0-S0260877405000415-main

    1/10

    Determination of suitable thin layer drying curve modelfor some vegetables and fruits

    Ebru Kavak Akpinar *

    Mechanical Engineering Department, Firat University, 23279 Elazig, Turkey

    Received 16 August 2004; accepted 7 January 2005

    Available online 25 February 2005

    Abstract

    This study presents a mathematical modeling of thin layer drying of potato, apple and pumpkin slices in a convective cyclone

    dryer. In order to estimate and select the appropriate drying curve equation, 13 different models, which are semi-theoretical

    and/or empirical, were applied to the experimental data and compared according to their coefficients of determination ( r,v2), which

    were predicted by non-linear regression analysis using the Statistica Computer Program. Moreover, the effects of drying air temper-

    ature, velocity and sample area on the model constants and coefficients were also studied by multiple regression analysis.

    Consequently, of all the drying models, a semi-theoretical MidilliKucuk model was selected as the best one, according to r and v2.

    2005 Elsevier Ltd. All rights reserved.

    Keywords: Drying; Thin layer; Mathematical modelling

    1. Introduction

    Drying of moist materials is a complicated process

    involving simultaneous, coupled heat and mass transfer

    phenomena, which occur inside the material being dried

    (Yilbas, Hussain, & Dincer, 2003).

    Thin layer drying mean to dry as one layer of sample

    particles or slices. Thin layer drying models that describe

    the drying phenomenon of agricultural materials mainly

    fall into three categories, namely theoretical, semi-theo-

    retical and empirical (Midilli, Kucuk, & Yapar, 2002;

    Panchariya, Popovic, & Sharma, 2002). The first takesinto account only internal resistance to moisture transfer

    while the other two consider only external resistance to

    moisture transfer between product and air (Bruce,

    1985; Ozdemir & Devres, 1999; Parti, 1993). The most

    widely investigated theoretical drying model has been

    Ficks second law of diffusion. Drying of many food

    products such as rice (Ece & Cihan, 1993) and hazelnut

    (Demirtas, Ayhan, & Kaygusuz, 1998) has been success-

    fully predicted using Ficks second law. Semi-theoretical

    models offer a compromise between theory and ease of

    use (Fortes & Okos, 1981). Simplifying general series

    solution of second law Ficks or modification of simpli-

    fied models generally derives semi-theoretical models.

    But they are only valid within the temperature, relative

    humidity, and airflow velocity and moisture content

    range for which they were developed. They require small

    time compared to theoretical thin layer models and donot need assumptions of geometry of a typical food, its

    mass diffusivity and conductivity (Parry, 1985). Among

    semi-theoretical thin layer drying models, the Newton

    model, Page model, the modified Page model, the

    Henderson and Pabis model, the logarithmic model,

    the two-term model, the two-term exponential, the diffu-

    sion approach model, the modified Henderson and Pabis

    model, the Verma et al. model and the MidilliKucuk

    model are used widely (Table 2). Empirical models derive

    0260-8774/$ - see front matter 2005 Elsevier Ltd. All rights reserved.

    doi:10.1016/j.jfoodeng.2005.01.007

    * Tel.: +90 424 237 5343; fax: +90 424 241 5526.

    E-mail addresses: [email protected], [email protected]

    www.elsevier.com/locate/jfoodeng

    Journal of Food Engineering 73 (2006) 7584

    mailto:[email protected]:[email protected]:[email protected]:[email protected]
  • 8/2/2019 1-s2.0-S0260877405000415-main

    2/10

    a direct relationship between average moisture content

    and drying time. They neglect fundamentals of the dry-

    ing process and their parameters have no physical mean-

    ing. Therefore they cannot give clear accurate view of the

    important processes occurring during drying althoughthey may describe the drying curve for the conditions

    of the experiments (Ozdemir & Devres, 1999). Among

    them, the Wang and Singh model (see Table 2) has been

    found application in the literature.

    Recently, there have been many studies on the thin

    layer drying and mathematical modeling of various veg-

    etables, fruits and agro-based products such as potato

    (Diamante & Munro, 1993), onion (Rapusas & Driscoll,

    1995; Sarsavadia, Sawhney, Pangavhane, & Singh,

    1999), hazelnut (Ozdemir & Devres, 1999), green pep-

    per, green bean and squash (Yaldiz & Ertekin, 2001),

    tea (Panchariya et al., 2002), green chilli (Hossain &Bala, 2002), banana (Dandamrongrak, Young, &

    Mason, 2002), pistachio (Midilli & Kucuk, 2003) red

    pepper (Akpinar, Bicer, & Yildiz, 2003). Moreover, the

    effects of some parameters related to the product or dry-

    ing conditions such as slice thickness, drying air temper-

    ature, relative humidity, etc. were investigated by many

    researchers (Hossain & Bala, 2002; Midilli & Kucuk,2003; Midilli et al., 2002; Ozdemir & Devres, 1999;

    Sarsavadia et al., 1999; Yaldiz & Ertekin, 2001).

    The main objective of this study was to determine

    and test the most appropriate thin layer drying model

    for understanding the drying behavior of potato, apple

    and pumpkin slices.

    2. Materials and methods

    2.1. Experimental set-up

    Fig. 1 illustrates the schematic diagram of the cyclone

    type dryer, developed for experimental work (Akpinar,

    Nomenclature

    A sample area, the area having the product be-

    fore drying process start, m2

    a, b, c,g, h, n empirical constants in the drying models

    k, k0, k1 empirical coefficients in the drying modelsn number constants

    N number of observations

    MR moisture ratio

    MRexp experimental moisture ratio

    MRpre predicted moisture ratio

    M local moisture content, %dry basis

    Mt mean moisture content at t, %dry basis

    Me mean equilibrium moisture content, %dry

    basis

    Mi initial moisture content, %dry basisr correlation coefficient

    r the diffusion path (m)

    t time, min

    T temperature, C

    V velocity, m/s

    v2 chi-square

    Fig. 1. Experimental set-up: 1Drying chamber; 21st tray; 32nd tray; 4digital balance; 5Observed windows; 6digital thermometer; 7

    the balance bar; 8control panel; 9thermocouples; 10digital thermometer and channel selector; 11rheostat; 12resistance; 13fan; 14wet

    and dry thermometers; 15adjustable flab; 16duct; 17the outlet of air from dryer.

    76 E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584

  • 8/2/2019 1-s2.0-S0260877405000415-main

    3/10

    2002). The system was introduced in the literature

    (Akpinar, Midilli, & Bicer, 2003). Briefly, it consists of

    fan, resistance and heating control systems, air-duct,

    drying chamber in cyclone type, and measurement

    instruments. The heating system consisted of an electric

    4000 W heater placed inside the duct. The rectangular

    duct included air fan and resistance was constructedfrom sheet iron in 1000 mm length, 200 mm width and

    250 mm height. The drying chamber was constructed

    from sheet iron in 600 mm diameter and 800 mm height

    cylinder. The inside and outside surfaces of the drying

    chamber was painted with a spray dye to prevent rust

    in the sheet iron surface. The drying chamber was con-

    structed in concentric form and 30 mm annulus was iso-

    lated by polystyrene. Both topside and bottom side of

    drying chamber was closed. Also, the covers made

    of the steel were isolated by polystyrene. This top

    cover was used to load or unload the chamber. In the

    measurements of temperatures, J type ironconstantan

    thermocouples were used with a manually controlled

    20-channel automatic digital thermometer (ELIMKO,

    6400, Turkey), with reading accuracy of 0.1 C. A

    thermo hygrometer (EXTECH, 444731, China) was

    used to measure humidity levels at various locations of

    the system. A 015 m/s range anemometer (LUTRON,

    AM-4201, Taiwan) measured the velocity of air passing

    through the system. Moisture loss was recorded at

    20 min intervals during drying for determination of dry-

    ing curves by a digital balance (BEL, Mark 3100, Italy)

    in the measurement range of 03100 g and an accuracy

    of 0.01 g.

    2.2. Procedure

    Before drying process, the products were peeled, po-

    tato and apple cut into slices of 12.5 12.5 25 mm

    and 8 8 18 mm (width thickness length) and

    pumpkin cut into slices of 5 mm thickness and 35 mm

    diameters with a mechanical cutter. The trays were

    loaded as thin layer. The first weight of the potato, apple

    and pumpkin was approximately 250, 125 and 200 g,

    respectively. The product slices were carefully and or-

    derly placed on the trays that are made of nylon so that

    the airflow could pass across the trays. The initial andfinal moisture contents of the products were determined

    at 80 C by using an Infrared Moisture Analyzer (MET-

    TLER, LJ16, Switzerland). After the dryer is reached at

    steady state conditions for operation temperatures, the

    samples are put on the trays of dryer and dried there.

    Drying experiments were carried out at 60, 70, and

    80 C drying air temperatures and 1, 1.5 m/s drying air

    velocities. The velocities and temperatures were mea-

    sured in the centre of drying chamber. In the velocity

    measurements, the values of the velocity in the centre

    of the drying chamber were taken into account. The tan-

    gential airflow was across the layer during drying pro-

    cess. Drying was continued until the final moisture

    content of the potato, apple and pumpkin reached to

    approximately 10% (wb), 13% (wb), 6% (wb), respec-

    tively. During the experiments, ambient temperature

    and relative humidity, inlet and outlet temperatures of

    drying air in the duct and dryer chamber were recorded.

    Drying air was tangentially entered in drying chamber.

    In this way, the samples were dried in swirl flow in place

    of uniform flow. The flow diagram of the thin layer dry-

    ing process is presented in Fig. 2.

    2.3. Experimental uncertainty

    Errors and uncertainties in the experiments can arise

    from instrument selection, condition, calibration, envi-

    ronment, observation, and reading, and test planning

    (Akpinar, 2002; Akpinar et al., 2003). In drying

    Fig. 2. The flow diagram of thin layer drying process of products.

    E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584 77

  • 8/2/2019 1-s2.0-S0260877405000415-main

    4/10

    experiments of the apple slices, the temperatures, veloc-

    ity of drying air, weight losses were measured with

    appropriate instruments. During the measurements of

    the parameters, the uncertainties occurred were pre-

    sented in Table 1.

    3. Mathematical modeling of drying curves

    It has been accepted that drying phenomenon of the

    biological products during the falling rate period is con-

    trolled by the mechanism of liquid and/or vapour diffu-

    sion. This behavior suggested strongly an internal mass

    transfer type drying with moisture diffusion as the con-

    trolling phenomena. Assuming that the resistance to

    moisture flow is uniformly distributed throughout the

    interior of the homogeneous isotropic material, the dif-

    fusion coefficient, D is independent of the local moisture

    content and if the volume shrinkage is negligible, Fickssecond law can be derived as follows (Khraisheh,

    Cooper, & Magee, 1997):

    dM

    dt D

    d2M

    dr21

    where, M is the local moisture content (kg water/kg dry

    solids), r is the diffusion path (m), t is the time (s) and D

    is the moisture dependent diffusivity (m2/s).

    Crank (1975) gave the analytical solutions of Eq. (1)

    for various regularly shaped bodies such as rectangular,

    cylindrical and spherical.

    With the appropriate initial and boundary conditions

    t 0; 0 < r< L; M Mi

    t> 0; r 0;dM

    dt 0

    t> 0;

    r

    L;

    M

    Me

    The solution of Eq. (1) for a slice, for a constant diffu-

    sivity, D, in terms of infinite series is given in literature.

    MR Mt MeMi Me

    8

    p2

    X1n1

    1

    2n 12exp

    2n 12p2Dt

    4L2

    " #2

    where, MR is the fractional moisture ratio, Mi is the ini-

    tial moisture, Mt is the mean moisture at time t, Me is

    the mean moisture at equilibrium and L is the half thick-

    ness of the slice for drying from both sides or the thick-

    ness of the slice for drying from one sides.

    For slices shapes, the first boundary condition states

    that moisture is initially uniformly distributed through-

    out the sample. The second implies that the mass trans-

    fer is symmetric with respect to the centre of the slice.

    The third conditions states that the surface moisture

    content of the samples instantaneously reaches equilib-

    rium with the conditions of surroundings air.

    The theoretical and the semi-theoretical models are

    summarized in Table 2. Semi-theoretical thin layer dry-

    ing models are generally derived by simplifying general

    series solution of Ficks second law. For example, the

    Henderson and Pabis model is the first term of a generalseries solution of Ficks second law (Doymaz, 2005). For

    mathematical modeling, the thin layer drying models in

    Table 2 were tested to select the best model for describ-

    ing the drying curve equation of potato, apple and

    pumpkin slices during drying process by the convective

    type cyclone dryer (Akpinar et al., 2003; Diamante &

    Munro, 1993; Hossain & Bala, 2002; Midilli & Kucuk,

    Table 1

    Uncertainties of the parameters during drying of products

    Parameter Unit Comment

    Uncertainty in the temperature

    measurement

    Fan inlet temperature C 0.3800.576

    Heaters outlet temperature C 0.576

    Cyclone inlet temperature C 0.380Cyclone outlet temperature C 0.380

    Centre temperature of product C 0.380

    Temperature between of trays C 0.380

    Ambient air temperature C 0.380

    Inlet of fan with dry and wet

    thermometers

    C 0.5590.707

    Uncertainty in the time measurement

    Mass loss values min 0.1

    Temperature values min 0.1

    Uncertainty in the mass loss measurement g 0.5

    Uncertainty in the air velocity measurement ms1 0.14

    Uncertainty of the measurement of relative

    humidity of air

    RH 0.1

    Uncertainty in the measurement of moisture

    quantity

    g 0.001

    Uncertainty in reading values of table

    (q, cp. . .)

    % 0.10.2

    Table 2

    Thin layer drying curve models

    Model name Model

    Newton MR = exp(kt)Page MR = exp(ktn)Modified Page MR = expb(kt)ncModified Page MR = expb(kt)ncHenderson and Pabis MR = a exp(kt)Logarithmic MR = a exp(kt) + cTwo term MR = a exp(k0t) + b exp(k1t)Two-term exponential MR = a exp(kt) + (1a)exp(kat)Wang and Singh MR = 1 + at + bt2

    Di ffusion approac h MR = a exp(kt) + (1a)exp(kbt)Modified Henderson

    and Pabis

    MR = a exp(kt) + b exp(gt) + c exp(ht)

    Verma et al. MR = a exp(kt) + (1a)exp(gt)MidilliKucuk MR = a exp(ktn) + bt

    78 E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584

  • 8/2/2019 1-s2.0-S0260877405000415-main

    5/10

    2003; Midilli et al., 2002; Mujumdar, 1995; Ozdemir &

    Devres, 1999; Yaldiz & Ertekin, 2001). The regression

    analysis was performed using Statistica Computer

    Program. The correlation coefficient (r) was primary cri-

    terion for selecting the best equation to describe the dry-

    ing curve equation (Guarte, 1996; Midilli et al., 2002). In

    addition to r, the reduced chi-square (v2

    ) was used todetermine the best fit (Midilli & Kucuk, 2003). This

    parameter can be calculated as following:

    v2

    Pn

    i1MRexp;i MRpre;i2

    N n3

    where, MRexp,i is the ith experimentally observed mois-

    ture ratio, MRpre,i the ith predicted moisture ratio, Nthe

    number of observations and n is the number constants

    (Midilli et al., 2002; Sarsavadia et al., 1999).

    Modeling the drying behavior of different agricultural

    products often requires the statistical methods of regres-

    sion and correlation analysis. Linear and non-linear

    regression models are important tools to find the rela-tionship between different variables, especially, for

    which no established empirical relationship exists. In

    this study, the constants and coefficients of the best fit-

    ting model involving the drying variables such as tem-

    perature, velocity of the drying air and product size

    were determined. The effects of these variables on the

    constants and coefficients of drying expression were also

    investigated by multiple linear regression analysis.

    4. Results and discussion

    Potato slices of 83% (wb), apple slices of 87% (wb),

    pumpkin slices of 93% (wb) average initial moisture con-

    tent were dried to 10% (wb), 13% (wb), 6% (wb), respec-

    tively, at temperatures of 60, 70 and 80 C in the

    velocities of drying air of 1 and 1.5 m/s by using a con-

    vective dryer. The final moisture contents represent

    moisture equilibrium between the sample and drying

    air under dryer conditions, beyond which any changes

    in the mass of sample could not occur.

    The moisture content data at the different drying con-

    ditions were converted to the more useful moisture ratio

    expression and then curve fitting computations with the

    drying time were carried on the 13 drying models evalu-

    ated by the previous workers (see Table 2). The statisti-cal analyses results applied to these models at drying

    process at 80 C drying air temperature and 1.5 m/s dry-

    ing air velocity are given in Table 3 for potato slices,

    apple slices and pumpkin slices. The best model describ-

    ing the thin layer-drying characteristic was chosen as the

    one with the highest r-value and the lowest v2 values.

    From Table 3, it was determined that r = 0.99989,

    v2 = 1.79 105 for potato slices, r = 0.99996, v2 =

    1.00 105 for apple slices, r = 0.99967, v2 = 8.38

    105 for pumpkin slices by the MidilliKucuk model.

    The results have shown that the v2 values of the

    MidilliKucuk model lower than the values determined

    Table 3

    Values of the drying constants and coefficients of different models determined through regression method for potato, apple and pumpkin slices

    T= 80 C, V= 1.5 m/s; 1st tray Potato (12.5 12.5 25) Apple (12.5 12.5 25) Pumpkin

    Model name R v2 R v2 R v2

    Newton 0.99871 1.88 104 0.99875 2.35 104 0.98973 2.14 103

    Page 0.99942 8.92 105 0.99930 1.43 104 0.99930 1.55 104

    Modified Page 0.99942 8.92 105 0.99930 1.43 104 0.99930 1.55 104

    Modified Page 0.99871 1.97 104 0.99876 2.53 104 0.98973 2.26 103

    Henderson and Pabis 0.99915 1.30 104 0.99885 2.34 104 0.99235 1.69 103

    Logarithmic 0.99917 1.34 104 0.99993 1.6 105 0.99757 5.72 104

    Two term 0.99970 5 105 0.99885 2.76 104 0.99235 1.91 103

    Two-term exponential 0.99970 5.56 105 0.99869 2.68 104 0.98952 2.31 103

    Wang and Singh 0.95661 6.63 103 0.98977 2.07 103 0.99902 2.17 104

    Diffusion approach 0.99970 4.77 105 0.99941 1.31 104 0.99912 2.08 104

    Modi fied Henderson and Pabis 0.99970 5.56 105 0.99885 3.38 104 0.99234 2.21 103

    Verma et al. 0.99969 4.78 105 0.99940 1.31 104 0.99911 2.09 104

    Midilli and Kucuk 0.99989 1.79 105 0.99996 1.00 105 0.99967 8.38 105

    T=60C Potato

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800

    Drying time (min)

    V = 1.5 m/s, 8x8x18 mm, 1st trayV = 1.5 m/s, 8x8x18 mm, 2nd trayV = 1.5 m/s, 12.5x12.5x25 mm, 1st trayV = 1.5 m/s, 12.5x12.5x25 mm, 2nd trayV = 1 m/s, 8x8x18 mm, 1st trayV = 1 m/s, 8x8x18 mm, 2nd trayV = 1 m/s, 12.5x12.5x25 mm,1st trayV = 1m/s, 12.5x12.5x25 mm, 2nd trayMidilli-Kucuk modelDiffusion approach model

    MR=(M

    t-M

    e)/(M

    i-M

    e)

    Fig. 3. Variation of the experimental and predicted moisture ratio by

    the MidilliKucuk model and diffusion approach model with drying

    time at 60 C of drying air for potato slices.

    E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584 79

  • 8/2/2019 1-s2.0-S0260877405000415-main

    6/10

    by other models. It was noticed that the MidilliKucuk

    model gave the highest r and the lowest v2 for all drying

    conditions. Figs. 311 present the variations of moisture

    ratio versus drying time for potato, apple and pumpkin

    slices dried at the different drying air temperatures,

    velocities and sample size. Additionally, Figs. 311 show

    the comparison of experimental and predicted moistureratio by the MidilliKucuk model and the model has

    correlation coefficient and chi-square, which is near to

    this model. The results of non-linear regression analyses

    and of statistical analyses applied to the MidilliKucuk

    model for all drying conditions have shown in Table 4

    for potato slices, Table 5 for apple slices, Table 6 for

    pumpkin slices. Generally, r-values obtained by using

    this model were varied between 0.999770.99995 for

    potato slices, 0.999650.99997 for apple slices and

    0.999400.99985 for pumpkin slices (see Tables 46).

    As shown in Figs. 311, the MidilliKucuk model

    showed good agreement with the experimental data

    and gave the best results for potato, apple and pumpkin

    slices according to r and v2. Therefore, the Midilli

    Kucuk was selected to represent the thin layer-drying

    behavior of these agricultural products according to

    the highest r and the lowest v2. Consequently, it can

    be said that the MidilliKucuk model could sufficiently

    define the thin layer drying of potato, apple and pump-

    kin slices.

    T=70C Potato

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500 550 600 650 700

    Drying time (min)

    V=1.5 m/s, 8x8x18 mm, 1st trayV=1.5 m/s, 8x8x18 mm, 2nd trayV=1.5 m/s, 12.5x12. 5x25 mm, 1st trayV=1.5 m/s, 12.5x12. 5x25 mm, 2nd trayV=1 m/s, 8x8x18 mm, 1st trayV=1 m/s, 8x8x18 mm, 2nd trayV=1 m/s, 12.5x12.5x25 mm, 1st trayV=1 m/s, 12.5x12.5x25 mm, 2nd trayMidilli-Kucuk modelDiffusion approach model

    MR=(Mt-

    Me

    )/(Mi-

    Me)

    Fig. 4. Variation of the experimental and predicted moisture ratios bythe MidilliKucuk model and diffusion approach model with drying

    time at 70 C of drying air for potato slices.

    T=80 C Potato

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500 550 600 650

    Drying time (min)

    V=1.5 m/s, 8x8x18 mm, 1st trayV=1.5 m/s, 8x8x18 mm, 2nd trayV=1.5 m/s, 12.5x12.5 x25 mm, 1st trayV=1.5 m/s, 12.5x12.5 x25 mm, 2nd trayV=1 m/s, 8x8x18 mm, 1st trayV=1 m/s, 8x8x18 mm, 2nd trayV=1 m/s, 12.5x12.5x25 mm, 1st trayV=1 m/s, 12.5x12.5x25 mm, 2nd trayMidilli-Kucuk modelDiffusion approach model

    MR=(M

    t-M

    e)/(M

    i-M

    e)

    Fig. 5. Variation of the experimental and predicted moisture ratios by

    the MidilliKucuk model and diffusion approach model with drying

    time at 80 C of drying air for potato slices.

    T=60C Apple

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500 550

    Drying time (min)

    V=1.5 m/s, 8x8x18 mm, 1st trayV=1.5 m/s, 8x8x18 mm, 2nd trayV=1.5 m/s, 12.5x12.5x 25 mm, 1st trayV=1.5 m/s, 12.5x12.5 x25 mm, 2nd trayV=1 m/s, 8x8x18 mm, 1st t rayV=1 m/s, 8x8x1 8 mm, 2n d trayV=1 m/s, 12.5x12.5x25 mm, 1st trayV=1 m/s, 12.5x12.5x25 mm, 2nd trayMidilli-Kucuk modelLogarithmic model

    MR=(Mt-

    Me)

    /(Mi-

    Me

    )

    Fig. 6. Variation of the experimental and predicted moisture ratios by

    the MidilliKucuk model and Logarithmic model with drying time at

    60 C of drying air for apple slices.

    T=70 C Apple

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500

    Drying time (min)

    V=1.5 m/s, 8x8x18 mm, 1st trayV=1.5 m/s, 8x8x18 mm, 2nd trayV=1.5 m/s, 12.5x12.5x2 5 mm, 1st trayV=1.5 m/s, 12.5x12.5x2 5 mm, 2nd trayV=1 m/s, 8x8x18 mm, 1st trayV=1 m/s, 8x8x18 mm, 2nd trayV=1 m/s, 12.5x12.5x25 mm, 1st trayV=1 m/s, 12.5x12.5x25 mm, 2nd trayMidilli-Kucuk modelLogarithmic

    MR=(M

    t-M

    e)/(M

    i-M

    e)

    Fig. 7. Variation of the experimental and predicted moisture ratios by

    the MidilliKucuk model and Logarithmic model with drying time at

    70 C of drying air for apple slices.

    80 E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584

  • 8/2/2019 1-s2.0-S0260877405000415-main

    7/10

    The fitting procedure indicated that the mentioned

    results of the MidilliKucuk model could be used to

    model the thin layer drying behavior of these agricul-

    tural products, but it did not indicate the effect of dryingconditions. To take into account the effect of the drying

    variables on the MidilliKucuk model constants a, k, n

    and b were regressed against those of drying air temper-

    ature, velocity and sample area using multiple regression

    analysis. All possible combinations of the different dry-

    ing variables were tested and included in the regression

    analysis. Based on the multiple regression analysis, the

    accepted model, the constants and coefficients were as

    follows:

    MRa; k; b; t Mt MeMi Me

    a expktn bt 4

    where

    for potato,

    a 0:986173 0:000069 T 0:005702 V 0:098206 A 5k 0:015582 0:000156 T 0:013467 V 0:266761 A 6

    n 1:218379 0:000802 T 0:162776 V 138:525 A 7

    b 0:0000085 0:00000029 T 0:0000393 V 0:0203022 A

    8

    for apple,

    a 1:004084 0:000073 T 0:001960 V 3:944759 A 9

    k 0:006391 0:000065 T 0:009775 V 1:576723 A 10

    n 1:187734 0:002467 T 0:128878 V 202:536 A 11

    b 0:000082 0:000002 T 0:000041 V 0:041667 A 12

    T=80 C Apple

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400

    Drying time (min)

    V=1.5 m/s, 8x8x18 mm, 1st trayV=1.5 m/s, 8x8x18 mm, 2nd trayV=1.5 m/s, 12.5x12.5x 25 mm, 1st trayV=1.5 m/s, 12.5x12.5 x25 mm, 2nd trayV=1 m/s, 8x8x18 mm, 1st t rayV=1 m/s, 8x8x1 8 mm, 2 nd trayV=1 m/s, 12.5x12.5x25 mm, 1st trayV=1 m/s, 12.5x12.5x25 mm, 2nd trayMidilli-Kucuk modelLogarithmic model

    MR=(M

    t-M

    e)/(M

    i-M

    e)

    Fig. 8. Variation of the experimental and predicted moisture ratios by

    the MidilliKucuk model and Logarithmic model with drying time at

    80 C of drying air for apple slices.

    T=60 C Pumpkin

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800

    Drying time (min)

    V=1.5 m/s, 1st tray

    V=1.5 m/s, 2nd tray

    V=1 m/s, 1st tray

    V=1 m/s, 2nd tray

    Midilli-Kucuk model

    Page model

    MR=(M

    t-M

    e)/(M

    i-M

    e)

    Fig. 9. Variation of the experimental and predicted moisture ratio by

    the MidilliKucuk model and Page model with drying time at 60 C of

    drying air for pumpkin slices.

    V=1.5 m/s, 1st tray

    V=1.5 m/s, 2nd tray

    V=1 m/s, 1st tray

    V=1 m/s, 2nd tray

    Midilli-Kucuk model

    Page model

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500 550 600

    Drying time (min)

    T=70C Pumpkin

    MR=(Mt-

    Me)/(Mi-

    Me

    )

    Fig. 10. Variation of the experimental and predicted moisture ratio by

    the MidilliKucuk model and Page model with drying time at 70 C of

    drying air for pumpkin slices.

    T=80C Pumpkin

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    0.8

    0.9

    1

    0 50 100 150 200 250 300 350 400 450 500

    Drying time (min)

    V=1.5 m/s, 1st tray

    V=1.5 m/s, 2nd tray

    V=1 m/s, 1st tray

    V=1 m/s, 2nd tray

    Midilli-Kucuk model

    Page model

    MR=(M

    t-M

    e)/(M

    i-M

    e)

    Fig. 11. Variation of the experimental and predicted moisture ratios

    by the MidilliKucuk model and Page model with drying time at 80 C

    of drying air for pumpkin slices.

    E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584 81

  • 8/2/2019 1-s2.0-S0260877405000415-main

    8/10

    Table 4

    Values of the drying constant and coefficients of the MidilliKucuk model determined through regression method for potato slices at all drying

    conditions

    Drying air

    temperature T, C

    Air flow

    rate V, m/s

    Sample

    area A, m2Tray no. a k n b r v2

    80 1.5 0.000544 1 0.9989 0.0169 0.9804 0. 000052 0.99994 1. 17 105

    70 1.5 0.000544 1 0.9995 0.0166 0.9457 0. 000048 0.99989 2. 05 105

    60 1.5 0.000544 1 0.9989 0.0149 0.9351 0. 000038 0.99994 1. 01 105

    80 1.5 0.001250 1 0.9973 0.0170 0.8893 0. 000051 0.99989 1. 79 105

    70 1.5 0.001250 1 1.0069 0.0176 0.8484 0. 000049 0.99980 3. 14 105

    60 1.5 0.001250 1 0.9957 0.0134 0.8459 0. 000075 0.99980 2. 93 105

    80 1 0.000544 1 1.0001 0.0115 1.0214 0. 000041 0.99993 1. 25 105

    70 1 0.000544 1 0.9991 0.0080 1.0594 0. 000018 0.99993 1. 03 105

    60 1 0.000544 1 0.9968 0.0085 1.0109 0. 000022 0.99995 7. 65 106

    80 1 0.001250 1 0.9966 0.0120 0.9153 0. 000041 0.99993 9. 73 106

    70 1 0.001250 1 0.9971 0.0117 0.8992 0. 000038 0.99992 1. 19 105

    60 1 0.001250 1 0.9963 0.0084 0.9283 0. 000034 0.99995 6. 33 106

    80 1.5 0.000544 2 0.9996 0.0171 0.9686 0. 000016 0.99994 9. 05 106

    70 1.5 0.000544 2 0.9985 0.0149 0.9559 0. 000049 0.99984 2. 98 105

    60 1.5 0.000544 2 0.9967 0.0161 0.9095 0. 000030 0.99992 1. 15 105

    80 1.5 0.001250 2 0.9977 0.0174 0.8756 0. 000045 0.99985 1. 86 105

    70 1.5 0.001250 2 1.0047 0.0158 0.8564 0. 000051 0.99987 1. 98 105

    60 1.5 0.001250 2 1.0016 0.0119 0.8623 0. 000074 0.99977 3. 53 105

    80 1 0.000544 2 0.9976 0.0100 1.0341 0. 000017 0.99992 1. 28 105

    70 1 0.000544 2 0.9991 0.0080 1.0594 0. 000018 0.99995 8. 63 106

    60 1 0.000544 2 0.9936 0.0058 1.0678 0. 000025 0.99990 1. 72 105

    80 1 0.001250 2 0.9963 0.0093 0.9568 0. 000022 0.99994 8. 68 106

    70 1 0.001250 2 0.9957 0.0087 0.9438 0. 000032 0.99990 1. 53 105

    60 1 0.001250 2 0.9934 0.0071 0.9532 0. 000034 0.99991 1. 29 105

    Table 5

    Values of the drying constant and coefficients of the MidilliKucuk model determined through regression method for apple slices at all drying

    conditions

    Drying air

    temperature T, C

    Air flow

    rate V, m/s

    Sample

    area A, m2Tray no. a k n b r v2

    80 1.5 0.000544 1 0.9987 0.0167 1.0520 0.000185 0.99990 3.50 105

    70 1.5 0.000544 1 0.9970 0.0140 1.0325 0.000156 0.99974 7.63 105

    60 1.5 0.000544 1 0.9988 0.0164 0.9681 0.000130 0.99995 1.11 105

    80 1.5 0.001250 1 1.0014 0.0126 0.9893 0.000137 0.99996 1.00 105

    70 1.5 0.001250 1 1.0028 0.0147 0.9348 0.000051 0.99987 2.42 105

    60 1.5 0.001250 1 1.0003 0.0157 0.8898 0.000046 0.99990 1.67 105

    80 1 0.000544 1 0.9991 0.0112 1.1375 0.000078 0.99997 9.35 105

    70 1 0.000544 1 0.9990 0.0076 1.1481 0.000119 0.99991 2.63 105

    60 1 0.000544 1 0.9996 0.0058 1.1445 0.000057 0.99997 7.29 106

    80 1 0.001250 1 1.0011 0.0108 0.9859 0.000098 0.99994 1.20 105

    70 1 0.001250 1 0.9988 0.0125 0.9261 0.000108 0.99993 1.30 105

    60 1 0.001250 1 1.0076 0.0115 0.9145 0.000056 0.99982 3.00 105

    80 1.5 0.000544 2 0.9971 0.0164 1.0546 0.000130 0.99965 1.16 104

    70 1.5 0.000544 2 0.9993 0.0121 1.0514 0.000123 0.99994 1.71 105

    60 1.5 0.000544 2 0.9982 0.0118 1.0220 0.000072 0.99991 2.14 105

    80 1.5 0.001250 2 0.9998 0.0116 1.0021 0.000110 0.99992 1.86 105

    70 1.5 0.001250 2 1.0018 0.0144 0.9277 0.000076 0.99982 3.40 105

    60 1.5 0.001250 2 0.9990 0.0140 0.9040 0.000055 0.99983 2.80 105

    80 1 0.000544 2 0.9999 0.0108 1.1289 0.000071 0.99998 4.30 106

    70 1 0.000544 2 0.9983 0.0062 1.1746 0.000143 0.99984 4.68 105

    60 1 0.000544 2 0.9984 0.0054 1.1558 0.000068 0.99997 8.18 106

    80 1 0.001250 2 0.9998 0.0100 0.9896 0.000103 0.99994 1.24 105

    70 1 0.001250 2 0.9973 0.0107 0.9496 0.000073 0.99989 1.99 105

    60 1 0.001250 2 1.0071 0.0094 0.9435 0.000059 0.99987 2.33 105

    82 E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584

  • 8/2/2019 1-s2.0-S0260877405000415-main

    9/10

    for pumpkin,

    a 0:966467 0:000184 T 0:007014 V 13k 0:005645 0:000095 T 0:003791 V 14

    n 0:572175 0:009074 T 0:064652 V 15

    b 0:000050 0:000001 T 0:000024 V 16

    These expressions can be used to estimate the mois-

    ture ratio of potato, apple and pumpkin slices at any

    time during the drying process with a great accuracy.

    The consistency of the model and relationship between

    the coefficients and drying variables evident with

    rpotato 0:9984; v2

    potato

    2:26 104 and

    rapple 0:9976; v2apple 4:03 10

    4 and

    rpumpkin 0:9955; v2pumpkin 7:76 10

    4

    5. Conclusion

    In order to explain the drying behavior and develop

    the mathematical modeling of agricultural products as

    potato, apple and pumpkin, 13 models in the literature

    were applied. Among these models, in each of three

    products, the MidilliKucuk model gave the best results

    and showed good agreement with the experimental data

    obtained from the experiments including the thin layer

    drying process. When the effects of drying air tempera-

    ture, velocity and sample area on the constants and

    coefficients of the MidilliKucuk model were examined,

    the resulting model gave an r of 0.9984 and v2 of

    2.26 104 for potato slices, and an r of 0.9976 and v2

    of 4.03 104 for apple slices, and an r of 0.9955 and

    v2 of 7.76 104 for pumpkin slices. According to re-

    sults, it can be said that the MidilliKucuk model ade-

    quately described the drying behavior of potato, apple

    and pumpkin slices in the drying process at a tempera-

    ture range 6080 C and a velocity range 11.5 m/s of

    drying air.

    Acknowledgement

    The author thanks Prof. Ibrahim Dincer from the

    University of Ontario Institute of Technology and Dr.

    Adnan Midilli from Nigde University, and Firat Univer-

    sity Research Foundation (FUNAF) financial support,

    under project number 357.

    References

    Akpinar, E. K. (2002). The development of a cyclone type dryer for

    agricultural products. PhD Thesis. Firat University, Elazig,

    Turkey.

    Akpinar, E. K., Bicer, Y., & Yildiz, C. (2003). Thin layer drying of red

    pepper. Journal of Food Engineering, 59, 99104.

    Akpinar, E., Midilli, A., & Bicer, Y. (2003). Single layer drying

    behavior of potato slices in a convective cyclone dryer and

    mathematical modeling. Energy Conversion and Management, 44,

    16891705.

    Bruce, D. M. (1985). Exposed-layer barley drying, three models fitted

    to new data up to 150 C. Journal of Agricultural Engineering

    Research, 32, 337347.

    Crank, J. (1975). The mathematics of diffusion. Oxford: Clarendon

    Press.

    Dandamrongrak, R., Young, G., & Mason, R. (2002). Evaluation ofvarious pre-treatments for the dehydration of banana and selection

    of suitable drying models. Journal of Food Engineering, 55,

    139146.

    Demirtas, C., Ayhan, T., & Kaygusuz, K. (1998). Drying behaviour of

    hazelnuts. Journal of the Science of Food and Agriculture, 76,

    559564.

    Diamante, L. M., & Munro, P. A. (1993). Mathematical modeling of

    the thin layer solar drying of sweet potato slices. Solar Energy, 51,

    271276.

    Doymaz, I. (2005). Drying behaviour of green beans. Journal of Food

    Engineering, 69(2), 161165.

    Ece, M. C., & Cihan, A. (1993). A liquid diffusion model for drying

    rough rice. Transactions of American Society of Agricultural

    Engineers, 36, 837840.

    Table 6

    Values of the drying constant and coefficients of the MidilliKucuk model determined through regression method for pumpkin slices at all drying

    conditions

    Drying air

    temperature T, C

    Air flow

    rate V, m/s

    Tray no. a k n b r v2

    80 1.5 1 0.9798 0.0021 1.3114 0.000071 0.99967 8.38 105

    70 1.5 1 0.9932 0.0056 1.0561 0.000100 0.99961 8.13 105

    60 1.5 1 0.9936 0.0067 0.9786 0.000054 0.99965 6.04 105

    80 1 1 0.9979 0.0037 1.1320 0.000061 0.99985 3.18 105

    70 1 1 0.9883 0.0031 1.1269 0.000084 0.99940 1.25 104

    60 1 1 0.9698 0.0021 1.1793 0.000017 0.99956 7.92 105

    80 1.5 2 0.9878 0.0021 1.2951 0.000068 0.99967 8.38 105

    70 1.5 2 0.9928 0.0056 1.0404 0.000098 0.99972 5.54 105

    60 1.5 2 0.9922 0.0062 0.9804 0.000058 0.99975 4.25 105

    80 1 2 0.9901 0.0024 1.2057 0.000081 0.99976 5.42 105

    70 1 2 0.9867 0.0026 1.1321 0.000084 0.99967 6.74 105

    60 1 2 0.9854 0.0030 1.0800 0.000050 0.99967 6.05 105

    E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584 83

  • 8/2/2019 1-s2.0-S0260877405000415-main

    10/10

    Fortes, M., & Okos, M. R. (1981). Non-equilibrium thermodynam-

    ics approach to heat and mass transfer in corn kernels.Transactions

    of American Society of Agricultural Engineers, 22, 761769.

    Guarte, R. C. (1996). Modeling the drying behavior of copra and

    development of a natural convection dryer for production of high

    quality copra in the Philippines. PhD Dissertation. 287, Hohen-

    heim University, Stuttgart, Germany.

    Hossain, M. A., & Bala, B. K. (2002). Thin-layer drying characteristics

    for green chilli. Drying Technology, 20(2), 489505.

    Khraisheh, M. A. M., Cooper, T. J. R., & Magee, T. R. A. (1997). The

    transport mechanism of moisture during air drying processes.

    Trans IChemE, 75(C), 3439.

    Midilli, A., & Kucuk, H. (2003). Mathematical modelling of thin layer

    drying of pistachio by using solar energy. Energy Conversion and

    Management, 44, 11111122.

    Midilli, A., Kucuk, H., & Yapar, Z. A. (2002). New model for single-

    layer drying. Drying Technology, 20(7), 15031513.

    Mujumdar, A. S. (1995). Handbook of industrial drying (Vol. I). New

    York: Marcel Dekker, pp. 71111.

    Ozdemir, M., & Devres, Y. O. (1999). The thin layer drying

    characteristics of hazelnuts during roasting. Journal of Food

    Engineering, 42, 225233.

    Panchariya, P. C., Popovic, D., & Sharma, A. L. (2002). Thin-layer

    modelling of black tea drying process. Journal of Food Engineering,

    52, 349357.

    Parry, J. L. (1985). Mathematical modelling and computer simulation

    of heat and mass transfer in agricultural grain drying. Journal of

    Agricultural Engineering Research, 54, 339352.

    Parti, M. (1993). Selection of mathematical models for drying grain in

    thin layers. Journal of Agricultural Engineering Research, 54, 339

    352.

    Rapusas, R. S., & Driscoll, R. H. (1995). The thin-layer drying

    characteristics of white onion slices. Drying Technology, 13(89),

    19051931.

    Sarsavadia, P. N., Sawhney, R. L., Pangavhane, D. R., & Singh, S. P.

    (1999). Drying behavior of brined onion slices. Journal of Food

    Engineering, 40, 219226.

    Yaldiz, O., & Ertekin, C. (2001). Thin layer solar drying of some

    vegetables. Drying Technology, 19, 583596.

    Yilbas, B. S., Hussain, M. M., & Dincer, I. (2003). Heat and moisture

    diffusion in slab products to convective boundary condition. Heat

    and Mass Transfer, 39, 471476.

    84 E.K. Akpinar / Journal of Food Engineering 73 (2006) 7584