Design of Large-Scale Microwave Cavity for Uniform and Efficient Plastic Heating

To reduce carbon emissions during heating in the manufacturing processes, microwave technology has attracted significant attention. Microwaves have considerable advantages over traditional heating methods, including more rapid heating, lower thermal damage, and eco-friendly processes. To apply microwaves to the manufacturing process, uniform and efficient heating is required. We analyzed the effect of various design parameters for uniform and efficient heating by changing the cavity heights, application of the reflector, and number and positions of waveguides. We conducted a numerical simulation and verified the findings by experiments. The results showed that a slight change in the cavity height altered the electromagnetic field distribution and heating parameters, such as the coefficient of variance and power absorption efficiency. With reflectors installed, 66% of cases exhibited better comprehensive evaluation coefficient (CEC) with consideration of uniform heating and power absorption. The spherical reflector showed 81% of cases, better than those of the ordinary model without a reflector. Furthermore, when double waveguides were installed, the average coefficient of variance (COV) was improved by 22%, and power absorption efficiency was increased by 53% compared to the single waveguide case. When the power applied to the waveguides was doubled, the average COV values improved by 18%. This large-scale analysis will be helpful in applying microwaves to actual industrial sites.


Introduction
Fossil fuels represent approximately 80% of the total energy sources as of 2020 [1]. Besides energy, chemicals produced from fossil fuels are used daily. They are also employed in various fields, such as petroleum, construction, and electricity [2]. However, the processes of manufacturing chemical products and the use of fossil fuels produce emissions of large amounts of carbon dioxide (CO 2 ) (33.5 Gt of CO 2 annually), which is a greenhouse gas. The emission of CO 2 into the atmosphere causes global warming by rapidly increasing the temperature of the earth. There have been active discussions among countries around the world to reduce the use of fossil fuels. To manage CO 2 emissions from the use of fossil fuels, clean renewable energy sources, such as sunlight and wind power, are being used instead of conventional fossil fuels. Moreover, technologies that reduce CO 2 emissions in the process of converting fossil fuels into different energy sources, such as electricity, are being employed. Although clean energy production and various efforts to reduce CO 2 emissions are improving, the consumption of fossil-fuel-based plastic products is rapidly increasing each year, compromising the process of low-carbon transitioning. Therefore, there is a need to reduce the large amount of CO 2 emissions from plastic-producing processes. To address this problem, many countries have implemented carbon reduction policies that mainly Thus, it is necessary to increase heating uniformity and energy efficiency simultaneously in industrial sites.
In this study, the hopper model, a large-capacity equipment model that is applicable to actual industrial sites, was designed. The results of a numerical analysis were compared with the experimental results. Various cavity heights, the application of the reflector, and the number and positions of waveguides were considered. The electromagnetic field distribution and heating effect according to the design parameters were calculated through the simulation to analyze the influence of each parameter. The results of this study can be applied to actual industrial sites and contribute to realizing a sustainable society. To process a material using microwave heating, it is essential to predict the temperature profile of the target material [23]. However, it is difficult to visually examine the temperature profile and electromagnetic field distribution inside a closed metal cavity during an actual test. Therefore, the temperature profile and electromagnetic field distribution were predicted through computer-based finite element analysis. Commercial finite element analysis software (COMSOL Multiphysics ® V5.6, COMSOL Inc., Stockholm, Sweden), which is widely used for microwave heating analysis, was used for the numerical simulation. The simulation was conducted under the frequency-transient condition using the heat transfer, electromagnetic, and transport diluted species modules. Figure 1 shows the relationship and interaction among the modules. The microwaves generated in the electromagnetic module form an electromagnetic field distribution inside the cavity based on Maxwell's equations. The heat transfer module is associated with the electromagnetic module, and the electromagnetic field loss is regarded as a heat source in the heat transfer part. Based on the calculated heat source, the temperature and moisture distribution of the target material were derived from the heat transfer and transport diluted species modules. All these processes occurred according to the energy-mass conservation law.

Simulation and Experiments
high heating efficiency is also the major interest [28,29]. However, focusing only on heating efficiency and not considering uniform heating can lead to a serious result. For example, overheating regions deteriorate the performance of the reactor and the quality of the material [17]. Thus, it is necessary to increase heating uniformity and energy efficiency simultaneously in industrial sites.
In this study, the hopper model, a large-capacity equipment model that is applicable to actual industrial sites, was designed. The results of a numerical analysis were compared with the experimental results. Various cavity heights, the application of the reflector, and the number and positions of waveguides were considered. The electromagnetic field distribution and heating effect according to the design parameters were calculated through the simulation to analyze the influence of each parameter. The results of this study can be applied to actual industrial sites and contribute to realizing a sustainable society.

