Chilean Journal of Agricultural Research 69(2):171-178 (April-June 2009)

**RESEARCH**

**Kinetic Study of Convective Drying of Blueberry Variety O'neil ( Vaccinium corymbosum L.) Estudio de la Cinética del Secado Convectivo de Arándano Variedad O'Neil (Vaccinium corymbosum L.) A. Vega-Gálvez^{1}*, R. Lemus-Mondaca^{1}, C. Tello-Ireland^{1}, M. Miranda^{1}, and F. Yagnam^{1 }**

^{1}Universidad de La Serena, Departamento de Ingeniería de Alimentos, Av. Raúl Bitrán s/n. Box 599, La Serena, Chile.*Corresponding author (avegag@userena.cl) ]]> Received: 04 January 2008. Accepted: 14 April 2008.

**ABSTRACT **The aim of this research was to study and to model the drying kinetics of the blueberry

*Vaccinium corymbosum*L.) at three temperatures (60, 70 and 80 ºC) with an airflow of 2.0 ± 0.2 m s

^{-1}. Modeling of the desorption isotherm was carried out with the GAB (Guggenheim, Anderson and de Boer) equation, showing a good fit to experimental moisture data, giving as a result a monolayer moisture level of 0.084 g water g

^{-1}dm. Newton, Henderson-Pabis, Page, Modified Page and Logarithmic mathematical models were applied in the study and in the modeling of the drying kinetics of this fruit. Kinetic parameters

*k*of each model showed dependence on temperature, and were evaluated by an Arrhenius-type equation, with an activation energy of between 36.2 and 54.5 kJ mol

^{-1}. Logarithmic and Modified Page models gave the best fits for each drying curve, based on the statistical test determination coefficient, sum square error, root mean sum errors and Chi-square. In consequence, both models are excellent tools for estimating the drying time of this product.

**Key words:**blueberry, GAB, drying, modeling, statistical tests.

**RESUMEN **El objetivo de esta investigación fue estudiar y modelar la cinética de secado del arándano (

*Vaccinium corymbosum*L.) a tres temperaturas (60, 70 y 80 ºC) con un flujo de aire de 2,0 ± 0,2 m s

^{-1}. El modelado de la isoterma de desorción se llevó a cabo con la ecuación de GAB (Guggenheim, Anderson y de Boer), mostrando un buen ajuste sobre los datos experimentales de humedad, dando como resultado una humedad de la monocapa de 0,084 g de agua g

^{-1}ms. Se aplicaron los modelos matemáticos de Newton, Henderson-Pabis, Page, Page modificado y Logarítmico para el modelado de la cinética de secado de esta fruta. Los parámetros cinéticos

*k*de cada modelo presentaron dependencia con la temperatura, evaluadas por una ecuación de tipo Arrhenius, con una energía de activación entre 36,2-54,5 kJ mol

^{-1}. Los modelos Logarítmico y Page modificado obtuvieron el mejor ajuste para cada curva de secado, basado en las pruebas estadísticas como coeficiente de determinación, suma de errores cuadrados, raíz media de los errores cuadrados y Chi-cuadrado. En consecuencia, ambos modelos son excelentes herramientas para estimar el tiempo de secado de este producto.

]]>

**Palabras clave:**arándanos, GAB, secado, modelado, pruebas estadísticas.

**INTRODUCTION **The blueberry variety O’Neil belongs to the genus

*Vaccinium*, of the family

*Ericaceae*, native of North America. It is in the product grouping of berries that includes strawberry, blackberry, raspberry, and others (USHBC, 2007). It was introduced in Chile in 1990 with very good adaptive results, with the result that Chile is now considered as the main producing country of this product in the Southern Hemisphere and the third largest producer worldwide. Of the total production, 85 to 90% is exported, mainly as fresh-cooled, and a small amount is shipped frozen (Chilealimentos, 2007). The production that does not meet whole-berry quality standards is generally converted to clarified and/or concentrated juice, and other technological alternatives are being explored in order to successfully market second quality material, including convection drying with previous pretreatment (osmo-drying, high pressure drying, etc.).

