versión On-line ISSN 0718-5839
Chilean J. Agric. Res. vol.72 no.2 Chillán jun. 2012
Chilean Journal of Agricultural Research 72(2) April - June 2012
Tomato Infrared Drying: Modeling and some Coefficients of the Dehydration Process
Secado de Tomate por Infrarrojo: Modelación y algunos Coeficientes del Proceso de Deshidratación.
Paulo Cesar Corrêa1, Gabriel Henrique Horta de Oliveira2*, Fernanda Machado Baptestini1, Mayra Darliane Martins Silva Diniz1, and Aline Almeida da Paixão1
1Universidade Federal de Viçosa, Agricultural Engineering Department, Campus UFV, Av. PH Rolfs, ssn, 36570-000, Viçosa, Minas Gerais, Brasil.
2Instituto Federal de Brasilia, Campus Gama, Praça 2, s/n, 72405980, Brasília-DF, Brasil."Corresponding author (email@example.com).
Infrared dehydration is more advantageous than the convective system under similar conditions, and studying this process is important to further develop equipment and procedures. Thus, the aim of this work was to study the dehydration process of tomato (Solanum lycopersicum L.) slices obtained by infrared drying at three different maturity stages throughout two different procedures: Firstly, the drying model was determined by applying the mass and energy balances under wet bulb temperature for the constant drying rate period and secondly, the mass effective diffusion coefficient was determined throughout the experimental data and the theory of diffusion of the liquid phase for the decreasing drying rate period. Tomato fruits cv. Santa Cruz were used. Three maturity stages were selected: green (stage 1), orange (stage 2), and red (stage 3). Mathematical models frequently used to represent drying of agricultural products were fitted to the experimental data of tomato drying. The effective diffusion coefficient was obtained by adjusting the liquid diffusion mathematical model to the experimental data of the descending period of dehydration. The two-term model was the best one to represent the tomato dehydration process. The critical moisture content for tomato dehydration was 2.97 kgw kgdm-1. There is an initial dehydration period in which the drying rate reaches its maximum (approximately 1.05 kgw kgdm-1, about 3 min). Three different methods were used to obtain values of the effective diffusion coefficient, including the finite element method, which had the lowest values for the least square sum of deviation 1.00 x 10-7 m2 s-1. The global coefficient of heat transfer was 12.45 W m-2K-1, and the global coefficient of mass transfer was 0.0105 m s-1.
Key words: Critical moisture content, diffusion, finite element method, heat and mass transfer, mathematical modeling, Solanum lycopersicum.
La deshidratación por infrarrojo es más ventajosa que el sistema convectivo, en condiciones similares, y el estudio de este proceso es importante para el desarrollo de equipos y procedimientos. Por tanto, el presente trabajo tuvo como objetivo el estudio del proceso de deshidratación en tres diferentes estados de maduración de rebanadas de tomate (Solanum lycopersicum L.) por medio del secado por infrarrojo, usando dos métodos: aproximación del período constante de deshidratación a la teoría de la transferencia de calor y masa para el termómetro de bulbo húmedo; y al período decreciente de la deshidratación a la teoría de difusión de líquidos. Se utilizaron frutos de tomate cv. Santa Cruz. Tres estados de maduración fueron obtenidos: verde (estado 1), naranja (estado 2) y rojo (estado 3). Modelos matemáticos usualmente utilizados para representar el secado de productos agrícolas fueron ajustados a los datos experimentales del secado de tomate. El coeficiente de difusión efectivo fue obtenido por medio del ajuste del modelo matemático de difusión líquida a los datos experimentales del período decreciente de deshidratación. El modelo Dos Términos representó mejor el proceso del secado de rebanadas de tomate. La humedad crítica para la deshidratación de tomate fue 2.97 kgw kgdm-1. Existe un período inicial de deshidratación en que la tasa de secado alcanzó el máximo (aproximadamente 1.05 kgw kgdm-1, alrededor de 3 min). Tres diferentes métodos fueron utilizados para obtener los valores del coeficiente efectivo de difusión, incluyendo el método de elementos finitos, en que presentó los menores valores de la suma de desviaciones mínimas al cuadrado, 1.00 x 10-7 m2 s-1. El coeficiente global de transferencia de calor fue 12.45 W m-2 K-1, y el coeficiente global de transferencia de masa fue 0.0105 m s-1.
Palabras clave: humedad crítica, difusión, método de elementos finitos, transferencia de calor y masa, modelación matemática, Solanum lycopersicum.
Tomato (Solanum lycopersicum L.) is one of the vegetables most consumed in the world. It is a source of vitamins A and C and mineral salts, such as K and Mg. Tomato is a climacteric fruit with a short maturation and color change period, firmness, aroma, and flavor throughout shelf-life (Oliveira, 2010).
Dehydration has recently been used to increase shelf-life of several fruits and vegetables (Vasques et al., 2006). Conservation through dehydration is based on the fact that microorganisms, enzymes, and the entire metabolic system of biological materials require water for their activities. A reduced amount of water decreases water activity and chemical reactions within the product, as well as the development of microorganisms (Christensen and Kaufmann, 1974). Therefore, the main objective of drying is reducing moisture content and water activity (Doymaz, 2005a; Corrêa et al., 2010a; Oliveira et al., 2010) since these parameters control product respiration rate and incidence of microorganisms and insects (Corrêa et al., 2010b; Goneli et al, 2010a; 2010b).
Infrared dehydration is more advantageous than the convective system under similar conditions (Sharma et al., 2005). When infrared is used to warm up or dehydrate wet materials, radiation encounters the exposed product, penetrates it, and radiation energy is converted into heat. The penetration depth of radiation depends on the material's properties and the radiation wavelength. When the product is exposed to radiation, it is intensively heated and the temperature gradient within the product is reduced in a short period of time.
Depending on the material to be dried, water can move through its interior by different mechanisms. In hygroscopic capillary-porous products, such as most agricultural products, the possible mechanisms of moisture movement are liquid diffusion, capillary diffusion, surface diffusion, hydrodynamic flow, vapor diffusion, and thermal diffusion (Brooker et al., 1992). The effective diffusion coefficient has been studied by different researchers in different products (Campos et al., 2009; Corrêa et al., 2010c); this parameter illustrates the moisture movement, intrinsic to the product, towards the outside due to drying.
Heat and mass transfer is intrinsically involved in reducing grain moisture (Hall, 1980) by drying. Simulating the behavior of each product during dehydration is an important parameter to develop and maximize grain drying equipment (Sacilik et al., 2006; Corrêa et al., 2010c), as well as selecting mathematical models that will represent moisture loss during drying. Moreover, transport phenomena of food and other important biological materials is a significant link between the processing of these materials and product quality and safety (Welti-Chanes et al., 2005).
Given the importance of studies on drying, the objective of this study was to model the dehydration process of tomato fruits at three different maturity stages, and also obtain the effective diffusion coefficient and global values of heat and mass transfer.
MATERIALS AND METHODS
The present work was conducted in the Laboratory of Physical Properties and Quality Evaluation of Agricultural Products at the National Grain Storage Training Center - CENTREINAR, Federal University of Viçosa, Viçosa, Minas Gerais, Brazil.
Tomato fruits cv. Santa Cruz were acquired at the local market. Three maturity stages were selected: green (stage 1), orange (stage 2), and red (stage 3). Each maturity stage had five replicates. Tomatoes were sliced approximately 5 mm thick and weighed on a digital scale with 0.001 g precision.
The study used a dryer with an infrared radiation source (model IV 2500, Gehaka, Sao Paulo, Brazil) for dehydration. This equipment contains a scale with 0.001 g precision and an automatic data acquisition system. The radiation source is located at a fixed distance of 15 mm from the product with an infrared power of 300 W. Dehydration was carried out at 60 °C. The moisture ratio of the tomato during drying was obtained by Equation :
where MR is the moisture ratio of the product, dimensionless; Ut* is the moisture content of the product at a certain drying time, decimal d.b.; Ui* is the initial moisture content of the product, decimal d.b.; and Ue* is the equilibrium moisture content of the product, decimal d.b.
The mathematical models frequently applied to represent drying of agricultural products were fitted to the experimental data of tomato drying. These models are shown in Table 1.
The goodness of fit of the mathematical models was analyzed through non-linear regression and the Gauss Newton method with Statistica 8.0 (StatSoft, Tulsa, Oklahoma, USA). The fitness degree of each model was based on the determination coefficient (R2), mean relative percent deviation (P), and standard error (SE). The values of P and SE for each model were calculated by Equations  and :
where P is the mean relative percent deviation (%); SE is the standard error, decimal d.b.; Y is the observed value, dimensionless; Y is the value estimated by the model, dimensionless; n is the number of observed data; and DF is the degrees of freedom of the model (number of observed data minus the number of model parameters).
The theories that govern the constant rate period of dehydration of agricultural products, according to Brooker et al. (1992), can be approximated to the heat and mass transfer balance studies for wet bulb temperature. Thus, the constant rate of dehydration can be represented by Equation :
where hc is the global coefficient of heat transfer (W m-2 K-1); hm is the global coefficient of mass transfer (m s-1); L is the latent heat of vaporization (J kg-1); Rv is the universal constant of water vapor (461.91 J kg-1 K-1); A is the superficial area of tomato slices (m2); Pvbu is the vapor pressure for the wet bulb temperature (Pa); Pv∞ is the vapor pressure for dehydration temperature (Pa); Tbu is the wet bulb temperature (K); Tc is the dehydration temperature (K); and dM/dt is the constant rate of dehydration (kg s-1).
Coefficients were calculated with the following parameters: plate geometric shape type, isothermal process, and neglecting volumetric shrinkage of the product.
The effective diffusion coefficient was obtained by adjusting the liquid diffusion mathematical model (Equation ) to the experimental data of the descending period of dehydration. This equation is the analytical solution of Fick's second law when considering the parameters of Equation :
where L is the half thickness of the slab (m); Def is the effective diffusion coefficient (m2 s-1); and n is the number of equation parameters.
Determining the coefficient (R2) for non-linear models is not a good decision-making tool. It is necessary to analyze different statistical parameters. Thus, for model adequacy in describing a phenomenon, P values must be taken into consideration; the lowest values of P indicate a higher representation of the model for the phenomenon under study. Analyzing SE also helps to select models. Lower SE values indicate higher model reliability. Therefore, by analyzing these three statistical parameters shown in Table 2, it was concluded that the two-term model is the best one to represent infrared drying of tomato fruit.
No significant difference was observed among tomato maturity stages at a 5% probability level by the Tukey test. As a result, Figure 1 represents the tomato drying curve. This figure shows the drying values observed and estimated by the two-term model at maturity stage 1. By analyzing Figure 1, the good adjustment of the two-term model to the experimental data can be observed, which satisfactorily represents tomato drying kinetics.
Thus, it is possible to estimate the moisture content of tomato during infrared drying at 60 °C by Equation :
MR= 0.8952 exp(-0.1036t) + 0.1503exp(-0.0064t) 
Figure 1 shows that some linearity occurs at approximately 10 min after the moisture ratio is reduced followed by an exponential reduction. The critical moisture content is defined at this point. This parameter divides the drying curve into two parts: before this point, it can be inferred that there is a constant falling rate period in which the velocities of moisture removal at the periphery and moisture replacement from the center to the periphery are considered to be equal; after this point of critical moisture content, a decrease in the moisture ratio with lower velocity is observed in which moisture replacement previously observed cannot supply the removal of moisture by infrared drying at the same velocity.
By analyzing the linearity of experimental data, a critical moisture content of 2.97 kgw kgdm-1 was found for tomato dehydration; this corroborates Brooker et al. (1992), who stated that the critical moisture content for the dehydration of biological products presenting a high initial moisture content ranges from 2.33 to 3.00 kgw kgdm-1.
Figure 2 reports drying rate vs. drying time of tomato fruits dried at 60 °C. It can be observed in Figure 2 that there is an initial period in which the drying rate reaches its maximum value (approximately 1.05 kgw kgdm-1, about 3 min). The product is being heated during this period of time. A decrease follows this increase. According to Celma et al. (2009), the initial evaporation occurs on the sample surface. However, this fact becomes less important over time. Therefore, moisture diffusion becomes the main factor affecting the drying rate. This trend was also observed for other agricultural products, such as okra (Abelmoschus esculentus (L.) Moench; Doymaz, 2005b) and olive cake (Akgun and Doymaz, 2005).
For long dehydration periods (MR < 0.6), the diffusion equation can be simplified by a one-term series (Senadeera et al., 2003). Therefore, assuming that De/is independent of moisture content, Equation  can be rewritten as follows:
The effective diffusion coefficient can thus be calculated by replacing the experimental data in Equation . The effective diffusion coefficient is determined by plotting ln MR vs. drying time, in which the angular coefficient of the linear regression is π2Def/4L2.
According to Karathanos et al. (1990), the effective diffusion coefficient of biological materials characterizes the transport of intrinsic mass in the form of water molecules, including molecular diffusion, liquid diffusion, gaseous diffusion, hydrodynamic flow, and other mechanisms. Values of Def for tomato fruits dried by infrared at maturity stages 1, 2, and 3 were 3.33, 2.78, and 3.05 x 10-7 m2 s-1, respectively
Figure 3 shows the plotting of ln MR vs. drying time of tomato fruits at maturity stage 1. Although the relationship between MR and drying time exhibited a high determination coefficient (R2 = 97.99%), it can be observed in Figure 3 that there is a nonlinear tendency of this relationship in the final stage of drying. Celma et al. (2009) reported that the effective diffusion coefficient depends on moisture content when this nonlinearity is observed. These authors reported a relationship between De/ and moisture content represented by Equations  and :
where F0 is the Fourier number, dimensionless.
Values were found by this method for the effective diffusion coefficient of 3.75, 3.42, and 3.80 x 10-7 m2 s-1 at maturity stages 1, 2, and 3, respectively. The values of these two methods were different and the ANSYS 5.3 (ANSYS, Canonsburg, Pennsylvania, USA) software system was used to detect which method showed the lowest error using the finite element method.
The finite element method (FEM), with its optimization routines, has been recently applied because it is easy to use for complex models and geometries (Gastón et al., 2002) . This methodology was used to determine the effective diffusion coefficient of wheat with spherical and ellipsoidal geometries (Gastón et al., 2002); to determine Def values of different apple tissues (Veraverbeke et al., 2003) ; and to obtain Def values during the soaking of rice by considering the expansion feature (Bakalis et al., 2009). Therefore, Def values were found by using the two abovementioned methods and a geometry representing the shape of tomato samples.
Moisture transfer within tomato fruit was close to the liquid diffusion theory without generating water inside the control volume and is represented by Equation . The initial condition is stated by Equation  by considering that water distribution at the beginning of the process is homogeneous. The boundary condition is represented by Equation  by considering the water vapor mass equilibrium between the product and the drying air:
where Ω is the control volume and r is the contour surface.
The equilibrium moisture content was calculated by the two-term equation. The geometry of the product was approximated to a sphere with an equivalent diameter of 70.55 mm. This diameter was calculated based on the mean of characteristic dimensions obtained with a 0.01 mm precision digital caliper. The analysis was carried out with a 45° angle semi-circle because of geometric symmetry.
A mesh with quadrangular elements and axis symmetry at the surface control was created. Geometry, mesh generation, and model implementation were carried out with the ANSYS 5.3 software system. Simulation results were compared with experimental data in each effective diffusion coefficient. The De/ values were selected with the least square sum of deviations (Equation ) between experimental and simulatedndata.
where RSD is the least square sum of deviations; Yj is the observed value; Yj is the average value predicted by the model; and n is the number of observed values.
Table 3 shows the values of the effective diffusion coefficient at each tomato maturity stage and the values of RSD obtained for each method. The RSD values ranged from 0.0129 to 0.0197 moisture content d.b. Despite these low values none of the methods were able to effectively predict the drying curve (Figure 4). Therefore, an optimization procedure of Def values was simulated with the ANSYS 5.3 software system. The Def value was 1.00 x 10-7 m2 s-1 (RSD = 0.0009 d.b.).
Table 4 shows a wide variation among products and among cultivars of the same product. This trend is probably due to differences in the chemical composition, physical structure, and geometry of the products; these values can vary with boundary conditions, drying methods, and prediction methods of Def values (Campos et al., 2009; Goneli et al., 2009; Corrêa et al., 2010c).
Mass transfer coefficients
The global coefficient of heat transfer (hc) was 12.45 W m-2 K-1 and the global coefficient of mass transfer (hm) 0.0105 m s-1. The values of hm in the present work corroborate the values reported by Touré and Kibangu-Nkembo (2004), who found 0.00974 m s-1 for natural dehydration of mango fruits with solar energy, 0.00475 m s-1 for cassava roots, and 0.00876 and 0.00943 m s-1 for two banana cultivars, respectively. Furthermore, the values of hc were expected since, according to Bird et al. (2004), the values of the global coefficient of heat transfer for free convection is between 3 and 20 W m-2 K-1.
Dehydration kinetics by infrared drying of tomato slices can be similar to the theoretical studies on heat and mass transfer mechanisms. The values found for the global coefficient of heat and mass transfer were 12.45 W m-2 K-1and 0.0105 m s-1, respectively. An effective diffusion coefficient of 1.00 x 10-7 m2 s-1 was found through an optimization procedure using the finite element method. The two-term model was the best to represent the tomato dehydration process, and the critical moisture content for tomato dehydration was 2.97 kgw kgdm-1. The drying rate reached its maximum value (1.05 kgw kgdm-1) at 3 min of dehydration.
The authors thank FAPEMIG and CNPq for funding, which was indispensable to carry out this work.
Akgun, N.A., and I. Doymaz. 2005. Modelling of olive cake thin-layer drying process. Journal of Food Engineering 68:455-461. [ Links ]
Bakalis, S., A. Kyritsi, V.T. Karathanos, and S. Yanniotis. 2009. Modeling of rice hydration using finite elements. Journal of Food Engineering 94:321-325. [ Links ]
Bird, R.B., W.E. Stewart, and E.N. Lightfoot. 2004. Fenomenos de transporte. 2a ed. 838 p. LTC Rio de Janeiro, Rio de Janeiro, Brazil. [ Links ]
Brooker, D.B., F.W. Bakker-Arkema, and C.W. Hall.1992. Drying and storage of grains and oilseeds. The AVI Publishing Company, Westport, Connecticut, USA. [ Links ]
Caglar, A., I.T. Togrul, and H. Togrul. 2009. Moisture and thermal diffusivity of seedless grape under infrared drying. Food and Bioproduct Process 87:292-300. [ Links ]
Campos, S.C., P.C. Corrêa, F.M. Botelho, B.L. Nogueira, G.H.H. Oliveira, and A.A. Paixao. 2009. Avaliacao da secagem dos graos de café de diferentes etapas do processamento via úmida. Revista Brasileira de Armazenamento. Especial Café. p. 88-98.
Celma, A.R., F. Cuadros, and F. Lópes-Rodriguez. 2009. Characterisation of industrial tomato by-products from infrared drying process. Food and Bioproducts Processing 89:282-291. [ Links ]
Christensen, C.M., and H.H. Kaufmann. 1974. Storage of cereal grain and their products. p. 158-192. American Association of Cereal Chemists, St. Paul, Minnesota, USA. [ Links ]
Corrêa, P.C., A.L.D. Goneli, P.C. Afonso Júnior, G.H.H. Oliveira, and D.S.M. Valente. 2010b. Moisture sorption isotherms and isosteric heat of sorption of coffee in different processing levels. International Journal of Food Science and Technology 45:2016-2022. [ Links ]
Corrêa, P.C., G.H.H. Oliveira, F.M. Botelho, A.L.D. Goneli, and F.M. Carvalho. 2010c. Modelagem matemática e determinacao das propriedades termodinàmicas do café (Coffea arabica L.) durante o processo de secagem. Revista Ceres 57:595-601.
Corrêa, P.C., G.H.H. Oliveira, A.P.L. Rodrigues, S.C. Campos, and F.M. Botelho. 2010a. Hygroscopic equilibrium and physical properties evaluation affected by parchment presence of coffee grain. Spanish Journalof Agricultural Research 8:694-702. [ Links ]
Diamente, L.M., and P.A. Munro. 1993. Mathematical modelling of the thin layer solar drying of sweet potato slices. Solar Energy 51:271-276. [ Links ]
Doymaz, I. 2005a. Sun drying of figs: An experimental study. Journal of Food Engineering 71:403-407. [ Links ]
Doymaz, I. 2005b. Drying characteristics and kinetics of okra. Journal of Food Engineering 69:275-279. [ Links ]
Gastón, A.L., R.M. Abalone, and S.A. Giner. 2002. Wheat drying kinetics. Diffusivities for sphere and ellipsoid by finite elements. Journal of Food Engineering 52:313-322. [ Links ]
Goneli, A.L.D., P.C. Corrêa, P.C. Afonso Júnior, and G.H.H. Oliveira. 2009. Cinética de secagem dos graos de café descascados em camada delgada. Revista Brasileirade Armazenamento. Especial Café. p. 64-73. [ Links ]
Goneli, A.L.D., P.C. Corrêa, G.H.H. Oliveira, and F.M. Botelho. 2010a. Water desorption and thermodynamic properties of okra seeds. Transactions of ASABE 53:191-197. [ Links ]
Goneli, A.L.D., P.C. Corrêa, G.H.H. Oliveira, C.F. Gomes, and F.M. Botelho. 2010b. Water sorption isotherms and thermodynamic properties of pearl millet grain. International Journal of Food Science and Technology 45:828-838. [ Links ]
Hall, C.W. 1980. Drying and storage of agricultural crops. 381 p. AVI Publishing Company, Westport, Connecticut, USA. [ Links ]
Henderson, S.M. 1974. Progress in developing the thin-layer drying equation. Transactions of the ASAE 17:1167-1168. [ Links ]
Henderson, S.M., and S. Pabis. 1961. Grain drying theory. I. Temperature effect on drying coefficient. Journal of Agriculture Engineering Research 6:169-174. [ Links ]
Karathanos, V.T., G. Villalobos, and G.D. Saravacos. 1990. Comparison of two methods of estimation of the effective moisture diffusivity from drying data. Journal of Food Science 55:218-223. [ Links ]
Midilli, A., H. Kucuk, and Z. Yapar. 2002. A new model for single- layer drying. Drying Technology 20:1503-1513. [ Links ]
Mota, C.L., C. Luciano, A. Dias, M.J. Barroca, and R.P.F. Guiné. 2010. Convective drying of onion: kinetics and nutritional evaluation. Food and Bioproducts Processing 88:115-123. [ Links ]
Oliveira, G.H.H. 2010. Características físico-químicas do tomate submetido a esforco controlado de compressao. Tesis. Universidade Federal de Vicosa, Vicosa, Minas Gerais, Brazil.
Oliveira, G.H.H., P.C. Corrèa, E.F. Araújo, D.S.M. Valente, and F.M. Botelho. 2010. Desorption isotherms and thermodynamic properties of sweet corn cultivars (Zea mays L.) International Journal of Food Science and Technology 45:546-554. [ Links ]
Sacilik, K., A.K. Elicin, and G. Unal. 2006. Drying kinetics of Uryani plum in a convective hot-air dryer. Journal of Food Engineering 76:362-368. [ Links ]
Senadeera, W., B.R. Bhandari, G. Young, and B. Wijesinghe. 2003. Influence of shapes of selected vegetable materials on drying kinetics during fluidized bed drying. Journal of Food Engineering 58:277-283. [ Links ]
Sharma, G.P., Verma, R.C., and P.B. Pathare. 2005. Thin-layer infrared radiation drying of onion slices. Journal of Food Engineering 67:361-366. [ Links ]
Touré, S., and S. Kibangu-Nkembo. 2004. Comparative study of natural solar drying of cassava, banana and mango. Renewable Energy 29:975-990. [ Links ]
Vasques, A.R., S.L. Bertoli, R.C.S.C. Valle, and J.A.B. Valle. 2006. Avaliacao sensorial e determinacao de vida-de-prateleira de macas desidratadas. Ciencia e Tecnologia de Alimentos 26:759-765.
Veraverbeke, E.A., P. Verboven, N. Scheerlinck, M.L. Hoang, and B.M. Nicolai'. 2003. Determination of the diffusion coefficient of tissue, cuticle, cutin and wax of apple. Journal of Food Engineering 58:285-294. [ Links ]
Verma, L.R., R.A. Bucklin, J.B. Endan, and F.T. Wratten. 1985. Effects of drying air parameters on rice dryingmodels. Transaction of ASAE 28:296-301. [ Links ]
Welti-Chanes, J., F. Vergara-Balderas, and D. Bermúdez-Aguirre. 2005. Transport phenomena in food engineering: Basic concepts and advances. Journal of Food Engineering 67:113-128. [ Links ]
Yagcioglu, A., A. Degirmencioglu, and F. Cagatay. 1999. Drying characteristic of laurel leaves under different conditions. p. 565569. In Proceedings of 7th International Congress on Agricultural Mechanization and Energy, Adana, Turkey. [ Links ]
Yaldiz, O., and C. Ertekin. 2001. Thin layer solar drying some different vegetables. Drying Technology 19:583-597. [ Links ]
Received: 5 August 2011.
Accepted: 18 May 2012.