The Regularity of Solution for Weakly Coupled System Derived by Microwave Heating Model

: In this paper, we study the regularity of the weak solution of the coupled system derived from the microwave heating model with frequency variable. We ﬁrst show that the weak solution E of the system is Hölder continuous near the boundary of S = ∂ Ω . The main idea of the proof is based on the estimation of linear degenerate system in Campanato space. Then we show that the solution u of the heat conduction equation is Hölder continuous with exponent α 2 . Finally, under the appropriate conditions we show that the coupled system with microwave heating has a weak solution. Moreover the regularity of the weak solution is studied.


Introduction
Microwave heating is a kind of heating method that can be converted into heat energy by absorbing microwave energy from an object. As a new type of energy carrier, it has been widely used in many fields such as food heating and thawing, biological sterilization, chemical synthesis and metal smelting (see [1][2][3][4][5]). Compared with the traditional heating mode, microwave heating has the characteristics of short heating time, high efficiency and small thermal inertia. Hence, the automatic control skill is able to be used in microwave heating process. However, due to the internal heat energy generated by the microwave, the heating of some heated stuff may be uneven or even "runaway heating" (see [6][7][8][9][10]). Therefore, it is a meaningful problem to realize uniform heating and minimum energy dissipation in microwave heating process.
The heat produced by microwaves causes the temperature of the heated object to be unevenly distributed. This is partly because various physical parameters such as electric permittivity, electric conductivity and magnetic permeability strongly depend on the temperature. Understanding the complicated dynamic interactions between electric field, magnetic field and temperature is of great important to the system modeling and design process.
If the variable coefficients of Maxwell's equations are not smooth, it will be difficult to obtain the regularity of the solution for Maxwell's equations. This also leads to the heat source of heat conduction equation is only belong to L 1 (Q T ). Therefore, it is difficult to prove the existence and regularity of solutions. Furthermore, the necessary optimality conditions for optimal control cannot be derived. Under the condition of the coefficient has better smoothness, Yin [9] proved the existence of the solution of the coupled system. But the smoothness of the coefficient is unattainable in the practical problems, even is discontinuity. The study of weak solution to the corresponding coupling system has become a significant problem. At the same time, in order to study necessary optimality conditions of the optimal control problem with frequency variable, we need to study the regularity of the weak solution.
In this paper, we prove that the weak solution of system (4) is continuous. When γ ω is bounded, the real part of γ ω has a positive lower bound and ξ is Lipschitz continuous. Assume that the regularity of the coefficients γ ω and ξ are minimal, then this regularity is optimal. The basic idea is to use Campanato type of estimates(see [11,12]).
We have derived weak coupling mathematical model of microwave heating system with frequency variation in [13]. In paper [13], the existence of global weak solutions is established by assuming suitable conditions. By applying Lax-Milgram theorem and monotone operator theory, the existence and uniqueness of solutions for weakly coupled systems are proved in [13]. However, the regularity of the weak solution has not been studied for the coupled system in [13]. Particularly, we consider the uniformity and boundedness of temperature. In this paper. We study that the temperature is Hölder continuous by appropriate assumptions.The regularity results provide the basis about the derivation of optimal conditions for uniform microwave heating This paper is organized as follows. In Section 2, for completeness we introduce the mathematical model of microwave heating with frequency variable from [13]. In Section 3, we impose some basic assumptions which ensure the well-posedness of the underlying system. Existence of a weak solution is also established from [13]. In Section 4, some important estimates are derived. In Section 5, regularity of the weak solution is investigated.

Weak Coupling Mathematical Model of Microwave Heating System
We consider that the electromagnetic field does not change with temperature during microwave heating. In particular, the distribution of dielectric, permeability, and electric conductivity functions are temperature independent. The frequency of microwave is dependent only on its position x, i.e., ω = ω(x), and has nothing to do with time.
In this paper, we assume that Ω is a bounded simply-connected domain in three-dimensional space R 3 . The boundary of Ω is C 1 continuous. Let Q T = Ω × (0, T] and S T = ∂Ω × (0, T]. Hereafter, a bold letter represents a vector in R 3 .
According to literature [13], the microwave heating process can be described by the weakly coupled system of Maxwell equations and heat conduction equations and the corresponding initial boundary value conditions as follows: where the unknown functions E(x, t) and u(x, t) respectively represent the electric field intensity and temperature at time t at x ∈ Ω, and n is the outward unit normal vector on S = ∂Ω and u n = ∇u · n is the outward the normal derivative on S = ∂Ω. The G(·) is the time-harmonic electromagnetic field generated by the external photoelectric device.
where functions ε, µ and σ represent the dielectric, permeability and electric conductivity respectively.

The Well-Posedness of Solution of Weak Coupling System
V(·) is the vector valued function in three dimensions.

