# Modeling the Temperature Field of the Microwave Heating Asphalt Mixture with the Utilization of a Mathematical Model

Authors: Yifei Zhang; Qingxian Wang; Bofu Guo
Journal: Applied Bionics and Biomechanics
 (2022)

---

## Abstract

Microwave heating asphalt mixture is an auspicious method since it requires a minimum level of outside intervention. It has been a recent method with some certain advantages over conventional ones. To obtain the temperature field of microwave heating asphalt mixture utilizing an oblique horn antenna, the aperture field of the oblique horn antenna was deduced based on the Maxwell equations. The three-dimensional mathematical model of a microwave heating asphalt mixture was established by introducing a position relation of the horn caliber called an oblique horn antenna, and the calculation results were simplified and verified by experiments. The results showed that the calculated values of the mathematical model tally with the experimental outcomes. The model is found to be accurate and reliable.


---

## Body

## 1. Introduction

With ever-increasing road mileage around the globe, the focus of pavement construction has shifted from new pavement construction to pavement maintenance and rehabilitation [1, 2]. When compared with methods of conventional asphalt pavement heating such as hot air and infrared, microwave heating has the advantages of less pollution, large penetration into the depths, a small temperature gradient, no thermal inertia, and 100% recycling capability of an old asphalt mixture [3, 4]. Thus, microwave heat time (cycles) has a limited influence on asphalt aging within a controlled temperature range [5], with no air emissions or liquid pollutants [6]. Moreover, microwave heating can ensure rapid, uniform, and deep penetration [7].The amount of electricity used by microwave devices is much less than those required to produce a similar effect by electromagnetic induction [8]. Therefore, microwave heating has been a more and more widely employed technique in the maintenance of asphalt roads. The widely benefitted features of this method can be found in [6, 9, 10].The device used for mixing microwaved heated asphalt is mainly composed of a magnetron, a power supply, a horn antenna, and a cooling device. The size of the horn antenna is closely related to microwave frequency. To reduce the volume of microwave heating equipment and improve the microwave power density per unit area of an asphalt road, an oblique horn antenna is utilized in microwave maintenance equipment in practical implementations. The scholars in [11–17] researched the aperture field of the straight horn antenna according to the Maxwell equations. The propagation direction of the microwave is perpendicular to the aperture of the horn, and the microwave power is the largest at the center of the aperture of the horn antenna. Thus, the mathematical models of the temperature field are established based on the above aperture field. While the highest temperature is located directly below the horn antenna, the highest temperature is located below the center of the horn antenna at the same depth. However, it is found that the high-temperature area after heating with the oblique horn antenna is not located at the center of the horn aperture but deviates to the bevel direction of the horn antenna; in other words, the aperture field of the oblique horn antenna is found to be different from that of the straight horn antenna and cannot be replaced by the aperture field of the straight horn antenna. Therefore, more research is needed to be conducted in this direction.


## 2. The Aperture Field of an Oblique Horn Antenna Port

