Ellipsoidal Thermal Concentrator and Cloak with Transformation Media

Yong Gao(高勇)$^{\hyperlink{s}{*}}$

doi:10.1088/0256-307X/38/2/020501

⇦Chinese Physics Letters, 2021, Vol. 38, No. 2, Article code 020501⇨ Ellipsoidal Thermal Concentrator and Cloak with Transformation Media Yong Gao (高勇)* Affiliations Department of Physics, Research Center of Resource Recycling Science and Engineering, Shanghai Polytechnic University, Shanghai 201209, China Received 10 October 2020; accepted 9 November 2020; published online 27 January 2021 Supported by the Gaoyuan Discipline of Shanghai-Environmental Science and Engineering (Resource Recycling Science and Engineering).
^*Corresponding author. Email: gaoyong@sspu.edu.cn Citation Text: Gao Y 2021 Chin. Phys. Lett. 38 020501 Abstract Thermal concentrators and cloaks with ellipsoidal shapes are designed by utilizing the transformation thermotics method and finite element simulations. The thermal conductivities for the concentrator and cloak are directly derive in Cartesian coordinates. The simulation results show that the ellipsoidal thermal concentrator can focus heat flux into a central region and that the ellipsoidal thermal cloak can guide heat flux around the cloaked region smoothly without disturbing the external temperature distribution and heat flux. The present method can be extended to design arbitrarily shaped thermal metadevices with novel properties. DOI:10.1088/0256-307X/38/2/020501 © 2021 Chinese Physics Society Article Text In recent years, the developed technique of transformation optics (TO)[1–6] has become a powerful means to control electromagnetic waves in almost arbitrary ways. The TO method is based on the form invariant of Maxwell's equations under coordinate transformations.[1,7] However, this approach is not limited only to electromagnetic waves. For example, the TO method has been extended to other fields including acoustics,[8] thermotics,[9–27] elastic waves,[28] matter waves,[29] fluid flow,[30] etc. The fundamental reason for such an extension is that the governing equations describing these systems are form invariant under coordinate transformations. Over the last few decades, increasing attention has been paid to renewable energy. In addition to solar energy, another important energy source is thermal energy, which can be collected from environment. Therefore, it is particularly important to manipulate heat transfer at will. To date, studies on thermal concentrators have been limited to cylinders[12,13] and spheres.[31] In the present work, we investigate an ellipsoidal thermal concentrator and cloak for heat flux using the transformation thermotics method. For convenience, we directly derive the thermal conductivity tensors for the concentrator and cloak in Cartesian coordinates. With the derived material parameters, the finite element simulation results are analyzed to display the corresponding functions. For the designed ellipsoidal thermal concentrator, the heat flux can be concentrated in the central region. However, for the ellipsoidal thermal cloak, the heat flux can be excluded from an object without perturbing the exterior heat flux. For the designed devices, the heat flux can be guided back to its original path. Based on the transformation thermotics method and thermal metamaterials, we can design thermal concentrator and cloak with different shapes, and furthermore, develop various thermal devices with novel functionalities, such as thermal diodes.[21] The extension of the TO method to the thermal field could open a new field of study since it promises unprecedented manipulation of the heat flux, which may lead to new phenomena and applications. Let us start by investigating the heat conduction equation. At the equilibrium, the stationary heat conduction equation in the absence of sources is given by $$ {\bf \nabla}\cdot (\kappa {\bf \nabla} T)=0,~~ \tag {1} $$ where $T$ is the temperature state variable, and $\kappa$ is the material thermal conductivity. Because the heat conduction equation is invariant in its form under the transformation from an original coordinate to a transformed coordinate,[7] in the new coordinate system, we must use the renormalized value of the thermal conductivity: $$ \kappa ^{i^{\prime }j^{\prime }}=|\mathrm{{\det}}(\varLambda _{i}^{i^{\prime }})|^{-1}\varLambda _{i}^{i^{\prime }}\varLambda _{j}^{j^{\prime }}\kappa ^{ij},~~ \tag {2} $$ where the Roman indices extend from 1 to 3 for the three spatial coordinates, and $\varLambda _{i}^{i^{\prime }}$ is the transformation operator. The bases denoted by the indices $i$ and $i^{\prime}$ are coordinate bases. The transformation operator is derived as $$ \varLambda _{i}^{i^{\prime }}=\frac{\partial x^{i^{\prime}}}{\partial x^{i}}.~~ \tag {3} $$ The coordinate transformation method is valid for any complicated transformation. As an example, we investigate an ellipsoidal concentrator. The coordinate transformation used in the studies of spherical cloaks[1,3] can be readily generalized to ellipsoids. For an ellipsoid with its semi-axes lying along the $x$, $y$, and $z$ directions with lengths $\alpha_{1}c$, $\alpha_{2}c$, and $\alpha_{3}c$, respectively, the ellipsoidal equation can be expressed as $$ \Big(\frac{x^{1}}{\alpha_{1} c}\Big)^{2}+\Big(\frac{x^{2}}{\alpha_{2} c}\Big)^{2}+\Big(\frac{x^{3}}{\alpha_{3} c}\Big)^{2}=1,~~ \tag {4} $$ where $\alpha_{1}$, $\alpha_{2}$, and $\alpha_{3}$ are positive numbers describing the aspect ratio of the ellipsoid and do not transform as a vector. For a prolate sphere, we choose $\alpha_{1} = \alpha_{2} < \alpha_{3}$. For an oblate sphere, we choose $\alpha_{1} < \alpha_{2} = \alpha_{3}$. For the design, we obtain three ellipsoidal equations by substituting $c=R_{1}$, $c=R_{2}$ and $c=R_{3}$ into Eq. (4). We redefine $r$ as a scaled distance:[32] $$ r=\Big[\sum_{i=1}^{3}\Big(\frac{x^{i}}{\alpha_{i}}\Big)^2\Big]^{1/2}.~~ \tag {5} $$ For the ellipsoidal thermal concentrator, we adopt the following linear transformations: $$\begin{alignat}{1} &r^{\prime}=\frac{R_{1}}{R_{2}}r, ~~~~~~0 \leqslant r \leqslant R_{2}, \\ &r^{\prime}=\frac{R_{3}-R_{1}}{R_{3}-R_{2}}r+\frac{R_{1}-R_{2}}{R_{3}-R_{2}}R_{3},~~R_{2} < r \leqslant R_{3}.~~ \tag {6} \end{alignat} $$ Now, since our transformation is radially symmetric, the unit vectors in the transformed space and in the original space must be equal to each other:[3] $$ \frac{x^{i^{\prime}}}{r^{\prime}}=\frac{x^{i}}{r}\delta_{i}^{i^{\prime }}.~~ \tag {7} $$ The position vector in the transformed space in terms of only the components in the original space for the region of the ellipsoid with semi-axes $\alpha_{1}R_{2}$, $\alpha_{2}R_{2}$, and $\alpha_{3}R_{2}$ can be written as $$ x^{i^{\prime}}=\frac{R_{1}}{R_{2}}x^{i}\delta_{i}^{i^{\prime }}.~~ \tag {8} $$ Thus, we can obtain the transformation matrix $$ \varLambda_{j}^{i^{\prime }}=\frac{r^{\prime }}{r}\delta_{j}^{i^{\prime }}.~~ \tag {9} $$ For the ellipsoidal shell with lengths of outer semi-axes $\alpha_{1}R_{3}$, $\alpha_{2}R_{3}$, and $\alpha_{3}R_{3}$, and lengths of inner semi-axes $\alpha_{1}R_{2},\alpha_{2}R_{2}$, and $\alpha_{3}R_{2}$, the position vector can be written as $$ x^{i^{\prime}}=\frac{R_{3}-R_{1}}{R_{3}-R_{2}}x^{i}\delta_{i}^{i^{\prime }}+\beta\frac{x^{i}}{r}\delta _{i}^{i^{\prime }}.~~ \tag {10} $$ Then, we can obtain the transformation matrix $$ \varLambda_{j}^{i^{\prime }}=\frac{r^{\prime }}{r}\delta_{j}^{i^{\prime }}-\frac{\beta x^{i}x^{k}\delta _{i}^{i^{\prime }}\delta _{k j}}{\alpha_{j}^{2}r^{3}},~~ \tag {11} $$ where $\beta$ can be written as $$ \beta=\frac{(R_{1}-R_{2})R_{3}}{R_{3}-R_{2}}.~~ \tag {12} $$ The transformation in Eq. (8) maps the field in the region of the ellipsoid with the semi-axes $\alpha_{1}R_{2}$, $\alpha_{2}R_{2}$, and $\alpha_{3}R_{2}$ onto the region of the ellipsoid with the semi-axes $\alpha_{1}R_{1}$, $\alpha_{2}R_{1}$, and $\alpha_{3}R_{1}$. Equation (10) transforms an ellipsoidal shell with lengths of outer semi-axes $\alpha_{1}R_{3}$, $\alpha_{2}R_{3},$ and $\alpha_{3}R_{3}$ and lengths of inner semi-axes $\alpha_{1}R_{2}$, $\alpha_{2}R_{2}$, and $\alpha_{3}R_{2}$ into an ellipsoidal shell with the same lengths of outer semi-axes and lengths of inner semi-axes $\alpha_{1}R_{1}$, $\alpha_{2}R_{1}$ and $\alpha_{3}R_{1}$. Using Eq. 2, we can obtain the material parameters in the region of the ellipsoid with the semi-axes $\alpha_{1}R_{1}$, $\alpha_{2}R_{1}$, and $\alpha_{3}R_{1}$: $$ \overleftrightarrow{\kappa}=\begin{pmatrix} {R_{2}}/{R_{1}} & 0 & 0 \\ 0 & {R_{2}}/{R_{1}} & 0 \\ 0 & 0 & {R_{2}}/{R_{1}} \end{pmatrix}\kappa_{0},~~ \tag {13} $$ where $\kappa_{0}$ is the thermal conductivity in the original coordinates. The material parameters in an ellipsoidal shell with lengths of outer semi-axes $\alpha_{1}R_{3}$, $\alpha_{2}R_{3}$, and $\alpha_{3}R_{3}$ and lengths of inner semi-axes $\alpha_{1}R_{1}$, $\alpha_{2}R_{1}$, and $\alpha_{3}R_{1}$ can be expressed as $$\begin{align} \kappa^{ij}={}&\frac{R_{2}-R_{3}}{R_{1}-R_{3}}\Big[\delta ^{ij}-x^{i}x^{j}\Big(\frac{\beta}{r^{3}}\Big(\frac{1}{\alpha_{i}^{2}}+\frac{1} {\alpha_{j}^{2}}\Big)\\ &-\frac{{\beta}^{2}}{\tilde{r}^{2}r^{4}}\sum_{k=1}^{3}\frac {(\tilde{x}^{k})^{2}}{\alpha_{k}^{4}}\Big)\Big]\kappa_{0} ,~~ \tag {14} \end{align} $$ where $\tilde{r}$ and $\tilde{x}^{k}$ read $$\begin{align} \tilde{r}={}&\frac{(r-\beta)(R_{3}-R_{2})}{R_{3}-R_{1}},~~ \tag {15} \end{align} $$ $$\begin{align} \tilde{x}^{k}={}&\frac{x^{k}}{\frac{R_{3}-R_{1}}{R_{3}-R_{2}}+\frac{\beta}{r}}.~~ \tag {16} \end{align} $$ Next, we investigate the ellipsoidal thermal cloak. For the design, we adopt the linear transformation $$ r^{\prime}=\frac{b-a}{b}r+a,~~~(0 < a < b).~~ \tag {17} $$ Based on Eqs. (7) and (17), we can obtain the components of the position vector in the transformed space in terms of only the components in the original space: $$ x^{i^{\prime }}=\frac{b-a}{b}x^{i}\delta_{i}^{i^{\prime }}+a\frac{x^{i}}{r}\delta _{i}^{i^{\prime }}.~~ \tag {18} $$ Therefore, the coordinate transformation given by Eqs. (17) and (18) transforms a closed ellipsoid with $\alpha_{1}b$, $\alpha_{2}b$, and $\alpha_{3}b$ into an ellipsoidal shell with the same lengths of outer semi-axes and lengths of inner semi-axes $\alpha_{1}a$, $\alpha_{2}a$, and $\alpha_{3}a$. By using the same method discussed above, the material property in the ellipsoidal shell is then given by $$\begin{align} \kappa^{ij}={}&\frac{b}{b-a}\Big\{ \delta ^{ij}-x^{i}x^{j}\Big[\frac{a}{r^{3}}\Big(\frac{1}{\alpha_{i}^{2}}+\frac{1} {\alpha_{j}^{2}}\Big)\\ &-\frac{a^{2}}{\tilde{r}^{2}r^{4}}\sum_{k=1}^{3}\frac{(\tilde{x}^{k})^{2}} {\alpha_{k}^{4}}\Big]\Big\}\kappa_{0} .~~ \tag {19} \end{align} $$ Here, $\tilde{r}$ and $\tilde{x}^{k}$ are given by $$ \tilde{r}=\frac{r-a}{b-a}b,~~ \tag {20} $$ $$ \tilde{x}^{k}=\frac{x^{k}}{\frac{b-a}{b}+\frac{a}{r}}.~~ \tag {21} $$ We have obtained the thermal conductivity tensor of an ellipsoidal thermal cloak. In the special case of spheres, $\alpha_{i} = 1$, the coefficient of Eq. (19) reduces to the simple result, which is exactly the same form as Eq. (20) in Ref. [3].