Simulation Domain and Modules
To process a material using microwave heating, it is essential to predict the temperature profile of the target material [23]. However, it is difficult to visually examine the temperature profile and electromagnetic field distribution inside a closed metal cavity during an actual test. Therefore, the temperature profile and electromagnetic field distribution were predicted through computer-based finite element analysis. Commercial finite element analysis software (COMSOL Multiphysics ®®® V5.6, COMSOL Inc., Stockholm, Sweden), which is widely used for microwave heating analysis, was used for the numerical simulation. The simulation was conducted under the frequency-transient condition using the heat transfer, electromagnetic, and transport diluted species modules. Figure 1 shows the relationship and interaction among the modules. The microwaves generated in the electromagnetic module form an electromagnetic field distribution inside the cavity based on Maxwell's equations. The heat transfer module is associated with the electromagnetic module, and the electromagnetic field loss is regarded as a heat source in the heat transfer part. Based on the calculated heat source, the temperature and moisture distribution of the target material were derived from the heat transfer and transport diluted species modules. All these processes occurred according to the energy-mass conservation law.  with varying positions and several waveguides. The dimensions of the cavity and the design parameters to be applied in the simulation are summarized in Figure 2a,b. The cavity was cylindrical and had a diameter and height of 430 and 461 mm, respectively. The lower hopper had a height of 272 mm, an upper diameter of 430 mm, and a lower diameter of 58 mm. The height of the cavity was adjusted based on the wavelength of the 2.45 GHz microwave. The cavity walls and waveguides were made of copper, which reflects microwaves. In addition, 25 kg of polycarbonate was placed inside the hopper. Table 1 summarizes the material properties of the cavity and the target material required for analysis. Figure 2a shows the cavity model with the hopper structure used in the experiment and analysis. A 2.45 GHz magnetron (frequency spectrum from 2448 to 2458 MHz) was attached to the cavity side face, which operated in the transverse electric 10 mode (TE10). Microwaves with a power of 1 kW were generated by the magnetron and propagated into the cavity through rectangular waveguides. The attached waveguides had dimensions of WR-340 (a width of 86.4 mm and a height of 43.2 mm), and the simulation was performed with varying positions and several waveguides. The dimensions of the cavity and the design parameters to be applied in the simulation are summarized in Figure 2a,b. The cavity was cylindrical and had a diameter and height of 430 and 461 mm, respectively. The lower hopper had a height of 272 mm, an upper diameter of 430 mm, and a lower diameter of 58 mm. The height of the cavity was adjusted based on the wavelength of the 2.45 GHz microwave. The cavity walls and waveguides were made of copper, which reflects microwaves. In addition, 25 kg of polycarbonate was placed inside the hopper. Table 1 summarizes the material properties of the cavity and the target material required for analysis.  The following assumptions were applied to the analysis model to efficiently simulate complex phenomena during the microwave heating [30].

Assumption
The following assumptions were applied to the analysis model to efficiently simulate complex phenomena during the microwave heating [30].

1.
The cavity walls and waveguides are in the impedance condition; 2.
The cavity walls and waveguides are completely insulated; 3.
All materials are homogeneous; 4.
Chemical reactions are not considered in the heating process; 5.
The heat transfer model applies only to polycarbonates; 6.
Shrinkage of pellets is not considered; 7.
The influence of air flow is not considered; 8.
The dielectric constant and loss factor are independent of temperature and humidity; 9.
The frequency of the microwave is fixed at 2.45 GHz, which remains constant. In this simulation, the governing equation of the electric field wave is given by the following equation [31]: where E is the electric field vector, ω is angular frequency, σ is the conductivity, ε 0 8.854 × 10 −12 F/m is permittivity of free space, and k 0 is the wave number of free space, which is defined as [19] where c 0 is the speed of light in vacuum. The relative permeability of µ r can be defined as [32]: where µ is the magnetic constant, that is, a measure of the ability of the dielectric materials to store magnetic energy, µ represents the loss of magnetic field energy. The relative permittivity of ε r can be defined as [33]: where ε is the dielectric constant that represents the ability to store electromagnetic energy, and ε is the loss factor that is known to dissipate the absorbed electromagnetic energy, converting electromagnetic energy into heat. When electromagnetic waves interact with the material, the electromagnetic energy is converted into heat based on Fourier's energy balance equation within the area. The heat transfer in the hopper model can be expressed by Equation (5) [17].
where ρ is the density, C p is the specific heat capacity, u is the velocity field from Darcy's law, and T is the temperature. Heat flux q can be defined as [17]: where k is the thermal conductivity. The electromagnetic loss, Q e , is regarded as the heat source on the right side of Equation (7). Q e can be defined as [22]: where the resistive losses are [34]: and the magnetic losses are [34]: J in Equation (8) is the current density. In Equation (9), B is the magnetic flux density, and H is the magnetic field intensity. The existence of moisture plays an important role in heating. Water evaporation and heat convection can affect the temperature. When the temperature reaches a certain level, the phase and distribution of moisture will change, which can influence the heating effect; hence mass transfer processes should be considered. The mass transfer solved in polycarbonates is expressed as [4] ∇·J i + u·∇c i = R i (10) where c i denotes the concentration of species, u is the fluid velocity, R i is the molar production constant, and J i is the diffusive flux vector, which is given as [4]: where D i is the diffusion coefficient.

