ABSTRACT

Microwave heating of liquid foods in laminar flow through a circular tube has been modeled. In particular, skim milk as a Newtonian fluid and apple sauce and tomato sauce as non-Newtonian fluids have been considered. The temperature profiles have been obtained solving the motion and energy equations in transient regime and Maxwell’s equations in the frequency domain. Numerical resolution of Finite Element Method has been implemented in Comsol Multiphysics. The generation term due to the microwave heating has been evaluated according both to Lambert’s law and Poynting theorem. Finally, a comparison between the two methods has been made in order to check to what extent the results obtained with the simpler Lambert’s law approximation are comparable with those deriving from the exact solution of Maxwell equations. Dielectric properties are considered to be temperature dependent.

Keywords:

Continuous Microwave Heating, Lambert’s Law, Maxwell Equations, Poynting Theorem, Non-Newtonian Fluids, FEM

1. Introduction

Microwave heating has been utilized since the 1940s [1] in different fields such as polymer and ceramics industries [2] [3] and medicine [4] [5] . However, food processing is the largest consumer of microwave energy, that can be employed for cooking, thawing, tempering, drying, freeze drying, pasteurization, sterilization, baking, heating and re-heating. In microwave heating, electromagnetic field polarizes the molecules of dielectric materials and creates dipole moments that cause these molecules to rotate. The resulting molecular friction causes heat generation in the body. Due to intrinsic heat generation capability, microwave heating can provide prompt rise of temperature within the low thermal conductive products, especially in food items. On the other hand, the presence of hot and cold spots determines an uneven heating. This is primarily caused by the non-uniform distribution of microwave energy in the foodstuff, due to factors such as dielectric loss, penetration depth, thickness, shape and size of the sample.

Microwave heating of solid foods has been largely investigated in the last years, writing the energy equation as a conductive heat transfer with a generation term.

The latter has been modeled by many authors, using two different approaches to evaluate the effects of the microwave distribution: by solving the Maxwell’s equations [2] [6] [7] or by applying the Lambert’s law [8] [9] . The Lambert’s law is a simple power formulation that was believed to simulate temperature profiles for Cartesian geometries and for cylindrical geometries with high radius. Ayappa et al. [2] deduced the minimum value of the characteristic sample dimension to successfully apply the law for slabs; Oliveira and Franca [7] found that this value was higher for cylinders than for slabs. Finally, Romano et al. [10] correctly applied Lambert’s law considering the geometry of the sample, thus finding the same field of application for both the shapes and a power concentration along the central axis for a cylindrical domain, according with the experimental observations [9] .

As regards fluids, there are fewer studies and the most are about batch processes, conducted only by means of Maxwell equations solution [11] [12] .

As regards continuous microwave heating, the problem of temperature distribution analysis in liquid foods that flow in cylindrical ducts was faced by Romano and Apicella [13] , taking into account also the velocity field and different rheological characteristics. In that work, generation term due to microwaves in the energy balance was modeled with Lambert’s law while in the present work, problem is faced also by solving the electromagnetic problem and results are then compared to check at what extent Lambert’s law estimation produces valid and accurate outcomes respect to Maxwell equation solution.

Sheryl Barringer and others [14] compared Lambert’s law and Maxwell equations solution on agar gel made slabs, observing that only exact solution of Maxwell equation can appreciate electric field oscillations; moreover, Lambert’s law seems to underestimate the absorbed power and, therefore, the temperature profile. However, an approximately exponential decay from surface towards the inner of sample is followed in both cases.

In the present paper, an analysis on liquid foods and a cylindrical geometry is made. The results of the two different models hypothesis are compared in terms of absorbed power and temperature profile.

2. Mathematical Model

2.1. Physical System

The physical system to be modeled is a cylindrical horizontal tube with a length

L = 0.3 m and a radius R = 0.02 m (Figure 1) in which different liquid foods are heated by microwave irradiance in the radial direction. A laminar flow is realized by a pressure drop between the inlet and the outlet sections. Only the axial component of the velocity is different from zero and it is a function of the radial direction : $v_z=v_z\left(r\right)$. The temperature, even though the microwave penetration is only radial, is also function of the axial direction by the effect of the flow field: $T=T\left(r,z\right)$.