Fig. 1. Three-dimensional finite element simulation of a prolate spheroidal concentrator with parameters $\kappa_0=1$ W/mK, $\alpha_{1}=1$, $\alpha_{2}=1$, $\alpha_{3}=2$, $R_{1}=1$ cm, $R_{2}=2$ cm, and $R_{3}=3$ cm.

Fig. 2. The same as Fig. 1 but the slice in the $zy$ plane.

We are now in a position to perform 3D finite element simulations based on the commercial software COMSOL Multiphysics.[33] In the simulations, the boundary conditions are set to be 100 K and 300 K on the lower and upper boundaries of the cuboid simulation area, respectively, with thermal insulation for the other boundaries. The thermal conductivity of the background is set to be 1 W/mK. The color represents the temperature distribution, and the red lines denote the pathway of heat flux. For the ellipsoidal concentrator, we choose the parameters $\kappa_0=1$ W/mK, $\alpha_{1}=1, \alpha_{2}=1$, $\alpha_{3}=2$, $R_{1}=1$ cm, $R_{2}=2$ cm, and $R_{3}=3$ cm. Figure 1 shows the results of the 3D simulation of a prolate spheroidal concentrator. It is proven to be well suited for the confinement of thermal energy from outside, which may be used for thermally enhanced sensing.

