Local heating realization by reverse thermal cloak.

Transformation thermodynamics, as one of the important branches among the extensions of transformation optics, has attracted plentiful attentions and interests recently. The result of transformation thermodynamics, or called as "thermal cloak", can realize isothermal region and hide objects from heat. In this paper, we presented the concept of "reverse thermal cloak" to correspond to the thermal cloak and made a simple engineering definition to identify them. By full-wave simulations, we verified that the reverse thermal cloak can concentrate heat and realize local heating. The performance of local heating depends on the anisotropic dispersion of the cloaking layer's thermal conductivity. Three-dimensional finite element simulations demonstrated that the reverse thermal cloak can be used to heat up objects. Besides pre-engineered metamaterials, such reverse thermal cloak can even be realized with homogenous materials by alternating spoke-like structure or Hashin coated-sphere structure.

With the born of transformation optics in 200612, it grows prosperously with a lot of attention and has turned into a subject of great interest since it not only provides a general method in coordinate transformation but also provides an effective design tool for manipulating electromagnetic wave3. In the past few years, transformation optics has been well developed beyond the domain of electromagnetic wave to plentiful extensions in other fields, including acoustic wave45, liquid wave6, matter wave7, and elastic wave89, etc. The results of these transformations, or called as cloaks, can be realized by pre-engineered metamaterials or even homogeneous natural materials with special-designed structures. One common type of cloaks is the core-shell structure, which maps a concealment region into a surrounding shell region with inhomogeneous distribution of constitutive parameters101112. The core-shell structure, like the typical Hashin coated sphere structure, was initiatively described by Hashin913 as a model of composite comprised of spherical grains of one phase embedded in a matrix of another phase. Hashin pointed out how to design the phase materials and the corresponding fractions so that inserting the coated sphere into the matrix would not disturb the uniform electromagnetic field outside14.

Analogy to transformation optics, Guenneau et al.15 first extended the transformation methodology to the area of thermodynamics and proposed the concept of thermal cloak to channel the heat flow to realize a local isothermal region where the temperature gradient vanishes. Further studies16171819 demonstrated that thermal cloak can be realized with homogenous and finite thermal conductivity by simulations and experiments. More specifically, Han et al.16 presented that thermal cloak can be realized by isotropic natural materials; Schittny et al.17 designed and manufactured a thermal cloak comprised of copper rings interlaid with plastic. From the view of guiding heat flux, Narayana and Sato18 did excellent and systematic work about controlling the heat flow by the engineering materials. By shielding, concentrating, and inverting the heat current, they experimentally showed that artificially engineered thermal materials can yield characteristics that go far beyond the conventional media.

The concept of thermal cloak, in fact, is derived from the optical invisible cloak, and it is aimed at hiding object from heat or rendering heat invisible in certain region. The thermal cloak can be achieved by anisotropic thermal conductivity distribution where the azimuthal component is much larger than the radial component (>). Inspired by the previous work, one possible idea turns up that could we apply the thermal cloak reversely (>) and realize a local heating region where the temperature gradient is large enough to heat up objects?

By utilizing the transformation thermodynamics and tailoring the dispersion of thermal conductivity, here we demonstrated that the reverse application of thermal cloak can efficiently realize a local heating region where the temperature gradient is more considerable than that of the outside region. Steady full-wave simulations show that the temperature profile in the central region strongly depends on the anisotropy of thermal conductivity. Transient simulations illustrate that heat diffuses faster in the left part of the reverse thermal cloak and slower in the right part, thus large temperature gradient is formed in the central region. By modeling a substrate consisting of reverse thermal cloak and assembling a cubic probe onto the central region, the corresponding temperature fields on the cubic probe show that local heating effect was enhanced with reverse thermal cloak. Two possible methods were also presented to realize such reverse thermal cloak.

Results

With transformation thermodynamics, we derived that the radial thermal conductivity and the azimuthal thermal conductivity are the inverse of each other, as shown in Eq. (1).