The aim is to check the temperature profile and to compare results obtained with Lambert’s law and Maxwell equations resolution.

2.2. Balance Equations

The mathematical model consists of the following three differential equations in vector form with their boundary conditions [15] .

$\frac{D}{Dt}\rho =-\rho \left(\underset{\_}{\nabla }\cdot \underset{\_}{v}\right)$ (continuity equation) (1)

$\rho \frac{D}{Dt}\underset{\_}{v}=-\underset{\_}{\nabla }p-\left[\underset{\_}{\nabla }\cdot \underset{\_}{\underset{\_}{\tau }}\right]+\rho \underset{\_}{g}$ (momentum equation) (2)

$\rho c_p\frac{D}{Dt}T=k{\nabla }^{2}T+Q$ (energy equation) (3)

In cylindrical coordinate system, Equations (1) and (2) for the specific problem become:

$\frac{\partial v_z}{\partial z}=0$ (4)

$\rho \frac{\partial v_z}{\partial t}=-\frac{\partial }{\partial r}\left(r{\tau }_{rz}\right)+\frac{P_{in}-P_{out}}{L}$ (5)

with the following initial and boundary conditions.

As initial condition, the fluid is stationary.

$@t=0,v_z=0,\forall r,z$ (i)

As boundary conditions, we assume no slipping at the wall

$@r=R,v_z=0,\forall z,\forall t>0$ (ii)

And we consider the symmetry on the axis.

$@r=0,{\tau }_{rz}=0,\forall z,\forall t>0$ (iii)

Similarly, Equation (3) turns out to be

$\rho c_p\left(\frac{\partial T}{\partial t}+v_z\frac{\partial T}{\partial z}\right)=k\left[\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial T}{\partial r}\right)+\frac{{\partial }^{2}T}{\partial z^2}\right]+Q$ (6)

With the following initial and boundary conditions.

As initial condition, the temperature of the fluid is assumed to be uniform

$@t=0,T=T_0,\forall r,z$ (iv)

As boundary conditions for the radial direction we consider the symmetry respect the z-axis and the heat exchange at the wall by external convection.

$@r=0,\frac{\partial T}{\partial r}=0,\forall z,t$ (v)

$@r=R,k\frac{\partial T}{\partial r}=-hc\ast \left(T-T_{air}\right),\forall z,t$ (vi)

With initial and boundary conditions for the radial direction, so defined, the system results to be axial-symmetric.

For the z direction uniform temperature in the inlet and Danckwerts’s condition in the outlet section are imposed

$@z=0,T=T_{in},\forall r,t$ (vii)

$@z=L,\frac{\partial T}{\partial z}=0,\forall r,t$ (viii)

Initial temperature $T_0$ and input temperature $T_{in}$ are both equal to environmental temperature
$T_{air}$.

Heat generation due to microwaves can be modeled as the QLR predicted by Lambert’s law or as the QMR provided by the Poynting theorem after the resolution of Maxwell equations.

Lambert’s law forecasts an exponential decay inside the sample, Romano et al. (2005)

$QLR=2\alpha \cdot \left(\frac{R}{r}\right)\cdot Q_1\cdot \left[{\text{e}}^{-2\alpha \left(R-r\right)}+{\text{e}}^{-2\alpha \left(R+r\right)}\right]$ (7)

where

$Q_1=\frac{Q_0}{2\text{π}R\left(R+L\right)}$

$\alpha =\frac{2\text{π}}{{\lambda }_0}\sqrt{\frac{{\epsilon }^{\prime }\left\{{\left[1+{\left(\frac{{\epsilon }^{″}}{{\epsilon }^{\prime }}\right)}^2\right]}^{\frac{1}{2}}-1\right\}}{2}}$

$\frac{{\epsilon }^{″}}{{\epsilon }^{\prime }}$ is the loss tangent, known as tanδ.

The attenuation factor α for each fluid has been considered as a temperature function, calculated by interpolation starting from graphic relationships for dielectric constant and loss tangent versus temperature in a range 10˚C - 90˚C [16] .