Blueberries are almost spherical, depending on the species and cultivation conditions, and between 0.7 and 1.5 cm in diameter, with a dark blue color. They contain up to 100 very small seeds in a central core (Stückrath and Petzold, 2007). In Chile, blueberries mature between December and late January, depending on the cultivation zone, and extending over a period of 4 to 5 wk. These berries are a rich source of antioxidant phytonutrients, which are believed to be associated with the antioxidant activity of anthocyanin pigments, flavonoids, and other phenolic compounds containing about three times the total phenolic compounds found in other berries (Skrede

*et al*., 2000). Blueberries have become very popular with consumers because of the research findings that associate their consumption with improvements in human health (Nindo

*et al*., 2007). However, seasonality, market accessibility and costs commonly limit the consumption of fresh fruits and vegetables, so dehydrated products are preferred over fresh products for their convenience, availability and shelf life (Azzouz

*et al*., 2002).

Hot air drying is one of most widely used methods for food preservation. The advantage of dehydrated foods is that decreased moisture content reduces thermodynamic water activity, thus preventing the growth of microorganisms that cause spoilage reactions (Babalis and Belessiotis, 2004). The optimization of dehydration processes in the agro-food industry has led to choosing the technological variables involved in the process itself. Good experimental designs, along with statistical programs, help to obtain a higher yield from the operational and capital points of view. This must be done considering a product of the highest quality, depending on the present variables, one of the most important being the time of the drying process. Prior to the study of the drying of any food, it is necessary to recognize and evaluate the moisture sorption isotherms, as these mathematically describe the relation between the water activity and the equilibrium moisture content of the product under study (Vega-Gálvez

*et al*., 2008). The use of mathematical equations is necessary for the simulation of the kinetics of matter transfer (water) that occurs during this unit operation (Vega

*et al*., 2007). There are several empirical equations used to simulate the drying process that are exceptionally practical for the study, the modeling of kinetics and the process optimization, as well as for dryer design (Senadeera

*et al*., 2003). Notably among these equations are those proposed by Newton, Henderson-Pabis, Page and Modified Page, etc. (Doymaz, 2004; Akpinar and Bicer, 2006).

The aim of the present study was to determine and model the drying kinetics and desorption isotherm of the blueberry variety O’Neil, using mathematical equations for both phenomena, and to evaluate the influence of drying air temperature on the kinetic parameters. ]]>

**MATERIALS AND METHODS**

**Raw material and proximate analysis**

Blueberries of the variety O’Neil (

*Vaccinium corymbosum*L.) were cultivated and purchased in the province of Salamanca, Chile. Samples were selected to provide a homogeneous group, based on their date of harvest, color, size, and freshness according to visual analysis. The moisture content was determined according to AOAC methodology Nº. 934.06 (AOAC, 1990), using a vacuum oven (Gallenkamp, OVL570, Leicester, UK) at 70 ºC for 72 h, and an analytical balance (CHYO, Jex120, Kyoto, Japan) with an accuracy of ± 0.0001 g. Crude protein content was determined using the Kjeldahl method with a conversion factor of 6.25. Lipid content was analyzed gravimetrically following Soxhlet extraction. Crude fiber was estimated by acid/alkaline hydrolysis of insoluble residue. Crude ash was estimated by incineration in a muffle furnace at 550 ºC. Acidity was determined by the adapted AOAC methodology Nº. 942.15A (AOAC, 1990), pH was measured using a potentiometer (Extech Instruments, Microcomputer pH-Vision 246072, Waltham, Massachusetts, USA), and sugar content was measured using an Abbe refractometer (ATAGO, 1-T, Tokyo, Japan). All the analyses were made in triplicate and expressed in g 100 g

^{-1}sample.

**Isotherm experiments**

Desorption isotherms were measured at 60 ºC. A known mass of sample (in triplicate) was allowed to come to equilibrium with the atmosphere (relative humidity) inside a hermetically sealed flask, which contained a glass dish with a saturated salt solution of known water activity. This standard gravimetric method was recommended by The European Cooperative Project COST 90, which deals with the physical properties of foods (Spiess and Wolf, 1983). The weight of the samples was taken every 15 days until reaching constant weight (equilibrium condition). The salts used to obtain a range of water activity of 0.10 to 0.95 included LiCl, KC

_{2}H

_{3}O

_{2}, MgCl

_{2}, K

_{2}CO

_{3}, NaNO

_{2}, KI, NaCl and KNO