The straight horn antenna and oblique horn antenna are shown in Figure1 [18, 19]. The magnetron is installed on the side of the straight horn antenna. Because the magnetron size is large, it can only be placed after increasing the height of the horn antenna and the size of the bell aperture or increasing the distance between adjacent horn antennas is necessary, which will inevitably lead to a lower power density and slower heating. While the magnetron is installed on the top of the oblique horn antenna with a compact structure and a lower height, the adjacent horn antennas can be contacted. The density of the output power of the horn is larger and so is higher the heating efficiency.Figure 1
The two kinds of the horn antennas. (a) A straight horn antenna. (b) An oblique horn antenna.
(a)(b)In the oblique horn antenna shown in Figure2, a straight horn antenna is constructed by extending the oblique edge M’m and right angle NF to intersect at point O, and the included angle is β. The angle bisector of β and the antenna bottom edge M’n intersect at the O” point, and the vertical line NO’ is made from the point N and extended, which is compared with the oblique edge at the point M. The constructed △OMN is the standard H-plane horn antenna, and the microwave can be considered to be transmitted from point O and propagated along the direction of OO’.Figure 2
An oblique horn antenna.Equation (1) presents the well-known electromagnetic field equations where E⟶ denotes electric field, H⟶ represents magnetic intensity, ∇ notation and ∇x notation denote gradient and curl operations, respectively, j denotes the total current per unit area, ρ denotes the total charge per unit area, and μ and ε denote the spaces related to permittable and permeable.
(1)∇×E⟶=−jωμH⟶,∇×H⟶=jωεE⟶.Then, the cylindrical form,Hy,Hφ,Hρ and Ey,Eφ,Eρ, of this equation is rewritten by
(2)1ρ∂ρEϕ∂ρ−1ρ∂Eρ∂ϕ=−jωμHy1ρ∂Ey∂ϕ−∂Eϕ∂y=−jωμHρ∂Eρ∂y−∂Ey∂ρ=−jωμHϕ,1ρ∂ρHϕ∂ρ−1ρ∂Hρ∂ϕ=jωεEy1ρ∂Hy∂ϕ−∂Hϕ∂y=jωεEρ∂Hρ∂y−∂Hy∂ρ=jωεEϕ.Under the single-mode excitation of the TE10 wave, the expression of the field quantity in the H-plane sector horn can be obtained by resolving the above equation system. Thus,
(3)Ey=CcosnϕHn2kρHϕ=−jkωμ0CcosnϕddkρHn2kρHρ=−jnωμ0ρCsinnϕHn2kρ,where n=π/2φH and k=ωμ0ε0. Hn2kρ is called the second kind of the Hankel function with kρ the number of cases, which represents the outward wave propagating from the top of the horn. When the value kρ grows large, the Hankel function has the following asymptotic formula expressed by
(4)Hn2kρ=2πkρe−jkp−2n+1/4π.Far away from the top of the horn, theHρ is much smaller than the Hφ and can be ignored. Therefore, there exist only two electromagnetic field components, namely, Ey and Hφ. Thus,
(5)Ey=C2πkρcosnϕe−jkρ−n+1/2π/2Hϕ=−Ey120πHρ=Hy=Eϕ=Eρ=0.Suppose that the amplitude of the electric field is cosine distributed along theY direction, according to cosnϕ=cosπ/2ϕ/ϕH, where φH is the half of the opening angle of the H-plane fan-shaped horn. Namely,
(6)ϕH=tan−1D1/2LH≈D12LH,ϕ=tan−1XSLH≈XSLH.Therefore,cosnϕ=cosπXS/D1.Besides, Equation (1) can be expressed in the form of cylindrical coordinates shown in Figure 3.Figure 3
The position relation of the horn caliber.In addition, the phase of each point of the aperture field is different according to the characteristics of surface waves. Figure3 depicts that XS is away from the midpoint O of the aperture at point M, and the phase of the field lags the phase at the center point O, which is caused by the travel difference.
(7)MN=LH2+XS2−LH≈12XS2LH.Therefore, the phase lag is represented by(8)ϕS=k⋅MN=πλ⋅XS2LH.Substituting the expression of the electric field intensity into the H-plane of the horn, we attain(9)Ey=E0cosπXsD1e−jπ/λXs2/LHHx=−Ey120π.For the internal field distribution of the E-plane sector horn, the aperture field of the E-plane horn can be obtained by the same method by analyzing the H-plane horn expressed by(10)Ey=E0cosπXsD1e−jπ/λYs2/LEHx=−Ey120π.Currently, the amplitude distribution law of the aperture field is the same as that of the H-plane horn. Under the excitation of the TE10 wave, it is distributed in cosine form along the X direction and evenly along the Y direction. The phase of the aperture is different, and there is also a phase difference that, however, occurs in the Y direction.The relationship betweenD1 and D is expressed by
(11)D1sinα=Dsinπ−π−β/2,Namely,D1=Dsinα/sinβ/2.The amplitude of the aperture field of the pyramid horn can be considered the same as that of the fan-shaped horn on the H-plane and distributed according to the cosine form. On the E-plane, it can be considered the same as the fan-shaped horn. According to the uniform distribution, the phase of the aperture field is distributedπ/λXS2/LH+YS2/LE. Therefore, the aperture field of the pyramid horn can be expressed by
(12)E=E0cosπXssinβ/2Dsinαe−jπ/λXs2/LH+Ys2/LE,where both LH and LE are the heights of the H-plane and E-plane of the constructed pyramid horn, respectively.