Boundary Condition Electromagnetic Boundary Conditions
Boundary conditions of impedance are applied on the cavity wall and waveguides where the electromagnetic field is known to penetrate only a short distance outside the boundary [34]: where E s is the source electric field and n is the normal vector toward the exterior of the target boundary.
In this study, the rectangular waveguides were excited by the frequency of the 2.45 GHz magnetron that operated in the TE10 mode.
The port boundary conditions are expressed as [22]: where β is the propagation constant, v is the frequency of microwave, and v c is the cutoff frequency.

Thermal Condition
The thermal insulation boundary condition means that there is no heat flux, which can be expressed by the following equation [4]:

Mass Condition
Convection is applied to the target and can be expressed as [4] − n·(−D∇c) = J 0 (15) where D is the moisture diffusion coefficient in the sample, c b is the air moisture concentration, k refers to the mass transfer coefficient, and J 0 is the bottom face of the sample and is defined under the no-flux boundary condition.

Mesh
The mesh size is an important element in the simulation as it affects the convergence and accuracy of the analysis. For example, when the mesh size decreases by half, the computation time increases 16 times, and the memory usage increases eightfold [19]. According to the results of a previous study on numerical simulations, at least five mesh elements per wavelength are recommended for proper mesh size in microwave heating [35]. The mesh element quality (MEQ), one of the most commonly used methods, was evaluated by applying five or more mesh elements per wavelength for accurate simulation. In COMSOL Multiphysics, MEQ can be calculated using the following equation [36,37].
where A is the area and h 1 , h 2 , and h 3 are the side lengths of the triangle. 227,146 meshes were used. As shown in Figure 3, the average mesh quality was 0.9, which exceeded the criterion (0.65) of previous studies.

Mesh
The mesh size is an important element in the simulation as it affects the convergence and accuracy of the analysis. For example, when the mesh size decreases by half, the computation time increases 16 times, and the memory usage increases eightfold [19]. According to the results of a previous study on numerical simulations, at least five mesh elements per wavelength are recommended for proper mesh size in microwave heating [35]. The mesh element quality (MEQ), one of the most commonly used methods, was evaluated by applying five or more mesh elements per wavelength for accurate simulation. In COM-SOL Multiphysics, MEQ can be calculated using the following equation [36,37].
where is the area and ℎ , ℎ , and ℎ are the side lengths of the triangle. 227,146 meshes were used. As shown in Figure 3, the average mesh quality was 0.9, which exceeded the criterion (0.65) of previous studies.  Figure 4a shows the experimental setup, and its specifications are the same as those mentioned in the simulation. As shown in Figure 4b, only one 1 kW magnetron was attached to the 221 mm position on the side, and another magnetron was attached to the upper part at a position 120 mm away from the center. The densely filled polycarbonate pellets ensure reliable contact between thermocouples and target materials. The roundended thermocouples were grounded to prevent spark.  Figure 4a shows the experimental setup, and its specifications are the same as those mentioned in the simulation. As shown in Figure 4b, only one 1 kW magnetron was attached to the 221 mm position on the side, and another magnetron was attached to the upper part at a position 120 mm away from the center. The densely filled polycarbonate pellets ensure reliable contact between thermocouples and target materials. The roundended thermocouples were grounded to prevent spark.

Experimental Procedure
Polycarbonate was exposed to a humid environment for several hours to create conditions similar to actual industrial environments. To prevent the difference in heating effect caused by the moisture content difference, the moisture content was maintained at 0.11% by monitoring it using a moisture meter (MS-70, AND, Tokyo, Japan). After inserting 25 kg of polycarbonate pellets into the experimental setup, the magnetron was operated for 600 s. The temperature of the polycarbonate was measured every 2 s using a thermocouple. The environment temperature was maintained 25 °C using Air-conditioning system.