_{2}(Lim

*et al*., 1995). Thymol was added separately in a Petri dish to the recipients containing saturated salt solutions with a relative humidity higher than 75% in order to avoid microbial growth, especially mould (Vega

*et al*., 2007). Once equilibrium was reached, the moisture content of the samples was determined in triplicate.

The relationship between equilibrium moisture content and water activity of blueberries was expressed by the equation proposed by Guggenheim, Anderson and de Boer, commonly termed GAB (Equation [1], Quirijns

*et al.*, 2005; Blahovec and Yanniotis, 2008). GAB is commonly used in moisture sorption experiments in different foods, because of its important parameters and physiochemical description, as monolayer moisture (X

_{m}) and the parameters C

_{o}and K

_{o}(Yu

*et al*., 1999; Timmermann

*et al*., 2001). These parameters were obtained by nonlinear regression analysis, using the EXCEL

^{®}program of Microsoft

^{®}Windows

^{®}XP (Redmond, Washington, USA). ]]>

The criteria to evaluate the fit quality of GAB model were the statisticals determination coefficient (r

^{2}) and mean percent square error (%E) (Equation [2]) (Vega

*et al*., 2007), where X

_{we}is the equilibrium moisture content (g water g

^{-1}dm), a

_{w}is the water activity (dimensionless), X

_{m}monolayer moisture (g water g

^{-1}dm), C

_{o}and K

_{o}are constants of GAB model (dimensionless), X

_{ei}is the experimental moisture content (g water g

^{-1}dm), X

_{ci}is the calculated moisture content (g water g

^{-1}dm), N is the number of data values and

*i*is the number of terms.

** Drying experiments **The drying experiment was carried out using a convective dryer designed and built in the Faculty of Engineering of Universidad de La Serena (Figure 1), La Serena. Three temperatures were used in the study of the drying kinetics (60, 70 and 80 ºC). Drying air velocity was held constant at 2.0 ± 0.2 m s

^{-1}and measured with an omnidirectional anemometer (Extech Instrument Inc., 451112, Waltham, Massachusetts, USA). The inlet relative humidity was 62.0 ± 5.2%, measured by an ambient digital hygro-thermometer (Extech Instrument Inc., 445703, Waltham, Massachusetts, USA). All the drying experiments were carried out in triplicate, using a sample mass of 100.0 ± 2.4 g and the charge density was 8.5 ± 0.3 kg m

^{-2}. Prior to the drying of the samples, they were pretreated with a Pectinex

^{®}solution, concentration 0.8% (Novo Nordisk Ferment Ltd., Flawil, Switzerland). The mass was measured on an analytical balance (Ohaus, SP402, New Jersey, USA) with an accuracy of ± 0.01 g at defined time intervals, connected by a interface system (Ohaus, RS232, Pine Brook, New Jersey, USA) to a PC, which recorded and stored the data. The experiments were finished at the point of reaching constant weight (equilibrium condition). The dried samples were packaged in polypropylene bags.

**Figure 1. Schematic diagram of drying equipment. ]]>
**

**Modeling of drying kinetics **Drying kinetic was modeled by means of five empirical equations widely used in most organic and biological materials. These equations use a relationship termed moisture ratio (MR) as a dependant variable (Eq. 3), relating the gradient of the sample moisture in real time (X

_{wt}) with initial moisture (X

_{wo}) and equilibrium moisture (X

_{we}) (Babalis and Belessiotis, 2004). These are Newton (Eq. 4), Henderson-Pabis (Eq. 5), Page (Eq. 6), Modified Page (Eq. 7) and Logarithmic (Eq. 8). In this research, the shrinkage and external resistance were assumed as negligible. Where

*k*is kinetic parameters (min

^{-1}),

*n*and

*a*are empirical parameters (dimensionless) and t is drying time (min).

In order to observe any influence of drying temperature on the kinetic parameters

*k*, an Arrhenius-type equation was applied (Equation [9]), from which the activation energy is obtained (E

_{a}kJ mol

^{-1}), which shows sensitivity of the parameter to temperature (Simal

*et al*., 1996; Simal

*et al*., 2005). Activation energy can be determined by the graphic representation between ln

*k*versus T

^{-1}(K).

**
]]>
**

**Statistical analysis of the models**