Fig. 3. The same as Fig. 1 but for oblate spheroidal concentrator with parameters $\alpha_{1}=1$, $\alpha_{2}=2$, and $\alpha_{3}=2$, and the slice is in the $zy$ plane.

Fig. 4. Three-dimensional simulation of a prolate spherical cloak with $\alpha_{1}=1$, $\alpha_{2}=1$, $\alpha_{3}=2$, $a=1$ cm, and $b=2$ cm.

Figure 2 illustrates the results of the 3D simulation of a prolate spheroidal concentrator with $\alpha_{1}=1$, $\alpha_{2}=1$, and $\alpha_{3}=2$. The slice is in the $zy$ plane. Figure 3 illustrates the results of the 3D simulation of an oblate spheroidal concentrator with $\alpha_{1}=1$, $\alpha_{2}=2$, and $\alpha_{3}=2$. The slice is in the $zy$ plane. The proposed ellipsoidal thermal concentrator is proven to be well suited for the confinement of thermal energy from outside. Figure 4 illustrates the results of the 3D simulation of a prolate spheroidal cloak with $\alpha_{1}=1$, $\alpha_{2}=1$, and $\alpha_{3}=2$. The heat flux is smoothly bent around the center of the ellipsoidal cloak. Outside the cloak, the heat flux returns to the original direction as if nothing happens.

