Degree of Coupling in Microwave-Heating Polar-Molecule Reactions

Microwave-assisted chemical reactions have been widely used, but the overheating effect limits further applications. The aim of this paper is to investigate the coupling degree of the electromagnetic field and thermal field in microwave-heating chemical reactions whose polarization changes as the reactions proceed. First, the entropy-balance equation of microwave-heating polar-molecule reactions is obtained. Then, the coupling degree of the electromagnetic field and the thermal field in microwave-heating polar-molecule reactions is derived, according to the entropy-balance equation. Finally, the effects of reaction processes on the degree of coupling are discussed. When the time scale of the component-concentration variation is much greater than the wave period during the chemical processes, the degree of coupling is sufficiently small, and the electric and thermal fields are considered as weakly coupled. On the other hand, the degree of coupling may change during the reactions. If the absolute value of the coupling degree becomes larger, due to the change in component concentration, this will lead to a transition from weak coupling to strong coupling.


Introduction
Microwave heating has been used increasingly in the chemical industry. Compared with traditional heating, microwave heating has many advantages, such as acceleration of reaction rate, reduction in energy consumption and improvement of physical and mechanical properties of products [1][2][3]. Unfortunately, problems associated with the scaling-up of laboratory units to industrial capacities, such as hot spots and overheating effect, have limited the application of microwaves in the chemical industry [4,5]. To overcome these problems, it is essential to investigate the thermal analysis of microwaveheating chemical reactions, based on the multiphysics calculation.
Solving a fully time-transient coupled electromagnetic-thermal problem is very expensive, as the solution to the final set of equations is required at each time-step [6]. In the majority of industrial-heating problems, the electromagnetic dynamics are much faster than the thermal dynamics. On the time scale, the electromagnetic field is weakly coupled with a time-transient non-linear thermal field, and a weakly coupled electromagnetic-thermal model is adequate [7,8]. Therefore, in terms of the multiphysics calculation, it is possible to implement an iterative method based on a weak-coupling procedure. However, in a reaction process, certain components in the chemical reactions may be sensitive to the thermodynamic state, and chemical reactions may have a significant influence on the thermodynamic state [9]. Meanwhile, most chemical reactions are non-stationary media, which exhibit the time-varying characteristics because the component concentrations change with the reaction time [10]. Consequently, the degree of physical interaction between the electromagnetic field and the thermal field may be changed profoundly during the reaction processes. More importantly, the degree of coupling may affect the accuracy of the solution of an iterative method [11]. Unfortunately, little attention to detail and implications has been given to the multiphysics coupling between the electromagnetic field and the thermal field in microwave-heating chemical reactions.
In this paper, we investigate the degree of coupling in microwave-heating polarmolecule reactions. In the next section, the entropy-balance equation for microwave heating a polarized system is reviewed and the entropy-balance equation for microwave heating reactions is derived according to the energy-balance equation. Then, based on the entropy-source strength, the coupling degree of the electromagnetic field and thermal field can be given in terms of the phenomenological equations. Finally, we discuss the effects of the component-concentration variation on the coupling degree of the electromagnetic field and the thermal field in microwave-heating reactions.

In the Polarized System
For a non-reacting multicomponent-polarized-system with n components inside the microwave fields, the balance equation of internal energy can be described by the first law of thermodynamic, as [12,13]: where u is the specific internal energy, q is the heat added to a mass element of unit mass, p is the equilibrium pressure, v is the specific volume, and ρ is the mass density of the material; µ j and c j are the chemical potential and the weight fraction of the jth component, respectively. E and B, respectively, are the instantaneous electric and magnetic fields inside the system, which induce the local electric polarization, P, and the local magnetic polarization, M. Although microwave heating is not at global equilibrium, it is assumed that there exists a state of local equilibrium in any small-mass element when the deviations from the equilibrium are not too large. For the hypothesis of local equilibrium, the local entropy is the same function as dq = Tds in equilibrium, in which s is the entropy per unit mass and T is the absolute temperature. E eq and B eq represent the instantaneous values of the electric and magnetic field at local equilibrium, respectively. Equation (1) can be expressed as Introducing Equation (2) into Equation (1), the entropy-balance equation in the polarized system can be expressed as