The fit quality of the experimental data to the desorption isotherm and all the models proposed for drying kinetics were evaluated using the determination coefficient (r

^{2}), sum squared errors (SSE, Equation [10]), root mean sum errors (RMSE, Equation [11]) and Chi-square (χ

^{2}, Equation [12]) statisticals. The values closest to 1.0 for r

^{2}, and those closest to zero for SSE, RMSE and

*X*

^{2}, are commonly considered as optimum criteria to evaluate the fit quality of the models used (Doymaz, 2004; Akpinar and Bicer, 2006; Doymaz, 2007; Vega

*et al*., 2007). The statistical evaluations were made on the predictions of equilibrium moisture content by the desorption isotherm, as well as on the drying kinetics. Where MR

_{ej}is the experimental moisture ratio (dimensionless), MR

_{cj}is the calculated moisture ratio (dimensionless), z is the number of constants and j is the number of terms.

**RESULTS AND DISCUSSION **

**Experimental analysis and desorption isotherm**

The proximate analysis of blueberry variety O’Neil gave an initial moisture content of 78.13 ± 0.05 g 100 g

^{-1}sample (a

_{w}= 0.99 ± 0.01); crude protein (nitrogen x 6.25) of 0.62 ± 0.12 g 100 g

^{-1}sample; total lipids of 0.40 ± 0.04 g 100 g

^{-1}sample; crude fiber of 6.51 ± 0.27 g 100 g

^{-1}sample; ash of 1.20 ± 0.03 g 100 g

^{-1}sample; non-nitrogen extract (by difference) of 12.90 g 100 g

^{-1}sample; acidity of 2.21 ± 0.12% (monohydrated citric acid); pH 2.72 ± 0.09, and soluble solids of 15.01 ± 0.07 ºBrix. ]]>

Figure 2 shows the experimental data of the equilibrium moisture contents of the blueberry desorption isotherm at 60 ºC, where equilibrium moisture content increases as water activity increases from 0.10 to 0.95. The tendency of this isotherm could correspond to that of type II isotherms, according to Van der Waals' classification (Brunauer

*et al*., 1938). Various authors working with other types of foods have observed this behavior (Kiranoudis

*et al*., 1993; Timmermann

*et al*., 2001; Kaymak-Ertekin and Gedik, 2004). The type II isotherm appears when the bonding energy between the water and the primary layer is less than that occurring among water molecules (Lomauro

*et al*., 1985). It was observed the good fit on the experimental data obtained by GAB (r

^{2}= 0.97; %E = 8.64) for a whole range of water activity (Figure 2). The sorption parameters obtained were X

_{m}= 0.084 g water g

^{-1}dm, C

_{o}= 13.319, and K

_{o}= 0.933. Similar results of monolayer moisture have been observed in raisins, figs, plums, potatoes, onions and tomatoes with 0.212-0.087 g water g

^{-1}dm for 30-60 ºC (Kiranoudis

*et al*., 1993); in red and green peppers with 0.113-0.038 g water g

^{-1}dm for 30-60 ºC (Kaymak-Ertekin and Sultanoğlu, 2001); in grapes, apples, potatoes and apricots with 0.220-0.095 g water g

^{-1}dm for 30-60 ºC (Kaymak-Ertekin and Gedik, 2004). Monolayer moisture (X

_{m}) is an important parameter since it has a physicochemical behavior that represents the first layer of water molecules, which can thermodynamically interact with other food compounds (Lim

*et al*., 1995; Yu

*et al*., 1999).

**Figure 2. Desorption isotherm of blueberries modeled by the GAB equation at 60 ºC. **

**Behavior of drying curves **For all drying experiments, an average outlet temperature of drying air of 60.0 ± 4.2 ºC was obtained. For this reason, the desorption isotherm modeled by the GAB equation was used to estimate the equilibrium moisture content to each temperature, with X

_{we}= 0.031 g water g

^{-1}dm for 60 ºC; X

_{we}= 0.027 g water g

^{-1}dm for 70 ºC and X

_{we}= 0.026 g water g

^{-1}dm for 80 ºC. In general, for all drying experiments, the equilibrium moisture content was lower than 18%, which gives product stability from the commercial and hygienic points of view (Karathanos and Belessiotis, 1999). Increasing the drying temperature decreased drying time (Figure 3). The shortest time (500 min) was obtained at 80 ºC in comparison to drying at 70 and 60 ºC, which required times of 800 and 1400 min, respectively. Other research on fruit and vegetable drying present the same tendency and behavior of the drying curves, such as Azzouz