By defining , where C is a constant, a cloak can be realized with constant thermal conductivity. For thermal cloaks, the constant C is usually less than 116, implying the radial conductivity is much smaller than the azimuthal conductivity . The performance characteristics of thermal cloak are totally attributed to its anisotropic distribution of thermal conductivity. In this paper, we proposed the concept of “reverse thermal cloak” where is much larger than , i.e. the constant C is larger than 1.

With full-wave simulations based on finite element method (FEM), the steady-state performances of the cloaks were examined. The large commercial software package COMSOL Multiphysics was used to conduct the simulations. For comparison, a homogeneous plane with constant thermal conductivity was also simulated. In the simulations, the left and right boundaries were kept at constant temperature, i.e. 400 K and 300 K, respectively. The top and bottom boundaries are thermally insulated. The interior and exterior radii of the annular thermal cloak (reverse thermal cloak) are 1 cm and 2 cm, respectively. The dimensions of the whole domain are 10×10 cm2. The thermal conductivities of the core region and the substrate are kept as 0.1 W/(m·K) and 1.0 W/(m·K), respectively.

Fig. 1 shows the temperature profiles and corresponding temperature curves of the section lines on a homogenous plane, a thermal cloak and a reverse thermal cloak, respectively. From Fig. 1(a), the temperature profile of the homogenous plane is uniform and the temperature gradient is constant throughout the plane. From Fig. 1(b), the temperature profiles outside the thermal cloak are all the same with Fig. 1(a), but the heat tends to be away from the inner core region and travels close to the outside edge of the thermal cloak (C = = 0.1 W/(m·K), = 10 W/(m·K)). From both the temperature contour and the temperature gradient of the section line, we can see that thermal cloak can realize an isothermal region in the central domain. This implies that any object placed in the interior region appears to be concealed since the outside heat flow looks like as if there is no disturbance at all in the central region (the temperature gradient outside the thermal cloak is the same with that in Fig. 1(a)). This is the reason and origin of the definition of “thermal cloak”. Fig. 1(c) shows that temperature field of the reverse thermal cloak in which we adjust the constant C to 10 (C = = 10 W/(m·K), = 0.1 W/(m·K)). We can see that the temperature profiles bend toward the central region and the temperature gradient in the central region is enhanced greatly. More vivid comparison between the thermal cloak and the reverse thermal cloak can be found in Fig. 2. The reverse thermal cloak tends to concentrate the heat and have the most heat diffuse from it, resulting in local heating and large temperature gradient. Isothermal region forms in the thermal cloak while a sharp temperature drop arises in the reverse thermal cloak.

Figure 1

Temperature profile comparison among a homogenous plane, an annular thermal cloak and an annular reverse thermal cloak with constant boundary temperatures at 400 K and 300 K respectively: (a) a plane and the corresponding temperature curve of section line at y = 0; (b) a thermal cloak and the corresponding temperature curve of section line at y = 0, (C = = 0.1, = 10); (c) a reverse thermal cloak and the corresponding temperature curve of section line at y = 0 (C = = 10, = 0.1).

Isotherms are illustrated in green color in panel. Cyan dash lines are the section lines.

Figure 2

Visual illustration of the temperature profiles of (a) a thermal cloak and (b) a reverse thermal cloak.

To further demonstrate the local heating phenomenon, we plotted the temperature curves on the cut plane (x = −0.5 cm) in Fig. 3. Obvious local high temperature can be observed in the central region. We can also see that when the constant C increases from 5 to 20, the local temperature in the central region is enhanced. The larger the constant C is, the larger the thermal conductivity anisotropy is. The trend in Fig. 3 implies that large anisotropy of the thermal conductivity distribution leads to enhanced local heating effect.

Figure 3

Temperature profiles of a cut plane (x = −0.5 cm) on an annular reverse thermal cloak with changing the constant C from 5 to 20.

