Comparison of efficiency between two different numerical modeling methods to predict tomato paste temperature during pasteurization process.

Document Type : Research Paper

Author

Department of Food Science and Technology, Faculty of Agriculture, Jahrom University, Jahrom, Fars, Iran, P.O. Box 74137-66171

Abstract

Introduction: Tomato (Solanum lycopersicum) is the second most important horticultural crops next to potato, with an estimated world production of over 120 million tons per year. Iran has always been among the top tomato producing countries in the world due to the diverse topography and climactic conditions prevailing in different parts of Iran. Tomato paste is one of the most important tomato products originates from tomato juice in which water is removed by evaporation and thermal processing plays a key role in the successful production process (Singh and Headman, 2014). Since it can guarantee safety and extend shelf life of the product. Although several emerging technologies such as ohmic heating, microwave heating, and non-thermal processing techniques such as pulsed electric fields and high pressure processing, have been developed for food preservation, Conventional thermal processing (pasteurization and sterilization) is still widely used in the Iranian food industry. However, the challenges of accurately determining both optimal operating conditions and developing a control system for industrial pasteurization process to prevent either under- or overprocessing are significant (Chen and Ramaswamy, 2007). Mathematical modeling and simulation is one of the most commonly used methods to gain a better understanding the process. Modeling can provide insights into complex processes, shorten the design cycle and optimize the process as a function of various variables at lower cost and time. Different mathematical methods for solving heat conduction problems have been proposed, but numerical methods are more useful, especially when such problems cannot be handled by the exact analysis because of nonlinearities, complex geometries, and complicated boundary conditions (Incropera and De Witt, 1990). Among the numerical methods developed so far, finite difference and finite element techniques have been widely used to analyze heat transfer phenomena in cylindrical cans of food. Although many researchers have tried to develop mathematical models based on numerical methods for predicting temperatures in tomato products during thermal processing, such as (Bichier et al., 1995), (Nicolai et al., 1998), (Tattiyakul et al., 2002) and (Plazl et al., 2006), there is not much research available on the subject of thermal processing of tomato paste in Iran.
In this study, the efficiency of two different modeling approaches (finite difference vs. finite element) for predicting temperature of tomato paste during pasteurization process were compared. The developed model allows to reduce energy consumption during operations while high-quality products are also produced in a short time.
Material and method: Experiments were run with batches of 400 g of tomato paste (pH = 4.1 and 28°Brix) in cylindrical cans (211×400) and hot water was used as the heating medium. The chemical analysis of tomato paste sample (Brix, pH salt, moisture, fat and protein) was performed in the first step and the thermal properties of the tomato paste product including thermal conductivity, specific heat and density were determined based on the sample chemical composition and structures model.
Temperature changes at various positions in the container were checked with a data logger (Testo, Germany) coupled with computer and thermocouples type-K (at 2 min intervals). The cylindrical can was immersed in a vertical position in the water bath and the temperature recording was started. After finishing the heating time, the can was cooled in another water bath (20ºC). The data were used to validate the developed model.
In the next step, 2D heat transfer model was developed in a cylindrical can by using the numerical solution of the Fourier second law with two different methods. 1) Finite difference (explicit scheme) and finite element methods. Computer simulation is done using MATLAB R2009a software (Math works, Inc., Natick, MA, USA) and COMSOL Multiphysic, Ver. 4.0. Finally, in order to evaluate the best model, two criteria, coefficient of determination (R2) and root mean squared error (RMSE) were used.
Result and discussion: The results showed that, by placing the sample in the bath, the surface temperature rises rapidly, while the temperature in the center is much slower. In addition, as can be expected, increasing hot water temperature enhanced the heating rate considerably due to the larger temperature gradient between the center and surface of the can at the higher temperatures. The models have been verified by comparing results with two analytical solutions and validated against experimental data. The statistical analysis results showed that the finite element model developed by COMSOL software can predicted temperature more accurate than finite difference model and may be more useful. The reason for this difference between the results of two numerical methods can be attributed to the consideration of a layer of air-steam mixture on the top of the can (head space) in the finite element method which increases the accuracy of the model in temperature prediction. After validation, the developed model was used to determine the cold spot location of the tomato paste can. Results also showed that the cold point was a stationary point and located at the radial center at a height of 60% of the can height from the bottom (Tattiyakul et al., 2002). Two simulations were conducted at two different head space volume (6 and 12% of total can height) to determine the importance of head space volume on cold point location. Results showed that there was no significant difference in the location and temperature of the cold spot in two simulations (Khakbaz Heshmati et al., 2014).
Conclusion: In this study, the pasteurization process of tomato paste (Brix=28) is investigated by two different numerical methods (finite difference& finite element). The results were compared with experimental data and it was found that the predicted temperature by finite element model is more accurate than finite difference method. Although it is generally believed that the coldest point for a solid product will be at the geometric center of the cylinder, our results indicated that the slowest cooling point was located at a height of 60% of the can height from the bottom. The developed model can predict temperature in tomato paste with different degree of concentration (brix) or different thermal processing conditions and with minor modifications, the model may be used to design and control the process of industrial pasteurization for various solid products. In addition, the results of this study is expected to be a significant contribution for further optimization studies.

