Next Article in Journal
High-Sensitivity in Dielectrophoresis Separations
Next Article in Special Issue
Resistance Reduction of Conductive Patterns Printed on Textile by Curing Shrinkage of Passivation Layers
Previous Article in Journal
Editorial for the Special Issue on Point-of-Care Devices
Previous Article in Special Issue
Recent Developments of Flexible and Stretchable Electrochemical Biosensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Recent Advances on Thermal Management of Flexible Inorganic Electronics

1
Institute of Solid Mechanics, Beihang University (BUAA), Beijing 100191, China
2
State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China
3
Institute of Structure and Environment Engineering, Beijing 100076, China
4
Department of Engineering Mechanics, Soft Matter Research Center, and Key Laboratory of Soft Machines and Smart Devices of Zhejiang Province, Zhejiang University, Hangzhou 310027, China
*
Author to whom correspondence should be addressed.
Micromachines 2020, 11(4), 390; https://doi.org/10.3390/mi11040390
Submission received: 28 February 2020 / Revised: 27 March 2020 / Accepted: 5 April 2020 / Published: 9 April 2020

Abstract

:
Flexible inorganic electronic devices (FIEDs) consisting of functional inorganic components on a soft polymer substrate have enabled many novel applications such as epidermal electronics and wearable electronics, which cannot be realized through conventional rigid electronics. The low thermal dissipation capacity of the soft polymer substrate of FIEDs demands proper thermal management to reduce the undesired thermal influences. The biointegrated applications of FIEDs pose even more stringent requirements on thermal management due to the sensitive nature of biological tissues to temperature. In this review, we take microscale inorganic light-emitting diodes (μ-ILEDs) as an example of functional components to summarize the recent advances on thermal management of FIEDs including thermal analysis, thermo-mechanical analysis and thermal designs of FIEDs with and without biological tissues. These results are very helpful to understand the underlying heat transfer mechanism and provide design guidelines to optimize FIEDs in practical applications.

1. Introduction

Flexible inorganic electronics, which offer the high performance of conventional rigid electronics but with the ability to be deformed like rubber, have enabled many novel devices, with some examples illustrated in Figure 1, including flexible microscale inorganic light-emitting diodes (μ-ILEDs) [1], epidermal electronics [2], flexible and stretchable 3ω sensors [3], skin-mounted power management systems [4], and soft robotics [5]. The challenge to develop flexible inorganic electronic devices (FIEDs) lies in the mismatch between the intrinsic brittle inorganic electronic material (e.g., silicon) and the flexibility requirement in applications. One effective strategy to overcome this mismatch is to integrate the functional inorganic electronic material in a stretchable format [6,7,8,9,10,11,12] at strategic locations on a soft polymer substrate by advanced transfer printing techniques [13,14,15,16]. Due to the unique mechanical properties of FIEDs, they can have conformal contacts with biological tissues under complex deformations, which enable the real-time monitoring of human vital signs for early diagnosis [17]. In the past decade, researchers have demonstrated several flexible real-time health monitoring devices to measure human vital signs including body temperature [4,18], blood glucose [19,20,21], blood oxygen [22], skin hydration [3], blood flow [2], brain electrophysiological signal [23,24], electrocardiography (ECG) [25], etc.
Microscale inorganic light-emitting diodes, as a typical inorganic electronic device, show excellent interactions with biological tissues as excitation sources. For example, μ-ILEDs can be an alternative excitation source to an optoelectronic tweezer via noncontact manipulation methods to investigate behaviors of various biological samples such as cells [26,27,28,29,30,31,32], proteins [33], DNA molecules [34,35] and particles [36,37]. In addition to the above applications, μ-ILEDs can also be adopted in optogenetics to control, affect and readout the neural activities of biological creatures, especially for stimulation in brain [38,39,40,41], retina [42,43,44], and audition [45,46,47,48,49]. However, a large amount of heat energy will be generated by μ-ILEDs when they are in service.
The poor heat dissipation capacity of soft polymer substrate with a thermal conductivity on the order of 0.1 W/m/K in FIEDs, which is much lower than that of the conventional silicon substrate (~100 W/m/K), may induce a high temperature rise in the functional component and cause device performance degradation [50]. When FIEDs are integrated with biological tissues, the sensitive nature of tissues to temperature poses more stringent requirements on the temperature control in practical applications. Moreover, the thermo-mechanical coupling due to the high temperature rise could affect the device reliability [51,52]. Therefore, thermal management of FIEDs becomes one of the most important challenges in the applications of FIEDs, especially in those involving biological tissues. Thermal models will help to establish design guidelines to reduce the adverse thermal effects [53]. This paper aims to review the recent advances in thermal management of FIEDs with the focus on thermal analysis, thermo-mechanical analysis and thermal designs of FIEDs with and without biological tissues, through discussions of analytical models, finite element analysis and experiments.

2. Thermal Management of FIEDs

2.1. Thermal Analysis of FIEDs