Fig. 5. The same as Fig. 4 but the slice in the $zy$ plane.

Fig. 6. The same as Fig. 4 but for oblate spherical cloak with parameters: $\alpha_{1}=1$, $\alpha_{2}=2$, and $\alpha_{3}=2$, and for the slice in the $zy$ plane.

Figure 5 illustrates the results of the 3D simulation of a prolate spherical cloak with $\alpha_{1}=1$, $\alpha_{2}=1$, and $\alpha_{3}=2$. The slice is in the $zy$ plane. Figure 6 illustrates the results of the 3D simulation of an oblate spherical cloak with $\alpha_{1}=1$, $\alpha_{2}=2$, and $\alpha_{3}=2$. The slice is in the $zy$ plane. The heat flux is smoothly bent around the center of the ellipsoidal cloak. In summary, we have theoretically and numerically demonstrated transformation media based devices for controlling heat flux. An ellipsoidal thermal concentrator and cloak are designed, and the corresponding thermal conductivities are derived explicitly using the transformation thermotics method. With the derived material parameters, 3D finite element simulations have been performed to confirm the theoretical predictions. The numerical results show that the ellipsoidal concentrator can concentrate heat flux in a small region, which can be used to collect thermal energy from the environment. The ellipsoidal cloak can mold the dynamic flow of heat around an object as if no object is there. Such thermal cloaks can temporarily protect sensitive regions from excessive heating. Our theory can be extended to investigate irregularly shaped concentrators and cloaks, such as rounded cuboids and rounded cylinders. In principle, we can design thermal concentrator and cloak with arbitrary geometries and can also suggest some thermal devices with novel properties. Based on thermal metamaterials, the developed technique of transformation thermotics provides a powerful means to precisely control the heat flux in almost arbitrary ways. The devices suggested in our work may have potential applications in devices such as solar thermal panels and thermal-electric devices. The devices could be probably achieved using composite materials or a layered structure with the designed material parameters. These kinds of artificially engineered thermal metamaterials can yield characteristics that far exceed those available with natural materials, and such thermal metamaterials should lead to unprecedented heat flux control and intriguing energy-related applications. References Controlling Electromagnetic FieldsMetamaterial Electromagnetic Cloak at Microwave FrequenciesCalculation of material properties and ray tracing in transformation mediaDesign of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s equationsTransformation optics and metamaterialsRecent advances in transformation opticsOn cloaking for elasticity and physical equations with a transformation invariant formAcoustic cloaking in three dimensions using acoustic metamaterialsShaped graded materials with an apparent negative thermal conductivityA bifunctional cloak using transformation mediaAnisotropic conductivity rotates heat fluxes in transient regimesHeat Flux Manipulation with Engineered Thermal MaterialsTransformation thermodynamics: cloaking and concentrating heat fluxExperiments on Transformation Thermodynamics: Molding the Flow of HeatStructured thermal surface for radiative camouflageThermal Illusion of Porous Media with Convection-Diffusion Process: Transparency, Concentrating, and CloakingChameleonlike metashells in microfluidics: A passive approach to adaptive responsesNegative Thermal Transport in Conduction and AdvectionFull-Parameter Omnidirectional Thermal Metadevices of Anisotropic GeometryTemperature-Dependent Transformation Thermotics: From Switchable Thermal Cloaks to Macroscopic Thermal DiodesMulti-Physics Bi-Functional Intelligent Meta-Device Based on the Shape Memory AlloysThermal cloak-concentratorDesign of layered structure for thermal cloak with complex shapeThermal Metamaterial: Fundamental, Application, and OutlookActive Thermal Wave CloakTransformation Omnithermotics: Simultaneous Manipulation of Three Basic Modes of Heat TransferUltrabroadband Elastic Cloaking in Thin PlatesCloaking of Matter WavesFluid Flow Control with Transformation MediaTheoretical realization of an ultra-efficient thermal-energy harvesting cell made of natural materialsInvisibility cloaks for irregular particles using coordinate transformations