Maxwell equations are a system of four vectory equations [17] [18] :

$\nabla \times E=-\frac{\partial B}{\partial t}$ Faraday’s law (8)

$\nabla \times H=\frac{\partial D}{\partial t}+J+\sigma E$ Ampère-Maxwell law (9)

$\nabla \cdot D=\rho$ Gauss’s law for electric field (10)

$\nabla \cdot B=0$ Gauss’s law for magnetic field (11)

Keeping in mind also the following constitutive equations:

$D=\epsilon E$ (12)

$B=\mu H$ (13)

Making the rotor of first and second member of (8)

$\nabla \times \left(\nabla \times E\right)=\nabla \times \left(-\mu \frac{\partial H}{\partial t}\right)$

Hence

$\nabla \left(\nabla \cdot E\right)-{\nabla }^{2}E=-\mu \frac{\partial }{\partial t}\left(\nabla \times H\right)$

Using Equations (9) and (10)

$\nabla \frac{\rho }{\epsilon }-{\nabla }^{2}E=-\mu \epsilon \frac{{\partial }^{2}E}{\partial t^2}-\mu \frac{\partial J}{\partial t}-\mu \sigma \frac{\partial E}{\partial t}$

${\nabla }^{2}E=\frac{1}{\epsilon }\nabla \rho +\mu \sigma \frac{\partial E}{\partial t}+\mu \epsilon \frac{{\partial }^{2}E}{\partial t^2}$

For source-free region ( $\rho =0$ ) and media with negligible electrical conductivity ( $\sigma \approx 0$ ) it becomes:

${\nabla }^{2}E-\frac{1}{c^2}\frac{{\partial }^{2}E}{\partial t^2}=0$ (14)

that is the wave equation in which $c=\sqrt{\frac{1}{\mu \epsilon }}$ is the speed of light (299,792,458 m/s).

It can be noticed that solution of the latter equation, that is second order respect to time, is equivalent to solution of the original system made by two first order equations. Passing from time domain to frequency domain:

${\nabla }^{2}E+k^{2}E=0$ (15)

where

$k=\omega \sqrt{\mu \epsilon }=\frac{\omega }{c}=\frac{2\text{π}f}{c}=\frac{\text{2π}}{\lambda }$

The equation assumes the following form:

${\nabla }^{2}E+{\omega }^2\mu \epsilon E=0$ (16)

Anyway, once solved for electric field, magnetic field H can be easily found using Equation (9).

After that electric and magnetic fields have been calculated, one can also determine the volumetric unit power generated by the correspondent wave thanks to the so called Poynting theorem.

Poynting vector $S[=]\frac{\text{W}}{{\text{m}}^{\text{2}}}$ is defined as the vector product $E\times H$ and represents the directional energy flux of an electromagnetic field.

Poynting’s theorem is a statement of conservation of energy for the electromagnetic field. In differential form, defining energy density as e:

$-\frac{\partial e}{\partial t}=\nabla \cdot S+J\cdot E[=]\frac{\text{W}}{{\text{m}}^{\text{3}}}$ (17)

$-\frac{\partial e}{\partial t}=\nabla \cdot \left(E\times H\right)+J\cdot E$

$-\frac{\partial e}{\partial t}=H\cdot \nabla \times E-E\cdot \nabla \times H+J\cdot E$

Using (8), (9), (10) and (11)

$\frac{\partial e}{\partial t}=\mu H\frac{\partial H}{\partial t}+\epsilon E\frac{\partial E}{\partial t}$

$\frac{\partial e}{\partial t}=\frac{1}{2}\frac{\partial }{\partial t}\left(\mu {\left|H\right|}^2+\epsilon {\left|E\right|}^2\right)$ (18)

In dielectric materials, the magnetic permeability is usually small and the first term of the equation second side can be neglected; so, shifting to the frequency domain, considering that is $\epsilon ={\epsilon }^{\prime }-j{\epsilon }^{″}$ and taking only the real component:

$\frac{\partial e}{\partial t}=P=\frac{1}{2}\omega {\epsilon }_0{\epsilon }^{″}{\left|E\right|}^2$ (19)