## 3. The Mathematical Model of the Temperature Field for the Microwave Heating

In this research, the discrete element method [20] was adopted. Asphalt is assumed to be heated as an isotropic material. Utilizing the conduction analysis of an infinitesimally small (differential) control volume (DCV) [21], dx⋅dy⋅dz in the asphalt pavement is shown in Figure 4 [22].Figure 4
The energy flow of the microelement.Equation (13) is obtained from the law of the conservation of energy.
(13)Ėin+Ėg−Ėout=Ėst,where Ėin denotes the rate of the energy transfer into a control volume, Ėout denotes the rate of the energy transfer out of a control volume, Ėst represents the rate of the energy change in a control volume, q denotes heat transfer rate, ρ is its density, c is its specific heat capacity, T denotes temperature, t is time, and x, y, and z are rectangular coordinates, respectively. Thus,
(14)Ėg=q̇dxdydz,(15)Ėst=ρc∂T∂tdxdydz,(16)Ėin−Ėout=∂qx∂xdx+∂qy∂ydy+∂qz∂zdz,(17)qx=−k∂T∂xdydz,(18)qy=−k∂T∂ydxdz,(19)qz=−k∂T∂zdxdy.Substituting Equations (14) through (19) simultaneously into Equation (13), it leads to the heat diffusion equation denoted by
(20)k∂2T∂x2+∂2T∂y2+∂2T∂z2+q̇=ρc∂T∂t.In the process of a microwave heating asphalt mixture, if the internal heat source is the microwave power per unit volume, then the expression is defined by(21)∂2T∂x2+∂2T∂y2+∂2T∂z2+P∗,(22)k∂2T∂x2+∂2T∂y2+∂2T∂z2+P∗=ρc∂T∂t.Besides, Equations (21) and (22) denote the temperature field model of a microwave-heated asphalt mixture. Thus, P∗ is represented by Equation (23) and Pd is represented by Equation (24).
(23)P∗=∂Pz∂z=−αP0e−αz,(24)Pd=πfε0ε′tgδE2.E denotes the distribution of the calculated aperture field.


## 4. The Results and Discussion

### 4.1. Simulation

The simulation model of a microwave heating asphalt mixture is established, as shown in Figure5. Firstly, the microwave field is simulated in the Microwave Studio of the CST software. The simulation parameters are set as follows: the size of the horn antenna aperture is 120×90mm, the relative dielectric constant of asphalt mixture is 8.5, the tangent of loss angle is 0.03, the specific heat capacity is 0.9 kJ/kg °C, the thermal conductivity is 0.55 W/k · m, and the density is 2020 kg/m3. Then, the data are imported into multiple physical fields to calculate the temperature field. The simulation parameters are microwave output power of 1.3 kW, the initial temperature of 10°C, and heating time of 14 min.Figure 5
The simulation model.


### 4.2. Experiment