[1]	Pendry J B, Schurig D and Smith D R 2006 Science 312 1780
[2]	Schurig D, Mock J J, Justice B J, Cummer S A, Pendry J B, Starr A F and Smith D R 2006 Science 314 977
[3]	Schurig D, Pendry J B and Smith D R 2006 Opt. Express 14 9794
[4]	Rahm M, Schurig D, Roberts D A, Cummer S A, Smith D R and Pendry J B 2008 Photon. Nanostruct. Fundam. Appl. 6 87
[5]	Chen H Y, Chan C T and Sheng P 2010 Nat. Mater. 9 387
[6]	Liu Y M and Zhang X 2012 Nanoscale 4 5277
[7]	Milton G W, Briane M and Willis J R 2006 New J. Phys. 8 248
[8]	Chen H Y and Chan C T 2007 Appl. Phys. Lett. 91 183518
[9]	Fan C Z, Gao Y and Huang J P 2008 Appl. Phys. Lett. 92 251907
[10]	Li J Y, Gao Y and Huang J P 2010 J. Appl. Phys. 108 074504
[11]	Guenneau S and Amra C 2013 Opt. Express 21 6578
[12]	Narayana S and Sato Y 2012 Phys. Rev. Lett. 108 214303
[13]	Guenneau S, Amra C and Veynante D 2012 Opt. Express 20 8207
[14]	Schittny R, Kadic M, Guenneau S and Wegener M 2013 Phys. Rev. Lett. 110 195901
[15]	Li Y, Bai X, Yang T Z, Luo H L and Qiu C W 2018 Nat. Commun. 9 273
[16]	Yang F B, Xu L J and Huang J P 2019 ES Energy & Environ. 6 45
[17]	Xu L J and Huang J P 2020 Sci. Chin. Phys. Mech. & Astron. 63 228711
[18]	Xu L J and Huang J P 2020 Chin. Phys. Lett. 37 080502
[19]	Huang J P 2020 Theoretical Thermotics: Transformation Thermotics and Extended Theories for Thermal Metamaterials (Berlin: Springer)
[20]	Han T C, Yang P, Li Y, Lei D Y, Li B W, Hippalgaonkar K and Qiu C W 2018 Adv. Mater. 30 1804019
[21]	Li Y, Shen X Y, Wu Z H, Huang J Y, Chen Y X, Ni Y S and Huang J P 2015 Phys. Rev. Lett. 115 195503
[22]	Jiang C R, Fang C C and Shen X Y 2019 Crystals 9 438
[23]	Shen X Y, Li Y, Jiang C R, Ni Y S and Huang J P 2016 Appl. Phys. Lett. 109 031907
[24]	Yuan X B, Lin G C and Wang Y S 2016 Mod. Phys. Lett. B 30 1650256
[25]	Wang J, Dai G L and Huang J P 2020 iScience 23 101637
[26]	Xu L J and Huang J P 2020 Chin. Phys. Lett. 37 120501
[27]	Xu L J, Yang S, Dai G L and Huang J P 2020 ES Energy & Environ. 7 65
[28]	Farhat M, Guenneau S and Enoch S 2009 Phys. Rev. Lett. 103 024301
[29]	Zhang S, Genov D A, Sun C and Zhang X 2008 Phys. Rev. Lett. 100 123002
[30]	Urzhumov Y A and Smith D R 2011 Phys. Rev. Lett. 107 074501
[31]	Han T C, Zhao J, Yuan T, Lei D, Li B and Qiu C W 2013 Energy & Environ. Sci. 6 3537
[32]	You Y, Kattawar G W, Zhai P W and Yang P 2008 Opt. Express 16 6134
[33]	http://www.comsol.com/