Transient simulations of the reverse thermal cloak were carried out and the heat diffusion processes are shown in Fig. 4 (C = = 10 W/(m·K), = 0.1 W/(m·K)). In the beginning, heat has not conducted to the right part yet and the isotherms squeeze in the left part. With time elapse, heat conducts to the right part gradually and the isotherms stretch to the whole domain. At t = 60 second, the thermal system almost get its equilibrium state and the corresponding isotherms are symmetric in the left and right parts. No matter in the transient state (see Fig. 4) or in the steady state (see Fig. 1 and Fig. 2), we can discover that the isotherms bend toward the central region in both the left and right parts. This implied that in the reverse thermal cloak, heat conducts faster in the left part and slower in the right part. As a result, compact isotherms and large temperature gradient form in the central region.

Figure 4

Transient heat conduction through an annular reverse thermal cloak with the left and right boundaries at constant temperature.

Snapshots of temperature profile with time elapse from t = 5 s to t = 60 s. Isotherms are illustrated with green color in panel.

From above analyses, we can observe the local heating effect and there exists a large temperature gradient in the central region of the reverse thermal cloak. Since the driving force of heat diffuse is the temperature gradient, one idea is that we could utilize the large temperature gradient to heat up objects. To verify our idea, we built a three-dimensional model in which we modeled a substrate consisting of a reverse thermal cloak (C = = 10 W/(m·K), = 0.1 W/(m·K)) and then mounted a cubic probe on the central region, as shown in Fig. 5(a). The dimensions of the substrate are 10×10×0.5 cm3. The dimensions of the cubic probe are 1×1×0.5 cm3. The interior and exterior radii of the annular reverse thermal cloak are 1 cm and 2 cm, respectively. For comparison, we also built the same-dimensional model except that the substrate is a homogeneous plane without reverse thermal cloak. After steady-state simulations, we selected the temperature contours on the front surface of the cubic probe (on the plane of x = −0.5 cm) to demonstrate our idea. The temperature comparisons are shown in Fig. 5(b) and Fig. 5(c). We can see that the maximum temperature of the cubic probe is 370 K when the substrate is with a reverse thermal cloak; the maximum temperature is 357 K when the substrate is without a reverse thermal cloak. Thus the reverse thermal cloak can be utilized to heat up objects.

Figure 5

(a) 3-D model of a substrate consisting of a reverse thermal cloak and a cubic probe, (b) temperature contours of the cubic probe on a reverse thermal cloak, and (c) temperature contours of the cubic probe on a homogenous plane with the same dimensions.

To realize such a reverse thermal cloak whose radial thermal conductivity is smaller than its azimuthal thermal conductivity, we designed a composite structure made by alternating spoke-like layers of two homogeneous materials. As shown in Fig. 6, the structure consists of two isotropic materials, i.e. material A and material B. The thermal conductivities of these two materials are 10 W/(m·K) and 0.1 W/(m·K), respectively. All the boundary conditions are kept the same. The simulated temperature profile is similar to above simulations where compact isotherms and local heating effect form in the central region. This implies that the reverse thermal cloak not only can be realized by pre-engineered metamaterials, but also can be realized by homogeneous materials with alternating spoke-like structures.

Figure 6

Simulation of a reverse thermal cloak consisting of alternating spoke-like layers of two homogeneous materials (A and B): (a) the temperature profile and (b) its structure.

Another method to realize such kind of reverse thermal cloak with homogeneous materials can ascend to the Hashin coated sphere. According to the two-dimensional Hashin-Shtrikman formula14, the effective thermal conductivity can be calculated asWithwhere , are the thermal conductivities of the core region and the shell region, respectively. f and f are the two-dimensional volume fraction, r and r are the exterior radii of the circular core and the cloaking shell. When equals to the thermal conductivity of the outside matrix , reverse thermal cloak can be realized.

According to dimensions and thermal conductivities mentioned above, the thermal conductivity of the cloaking shell can be calculated as 0.918 W/(m·K) in this paper. With the same method, we simulated the temperature field and the result was plotted in Fig. 7. It is seen that the green isotherms and the white streamlines are orthometric to each other in panel. The isotherms bend to the core region and a large temperature gradient forms there. Therefore, with the Hashin-Shtrikman formula, we can design the dimensions and thermal conductivities to realize reverse thermal cloaks. But when comparing Fig. 7 with Fig. 4 and Fig. 6, it is also seen that the thermal gradient in Hashin coated sphere is not as compact as that realized by pre-engineered metamaterials, nor homogeneous materials with alternating spoke-like structures. Thus it is perceived that the local heating effect in Hashin coated sphere is not as strong as the other two, and the heat-concentration ability of Hashin coated sphere is also limited.