Definition 1.
(see [14]) H(curl, Ω), H(div, Ω) and H(curl, div, Ω) space H(curl, Ω) is the Hilbert space with the following inner product: where F * is the complex conjugate of F (see [9]). The norm is (derived from inner product) We impose some basic assumptions which ensure the well-posedness of the underlying system and proof of regularity.
The following inequality conditions are satisfied: for all x ∈ Ω, u ∈ R with positive constants k 0 and k 1 .
(a) The functions a 1 , a 2 are bounded. Assume that a 1 , a 2 are Lipschity continuous with respect to ω−variables. At the same time there is for all x ∈ Ω and some constant γ 0 > 0. (c) The function G is given and defined on S T with the extended function G ∈ H(curl, Ω). G is Hölder continuous on Ω, and we have G(·, t) H(curl,Ω) ≤ c 0 G(·, t) L 2 (S) for all t, where c 0 is a constant that depends only on Ω. (d) Assumed that the function ω : Ω → R is measurable. Let where ω 0 and ω 1 are positive constants.
Since the controlled system we are considering is a weakly coupled system. Its heating process can be described by Maxwell equations determining the distribution of its electric field intensity.
The intensity of the electric field produces an internal heat source in the heating process, causing the temperature of the heated material to rise. In this way, the material can be heated.
Therefore, in mathematical analysis, we can first study the existence and uniqueness of Maxwell equations. Then we consider the existence and uniqueness of the solution for the initial boundary value problem of heat conduction equation under a given electric field intensity distribution. We consider the system: From Theorem 3.2. of literature [13], the system (4) exists unique solution E ∈ L ∞ (0, T; H 0 (curl, Ω)) for given ω(·) ∈ U ad .
Next, we study the existence and uniqueness of the solution to the heat conduction equation under a given electric field intensity distribution. Consider the following system: Similarly, from Theorem 3.5. of literature [13], there exists a unique solution u ∈ W[0, T] to the problem (5) for any given ω(·) ∈ U ad and its corresponding solution E ∈ L ∞ (0, T; H(curl, Ω)) to (4).
In conclusion, the existence and uniqueness theorem (see [13] Theorem 3.6) of solutions for the coupled system (1).
For weakly coupled systems, we first consider the regularity of the solution to the linear Maxwell Equations (4). Let

Lemma 5. Under the assumption H(3), let q
for any ω ∈ U ad . There exists positive constant α 0 and β 0 such that and there exists a positive constant B 0 such that for F ∈ L 2 (Ω).
The System (4) is degenerate. The classical regularity theory for elliptic system is not valid here. We will prove Hölder continuous of the solution for System (4). The following Lemma plays a key role in the following proof.
in Ω, n · V = 0, on S, and where C 3 depends only on Ω. (4) is Hölder continuous inΩ for any ω ∈ U ad , and there exists a constant C 4 > 0 such that

Theorem 7. Under the assumption H(3), then the weak solution E of the System
where C 4 depends only on known data.
From Lemma 6, we know that there exists a potential vector field V(x) such that Moreover where C depends only on Ω and the upper bound of ξ. Substitute V(x) into Equations (12), and we have From the boundary condition n ×Ẽ = 0, we have For any smooth function Ψ ∈ C ∞ (Ω), we can use Gauss's divergence theorem to get Hence, n · (∇ ×Ẽ) = 0, x ∈ ∂Ω.
Since Ω is simply-connected, there exists the potential function Ψ(x). The following equation holds Therefore, By calculation that is ω(x) has a positive lower bound, we use L 2 (Ω)-theory of elliptic equations to get where C depends only on known data. According to the estimate of Campanato space L 2,µ (Ω), we see that there exist µ ∈ (1, 2) such that which is bounded from the estimate for V(x). From the Equations (15), one can see that On the other hand, from (13) we know From the boundedness of ξ(x), we have for ∀x ∈ Ω and ρ > 0, where B ρ is B ρ ∩ Ω andẼ ∈ H 1 (Ω). From Lemma 4, we know thatẼ ∈ L 2,2 (Ω) and From what has been discussed above, we get where C only depends on the known data constant, i.e., From Lemma 4, the space dimension n = 3 and µ 0 = 2 + µ > 3. We see that where α = µ−n 2 > 0. HenceẼ is Hölder continuous. LetẼ = E − G, we have E =Ẽ + G. From the hypothesis of G in assumption H(3), we know that E ∈ C α (Ω).
Next we study the regularity of the heat conduction equation. Let for B(x 0 , ρ). One has a similar result to Lemma 4 with dimension n replaced by n + 1.
We consider the following parabolic equation The following basic assumptions are needed. H(4).
(a) Let a ij (x, t) be measurable in Q T and satisfy the following conditions for ellipticity: (c) Let for some µ > 0.

Lemma 8.
(see [17]) Under the assumption H(4), the solution u ∈ C([0, T]; H 1 (Ω)) of the general parabolic equation satisfies the following estimate: where C is a constant depending only on known data a 1 , a 2 , A 3 and Q T . In particular, there is for µ ∈ (n, n + 2), where α = µ−n 2 .
In particular, Proof. From Theorem 3.6. of literature [13], we know that there exists the solution (E, u) for coupled system (1) with We assume that G(x, t) = 0 on ∂Ω. Otherwise, we definê From Theorem 7, we have E ∈ C α (Ω), for any given t ∈ [0, T]. Moreover Now we think that k(x, u) is independent of u. Let k(x, u) =k(x), we use L 2,µ theory of parabolic equation and Lemma 8. We have ∇u L 2,2+µ 0 (Q T ) ≤ C ∇ × E L 2,2+µ (Q T ) .
Since k(x, u) is smooth, we can use Schauder's theorem to get the following result u ∈ C 2+α,1+ α 2 (Q T ).
Remark 11. The necessary optimality conditions for the optimal control problem in the microwave heating process will be carried out in a separate paper.

Conclusions
In this paper, we study the regularity of the weak solution of the coupled system derived from the microwave heating model with frequency variable. The distribution of dielectric, permeability, and electric conductivity functions are temperature independent. The frequency of microwave is dependent only on its position x, i.e., ω = ω(x), and has nothing to do with time. We established that the temperature is Hölder continuous, even if electric conductivity has a jump discontinuity with respect to the temperature change. This paper presents a method that is based on the estimation of linear degenerate system in the Campanato space. We use this idea to deal with time harmonic Maxwell equations with rough coefficients. On the other hand, by using the regularity results of the coupled system, the necessary optimality conditions for uniform microwave heating in three directions are derived.