Original Research ARTICLE
Nonlinear Heat Radiation Induces Thermal Rectifier in Asymmetric Holey Composites
- Center for Phononics and Thermal Energy Science, China-EU Joint Center for Nanophononics, Shanghai Key Laboratory of Special Artificial Microstructure Materials and Technology, School of Physics Science and Engineering, Tongji University, Shanghai, China
Asymmetric design and nonlinearity are essential to give rise to thermal rectification. In addition to the linear Fourier’s law of thermal conduction, we introduce the nonlinear thermal radiation following Stefan–Boltzmann law to asymmetric holey composites to realize the thermal rectifier. The thermal rectification results from the competition between the linear thermal conduction and the nonlinear thermal radiation, which we explain as a thermal resistor network with the linear and nonlinear components connected in series-parallel. All the results are confirmed by an analytical model that is composed of one central nonlinear thermal radiation region sandwiched by two linear conduction regions with different conductances. We can get the optimal rectification effect around the room temperature when the temperature difference is fixed, which allows the application of the designed materials to energy-saving buildings.
In analogy to the familiar electric diode in modern electronics, thermal rectifier is a device that can achieve non-reciprocal heat transport. Theoretical models for thermal rectifiers are proposed relying on various mechanisms (Roberts and Walker, 2011), including different temperature-dependent thermal properties between dissimilar materials at a contact (Stevenson et al., 1991; Dames, 2009; Kobayashi et al., 2009), asymmetric design of nanostructure (Wu and Li, 2008; Yang et al., 2009), and quantum thermal systems, such as quantum dots and other quantum heterojunctions (Scheibner et al., 2007; Ruokola et al., 2009; Wu and Segal, 2009; Ren and Zhu, 2013a). Similar concepts are even extended to the spin Seebeck diode and transistor (Ren, 2013; Ren and Zhu, 2013b; Ren et al., 2013), and the multiferroic thermal diode (Chotorlishvili et al., 2015), multiferroic switch and memory (Li et al., 2014). All these device concepts lay the foundation for the emerging field: phononics, where scientists are devoted to control heat flow and thermal information carried by phonons and beyond with electronic analogs (Li et al., 2012; Maldovan, 2013; Ren and Li, 2015).
It has been well established that asymmetric design and nonlinearity are essential to give rise to thermal rectification (Roberts and Walker, 2011; Li et al., 2012). Thus, there have been more theoretical studies of thermal rectification focusing on nonlinear systems (Wehmeyer et al., 2017). Nonlinearity can be provided in different ways to design a thermal rectifier, including one-dimensional nonlinear lattice (Terraneo et al., 2002; Li et al., 2004; Hu and Yang, 2005), nonlinear electron gas dispersion in conductors (Segal, 2008), nonlinear Schrödinger equation (Li and Ren, 2014), nonlinear near-field radiation (Otey et al., 2010; Basu and Francoeur, 2011), and even a vapor generation rectifier (Ding et al., 2017).
Punching holes in an isotropic material can induce a new heat transfer mechanism—radiation or convection (if the hole is filled with fluid or gas) in addition to the thermal conduction of the material. Thus, some interesting thermal work of materials with holes has been studied for thermal buckling analysis (Zhao et al., 2009), thermal rectifier (Hu et al., 2017) and thermal cloaking (Maldovan, 2013). Yet, few researches of materials with holes consider the heat radiation process on the surface of holes as a nonlinear procedure required in thermal rectification.
In this paper, we design a simple model of thermal rectifier, in which required nonlinearity is met by thermal radiation which is proportional to the fourth power of temperature within the framework of Stefan–Boltzmann Law. The rectification design comes from the asymmetry of the structure and the competition between nonlinear thermal radiation and linear thermal conduction.
We discuss how asymmetry and nonlinearity works to achieve thermal rectification in our model, with results of simulation and analytical explanation. In our simulation models, we consider two kinds of heat transfer mechanisms: (i) thermal conduction that is linear as we consider constant thermal conductivity and (ii) thermal radiation that is proportional to the fourth power of temperature. We report the simulation results of our holey diodes that have non-reciprocal heat transfer due to the nonlinearity of the thermal radiation and the asymmetric geometry. We obtain the clear thermal rectification, and find that there is a best average temperature, where the rectification reaches maximum, when the temperature difference is fixed. All the simulation results can be analytically explained by a simple thermal circuit model, which is composed of a central thermal radiation region sandwiched by two conduction regions with different thermal conductances.
Since there are already thermal diodes achieved with different temperature-dependent conductivity (Dames, 2009) or emissivity (Roberts and Walker, 2011), in our study we consider the bulk material itself to be homogenous and isotropic. In the homogenous material, the local heat flux density is expressed as , according to Fourier’s Law, where J is the heat flux, κ is the material’s thermal conductivity and on the vacuum interface of the holes, heat transfer follows the Stefan–Boltzmann Law , where ε is the emissivity and σ is the Stefan–Boltzmann constant. For simplification, we consider no medium inside the holes so that there is no need to consider the convection and conduction mechanism of the medium. It has been proven that (Li and Ren, 2014) combination of nonlinearity and asymmetry is an effective method to give rise to rectification. It is, therefore, reasonable to infer that holey structures could be able to transfer heat non-reciprocally and the analytical explanation will be demonstrated later.
2.1. Simulation Results for a Feasible Holey Diode Model
To achieve heat rectifier effect, we design a structure with dense periodic organization of circle holes in a pattern, only on a half of the structure itself, as shown in Figure 1A,B. This system is too complicated to be tackled analytically. Therefore, we perform two dimensional finite-element simulations, to investigate the heat transport across the structure. The temperature of the left and the right boundaries of the diode are fixed at T0 + ΔT and T0 − ΔT, respectively, where T0 represents the average temperature of the structure approximately and 2ΔT the temperature difference of two terminals. In the reverse temperature biased scenario (ΔT = −ΔT), the temperatures of the boundaries are exchanged.
Figure 1. Simulation model of holey concrete structure. (A) Sketch of the thermal diode model with periodic boundary condition in the vertical direction. The width of the structure is 83.4 cm. For the hole patterns, the hole diameter is 7.4 cm, and the horizontal distance between two centers of the adjacent holes is 12.95 cm. The gray area is the body of the concrete structure, and the blue area is the background that we consider as vacuum. The thermal conductivity of the concrete is κ = 1.7 W/(m⋅K) and emissivity is set as ε = 1. (B) Schematic graph depicting the infinite periodic holey thermal diode. We consider the model infinite in x and z axis direction. (C) Heat flux in forward (J+) and reverse (J−) temperature biased system and rectification factor R versus ΔT for thermal rectifier model in (A) and (B). T0 is fixed at 318 K. (D) Heat flux in forward and reverse temperature biased system and rectification factor versus T0. ΔT is fixed at 25 K.
If the magnitudes of thermal energy transfer in the two directions are different, we can safely say that there is a thermal rectification. Since no energy escapes from the system in the two-dimension and infinite periodic structure we proposed, we can calculate the total heat flux passing through the whole structure by integrating the normal heat flux along the boundary either the left one or the right one. We define J+ as heat flux flowing from left to right when the value of ΔT is positive in the setting of the boundary temperature, while J− with negative ΔT denotes the right-to-left heat flux in the reversed temperature biased scenario. To quantify the magnitude of thermal rectification, we define a rectification factor .
Heat flux in the forward and reverse bias as well as the corresponding rectification factor as a function of ΔT are plotted in Figure 1C. As shown in the plot, the heat transfer in the forward bias is always stronger than in the reverse bias. J+, J−, and rectification factor R are shown to increase with the increasing ΔT, in an approximately linear relationship. Moreover, it is shown in Figure 1D that although J+ and J− increase with the average temperature T0, R does not monotonically increase with T0, and reaches its maximum at about 330 K in this settings.
The thermal rectification is achieved by satisfying the required asymmetry and nonlinearity. Therefore, if the holes in our proposed model are punched exactly in the center of the structure, rectification would not appear, since the geometry of the structure is totally symmetric. To achieve a strong rectification effect, we design the model in an asymmetric way. The asymmetric setup is constructed by using two materials with different constant thermal conductivity, where material on the left holey side has thermal conductivity κL and the right complete side κR. Both the results of the simulation and the sketch of the new model are shown in Figure 2. Note the left and right materials can be identical so that κL = κR is allowed.
Figure 2. (A) Sketch of the two-material junction model, which have different thermal conductivity in left half (κL) and the right half (κR). Parameters and geometries are the same as the previous ones. (B) Schematic graph depicting the equivalent thermal circuits of the new model. Note that thermal resistance R = 1/G, where G is the corresponding thermal conductance. (C) Thermal rectification factor R versus ΔT of the model with different combination of thermal conductivity of left and right with [κL, κR] in the legend. Unit of all thermal conductivity mentioned is W/(m⋅K).
As shown in Figure 2C, we keep one of the κL and κR in the new model equal to the thermal conductivity of the single material model, 1.7 W/(m⋅K), and compare the new model with the previous one. The simulation results show that if κL = 1.7 W/(m⋅K) and κR = 1 or 0.7 W/(m⋅K), the rectification factor is improved, since not only the geometry but also the thermal properties are asymmetric in this model. However, if κR = 1.7 W/(m⋅K) and κL = 10 or 5 W/(m⋅K), even though with stronger asymmetry, the rectification effect becomes weaker.
To understand the phenomenon in a more clear way, we introduce the equivalent thermal circuit of the model in Figure 2B. GR is the conductance of the right half of the structure, which can be easily calculated by GR = κR/0.417 m, since the right half is totally homogenous. However, the left half is more complicated. GL is the equivalent conductance of the conduction mechanisms on the left. Moreover, Grad denotes the equivalent conductance contributed by the radiation of holes on the left. Due to the circle geometry of the holes, GL can not be easily expressed like GR. But, It is evident that GL ∝ κL and Grad has nonlinear positive relationship with temperature. Therefore, when κL is larger, so does GL, the resistance of conduction would be smaller guiding more heat flux flow through it rather than through radiation. Since little flux through radiation means little nonlinearity, although in this case the model has stronger asymmetry, the nonlinearity is weaker leading to the weaker diode effect. In the proposed model, thermal resistance of conduction and radiation mechanisms is in series and in parallel. In the holey side, heat can transfer in two channels, and competition between the channels weaken the required nonlinearity to give rise to rectification. To improve the rectification effects, we construct a more effective model, where thermal resistances of radiation and conduction are only in series, in the next subsection. This case also renders us an analytical explanation.
2.2. Two Infinite Plate Model with Analytical Explanation
As shown in the Figure 3, we analyze the heat transfer between two infinite parallel plates of different conductance in the simplified model. The thermal conductance instead of thermal conductivity are considered here, and can be calculated through Gi = κi/wi, where κi is the thermal conductivity and wi represents the thickness of two infinite plates, with i = 1, 2. In this case, the combined processes of thermal conduction and thermal radiation form an equivalent thermal circuit in series. Since all the heat flux leaving one plate will arrive at the other plate, the heat flux through the three regions (I, II, III) would be the same:
Figure 3. (A) Sketch of the two infinite plate model. This is a one-dimensional model with three separated region. Temperature of two outer surfaces is set T0 + ΔT and T0 − ΔT, respectively. w1 and w2 denotes the thickness of two infinite plate. ε is the emissivity of the inner surface. T1 and T2 are measured as temperature of two inner surface. (B) Schematic graph depicting the equivalent thermal circuits of the plate-vaccum-plate sandwich model. Thermal resistance of conductance of region I and region III, and the equivalent thermal resistance of radiation in region II are in series. (C) w1 versus rectification factor R. T0 and ΔT are fixed at 318 and 25 K. (D) G2 versus rectification factor R in the two-part model. T0 and ΔT are fixed at 318 and 25 K.
Note the Stefan-Boltzmann constant σ is at −8 order of magnitude in SI units, which makes the nonlinear radiation contribution rather weak. If we consider the situation where ΔT as well as J is small, we can obtain the expressions of J+ and J− up to the second order of thermal bias, as following:
The first term of equation (2) is a linear contribution that has a clear physical picture: with linear response, region II is approximated as a homogenous material with conductance , so that the equivalent conductance of the model is composed of thermal resistances of three regions in series: , with the temperature bias 2ΔT. For a temperature difference 2ΔT smaller than T0, the second order term is small compared with the first term. Therefore, the heat flux J+ and J− is approximately proportional to the temperature difference ΔT, which also explain the nearly linear relationship of J+ and J− versus ΔT shown in the Figure 1C. Equation (2) also shows that the dominant first term of heat flux, although linear respect to thermal bias, has a nonlinear positive relationship with T0, which well explains the nonlinear positive relationship and gradually saturated dependence on T0, revealed in Figure 1D.
The second term in the heat flux expression equation (2) is very important to heat rectification. From the expression of heat flux, we can obtain the rectification factor up to the first order of thermal bias:
Therefore, the rectification coefficient is also linear to the ΔT as showed in Figure 1C and Figure 2C, which will vanish if the system becomes symmetric G1 = G2. Through the analytical result, we can understand the nontrivial relationship between R versus T0, as shown in Figure 1D. Clearly from equation (2), when T0 increases from low temperature, R increases while when T0 is already at high temperature, R decreases as T0 increases. With simple differential method, we can calculate that the denominator of the R reaches its minimum value when , where R has a maximum value when changing the average temperature T0. This behavior coincides with what one observed in Figure 1D.
We also study the dependence of R on the width w1 in the plate–vacuum–plate sandwich model, to see the effect of competition between linear conduction and nonlinear radiation on the rectification. In Figure 3C, we keep w2 fixed at 0.3 m, and both κ1 and κ2 are equal to 1.7 W/(m⋅K). By changing w1 only, we can see the importance of asymmetry in giving rise to the rectification that can be reversed. There is an exact vanishing point of R at w1 = w2, which means G1 = G2 in this case. If we increase w1 from the initial value 0.3 m, |R| first increases with w1, since the asymmetry of the system increases. However, while increasing w1 further, the left conductance G1 gets smaller, so the temperature difference at two inner surfaces (T1–T2) gets smaller. Note that T1–T2 determines the flux (so as the total flux) through radiation within the nonlinear framework of Stefan–Boltzmann’s Law, so when w1 goes even larger, the radiative flux of nonlinear nature goes smaller so that the rectification R falls down. When we decrease w1 from the initial value 0.3 m, R increases quickly, and reaches a maximum at about 12%. This is reasonable since when w1 becomes smaller and smaller, approaching to 0, G1 becomes larger and gets close to infinity. The strongest asymmetry is achieved in this way, so that R can reach its maximum limit in this case when G1 → ∞.
To further maximize the rectification factor, we may consider a two-segment model with only region II and III where the left plate region I is removed. Temperature at the left boundary is imposed at the border of the vacuum part (region II). In Figure 3D, we show the numerical rectification factor R vs G2 with fixed T0 and ΔT in the two-segment model. As we discussed above, the thermal rectification results from the competition between nonlinear radiation and linear conduction. When there is only one mechanism—either radiation or conduction, rectification can not appear since the system is symmetrical even with nonlinearity. Only if G2 and Grad are comparable, which means the nonlinear and linear mechanisms are comparable to make the system most asymmetric, rectification effect becomes strongest. In other word, R reaches its maximum around G2 ≈ Grad. As shown in the case of Figure 3D, , and R indeed reaches the maximum when G2 is about 7, comparable with Grad.
3. Conclusion and Discussion
We have designed the asymmetric holey composites as a thermal rectifier, which has non-reciprocal heat transfer. The non-reciprocal heat transfer is attributed to the thermal radiation of nonlinear nature and the asymmetric geometric of the system. We have simulated and then confirmed from the theoretical analysis that the thermal rectification is approximately linear to the temperature difference of the system with the average temperature fixed and the maximal thermal rectification can be obtained by tuning the average temperature when the temperature difference is fixed. There have been considerable interests in realizing thermal rectifiers for its ability in thermal management and control for heat energy in the recent years. Most researches focus on achieving thermal rectifiers in nanoscale or even smaller, for its numerous potential applications in micro/nano-electronic cooling, thermal memory and computations. We have designed our thermal rectification model at macro-scale by inducing the nonlinearity with holes in the design, which is a simple model for thermal rectifiers. This means that the setup can be easily built, which might have potential applications in materials of architecture for diverse thermal management problems, e.g., energy-saving buildings, solar cells, and efficient refrigerator.
Furthermore, the rectification strength of the holely diodes can be further improved by constructing with two materials that with different temperature-dependent thermal conductivity of opposite trends [similar to the one in Dames (2009)]. In some way, the crucial nonlinearity induces the effective temperature-dependent conductance. For the nonlinear radiation, the conductance , which under asymmetric geometry and bias reversing will effectively change the values. Therefore, a rectification action emerges. Beyond the nonlinear thermal radiation, the nonlinearity can also be introduced by the ballistic transport at quantum regime, since the phonon population is described by the Bose–Einstein distribution that has nonlinear dependence on temperature. Explicitly, the ballistic quantum thermal transport is described by the well-known Landauer-Büttiker formula:
with the phononic transmission coefficient, NL(ω) − NR(ω) the phonon population difference between two sides. When combining this ballistic transport with diffusive transport described by the linear Fourier’s law and with asymmetric geometry of the system, the rectification emerges due to the nonlinear Bose–Einstein distribution. However, the rectification disappears in the high temperature classical regime since . As such, equation (4) reduces to the one with linear temperature dependence: . Therefore, it is readily to prove that when combing with the Fourier’s law, even with asymmetric geometry, the system can not have rectification. Both nonlinearity and asymmetry are necessary to construct the thermal rectifier (Li and Ren, 2014).
GW and WZ carried out the simulation, performed the theoretical analysis, and contributed equally to this work. GW, WZ, HC, and JR, all contributed to discussion and the writing. JR conceived the project.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer, QC, and handling editor declared their shared affiliation.
This work is supported by the National Natural Science Foundation of China (Grant Nos. 11775159 and 61621001), the National Key Research Program of China (No. 2016YFA0301101), and the National Youth 1000 Talents Program in China.
Ding, H., Peng, G., Mo, S., Ma, D., Sharshir, S. W., and Yang, N. (2017). Ultra-fast vapor generation by a graphene nano-ratchet: a theoretical and simulation study. Nanoscale 9, 19066–19072. doi:10.1039/C7NR05304E
Li, N., Ren, J., Wang, L., Zhang, G., Hanggi, P., and Li, B. (2012). Phononics: manipulating heat flow with electronic analogs and beyond. Rev. Mod. Phys. 84, 1045–1066. doi:10.1103/RevModPhys.84.1045
Ren, J., Fransson, J., and Zhu, J. X. (2013). Nanoscale spin Seebeck rectifier: controlling thermal spin transport across insulating magnetic junctions with localized spin. Phys. Rev. B 89, 106–114. doi:10.1103/PhysRevB.89.214407
Ren, J., and Zhu, J. X. (2013a). Heat diode effect and negative differential thermal conductance across nanoscale metal-dielectric interfaces. Phys. Rev. B Condens. Matter 87, 291–295. doi:10.1103/PhysRevB.87.241412
Ren, J., and Zhu, J. X. (2013b). Theory of asymmetric and negative differential magnon tunneling under temperature bias: towards a spin Seebeck diode and transistor. Phys. Rev. B Condens. Matter 88, 3169–3182. doi:10.1103/PhysRevB.88.094427
Wehmeyer, G., Yabuki, T., Monachon, C., Wu, J., and Dames, C. (2017). Thermal diodes, regulators, and switches: physical mechanisms and potential applications. Appl. Phys. Rev. 4, 041304. doi:10.1063/1.5001072
Keywords: thermal rectification, thermal radiation, non-reciprocal heat transfer, metamaterial, heat conduction, thermal diode, energy conservation
Citation: Zhu W, Wu G, Chen H and Ren J (2018) Nonlinear Heat Radiation Induces Thermal Rectifier in Asymmetric Holey Composites. Front. Energy Res. 6:9. doi: 10.3389/fenrg.2018.00009
Received: 22 December 2017; Accepted: 14 February 2018;
Published: 05 March 2018
Edited by:Nuo Yang, Huazhong University of Science and Technology, China
Reviewed by:Qiang Cheng, Huazhong University of Science and Technology, China
Dahai He, Xiamen University, China
Copyright: © 2018 Zhu, Wu, Chen and Ren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Ren, firstname.lastname@example.org
†These authors have contributed equally to this work.