3. Materials and Methods

To solve the previous partial differential equations, a Finite Elements Method (FEM) [19] has been used. To practically implement this solution, COMSOL Multiphysics^® 5.3 has been utilized and in Figure 2 a possible mesh of the system is

shown.

Three fluid foods have been considered: skim milk, with a Newtonian behavior and constant viscosity, apple sauce and tomato sauce as non Newtonian fluids, modeled with a power law having different fluid consistency coefficient and flow behavior index [20] , as reported in Table 1 together with the other physical and rheological properties.

As ε' and ε" are temperature functions, average values have been obtained by integrating in the entire domain and in the time (range 0 - 100 s). They are fundamental for microwave heating, because they determine respectively the energy absorbed and the fraction converted in heat power.

4. Results

In the next pages, velocity and temperature profiles and generation term modeled with both Lambert’s law and Poynting theorem are reported.

In all cases, the average velocity is the same (4 mm/s). To practically realize that, as the three fluid present different rheological properties, different pressure drops, reported in Table 2, have been calculated from integration of the momentum flux distribution for the flow through a circular tube for Newtonian and non-Newtonian fluids.

To ensure a comparability between the two heat generation terms, a proper value of electric field has been chosen as boundary condition in wave equation so that the volume integral over all the domain is the same:

${\displaystyle {\int }_{V}QLR\text{d}V}={\displaystyle {\int }_{V}QMR\text{d}V}$

In this way, the total heat supplied in the different cases is the same. Values of total power P for all cases and correspondent values of electric field in case of Poynting generationare reported in Table 3.

All the results showed above are obtained for the same fixed time of 100 seconds. It can be seen that skim milk, that is a Newtonian fluid, shows a parabolic shape (Figure 3), while the other two fluids present a flatter profile. In particular tomato sauce (Figure 5) is flatter than apple sauce (Figure 4), having a lower flow behavior index value n.

The three fluids show a very similar temperature distribution even though apple sauce values are smaller (Figure 6, Figure 7). Temperature profiles can be explained observing the generation terms (Figure 8). Barringer’s result of

underestimation of heating in case of Lambert’s law use is confirmed also in this case with a different geometry.

Actually, in order to better understand the phenomenon, it has been considered that the fluid elements of each of the three liquids have different velocity profiles along the radial direction (Figures 3-5) and so different residence times depending to the radius. For the fluid elements having higher velocities, near the tube axis, the residence times are lower; for those near the wall, the velocities are lower and so the residence times are longer. In skim milk this difference is more important, whereas in the other two cases the profile is quite flat, especially in the case of tomato sauce. Anyway, outlet section (z = 30 cm) has been fixed and temperature profile evaluated for each position based on the real residence time are overlapped in Figure 9.

It is evident that in the three cases the trend is always the same: in the first part near the axis in which Lambert’s laws overestimates Maxwell results in term of temperature. This is due to Lambert’s generation term that goes to infinite for r = 0 that in turn is due to the pre exponential factor R/r. In the larger middle part, Lambert’s law curve is instead always under Maxwell predictions, save to slightly raise in the last part near the wall.

In both cases temperature increases near the wall because of the longer residence times, but in Lambert’s law case the effect is more evident. After all, apart few points, a good agreement between the two models has been found.

5. Conclusions

In conclusion, it can be stated that Lambert’s law allows to estimate power generated during a microwave heating process without solving for electric field and so in a relatively simply way, but it presents some problems like the singularity of infinite power for r = 0 and an appreciable underestimation of the real power and so of the temperature distribution. Moreover as it suppose an exponential decay, it can’t forecast the oscillations of the electric field and so of Poynting generation.

On the other hand, solution of Maxwell equations allows to exactly know the real electric field and then, thanks to Poynting theorem, the actual generated power, although it is obviously more expensive in terms of time and resource of calculations.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Apicella, R. and Romano, V. (2019) Comparison between Lambert and Maxwell Approaches in the Modelling of Microwave Heating of Liquid Foods. Engineering, 11, 1-13. https://doi.org/10.4236/eng.2019.111001

References