The experimental scheme is defined as follows: the experimental devices are shown in Figure6. To reduce the influence of many temperatures measuring thermocouples in an asphalt mixture on the distribution of the microwave field, five thermocouples were buried at the depths of Z=3 cm, 6 cm, 9 cm, 12 cm, and 15 cm in the center of the material box. The temperature measuring points are shown in Figure 7. The temperatures of the different coordinates were measured by moving the horn antenna. After measuring the temperature field of one point, the asphalt mixture needs to be cooled down to room temperature, and then, a measurement is taken at another point. The coordinates of the temperature measuring points are as follows: (-50,-60), (-50,-30), (-50,0), (-50,30), and (-50,60); (-25,-60), (-25,-30), (-25,0), (-25,30), and (-25,60); (0,-60), (0,-30), (0,0), (0,30), (0,60); (25,-60), (25,-30), (25,0), (25,30), and (25,60); and (50,-60), (50,-30), (50,0), (50,30), and (50,60).Figure 6
A microwave heating unit.Figure 7
The location of the temperature measuring points.Maintaining and constructing asphalt pavements sustainably could be realized by microwave heating which is an auspicious method since it requires a minimum level of outside intervention. Some of the underlined issues are called overheating and nonuniform heating. Thus, both frequency and power are two attributes that are under investigation.


### 4.3. Discussion

Both Figure8 and Table 1 depict the temperature field at the same depth z that is symmetrical about the y=0 planes. Therefore, the experiment can be simplified further. At the same depth, 25 temperature measurement points in the experimental scheme can be simplified to 15 temperature measurement points: (-50,-60), (-50,-30), and (-50,0); (-25,-60), (-25,-30), and (-25,0); (0,-60), (0,-30), and (0,0); (25,-60), (25,-30), and (25,0); and (50,-60), (50,-30), and (50,0).Figure 8
The distribution of the plane temperature whenx=0.Table 1
The results of both simulation and computation.
zySimulation resultsCalculation results306090120150306090120150-6034.931.228.125.523.236.93130.119.120.7-3082.470.657.144.935.480.972.654.54634.60122.9100.976.155.941.7128.3100.974.754.939.93082.470.657.144.935.480.666.852.545.334.26034.931.228.125.523.23629.929.421.520.3As shown in Figure9, in the y=0 and y=30 planes, the simulation values tally with the experimental values highly. In the y=60 plane, the simulation values tally with the experimental values relatively. The experimental area is located at the edge of the horn antenna and belongs to the low-temperature areas heated by the microwave sidelobe. The experimental values are slightly different from the simulation values. The reasons may be (1) the asphalt mixture used in the experiment cannot be completely uniform, (2) the processing error of the horn antenna leads to the change of the microwave field distribution, and (3) the opening surface of the horn antenna is not parallel to the asphalt mixture surface.Figure 9
The plane temperature field whenY is set to different values. (a) The plane temperature field when Y=0. (b) The temperature field when Y=30. (c) The plane temperature field when Y=60.
(a)(b)(c)


## 4.1. Simulation

The simulation model of a microwave heating asphalt mixture is established, as shown in Figure5. Firstly, the microwave field is simulated in the Microwave Studio of the CST software. The simulation parameters are set as follows: the size of the horn antenna aperture is 120×90mm, the relative dielectric constant of asphalt mixture is 8.5, the tangent of loss angle is 0.03, the specific heat capacity is 0.9 kJ/kg °C, the thermal conductivity is 0.55 W/k · m, and the density is 2020 kg/m3. Then, the data are imported into multiple physical fields to calculate the temperature field. The simulation parameters are microwave output power of 1.3 kW, the initial temperature of 10°C, and heating time of 14 min.Figure 5
The simulation model.


## 4.2. Experiment