In Polar-Molecule Reactions
We consider an isotropic mixture consisting of n components, amongst which m chemical reactions are possible. At constant pressure and temperature, the heat conduction and diffusion caused by uniform pressure or temperature may be neglected. Furthermore, we restrict ourselves to the diffusion process where the jth component concentration is non-uniform over the system. Under these conditions, based on the conservation of mass in chemical reactions and the diffusion process, the rate of the mass change of the jth component can be expressed as [14] ρ J j is the diffusion flux of the jth component; ν jr J r is the production of the jth component per unit volume in the rth chemical reaction, where J r is the chemical-reaction rate of the rth reaction. In addition, ν jr is proportional to the stoichiometric coefficient with which the jth component appears in the rth chemical reaction. For chemical reactions and the diffusion process, with the addition of Equation (4), the rate of entropy change due to chemical reactions is then given as where the chemical affinity of the rth reaction, a r , is defined as [15,16]. Finally, by merging Equations (5) and (3), in microwave-heating chemical reactions, the explicit form of the entropy-balance equation can be expressed as

The Heat Flow and Polarization Current of Microwave Heating in Polar-Molecule Reactions
The differential form of the energy-conservation principle defines the relationship between the heat added per unit of mass, q, and the heat flow, J q , as Inserting Equation (7) into Equation (6), in microwave-heating chemical reactions, we get As the magnetic permeability of chemical reactions is close to that of free space, the magnetic polarization M is not considered. Consequently, we obtain ρ ∂s ∂t According to the second law of thermodynamics [14], where J s is the entropy flux per unit area and unit time, and Θ is the entropy-source strength.
Comparing Equation (9) with Equation (10), the expressions for the entropy flux and the entropy production are given as and Furthermore, if chemical reactions have uniform density and the diffusion flux J j is neglected, Equation (12) reduces to To further investigate the influence of the reaction process on the heat and polarization without loss of generality, we assume that chemical reactions are homogeneous isotropic media without free charges and currents, and only the jth component has a polarization P = φ j * [c j E]. φ j and c j are the rotational-diffusion function and concentration of the jth component, respectively; * represents the convolution; φ j = ε 0 χ exp(−t/τ)/τ, where χ is the electric susceptibility, ε 0 is the vacuum permittivity and τ is the dielectric relaxation time [10]. Accordingly, the differential form of polarization is ∂P/∂t = −(P − ε 0 χc j E)/τ, where ∂P/∂t represents the polarization current, J p , and we can get E eq = P/(ε 0 χc j ). Equation (13) becomes When uniform temperature is attained (∇T = 0) and the reactions are at equilibrium (the chemical affinity, a r , vanishes), we have Specifically, the term τ ε 0 χc j ∂P ∂t 2 on the right-hand side denotes the power converted into local entropy production by electro-thermal coupling, which corresponds to the transient power-loss-density of electromagnetic waves in the reactions reported by Liao et al. [17].
Based on the phenomenological equations and Curie-Prigogine principle, Equation (14) leads to the related phenomenological equations and Here, the phenomenological coefficients L qq , L qp , L pq , L pp and l rh , which are scalar quantities in isotropic media, represent the interactions between the various processes. Equations (16) and (17)  After comparing Equations (16) and (17) with Fourier's law J q = −k∇T and the differential form of the polarization expression, where k is the thermal conductivity, we get L qq = kT 2 (19) and Then, the coefficients L pq and L qp are expressed by [18] L pq = ∂J p and The cross-phenomenological coefficient represents the coupling of the heat flow and polarization current in isotopic media; namely, the polarization current arises from the thermal diffusion. According to the Onsager reciprocal relation, the coefficient L qp is equal to L pq , that is L qp = L pq .
The final expressions for the heat flow and polarization current are then expressed as and Equation (23) indicates that the heat flow is determined by both thermal conduction and polarization, and quantitatively describes a particular excess effect due to polarization loss. In fact, in the absence of a thermal gradient, that is ∇T = 0, the heat flow arises from the power loss of the electromagnetic waves. On the other hand, Equation (24) suggests that the polarization current is influenced not only by the instantaneous polarization and electric field, but also by the thermal gradient. Here, we suppose that the instantaneous electric field, E, is zero, and therefore P = 0. Thus, Equation (24) implies that the polarization current can flow in non-uniformly heated reactions even if the electric field, E, is zero. This can be explained by the directional motion of the polar molecules due to the thermal gradient, namely the thermal diffusion.

