Microwave Heating in Moist Materials

The definition of “microwaves” is somewhat arbitrary; however microwaves are usually considered to be electromagnetic waves in the frequency range from 300 MHz to 300 GHz. Before World War II, there is little evidence of work on radio frequency or microwave heating; however Kassner (1937b) mentions industrial applications of microwave energy in two of his patents on spark-gap microwave generators (Kassner, 1937a, 1937b, 1938). Unfortunately early studies in radio frequency heating concluded that microwave heating of food stuffs would be most unlikely because the calculated electric field strength required to heat biological materials would approach the breakdown voltage of air (Shaw & Galvin, 1949). A fortuitous discover by Spencer that microwave energy could heat food (Murray, 1958) lead to a series of patents (Spencer, 1947, 1949, 1952) and the development of microwave cooking equipment. The major advantages of microwave heating are its short start-up, precise control and volumetric heating (Ayappa et al., 1991); however microwave heating suffers from: uneven temperature distributions (Van Remmen et al., 1996; Brodie, 2008); unstable temperatures (Vriezinga, 1996, 1998, 1999; Vriezinga et al., 2002); and rapid moisture movement (Brodie, 2007b). This chapter will explore the theory and practice of microwave heating in moist materials.


Introduction
The definition of "microwaves" is somewhat arbitrary; however microwaves are usually considered to be electromagnetic waves in the frequency range from 300 MHz to 300 GHz. Before World War II, there is little evidence of work on radio frequency or microwave heating; however Kassner (1937b) mentions industrial applications of microwave energy in two of his patents on spark-gap microwave generators (Kassner, 1937a(Kassner, , 1937b(Kassner, , 1938. Unfortunately early studies in radio frequency heating concluded that microwave heating of food stuffs would be most unlikely because the calculated electric field strength required to heat biological materials would approach the breakdown voltage of air (Shaw & Galvin, 1949). A fortuitous discover by Spencer that microwave energy could heat food (Murray, 1958) lead to a series of patents (Spencer, 1947(Spencer, , 1949(Spencer, , 1952 and the development of microwave cooking equipment. The major advantages of microwave heating are its short start-up, precise control and volumetric heating (Ayappa et al., 1991); however microwave heating suffers from: uneven temperature distributions (Van Remmen et al., 1996;Brodie, 2008); unstable temperatures (Vriezinga, 1996(Vriezinga, , 1998(Vriezinga, , 1999Vriezinga et al., 2002); and rapid moisture movement (Brodie, 2007b). This chapter will explore the theory and practice of microwave heating in moist materials.

Theory
Microwave heating is governed by Maxwell's electromagnetic equations: Because the angular frequency ω is so large at microwave frequencies, equation (4) can be nicely approximated by: This can be rearranged to become: The form of the Laplacian operator (∇ 2 ) in equation (3) depends entirely on the geometry of the problem.

Simultaneous heat and moisture movement
Any realistic analysis of microwave heating in moist materials must account for simultaneous heat and moisture diffusion through the material. The coupling between heat and moisture transfer is well known but not very well understood (Chu & Lee, 1993). Henry www.intechopen.com (1948) first proposed the theory for diffusion of heat and moisture into a textile package. Crank (1979) presented a more thorough development of Henry's work. This theory has been rewritten and used by many authors (Chu & Lee, 1993;Vos et al., 1994;Fan et al., 2000;Wang et al., 2002;Fan, 2004;Brodie, 2007b). The amount of water vapour moving into a small section of porous material is the sum of any net increase in moisture content in the air space and the net increase in moisture content of the material's fibres (Crank, 1979). Therefore: Heat is evolved because of microwave heating and moisture absorption by the material; therefore the thermal diffusion equation, which includes a volumetric heat source, is: If a linear relationship between the moisture content of a material, the moisture vapour concentration in the air spaces in the material and the temperature is assumed (Henry, 1948;Crank, 1979), then: where σ and ω are constants of association between heat transfer and moisture vapour concentration. Substituting from equation (11) into equations (9) and (10) This can be expressed in a simpler form if Ω = pM v + nT: The constants of association, p and n, are calculated to satisfy: Equation (14) implies that the combined heat and moisture diffusion coefficient (γ) has two independent values. This is consistent with Henry's (1948) equation for simultaneous heat and moisture diffusion in textiles. Henry (1948) showed that a change in external temperature or humidity (or both), "results in a coupled diffusion of moisture and heat". Henry (1948) also states that "the diffusion constants www.intechopen.com appropriate to these two quantities are always such that one is greater and the other less than either of the diffusion constants which would be observed for the moisture or heat, were these not coupled by the interaction". The slower diffusion coefficient of the coupled system is less than either the isothermal diffusion constant for moisture or the constant vapour concentration coefficient for heat diffusion, whichever is less, but never by more than one half (Henry, 1948). The faster diffusion coefficient is always many times greater than either of these independent diffusion constants (Henry, 1948). Considerable evidence exists in literature for rapid heating and drying during microwave processing (Rozsa, 1995;Zielonka & Dolowy, 1998;Torgovnikov & Vinden, 2003); therefore it is reasonable to assume that the faster diffusion wave dominates microwave heating in moist materials. A slow heat and moisture diffusion wave should also exist; however observing this slow wave may be difficult and no evidence of its influence on microwave heating has been seen in literature so far.

Solutions to the microwave heating equation
Approximate solutions to equation (13) are sought under assumptions of uniform material properties for rectangular, cylindrical and spherical geometries. The aim is to readily identify the influence of geometry over microwave fields and temperature distributions. It is recognised that the derived equations will only provide approximate solutions to the diffusion problem, however the solutions are sufficiently accurate to provide some insight into the influence of geometry over temperature distributions during microwave heating.

Microwave heating in rectangular coordinates
Beginning with Maxwell's equations, Ayappa et al. (1991) demonstrate that the average power density (W m −2 ) at a distance z below the irradiated surface of a material block of thickness W is: The Laplace transformation of equation (13) is: This can now be solved using standard mathematical tables (Crank, 1979) Applying the first boundary condition to evaluate A yields: The analysis used to obtain equation (23) assumed that the incoming electromagnetic wave has a constant phase. In a multi-mode oven, this scenario is unlikely because the load is usually placed on a rotating turntable, or mode stirrers are employed to perturb the microwave fields. The net result is that the phase (δ) of the incoming wave varies considerably during the heating process. For simplicity, it has been assumed that δ regularly cycles through 2π radians during the heating process. Therefore: This renders Ψ = 0, allowing equation (23) to simplify to:  (23) and (29) predict that when a slab or block is heated in a microwave oven, subsurface heating occurs with the maximum temperature located slightly below the material surface. This agrees with information presented by Van Remmen et al., (1996) and Zielonka and Gierlik (1999). When the thickness of the material is small, the maximum temperature occurs in the center of the heated slab because the "subsurface" temperature peaks coincide with the center of the slab. However when the thickness of the material increased, two temperature peaks appear. This is a direct result of field attenuation inside the material. When the slab or block is fairly small, very little wave attenuation occurs; however, when the slab or block becomes much larger, the fields are significantly attenuated near the center of the object.

Microwave heating in cylindrical coordinates
Temperature profiles in small-diameter cylinders usually exhibit pronounced core heating (Ohlsson & Risman, 1978;Van Remmen et al., 1996); on the other hand, temperature profiles in large cylinders exhibit subsurface heating, with the peak temperature occurring slightly below the surface (Van Remmen et al., 1996). The same is also true when the microwave fields can not penetrate very far into the cylinder (Van Remmen et al., 1996).
The refractive index, at microwave frequencies, for many moist materials is very high (refractive index > 8). Therefore, by Snell's law, it is reasonable to assume that any wave propagates radially inside the cylinder irrespective of the incident angle of the external wave. It is also assumed that the microwave oven's turntable or mode stirrer ensures that the long-term illumination of the cylinder is uniform over all exposed surfaces. Therefore, in cylindrical coordinates, equation (3) becomes: The solution to equation (30) has the general form: For E to be bounded at the center of the cylinder, B must be equal to zero, and therefore E = AJ o (fr). Oliveira and Franca (2002) present the same solution for the electric field in a cylinder.
The Bessel function addition theorem is now applied to this complex field: When the complex argument is small, this simplifies to: Yousif and Melka (1997) further show that when the real part of the argument is small, as would be the case when heating relatively small cylindrical objects in domestic microwave ovens, Jx j y Iy j x Iy +≈ − . Therefore, except near the central axis of the cylinder, the real component of the electric field has the form E = AI o (βr). When the electric field in the cylinder is normalized to the field at the surface of the cylinder, equation (31) becomes: Therefore, the heat and moisture diffusion equation in cylindrical coordinates is: Now γ is usually very small (≈10 -7 ), so except in the case when r is infinitesimally close to zero, equation (37) can be approximated by: Therefore, the combined solution becomes: This can be solved using standard mathematical tables (Crank, 1979) To evaluate the full temperature/moisture distribution requires application of the following boundary conditions: For the sake of simplifying the mathematics somewhat, the second boundary condition is applied to: Because the approximation for the microwave power distribution used in equation (34) is not accurate near the center of the cylinder, the first boundary condition should be applied to: The analysis eventually yields: Following the same process as with the cylindrical analysis described above, the solution to equation (51) becomes: Vriezinga (1998) presents an expression for the radial component of the spherical field in terms of ordinary Bessel functions: Jf r j f r fr ; therefore, equation (52)  Because cylinders and spheres tend to focus the microwave energy into the center of the workload, there is a stronger tendency towards core heating than in the rectangular slab. However, as the radius of the cylinder or sphere increases, the attenuation of the microwave fields becomes more pronounced, and there is a shift from predominantly core heating to predominantly subsurface heating. This phenomenon is well documented (Ohlsson & Risman, 1978;Kritikos et al., 1981;Van Remmen et al., 1996;Romano et al., 2005). Even when the radius of the cylinder or sphere is fairly small, if the dielectric loss factor increases, there is a transition from predominantly core heating at lower loss factors to predominantly subsurface heating at higher loss factors. This was also observed by Van Remmen et al. (1996).

Multi-dimensional solutions
So far, this analysis has only considered single-dimensional heat transfer. Most real problems are multi-dimensional. The solution to multi-dimensional partial differential equations is a simple product of the single-dimensional solutions (Crank, 1979;Van Remmen et al., 1996). For example, the temperature distribution within a rectangular block, irradiated equally on all six faces, can be found by multiplying the temperature profiles associated with three slabs of thickness equal to the length, width, and height of the rectangular block. Fig. 1 (a) shows the predicted temperature distribution in a rectangular block of agar gel heated in a multi-mode microwave oven. This distribution agrees with thermal images and other data found in the literature (Van Remmen et al., 1996;Zielonka & Gierlik, 1999). Agar gel is a useful medium for identifying hot spots during microwave heating because it melts at the hot spots, but remains solid where the temperature of the object is lower than its melting temperature. The theoretical temperature distributions in the agar gel block compare very well with the melting patterns in a block of agar gel ( Fig. 1 (b)).

Length (m)
Width ( The temperature distribution within a cylinder, which is irradiated equally on all surfaces, can be found by multiplying the radial temperature profile for the cylinder by the temperature profile associated with a slab of thickness equal to the height of the cylinder. When the cylinder's radius and loss factor are small, the temperature distribution resembles a "dumbbell" with two temperature peaks along the longitudinal axis, as shown in fig. 2a. As the loss factor of the cylinder increases, the temperature profile is transformed such that there is an annulus of high temperature below the upper and lower circular surfaces, as shown in Fig. 3b. The same is true when the outer radius of the cylinder increases, as shown in Fig. 3c. The predicted temperature profiles shown in Fig. 3a and 3b are directly comparable with thermal images for cylindrical objects presented by Van Remmen et al. (1996). The temperature distribution within a sphere irradiated equally on all surfaces is purely radial. Fig. 4 shows the temperature distributions in spheres of various diameters and loss factors. Again, these predicted temperature profiles are directly comparable with thermal images for spherical objects presented by Van Remmen et al. (1996).

The effect of heating time on temperature and drying
As would be expected, increasing microwave heating time increases the temperature in moist objects (Fig. 5). However, the relationship between time and temperature at a particular location inside the heated material is often non-linear because of the complexity of synchronized heat and moisture diffusion, changes of the spatial distribution of heat with time, partial drying of the material and their combined effect on the absorption of microwave energy.   Heating time not only affects the magnitude of the temperature but also the spatial distribution of heat. Fig. 6 shows a heating sequence for a rectangular block, while Fig. 7 shows the heating sequence for a cylinder. The change in temperature distribution as heating time increases is due to internal thermal diffusion from the zones where the microwave energy is being absorbed by the material to the rest of the object. As may also be expected, microwave heating also dries moist materials. Conventional drying usually consists of a short constant rate drying period, where the rate of drying is controlled by surface evaporation, followed by a prolonged falling rate drying period (Fig.  8 a) during which the slow movement of moisture out of the object is controlled by the moisture diffusion properties of the material (Youngman et al., 1999). However, after a relatively short heat up period, microwave drying is usually quite linear with time ( Fig. 8 b). This is because of the strong coupling between heat and moisture movement through the material. Therefore the rate of microwave drying is not limited by the internal diffusion properties of the material, but by the rate at which moisture can be evaporated from the surface of the material. As the material dries, its dielectric properties will change significantly. These changes are too slow to affect the behaviour of individual cycles of the microwave fields; however they do affect microwave heating on the much longer thermal time scale.

The effect of changing dielectric properties on microwave heating
As an example of how the dielectric properties of moist materials change with moisture content, Ulaby and El-Rayes (1987) Plant material with high moisture content have higher dielectric constants (Fig. 9) and will therefore interact more with microwave fields. It has also been well documented that the dielectric properties of most other materials are temperature and moisture dependent (Hill & Marchant, 1996;Vriezinga, 1998;Nelson et al., 2001). In order to explore what can happen when a material dries out during microwave heating, a semi-analytical technique can be used to analyse the microwave heating problem.
As an example, the heating in a cylinder of plant material will be considered: www.intechopen.com Equation (47) was used in an iterative calculation, where the temperature and moisture content of a plant based material was calculated after 0.5 second of microwave heating. This new temperature and moisture content was used to calculate the dielectric properties of the material using equation (59). These new dielectric properties were then used to calculate the next rise in temperature after another 0.5 seconds of microwave heating and so on. Moisture loss was assumed to have a linear relationship with microwave heating time as shown in Fig. 8 b. The results of this process are presented in Fig. 10. The heating rate for 10 mm and 12 mm diameter cylinders is relatively constant with time; however there is a sudden temperature jump in the 15 mm diameter case when there is no change in the applied microwave power (Fig. 10). This sudden jump in temperature is the result of "thermal runaway". Thermal runaway often destroys or chars the treated material, as seen in Fig. 11.   Fig. 11. Charing, due to thermal runaway, in the core of a wood sample that was being dried in a microwave oven (Note: the location of the charring in comparison with the last image in Fig. 7)

Thermal runaway
Several authors (Vriezinga, 1998;Liu & Marchant, 2002) have suggested that the relationship between the applied microwave field strength and temperature follows a characteristic "S" shaped curve (Fig. 12). As microwave power is steadily increased, temperature rises steadily along the stable lower arm of the power curve shown in Fig. 12. In the critical power range, a small increase in electric field strength may shift the equilibrium temperature from the lower limb of the temperature curve to the higher limb, as indicated by the arrow in Fig. 12. The resulting change in temperature could be quite substantial and very rapid. This temperature jump gives rise to the potentially catastrophic phenomenon known as thermal runaway (Vriezinga, 1998). Thermal runaway is predicted by equations (23), (48) and (56). In each case, under the right combination of physical dimensions and dielectric properties of the material, the denominators of these equations can become very small thus forcing the temperature to become very large. Thermal runaway manifests itself as a sudden jump in temperature over a very short time. It is linked to the temperature and moisture dependence of the material's dielectric properties.

Microwave Field Intensity (V m
In some cases this can result in the material absorbing more microwave energy as it is heated, which in turn leads to faster heating which in turn leads to more energy absorption, etc. The net result is a sudden and uncontrolled jump in temperature, depending on the strength of the electromagnetic fields and the heating time. The transfer of microwave energy into a material depends on the transmission coefficient of the material's surface, which also depends on the dielectric properties of the material. In many moist materials, the dielectric constant decreases with temperature and moisture loss (Vriezinga, 1998). As the dielectric properties decrease with microwave heating, more energy enters the material due to increased transmission across the material's surface. This leads to faster heating, which leads to more energy transfer and so on. Thermal runaway can also be caused by resonance (Vriezinga, 1999) of the electromagnetic waves inside the object. As dielectric properties decrease with increasing temperature and moisture loss, the wavelength inside the material will increase in length. At a certain temperature and moisture, the electromagnetic wave will fit exactly within the dimensions of the heated object, causing resonance (Vriezinga, 1999). Exactly at that moment the temperature will suddenly rise. Because heating and drying rates are dependent on the microwave field intensity, thermal runaway is strongly dependent on the intensity of the microwave's electric field strength (Fig. 13).

Examples of using thermal runaway in industrial microwave processes
In most cases, thermal runaway is a problem becasue it usually leads to destruction of the material (Zielonka et al., 1997); however it has been very effectively used in some applications: Jerby et al. (2002) have developed a microwave drill that can drill holes through glass and ceramics by super-heating a very small section of the material using microwave induced thermal runaway. The system works by creating very intense microwave fields immediately in front of a needle like probe that extends into the material as the drilling process proceeds. Vinden and Torgovnikov (2000) have also shown that the controlled application of intense microwave fields to green timber can rapidly boil free water inside the wood cells (Vinden & Torgovnikov, 2000;Ximing et al., 2002) causing microscopic fractures in the wood cells. When exposure to these intense fields is carefully managed, there is very little change in the vissible appearance or strength of the wood; however there is a substantial change in wood density.

Time (s)
E field = 1000 E field = 1800 E field = 3000 Fig. 13. Theoretical influence of microwave field strength over thermal runaway in a cylinder of plant based material of 16 mm diameter Lawrence (2006) has demonstrated that this process can reduce the density of Eucalyptus obliqua wood by up to 12 %, as the wood cells are ruptured and the cross section of the wood is expanded due to internal steam explosions. This has been confirmed by other experimental work (Brodie, 2005(Brodie, , 2007a. This reduction in wood density leads to rapid drying in conventional systems with little loss in wood strength or quality. Microwave induced cell rupture has been linked to thermal runaway (Brodie, 2007a) in the timber.

Determining the electric field strength
The magnitude of the temperature/moisture vapour wave and the onset of phenomenon such as thermal runaway depend on the electric field strength at the surface of the heated object; therefore it is important to model the field strength in the immediate vicinity of the heated object. Analytical solutions to Maxwell's equations can be found for simple microwave applicators, especially when they are empty (Meredith, 1994); however modeling of the microwave field strength in a loaded applicator requires the solution of Maxwell's equations in three dimensions when complex boundary conditions are imposed onto the system. It is almost impossible to determine these fields using analytical techniques. Several techniques can be employed to simulate microwave fields in complex systems, including numerical techniques such as the Finite Difference Time Domain (FDTD) technique, proposed by Yee (1966). The FDTD method is a simple and elegant way to transform the differential form of Maxwell's equations from equation (1) into difference equations. Yee used an electric field grid (Fig. 14), which was offset both spatially and temporally from a magnetic field grid to calculate the present field distribution throughout the computational domain in terms of the past field distribution.
Equations (64) and (65) calculate the electric and magnetic fields in the x direction. Similar equations are required to calculate Ey, Ez, Hy and Hz. These equations are used in a leapfrog scheme to incrementally march the electric and magnetic fields forward in time; therefore this numerical technique is a simulation of the microwave field rather than a direct solution of the field equations in space and time.
The FDTD method was originally developed to solve electromagnetic transmission problems associated with the communications industry (Kopyt et al., 2003), but it is now commonly used to determine electromagnetic field distributions within industrial microwave applicators (Kopyt et al., 2003). Numerical simulations of the microwave fields can provide valuable insight into the behaviour of microwave heating systems. www.intechopen.com

Strengths of the FDTD method
Every modeling technique has some strengths and some weaknesses. FDTD is a very versatile modeling technique. It is fairly intuitive, so users can easily understand what to expect from a given model. FDTD is a time domain technique, and when a time-domain pulse (such as a Gaussian pulse) is used as an input to the computational model, a wide frequency range can be explored in a single simulation. This is extremely useful in applications where resonant frequencies are not known or when broadband performance is desirable.
Since FDTD is a time-domain technique which finds the E/H fields everywhere in the computational domain, it lends itself to providing animation displays (movies) of the E/H field movement throughout the model. FDTD also allows the user to specify the material properties at all points within the computational domain, which provides useful insight into the field distributions inside dielectric materials as they are processed using microwave systems. For example, when the FDTD technique is applied to a domestic microwave oven cavity in which various sized samples of water are included in the simulation space, it becomes clear that the microwave fields are significantly attenuated as the dimensions of the water sample increase (Fig. 15). This is becasue the water sample reduces the quality factor (a) (b) (c) (d) Fig. 15. FDTD simulations of the microwave fields inside a domestic microwave oven with various sized water loads (a) empty cavity, (b) a small water load, (c) a medium sized water load and (d) a large water load (Q) of the field resonance inside the chamber and effectively dampens the wave resonance. This field damping will affect the microwave heating time needed for a sample to reach a desired temperature. Consequently, small samples will heat much faster per unit volume than larger samples. Similarly, when the FDTD technique is applied to a domestic oven cavity in which a disk of wood with varying moisture content is placed, it becomes apparant that as the moisture content reduces, the microwave field amplitude inside the wood sample increases (Fig. 16). This supports the earlier discussion about wave resonance inside the heated material that may lead to thermal runaway. Fig. 16. FDTD simulations of the microwave fields inside a domestic microwave oven with a wooden disk of varying moisture content, measured as a percentage of oven dried wood.

Problems with the FDTD method
If aliasing of the final solution is to be avoided (Vandoren, 1982), the grid size must be small enough so that the electromagnetic field does not change significantly from one grid cell to the next. Similarly the time steps for each computational cycle must satisfy the Courant stability criterion (Yee, 1966;Taflove, 1988;Walker, 2001). This implies that the simulated time step for many problems may be no more than a few picoseconds. Therefore solving heating problems that may span many seconds or minutes of real time requires a significant number of computational iterations.

www.intechopen.com
For computational stability to be satisfied (Yee, 1966;Walker, 2001), the time step used in the program must satisfy: When establishing an appropriate grid for analysis, it must be remembered that the wavelength of the microwaves inside a dielectric material is: The values of Δx, Δy and Δz must be small compared to λ to ensure accurate simulations.
Since FDTD requires that the entire computational domain be gridded, and these grids must be small compared to the smallest wavelength and smaller than the smallest feature in the model, very large computational domains can be developed. This implies that very long calculation times are needed to simulate real microwave heating systems.
To put this into perspective, simulating the first 20 nanoseconds of electromagnetic activity inside a 335 mm by 335 mm by 205 mm microwave oven that has a cup (cylinder) of water in it, using a 4 mm grid, requires approximately 30 minutes of uninterrupted computational time on a Pentium 4 computer. Trying to simulate 10 minutes of real time would take 950 years on the same machine using the same code. FDTD finds the E/H fields everywhere in the computational domain. If the field values at some distance from the microwave source (like 10 meters away) are required, the computational domain will be excessively large. Far field extensions, which apply a Fast Fourier Transformation to the field distribution in a fixed plane of the FDTD computational domain are available, but this usually requires some post processing rather than direct calculations.
Most authors use FDTD and other numerical techniques to analyse "static systems". This implies that apart from the initial field excitation imposed onto the system, all other components are stationary. Commercial microwave applicators use mode stirrers (Metaxas & Meredith, 1983) or movement of the heated product relative to the microwave fields to deliberately perturb the fields and more evenly irradiate the load. Consequently, any attempt to evaluate the electric field, including numerical techniques, only provides an approximation of the field strength (Kopyt et al., 2003).
More recently, microwave heating has been used to modify the structure of wood (Vinden & Torgovnikov, 2000;Torgovnikov & Vinden, 2005) allowing it to dry much faster (Brodie, 2007a) and be impregnated with resins to make it stronger (Przewloka et al., 2007). Microwave energy has also been used to manage weeds and other soil born pathogens (Nelson, 1987;Bansal, 2006) and to treat various animal feeds to improve their digestibility (Sadeghi & Shawrang, 2006Brodie et al., 2010). Microwaves are commonly used to process ceramics (Takayama et al., 2005) and can be used to drill through ceramics and glass (Jerby et al., 2002).

Conclusion
The major advantages of microwave heating are its short startup, precise control, and volumetric heating (Ayappa et al., 1991); however non-uniform heating, rapid drying and unstable temperatures are commonly reported. This chapter has explored some of the reasons why these three features of microwave heating occur. Mathematical models and simulation techniques reveal that object geometry, size, moisture content, dielectric properties, and heating time all profoundly affect temperature/moisture distributions during microwave heating. Small slabs, cylinders, and spheres, with low dielectric loss, exhibit pronounced core heating, while increasing the size of the object or the dielectric loss of the material results in pronounced subsurface heating. = phase shift of microwave fields at the surface of a material ε* = complex electrical permittivity of the space through which the waves are propagating ε s = the dielectric properties of a material at very low frequencies γ = combined diffusivity for simultaneous heat and moisture transfer κ'

Nomenclature
= relative dielectric constant of the material κ" = dielectric loss factor of the material λ = wave length inside a material (m) λo = wave length in free space (m) μ = magnetic permeability of the space through which the waves are propagating θ A = phase angle associated with forward propagating wave in rectangular coordinate system θ Β = phase angle associated with reverse propagating wave in rectangular coordinate system ρ = composite material density (kg m -3 ) ρ c = charge density within the space through which the waves are propagating ρ s = density of the solid material (kg m -3 ) σ = constant of association relating moisture vapor concentration to moisture content in a solid σ c = electrical conductivity of the material (S m -1 ) τ = transmission coefficient for incoming microwave τ v = tortuosity factor ω = angular frequency (rad s -1 ) ℜ(z) = real part of the complex number z www.intechopen.com