Experimental Procedure
Polycarbonate was exposed to a humid environment for several hours to create conditions similar to actual industrial environments. To prevent the difference in heating effect caused by the moisture content difference, the moisture content was maintained at 0.11% by monitoring it using a moisture meter (MS-70, AND, Tokyo, Japan). After inserting 25 kg of polycarbonate pellets into the experimental setup, the magnetron was operated for 600 s. The temperature of the polycarbonate was measured every 2 s using a thermocouple. The environment temperature was maintained 25 • C using Air-conditioning system.

Experiment Validation
To verify the validity of the simulation model, the results of the experiment described above were compared with the simulation results. The atmospheric temperature at the time of the experiment, 25 • C, was set as the initial value, and the temperature values at the same positions as in the experimental setup were collected while 1 kW and 2 kW microwaves were radiated inside the cavity for 600 s. As shown in Figure 5, the root mean square error between the experimental values and simulation results were 0.791, 2.565 and 1.982 • C at the measuring depth of 17, 106 and 197 mm, respectively, when heating was performed using 1 kW for 10 min and 3.224, 6.408, 6.971 • C at the same points when 2 kW was supplied by installing two waveguides. The small errors between the simulation and experimental results confirmed that the numerical analysis model is reliable.

Parameters
Microwaves are generated from a magnetron and transmitted to the microwave cavity through the waveguide. In this instance, various design parameters, such as the cavity height, geometry of the reflector inside the cavity, positions and number of microwave port that reflect the microwaves in different directions, create various electromagnetic field distributions. The created electromagnetic field converts the energy of the microwaves into thermal energy through interactions with objects, thereby causing different heating results. As mentioned in the Introduction, reducing hot spots through uniform heating is an important element for microwave-based heating. The following coefficient of variance (COV) has been widely used to quantify the heating uniformity after microwave irradiation [38]: where Ti is the temperature at a specific point, Ta is the average temperature of the target material to be heated, and T0 is the initial temperature. A lower COV value implies that the target material is heated more uniformly. The power absorption efficiency (PAE, ), which is a value indicating the energy absorption efficiency, has also been widely em-

Parameters
Microwaves are generated from a magnetron and transmitted to the microwave cavity through the waveguide. In this instance, various design parameters, such as the cavity height, geometry of the reflector inside the cavity, positions and number of microwave port that reflect the microwaves in different directions, create various electromagnetic field distributions. The created electromagnetic field converts the energy of the microwaves into thermal energy through interactions with objects, thereby causing different heating results. As mentioned in the Introduction, reducing hot spots through uniform heating is an important element for microwave-based heating. The following coefficient of variance (COV) has been widely used to quantify the heating uniformity after microwave irradiation [38]: where T i is the temperature at a specific point, T a is the average temperature of the target material to be heated, and T 0 is the initial temperature. A lower COV value implies that the target material is heated more uniformly. The power absorption efficiency (PAE, η), which is a value indicating the energy absorption efficiency, has also been widely employed together with COV to understand the heating phenomenon [39].
Here, P i is the total energy supplied to the model, and P a is the absorbed power, which is calculated through the volumetric integration of the power loss density. A higher PAE value indicates that the supplied energy is more efficiently converted into heat. In this study, the positions and sizes of hot spots that occurred inside the cavity were confirmed through an electromagnetic field distribution, and the heating uniformity and efficiency were analyzed for various parameters using the COV and PAE values.