The Coupling Degree of the Electromagnetic Field and Thermal Field
Kedem and Caplan have defined the degree of coupling, λ qp , as which has been used as a basis for comparison of systems with various coupled flows [19,20]. By using Equations (19), (20) and (22), the degree of coupling can be expressed as Equation (26) shows that λ qp is directly proportional to the rate of J q variation, and inversely proportional to the rate of J p variation.
The reactions are proofed to exhibit the time-varying and dispersive characteristics. The time-varying characteristics originate from the changes in component concentration.
When the time scale of the component-concentration variation is greater than the wave period, adiabatic approximations can be used to neglect transients [10]. The reactions can be considered as general dispersive media, in which on the time scale, the rate of J q variation is much smaller than the rate of change in J p . For this situation, based on Equation (26), λ qp is usually sufficiently small. Therefore, the influence of the reaction processes on the degree of coupling can be ignored. Furthermore, the heat flow and polarization current are independent of each other, especially when λ qp is approximately close to zero. Thus, in the majority of microwave-heating problems, electric and thermal fields are considered as weakly coupled.
On the other hand, while the time scale of the component-concentration variation is smaller than the wave period, that is, the component concentration changes rapidly with time, transients become apparent, and frequency changes in the electromagnetic field can be observed [10]. The effect of time-varying characteristics on the polarization current J p is much more obvious, and the rate of change in J p changes rapidly with time. This means that the degree of coupling may vary during the reaction processes. Thus, the effects of the reaction processes on the degree of coupling are noticeable. Most importantly, if the rate of change of J p decreases in the reactions, the reduction leads to the increase in λ qp . The higher the absolute value of λ qp , the higher the dependence of the heat flow on the polarization current, thereby provoking a transition from weak coupling to strong coupling. More importantly, the degree of coupling may affect the accuracy of the numerical solution as well as the computational efficiency, which should be noticed carefully in the simulation of microwave-assisted chemical reactions.
All in all, the coupled conditions under different reaction processes exhibit various degrees of interaction between the electromagnetic and thermal fields. Generally, in the microwave-heating process, the time scale of the component-concentration variation is greater than the wave period. Hence, the electromagnetic-thermal models are considered as weakly coupling, and it is possible to implement an iterative method based on a weakcoupling procedure [21,22]. On the other hand, when the time scale of the componentconcentration variation is smaller than the wave period, the interaction changes profoundly. Classical examples are frequency conversions in the rapidly created plasma [23,24]. Thus, in terms of the multiphysics calculation, numerical calculations need to be converted into strong-coupling algorithms.

Conclusions
In this paper, we have derived the coupling degree of the electromagnetic field and the thermal field in microwave-heating polar-molecule reactions, based on the entropybalance equation. Because of the time-varying characteristics of the reactions, the effects of the reaction processes on the degree of coupling are discussed. If the time scale of the component-concentration variation in reaction processes is much greater than the wave period, the degree of coupling approaches zero and there is a weak coupling between the electric and thermal fields. On the other hand, when the time scale of the componentconcentration variation is smaller than the wave period, the degree of coupling may vary during the reaction processes. When the absolute value of the degree of coupling becomes larger, there may be a change from weak coupling to strong coupling.