Impact of Filled Materials on the Heating Uniformity and Safety of Microwave Heating Solid Stack Materials

Microwave heating of solid stack materials is common but bothered by problems of uneven heating and electric discharge phenomena. In this paper, a method introducing fluid materials with different relative permittivity is proposed to improve the heating uniformity and safety of solid stack materials. Simulations have been computed based on the finite element method (FEM) and validated by experiments. Simulation results show that the introducing of fluid materials with proper relative permittivity does improve the heating uniformity and safety. Fluid materials with the larger real part of relative permittivity could obviously lower the maximum modulus value of the electric field for about 23 times, and will lower the coefficient of variation (COV) in general, although in small ranges that it has fluctuated. Fluid materials with the larger imaginary part of relative permittivity, in a range from 0 to 0.3, can make a more efficient heating and it could lower the maximum modulus value of the electric field by 34 to 55% on the whole studied range. However, the larger imaginary part of relative permittivity will cause worse heating uniformity as the COV rises by 246.9% in the same process. The computed results are discussed and methods to reach uniform and safe heating through introducing fluid materials with proper relative permittivity are proposed.


Introduction
Unlike traditional heating methods that heat materials by radiation, convection and conduction, microwaves heat materials through the direct interaction with the inner polar molecules and charged particles of materials [1].Thus, characterized as efficient heating, volumetric heating and selective heating [2][3][4][5][6], microwave heating has been rapidly developed and applied in various areas.
However, there are also some natural drawbacks of microwave heating: (a) When heating materials by microwave, there exist "hot spots" and thermal runaway problems in the system temperature [7], which will cause degradation of the processing materials [8][9][10][11] and even cause burning and explosion of the microwave reactor and reactants [12,13].(b) When dealing with numerous solid materials, their stacking in the heating cavity can easily result in some sharp edges, tips or submicroscopic irregularities, which may lead to electric sparks or electric arcs [14,15].These electric discharges will affect the heating process and the product composition, and lead to safety problems when dealing with materials that have a low flashing point or in high temperature conditions [16,17].
For studies on microwave heating uniformity improvement, the common methods could be categorized into two aspects [18]: (1) To improve the uniformity of the electromagnetic field in the microwave cavity, such as optimization of microwave power feeding situations [19][20][21][22], introducing of mode stirrers [23,24], application of conductive matter [25,26] and optimization of heating chamber structure [27,28]; (2) to improve the uniformity of microwave energy absorption, such as introducing of a rotating turntable or conveyor belt [29,30], optimization of the shapes and sizes of processing materials [26,31,32].For studies on the electric discharge phenomena, related research has focused on microwave discharges in sintering of powdered metals.Attention is mainly paid to factors that influence the intensity and frequency of microwave discharges, such as the magnetron output power [14,33], the size, amount, conductivity, morphology and surface conditions of the metal [14,[33][34][35] and the dielectric properties of the surrounding medium [15,[35][36][37].The lack of comprehensive consideration of the heating uniformity and safety of microwave heating solid stack materials has hindered the further application of microwave heating.
In this paper, a method to improve the heating uniformity and safety of microwave heating solid stack materials by introducing fluid materials with different dielectric properties is proposed.In Section two, a general model is built and computed with coupled physics of electromagnetic and heat transfer in solid and liquid.Based on simulation conditions, corresponding physical experiments are designed to complete the validation work.In Section three, comparisons between experimental results and computational ones are performed.Following the analysis of the influence of fluid materials' complex relative permittivity on the heating process, solutions to uniformity and unsafety issues are discussed.

Geometry
To build the general model for analyzing, the simulation model should be characterized with a simple structure and be intuitive to show the difference of heating results.Hence, a standard BJ22 waveguide with one end shortened is used as the heating cavity.The processing materials, silicon carbide spheres, and the fluid materials are placed in the cavity in a quartz cup.The model structure is built in a multi-physics software, COMSOL Multiphysics (5.3, COMSOL Inc., Stockholm, Sweden), which is shown in Figure 1.
Processes 2018, 6, x FOR PEER REVIEW 2 of 13 problems when dealing with materials that have a low flashing point or in high temperature conditions [16,17].
For studies on microwave heating uniformity improvement, the common methods could be categorized into two aspects [18]: (1) To improve the uniformity of the electromagnetic field in the microwave cavity, such as optimization of microwave power feeding situations [19][20][21][22], introducing of mode stirrers [23,24], application of conductive matter [25,26] and optimization of heating chamber structure [27,28]; (2) to improve the uniformity of microwave energy absorption, such as introducing of a rotating turntable or conveyor belt [29,30], optimization of the shapes and sizes of processing materials [26,31,32].For studies on the electric discharge phenomena, related research has focused on microwave discharges in sintering of powdered metals.Attention is mainly paid to factors that influence the intensity and frequency of microwave discharges, such as the magnetron output power [14,33], the size, amount, conductivity, morphology and surface conditions of the metal [14,[33][34][35] and the dielectric properties of the surrounding medium [15,[35][36][37].The lack of comprehensive consideration of the heating uniformity and safety of microwave heating solid stack materials has hindered the further application of microwave heating.
In this paper, a method to improve the heating uniformity and safety of microwave heating solid stack materials by introducing fluid materials with different dielectric properties is proposed.In Section two, a general model is built and computed with coupled physics of electromagnetic and heat transfer in solid and liquid.Based on simulation conditions, corresponding physical experiments are designed to complete the validation work.In Section three, comparisons between experimental results and computational ones are performed.Following the analysis of the influence of fluid materials' complex relative permittivity on the heating process, solutions to uniformity and unsafety issues are discussed.

Geometry
To build the general model for analyzing, the simulation model should be characterized with a simple structure and be intuitive to show the difference of heating results.Hence, a standard BJ22 waveguide with one end shortened is used as the heating cavity.The processing materials, silicon carbide spheres, and the fluid materials are placed in the cavity in a quartz cup.The model structure is built in a multi-physics software, COMSOL Multiphysics (5.3, COMSOL Inc., Stockholm, Sweden), which is shown in Figure 1.

Governing Equations
The whole simulation is coupled of two physical processes: The microwave propagation process and heat transfer in solid and liquid process.

Governing Equations
The whole simulation is coupled of two physical processes: The microwave propagation process and heat transfer in solid and liquid process.
For the microwave propagation process, Maxwell's equations are used to solve the distribution of the electric field.The governing equations here are given as: where → H is the magnetic field intensity, → J is the ampere density, → E is the electric field intensity, t is the time, → B is the magnetic induction intensity, → D is the electric displacement vector and ρ e is the electric charge density.Equation (1) in the time harmonic field can be then written as the Helmholtz equation [38] ∇ where µ r is the relative permeability, k 0 is the wave number in free space, ε r is the relative permittivity, ε 0 is the permittivity of vacuum, ω is the angular frequency, σ is the electrical conductivity and µ 0 is the permeability of vacuum.Then, the electromagnetic power loss Q e of the processing materials can be gained from the computed electric field by the following equation [39,40] where ε is the real part of the relative permittivity, ε is the imaginary part of the relative permittivity and tan δ is the loss tangent of the processing materials.
For the heat transfer process, only the processing materials are involved to reduce the computation cost.The governing equation for heat transfer in solid and fluid is given as [41][42][43] where ρ is the material density, C p is the material heat capacity under atmospheric pressure, T is the temperature, Q is the heat source and k is the thermal conductivity.

Input Parameters and Boundary Conditions
To complete the simulation, necessary property parameters and boundary conditions are needed.Related input parameters are shown in Table 1.The thermal and dielectric properties of the processing materials are obtained from related literature [44][45][46][47] and modified by comparisons between experimental results and simulation results.The loss tangent of water (deionized), glycerol (analytical reagent, Yangzhou Feiyang chemical industry Co., Ltd, Yangzhou, China) and silicon carbide (99% purity, Hai Ning Zhijie pottery bearing Co., Ltd, Haining, China) is compared through the references [44][45][46] and the experiments at different temperatures, and the optimized data is chosen.It should be noted that the conductivity of silicon carbide is included by the loss tangent.The complex relative permittivity is defined as the processing materials to give a near-actual model, which is expressed as [48].For the microwave propagation process, one surface of the cavity is selected as the port where it feeds the electromagnetic energy at 2450 MHz in TE10 mode.Other surfaces are all defined as a perfect electric conductor, which can be expressed as the equation: where → n is the unit normal vector of the corresponding surface.For the heat transfer process, as only the processing materials are involved, all boundaries of the processing materials are defined as the thermal insulation boundary condition.The governing equation is given as where → q is the heat flux.

Mesh Size
Appropriate mesh size in the model simulation can provide accurate simulation results with higher computation efficiency.The space discretization errors could be down to a quarter when mesh size is halved, while the computation time will increase by almost 16 times [49].
To determine the appropriate mesh size in our model, normalized power absorption (NPA) of the processing materials has been employed to complete the mesh independent study [37].Here, the variation of NPA with mesh sizes is shown in Figure 2. The manual of software QuickWave (QWED, Warsaw, Poland) suggests there to be 12 cells per wavelength for mesh independent results while other researchers [50] suggest 10 cells per wavelength is enough.The mesh size used in this paper is defined as For the microwave propagation process, one surface of the cavity is selected as the port where it feeds the electromagnetic energy at 2450 MHz in TE10 mode.Other surfaces are all defined as a perfect electric conductor, which can be expressed as the equation: where n is the unit normal vector of the corresponding surface.
For the heat transfer process, as only the processing materials are involved, all boundaries of the processing materials are defined as the thermal insulation boundary condition.The governing equation is given as n q 0    (7) where q is the heat flux.

Mesh Size
Appropriate mesh size in the model simulation can provide accurate simulation results with higher computation efficiency.The space discretization errors could be down to a quarter when mesh size is halved, while the computation time will increase by almost 16 times [49].
To determine the appropriate mesh size in our model, normalized power absorption (NPA) of the processing materials has been employed to complete the mesh independent study [37].Here, the variation of NPA with mesh sizes is shown in Figure 2. The manual of software QuickWave (QWED, Warsaw, Poland) suggests there to be 12 cells per wavelength for mesh independent results while other researchers [50] suggest 10 cells per wavelength is enough.The mesh size used in this paper is defined as

Simulation Process
In the simulation, the electromagnetic field distribution is firstly calculated in the frequency domain, and then the dissipated power is calculated.Finally, the temperature rise of the material is updated through the heat transfer equation in the time domain.Since the electrical properties of samples do not vary with temperature, the electromagnetic field distribution will not change and is only calculated once.The whole heating process lasts 120 s, and the computation time of the model depends on the mesh size and the memory of the computer.In this simulation, the mesh sizes of the model with water and model with glycerol are 808,998 and 455,672, respectively.The computation time of each model is 7317 s and 3354 s in a computer with a memory of 128 GB.

Experiment System
To validate the computation results, physical experiments need to be completed.The whole experimental system is shown in Figures 3 and 4. Microwave power is generated from a microwave generator (WSPS-2450-1000M, Wattsine, Chengdu, China) to a waveguide coaxial connector (CAWG-26-N, Euler, Nanjing, China).A circulator connected with water load is employed to protect the source from the reflection power.With a dual directional coupler and a power meter (AV2433, the 41st Institute of CETC, Tsingtao, China) as power monitoring, microwave power is finally fed into the BJ22 waveguide connected with a slide short.Final heating temperature is gained by a fiber optic thermometer (FISO FOT-NS-967A, FISO Technologies, Quebec, QC, Canada).

Simulation Process
In the simulation, the electromagnetic field distribution is firstly calculated in the frequency domain, and then the dissipated power is calculated.Finally, the temperature rise of the material is updated through the heat transfer equation in the time domain.Since the electrical properties of samples do not vary with temperature, the electromagnetic field distribution will not change and is only calculated once.The whole heating process lasts 120 s, and the computation time of the model depends on the mesh size and the memory of the computer.In this simulation, the mesh sizes of the model with water and model with glycerol are 808,998 and 455,672, respectively.The computation time of each model is 7317 s and 3354 s in a computer with a memory of 128 GB.

Experiment System
To validate the computation results, physical experiments need to be completed.The whole experimental system is shown in Figures 3 and 4. Microwave power is generated from a microwave generator (WSPS-2450-1000M, Wattsine, Chengdu, China) to a waveguide coaxial connector (CAWG-26-N, Euler, Nanjing, China).A circulator connected with water load is employed to protect the source from the reflection power.With a dual directional coupler and a power meter (AV2433, the 41st Institute of CETC, Tsingtao, China) as power monitoring, microwave power is finally fed into the BJ22 waveguide connected with a slide short.Final heating temperature is gained by a fiber optic thermometer (FISO FOT-NS-967A, FISO Technologies, Quebec, QC, Canada).In the simulation, the electromagnetic field distribution is firstly calculated in the frequency domain, and then the dissipated power is calculated.Finally, the temperature rise of the material is updated through the heat transfer equation in the time domain.Since the electrical properties of samples do not vary with temperature, the electromagnetic field distribution will not change and is only calculated once.The whole heating process lasts 120 s, and the computation time of the model depends on the mesh size and the memory of the computer.In this simulation, the mesh sizes of the model with water and model with glycerol are 808,998 and 455,672, respectively.The computation time of each model is 7317 s and 3354 s in a computer with a memory of 128 GB.

Experiment System
To validate the computation results, physical experiments need to be completed.The whole experimental system is shown in Figures 3 and 4. Microwave power is generated from a microwave generator (WSPS-2450-1000M, Wattsine, Chengdu, China) to a waveguide coaxial connector (CAWG-26-N, Euler, Nanjing, China).A circulator connected with water load is employed to protect the source from the reflection power.With a dual directional coupler and a power meter (AV2433, the 41st Institute of CETC, Tsingtao, China) as power monitoring, microwave power is finally fed into the BJ22 waveguide connected with a slide short.Final heating temperature is gained by a fiber optic thermometer (FISO FOT-NS-967A, FISO Technologies, Quebec, QC, Canada).

Experimental Procedures
Physical experiments are carried out with 55 W microwave power at 2.45 GHz.For temperature measurements, several fibers of the fiber optic thermometer have been used and connected to the processing materials through the cut-off waveguide on the BJ22 waveguide.For power measurements, a power meter is employed to verify the incidence and reflection situation of the microwave power.To give more evident temperature rise and make the experiment easier, the length between the quartz cup and the slide short is set as λ e , which is defined as where λ is the wavelength of the microwave in vacuum, f c is the cut-off frequency of the BJ22 waveguide, f is the frequency of the microwave.The corresponding modification of the experimental length of the BJ22 waveguide is realized by the slide short.A simple test system, shown in Figure 5, is performed to adjust the position of the slide short.By combining the S 11 gained from the vector network analyzer (N5230A, Agilent Technologies Inc., Santa Clara, CA, USA), the position of the slide short is confirmed to match the simulation.

Experimental Procedures
Physical experiments are carried out with 55 W microwave power at 2.45 GHz.For temperature measurements, several fibers of the fiber optic thermometer have been used and connected to the processing materials through the cut-off waveguide on the BJ22 waveguide.For power measurements, a power meter is employed to verify the incidence and reflection situation of the microwave power.To give more evident temperature rise and make the experiment easier, the length between the quartz cup and the slide short is set as λe, which is defined as where λ is the wavelength of the microwave in vacuum, fc is the cut-off frequency of the BJ22 waveguide, f is the frequency of the microwave.The corresponding modification of the experimental length of the BJ22 waveguide is realized by the slide short.A simple test system, shown in Figure 5, is performed to adjust the position of the slide short.By combining the S11 gained from the vector network analyzer (N5230A, Agilent Technologies Inc., Santa Clara, CA, USA), the position of the slide short is confirmed to match the simulation.

Experiment Validation
To validate the simulation results, one experiment group with water as the polar solvent while another with glycerol as the non-polar solvent is employed.In the experiments, the silicon carbide spheres are actually placed in a two-layer structure as shown in Figure 6.The same adjustment has been done to the corresponding simulation model.

Experiment Validation
To validate the simulation results, one experiment group with water as the polar solvent while another with glycerol as the non-polar solvent is employed.In the experiments, the silicon carbide spheres are actually placed in a two-layer structure as shown in Figure 6.The same adjustment has been done to the corresponding simulation model.

Experimental Procedures
Physical experiments are carried out with 55 W microwave power at 2.45 GHz.For temperature measurements, several fibers of the fiber optic thermometer have been used and connected to the processing materials through the cut-off waveguide on the BJ22 waveguide.For power measurements, a power meter is employed to verify the incidence and reflection situation of the microwave power.To give more evident temperature rise and make the experiment easier, the length between the quartz cup and the slide short is set as λe, which is defined as where λ is the wavelength of the microwave in vacuum, fc is the cut-off frequency of the BJ22 waveguide, f is the frequency of the microwave.The corresponding modification of the experimental length of the BJ22 waveguide is realized by the slide short.A simple test system, shown in Figure 5, is performed to adjust the position of the slide short.By combining the S11 gained from the vector network analyzer (N5230A, Agilent Technologies Inc., Santa Clara, CA, USA), the position of the slide short is confirmed to match the simulation.

Experiment Validation
To validate the simulation results, one experiment group with water as the polar solvent while another with glycerol as the non-polar solvent is employed.In the experiments, the silicon carbide spheres are actually placed in a two-layer structure as shown in Figure 6.The same adjustment has been done to the corresponding simulation model.Two points of the processing materials have been chosen to perform the temperature variation comparisons between the experiment and the simulation (shown in Figure 7).The center bottom of the quartz cup is selected as the origin of the axes and directions of the axes are shown in Figure 6b.Comparisons show that simulation results match the experiential results well.Two points of the processing materials have been chosen to perform the temperature variation comparisons between the experiment and the simulation (shown in Figure 7).The center bottom of the quartz cup is selected as the origin of the axes and directions of the axes are shown in Figure 6b.Comparisons show that simulation results match the experiential results well.

Effect of Introducing Fluid Materials with Different Dielectric Properties
To give intuitive and convictive results of the influence of the fluid materials on the heating uniformity and safety, a more compact three-layer structure is applied to the simulated silicon carbide spheres and the computations are divided into two situations, where only ε′ of fluid materials is computed and where both ε′ and ε″ are taken into consideration, respectively.Since properties of other related materials stay the same, the other physical properties of the fluid materials have all been defined to be the same with water to avoid their distractions.Meanwhile, the microwave power is set as 100 W and the density of water is set as 1 kg/m 3 in order to get a higher temperature rise, and the initial temperature is set as 20 degrees centigrade.The whole heating process lasts 120 s.In the computation results, the reflection parameters S11 of the heating system, the average body temperature T of the solids, the coefficient of variation (COV) value of the solids' final temperature and the maximum modulus value of the electric field | E |max in the whole processing materials are analyzed.The COV of temperature can be expressed by

Effect of Introducing Fluid Materials with Different Dielectric Properties
To give intuitive and convictive results of the influence of the fluid materials on the heating uniformity and safety, a more compact three-layer structure is applied to the simulated silicon carbide spheres and the computations are divided into two situations, where only ε of fluid materials is computed and where both ε and ε" are taken into consideration, respectively.Since properties of other related materials stay the same, the other physical properties of the fluid materials have all been defined to be the same with water to avoid their distractions.Meanwhile, the microwave power is set as 100 W and the density of water is set as 1 kg/m 3 in order to get a higher temperature rise, and the initial temperature is set as 20 degrees centigrade.The whole heating process lasts 120 s.In the computation results, the reflection parameters S 11 of the heating system, the average body temperature T of the solids, the coefficient of variation (COV) value of the solids' final temperature and the maximum modulus value of the electric field | → E| max in the whole processing materials are analyzed.The COV of temperature can be expressed by Processes 2018, 6, 220 where N is the total number of the temperature points, T i is the temperature at a point, T is the average temperature and T 0 is the initial temperature.Generally speaking, the smaller the COV of the temperature, the smaller the dispersion degree of the temperature distribution, which means the better heating uniformity.
For the first situation, S 11 , T, COV and | → E| max along with ε are shown in Figure 8.It is worthy to note that the reflection parameter is calculated by the rate of reflection power and incident power.Computation results of the system are firstly characterized by the reflection coefficient S 11 , namely the power absorbed by the processing materials.As shown in Figure 8a, the increasing of ε has obvious but nonlinear effects on the S 11 and will thus decide the corresponding T through its impact on the microwave feeding condition.
where N is the total number of the temperature points, i T is the temperature at a point, T is the average temperature and 0 T is the initial temperature.Generally speaking, the smaller the COV of the temperature, the smaller the dispersion degree of the temperature distribution, which means the better heating uniformity.
For the first situation, S11, T , COV and | E |max along with ε′ are shown in Figure 8.It is worthy to note that the reflection parameter is calculated by the rate of reflection power and incident power.Computation results of the system are firstly characterized by the reflection coefficient S11, namely the power absorbed by the processing materials.As shown in Figure 8a, the increasing of ε′ has obvious but nonlinear effects on the S11 and will thus decide the corresponding T through its impact on the microwave feeding condition.
To the heating uniformity of the processing solids, the COV variation shows that it will mainly decrease with the rise of ε′, e.g., it decreases from 0.1931 with ε′ = 1 to 0.0782 with ε′ = 85.However, in small ranges, it should be analyzed more precisely as the COV has fluctuated.Combined with the analyzed T condition, areas where uniform heating with high efficiency could be gained just like the marked areas in Figure 8b.To the heating uniformity of the processing solids, the COV variation shows that it will mainly decrease with the rise of ε , e.g., it decreases from 0.1931 with ε = 1 to 0.0782 with ε = 85.However, in small ranges, it should be analyzed more precisely as the COV has fluctuated.Combined with the analyzed T condition, areas where uniform heating with high efficiency could be gained just like the marked areas in Figure 8b.
For the heating safety of the processing materials, the increase of ε in a certain region (from ε = 1 to ε = 10 as marked in Figure 8c) has a remarkable effect on lowering the | → E| max as it changes from 478,870 V/m to 19,982 V/m, which is reduced by about 24 times.Out of this range, the influence of the increasing ε stays weak.The differences could be further described by the mean square error as it is 191,152.0V/m in the former range while it is 3109.9V/m in the latter, which is reduced by about 60.5 times.Compared with normal heating solids with air surroundings, namely ε = 1, introduction of fluid materials with proper ε shows a more convenient way to achieve uniform and safe microwave heating.
After computations of ε only, the influences of ε" have also been studied.The related computation results shown in Table 2 and Figure 9, and the increasing of ε" from 0 to 0.3 will lower the S 11 and thus bring more efficient heating.While the S 11 and T almost stay the same out of this range.
To the heating uniformity and safety, introducing fluid materials with higher ε" could reduce the | → E| max , as the increasing of ε" from 0.1 to 0.9 has reduced the | → E| max by about 34% to 55%.However, higher ε" will cause worse heating uniformity as shown in Figure 9b.For the heating safety of the processing materials, the increase of ε′ in a region (from ε′ = 1 to ε′ = 10 as marked in Figure 8c) has a remarkable effect on lowering the | E |max as it changes from 478,870 V/m to 19,982 V/m, which is reduced by about 24 times.Out of this range, the influence of the increasing ε′ stays weak.The differences could be further described by the mean square error as it is 191,152.0V/m in the former range while it is 3109.9V/m in the latter, which is reduced by about 60.5 times.Compared with normal heating solids with air surroundings, namely ε′ = 1, the introduction of fluid materials with proper ε′ shows a more convenient way to achieve uniform and safe microwave heating.
After computations of ε′ only, the influences of ε″ have also been studied.The related computation results shown in Table 2 and Figure 9, and the increasing of ε″ from 0 to 0.3 will lower the S11 and thus bring more efficient heating.While the S11 and T almost stay the same out of this range.
To the heating uniformity and safety, introducing fluid materials with higher ε″ could reduce the | E |max, as the increasing of ε″ from 0.1 to 0.9 has reduced the | E |max by about 34% to 55%.However, higher ε″ will cause worse heating uniformity as shown in Figure 9b.As the computation results analyzed above, the introduction of fluid materials with proper relative permittivity does improve the heating uniformity and safety of microwave heating of common stack solid materials.For the real part of relative permittivity of fluid materials, the computed E |max in simulation models is always gained in the interface between the solids and the fluid materials.Combined with the cure trend shown Figure 8b and the relative permittivity of the solids shown in Table 1, it is deduced that the maximum modulus value of the electric field | E |max stays huge when the ε′ of fluid materials is much different with the one of solids.A comparison of the variation of | E |max on fluid materials' ε′ with solids of different relative permittivity is performed in Figure 10 and the simulation results have agreed with our deduction.When dealing with certain solids, introducing fluid materials with higher ε′ can efficiently lower the maximum modulus value of the electric field in order to make a safer heating process.The higher ε′ is also good for improving heating uniformity.On the other hand, the increasing of ε″ in a certain range could obviously improve the heating efficiency and in the whole range its increase will lower the maximum modulus value of the electric field, but bring worse heating uniformity.Since the involved factors are complicated, a proper choice of ε″ to balance the requirements between heating uniformity and safety should be made in the practical

Conclusions
In this work, a general model has been built to study the influence of introducing fluid materials in the heating uniformity and safety of microwave heating solid stack materials.Simulations are computed by COMSOL Multiphysics based on the finite element method.With the mesh independence study and adjustment corresponding to physical experiments, the model has been verified and its simulation results agree well with validation experiments.
The real part of relative permittivity of the introduced fluid materials could significantly lower the maximum modulus value of the electric field in a range of ε′ = 1 to ε′ = 10, but out of this range its influence stays low.The increase of ε′ on the whole studied range will improve the heating uniformity, but in small ranges it has fluctuated.The increase of the imaginary part of relative permittivity from 0 to 0.3 could obviously improve the heating efficiency.On the whole studied range, higher ε″ could lower the maximum modulus value of the electric field and give a safer heating, but it will make the heating uniformity worse.Combined with computational results, the general rules of the influence of fluid materials' relative permittivity on the heating uniformity and safety may be useful for the practical heating process.

Figure 2 .
Figure 2. Normalized power absorption (NPA) variation of heating computations with different mesh sizes.

Figure 2 .
Figure 2. Normalized power absorption (NPA) variation of heating computations with different mesh sizes.

Figure 3 .
Figure 3. Schematic of the laboratory microwave heating system.

Figure 3 .
Figure 3. Schematic of the laboratory microwave heating system.

Figure 3 .
Figure 3. Schematic of the laboratory microwave heating system.

Figure 4 .
Figure 4.The laboratory microwave heating system.

Figure 4 .
Figure 4.The laboratory microwave heating system.

Figure 5 .
Figure 5. Schematic of the simple system for the slide short position confirmation.

Figure 5 .
Figure 5. Schematic of the simple system for the slide short position confirmation.

Figure 4 .
Figure 4.The laboratory microwave heating system.

Figure 5 .
Figure 5. Schematic of the simple system for the slide short position confirmation.

Figure 8 .
Figure 8. Influence of ε′ on the final parameters at time 120 s: (a) Variations of the S11 and the T ; (b) variations of the coefficient of variation (COV) and T ; (c) variations of the | E |max.

Figure 8 .
Figure 8. Influence of ε on the final parameters at time 120 s: (a) Variations of the S 11 and the T; (b) variations of the coefficient of variation (COV) and T; (c) variations of the |

Figure 9 .
Figure 9. Influence of ε″ on the final parameters with ε′ = 85 at time 120 s: (a) Variations of the S11 and the T ; (b) variations of the COV and the | E |max.

Figure 9 .
Figure 9. Influence of ε" on the final parameters with ε = 85 at time 120 s: (a) Variations of the S 11 and the T; (b) variations of the COV and the |

Figure 10 .
Figure 10.The influence of ε′ on the | E |max.

Table 1 .
Summary of material properties applied in the model.

Table 1 .
Summary of material properties applied in the model.

Table 2 .
Comparisons of the COV and |

Table 2 .
Comparisons of the COV and | E |max while ε″ is considered.