Effect of Cavity Height
Since the distribution of the electromagnetic field changes considerably depending on the geometry of the cavity, the analysis was conducted according to the cavity height first. The cavity height was increased by 12.2 mm, which corresponded to a tenth of the wavelength of the 2.45 GHz microwave from 461 mm. The total increase was 122 mm. Figure 6 shows the electromagnetic field distribution according to height. Finding the pattern of the electromagnetic field distribution according to height using the electromagnetic field distribution map alone was challenging, as was predicting the distribution of hot spots and cold spots. For example, in the case of the cavity whose height was increased by 10% of the wavelength, cold spots were formed in the central part, whereas hot spots were formed in parts close to the surface. When the cavity height was increased by 50% of the wavelength, cold spots were formed toward the surface and wide hot spots were formed in the central part. This showed significant changes in the electromagnetic field despite a small change in geometry. Moreover, it showed that it is almost impossible to intuitively predict the electromagnetic field distribution.
To examine the heating effect according to height, the COV and PAE are summarized in Figure 7. When the height was increased by 20% of the wavelength, the largest COV value of 0.852 was observed, and the uniformity value became approximately two times higher than that when the height was increased by 80%, of the wavelength and the lowest COV value was observed. The heating efficiency also significantly changed depending on height, similar to the heating uniformity. When the height was increased by 50% of the wavelength, 141 W of overall power was absorbed. This occurred when the PAE value was 0.141 (14.1%), and the absorption rate was approximately 11 times higher than when the height was increased by 60% of the wavelength and the lowest PAE value was observed. A comparison between the PAE value and the electromagnetic field distribution summarized in Figure 6 reveals that the PAE value was also high when a strong electromagnetic field was generated. The positions and number of hot spots in the target material changed depending on the microwave cavity height owing to changes in the electromagnetic field distribution and the amount of electromagnetic energy passing through the polycarbonate. This resulted in differences in the temperature distribution and energy absorption efficiency. It indicated that the cavity height is a design parameter that significantly affects the heating effect and uniformity.
by 10% of the wavelength, cold spots were formed in the central part, whereas hot spots were formed in parts close to the surface. When the cavity height was increased by 50% of the wavelength, cold spots were formed toward the surface and wide hot spots were formed in the central part. This showed significant changes in the electromagnetic field despite a small change in geometry. Moreover, it showed that it is almost impossible to intuitively predict the electromagnetic field distribution. To examine the heating effect according to height, the COV and PAE are summarized in Figure 7. When the height was increased by 20% of the wavelength, the largest COV value of 0.852 was observed, and the uniformity value became approximately two times higher than that when the height was increased by 80%, of the wavelength and the lowest COV value was observed. The heating efficiency also significantly changed depending on height, similar to the heating uniformity. When the height was increased by 50% of the wavelength, 141 W of overall power was absorbed. This occurred when the PAE value was 0.141 (14.1%), and the absorption rate was approximately 11 times higher than when the height was increased by 60% of the wavelength and the lowest PAE value was observed. A comparison between the PAE value and the electromagnetic field distribution summarized in Figure 6 reveals that the PAE value was also high when a strong electromagnetic field was generated. The positions and number of hot spots in the target material changed depending on the microwave cavity height owing to changes in the electromagnetic field distribution and the amount of electromagnetic energy passing through the polycarbonate. This resulted in differences in the temperature distribution and energy absorption efficiency. It indicated that the cavity height is a design parameter that significantly affects the heating effect and uniformity.

Effect of a Reflector
When an object is irradiated with microwaves, three reactions occur; that is, some waves are reflected, some pass through, and some are absorbed. Transparency and absorption are intrinsic characteristics of the material and vary depending on the wavelength of the irradiated microwaves. Reflection could be related to the target material, although it mainly occurs in the cavity because the reflectivity of the target material is relatively low. Thus, the heating pattern can be adjusted by changing the geometry of the cavity because it changes the energy distribution of microwaves. In this regard, research has been widely conducted to achieve more uniform heating than before the use of inter-

Effect of a Reflector
When an object is irradiated with microwaves, three reactions occur; that is, some waves are reflected, some pass through, and some are absorbed. Transparency and absorption are intrinsic characteristics of the material and vary depending on the wavelength of the irradiated microwaves. Reflection could be related to the target material, although it mainly occurs in the cavity because the reflectivity of the target material is relatively low.
Thus, the heating pattern can be adjusted by changing the geometry of the cavity because it changes the energy distribution of microwaves. In this regard, research has been widely conducted to achieve more uniform heating than before the use of internal modification by adjusting the reflection of electromagnetic waves through a change in the cavity geometry or by installing objects, such as a rotating body or a conveyor belt, inside the cavity [40]. In this study, spherical, cylindrical, and conical reflectors that reflect microwaves were installed inside the microwave cavity, and changes in the electromagnetic field and the heating effect were analyzed. To evaluate the performance of the cavity in various directions, the comprehensive evaluation coefficient (CEC), which considers the heating efficiency and heating uniformity, as well as the PAE and COV, is defined as follows: A higher CEC value indicates a higher energy absorption rate and more uniform heating. Therefore, CEC can be used as a comprehensive parameter [41].

