Microwave Modeling and Validation in Food Thawing Applications

Developing temperature fields in frozen cheese sauce undergoing microwave heating were simulated and measured. Two scenarios were investigated: a centric and offset placement on the rotating turntable. Numerical modeling was performed using a dedicated electromagnetic Finite Difference Time Domain (FDTD) module that was two-way coupled to the PHYSICA multiphysics package. Two meshes were used: the food material and container were meshed for the heat transfer and the microwave oven cavity and waveguide were meshed for the microwave field. Power densities obtained on the structured FDTD mesh were mapped onto the unstructured finite volume method mesh for each time-step/turntable position. On heating for each specified time-step the temperature field was mapped back onto the FDTD mesh and the electromagnetic properties were updated accordingly. Changes in thermal/electric properties associated with the phase transition were fully accounted for as well as heat losses from product to cavity. Detailed comparisons were carried out for the centric and offset placements, comparing experimental temperature profiles during microwave thawing with those obtained by numerical simulation.

[1]  G. Mur The Modeling of Singularities in the Finite-Difference Approximation of the Time-Domain Electromagnetic-Field Equations , 1981 .

[2]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[3]  G. Mur Absorbing Boundary Conditions for the Finite-Difference Approximation of the Time-Domain Electromagnetic-Field Equations , 1981, IEEE Transactions on Electromagnetic Compatibility.

[4]  L. Barratt,et al.  Experimental validation of a combined electromagnetic and thermal FDTD model of a microwave heating process , 1995 .

[5]  A K Datta,et al.  Coupled Electromagnetic and Termal Modeling of Microwave Oven Heating of Foods , 2000, The Journal of microwave power and electromagnetic energy : a publication of the International Microwave Power Institute.

[6]  E. Süli,et al.  A convergence analysis of Yee's scheme on nonuniform grids , 1994 .

[7]  Koulis Pericleous,et al.  A new computational approach to microwave heating of two‐phase porous materials , 2004 .

[8]  Kazuo Aoki,et al.  The characteristics of microwave melting of frozen packed beds using a rectangular waveguide , 2002 .

[9]  M. Cross,et al.  An explicit scheme for coupling temperature and concentration fields in solidification models , 2004 .

[10]  J. Benek,et al.  A flexible grid embedding technique with application to the Euler equations , 1983 .

[11]  P. Risman,et al.  Coupled electromagnetic and thermal modeling of microwave oven heating of foods. , 2000, The Journal of microwave power and electromagnetic energy : a publication of the International Microwave Power Institute.

[12]  A NUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE MICROWAVE HEATING OF POLYMER MATERIALS INSIDE A RIDGE WAVEGUIDE , 1996 .

[13]  K. Ayappa,et al.  Analysis of microwave heating of materials with temperature-dependent properties , 1991 .

[14]  M. Celuch,et al.  Towards a multiphysics simulation system for microwave power phenomena , 2005, 2005 Asia-Pacific Microwave Conference Proceedings.

[15]  A. Datta,et al.  Heating Concentrations of Microwaves in Spherical and Cylindrical Foods: Part One: in Planes Waves , 2005 .

[16]  V. Voller,et al.  The modelling of heat, mass and solute transport in solidification systems , 1989 .

[17]  K. Yee Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[18]  A. C. Metaxas,et al.  Industrial Microwave Heating , 1988 .

[19]  G. Georghiou,et al.  Numerical Modeling for Rf and Microwave Heating Using Fe Comes of Age , 2001, The Journal of microwave power and electromagnetic energy : a publication of the International Microwave Power Institute.

[20]  Joseph O'Rourke,et al.  Computational Geometry in C. , 1995 .

[21]  N. Sakai,et al.  Power and temperature distribution during microwave thawing, simulated by using Maxwell's equations and Lambert's law , 2005 .

[22]  V. Voller,et al.  A fixed grid numerical modelling methodology for convection-diffusion mushy region phase-change problems , 1987 .

[23]  J. O´Rourke,et al.  Computational Geometry in C: Arrangements , 1998 .

[24]  Vaughan R Voller,et al.  Modelling the mushy region in a binary alloy , 1990 .

[25]  C. Balanis Advanced Engineering Electromagnetics , 1989 .

[26]  P. Kopyt,et al.  FDTD modeling and experimental verification of electromagnetic power dissipated in domestic microwave ovens , 2003 .

[27]  Vaughan R Voller,et al.  Towards a general numerical scheme for solidification systems , 1997 .

[28]  Ashim K. Datta,et al.  Heating Concentrations of Microwaves in Spherical and Cylindrical Foods: Part Two: in a Cavity , 2005 .

[29]  Peter Monk,et al.  Phase-accuracy comparisons and improved far-field estimates for 3-D edge elements on tetrahedral meshes , 2001 .