The experimental scheme is defined as follows: the experimental devices are shown in Figure6. To reduce the influence of many temperatures measuring thermocouples in an asphalt mixture on the distribution of the microwave field, five thermocouples were buried at the depths of Z=3 cm, 6 cm, 9 cm, 12 cm, and 15 cm in the center of the material box. The temperature measuring points are shown in Figure 7. The temperatures of the different coordinates were measured by moving the horn antenna. After measuring the temperature field of one point, the asphalt mixture needs to be cooled down to room temperature, and then, a measurement is taken at another point. The coordinates of the temperature measuring points are as follows: (-50,-60), (-50,-30), (-50,0), (-50,30), and (-50,60); (-25,-60), (-25,-30), (-25,0), (-25,30), and (-25,60); (0,-60), (0,-30), (0,0), (0,30), (0,60); (25,-60), (25,-30), (25,0), (25,30), and (25,60); and (50,-60), (50,-30), (50,0), (50,30), and (50,60).Figure 6
A microwave heating unit.Figure 7
The location of the temperature measuring points.Maintaining and constructing asphalt pavements sustainably could be realized by microwave heating which is an auspicious method since it requires a minimum level of outside intervention. Some of the underlined issues are called overheating and nonuniform heating. Thus, both frequency and power are two attributes that are under investigation.


## 4.3. Discussion

Both Figure8 and Table 1 depict the temperature field at the same depth z that is symmetrical about the y=0 planes. Therefore, the experiment can be simplified further. At the same depth, 25 temperature measurement points in the experimental scheme can be simplified to 15 temperature measurement points: (-50,-60), (-50,-30), and (-50,0); (-25,-60), (-25,-30), and (-25,0); (0,-60), (0,-30), and (0,0); (25,-60), (25,-30), and (25,0); and (50,-60), (50,-30), and (50,0).Figure 8
The distribution of the plane temperature whenx=0.Table 1
The results of both simulation and computation.
zySimulation resultsCalculation results306090120150306090120150-6034.931.228.125.523.236.93130.119.120.7-3082.470.657.144.935.480.972.654.54634.60122.9100.976.155.941.7128.3100.974.754.939.93082.470.657.144.935.480.666.852.545.334.26034.931.228.125.523.23629.929.421.520.3As shown in Figure9, in the y=0 and y=30 planes, the simulation values tally with the experimental values highly. In the y=60 plane, the simulation values tally with the experimental values relatively. The experimental area is located at the edge of the horn antenna and belongs to the low-temperature areas heated by the microwave sidelobe. The experimental values are slightly different from the simulation values. The reasons may be (1) the asphalt mixture used in the experiment cannot be completely uniform, (2) the processing error of the horn antenna leads to the change of the microwave field distribution, and (3) the opening surface of the horn antenna is not parallel to the asphalt mixture surface.Figure 9
The plane temperature field whenY is set to different values. (a) The plane temperature field when Y=0. (b) The temperature field when Y=30. (c) The plane temperature field when Y=60.
(a)(b)(c)


## 5. Conclusion

The aperture field of an oblique horn antenna is different from that of a straight horn antenna. The main reason is that the oblique horn antenna has an asymmetric structure, and there is an included angle between the microwave propagation direction and the aperture, while the microwave propagation direction of the straight horn antenna is vertical and the aperture outward.The distribution of the temperature field of an oblique horn antenna is different from that of a straight horn antenna. At different depths of asphalt mixture, the highest point of the temperature field of the oblique horn antenna is not directly below the center of the horn surface, and it becomes closer and closer to the oblique edge with the increase of depthz. At different depths, the highest point of the temperature field of the straight horn antenna is directly below the center of the horn surface.On the other hand, there is no one-to-one correspondence between the experimental outcomes and the simulation results. The potential difference could stem from those issues as follows: (1) the asphalt mixture used in the experiment cannot be completely uniform, (2) the processing error of the horn antenna leads to the change of the microwave field distribution, and (3) the opening surface of the horn antenna is not parallel to the asphalt mixture surface.The horn antenna has symmetry, and the temperature field also has symmetry too. When measuring the temperature field, the measuring points can be reduced according to the actual situation to improve efficiency.Consequently, the model owns high precision and can be used to guide the design of microwave heating equipment.



---

*Source: 9356893-2022-06-29.xml*