Effect of Reflector Shape
Conical, spherical, and cylindrical reflectors were placed inside the cavity and on the opposite side of the magnetron. Their dimensions and geometry are summarized in Figure 2.
To analyze the influence of the installed reflectors on the internal environment and heating effect, a simulation was performed by applying three types of reflectors and changing their cavity heights. As can be observed from the electromagnetic field distribution in Figure 8, more notable changes in the electromagnetic field were observed when the reflectors were installed than when only the height of the cavity was changed. The highest intensity of the electromagnetic field with no height change was 4.82 × 10 4 V/m for the existing model; however, it increased by approximately 60% to 7.84 × 10 4 V/m when a conical reflector was installed. The spherical and cylindrical reflectors also increased the highest intensity of the electromagnetic field to 4.88 × 10 4 and 5.82 × 10 4 V/m, respectively. To analyze the influence of the installed reflectors on the internal environment and heating effect, a simulation was performed by applying three types of reflectors and changing their cavity heights. As can be observed from the electromagnetic field distribution in Figure 8, more notable changes in the electromagnetic field were observed when the reflectors were installed than when only the height of the cavity was changed. The highest intensity of the electromagnetic field with no height change was 4.82 × 10 4 V/m for the existing model; however, it increased by approximately 60% to 7.84 × 10 4 V/m when a conical reflector was installed. The spherical and cylindrical reflectors also increased the highest intensity of the electromagnetic field to 4.88 × 10 4 and 5.82 × 10 4 V/m, respectively. The simulation was performed for 44 cases while varying the reflector geometry (4 cases) and cavity height (11 cases); the COV and PAE values are summarized in Table 2. Whether more uniform heating occurs in the presence of reflectors was analyzed based on the case with no reflector. Among the 33 cases with reflectors, lower COV values were observed in the presence of reflectors in 19 cases (57%), confirming more uniform heating than in the models with no reflector. In addition, the PAE value was higher in the presence of reflectors in 21 cases (63%) out of the 33 cases, showing that the installation of reflectors can increase the energy absorption. In particular, for the cavity whose height increased by 20% of the wavelength compared to the ordinary model, the energy absorption was approximately 3% when no reflector was installed, while it increased at least twofold when any of the reflectors was installed.  The simulation was performed for 44 cases while varying the reflector geometry (4 cases) and cavity height (11 cases); the COV and PAE values are summarized in Table 2. Whether more uniform heating occurs in the presence of reflectors was analyzed based on the case with no reflector. Among the 33 cases with reflectors, lower COV values were observed in the presence of reflectors in 19 cases (57%), confirming more uniform heating than in the models with no reflector. In addition, the PAE value was higher in the presence of reflectors in 21 cases (63%) out of the 33 cases, showing that the installation of reflectors can increase the energy absorption. In particular, for the cavity whose height increased by 20% of the wavelength compared to the ordinary model, the energy absorption was approximately 3% when no reflector was installed, while it increased at least twofold when any of the reflectors was installed.  Figure 9 shows the CEC values for 44 cases. For the CEC value that reflected both the energy efficiency and heating uniformity, 66% (filled marks in Figure 9) exhibited better values than the ordinary model with no reflector. When the reflectors were compared, the spherical and conical reflectors exhibited better performance than the cylindrical reflector. When only the spherical and conical reflectors were considered, 16 cases out of 22 (72%) exhibited higher CEC values than the model with no reflectors. In the case of the spherical reflector, 9 cases out of the 11 (81%) showed better performance. This appears to be because the spherical and conical reflectors reflected electromagnetic waves in three dimensions, whereas the cylindrical reflector reflected them in two dimensions. The installation of reflectors was found to be an important factor for cavity design as it changed the electromagnetic field distribution more notably even though it did not always lead to more uniform and higher heating efficiency.

Effect of Double Waveguides
In the preceding chapters, changes in the electromagnetic field and heat distribution due to changes in cavity geometry were analyzed. Changing the cavity geometry and the method of using multiple waveguides have also been proposed as measures to achieve uniform and high-efficiency heating [23]. However, most studies have focused on the application of multiple waveguides to the fixed cavity. However, studies on the influence of detailed cavity design parameters, such as the position, port direction, and power, on the electromagnetic field distribution and heating effect are insufficient and remain on a laboratory scale. Therefore, more waveguides were installed in the cavity of various geometries, and the temperature uniformity and energy efficiency of the target material during heating were analyzed.
When only the spherical and conical reflectors were considered, 16 cases out of 22 (72%) exhibited higher CEC values than the model with no reflectors. In the case of the spherical reflector, 9 cases out of the 11 (81%) showed better performance. This appears to be because the spherical and conical reflectors reflected electromagnetic waves in three dimensions, whereas the cylindrical reflector reflected them in two dimensions. The installation of reflectors was found to be an important factor for cavity design as it changed the electromagnetic field distribution more notably even though it did not always lead to more uniform and higher heating efficiency.

Effect of Double Waveguides, Position, and Direction
To analyze the influence of the double waveguides, another waveguide was installed in the upper part in addition to the existing waveguide position. The waveguides were vertically arranged on the side and in the upper part of the cavity to reduce the resonance and offset between the electromagnetic waves from them. Figure 10 summarizes the geometry and installation positions of the waveguides. In the case of the waveguide placed in the upper part, eight cases were simulated based on four positions and two directions. The two waveguides had dimensions of WR-340, and each waveguide had a power of 0.5 kW. Accordingly, a total of 1 kW of energy could be emitted as in the single waveguide model.