*et al*. (2002); Krokida

*et al*. (2003); Babalis and Belessiotis (2004); Simal

*et al*. (2005); Akpinar and Bicer (2006); Vega

*et al*. (2007); and Doymaz (2007).

**Figure 3. Drying curves for blueberries at different drying temperatures (60, 70, and 80 ºC). ]]>
**

**Modeling of the drying process **Drying curves (Figure 3) showed a clear exponential tendency and only a falling rate period was observed, which is very common in fruit and vegetable drying processes (Toğrul and Pehlivan, 2003; Kingsly

*et al*., 2007). Under these conditions, the use of the five empirical models is suggested, as mentioned before. Table 1 shows the values for the kinetic parameters

*k*of the five models for each drying temperature evaluated. A value of p < 0.05 for a confidence level of 95% was obtained by ANOVA, using Statgraphics Plus

^{®}5.1 software (Statistical Graphics Corp., Herndon, Virginia, USA), suggesting there are statistically significant differences among these kinetic parameters with respect to temperature. In order to prove the dependence of these parameters on the drying temperature, an Arrhenius-type equation (Figure 4) was applied, showing r

^{2}≥ 0.95. The activation energy obtained for each kinetic parameter was 51.05 kJ mol

^{-1}(Newton), 54.45 kJ mol

^{-1}(Henderson-Pabis), 36.20 kJ mol

^{-1}(Page), 46.39 kJ mol

^{-1}(Modified Page), and 48.34 kJ mol

^{-1}(Logarithmic). Several authors have presented very similar activation energy values, including Azzouz

*et al*. (2002); Babalis and Belessiotis (2004); Simal

*et al*. (2005); Akpinar and Bicer (2006), Doymaz (2007) and Vega

*et al*. (2007). Kinetic parameter

*k*(min

^{-1}) also showed a clear tendency of increasing as the working temperature increased (Table 1). Similar results were obtained by Karathanos and Belessiotis (1999) working with figs, plums and raisins; Toğrul and Pehlivan (2003) with apricots; and Akpinar and Bicer (2006) in strawberry drying.

**Table 1. Values of kinetic parameters k (x10^{-1} min^{-1}) for each drying curve. **

**Figure 4. Relationship between the kinetic parameters k versus the inverse of temperature (K) modeled by means of five empirical equations: Newton, Henderson-Pabis (H-P), Page, Modified Page (mod. Page) and Logarithmic (Log). **

**The empirical parameters**

*n*of Henderson-Pabis, Page, Modified Page and Logarithmic, as well as the empirical parameter

*a*, did not show statistically significant differences (p value > 0.05), suggesting they probably depend more on the characteristics of the tissue and the drying air flow (Akpinar and Bicer 2006; Menges and Ertekin 2006). Azzouz

*et al*. (2002), working with grapes, concluded that the parameter

*n*was in function of air flow rate and that the parameter

*k*of Page depended on the temperature and the initial moisture of the product. Karathanos and Belessiotis (1999), working on skinned and non-skinned fruit dehydration, proposed that parameter

*n*increased with the existence of the outer skin depending on its thickness and the kind of product to be dried.

**Statistical analyses of models**

Logarithmic and Modified Page models provided good fits to the experimental data for a whole drying process (Figure 5). All the models presented high values for the determination coefficient (r

^{2}≥ 0.95) at three temperatures (60, 70 and 80 ºC). Furthermore, the lowest SSE, RMSE and c

^{2}values were selected as optimal criteria in order to evaluate the fitting quality of the five models proposed. A good fit was observed based on this evaluation, since low SSE (< 0.0027), RMSE (< 0.0495) and c

^{2}(< 0.0029) values were obtained by Logarithmic and Modified Page models; followed by Henderson-Pabis, Page and Newton. Other authors have also obtained good results when applying these models in drying kinetics of other food and foodstuff (Krokida

*et al*., 2003; Toğrul and Pehlivan, 2003; Doymaz, 2004; Simal

*et al*., 2005; Menges and Ertekin, 2006; Doymaz, 2007).

**Figure 5. Graphic representation of the statistic test for each model. **