Figure 7

Simulation of another reverse thermal cloak comprised of homogenous materials.

The thermal conductivity of the reverse thermal cloak is 0.918 W/(m·K), which is designed according to the Hashin-Shtrikman theory. The isotherms and streamlines are also illustrated with green and white colors in panel, respectively.

Discussion

In summary, we demonstrated in this study that reverse thermal cloak can be applied to realize local heating. Unlike the thermal cloak can realize isothermal region and render heat invisible from the outside temperature field, the reverse thermal cloak can bend the isotherms to the central region and realize large temperature gradient there. Due to the large temperature gradient, local heating phenomenon was observed. Such “local heating” phenomenon can be enhanced by re-distributing the dispersion of thermal conductivity in the domain of the reverse thermal cloak. Here, a way to realize local heating was proposed by the reverse thermal cloak. Like the fabrication of thermal cloak, the reverse thermal cloak can also be fabricated with pre-engineered metamaterial or even with homogenous materials by spoke-like structures or Hashin coated-sphere structure. With the reverse thermal cloak, we can realize the local heating in many fields. For instance, in the electronic packaging, with local heating, we may avoid heating some temperature-sensitive components or causing printed circuit board (PCB) warping.

It should be pointed out that the reverse thermal cloak was proposed to correspond to the thermal cloak during the research of thermal cloak. In this paper, we simply defined the difference between these two cloaks as that: for thermal cloak, the constant C ( = radial thermal conductivity ) is less than 1; for the reverse thermal cloak, the constant C is larger than 1. The essence of both thermal cloak and reverse thermal cloak is the anisotropic dispersion of thermal conductivity. When compared with the thermal concentrator1518, the present reverse thermal cloak exhibits similar characteristics but in a much simpler engineering way. The thermal concentrator has specific formula to identify the dispersion of thermal conductivity15, while in this paper we extended the limitation of thermal conductivity and gave a general formula in Eq. (1). Moreover, we paid more attentions to the application of reverse thermal cloak, rather than just staying in the theoretical or laboratory stages. What should be emphasized is that the present reverse thermal cloak can also be applied to other shape cloaks, like square, rectangular, elliptic, etc.

Methods

When a temperature gradient exists in a body, experience has shown that there is an energy transfer from the high-temperature region to the low-temperature region. For a steady state without heat source, the control equation can be written as , where is the thermal conductivity and is the temperature. Similar to the invariance of Maxwell equations in the transformation optics, Guenneau et al.15 pointed out that the heat conduction equations can also be written in a form-invariant manner under coordinate transformations. Thus, in the transformed space, the heat conduction equation can be re-written asWithwhere is the Jacobian matrix with components defined as . The Jacobian matrix characterizes the geometrical variation in the original space and the transformed space . Eqs. (2) and (3) construct the basis of transformation thermodynamics and we can design and manipulate the flow of heat by an arbitrary coordinate transformation.

One seminal typical design of thermal cloak was consulted to the invisible cloak proposed by Pendry et al.110 With the simple coordinate transformation in Eq. (4), we can squeeze the heat flow from a central region () in the original space into a shell region () in the transformed space , while the rest of heat flow (in the region ) are maintained.

where R and R are the interior and exterior radii of the annular thermal cloak. Substituting Eq. (6) into Eqs. (4) and (5) yields Eq. (1).

Author Contributions

R.H., X.L.W., J.Y.H. and X.B.L. contributed in the discussion. R.H. and J.Y.H. participated in the simulation. R.H. prepared the manuscript. R.H. and X.B.L. conceived the idea and revised the manuscript. X.B.L. supervised the project.