Keywords


استاندارد ملی ایران 761 ، کنسرو رب گوجه فرنگی- ویژگی ها و روش های آزمون، موسسه استاندارد و تحقیقات صنعتی ایران.
بی نام، 1398. آمار نامه صادرات کالا. اداره کل آمار و اطلاعات سازمان توسعه تجارت ایران.
حسینی، ز، 1369. روشهای متداول در تجزیه مواد غذایی. انتشارات دانشگاه شیراز،150-200.
دلوی، م و همدمی، ن. 1389. مدل سازی عددی انتقال حرارت در پنیر سفید فراپالایش شده. نشریه پژوهش های صنایع غذایی، 3(2)، 45-60.
شهابی قهفرخی، ا و بهرامی، ب. مدل سازی انتقال حرارت در طی فرایند کنسرو رب گوجه فرنگی در داخل اتوکلاو. دومین همایش و نمایشگاه بزرگ صنایع غذایی ایران. دانشگاه صنعتی اصفهان. اصفهان،90-96.
Barreiro, J.A, Milano, M., Sandoval, A.J, 1997. Kinetics of colour change of double concentrated tomato paste during thermal treatment. Journal of Food Engineering 33(3), 359-371.
Bichier, J.G, Teixeira, A.A, Balaban, M.O, Heyliger, T.L, 1995. Thermal process simulation of canned foods under mechanical agitation1. Journal of Food Process Engineering 18(1), 17-40.
Chen, C.R, Ramaswamy, H.S, 2007. Visual Basics computer simulation package for thermal process calculations. Chemical Engineering and Processing: Process Intensification 46(7), 603-613.
Dalvi, M, Hamdami, N, 2011. Characterization of thermophysical properties of Iranian ultrafiltrated white cheese: measurement and modeling. Journal of agricultural science and technology (JAST) 13(1), 67-78.
Drusas, AE, Saravacos, G.D, 1985. Thermal conductivity of tomato paste. Journal of Food Engineering 4(3), 157-168.
Holdsworth, S.D, Simpson, R, 2016. Thermal processing of packaged foods. Springer International Publishing, New York.
Incropera, FP, De Witt, D.P, 1990. Fundamentals of Heat and Mass Transfer. John Wiley and Sons Inc., , New York.
Khakbaz Heshmati, M, Shahedi, M., Hamdami, N, Hejazi, M.A, Motalebi, A.A, Nasirpour, A, 2014. Mathematical Modeling of Heat Transfer and Sterilizing Value Evaluation during Caviar Pasteurization. Journal of Agricultural Science and Technology 16(4), 827-839.
Nicolai, B, De Baerdemaeker, J, 1992. Simulation of heat transfer in foods with stochastic initial and boundary conditions. Food Bioprod Process, Trans IChemE 70(C), 78–82.
Nicolai, B. De Baerdemaeker, J, 1997. Finite element perturbation analysis of non-linear heat conduction problems with random field parameters. International Journal of Numerical Methods and Heat Fluid Flow 7(5), 525-544.
Pérez-Tejeda, G, Vergara-Balderas, F.T, López-Malo, A, Rojas-Laguna, R., Abraham-Juárez, M.d.R, Sosa-Morales, M.E., 2016. Pasteurization treatments for tomato puree using conventional or microwave processes. Journal of Microwave Power and Electromagnetic Energy 50(1), 35-42.
Plazl, I, Lakner, M, Koloini, T, 2006. Modeling of temperature distributions in canned tomato based dip during industrial pasteurization. Journal of Food Engineering 75(3), 400-406.
Sandoval, A.J, Barreiro, JA, Mendoza, S, 1992. Thermal Resistance of Bacillus coagulans in Double Concentrated Tomato Paste. Journal of Food Science 57(6), 1369-1370.
Singh, R.P, Heldman, D.R., 2014. Chapter 8 - Evaporation, in: Singh, R.P., Heldman, D.R. (Eds.), Introduction to Food Engineering (Fifth Edition). Academic Press, San Diego, pp. 565-592.
Tattiyakul, J, Rao, MA, Datta, A.K, 2002. Heat Transfer to Three Canned Fluids of Different Thermo-Rheological Behaviour Under Intermittent Agitation. Food and Bioproducts Processing 80(1), 20-27.
Toledo, R.T, 2007. Fundamentals of Food Process Engineering. Springer, New York,.
Uyar, R, Erdogdu, F, 2012. Numerical Evaluation of Spherical Geometry Approximation for Heating and Cooling of Irregular Shaped Food Products. Journal of Food Science 77(7), E166-E175.