×

Estimating the temperature evolution of foodstuffs during freezing with a 3D meshless numerical method. (English) Zbl 1403.80039

Summary: Freezing processes are characterised by sharp changes in specific heat capacity and thermal conductivity for temperatures close to the freezing point. This leads to strong nonlinearities in the governing PDE that may be difficult to resolve using traditional numerical methods. In this work we present a meshless numerical method, based on a local Hermite radial basis function collocation approach in finite differencing mode, to allow the solution of freezing problems. By introducing a Kirchhoff transformation and solving the governing equations in Kirchhoff space, the strength of nonlinearity is reduced while preserving the structure of the heat equation. In combination with the high-resolution meshless numerical method, this allows efficient and stable numerical solutions to be obtained for freezing problems using 3D unstructured datasets. We demonstrate the proposed numerical formulation for the freezing of foodstuffs. By approximating the shape of the thermal conductivity and heat capacity curves in a piecewise linear fashion the temperature-dependent material curves may be described using eight independent parameters. We consider the optimisation of these parameters to match experimental data for the freezing of a hemispherical sample of mashed potato by using a simple manual procedure that requires a minimal number of simulations to be performed. Working in this way, a good approximation is obtained to the temperature profile throughout the sample without introducing instability into the numerical results.

MSC:

80M22 Spectral, collocation and related (meshless) methods applied to problems in thermodynamics and heat transfer
65M70 Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs
80A22 Stefan problems, phase changes, etc.
Full Text: DOI

References:

[1] Becker, B. R.; Fricke, B. A., Hydrocooling time estimation, Int Commun Heat Mass Transf, 29, 165-174, (2002)
[2] Hu, Z.; Sun, D. W., Cfd simulation of heat and moisture transfer for predicting cooling rate and weight loss of cooked ham during air-blast chilling process, J Food Eng, 46, 189-197, (2000)
[3] Puri, V. M.; Anantheswaran, R. C., The finite-element method in food processing—a review, J Food Eng, 19, 247-274, (1993)
[4] Pham, Tuan, Modelling heat and mass transfer in frozen foodsa review, Int J Refrig, 29, 876-888, (2006)
[5] Lin, S. H., Transient heat conduction with temperature-dependent thermal conductivity by the orthogonal collocation method, Lett Heat Mass Transf, 5, 29-39, (1978)
[6] LaRocca, A.; Hernandez Rosales, A.; Power, H., Radial basis function Hermite collocation approach for the solution of time dependent convection-diffusion problems, Eng Anal Bound Elem, 29, 4, 359-370, (2005) · Zbl 1182.76928
[7] LaRocca, A.; Power, H.; LaRocca, V.; Morale, M., A meshless approach based upon radial basis function Hermite collocation method for predicting the cooling and the freezing times of foods, Comput Mater Contin, 2, 239-250, (2005)
[8] Driscoll, T. A.; Fornberg, B., Interpolation in the limit of increasingly flat radial basis functions, Comput Math Appl, 43, 413-422, (2002) · Zbl 1006.65013
[9] Robert Schaback. Multivariate interpolation and approximation by translates of a basis function. In: Approximation theory VIII, 1995. · Zbl 1139.41301
[10] Brown, D., On approximate cardinal preconditioning methods for solving PDEs with radial basis functions, Eng Anal Bound Elem, 29, 343-353, (2005) · Zbl 1182.65174
[11] Beatson, R. K.; Cherrie, J. B.; Mouat, C. T., Fast Fitting of radial basis functionsmethods based on preconditioned GMRES iteration, Adv Comput Math, 11, 253-270, (1999) · Zbl 0940.65011
[12] Ling, L.; Kansa, E. J., A least-squares preconditioner for radial basis functions collocation methods, Adv Comput Math, 23, 31-54, (2005) · Zbl 1067.65136
[13] Baxter, B. J.C., Preconditioned conjugate gradients, radial basis functions, and Toeplitz matrices, Comput Math Appl, 43, 305-318, (2002) · Zbl 1002.65018
[14] Ling L, Schaback R. On adaptive unsymmetric meshless collocation. In: Atluri SN, Tadeu AJB, editors. Proceedings of the 2004 international conference on computational and experimental engineering and sciences. Madeira, Portugal: Tech. Science Press; 2004. · Zbl 1167.65059
[15] Ling, L.; Opfer, R.; Schaback, R., Results on meshless collocation techniques, Eng Anal Bound Elem, 30, 247-253, (2006) · Zbl 1195.65177
[16] Hernandez-Rosales, A.; Power, H., Non-overlapping domain decomposition algorithm for the Hermite radial basis function meshless collocation approachapplications to convection diffusion problems, J Alg Technol, 33, 127-159, (2007)
[17] Zhou, X.; Hon, Y. C.; Li, J., Overlapping domain decomposition method by radial basis functions, Appl Numer Math, 44, 241-255, (2003) · Zbl 1013.65134
[18] Zhang, Yunxin, Reconstruct multiscale functions using different RBFs in different subdomains, Appl Math Comput, 189, 893-901, (2007) · Zbl 1127.65010
[19] Wong, A. S.M.; Hon, Y. C.; Li, T. S.; Chung, S. L.; Kansa, E. J., Multizone decomposition for simulation of time-dependent problems using the multiquadric scheme, Comput Math Appl, 37, 23-43, (1999) · Zbl 0951.76066
[20] Wu, Z., Compactly supported positive definite radial basis functions, Adv Comput Math, 4, 75-97, (1995)
[21] Wendland, H., Piecewise polynomial, positive definite and compactly supported radial basis functions of minimal degree, Adv Comput Math, 4, 389-396, (1995) · Zbl 0838.41014
[22] Schaback, R., On the efficiency of interpolation by radial basis functions, (LeMehaute, A.; Rabut, C.; Schumaker, L. L., Surface fitting and multiresolution methods, (1997), Vanderbilt University Press Nashville, TN), 309-318 · Zbl 0937.65013
[23] Fasshauer, G. E., Solving differential equations with radial basis functionsmultilevel methods and smoothing, Adv Comput Math, 11, 139-159, (1999) · Zbl 0940.65122
[24] Lee, C.; Liu, X.; Fan, S., Local multiquadric approximation for solving boundary value problems, Comput Mech, 30, 396-409, (2003) · Zbl 1035.65136
[25] Sarler, B.; Vertnik, R., Meshless explicit local radial basis function collocation methods for diffusion problems, Comput Math Appl, 51, 1269-1282, (2006) · Zbl 1168.41003
[26] Divo, E.; Kassab, A., An efficient localised radial basis function meshless method for fluid flow and conjugate heat transfer, J Heat Transf, 129, 124-136, (2007)
[27] Wright, G.; Fornberg, B., Scattered node compact finite difference-type formulas generated from radial basis functions, J Comput Phys, 212, 99-123, (2006) · Zbl 1089.65020
[28] Stevens, D.; Power, H.; Lees, M.; Morvan, H., The use of PDE centres in the local RBF Hermitian method for 3D convective-diffusion problems, J Comput Phys, 228, 4606-4624, (2009) · Zbl 1167.65447
[29] Stevens, D.; Power, H., A scalable meshless formulation based on RBF Hermitian interpolation for 3d nonlinear heat conduction problems, Comput Model Eng Sci, 55, 111-146, (2010) · Zbl 1231.65182
[30] Rohsenow, W.; Hartnett, J.; Ganic, E., Handbook of heat transfer fundamentals, (1985), McGraw-Hill Book Company New York
[31] Carlslaw, H.; Jaeger, J., Conduction of heat in solids, (1959), Oxford at the Clarendon Press Oxford
[32] Bagnall, K. R.; Muzychka, Y. S.; Wang, E. N., Application of the Kirchhoff transform to thermal spreading problems with convection boundary conditions, IEEE Trans Compon Packag Manuf Technol, 4, 408-420, (2014)
[33] Chantasiriwan, S., Multiquadric collocation method for time-dependent heat conduction problems with temperature-dependent thermal properties, J Heat Transf, 129, 109-113, (2007)
[34] Kim, S., A simple direct estimation of temperature-dependent thermal conductivity with Kirchhoff transformation, Int Commun Heat Mass Transf, 28, 537-544, (2001)
[35] Bonacina, C. G.; Comini, G.; Fasano, A.; Primicerio, M., On the estimation of thermophysical properties in non linear heat-conduction problems, Int J Heat Mass Transf, 17, 861-867, (1974)
[36] Wu, Z., Hermite-Birkhoff interpolation of scattered data by radial basis functions, Approx Theory Appl, 8, 1-11, (1992) · Zbl 0757.41009
[37] Schaback, R., Convergence of unsymmetric kernel-based meshless collocation methods, SIAM J Numer Anal, 45, 333-351, (2007) · Zbl 1165.65066
[38] Hon, Y. C.; Schaback, R., On unsymmetric collocation by radial basis functions, J Appl Math Comput, 119, 177-186, (2001) · Zbl 1026.65107
[39] Power, H.; Barraco, V., A comparison analysis between unsymmetric and symmetric radial basis function collocation methods for the numerical solution of partial differential equations, Comput Math Appl, 43, 551-583, (2002) · Zbl 0999.65135
[40] LaRocca, A.; Power, H., A double boundary collocation Hermitian approach for the solution of steady state convection diffusion problems, Comput Math Appl, 55, 1950-1960, (2007) · Zbl 1139.65077
[41] Madych, W.; Nelson, S., Multivariate interpolation and conditionally positive definite functions ii, Math Comput, 54, 211-230, (1990) · Zbl 0859.41004
[42] Stevens, D.; Power, H., A scalable and implicit meshless RBF method for the 3D unsteady nonlinear Richards equation, with single and multi-zone domains, Int J Numer Methods Eng, 85, 135-163, (2011) · Zbl 1217.76056
[43] Suzuki, S.; Kashiwagi, J.; Nakagawa, S.; Souma, K., Mechanism of hysteresis in thermal conductivity of frozen solid between freezing and thawing processes, Trans Jpn Soc Irrig Drain Reclam Eng, 218, 223-231, (2002)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.