Effect of Double Waveguides
In the preceding chapters, changes in the electromagnetic field and heat distribution due to changes in cavity geometry were analyzed. Changing the cavity geometry and the method of using multiple waveguides have also been proposed as measures to achieve uniform and high-efficiency heating [23]. However, most studies have focused on the application of multiple waveguides to the fixed cavity. However, studies on the influence of detailed cavity design parameters, such as the position, port direction, and power, on the electromagnetic field distribution and heating effect are insufficient and remain on a laboratory scale. Therefore, more waveguides were installed in the cavity of various geometries, and the temperature uniformity and energy efficiency of the target material during heating were analyzed.

Effect of Double Waveguides, Position, and Direction
To analyze the influence of the double waveguides, another waveguide was installed in the upper part in addition to the existing waveguide position. The waveguides were vertically arranged on the side and in the upper part of the cavity to reduce the resonance and offset between the electromagnetic waves from them. Figure 10 summarizes the geometry and installation positions of the waveguides. In the case of the waveguide placed in the upper part, eight cases were simulated based on four positions and two directions. The two waveguides had dimensions of WR-340, and each waveguide had a power of 0.5 kW. Accordingly, a total of 1 kW of energy could be emitted as in the single waveguide model. As shown in the surface electromagnetic field distribution in Figure 11, a hot spot occurred at the top center of the target material when double waveguides were installed, and it was wider compared with the existing single waveguide model. When double waveguides were installed, the intensity of the electromagnetic field was found to be higher than the maximum intensity of the electromagnetic field in the existing model (2.54 × 10 4 V/m), irrespective of the parameters. As shown in the surface electromagnetic field distribution in Figure 11, a hot spot occurred at the top center of the target material when double waveguides were installed, and it was wider compared with the existing single waveguide model. When double waveguides were installed, the intensity of the electromagnetic field was found to be higher than the maximum intensity of the electromagnetic field in the existing model (2.54 × 10 4 V/m), irrespective of the parameters.  Table 3 shows the COV and PAE for the eight positions and directions of the waveguides. First, in the case of Direction 1, the most uniform heating was observed when the waveguide was located at position D in Figure 11, where the COV value was found to be 0.659. All four positions were more uniformly heated than when the single waveguide was used, and the average COV value was 0.710, exhibiting 22% higher uniformity than when one waveguide port was used. The highest energy absorption was also observed at position D, similar to the COV, with the PAE value being 0.116. The PAE value was also higher at all positions than when the single waveguide was used. The average value of the four positions was 0.08, which was 22% higher than that of the ordinary model.
In the case of Direction 2, when the waveguide was located at position D, the most uniform heating was achieved among the eight cases, as the COV value was 0.566. The average COV value for the four positions was 0.701, which was higher than that when the single waveguide was installed. The average PAE value for the four positions was also 0.059, which was excellent compared to the case when the single waveguide was installed.
Consequently, when double waveguides were used, more excellent heating uniformity and energy absorption were observed in all eight cases. The use of the double waveguides made it possible to reach areas that were difficult to reach when microwaves were radiated from the side; by radiating microwaves in the vertical direction as well, more uniform heating occurred. The target material was uniformly heated, and the energy absorption increased compared to when the single waveguide was used because the cases of generating hot spots due to the concentration of electromagnetic waves caused by interference and cold spots due to the offset between the microwaves decreased.   Table 3 shows the COV and PAE for the eight positions and directions of the waveguides. First, in the case of Direction 1, the most uniform heating was observed when the waveguide was located at position D in Figure 11, where the COV value was found to be 0.659. All four positions were more uniformly heated than when the single waveguide was used, and the average COV value was 0.710, exhibiting 22% higher uniformity than when one waveguide port was used. The highest energy absorption was also observed at position D, similar to the COV, with the PAE value being 0.116. The PAE value was also higher at all positions than when the single waveguide was used. The average value of the four positions was 0.08, which was 22% higher than that of the ordinary model.
In the case of Direction 2, when the waveguide was located at position D, the most uniform heating was achieved among the eight cases, as the COV value was 0.566. The average COV value for the four positions was 0.701, which was higher than that when the single waveguide was installed. The average PAE value for the four positions was also 0.059, which was excellent compared to the case when the single waveguide was installed. Consequently, when double waveguides were used, more excellent heating uniformity and energy absorption were observed in all eight cases. The use of the double waveguides made it possible to reach areas that were difficult to reach when microwaves were radiated from the side; by radiating microwaves in the vertical direction as well, more uniform heating occurred. The target material was uniformly heated, and the energy absorption increased compared to when the single waveguide was used because the cases of generating hot spots due to the concentration of electromagnetic waves caused by interference and cold spots due to the offset between the microwaves decreased.