Figure 6 shows the experimental and calculated values for the drying curves represented by MR vs. time for the two best models found in this investigation (Logarithmic and Modiﬁed Page). Both models gave similar results over the whole drying process from the beginning to the end, considering the middle stage of the drying process. The good estimations obtained by these two models represent results that are not usually given by many models, since the middle stage of the drying process is the segment where most of the water is removed from the food. Thus, a good simulation is required for the estimation of the optimum drying time.

]]>**Figure 6. Experimental and calculated drying curves for a) Logarithmic (Log) and b) Modified Page (Mod. Page) at the different working temperatures (ºC).**

**CONCLUSIONS **The results of this study showed that the GAB model provided a good fit to the experimental data of desorption isotherm, with 8.64 %E and r

^{2}= 0.97. Drying of blueberries presents a clear dependence on drying air temperature, showing only a falling rate period, and reaching an average equilibrium moisture close to 0.03 g water g

^{-1}dm. All models used to describe the dehydration kinetics were useful. Nevertheless, Logarithmic and Modified Page models gave the best fit quality to drying experimental data at the three temperatures used, based on the statistical tests used for evaluation. All the kinetic parameters

*k*were dependant on the drying temperature, giving an activation energy of 48.34 and 46.39 kJ mol

^{-1}for Logarithmic and Modified Page, respectively. In consequence, both models are excellent tools for estimating the drying time of this product.

**ACKNOWLEDGEMENTS**

]]>The authors wish to thank the Research Office of the Universidad de La Serena (

]]>

*DIULS*) and Ricardo León Alcarraz of Agrícola Estero Camisas for the financial support given to the publication of this research.

**LITERATURE CITED**

AOAC. 1990. Official method of analysis. 15

^{th}ed. Association of Official Analytical Chemists (nº 934.06), Washington D.C., USA. [ Links ]

Akpinar, E.K., and Y. Bicer. 2006. Mathematical modeling experimental study on thin layer drying of strawberry. Int. J. Food Eng. 2:5. [ Links ]

Azzouz, S., A. Guizani, W. Jomaa, and A. Belghith. 2002. Moisture diffusivity drying kinetic equation of convective drying of grapes. J. Food Eng. 55:323-330. [ Links ]

]]>

Babalis, S.J., and V.G. Belessiotis. 2004. Influence of the drying conditions on the drying constants moisture diffusivity during the thin-layer drying of figs. J. Food Eng. 65:449-458. [ Links ]

Blahovec, J., and S. Yanniotis. 2008. GAB generalized equation for sorption phenomena. Food Bioprocess Technol. 1:82-90. [ Links ]

Brunauer, S., P. Emmett, and E. Teller. 1938. Adsorption of gases in multimolecular layers. J. Am. Chem. Soc. 60:309-319. [ Links ]

Chilealimentos. 2007. Asociación de empresas de alimentos de Chile. Productos y empresas. Ventajas de Chile. Available at http://www.chilealimentos.com. Accessed 10 June. [ Links ]

Doymaz, I. 2004. Drying kinetics of white mulberry. J. Food Eng. 61:341-346. [ Links ]

]]>

Doymaz, I. 2007. Influence of pretreatment solution on the drying of sour cherry. J. Food Eng. 78:591-596. [ Links ]

Karathanos, V.T., and V.G. Belessiotis. 1999. Application of a thin-layer equation to drying data of fresh semi-dried fruits. J. Agric. Eng. Res. 74:355-361. [ Links ]

Kaymak-Ertekin, F., and A. Gedik. 2004. Sorption isotherms and isosteric heat of sorption for grapes, apricots, apples and potatoes. LWT - Food Sci. Technol. 37:429-438. [ Links ]

Kaymak-Ertekin, F., and M. Sultanoğlu. 2001. Moisture sorption isotherm characteristics of peppers. J. Food Eng. 47:225-231. [ Links ]

Kingsly, R.P., R.K. Goyal, M.R. Manikantan, and S.M. Ilyas. 2007. Effects of pretreatments drying air temperature on drying behaviour of peach slice. Int. J. Food Sci. Technol. 42:65-69. [ Links ]

]]>

Kiranoudis, C.T., Z.B. Maroulis, E. Tsami, and D. Marino-Kouris. 1993. Equilibrium moisture content and heat of desorption of some vegetables. J. Food Eng. 20:55-74. [ Links ]