We take micro-scale inorganic light-emitting diodes (μ-ILEDs), which serve as the functional components of FIEDs, to induce the undesired heating, as an example to review the recent advances on thermal analysis of FIEDs. Kim et al. [1] reported an unusual strategy using anisotropic etching and microscale device assembly to fabricate μ-ILEDs on unconventional substrates. Figure 2a shows the schematic diagram of the μ-ILED structure with the μ-ILED encapsulated by benzocyclobutene (BCB) and metal layers (i.e., interconnection) on a substrate. Due to the symmetric conditions, a quarter geometry of the μ-ILED layout with key geometrical dimensions labelled is shown in Figure 2b. Cui et al. [54] established a three-dimensional analytical model to study the thermal properties of rectangular μ-ILEDs with the length of 2a × 2b and a thickness of HLED. The steady state heat conduction equation can be expressed as
2 Δ T x 2 + 2 Δ T y 2 + 2 Δ T z 2 = 0
where Δ T = T ( x , y , z ) T 0 is the temperature increase from the ambient temperature T 0 . The μ-ILED is modeled as a planar heat source since its thickness is much smaller than its in-plane dimension. The boundary conditions include the heat convection on the top surface, the constant temperature on the bottom surface, and the heat generation in the μ-ILED region. The heat flux and temperature are continuous at the interfaces. The temperature increase of μ-ILED is obtained as
Δ T LED = 2 P a b [ 2 a b + ( a + b ) H LED ] π 2 0 + 0 + ( E + F ) sin 2 ( α a ) sin 2 ( β b ) G α 2 β 2 α 2 + β 2 d α d β ( ξ )
where P is the heat generation power; E, F, G are
{ E = ( 1 + k m k B ) k m α 2 + β 2 + h k m α 2 + β 2 h e 2 ( H m + H B ) α 2 + β 2 + ( 1 k m k B ) e 2 H B α 2 + β 2 F = ( 1 k m k B ) k m α 2 + β 2 + h k m α 2 + β 2 h e 2 H m α 2 + β 2 + ( 1 + k m k B ) G = k B ( E F ) + k s ( E + F ) coth ( H s α 2 + β 2 )
where k, H and h are the thermal conductivity, thickness and coefficient of heat convection, respectively; the subscripts “m”, “B”, and “s” stand for interconnection, BCB and substrate, respectively.
To clearly understand the effects of the parameters (such as geometrical parameters, material parameters and loading parameters) on the μ-ILED temperature, Equation (2) was further simplified based on a few reasonable assumptions. First, the effect of the heat convection boundary on the top surface is negligible. Second, the substrate is regarded as a semi-infinite solid because the substrate thickness is much larger than other dimensions. Third, e 2 H m α 2 + β 2 2 H m α 2 + β 2 , and e 2 H B α 2 + β 2 2 H B α 2 + β 2 . After adopting the above assumptions, Equation (2) becomes
Δ T LED 2 P a b [ 2 a b + ( a + b ) H LED ] π 2 0 + 0 + sin 2 ( α a ) sin 2 ( β b ) α 2 β 2 α 2 + β 2 k m H m α 2 + β 2 1 + k m k B H m H B ( α 2 + β 2 ) + k s d α d β
Figure 2c compares the temperature increase of the μ-ILED with the accurate analytical prediction from Equation (2), the approximate prediction from Equation (4), and finite element analysis (FEA) for μ-ILEDs with 100 μm × 2b × 5 μm under the heat generation power of 12 mW. The good agreement among them validates the treatments in analytical models. The temperature increase on the top surface of metal layer is obtained as
Δ T surface = 4 P π 2 a b 0 + 0 + k m H m α 2 + β 2 e ( H m + H B ) α 2 + β 2 k m H m α 2 + β 2 h sin ( α a ) sin ( β b ) sin ( α x ) sin ( β y ) G α β α 2 + β 2 d α d β
Again, the good agreement with an error less than 5% of the surface temperature between the analytical prediction from Equation (5) and FEA with power of 12 mW and μ-ILED size of 100 × 200 × 5 μm in Figure 2d validates the analytical model.
The results for a single μ-ILED in Equations (2)–(5) can be adopted to find the temperature of μ-ILED arrays based on the superposition method. Cui et al. [54] took an array of two μ-ILEDs with one’s dimension of 100 × 100 µm (2a × 2b), the other’s dimension of 100 × 200 µm (2a × 2b1) and the spacing of 200 µm (b0), as shown in Figure 2e, as an example to illustrate this approach. The temperature increase at any point (x, y, z) for the array of two μ-ILEDs is given by
Δ T array ( x , y , z ) = Δ T 1 ( x , y , z ) + Δ T 2 ( x , y , z )
As shown in Figure 2f, the analytical prediction of surface temperature agrees well with FEA, which validates the accuracy of the analytical model. Besides the surface temperature, the temperature increase of μ-ILED array can also be determined with any spacing for an array with finite or even infinite number of μ-ILEDs.
In addition to the 3D thermal model for rectangular heat sources, Lu et al. [55] established an analytical axisymmetric model to study the thermal properties of square μ-ILEDs with an effective radius of r 0 = L / π , with L as the length of the μ-ILED, which compares well with FEA and experiments. Yin et al. [56] studied the thermal properties of different shaped serpentine metal heater devices analytically, which can be used to evaluate the heating effects of complex shaped flexible heaters.
Pulsed operation is an effective routine for thermal management by further reducing the adverse thermal influences. Kim et al. [50] fabricated μ-ILEDs, which are located on top of a PI layer and hydrogel substrate and encapsulated with an SU8 layer (Figure 3a) to demonstrate the pulsed strategy. The key geometrical dimensions are labelled in Figure 3b, which shows a quarter geometry of the μ-ILED structure. The applied heat generation power (Figure 3c) is denoted by Q ( t ) = Q 0 U ( t ) , where Q0 is the peak power and U(t) as a unit pulsed power. Cui et al. [57] established an analytical model to study the thermal properties of μ-ILEDs in pulsed operation. The duty cycle is defined by D = τ / t 0 , where τ is the pulse duration and t0 is the period.
Figure 3d compares the μ-ILED temperature increase after saturation under pulsed power with the peak power of 30 mW, the period of 50 ms, the duty cycle of 50% and the μ-ILED dimension of 100 × 200 × 6.5 µm (2a × 2b × hLED) from analytical prediction by Equation (9) and FEA. The good agreement validates the analytical model. In order to study the effect of duty cycle, the maximum and minimum temperature increases of μ-ILED are shown in Figure 3e, which indicates that the duty cycle has a significant effect on the temperature of a μ-ILED. A smaller duty cycle will yield a lower temperature increase of the μ-ILED. Furthermore, Cui et al. [57] established a scaling law for the maximum temperature increase of μ-ILED as
Δ T ¯ L E D max = ( 1 + γ 2 ) 3 γ 2 ( 1 + 1 + γ 2 γ h ¯ L E D ) G ( χ , h ¯ S U 8 , D , γ , σ )
where the normalized maximum temperature Δ T ¯ L E D max = k S U 8 a Δ T L E D max / Q 0 , χ = k h y d r o g e l c h y d r o g e l ρ h y d r o g e l / ( k S U 8 c S U 8 ρ S U 8 ) , γ = b / a , h ¯ L E D = h L E D / a , h ¯ S U 8 = h S U 8 / a , σ = a 2 / ( λ S U 8 t 0 ) , and G is a non-dimensional function [57]. Figure 3f shows the normalized maximum temperature as a function of duty cycle under various values of σ with χ = 5.55 , γ = 2 , h ¯ L E D = 0.13 , and h ¯ S U 8 = 0.14 . It is noted that small D or large σ are effective to reduce the maximum temperature increase of the μ-ILED. More discussions on the influences of other non-dimensional parameters can be found in Cui et al.’s work [57].
On the other side, many researchers theoretically and experimentally investigated the heat exchange in polymers, which aims to replace metallic heat exchanges for the μ-ILED devices [58,59,60]. A series of new composites were fabricated to enhance the thermal properties of polymers by adding particles [61,62,63], fibres [64,65], carbon nanotubes [66,67,68] and so on.

2.2. Thermo-Mechanical Analysis of FIEDs

In addition to effects of temperature increase, the heating induced thermal stress (or strain) is not negligible since the thermal expansion coefficient of sot polymer substrate is much larger than that of the stiff functional inorganic component, which results in noticeable mismatch at the interface between the stiff functional component and the substrate. This mismatch may lead to interfacial delamination, thus causing the failure of FIEDs. Moreover, the non-negligible thermal stress may cause mechanical discomfort when FIEDs are integrated with biological tissues. Therefore, thermo-mechanical analysis of thermal management of FIEDs should be fully explored.
Figure 4a shows a typical structure of an FIED with a functional component encapsulated by an encapsulation layer on a soft polymer substrate. Zhang et al. [69] developed a three-dimensional steady state analytical model to study thermo-mechanical properties of FIEDs based on the transfer matrix method through the Fourier integral transform in space domain and the Laplace integral transform in time domain. According to the linear thermoelasticity theory, the thermo-mechanical equilibrium equations can be expressed as
{ 2 u + 1 1 2 μ e x 2 α ( 1 + μ ) 1 2 μ θ x = 0 2 v + 1 1 2 μ e y 2 α ( 1 + μ ) 1 2 μ θ y = 0 2 w + 1 1 2 μ e z 2 α ( 1 + μ ) 1 2 μ θ z = 0 2 θ x 2 + 2 θ y 2 + 2 θ z 2 = ρ c k θ t
where θ, ρ, c, t are temperature increase, density, heat capacity and time, respectively; u, v, w are the displacements along x, y and z directions, respectively; and e = u / x + v / y + w / z is the volumetric strain. After solving the above equilibrium equations combined with boundary and continuous conditions, the stresses, displacements, temperatures and heat fluxes at any point in transformation space are related the counterparts on the bottom surface and can be obtained by
{ σ ˜ z ( ξ , η , z , s ) τ ˜ z x ( ξ , η , z , s ) τ ˜ y z ( ξ , η , z , s ) u ˜ z ( ξ , η , z , s ) v ˜ z ( ξ , η , z , s ) w ˜ ( ξ , η , z , s ) θ ˜ ( ξ , η , z , s ) Q ˜ ( ξ , η , z , s ) } = [ Λ ] { σ ˜ z ( ξ , η , 0 , s ) τ ˜ z x ( ξ , η , 0 , s ) τ ˜ y z ( ξ , η , 0 , s ) u ˜ z ( ξ , η , 0 , s ) v ˜ z ( ξ , η , 0 , s ) w ˜ ( ξ , η , 0 , s ) θ ˜ ( ξ , η , 0 , s ) Q ˜ ( ξ , η , 0 , s ) }
where σ z , τ z x , τ y z , Q are the normal stress, shear stresses and heat flux, respectively; [ Λ ] is the transfer matrix; the superscript “∼” stands for the variables in the transformation space. The inverse transformation of Equation (9) gives the stresses, displacements and temperature.
In order to validate the analytical model, an experimental setup was established to acquire the strains on the top surface of the device via digital image correlation technology and the temperature via an infrared camera. The FIED with a heater embedded in polydimethylsiloxane (PDMS) is taken as an example to illustrate the analytical approach. Figure 4b compares the temperature increase at the center point on the top surface for a heater (2 × 1.25 mm) under a heat generation power of 291 mW. It is obvious that the temperature rise increases quickly at first and then slowly to the steady state. Figure 4c–e compares the strain contours on the top surface. The good agreement among the analytical predictions, FEA and experiments validates the accuracy of the analytical thermo-mechanical model.

3. Thermal Management of FIEDs Integrated with Biological Tissues

3.1. Thermal Analysis of FIEDs Integrated with Biological Tissues

The unique mechanical properties of FIEDs, which can sustain complex bending, twisting or stretching even with large deformations, enable their conformal contact with biological tissues for diagnostic and surgical functions. The biological tissues pose great challenges in thermal management. First, the biological tissue is very sensitive to temperature and even a few degrees in temperature increase may cause human discomfort or tissue lesion. Second, the heat transfer in biological tissue doesn’t obey Fourier’s law due to its strong biological features including blood perfusion, metabolism, and so on.
In the past decade, several analytical studies were performed to study the thermal properties of FIEDs integrated with biological tissues, especially human skin. Cui et al. [70] established a one-dimensional analytical model to predict the temperature in FIED/skin systems. In the model, the skin consists of four layers: stratum corneum layer, epidermis layer, dermis layer and fat layer, as shown in Figure 5a. It is known that blood perfusion only exists in the dermis layer and metabolism exists in all the layers of human skin. The heat transfer equations, including Fourier heat transfer equation in FIEDs and the Pennes bio-heat transfer equation in human skin tissues are given by
{ k i d 2 Δ T i d z 2 = 0 ( i = 1 , 2 ) k i d 2 Δ T i d z 2 + q met = 0 ( i = 3 , 4 , 6 ) k 5 d 2 Δ T 5 d z 2 ω b ρ b c b Δ T 5 + q met = 0
where the subscript i denotes the corresponding layer, k is the thermal conductivity, Δ T is the temperature increase from the ambient temperature; ωb, ρb and cb are the flow rate, density and heat capacity of blood, respectively, and qmet is the metabolic heat generation. Considering the boundary conditions with natural convection on the top surface and constant temperature on the bottom surface, and continuity conditions (Figure 5b), the temperature increase in the system can be derived analytically by solving Equation (10). Figure 5c shows the temperature increase in the system along the thickness direction. The results from the analytical model agree well with FEA. It is shown that the temperature increase varies nonlinearly in the dermis layer due to the effect of blood perfusion while temperature varies linearly in other layers.
Cui et al. [71] further extended the one-dimensional thermal model to a three-dimensional model to investigate the thermal properties of FIEDs involving a rectangular heating component (e.g., μ-ILED) integrated with human skin. The analytical model accounts for the coupling between the Fourier heat conduction in μ-ILED and the Pennes bio-heat transfer in skin. Figure 5d shows the schematic structure of the FIED/skin system with key dimensions labelled. The governing equations and boundary conditions are all non-homogeneous for such a system. The linear superposition principal is adopted to obtain the temperature increase, which is validated by FEA. Based on the thresholds of 39 and 43 °C for thermal pain and tissue injury, a thermal comfort map is obtained as shown in Figure 5e under the heat generation power of 5 mW, with the comfortable region, uncomfortable region and tissue injury region separated by the two white lines. It is shown that a thicker substrate with a higher thermal conductivity can reduce the maximum temperature increase to avoid human discomfort or tissue injury. Cui et al. [71] also investigated the thermal properties of FIED/skin system under a pulsed heat generation power to further reduce the temperature rise.
In order to consider the actual situations, the clothing effect [72], sweating effect [73], and interfacial thermal resistance effect [74] on thermal properties of the FIED/skin system were also investigated. Besides human skin, thermal properties of FIED integrated with other biological tissues were also studied. For example, Yin et al. [75] developed an axisymmetric cambered analytical model in a sphere coordinate system to investigate the thermal characteristics of FIEDs affixed on a large curvature myocardial surface via steady state analysis. The model, verified by FEA, can exactly predict the thermal properties of the heart or other large curvature structures. Effects of parameters, such as curvature, myocardium, substrate thickness and substrate material on the thermal properties of FIEDs were investigated systematically. Li et al. [76] established an analytical model to predict temperatures of four μ-ILEDs in the brain and yielded optimal designs to constrain the temperature rise below 1 °C [77].
Thermal models of FIED/skin systems can be used to characterize human vital signs and mechanical properties of human skin. For example, Webb et al. [2] developed an ultra-thin, soft, conformal temperature sensor and heater array that can measure blood flow. When the power is applied, the heaters induce an anisotropic temperature distribution, which can give the blood flow based on analytical thermal models. Tian et al. [3] proposed an effective three omega technology to obtain the skin hydration by comparing the analytical model, finite element analysis and experiments for the FIED/skin system under a transient periodical load. Li et al. [78] designed and fabricated flexible sensors on fingers/toenails to precisely study thermal transport characteristics of nail bed tissue.

3.2. Thermo-Mechanical Analysis of FIEDs Integrated with Biological Tissues

Li et al. [79] established an axisymmetric steady state thermo-mechanical model to obtain the temperature and stresses in the FIED/skin system. Figure 6a schematically shows the multilayer structure of FIED/skin system with the power density of Q0 loaded in a circular region with radius of r0. The thicknesses of substrate, epidermis, dermis and fat layers are denoted by hP, hE, hD and hF, respectively. The transfer matrix method was adopted to derive the temperature increase and thermal stresses in the system at any point. For example, Figure 6b shows the maximum principal stress distribution at the device/skin interface for a heater size of 2 mm. The analytical predictions agree well with FEA, which validates the accuracy of the analytical thermo-mechanical model. The maximum principal stress reaches an almost unchanged maximum value within the heating region (r < 2 mm) and then decreases quickly as the distance to the heater increases. The distribution of temperature increase at the device/skin interface shows a similar trend as that of the maximum principal stress [79].
The effects of size of the heating component and substrate thickness on the maximum temperature increase and the maximum principal stress are shown in Figure 6c,d, respectively. The black line in Figure 6c gives the critical design parameters, which yield the maximum temperature increase of 6 °C, to separate the thermal comfortable and uncomfortable regions. The region above the black line is the thermal comfortable region while the one below is the thermal uncomfortable region. The black and white lines in Figure 6d give the design parameters, which yield the principal stress with the maximum absolute value of 20 kPa, to separate the mechanical comfortable and uncomfortable regions. The region above the black and white lines is the mechanical comfortable region and the ones below are mechanical uncomfortable regions. These results are very helpful to provide design guidelines to reduce the undesired thermal responses in the FIED/skin system.
To further investigate the heater shape effect on the thermo-mechanical properties of the FIED/skin system, Zhang et al. [80] established a three-dimensional analytical steady state model to study the thermo-mechanical behaviors of the FIED/skin system with rectangular heating components based on the transfer matrix method. The transfer equation for each layer is derived analytically. Both cases of one heater and multiple heaters were systematically investigated.
Besides the steady state thermo-mechanical analysis, a transient thermo-mechanical model is also important considering that the perceptions of thermoreceptors and mechanoreceptors, as uncomfortable feelings are related to the response of thermoreceptors, the rate of thermal strain, and the thermal stress with thresholds of 19.8 imp/s, 0.21%/s and 20 kPa, respectively [81]. Zhang et al. [81] established a three-dimensional transient analytical model to predict the temperature increase and thermal stresses in the FIED/skin system with the structure layout shown in Figure 7a,b. Both cases, under a constant heat generation power and a pulsed power (Figure 7c), were systematically studied. Figure 7d–f shows some typical results under a pulsed heat generation power for the FIED/skin system with a heating component size of 1 mm × 1 mm and a thickness of PDMS encapsulation of 0.6 mm. The pulsed power density used for calculations has a magnitude of 1mW/mm2, pulse duration time of 0.5 s and period of 1 s. The maximum principal strain and strain rate are given in Figure 7d,e, respectively. The thermoreceptors and mechanoreceptors in the human body help humans feel temperature, strain or stress. Once the responses of these receptors exceed the thresholds, humans may have an uncomfortable feeling. The red dotted line in Figure 7e corresponds to the response threshold of strain rate (0.21%). It is noted that the maximum principal strain rate is below 0.21%, which indicates that the strain rate is not large enough to induce mechanical discomfort under these specific loading conditions. Figure 7f shows the influences of thicknesses of encapsulation and substrate on the peak value of the maximum rate of the principal strain. These results are very useful to determine whether the strain rate can induce mechanical discomfort. In addition to the maximum rate of the principal strain, the maximum thermoreceptor response and the maximum principal stress were also systematically investigated to study whether the temperature and thermal stress induce human discomfort in Zhang et al.’s work [81].

4. Thermal Designs of FIEDs

In practical applications of FIEDs integrated with biological tissues, it is important to prevent the undesired heating induced by FIEDs to dissipate to biological tissues through the substrate of FIEDs. Several active thermal designs [82,83,84,85,86] have been reported to help thermal management of the FIED/skin system. This section will briefly review the recently developed orthotropic substrate design [83] and the functional soft composite design [85].

4.1. Orthotropic Substrate Design

In order to prevent heating energy to dissipate from FIEDs to human skin, Li et al. [83] reported an orthotropic substrate design consisting of layered metal and polymer materials shown in Figure 8a, which can control the heat flow direction to achieve the goal of thermal management. The in-plane and off-plane thermal conductivities can be obtained as
k in - plane = k 1 + k 2 2 , k off - plane = 2 k 1 k 2 k 1 + k 2
where k1 and k2 are the thermal conductivities of metal and polymer materials, respectively. For a large difference in k1 and k2, Equation (11) yields k in - plane k off - plane , which indicates that the heat dissipation along the in-plane directions is much larger than that along the off-plane direction. This unique orthotropic feature offers an advantage in thermal management of FIEDs, especially in bio-integrated applications where the heat dissipation along the off-plane direction should be minimized.
The μ-ILED is taken as the heating component in FIEDs to demonstrate the performance of orthotropic substrate design and to identify the underlying heat transfer mechanism. Li et al. [83] established a three-dimensional analytical model to study the thermal properties of FIEDs with an orthotropic substrate. Figure 8b shows the influence of the orthotropic substrate design on the temperature of the heating component. The analytical prediction agrees well with FEA. The μ-ILED decreases quickly as the in-plane thermal conductivity of orthotropic substrate increases. For the thermal conductivity ratio k in - plane / k off - plane of 60, the μ-ILED temperature has a significant decrease from 21.5 to 3 °C. Figure 8c,d shows the temperature distributions on the top surface and cross-sectional plane for the cases of homogeneous PDMS substrate and orthotropic PDMS/aluminum substrate, respectively. Indeed, the orthotropic substrate changes the heat dissipation path with much more heat dissipating along the in-plane directions, as indicated by the black lines and red arrows. Li et al. [84] further studied the thermal properties of FIEDs on an orthotropic substrate for biointegrated applications with influences of thermal conductivity of orthotropic substrate, substrate thickness, and loading parameters fully investigated. It was shown that the orthotropic substrate design can reduce the temperature rise of skin significantly.

4.2. Functional Soft Composite Design

Shi et al. [85] proposed a functional soft composite to serve as the thermal protecting substrate for wearable electronics, with abilities to manipulate the heat flow and efficiently absorb the excessive heat energy. The thermal protecting substrate features embedded phase change materials with a thin metal film on the top in a soft polymer as shown in Figure 9a. The thin metal film helps to spread the heat energy along in-plane directions while the phase change material helps to absorb the heat energy once the temperature exceeds the transition temperature. A wearable electronics with soft heater as the heating component, copper as the metal film, paraffin as the phase change material and PDMS as the substrate material was developed as shown in Figure 9b. Figure 9c shows the optical image of a wearable device under bending deformation, which indicates its ability to have a conformal contact with skin.
Figure 9d compares the temperature contour on the bottom surface of the substrate between wearable electronic devices on a conventional PDMS substrate and a thermal protecting substrate from experiments, which can obviously show the high efficiency of the functional thermal protecting substrate in thermal management. In order to reveal the underlying physics associated with the thermal protection, a three-dimensional finite element model was established. Figure 9e shows the temperature distributions on the bottom surface of FEA, which agree well with experiments. Figure 9f compares the maximum temperature on the bottom surface with a heating duration time of 40 s. The good agreement between experiments and FEA again validates FEA. Figure 9g shows the temperature versus time to reveal the thermal protection mechanism. For wearable electronic devices with a conventional substrate, the temperature increases monotonically to steady state, while the temperature response can be described as four stages for wearable electronic devices with a thermal protecting substrate. At stage I, the temperature is reduced compared to the case of conventional substrate, mainly due to the heat flow manipulation of Cu film since the temperature is below the transition temperature of paraffin. At stage II, the temperature remains unchanged due to the occurrence of phase change of paraffin from sold state to liquid state and the excessive heat energy is absorbed during the phase change process. At stage III, the temperature increases monotonically to steady state, i.e., stage IV. It is shown that this thermal protecting substrate design can reduce the maximum temperature increase over 85% with optimal parameters.
The soft functional composite design in Figure 9a can only provide limited stretchability due to the existence of a stiff metal film. To further improve the stretchablity, Shi et al. [86] developed a composite with paraffin uniformly mixed with PDMS, which can be bent, twisted and stretched easily without losing the ability to absorb the excessive heat energy.
Other methods were also developed to effectively reduce the temperature through novel designs. For example, Kotagama et al. [87] adopted a thermally conductive composite and included liquid-cooled tubes to enhance the electronic cooling. Liquid metal inclusions were also embedded into the polymer substrate to achieve both flexibility and high thermal conductivity [88,89].

5. Conclusions

This paper overviews the recent advances on thermal management of FIEDs including thermal analysis, thermo-mechanical analysis and thermal designs of FIEDs with and without biological tissues to provide design guidelines for optimization. Despite these advances on thermal management of FIEDs, there are still many open challenges and opportunities for future research. For example, for FIEDs with high heat generation powers such as a flexible inorganic μ-ILED display with high resolution, active thermal design with associated thermal modeling to avoid the undesired heating is required. In addition, for FIEDs integrated with biological tissues, the thermal pain sensation based on a holistic neurodynamics model is needed for accurate modeling.

Author Contributions

Conceptualization, J.S. and Y.L.; resources, J.C.; writing—original draft preparation, J.C., and Y.L.; writing—review and editing, J.S. and S.Z.; supervision, J.S.; project administration, J.S.; funding acquisition, J.S. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos. 11772030 and 11872331), the Aeronautical Science Foundation of the PR China and Opening fund of State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology (GZ19117).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, H.S.; Brueckner, E.; Song, J.; Li, Y.; Kim, S.; Lu, C.; Sulkin, J.; Choquette, K.; Huang, Y.; Nuzzo, R.G.; et al. Unusual strategies for using indium gallium nitride grown on silicon (111) for solid-state lighting. Proc. Natl. Acad. Sci. USA 2011, 108, 10072–10077. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Webb, R.C.; Bonifas, A.P.; Behnaz, A.; Zhang, Y.H.; Yu, K.J.; Cheng, H.Y.; Shi, M.X.; Bian, Z.G.; Liu, Z.J.; Kim, Y.S.; et al. Ultrathin conformal devices for precise and continuous thermal characterization of human skin. Nat. Mater. 2013, 12, 938–944. [Google Scholar] [CrossRef] [PubMed]
  3. Tian, L.M.; Li, Y.H.; Webb, R.C.; Krishnan, S.; Bian, Z.G.; Song, J.Z.; Ning, X.; Crawford, K.; Kurniawan, J.; Bonifas, A.; et al. Flexible and stretchable 3 omega sensors for thermal characterization of human skin. Adv. Funct. Mater. 2017, 27, 1701282. [Google Scholar] [CrossRef]
  4. Lee, J.W.; Xu, R.; Lee, S.; Jang, K.I.; Yang, Y.; Banks, A.; Yu, K.J.; Kim, J.; Xu, S.; Ma, S.; et al. Soft, thin skin-mounted power management systems and their use in wireless thermography. Proc. Natl. Acad. Sci. USA 2016, 113, 6131–6136. [Google Scholar] [CrossRef] [Green Version]
  5. Wang, C.; Sim, K.; Chen, J.; Kim, H.; Rao, Z.; Li, Y.; Chen, W.; Song, J.; Verduzco, R.; Yu, C. Soft ultrathin electronics innervated adaptive fully soft robots. Adv. Mater. 2018, 30, 1706695. [Google Scholar] [CrossRef] [Green Version]
  6. Cheng, X.; Zhang, Y.H. Micro/nanoscale 3D assembly by rolling, folding, curving, and buckling approaches. Adv. Mater. 2019, 31, 1901895. [Google Scholar] [CrossRef]
  7. Song, J. Mechanics of stretchable electronics. Curr. Opin. Solid. St. M. 2015, 19, 160. [Google Scholar] [CrossRef]
  8. Cai, M.; Nie, S.; Du, Y.; Wang, C.; Song, J. Soft elastomers with programmable stiffness as strain-isolating substrates for stretchable electronics. Acs Appl. Mater. Interfaces 2019, 11, 14340–14346. [Google Scholar] [CrossRef]
  9. Yetisen, A.K.; Martinez-Hurtado, J.L.; Unal, B.; Khademhosseini, A.; Butt, H. Wearables in Medicine. Adv. Mater. 2018, 30, 1706910. [Google Scholar] [CrossRef]
  10. Sekitani, T.; Someya, T. Stretchable organic integrated circuits for large-area electronic skin surfaces. MRS Bull. 2012, 37, 236–245. [Google Scholar] [CrossRef]
  11. Qazi, R.; Kim, C.Y.; Byun, S.H.; Jeong, J.W. Microscale Inorganic LED Based Wireless Neural Systems for Chronic in vivo Optogenetics. Front. Neurosci.-SWITZ 2018, 12, 764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Linghu, C.; Zhang, S.; Wang, C.; Song, J. Transfer printing techniques for flexible and stretchable inorganic electronics. NPJ Flexible Electron. 2018, 2, 26. [Google Scholar] [CrossRef]
  13. Linghu, C.; Wang, C.; Cen, N.; Wu, J.; Lai, Z.; Song, J. Rapidly tunable and highly reversible bio-inspired dry adhesion for transfer printing in air and a vacuum. Soft Matter 2019, 15, 30–37. [Google Scholar] [CrossRef] [PubMed]
  14. Huang, Y.; Wu, H.; Xiao, L.; Duan, Y.; Zhu, H.; Bian, J.; Ye, D.; Yin, Z. Assembly and application of 3D conformal electronics on curvilinear surfaces. Mater. Horizons 2019, 6, 642–683. [Google Scholar] [CrossRef]
  15. Canals, J.; Franch, N.; Alonso, O.; Vila, A.; Dieguez, A. A Point-of-Care Device for Molecular Diagnosis Based on CMOS SPAD Detectors with Integrated Microfluidics. Sensors 2019, 19, 445. [Google Scholar] [CrossRef] [Green Version]
  16. Linghu, C.; Zhang, S.; Wang, C.; Yu, K.; Li, C.; Zeng, Y.; Zhu, H.; Jin, X.; You, Z.; Song, J. Universal SMP gripper with massive and selective capabilities for multi-scaled, arbitrarily shaped objects. Sci. Adv. 2020, 6, eaay5120. [Google Scholar] [CrossRef] [Green Version]
  17. Ma, Y.J.; Zhang, Y.C.; Cai, S.S.; Han, Z.Y.; Liu, X.; Wang, F.L.; Cao, Y.; Wang, Z.H.; Li, H.F.; Chen, Y.H.; et al. Flexible hybrid electronics for digital healthcare. Adv. Mater. 2019, 1902062. [Google Scholar] [CrossRef]
  18. Qiao, Y.C.; Wang, Y.F.; Tian, H.; Li, M.R.; Jian, J.M.; Wei, Y.H.; Tian, Y.; Wang, D.Y.; Pang, Y.; Geng, X.S.; et al. Multilayer graphene epidermal electronic skin. ACS Nano 2018, 12, 8839–8846. [Google Scholar] [CrossRef]
  19. Gao, W.; Emaminejad, S.; Nyein, H.Y.Y.; Challa, S.; Chen, K.; Peck, A.; Fahad, H.M.; Ota, H.; Shiraki, H.; Kirriya, D.; et al. Fully integrated wearable sensor arrays for multiplexed in situ perspiration analysis. Nature 2016, 529, 509–514. [Google Scholar] [CrossRef] [Green Version]
  20. Lee, H.; Choi, T.K.; Lee, Y.B.; Cho, H.R.; Ghaffari, R.; Wang, L.; Cho, H.J.; Chung, T.D.; Lu, N.S.; Hyeon, T.; et al. A graphene-based electrochemical device with thermoresponsive microneedles for diabetes monitoring and therapy. Nat. Nanotechnol. 2016, 11, 566–572. [Google Scholar] [CrossRef]
  21. Chen, Y.; Lu, S.; Zhang, S.; Li, Y.; Qu, Z.; Chen, Y.; Lu, B.; Wang, X.; Feng, X. Skin-like biosensor system via electrochemical channels for noninvasive blood glucose monitoring. Sci. Adv. 2017, 3, 1701629. [Google Scholar] [CrossRef] [Green Version]
  22. Li, H.C.; Xu, Y.; Li, X.M.; Chen, Y.; Jiang, Y.; Zhang, C.X.; Lu, B.W.; Wang, J.; Ma, Y.J.; Chen, Y.H.; et al. Epidermal inorganic optoelectronics for blood oxygen measurement. Adv. Healthc. Mater. 2017, 6, 1601013. [Google Scholar] [CrossRef] [PubMed]
  23. Yan, Z.C.; Pan, T.S.; Xue, M.M.; Chen, C.Y.; Cui, Y.; Yao, G.; Huang, L.; Liao, F.Y.; Jing, W.; Zhang, H.L.; et al. Thermal release transfer printing for stretchable conformal bioelectronics. Adv. Sci. 2017, 4, 1700251. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, S.; Wang, C.; Gao, H.; Yu, C.; Yan, Q.; Lu, Y.; Tao, Z.; Linghu, C.; Chen, Z.; Xu, K.; et al. A removable insertion shuttle for ultraflexible neural probe implantation with stable chronic brain electrophysiological recording. Adv. Mater. Interfaces 2020, 1901775. [Google Scholar] [CrossRef]
  25. Kim, D.H.; Lu, N.S.; Ma, R.; Kim, Y.S.; Kim, R.H.; Wang, S.D.; Wu, J.; Won, S.M.; Tao, H.; Islam, A.; et al. Epidermal electronics. Science 2011, 333, 838–843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lu, Y.; Huang, Y.; Yeh, J.A.; Lee, C. Controllability of non-contact cell manipulation by image dielectrophoresis (iDEP). Opt. Quantum Electron. 2005, 37, 1385–1395. [Google Scholar] [CrossRef]
  27. Wasisto, H.S.; Prades, J.D.; Gulink, J.; Waag, A. Beyond solid-state lighting: Miniaturization, hybrid integration, and applications of GaN nano- and micro-LEDs. Appl. Phys. Rev. 2019, 6, 041315. [Google Scholar] [CrossRef] [Green Version]
  28. Hsu, H.; Ohta, A.T.; Chiou, P.Y.; Jamshidi, A.; Neale, S.L.; Wu, M.C. Phototransistor-based optoelectronic tweezers for dynamic cell manipulation in cell culture media. Lab Chip 2010, 10, 165–172. [Google Scholar] [CrossRef]
  29. Ohta, A.T.; Chiou, P.Y.; Phan, H.L.; Sherwood, S.W.; Yang, J.M.; Lau, A.N.K.; Hsu, H.Y.; Jamshidi, A.; Wu, M.C. Optically controlled cell discrimination and trapping using optoelectronic tweezers. IEEE J. Sel. Top. Quantum Electron. 2007, 13, 235–243. [Google Scholar] [CrossRef] [Green Version]
  30. Ohta, A.T.; Chiou, P.Y.; Han, T.H.; Liao, J.C.; Bhardwaj, U.; McCabe, E.R.B.; Yu, F.Q.; Sun, R.; Wu, M.C. Dynamic cell and microparticle control via optoelectronic tweezers. J. Microelectromech. Syst. 2007, 16, 491–499. [Google Scholar] [CrossRef]
  31. Neale, S.L.; Ohta, A.T.; Hsu, H.Y.; Valley, J.K.; Jamshidi, A.; Wu, M.C. Trap profiles of projector based optoelectronic tweezers (OET) with HeLa cells. Opt. Express 2009, 17, 5231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Scholz, G.; Mariana, S.; Dharmawan, A.B.; Syamsu, I.; Hormann, P.; Reuse, C.; Hartmann, J.; Hiller, K.; Prades, J.D.; Wasisto, H.S.; et al. Continuous Live-Cell Culture Imaging and Single-Cell Tracking by Computational Lensfree LED Microscopy. Sensors 2019, 19, 1234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Hwang, H.; Park, J.K. Measurement of molecular diffusion based on optoelectrofluidic fluorescence microscopy. Anal. Chem. 2009, 81, 9163–9167. [Google Scholar] [CrossRef] [PubMed]
  34. Lin, Y.H.; Chang, C.M.; Lee, G.B. Manipulation of single DNA molecules by using optically projected images. Opt. Express 2009, 17, 15318. [Google Scholar] [CrossRef] [PubMed]
  35. Hoeb, M.; Raedler, J.O.; Klein, S.; Stutzmann, M.; Brandt, M.S. Light-induced dielectrophoretic manipulation of DNA. Biophys. J. 2007, 93, 1032–1038. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Chiou, P.Y.; Ohta, A.T.; Wu, M.C. Massively parallel manipulation of single cells and microparticles using optical images. Nature 2005, 436, 370–372. [Google Scholar] [CrossRef]
  37. Huang, K.W.; Su, T.W.; Ozcan, A.; Chiou, P.Y. Optoelectronic tweezers integrated with lensfree holographic microscopy for wide-field interactive cell and particle manipulation on a chip. Lab Chip 2013, 13, 2278. [Google Scholar] [CrossRef]
  38. McCall, J.G.; Kim, T.I.; Shin, G.; Huang, X.; Jung, Y.H.; Al-Hasani, R.; Omenetto, F.G.; Bruchas, M.R.; Rogers, J.A. Fabrication and application of flexible, multimodal lightemitting devices for wireless optogenetics. Nat. Protoc. 2013, 8, 2413–2428. [Google Scholar] [CrossRef] [Green Version]
  39. Jeong, J.W.; McCall, J.G.; Shin, G.; Zhang, Y.H.; Al-Hasani, R.; Kim, M.; Li, S.; Sim, J.Y.; Jang, K.I.; Shi, Y.; et al. Wireless optofluidic systems for programmable in vivo pharmacology and optogenetics. Cell 2015, 162, 662–674. [Google Scholar] [CrossRef] [Green Version]
  40. Pisanello, F.; Sileo, L.; Vittorio, M.D. Micro- and nanotechnologies for optical neural interfaces. Front. Neurosci. 2016, 10, 468. [Google Scholar] [CrossRef] [Green Version]
  41. Yang, F.; Gao, C.; Wang, P.; Zhang, G.J.; Chen, Z. Fish-on-a-chip: Microfluidics for zebrafish research. Lab Chip 2016, 16, 1106–1125. [Google Scholar] [CrossRef] [PubMed]
  42. Al-Atabany, W.; McGovern, B.; Mehran, K.; Berlinguer-Palmini, R.; Degenaar, P. A processing platform for optoelectronic/optogenetic retinal prosthesis. IEEE Trans. Biomed. Eng. 2013, 60, 781–791. [Google Scholar] [CrossRef] [PubMed]
  43. Barrett, J.M.; Berlinguer-Palmini, R.; Degenaar, P. Optogenetic approaches to retinal prosthesis. Vis. Neurosci. 2014, 31, 345–354. [Google Scholar] [CrossRef] [PubMed]
  44. Degenaar, P.; Grossman, N.; Memon, M.A.; Burrone, J.; Dawson, M.; Drakakis, E.; Neil, M.; Nikolic, K. Optobionic vision—A new genetically enhanced light on retinal prosthesis. J. Neural Eng. 2009, 6, 35007. [Google Scholar] [CrossRef]
  45. Moser, T. Optogenetic stimulation of the auditory pathway for research and future prosthetics. Curr. Opin. Neurobiol. 2015, 34, 29–36. [Google Scholar] [CrossRef]
  46. Gossler, C.; Bierbrauer, C.; Moser, R.; Kunzer, M.; Holc, K.; Pletschen, W.; Kohler, K.; Wagner, J.; Schwaerzle, M.; Ruther, P.; et al. GaN-based micro-LED arrays on flexible substrates for optical cochlear implants. J. Phys. D: Appl. Phys. 2014, 47, 205401. [Google Scholar] [CrossRef]
  47. Hernandez, V.H.; Gehrt, A.; Reuter, K.; Jing, Z.; Jeschke, M.; Schulz, A.M.; Hoch, G.; Bartels, M.; Vogt, G.; Garnham, C.W.; et al. Optogenetic stimulation of the auditory pathway. J. Clin. Invest. 2014, 124, 1114–1129. [Google Scholar] [CrossRef] [Green Version]
  48. Jeschke, M.; Moser, T. Considering optogenetic stimulation for cochlear implants. Hear. Res. 2015, 322, 224–234. [Google Scholar] [CrossRef]
  49. Klein, E.; Gossler, C.; Paul, O.; Ruther, P. High-density lLED-based optical cochlear implant with improved thermomechanical behavior. Front. Neurosci. 2018, 12, 659. [Google Scholar] [CrossRef]
  50. Kim, T.I.; Jung, Y.H.; Song, J.Z.; Kim, D.; Li, Y.H.; Kim, H.S.; Song, I.S.; Wierer, J.J.; Pao, H.A.; Huang, Y.G.; et al. High-Efficiency, microscale GaN light-emitting diodes and their thermal properties on unusual substrates. Small 2012, 8, 1643–1649. [Google Scholar] [CrossRef]
  51. Kim, C.; Kim, C.H. Universal Testing Apparatus Implementing Various Repetitive Mechanical Deformations to Evaluate the Reliability of Flexible Electronic Devices. Micromachines 2018, 9, 492. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Baek, K.; Jang, K.; Lee, Y.J.; Ryu, K.; Choi, W.; Kim, D.; Yi, J. Improvement of mechanical reliability by patterned silver/Indium-Tin-Oxide structure for flexible electronic devices. Thin Solid Films 2013, 531, 349–353. [Google Scholar] [CrossRef]
  53. Song, J.; Feng, X.; Huang, Y. Mechanics and thermal management of stretchable inorganic electronics. Nat. Sci. Rev. 2016, 3, 128–143. [Google Scholar] [CrossRef]
  54. Cui, Y.; Li, Y.H.; Xing, Y.F.; Ji, Q.G.; Song, J.Z. Thermal design of rectangular microscale inorganic light-emitting diodes. Appl. Therm. Eng. 2017, 122, 653–660. [Google Scholar] [CrossRef]
  55. Lu, C.F.; Li, Y.H.; Song, J.Z.; Kim, H.S.; Brueckner, E.; Fang, B.; Hwang, K.C.; Huang, Y.G.; Nuzzo, R.G.; Rogers, J.A. A thermal analysis of the operation of microscale, inorganic light-emitting diodes. P. Roy. Soc. A-Math. Phy. 2012, 468, 3215–3223. [Google Scholar] [CrossRef]
  56. Yin, Y.F.; Li, M.; Yuan, W.; Chen, X.L.; Li, Y.H. A widely adaptable analytical method for thermal analysis of flexible electronics with complex heat source structures. P. Roy. Soc. A-Math. Phy. 2019, 475, 20190402. [Google Scholar] [CrossRef] [Green Version]
  57. Cui, Y.; Bian, Z.G.; Li, Y.H.; Xing, Y.F.; Song, J.Z. 3D thermal analysis of rectangular microscale inorganic light-emitting diodes in a pulsed operation. J. Phys. D: Appl. Phys. 2016, 49, 405101. [Google Scholar] [CrossRef]
  58. Chen, X.J.; Su, Y.H.; Reay, D.; Riffat, S. Recent research developments in polymer heat exchangers—A review. Renew. Sust. Energ. Rev. 2016, 60, 1367–1386. [Google Scholar] [CrossRef]
  59. Chen, H.Y.; Ginzburg, V.V.; Yang, J.; Yang, Y.F.; Liu, W.; Huang, Y.; Du, L.B.; Chen, B. Thermal conductivity of polymer-based composites: Fundamentals and applications. Prog. Polym. Sci. 2016, 59, 41–85. [Google Scholar] [CrossRef]
  60. Cevallos, J.G.; Bergles, A.E.; Bar-Cohen, A.; Rodgers, P.; Gupta, S.K. Polymer Heat Exchangers-History, Opportunities, and Challenges. Heat Transfer Eng. 2012, 33, 1075–1093. [Google Scholar] [CrossRef]
  61. Mamunya, Y.P.; Davydenko, V.V.; Pissis, P.; Lebedev, E. Electrical and thermal conductivity of polymers filled with metal powders. Eur. Polym. J. 2002, 38, 1887–1897. [Google Scholar] [CrossRef]
  62. Gu, J.W.; Zhang, Q.Y.; Dang, J.; Zhang, J.P.; Yang, Z.Y. Thermal Conductivity and Mechanical Properties of Aluminum Nitride Filled Linear Low-Density Polyethylene Composites. Polym. Eng. Sci. 2009, 49, 1030–1034. [Google Scholar] [CrossRef]
  63. Leung, S.N.; Khan, M.O.; Chan, E.; Naguib, H.; Dawson, F.; Adinkrah, V.; Lakatos-Hayward, L. Analytical modeling and characterization of heat transfer in thermally conductive polymer composites filled with spherical particulates. Compos. Part B-Eng. 2013, 45, 43–49. [Google Scholar] [CrossRef]
  64. Chen, Y.M.; Ting, J.M. Ultra high thermal conductivity polymer composites. Carbon 2002, 40, 359–362. [Google Scholar] [CrossRef]
  65. Chen, L.; Sun, Y.Y.; Lin, J.; Du, X.Z.; Wei, G.S.; He, S.J.; Nazarenko, S. Modeling and analysis of synergistic effect in thermal conductivity enhancement of polymer composites with hybrid filler. Int. J. Heat Mass Tran. 2015, 81, 457–464. [Google Scholar] [CrossRef]
  66. Kim, P.; Shi, L.; Majumdar, A.; McEuen, P.L. Mesoscopic thermal transport and energy dissipation in carbon nanotubes. Phys. B: Condens. Matter 2002, 323, 67–70. [Google Scholar] [CrossRef]
  67. Xu, Y.S.; Ray, G.; Abdel-Magid, B. Thermal behavior of single-walled carbon nanotube polymer-matrix composites. Compos. Part A: Appl. Sci. Manuf. 2006, 37, 114–121. [Google Scholar] [CrossRef]
  68. Wang, Q.; Han, X.H.; Sommers, A.; Park, Y.; T’Joen, C.; Jacobi, A. A review on application of carbonaceous materials and carbon matrix composites for heat exchangers and heat sinks. Int. J. Refrig. 2012, 35, 7–26. [Google Scholar] [CrossRef]
  69. Zhang, J.P.; Li, Y.H.; Xing, Y.F. Theoretical and experimental investigations of transient thermo-mechanical analysis on flexible electronic devices. Int. J. Mech. Sci. 2019, 160, 192–199. [Google Scholar] [CrossRef]
  70. Cui, Y.; Li, Y.H.; Xing, Y.F.; Yang, T.Z.; Song, J.Z. One-dimensional thermal analysis of the flexible electronic devices integrated with human skin. Micromachines 2016, 7, 210. [Google Scholar] [CrossRef] [Green Version]
  71. Cui, Y.; Li, Y.H.; Xing, Y.F.; Yang, T.Z.; Song, J.Z. Three-dimensional thermal analysis of rectangular micro-scale inorganic light-emitting diodes integrated with human skin. Int. J. Therm. Sci. 2018, 127, 321–328. [Google Scholar] [CrossRef]
  72. Yin, Y.F.; Cui, Y.; Li, Y.H.; Xing, Y.F. Thermal management of flexible wearable electronic devices integrated with human skin considering clothing effect. Appl. Therm. Eng. 2018, 144, 504–511. [Google Scholar] [CrossRef]
  73. Nie, S.; Zhang, C.; Song, J. Thermal management of epidermal electronic devices/skin system considering insensible sweating. Sci. Rep.-UK 2018, 8, 14121. [Google Scholar] [CrossRef] [PubMed]
  74. Li, Y.; Zhang, J.; Xing, Y.; Song, J. Thermal analysis of epidermal electronic devices integrated with human skin considering the effects of interfacial thermal resistance. AIP Adv. 2018, 8, 055102. [Google Scholar] [CrossRef] [Green Version]
  75. Yin, Y.F.; Li, Y.H.; Li, M. Thermal analysis of the flexible electronics affixed on large curvature myocardium surface. Int. J. Heat Mass Tran. 2020, 147, 118983. [Google Scholar] [CrossRef]
  76. Li, Y.H.; Shi, X.T.; Song, J.H.; Lu, C.F.; Kim, T.I.; McCall, J.G.; Bruchas, M.R.; Rogers, J.A.; Huang, Y.G. Thermal analysis of injectable, cellular-scale optoelectronics with pulsed power. P. Roy. Soc. A-Math. Phy. 2013, 469, 2156. [Google Scholar] [CrossRef] [PubMed]
  77. Kim, T.I.; McCall, J.G.; Jung, Y.H.; Huang, X.; Siuda, E.R.; Li, Y.H.; Song, J.Z.; Song, Y.M.; Pao, H.A.; Kim, R.H.; et al. Injectable, cellular-scale optoelectronics with applications for wireless optogenetics. Science 2013, 340, 6129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Li, Y.J.; Ma, Y.J.; Wei, C.; Luan, H.W.; Xu, S.; Han, M.D.; Zhao, H.B.; Liang, C.M.; Yang, Q.S.; Yang, Y.Y.; et al. Thin, Millimeter scale fingernail sensors for thermal characterization of nail bed tissue. Adv. Funct. Mater. 2018, 28, 1801380. [Google Scholar] [CrossRef]
  79. Li, Y.H.; Zhang, J.P.; Xing, Y.F.; Song, J.Z. Thermomechanical analysis of epidermal electronic devices integrated with human skin. J. Appl. Mech.-T ASME 2017, 84, 111004. [Google Scholar] [CrossRef]
  80. Zhang, J.; Li, Y.; Xing, Y.; Song, J. Three-dimensional thermomechanical analysis of epidermal electronic devices on human skin. Int. J. Solids Struct. 2019, 167, 48–57. [Google Scholar] [CrossRef]
  81. Zhang, J.P.; Li, Y.H.; Xing, Y.F.; Song, J.Z. Transient thermomechanical analysis of epidermal electronic devices on human skin. Mech. Mater. 2019, 137, 103097. [Google Scholar] [CrossRef]
  82. Jung, H.H.; Song, J.; Nie, S.; Jung, H.N.; Kim, M.S.; Jeong, J.W.; Song, Y.M.; Song, J.; Jang, K.I. Thin metallic heat sink for interfacial thermal management in biointegrated optoelectronic devices. Adv. Mater. Technol. 2018, 3, 1800159. [Google Scholar] [CrossRef]
  83. Li, Y.; Chen, J.; Xing, Y.; Song, J. Analytical investigations on the thermal properties of microscale inorganic light-emitting diodes on an orthotropic substrate. AIP Adv. 2017, 7, 035208. [Google Scholar] [CrossRef] [Green Version]
  84. Li, Y.H.; Chen, J.; Xing, Y.F.; Song, J.Z. Thermal management of micro-scale inorganic light-emittng diodes on an orthotropic substrate for biointegrated applications. Sci. Rep.-UK 2017, 7, 6638. [Google Scholar] [CrossRef] [Green Version]
  85. Shi, Y.L.; Wang, C.J.; Yin, Y.F.; Li, Y.H.; Xing, Y.F.; Song, J.Z. Functional soft composites as thermal protecting substrates for wearable electronics. Adv. Funct. Mater. 2019, 29, 1905470. [Google Scholar] [CrossRef]
  86. Shi, Y.L.; Hu, M.; Xing, Y.F.; Li, Y.H. Temperature-dependent thermal and mechanical properties of flexible functional PDMS/paraffin composites. Mater. Design 2020, 185, 108219. [Google Scholar] [CrossRef]
  87. Kotagama, P.; Phadnis, A.; Manning, K.C.; Rykaczewski, K. Rational Design of Soft, Thermally Conductive Composite Liquid-Cooled Tubes for Enhanced Personal, Robotics, and Wearable Electronics Cooling. Adv. Mater. Technol. 2019, 4, 1800690. [Google Scholar] [CrossRef]
  88. Jeong, S.H.; Chen, S.; Huo, J.X.; Gamstedt, E.K.; Liu, J.H.; Zhang, S.L.; Zhang, Z.B.; Hjort, K.; Wu, Z.G. Mechanically Stretchable and Electrically Insulating Thermal Elastomer Composite by Liquid Alloy Droplet Embedment. Sci. Rep. 2015, 5, 18257. [Google Scholar] [CrossRef]
  89. Bartlett, M.D.; Kazem, N.; Powell-Palm, M.J.; Huang, X.N.; Sun, W.H.; Malen, J.A.; Majidi, C. High thermal conductivity in soft elastomers with elongated liquid metal inclusions. Proc. Natl. Acad. Sci. USA 2017, 114, 2143–2148. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Examples of flexible inorganic electronic devices. (a) A flexible inorganic microscale inorganic light-emitting diode (μ-ILED) display [1]. (b) An epidermal temperature sensor array in a twisting motion [2]. (c) A flexible and stretchable 3ω sensor [3]. (d) A soft, thin skin-mounted power management system [4]. (e) A soft ultrathin electronics enabled for fully soft robots [5].
Figure 1. Examples of flexible inorganic electronic devices. (a) A flexible inorganic microscale inorganic light-emitting diode (μ-ILED) display [1]. (b) An epidermal temperature sensor array in a twisting motion [2]. (c) A flexible and stretchable 3ω sensor [3]. (d) A soft, thin skin-mounted power management system [4]. (e) A soft ultrathin electronics enabled for fully soft robots [5].
Micromachines 11 00390 g001
Figure 2. Schematic illustration of (a) the μ-ILED layout and (b) a quarter geometry of the analytical modeled μ-ILED system [54]. (c) The shape influence on the temperature increase of μ-ILEDs with 100 μm × 2b × 5 μm under the input power of 12 mW [54]. (d) The distribution of the surface temperature increase of a single rectangular μ-ILED under the heat generation power 12 mW [54]. (e) Schematic diagram and (f) the distribution of surface temperature increase of a μ-ILED array consisting of two μ-ILEDs [54].
Figure 2. Schematic illustration of (a) the μ-ILED layout and (b) a quarter geometry of the analytical modeled μ-ILED system [54]. (c) The shape influence on the temperature increase of μ-ILEDs with 100 μm × 2b × 5 μm under the input power of 12 mW [54]. (d) The distribution of the surface temperature increase of a single rectangular μ-ILED under the heat generation power 12 mW [54]. (e) Schematic diagram and (f) the distribution of surface temperature increase of a μ-ILED array consisting of two μ-ILEDs [54].
Micromachines 11 00390 g002
Figure 3. Schematic illustration of (a) the μ-ILED layout and (b) a quarter geometry of the analytic modeled μ-ILED system [57]. (c) The pulsed power Q(t) with Q0 as the peak power [57]. (d) The μ-ILED temperature increase versus time [57] after saturation [57]. (e) The maximum and the minimum temperature increases versus duty cycle after saturation [57]. (f) The maximum normalized temperature increase of μ-ILED versus duty cycle with various values of σ [57].
Figure 3. Schematic illustration of (a) the μ-ILED layout and (b) a quarter geometry of the analytic modeled μ-ILED system [57]. (c) The pulsed power Q(t) with Q0 as the peak power [57]. (d) The μ-ILED temperature increase versus time [57] after saturation [57]. (e) The maximum and the minimum temperature increases versus duty cycle after saturation [57]. (f) The maximum normalized temperature increase of μ-ILED versus duty cycle with various values of σ [57].
Micromachines 11 00390 g003
Figure 4. (a) Schematic diagram of a flexible electronic device consisting of a functional component on a substrate encapsulated by an encapsulation layer [69]. (b) The temperature increase versus time at the center point on the top surface. [69]. The comparisons of (c) x-stain, (d) y-strain and (e) shear strain from experiments, finite element analysis (FEA) and analytical model [69].
Figure 4. (a) Schematic diagram of a flexible electronic device consisting of a functional component on a substrate encapsulated by an encapsulation layer [69]. (b) The temperature increase versus time at the center point on the top surface. [69]. The comparisons of (c) x-stain, (d) y-strain and (e) shear strain from experiments, finite element analysis (FEA) and analytical model [69].
Micromachines 11 00390 g004
Figure 5. Schematic illustration of (a) the skin tissue structure and (b) the geometry of the analytical modeled system [70]. (c) The comparison of temperature increase along the thickness direction between the analytical prediction and FEA [70]. (d) The three-dimensional analytical model showing a quarter geometry of the μ-ILED device integrated with skin tissue [71]. (e) The dependence of the thermal comfort map on the thickness and thermal conductivity of the substrate [71].
Figure 5. Schematic illustration of (a) the skin tissue structure and (b) the geometry of the analytical modeled system [70]. (c) The comparison of temperature increase along the thickness direction between the analytical prediction and FEA [70]. (d) The three-dimensional analytical model showing a quarter geometry of the μ-ILED device integrated with skin tissue [71]. (e) The dependence of the thermal comfort map on the thickness and thermal conductivity of the substrate [71].
Micromachines 11 00390 g005
Figure 6. (a) Schematic diagram of the cross-sectional structure for the flexible inorganic electronic device (FIED)/skin system [79]. (b) The distribution of the maximum principal stress at the device/skin interface along the radial direction [79]. (c) The maximum temperature increase and (d) the maximum principal stress at the device/skin interface vary with the size of heating component and the substrate thickness [79].
Figure 6. (a) Schematic diagram of the cross-sectional structure for the flexible inorganic electronic device (FIED)/skin system [79]. (b) The distribution of the maximum principal stress at the device/skin interface along the radial direction [79]. (c) The maximum temperature increase and (d) the maximum principal stress at the device/skin interface vary with the size of heating component and the substrate thickness [79].
Micromachines 11 00390 g006
Figure 7. Schematic diagrams of (a) the multilayer structure of the device/skin system with one rectangular heating component and (b) the analytical modeled system with a quarter of the geometry [81]. (c) Periodic pulsed power density versus time [81]. (d) Maximum principal strain. (e) Maximum rate of principal strain as functions of time in the skin [81]. (f) Influences of thicknesses of encapsulation and substrate on the peak values of the maximum rate of the principal strain [81].
Figure 7. Schematic diagrams of (a) the multilayer structure of the device/skin system with one rectangular heating component and (b) the analytical modeled system with a quarter of the geometry [81]. (c) Periodic pulsed power density versus time [81]. (d) Maximum principal strain. (e) Maximum rate of principal strain as functions of time in the skin [81]. (f) Influences of thicknesses of encapsulation and substrate on the peak values of the maximum rate of the principal strain [81].
Micromachines 11 00390 g007
Figure 8. (a) Scheme diagram of the μ-ILED device on an orthotropic substrate consisting of two layered materials [83]. (b) The temperature increase of the μ-ILED versus the thermal conductivity ratio [83]. Distributions of temperature increase for the μ-ILED device on the (c) isotropic and (d) orthotropic substrate [83].
Figure 8. (a) Scheme diagram of the μ-ILED device on an orthotropic substrate consisting of two layered materials [83]. (b) The temperature increase of the μ-ILED versus the thermal conductivity ratio [83]. Distributions of temperature increase for the μ-ILED device on the (c) isotropic and (d) orthotropic substrate [83].
Micromachines 11 00390 g008
Figure 9. (a) Cross-sectional illustration of a skin-integrated wearable device with a functional soft composite substrate design [85]. (b) Schematic illustration of a wearable device consisting of a heat source on the functional soft composite substrate with the embedded paraffin as the phase change material and the Cu film on its top in PDMS [85]. (c) Optical image of a wearable device under bending deformation [85]. Temperature distributions on the bottom of the PDMS substrate and thermal protecting substrate from (d) experiments and (e) FEA after heated 40 s [85]. The maximum temperature on the bottom of the thermal protecting substrate versus time with (f) a heating duration of 40 s and (g) a long heating duration [85].
Figure 9. (a) Cross-sectional illustration of a skin-integrated wearable device with a functional soft composite substrate design [85]. (b) Schematic illustration of a wearable device consisting of a heat source on the functional soft composite substrate with the embedded paraffin as the phase change material and the Cu film on its top in PDMS [85]. (c) Optical image of a wearable device under bending deformation [85]. Temperature distributions on the bottom of the PDMS substrate and thermal protecting substrate from (d) experiments and (e) FEA after heated 40 s [85]. The maximum temperature on the bottom of the thermal protecting substrate versus time with (f) a heating duration of 40 s and (g) a long heating duration [85].
Micromachines 11 00390 g009

Share and Cite

MDPI and ACS Style

Li, Y.; Chen, J.; Zhao, S.; Song, J. Recent Advances on Thermal Management of Flexible Inorganic Electronics. Micromachines 2020, 11, 390. https://doi.org/10.3390/mi11040390

AMA Style

Li Y, Chen J, Zhao S, Song J. Recent Advances on Thermal Management of Flexible Inorganic Electronics. Micromachines. 2020; 11(4):390. https://doi.org/10.3390/mi11040390

Chicago/Turabian Style

Li, Yuhang, Jiayun Chen, Shuang Zhao, and Jizhou Song. 2020. "Recent Advances on Thermal Management of Flexible Inorganic Electronics" Micromachines 11, no. 4: 390. https://doi.org/10.3390/mi11040390

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop