## Food Science and Technology

##
*Print version* ISSN 0101-2061*On-line version* ISSN 1678-457X

### Ciênc. Tecnol. Aliment. vol.22 no.3 Campinas Sept./Dec. 2002

#### http://dx.doi.org/10.1590/S0101-20612002000300021

**Air blast freezing of fruit pulp models in commercial boxes: influence of preferential channels in the bed on freezing times estimating**

** Congelamento de modelos de polpas de frutas com ar forçado. Influência de canais preferenciais no leito sobre as estimativas dos tempos de congelamento**

**Jaime Vilela de Resende ^{I, *}; Vivaldo Silveira Jr.^{II}; Lincoln de C. Neves Filho^{II}**

^{I}LEST/FEMEC/UFU Av. João Naves de Àvila, 2160 Bloco 1M Campus Santa Mônica CEP 38400-089. Uberlândia/MG. E-mail: jvresende@mecanica.ufu.br

^{II}DEA/FEA/ Universidade Estadual de Campinas Caixa Postal 6121 Cid. Universitária Zeferino Vaz Distr. Barão Geraldo. CEP 13083-970 Campinas/SP Brazil. E-mails: vivaldo@cers.fea.unicamp.br – neveslin@ceres.fea.unicamp.br

**SUMMARY**

The freezing times of fruit pulp models packed and conditioned in multi-layered boxes were evaluated under conditions similar to those employed commercially. Estimating the freezing time is a difficult practice due to the presence of significant voids in the boxes, whose influence may be analyzed by means of various methods. In this study, a procedure for estimating freezing time by using the models described in the literature was compared with experimental measurements by collecting time/temperature data. The following results show that the airflow through packages is a significant parameter for freezing time estimation. When the presence of preferential channels was considered, the predicted freezing time in the models could be 10% lower than the experimental values, depending on the method. The isotherms traced as a function of the location of the samples inside the boxes showed the displacement of the thermal center in relation to the geometric center of the product.

**Keywords: **freezing times; frozen foods; heat transfer.

**RESUMO**

Os tempos de congelamento da solução modelo de polpas de frutas embaladas e acondicionadas em diferentes arranjos de camadas no interior das caixas do empilhamento foram avaliados sob várias condições de operação encontradas nas práticas comerciais. Nestes tipos de processos, as caixas contêm uma porosidade significativa e métodos para avaliar como esta afeta o tempo de congelamento, torna o problema de difícil resolução. Procedimentos para as estimativas dos tempos de congelamento, usando modelos encontrados na literatura, foram comparados com medidas experimentais através da aquisição de dados de tempo/temperatura. Os resultados mostraram que o escoamento do ar entre as embalagens é um parâmetro significativo nas estimativas dos tempos de congelamento. Quando é considerada a presença de canais preferenciais, dependendo do método utilizado, estes valores são subestimados em cerca de 10% quando comparados com os valores experimentais. As isotermas traçadas em função da localização das amostras no interior das caixas comprovam o deslocamento do centro térmico, em relação ao centro geométrico do produto.

**Palavras-chave: **tempo de congelamento; alimentos congelados; calor-transmissão.

**1 – INTRODUCTION**

In Brazil, fruit pulp falls into the important category of frozen foodstuffs that are used *in natura* and as raw material for ice cream, yogurt and jelly production. As fruit pulp is presented in diverse package types and sizes, a great increase in its production and consumption has been recorded over the last decade.

When fruits are frozen in their pulp form, the effect of the freezing process is verified by the alteration in the consistency of the pulp and through changes in its composition caused by chemical reactions during the subsequent storage. The consistency of the pulp and global appearance is better kept when the fruit pulp is frozen quickly.

In industrial situations, the freezing process of fruit pulp inside boxes with large amounts of product may be very complex. The initial temperature of the product before being submitted to the freezing process may be uniform or not. The homogeneity of the product is not guaranteed and there is also the possibility of the existence of void, air bags or bubbles filled with gas inside the product and the packages [2] . The product may have form and regular size, but this does not occur frequently. The external heat transfer can be constant through the entire product surface, but wide variations in heat transfer coefficients, mainly the convective, may occur in different positions [4].

The final quality characteristics of the pulp are reflected by the efficiency of the freezing process, and the time during which the product remains in the process is one of the main factors for its preservation. Industrial processing, in most cases, uses the same operational conditions for different types of products, and the freezing time is pre-established in accordance with the amount of product to be processed. The process is performed in air blast freezers of empirical form and does not posses systematic control of operational parameters which can directly contribute to the reduction of freezing time and can influence the economic aspects of the process and of quality and safety of the product.

Considering the complexity of this industrial process, the objective of this study is to evaluate the effect of the irregularities of the bed (preferential channels) and the amount and product configuration study on the prediction of freezing times of fruits pulp packed in polyethylene bags (100 grams) and conditioned in commercial plastic boxes in air blast freezing process under different operational conditions of the equipment.

**2 – METHODS AND MATERIALS**

**2.1 – Monitoring of the boxes containing the samples**

A system model consisting of K-Carrageenan 0.5% and sucrose 10% (mass/volume of water) was used to simulate the fruit pulp. To obtain the description of temperature records during the freezing process, 51 thermocouples (Type T, cooper-constantan) were installed through the entire experimental configuration; with 17 having been installed in each of the three boxes of the pile up in points chosen between the layers of the array of samples such that the transient behavior of the heat transfer during the process could be better translated (*Figure 1-A-B*).

Each thermocouple was installed in the geometric center of the packages (100g). The packages were filled with the product, and after that sealed and placed inside the boxes. In *Figure 1-B*, the distribution of the samples containing thermocouples inside the boxes is shown for arrays of 3 layers of the product. The position of each thermocouple in the cartesian coordinates system for arrays of 3 layers is presented in *Figure 1-C*.

**2.2 – Initial stabilization of the temperature of the product**

To provide the initial stabilization of the temperature of the product and the thermal homogeneity of the samples, the entire monitored system, consisting of three boxes containing thermocouples, was placed in a stabilization chamber with an internal room temperature of 25 and 5 ± 0.5^{o}C for a period of 48 hours. These temperatures are commonly used during the industrial process, when the fruit is frozen immediately after it arrived to the facilities or when it is first cooled and frozen later.

**2.3 – Freezing of the product**

After obtaining the stabilization of the equipment at the operating temperature, the three boxes containing the samples (3, 5 and 7 layers) with thermocouples, were piled up in the freezing chamber. To prevent the formation of preferential channels, expanded polystyrene boards were placed in the lateral sections of the boxes making it so that the air stream of the entire chamber flowed through the openings of the boxes containing the product. The monitoring of the temperature was performed at 1 minute intervals. The temperature of the air cooling was established in 3 levels, -30, -25 and -20^{o}C referring to levels 1, 2 and 3 respectively.

**2.4 – Thermal properties of the samples**

Data of the thermal properties, specific for food models, for the calculations of heat transfer problems were obtained by means of experimental measurements and compared to structural models, prediction equations and adjusted in terms of polynomial functions, divided in a range of temperatures covering the entire experimental scene. These data are available in RESENDE [10].

**2.5 – Freezing time prediction for samples packed in polyethylene bags and conditioned in boxes.**

For the freezing time prediction of the samples packed in polyethylene bags and conditioned in boxes, the use of the geometric factor was considered for rectangular bricks representing the product inside the boxes (equation 1).

**2.5.1 – Freezing time of an infinite slab **

One of the most recommended equations [2, 5] for predicting the freezing time of one-dimensional foodstuffs is presented in the method of Pham [8] as a basic equation, for which the factor number (E) can be applied. The freezing is calculated by the equation (2). Equation (2) considers the entire freezing process (pre-cooling, phase change and tempering).

where,

A | heat transfer surface area (m^{2}); |

Bi_{s} | Biot number (h_{s}/RK); |

h_{s} | surface heat transfer coefficient (W/m^{2} k); |

DH_{i} | change in relative enthalpy content (J/m^{3}), i = 1 (pre-cooling) and i = 2 (freezing and tempering); |

DT_{i} | change in temperature (K), i = 1 (pre-cooling) and i = 2 (freezing and tempering); |

V | volume (m^{3}). |

The values for the effective heat transfer coefficients used in the calculations were obtained by the average of the values from the Nusselt *versus* Reynolds correlations for the different locations (upper and lower layers) in the boxes represented by the equations:

For the superior layer: Nu = 0,01575Re^{0,81439 }(3)

For the inferior layer: Nu = 0,00754Re^{0,85461 }(4)

When L'/D_{H} __<__ 60, h = h_{2(eq)}[1+1,4(D_{H}/L'), where h_{2(eq)} is the value from the equation (3) and (4) and L' is the length of the duct.

These correlations were developed from the temperatures records of an aluminum transducer with the same dimensions of the samples, in the experimental determination of the airflow rate across the cross section of air blast freezing, and in the concept of the hydraulic diameter, that is established by the amount and number of layers of product inside the boxes [10].

**2.5.2 – Geometric factor calculation (E)**

The geometric factor (E) was calculated using two methods: i) method EHTD ("Equivalent Heat Transfer Dimensions") considered by Cleland, Earle [3] given by the equation:

E = 1 + W_{1} + W_{2 }(5)

The functions W1 and W2 depend on the Biot number (Bi_{s}) and on the dimensional ratios b_{1} and b_{2} in accordance with the equations:

and, ii) the geometric factor (E_{AN}) derived from the analytical solution of heat transfer problems with phase change [7].

where the values of Z_{n} e Z_{m} are the roots of equations:

Bi_{s} = Z_{n}tanZ_{n }(9)

and Bi_{s}´b_{1} = Z_{m}tanZ_{m }(10)

and the values of Z_{nm} from the equation (11):

The six first roots of equations (9) and (10) are tabled in Carslaw & Jaeger [1] as a function of the values of Bi_{s} and Bi_{s}´b_{1}.

When trying to find a solution, it was verified that in many configurations tested values for the dimensional ratios b_{1} e b_{2} > 3 were found. In accordance with Hossain [6], when this occurs or when the Biot number is under 3 (Bi_{s} < 3), most of the methods become limited and inexact, and, in these cases, the full analytical method must be used for the best prediction of the freezing time.

**2.6 – Geometric configurations **

Two types of configuration of the rectangular bricks that could describe the real environment of the piling up of the boxes containing the samples of the product in air blast freezing were analyzed in accordance with the following proposals (*Figure 2*): (A) the product inside the boxes was considered as a compact rectangular brick; (B) the product was considered as a group of five identical rectangular bricks evenly spaced inside the boxes with the largest dimensions parallel to the direction of the airflow.

The geometric characteristics of each configuration and the values of the parameters (b_{1}, b_{2}, V, A, V/A), used in these equations for the freezing prediction, are shown in *Table (1).*

**2.7 – Comparison between experimental and calculated results for the freezing times**

The comparison was made through the relation of the difference between predicted and experimental values in accordance with equation:

The standard deviation (SD) was calculated by the equation:

where

where is the average error and n is the size of the sample.

**3 – RESULTS AND DISCUSSION**

**3.1 – Curves of freezing **

The graph of *Figure 3*, shows the descriptions of the temperatures records during the freezing process by the reading of the 17 thermocouples located in box 2 (central), for an experiment with arrays of 3 layers, level 2 of the temperature of the air cooling and initial temperature of the product around 5^{o}C. It was found that the freezing of the samples located in the front portion of the boxes (coordinates X = 0,0361 m) referring to thermocouples that suffer the direct impact of the airflow originated from the evaporator, presents a faster fall temperature rate when compared to those located in the back portion of the boxes (coordinate X = 0,3249m). It can also be observed that there are notable differences among the temperature records for samples placed between the lower and the upper layers of the boxes and for samples placed in the central points of the different layers.

**3.2 – Thermal center location**

*Figure 4* shows the isotherms obtained from the non linear regression built using the temperatures readings from the 17 thermocouples (*Figure 1*) installed in different positions within box 2 at intervals of 30 minutes during the freezing process. *Figure 4* shows this advance on (XY) coordinates system which encloses the width and the length of the box and the movement of the front freezing occurs in a parallel direction to the airflow, and in *Figure 5* this advance is seen on the (XZ) plane that represents the perpendicular section to the airflow correspondent to the height and the length of the box. The points of readings are represented by the points shown on the grid and the isotherms are obtained from the averages of these temperatures in function of their locations in the plane. In each assay, 51 temperature records were obtained supplying data for all the piling up.