Effect of Double Waveguide Power
In the preceding section, the effect of double waveguides was confirmed. For their application in industries, the heating effect according to supplied power is also an important element. Therefore, in this section, the heating effect was compared and analyzed while the power per waveguide was increased from 0.5 to 1 kW. The electromagnetic field distribution map in Figure 12 shows that when the power of the waveguides was increased, the distribution of the electromagnetic field remained similar, whereas its intensity increased. Thus, the position of the hot spot also remained similar, although its size increased.

Effect of Double Waveguide Power
In the preceding section, the effect of double waveguides was confirmed. For their application in industries, the heating effect according to supplied power is also an important element. Therefore, in this section, the heating effect was compared and analyzed while the power per waveguide was increased from 0.5 to 1 kW. The electromagnetic field distribution map in Figure 12 shows that when the power of the waveguides was increased, the distribution of the electromagnetic field remained similar, whereas its intensity increased. Thus, the position of the hot spot also remained similar, although its size increased. When the COV and PAE values described in Table 4 were calculated, the COV value was found to decrease regardless of the position and direction, as the power per port increased from 0.5 to 1 kW. In particular, when the power was 1 kW in Direction 1, the average COV value for the four positions was 0.647, which was 117% lower than the case when the power was 0.5 kW. The energy efficiency remained similar as the input power doubled even though the energy absorption almost doubled.
Since the inside of the cavity was insulated, the average temperature also doubled as the energy supply doubled. As the positions of the waveguides remained unchanged, the temperature profile was maintained and only its intensity increased. The increased average temperature led to a relatively uniform temperature distribution because the absorbed heat could not be transferred to the outside of the cavity. This indicated that the boundary conditions, as well as the cavity geometry, are important design parameters for the cavity. The excellent insulation effect of the cavity can save energy and create a uniform temperature distribution by preventing the loss of thermal energy. When the COV and PAE values described in Table 4 were calculated, the COV value was found to decrease regardless of the position and direction, as the power per port increased from 0.5 to 1 kW. In particular, when the power was 1 kW in Direction 1, the average COV value for the four positions was 0.647, which was 117% lower than the case when the power was 0.5 kW. The energy efficiency remained similar as the input power doubled even though the energy absorption almost doubled. Since the inside of the cavity was insulated, the average temperature also doubled as the energy supply doubled. As the positions of the waveguides remained unchanged, the temperature profile was maintained and only its intensity increased. The increased average temperature led to a relatively uniform temperature distribution because the absorbed heat could not be transferred to the outside of the cavity. This indicated that the boundary conditions, as well as the cavity geometry, are important design parameters for the cavity. The excellent insulation effect of the cavity can save energy and create a uniform temperature distribution by preventing the loss of thermal energy.

Conclusions
In this study, the microwave heating effect was analyzed on a pilot-scale according to the design parameters, such as the cavity height, installation and type of reflectors, and number and position of waveguides. The numerical analysis results were verified by comparison with the experimental temperature profile, and the following conclusions were drawn: (1) Despite a slight change in the cavity height caused by changing the wavelength by one-tenth, the electromagnetic field distribution and parameters were considerably changed. The COV value, which means heating uniformity, and the PAE value, which means heating efficiency, were changed up to double and 11 times, respectively; (2) The installation of reflectors significantly changed the electromagnetic field distribution. When reflectors were installed, uniform heating was generally achieved, and the power absorption also increased. The spherical reflector had the highest efficiency; 9 out of 11 cases showed better CEC values with a comprehensive coefficient that reflects both PAE and COV. (3) When double waveguides were used, the target material was uniformly heated, and the energy absorption increased at the same time because electromagnetic waves reached areas that were difficult to reach using the single waveguide. So, the average COV was improved by 22%, and power absorption efficiency was increased by 53% compared to the single waveguide model. These results came from the reduction of interference phenomenon by the vertical arrangement of waveguides. (4) As the applied power increased, the form of the electromagnetic field distribution did not change, although the size of the hot spot increased. The average COV value was 0.647, which was 117% lower than the basic model.
The heat inside the cavity increased owing to the increased power supply, and the increase in the average temperature of the target material led to uniform heating. The above observations revealed that the cavity geometry, application of reflectors, arrangement of waveguides, and insulation conditions of the cavity are important parameters for cavity design. This study provides the way microwave heating applies to actual industrial sites. Based on the uniform and efficient heating shown in this study, considering the various factors such as chemical reactions, phase change, etc., in the future, it is expected to be applicable to processes such as drying, synthesis and decomposition which are widely used in industry.