Krokida, M.K., V.T. Karathanos, Z.B. Maroulis, and D. Marinos-Kouris. 2003. Drying kinetics of some vegetables. J. Food Eng. 59:391-403. [ Links ]

Lim, L.T., J. Tang, and J. He. 1995. Moisture sorption characteristics of freeze dried blueberries. J. Food Sci. 60:810-814. [ Links ]

Lomauro, C.J., A.S. Bakshi, and T.P. Labuza. 1985. Evaluation of food moisture sorption isotherm equations. Part I: fruit, vegetable and meat products. LWT – Food Sci. Technol. 18:111-117. [ Links ]

Menges, H.O., and C. Ertekin. 2006. Mathematical modeling of thin layer drying of golden apples. J. Food Eng. 77:119-125. [ Links ]

]]>

Nindo, C.I., J. Tanga, J.R. Powersb, and P.S. Takhar. 2007. Rheological properties of blueberry puree for processing applications. LWT – Food Sci. Technol. 40:292-299. [ Links ]

Quirijns, E.J., A.J.B. van Boxtel, W.K.P. van Loon, and G. van Straten. 2005. Sorption isotherms, GAB parameters and isosteric heat of sorption. J. Sci. Food Agric. 85:1805-1814. [ Links ]

Senadeera, W., B.R. Bhari, G. Young, and B. Wijesinghe. 2003. Influence of shapes of selected vegetable material son drying kinetics during fluidized bed drying. J. Food Eng. 58:277-283. [ Links ]

Simal, S., A. Femenia, M.C. Garau, and C. Roselló. 2005. Use of exponential, Page’s and difusional models to simulate the drying kinetics of kiwi fruits. J. Food Eng. 66:323-328. [ Links ]

Simal, S., A. Mulet, P.J. Catalá, J. Cañellas, and C. Roselló. 1996. Moving boundary model for simulating moisture movement in grapes. J. Food Sci. 61:157-160. [ Links ]

]]>

Spiess, W., and W. Wolf. 1983. The results of the COST 90 project on water activity. p. 65. Jowitt, R., F. Escher, M. Kent. B. McKenna, and M. Roques (eds.) Physical properties of foods. Elsevier Applied Science Publisher, London, England. [ Links ]

Skrede, G., R.E. Wrolstad, and R.W. Durst. 2000. Changes in anthocyanins and polyphenolics during juice processing of O’Neil blueberries (

*Vaccinium corymbosum*L.). J. Food Sci. 65:357-364. [ Links ]

Stückrath, R., and G. Petzold. 2007. Formulación de una pasta gelificada a partir del descarte de arándanos (

*Vaccinium corymbosum*). Inf. Tecnol. 18:53-60. [ Links ]

Timmermann, E.O., J. Chirife, and H.A. Iglesias. 2001. Water sorption isotherms of foods and foodstuffs: BET and GAB parameters?. J. Food Eng. 48:19-31. [ Links ]

Toğrul, I.T., and D. Pehlivan. 2003. Modeling of drying kinetics of single apricot. J. Food Eng. 58:23-32. [ Links ]

]]>

USHBC. 2007. Blueberries storage and processing. United States Highbush Blueberry Council, Folsom, California, USA. Available at http://www.blueberry.org/Storage&Processing1.pdf (Accessed 21 May). [ Links ]

Vega-Gálvez, A., R. Lemus-Mondaca, C. Bilbao-Sainz, F. Yagnam, and A. Rojas. 2008. Mass transfer kinetics during convective drying of red pepper var. Hungarian

*Capsicum*

*annuum*L.): Mathematical modeling and evaluation of kinetics parameters. J. Food Process Eng. 31:120-137. [ Links ]

Vega, A., A. Andrés, P. Fito, and R. Lemus. 2007. Mathematical modeling of hot-air drying kinetics of red bell pepper (var. Lamuyo). J. Food Eng. 79:1460-1466. [ Links ]

Yu, L., G. Mazza, and D.S. Jayas. 1999. Moisture sorption characteristics of freeze-dried, osmofreeze-dried, and osmo-air-dried cherries and blueberries. Am. Soc. Agric. Eng. 42:141-147. [ Links ]