Previous work [9, 10] found that there was a significant difference between the effective heat transfer coefficients for the samples placed in different positions (lower and upper) of the layer arrays inside the boxes. These comments suggest that in the different processing conditions it is not possible to obtain a uniform heat transfer in all surfaces of the product and the replacement of the thermal center in the configurations certainly will not occur coinciding with the geometric center of the considered forms.

*Figures 4 *and* 5* prove the replacement of the thermal center compared to the geometric center of the product during the freezing and show that the samples placed in the second half of the (XY) plane of the configurations are the last ones to be frozen and the front freezing is sped up in the direction of the airflow.

It is observed that the replacement of the front freezing occurs in all surfaces, and the velocity of advance in direction to the thermal center is more noticeable from the front surface and on the upper layers of the boxes containing the samples. Therefore in these cases, the shape of the product is such that the heat transfer occurs in more than one direction and the calculations of the freezing time, must be dealt using models for multidimensional forms, as considered previously.

The thermal center is defined as the last point of the foodstuffs to suffer changes in temperature from the freezing process. In all the experiments the final temperature of thermal center (T_{c}) was of -18^{o}C and the experimental freezing time was taken when all the thermocouples placed in all the configuration had shown lower temperatures immediately after this.

**3.3 – Evaluation of the freezing times for samples conditioned inside the boxes**

*Table 2* shows the experimental freezing times for all the configuration proposals when the initial temperatures (T_{i}) of the product was around 25^{o}C and 5^{o}C after stabilization. In *Table 3* the results of the methods for the calculation of the geometric factor (E) are presented using the two methods and for the two configuration proposals.

The different methods for the calculation of the geometric factor have been tested by its authors mainly using data obtained from simulation tests. In accordance with Salvadori, De Michaelis, Mascheroni [11], the correct way to verify the effectiveness of a predictive method is to verify the predictions that correspond to an available experimental data set. In this study, (*Table 2*) 36 experimental assays were performed, with varied geometric configurations of the product, initial temperatures and operational conditions of the equipment. *Tables 4* and *5* show the average errors and standard deviation of the freezing times using the methods tested for these data. The results are divided by the geometric factor method of calculation, by the level of air temperature and by the configuration type considered. *Figures 6* and *7* allow a better visualization for comparison of the results, showing the data of *Tables 4* and *5* through bar diagrams.

It is observed that for the configuration that considers the product as a compact rectangular brick, the freezing times predicted by the models are overestimated, in some cases, with upper average errors of 100% when compared with the experimental freezing times, independently of the method of calculation and of the geometric factor. This configuration does not consider the empty spaces between the samples, and, the superficial area of heat transfer applied to the calculations is reduced not describing the real situation of the passage of air cooled through the product inside the freezing chamber. For the other configuration the existence of preferential channels between the samples is considered and the average errors produced by the models, are sufficiently reduced.

The best results are observed when the configuration applied to the model is an array consisting of 5 rectangular bricks of the product evenly spaced with its largest dimensions parallel to the airflow direction. In this configuration, for arrays of 7 layers of samples, the predicted freezing times are underestimated, however the average errors when compared to the experimental results are reduced. For procedures with arrays of 5 and 3 layers of samples, this type of configuration applied to the models yielded good results with average errors of prediction 12% lower when compared to the experimental values. This model considers the passage of air through preferential channels parallel to the air flow direction formed between a group of samples and an increase in the surface of heat transfer because it includes in the lateral area of the individual samples inside the bed into the calculations.

The effectiveness of the prediction method can be noticed when the procedure applied to the prediction models was the one constituted by arrays of various bricks inside the boxes (configuration B) and the geometric factor (E_{AN}) was calculated using the equation (8) of Hossain, Cleland, Cleland [7] developed by the analytical solution of heat transfer problems in the transient state.

The method presented good values for the freezing time prediction when compared to the experimental results in all levels of air temperature, with results more compatible for experiments with arrays of 5 and 3 layers of samples inside of the boxes. For the freezing time prediction for boxes containing arrays of 7 layers, a great difference between the predicted and experimental values was observed, but when compared to the other applied conditions, it presented the lesser average prediction errors. The presence of these channels inside the boxes between the samples increase the heat transfer rate by increasing the air velocity, and, is probably, also influenced by edge effect and sudden changes of the air direction in these spaces. The first type of tested configuration (compact rectangular bricks) does not consider the existence of these channels.

**4 – CONCLUSIONS**

A method of prediction developed for regular forms must be used with certain care in situations where apparently the conformation of the product is similar. It is necessary to consider the presence of preferential channels, voids and irregularities that can exist in the product or in the shape during processing. In this study, it was found that differences in the behavior of the cooling rate exist between the boxes placed in different positions of the pile, and that the heat transfer occurs in all the directions and is not uniform in all the surfaces. When the details of the flow are not taken into account, average errors in the estimation of freezing times has been above 100% in some cases. On the other hand, the configuration that considers the existence of preferential channels in the bed yielded good results in all the tested levels of temperature, resulting in average errors of the predicted values below 10%. The heat transfer coefficients greatly affect the accuracy of the prediction methods, and in a previous work [10] , specific correlations from measures in this experimental procedure were developed. Therefore, a prediction method needs: to be tested and compared with experimental data, to have quality in the measurements, and, an analysis must be made to represent the configuration of the experimental environment to get compatible estimates with the real process applicable in the calculations of engineering.

**5 – ACKNOWLEDGEMENTS**

Acknowledgements are made to CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) for financially supporting this study.

**6 – REFERENCES**

[1] Carslaw, H. S.; Jaeger, J. C. Conduction of heat in solids. 2th Edition, Oxford University Press, 1959. [ Links ]

[2] Cleland, A. C. Food refrigeration process. Analysis, design and simulation. Elsevier Applied Science: London e New York. 284p., 1990. [ Links ]

[3] Cleland, A. C.; Earle; R. L. Freezing time prediction for foods – a simplified procedure. **International Journal of Refrigeration**, v.5, n.3, p.134-140, 1982. [ Links ]

[4] Cleland, A.C.; Ôzilgen, S. Thermal design calculations for food freezing equipment- past, present and future. **International Journal of Food Refrigeration**, v.21, n.5, p. 359-371, 1998. [ Links ]

[5] Delgado, A. E.; Da-Wen Sun Heat and mass transfer models for predicting freezing process – a review. **Journal of Food Engineering**, v.47, p.157-174, 2001. [ Links ]

[6] Hossain, M. M. A simple method of freezing time calculation for foodstuffs of various shapes. **Food Austrália**, v.47, n.3, p.109-112, 1995. [ Links ]

[7] Hossain, M. M., Cleland, D. J.; Cleland, A. C. Prediction of freezing and thawing times for foods of regular multi-dimensional shape by using an analytically derived geometric factor*. ***International Journal of Refrigeration**, v.15, n.4, p.227-234, 1992. [ Links ]

[8] Pham, Q. T. Simplified equation for predicting the freezing time of foodstuffs. **Journal of Food Technology**, v.21, p.209-219, 1986. [ Links ]

[9] Pham, Q.T.; Willix, J. Heat transfer coefficients in the air blast freezing of rows of cartons. **Proceedings 17t ^{h }Int. Congress of Refrigeration.**, pp.350-357,1987. [ Links ]

[10] Resende, J. V. Heat transfer analyses for the freezing of fruit pulps in commercial boxes. Campinas, 2001. (D.Sc. thesis in Food Engineering), Faculty of Food Engineering, State University of Campinas (UNICAMP). [ Links ]

[11] Salvadori, V. O., De Michaelis; Mascheroni, R. H. Prediction of freezing times for regular multi-dimensional foods using simple formulae**. Lebensmittel Wissenchaft und Technology***.* v.30, p.30-35, 1997. [ Links ]

* Recebido para publicação em 11/10/2001. Aceito para publicação em 08/03/2002. *

** A quem a correspondência deve ser enviada.*