Next Article in Journal
Gaussian Multiple Access Channels with One-Bit Quantizer at the Receiver ,
Previous Article in Journal
Multi-Objective Evolutionary Rule-Based Classification with Categorical Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constructal Optimization for Cooling a Non-Uniform Heat Generating Radial-Pattern Disc by Conduction

1
Institute of Thermal Science and Power Engineering, Naval University of Engineering, Wuhan 430033, China
2
Military Key Laboratory for Naval Ship Power Engineering, Naval University of Engineering, Wuhan 430033, China
3
College of Power Engineering, Naval University of Engineering, Wuhan 430033, China
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(9), 685; https://doi.org/10.3390/e20090685
Submission received: 4 August 2018 / Revised: 29 August 2018 / Accepted: 5 September 2018 / Published: 7 September 2018
(This article belongs to the Section Thermodynamics)

Abstract

:
A heat conduction model in a radial-pattern disc by considering non-uniform heat generation (NUHG) is established in this paper. A series of high conductivity channels (HCCs) are attached on the rim of the disc and extended to its center. Constructal optimizations of the discs with constant and variable cross-sectional HCCs are carried out, respectively, and their maximum temperature differences (MTDs) are minimized based on analytical method and finite element method. Besides, the influences of the NUHG coefficient, HCC number and width coefficient on the optimal results are studied. The results indicate that the deviation of the optimal constructs obtained from the analytical method and finite element method are comparatively slight. When the NUHG coefficient is equal to 10, the minimum MTD of the disc with 25 constant cross-sectional HCCs is specifically reduced by 48.8% compared to that with 10 HCCs. As a result, the heat conduction performance (HCP) of the disc can be efficiently improved by properly increasing the number of HCCs. The minimum MTD of the disc with variable cross-sectional HCC is decreased by 15.0% when the width coefficient is changed from 1 to 4. Therefore, the geometry of variable cross-sectional HCC can be applied in the constructal design of the disc to a better heat transfer performance. The constructal results obtained by investigating the non-uniform heat generating case in this paper can contribute to the design of practical electronic device to a better heat transfer performance.

1. Introduction

Nowadays, with the fast-developing electronic science and technology, the electronic components are faced with new trends of miniaturization, low power consumption, compound design and modularization, which are unprecedented opportunities and tremendous challenges for the electronic industry. As to miniature electronic components, the conventional heat transfer approaches such as convection or radiation are inevitably impossible to realize due to the space constraints. Accordingly, heat conduction has naturally been the prospective strategy for the effective cooling of electronic devices, and one way to enhance heat conduction is arranging high conductivity material. How to properly arrange the high conductivity material to reduce the hot spot temperature of the heat generating area is a hot topic in the context of heat conduction optimization [1,2,3].
Constructal theory [4,5,6] has emerged as the evolutionary design philosophy for finite size flow system evolution through time and it states that this evolution is based on a simple law: “for a finite size flow system to survive in time, it should evolve in such a way that it provides easier access to the flow that goes through it”. Invoking this law in the engineering field, one should bear in mind that this law is about the direction of evolution in time, and the design phenomenon is dynamic. Furthermore, the connotation of the time arrow design evolution has become the focus of concern by many researches and scholars [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]. Bejan [27] firstly conducted constructal optimization of a rectangular heat generating volume with the optimization objective of maximum temperature difference (MTD) and obtained the optimal distribution of constant high conductivity channels (HCCs). Neagu and Bejan [28,29] introduced a new geometric feature of variable cross-sectional high conductivity channels into the constructal method with the aim of further minimizing the thermal resistance. Thereafter, many scholars performed constructal optimizations of various heat generating bodies to improve their heat conduction performance (HCP), such as rectangular bodies [30,31,32,33], square bodies [34,35,36,37,38,39,40,41], triangular bodies [42,43,44,45], disc-shaped bodies with radial-pattern [46,47] and tree-shaped HCCs [48,49,50,51,52,53,54], respectively.
In the constructal designs of disc-shaped bodies, da Silva et al. [46] proposed a new heat conduction model in a radial-pattern disc in which the rectangular HCCs were mounted on the perimeter and protruded into the center. By taking the MTD as optimization objective, the constructal optimizations were carried out by analytical and numerical methods, respectively. The results showed that the number of HCCs and conductivity ratio had great effects on the optimal constructs of the disc. Sharifi et al. [47] carried out an optimization for the disc inserted with incomplete variable cross-sectional HCCs. The relevant results reported that the incomplete variable cross-sectional HCCs structure could further reduce the thermal resistance of the disc. In addition to the researches of radial-pattern disc, tree-shaped disc has been widely studied in the meantime. Rocha et al. [48] established a uniform heat generating model in a disc, and obtained the optimal arrangements of tree-shaped HCCs with the objective of MTD. Rocha et al. [49] further verified the results obtained in Ref. [46] by numerical calculations, and optimized the disc which included tree-shaped HCCs with loops. Xiao et al. [50] re-optimized the tree-shaped disc by releasing the constraint that the optimized last order constructs constituted the new order construct, and significantly reduced the MTD by 49.3%. Bahadormanesh and Salimpour [51] studied the disc-shaped body inserted with high-emissivity inserts, and obtained the best configurations of the disc by analytical and numerical methods, respectively. Sharfi et al. [52] investigated a solid disc using incomplete inserts of high thermal conductivity with the objective of thermal resistance minimization. Furthermore, the constructal studies of the heat conduction discs were extended from conventional scale to micro and nano scales. Chen et al. [53] re-optimized the disc-shaped body under the circumstances of micro and nano scales, and the results indicated that the constructs with size effects could improve the HCP with authority. Daneshi et al. [54] numerically calculated the disc-shaped body with conductive tree-shaped HCCs at micro- and nanoscale, and the results manifested that increasing the number of HCCs didn’t appear to necessarily improve the HCP of the disc.
However, as a matter of fact, the heat generation of practical electronic devices are usually non-uniformly distributed due to the working conditions and operating environment. In this context, considering the non-uniform heat generating phenomenon makes the research outcomes closer to real life problems in the field of electronics cooling, which is of great pragmatic significance. As mentioned above, most of the constructal researches in recent years were apparently carried out under the precondition of uniform heat generation, which was incompatible with the actual situation. Cetkin and Oliani [55] investigated rectangular constructs with non-uniform heat generation (NUHG), and optimized the shapes of HCCs on account of the minimum MTD. Besides, they analyzed the influences of HCC locations and shapes on the hot spot temperature in the premise of linear and local heat generation, respectively. Feng et al. [56,57] constructed the NUHG rectangular bodies inserted with constant and variable cross-sectional HCCs, respectively, and numerically obtained the optimal constructs with the criteria of minimum MTD [56] and minimum entransy dissipation rate [57], respectively. You et al. [58] performed constructal optimization of a NUHG triangular body, and theoretically and numerically optimized the first order constructs with constant and variable cross-sectional HCCs, respectively. The corresponding results showed that the minimum MTD of the constructs with variable cross-sectional HCCs was decreased by 12.6% compared with the one with constant cross-sectional HCCs, which thereby meant that the HCP could be further improved by adopting the variable cross-sectional HCCs structure in this regard.
Among the heat conduction constructal investigations, the disc model with NUHG is rarely considered. On the basis of [28,46,55,56,57,58], a conduction model in the radial-pattern disc considering NUHG phenomenon will be established in this paper. Constructal optimizations of the radial-pattern disc with constant cross-sectional HCCs will be carried out analytically and numerically, and the minimum MTDs derived by the two methods will be compared. In order to quest for a better heat transfer performance, discarding the constraint that the width of HCC is constant, variable cross-sectional HCC architecture [28] will be further utilized in the radial-pattern disc model in the meantime. Thereafter, the influences of heat generating distribution and HCC shapes on the optimal results will be analyzed.

2. Constructal Optimization of a Radial-Pattern Disc with Analytical Solution

A generic heat conduction model of a radial-pattern disc is schematically shown in Figure 1. As shown in this figure, the radius of the disc is R , and the heat generating rate (HGR) in the disc (thermal conductivity k 0 ) is q = q 0 f ( r ) , where q 0 refers to the heat generating constant and f ( r ) refers to the HGR function. The heat generated in the radial-pattern disc converges into N HCCs (length L and width D ) which are uniformly distributed in the circumferential direction. By this means, the heat flows out of the disc-shaped body from the ends of HCCs which are situated on the rim of the disc.
Given that the total amount of high conductivity material (thermal conductivity k p ) is fixed, the HCCs could extend radially from the perimeter to the center of the disc. In this manner, the area fraction ϕ is defined as the ratio of k p area to the disc area:
ϕ = D L N / ( π R 2 )
According to the distribution of the HGR and HCCs, the radial-pattern disc can be divided into N ( = 2 π / α ) identical sectorial elements (i.e., the number of HCCs in the disc), where α is the apex angle of each sectorial element. When the apex angle α ( = 2 π / N ) is specified, the HCP of the radial-pattern disc can vary with the length ratio λ ( = L / D ) of HCCs. Supposing that the apex angle α is sufficiently small, each sectorial element can be approximately viewed as an isosceles triangle, whose base and height are α R and R , respectively. In this way, as shown in Figure 2, the triangular element (TE) consists of two parts, i.e., a trapezoid area which contains an HCC and a wedge-shaped area near the apex of the TE.
The total heat current generated in the sectorial area per unit of thickness is:
q = q 0 α R 2 / 2
Clearly, the MTD in the TE appears to be T max T min , where T max denotes the hot spot temperature and T min denotes the HCC temperature on the rim. With the exception of the end of HCC ( T min ), the arc boundary is the same adiabatic condition as the radial boundaries, which are sketched with dashed lines. In the case that the TE’s apex angle α is sufficiently small, the hot spot may occur at the center of the disc (temperature T c ), or may probably occur on the TE’s base far from the HCC (temperature T r ). In this manner, the hot spot temperature can be inspected as:
T max = max ( T c , T r )
The dimensionless temperature difference in the TE is defined as:
Δ T ˜ = ( T T min ) / ( q 0 R 2 / k 0 )
Likewise, the structure parameters related to the TE are non-dimensionalized as:
( L ˜ , D ˜ , r ˜ ) = ( L , D , r ) / R
In line with Ref. [57], assuming that the average HGR of the whole disc is a constant, the HGR decreases gradually along the radius, so that the HGR function f ( r ) is written as f ( r ) = 1 + 0.1 p 0.15 p r ˜ , in which p is termed as the nonuniform heat generating (NUHG) coefficient.
The heat current q t i p generated in the wedge-shaped area (radius R L and apex angle α ) is calculated as:
q t i p = q 0 ( 10 + p L ˜ ) / 20 α ( R L ) 2
Moreover, the heat flux q r collected at the longitudinal direction of y axis is given as:
q r = q 0 ( 1 + 0.1 p 0.15 p r ˜ ) ( α r D )
According to [46,57], in the context of ϕ < < 1 , k 0 < < k p and sufficiently small α , the heat conduction direction in the k p material can be approximately viewed as parallel to the r axis, while that in the k 0 material can be viewed as parallel to the y axis. Based on the above hypotheses, the heat conduction differential equation in the k p material is:
d d r ( k p D d T d r ) = q r
The corresponding boundary conditions are invoked as:
k p D d T d r = q t i p , r = R L
T = T min ,   r = R
By solving Equations (8)–(10), the dimensionless temperature difference along the HCC is:
T 0 T min q 0 R 2 / k 0 = λ k ˜ { α [ p 80 L ˜ 3 + 1 2 + L ˜ 2 ( 1 6 p 30 ) + L ˜ ( 1 2 + 1 40 p ) ] + D ˜ L ˜ ( 1 2 p 20 L ˜ + p 40 ) }
where k ˜ denotes the ratio of k p to k 0 , and T 0 denotes the temperature of the HCC terminal near the center of the disc.
In addition, the heat conduction differential equation in the wedge-shaped area is:
d d r ( k 0 α r ˜ d T d r ) = α r ˜ q 0 ( 1 + 0.1 p 0.15 p r ˜ )
The corresponding boundary conditions are invoked as:
k 0 α r d T d r = q t i p ,   r = R L
T = T c ,   r = 0
Integrating Equations (12)–(14) yields the dimensionless temperature difference between T c and T 0
T c T 0 q 0 R 2 / k 0 = 1 120 ( 30 + 2 p L ˜ + p ) ( 1 L ˜ ) 2
Combining Equations (11) and (15) leads to the dimensionless temperature difference between T c and T min . Therefore, the dimensionless temperature difference T ˜ c can be defined as:
T ˜ c = ( T c T min ) / ( q 0 R 2 / k 0 )          = α 2 { λ k ˜ ( p 40 L ˜ 3 + 5 p 15 L ˜ 2 + p 20 20 L ˜ + 1 ) + 1 k ˜ α [ 1 p 10 L ˜ +                + p 40 ] L ˜ 2 + 1 α ( 30 + p 60 + p 30 L ˜ ) ( 1 + L ˜ 2 2 L ˜ ) }
From Equation (1), L ˜ can be inspected as:
L ˜ = ( α / 2 ϕ λ ) 1 / 2
In terms of the structure of the TE, the HCC is geometrically contained in it, which should meet the following constraint:
D < α ( R L )
As discussed previously, another way to express Equation (18) is:
[ 1 + 1 / ( α λ ) ] ( α / 2 ϕ λ ) 1 / 2 < 1
From Equation (3), the temperature difference T r T min should be compared with T c T min to determine the MTD ( T max T min ), i.e.:
T max T min = max ( T c T min , T r T min )
The heat conduction differential equation in the k 0 material of trapezoid area is:
d 2 T d y 2 + q 0 ( 1 + 0.1 p 0.15 p r ˜ ) k 0 = 0
The corresponding boundary conditions are invoked as:
d T d y = 0 ,   y = α r / 2
T = T min ,   y = 0
Integrating Equations (21)–(23), the temperature difference T r T min on the base of the TE can be quantified as:
T r T min = q 0 ( 1 0.05 p ) α 2 R 2 / ( 8 k 0 )
The corresponding T ˜ r is also specified as:
T ˜ r = ( T r T min ) / ( q 0 R 2 / k 0 ) = ( 1 0.05 p ) α 2 / 8
Hence, the dimensionless MTD Δ T ˜ of the heat conduction disc is the larger one between T ˜ c and T ˜ r , that is:
Δ T ˜ = max ( T ˜ c , T ˜ r )
In view of Equation (16), under given conditions that the apex angle α , area ratio ϕ and conductivity ratio k ˜ are fixed, T ˜ c can be described as a function of the dimensionless HCC length L ˜ . By this means, the optimal L ˜ ( L ˜ o p t ) can be derived by solving T ˜ c / ( L ˜ ) = 0 , which is implicitly given by:
1 k ˜ ϕ [ 15 p L ˜ o p t 4 + ( 160 32 p ) L ˜ o p t 3 + ( 360 + 18 p ) L ˜ o p t 2 + 240 L ˜ o p t ] + 1 k ˜ [ 18 p L ˜ o p t 2 + ( 120 + 3 p ) L ˜ o p t ] + 6 [ p L ˜ o p t 2 + ( 10 p ) L ˜ o p t 10 ] = 0
By Equations (16) and (17), T ˜ c , min is depicted as:
T ˜ c , min = 1 2 { 1 k ˜ ϕ [ p 20 L ˜ o p t 5 + 10 2 p 15 L ˜ o p t 4 + p 20 10 L ˜ o p t 3 + 2 L ˜ o p t 2 ] + 1 k ˜ [ p 10 L ˜ o p t 3 + p 40 40 L ˜ o p t 2 ] + p 30 L ˜ o p t 3 + 10 p 20 L ˜ o p t 2 L ˜ o p t + 30 + p 60 }
From Equation (27), the optimal HCC length ( L ˜ o p t ) is definitely independent of α , but merely related to ϕ and k ˜ . Figure 3 and Figure 4 depict the effects of the NUHG coefficient p and area ratio ϕ on the relationship of T ˜ c , min and k ˜ ϕ , respectively. From Figure 3, one can see that T ˜ c , min tends to be 0 when k ˜ ϕ > > 10 , i.e., T ˜ c , min is rather minor when k ˜ is far larger than 1, so that the HCP can be improved by increasing the conductivity ratio in this regard. Additionally, T ˜ c , min increases with the increase of p when k ˜ ϕ is given, and the HCP is weaken. From Figure 4, when k ˜ ϕ is far less than 1, T ˜ c , min is equal to 0.291. By using Equation (25), α needs to be 1.762 when T ˜ r , min = 0.291 , but this could be seen as a breach of the prior assumption of sufficiently small apex angle α . Hence, one can highlight that the hot spot is always situated at the center of the disc without exception.
Figure 5 and Figure 6 depict the effects of the NUHG coefficient p and area ratio ϕ on the relationship of L ˜ o p t and k ˜ ϕ , respectively. From Figure 5, when p is initially specified, L ˜ o p t increases in an S-shaped mode with the increase of k ˜ ϕ , and obviously appears to have a quick change between 1 and 10. When k ˜ ϕ < < 1 , L ˜ o p t tends to be 0, and the high conductivity material is nearly stacked on the rim of the disc in this condition. When k ˜ ϕ > > 10 , L ˜ o p t tends to be 1, and the HCC almost extends to the center of the disc. When k ˜ ϕ is given, L ˜ o p t increases as p increases. In this case, the larger the HGR near the center of the disc is, the longer the HCC is in order to deliver the generated heat and reduce the temperature difference efficiently. From Figure 6, when k ˜ ϕ is specified, L ˜ o p t increases as ϕ increases. This is because the HCC needs to be longer to deliver more heat current near the center of the disc.

3. Constructal Optimization of a Radial-Pattern Disc with Numerical Solution

When the presuppositions of ϕ < < 1 , k 0 < < k p and small α are untenable, the actual heat conduction in k 0 and k p materials are both no longer a simple one-dimensional one, and Equations (8), (12) and (21) are invalidated in this manner. This problem always exists in analytical solution, and one can solve it by using higher dimensional model with numerical solution. As shown in Figure 2, the two-dimensional heat conduction differential equations in the k 0 and k p materials are, respectively, given as:
2 T r 2 + 2 T y 2 + q 0 f ( r ) k 0 = 0
2 T r 2 + 2 T y 2 = 0
With a view of the symmetric geometry of the TE, there appears to be sufficient to merely investigate the area of y 0 . In this regard, the requisite boundary conditions are:
T = T min    r = R ,    0 < y D / 2
T / r = 0    r = R ,    D / 2 < y α R / 2
T / n = 0    { 0 < r R ,     y = 0 0 r < R ,     y = α r / 2
where n is termed as the normal vector to the corresponding border.
For simplification, it is assumed that there exists no thermal contact resistance at the interface between the k 0 and k p materials, so that the heat current continuity equation is:
k 0 ( T / n ) Disc = k p ( T / n ) HCC
The dimensionless temperature difference and geometry parameters are defined the same as those in Equations (4) and (5). From Equation (4), one can know that the dimensionless MTD Δ T ˜ 0 is
Δ T ˜ 0 = ( T max T min ) / ( q 0 R 2 / k 0 )
where Δ T ˜ 0 is the dimensionless MTD of the whole disc.
The PDE toolbox of Matlab is invoked to solve Equations (29)–(34), in which the finite element method is adopted. During the processes of numerical calculations, the target area is sequentially refined until the relative error of MTDs between Δ T ˜ 0 j and Δ T ˜ 0 j 1 obtained from the j th and j 1 th refinements satisfies:
| ( Δ T ˜ 0 j Δ T ˜ 0 j 1 ) / Δ T ˜ 0 j 1 | < 0.005
According to Equations (16) and (17), it clearly indicates that when the conductivity ratio k ˜ , area ratio ϕ and HCC number N ( = 2 π / α ) are specified, the dimensionless MTD Δ T ˜ 0 of the disc is only a function of the length ratio λ ( = L / D ). As a result, by taking the NUHG into consideration, the constructal optimization of the disc with constant and variable cross-sectional HCCs can be carried out with λ as the variable, respectively.

3.1. Constant Cross-Sectional HCCs

The constant cross-sectional HCCs model is the same as that shown in Figure 1 and Figure 2, and one can further optimize its HCP with using numerical solution as follows: Figure 7 depicts the effect of the NUHG coefficient p on the relationship of the dimensionless MTD Δ T ˜ 0 and length ratio λ . It is easy to see from Figure 7 that when N , ϕ , k ˜ and p are already given, there appears to be an optimal λ ( λ o p t ) leading to the minimum Δ T ˜ 0 ( Δ T ˜ 0 , min ), which means that choosing an appropriate λ can effectively reduce the MTD and improve the HCP to some extent. Furthermore, λ o p t increases with the increase of p , in this case that the HCC needs to protrude into the center of the disc where the heat is relatively concentrated, which aims to reduce the MTD. On the one hand, the nonuniform heat generating model can be processed back to the model mentioned in [46] in the condition of p = 0 . On the other hand, when p = 10 , k ˜ = 5000 , ϕ = 0.01 and N = 100 , Δ T ˜ 0 is equal to 0.2764 by using the finite element method, and comparatively Δ T ˜ 0 is equal to 0.2711 by using the analytical method in Equation (16). These preconditions mentioned above mainly meet the foundations of Equations (8), (12) and (21). Clearly, if given these prerequisites, the two-dimensional model can be simplified as a one-dimensional one. The relative error between these two results is merely 1.9%, which in the meantime endorses the validity of analytical method under the given assumptions. Besides, in case of p = 10 , k ˜ = 5 , ϕ = 0.1 and N = 25 , Δ T ˜ 0 is equal to 0.2686 and 0.3868 obtained from the finite element method and the analytical method, respectively, and the deviation between the two is severely 30.6%. In the context, the given conditions violate the establishing conditions of the analytical method, which means that one can accordingly choose a better way to calculate the MTD of the disc based on the given conditions.
Figure 8 and Figure 9 depict the effects of area ratio ϕ and conductivity ratio k ˜ on the relationship of the dimensionless MTD Δ T ˜ 0 and length ratio λ . It is easy to see from Figure 8 that when N , k ˜ and p are already given, the minimum dimensionless MTD Δ T ˜ 0 , min decreases with the increase of ϕ as well as the corresponding λ o p t does. This indicates that whereas more high conductivity material is available, the HCC seems to be thicker and shorter for a better HCP. From Figure 9, when N , ϕ and p are already given, both Δ T ˜ 0 , min and λ o p t decrease with the increase of k ˜ , which means that the HCP can be improved by increasing the conductivity ratio.
Figure 10 depicts the effect of the NUHG coefficient ( p ) on the relationship of the minimum dimensionless MTD Δ T ˜ 0 , min and HCC number N . From Figure 10, when p is equal to 0 and 15, Δ T ˜ 0 , min tends to be 0.0117 and 0.0159, respectively, with a large value of N . Apparently, the MTD can be reduced by increasing the HCC number on condition that N varies from 10 to 25. When p is equal to 10, Δ T ˜ 0 , min of the disc with 25 HCCs is specifically reduced by 48.8% compared to that with 10 HCCs. Instead, Δ T ˜ 0 , min remains constant when N varies from 25 to 40. Whereupon, one can speculate that it is costly to further increase the number of HCCs when N > 25 , and the HCP is not improved to some extent in this case.

3.2. Variable Cross-Sectional HCCs

In order to further improve the HCP of the disc, the model of a TE with variable cross-sectional [28,29] HCCs is established as shown in Figure 11. The minimum and maximum widths of the HCCs are D and m D , respectively, where m is termed as the width ratio. The area ratio ϕ of k p material is:
ϕ = D L N ( m + 1 ) / ( 2 π R 2 )
Similarly, the dimensionless L ˜ is:
L ˜ = [ α / ( m + 1 ) λ ϕ ] 1 / 2
When the width ratio m is specified, constructal optimizations of the disc with variable cross-sectional HCCs can be carried out similar to that in Section 3.1. By invoking the PDE toolbox, one can numerically obtain the MTD in the heat generating area according to Equations (37) and (38) and the numerical method given in Section 3.1.
Figure 12 depicts the effect of the width ratio m on the relationship of the minimum dimensionless MTD Δ T ˜ 0 , min and HCC number N with ϕ = 0.05 and p = 10 . From Figure 12, it is clear that Δ T ˜ 0 , min first decreases and then remains a constant with the increase of N . Otherwise, when N is given, Δ T ˜ 0 , min gradually decreases as m increases. According to [46], when the disc is uniformly heat generating and inserted with constant cross-sectional HCCs, a.k.a. p = 0 , Δ T ˜ 0 , min is equal to 0.0107 in terms of k ˜ = 500 , ϕ = 0.1 and N = 20 . While applying the variable cross-sectional structure in the disc under the same constraint conditions, Δ T ˜ 0 , min is equal to 0.0097 which is reduced by 9.3%. In these two cases, the optimal constructs of HCCs both have no geometrical connections with the center of the disc. From Figure 12, one can see that as the disc consists of 15 TEs and the heat is nonuniformly generated, Δ T ˜ 0 , min is equal to 0.0174 in case of m = 1 while Δ T ˜ 0 , min is equal to 0.0152 in case of m = 4 , that is, Δ T ˜ 0 , min is decreased by 12.6% when m changes from 1 to 4. In these two nonuniform heat generating cases, the optimal constructs of HCCs likewise have no geometrical connections with the center. As revealed in Figure 13a,b the isothermal lines are labelled with corresponding dimensionless temperature differences. In the premise of ϕ = 0.05 , p = 10 and N = 25 , for the disc with constant cross-sectional HCCs (i.e., m = 1 ), the optimal length ratio ( λ o p t ) turns out to be 151 resulting in Δ T ˜ 0 , min , which is equal to 0.0147. Besides, for the disc with variable cross-sectional HCCs (i.e., m = 4 ), λ o p t turns out to be 385 resulting in Δ T ˜ 0 , min , which is equal to 0.0125. It can be seen that Δ T ˜ 0 , min of the disc with m = 4 is decreased by 15.0% compared to that with m = 1 . When the disc is nonuniformly heat generating, heat is relatively concentrated near the heat sink, i.e., the ends of HCCs situated on the perimeter of the disc. It is beneficial to place more high conductivity material near the heat sink so as to efficiently transfer the heat current. Therefore, the geometry of variable cross-sectional HCCs takes advantage of distributing more high conductivity material near the region where heat is more concentrated, which can be applied in the constructal design of radial-pattern disc to improve its HCP.

4. Conclusions

With the high demand to efficient cooling systems for heat generating area, the heat conduction model of a radial-pattern disc considering the NUHG is established in this paper, which is closer to real problems in the field of electronics cooling. Constructal optimizations of the discs with constant and variable cross-sectional HCCs are implemented, and their MTDs are minimized. The influences of the NUHG coefficient p , HCC number N and conductivity ratio k ˜ on the optimal results are analyzed. The HCP comparisons obtained by analytical and numerical solutions as well as constant and variable cross-sectional HCCs are carried out. The results are given as follows:
(1)
For the disc with constant cross-sectional HCCs, the optimal length ratio λ o p t increases with the increase of p , in this case that the HCCs need to protrude into the center of the disc where the heat is relatively concentrated. The deviation of the optimal constructs obtained from analytical method and finite element method is comparatively slight. When p is equal to 10, Δ T ˜ 0 , min of the disc with 25 HCCs is specifically reduced by 48.8% compared to that with 10 HCCs. As a result, the HCP of the disc can be efficiently improved by properly increasing the number of HCCs. In addition, the approach of increasing k ˜ can also contribute to reduce the MTD and improve the HCP of the disc.
(2)
For the disc with variable cross-sectional HCCs, the Δ T ˜ 0 , min of the disc in case of m = 4 is decreased by 15.0% compared to that of m = 1 in the premise of ϕ = 0.05 , p = 10 and N = 25 . Therefore, the geometry of variable cross-sectional HCCs takes advantage of distributing more high conductivity material near the region where heat is more concentrated, which can be applied in the constructal design of radial-pattern disc to improve its HCP.
In this paper, the heat conduction problem of the radial pattern disc with a single high conductivity channel arranged in each element is studied. In the following research, the radial-pattern disc with more complex multistage tree-shaped high conductivity channels will be investigated wishing to further improve the heat transfer performance. Furthermore, the effects of the location where the heat sink is situated on the heat conduction performance can be studied. Additionally, the percentage of the surface/volume in which the temperature is above a certain value is worth taking into consideration to improve the existing work.

Author Contributions

J.Y., H.F., L.C., Z.X. commonly finished the manuscript. All authors have read and approved the final manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (Grant Nos. 51506220 and 51579244) and the Natural Science Foundation of Hubei Province (Grant No. 2016CFB504). The authors wish to thank the reviewers for their careful, unbiased and constructive suggestions, which led to this revised manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclatures

D Width [m]
f ( r ) Heat generating rate function
k 0 Thermal conductivity of heat generating area [W/(m·K)]
k p Thermal conductivity of high conductivity material [W/(m·K)]
L Length [m]
m Width ratio
N Number of high conductivity channels
p Nonuniform heat generating coefficient
q 0 Heat generating constant in per unit volume [W/m3]
q r Heat current collected at the longitudinal direction of y axis [W/m2]
q t i p Heat current q t i p generated in the wedge-shaped area [W/m]
R Radius [m]
r Horizontal axis
T Temperature [K]
y Vertical axis
Greek symbols
α Apex angle of each sectorial element
λ Length ratio L / D
ϕ Area ratio of high conductivity material area to whole disc area
Subscripts
cCenter of disc
maxMaximum
minMinimum
optOptimum
rRim of disc
superscript
~nondimensionalized

Abbreviations

The following abbreviations are used in this manuscript:
HCCHigh conductivity channel
HCPHeat conduction performance
HGRHeat generating rate
MTDMaximum temperature difference
NUHGNonuniform heat generation

References

  1. Li, B.T.; Hong, J.; Tian, X.Y. Generating optimal topologies for heat conduction by heat flow paths identification. Int. Commun. Heat Mass Transfer 2016, 75, 177–182. [Google Scholar] [CrossRef]
  2. Dbouk, T. A review about the engineering design of optimal heat transfer systems using topology optimization. Appl. Therm. Eng. 2017, 112, 841–854. [Google Scholar] [CrossRef]
  3. Manuel, M.C.E.; Lin, P.T. Design explorations of heat conductive pathways. Int. J. Heat Mass Transfer 2017, 104, 835–851. [Google Scholar] [CrossRef]
  4. Bejan, A. Shape and Structure, from Engineering to Nature; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  5. Bejan, A.; Lorente, S. Constructal theory of generation of configuration in nature and engineering. J. Appl. Phys. 2006, 100, 041301. [Google Scholar] [CrossRef]
  6. Bejan, A.; Lorente, S. Design with Constructal Theory; Wiley: Hoboken, NJ, USA, 2008. [Google Scholar]
  7. Ding, X.; Yamazaki, K. Constructal design of cooling channel in heat transfer system by utilizing optimality of branch systems in nature. J. Heat Transfer 2007, 129, 245–255. [Google Scholar] [CrossRef]
  8. Lorenzini, G.; Moretti, S.; Conti, A. Fin Shape Thermal Optimization Using Bejan’s Constructal Theory; Morgan & Claypool Publishers: San Francisco, CA, USA, 2011; pp. 1–219. [Google Scholar]
  9. Hajmohammadi, M.R.; Poozesh, S.; Hosseini, R. Radiation effect on constructal design analysis of a T-Y-shaped assembly of fins. J. Therm. Sci. Technol. 2012, 7, 677–692. [Google Scholar] [CrossRef]
  10. Chen, L.G. Progress in study on constructal theory and its application. Sci. China Technol. Sci. 2012, 55, 802–820. [Google Scholar] [CrossRef]
  11. Bejan, A.; Lorente, S. Constructal law of design and evolution: Physics, biology, technology, and society. J. Appl. Phys. 2013, 113, 151301. [Google Scholar] [CrossRef]
  12. Hajmohammadi, M.R.; Poozesh, S.; Campo, A.; Nourazar, S.S. Valuable reconsideration in the constructal design of cavities. Energy Convers. Mgmt. 2013, 66, 33–40. [Google Scholar] [CrossRef]
  13. Hajmohammadi, M.R.; Poozesh, S.; Rahmani, M.; Campo, A. Heat transfer improvement due to the imposition of non-uniform wall heating for in-tube laminar forced convection. Appl. Therm. Eng. 2013, 61, 268–277. [Google Scholar] [CrossRef]
  14. Yang, J.; Oh, S.R.; Liu, W. Optimization of shell-and-tube heat exchangers using a general design approach motivated by constructal theory. Int. J. Heat Mass Transfer 2014, 77, 1144–1154. [Google Scholar] [CrossRef]
  15. Bejan, A. Constructal law: optimization as design evolution. Trans. ASME J. Heat Transfer 2015, 137, 061003. [Google Scholar] [CrossRef]
  16. Xie, G.N.; Shen, H.; Wang, C.C. Parametric study on thermal performance of microchannel heat sinks with internal vertical Y-shaped bifurcations. Int. J. Heat Mass Transfer 2015, 90, 948–958. [Google Scholar] [CrossRef]
  17. Zhang, F.L.; Sunden, B.; Zhang, W.H.; Xie, G.N. Constructal parallel-flow and counter-flow microchannel heat sinks with bifurcations. Numer. Heat Transfer A-Appl. 2015, 68, 1087–1105. [Google Scholar] [CrossRef]
  18. Zhang, R.P.; Chen, Z.Y.; Xie, G.N.; Sunden, B. Numerical analysis of constructal water-cooled microchannel heat sinks with multiple bifurcations in the entrance region. Numer. Heat Transfer A-Appl. 2015, 67, 632–650. [Google Scholar] [CrossRef]
  19. Chen, L.G.; Feng, H.J.; Xie, Z.H. Generalized thermodynamic optimization for iron and steel production processes: Theoretical exploration and application cases. Entropy 2016, 18, 353. [Google Scholar] [CrossRef]
  20. Bejan, A.; Errera, M.R. Complexity, organization, evolution, and constructal law. J. Appl. Phys. 2016, 119, 074901. [Google Scholar] [CrossRef] [Green Version]
  21. Hajmohammadi, M.R. Introducing a ψ-shaped cavity for cooling a heat generating medium. Int. J. Therm. Sci. 2017, 121, 204–212. [Google Scholar] [CrossRef]
  22. Feng, H.J.; Chen, L.G.; Xie, Z.H. Multi-disciplinary, multi-objective and multi-scale constructal optimizations for heat and mass transfer processes performed in Naval University of Engineering, a review. Int. J. Heat Mass Transfer 2017, 115, 86–98. [Google Scholar] [CrossRef]
  23. Chen, L.G.; Xiao, Q.H.; Feng, H.J. Constructal optimizations for heat and mass transfers based on the entransy dissipation extremum principle, performed at the Naval University of Engineering; a review. Entropy 2018, 20, 74. [Google Scholar] [CrossRef]
  24. Miguel, A.F.; Rocha, L.A.O. Tree-Shaped Fluid Flow and Heat Transfer; Springer: New York, NY, USA, 2018. [Google Scholar]
  25. Chen, L.G.; Feng, H.J.; Xie, Z.H. Progress of constructal theory in China over the past decade. Int. J. Heat Mass Transfer 2017, 55. [Google Scholar] [CrossRef]
  26. Bejan, A. Constructal law, twenty years after. Proc. Rom. Acad. Ser. A Math. Phys. Tech. Sci. Inf. Sci. 2018, 18, 309–311. [Google Scholar]
  27. Bejan, A. Constructal-theory network of conducting path for cooling a heat generating volume. Int. J. Heat Mass Transfer 1997, 40, 799–816. [Google Scholar] [CrossRef]
  28. Neagu, M.; Bejan, A. Constructal-theory tree networks of “constant” thermal resistance. J. Appl. Phys. 1999, 86, 1136–1144. [Google Scholar] [CrossRef]
  29. Neagu, M.; Bejan, A. Three-dimensional tree constructs of “constant” thermal resistance. J. Appl. Phys. 1999, 86, 7107–7115. [Google Scholar] [CrossRef]
  30. Ghodoossi, L.; Egrican, N. Exact solution for cooling of electronics using constructal theory. J. Appl. Phys. 2003, 93, 4922–4929. [Google Scholar] [CrossRef]
  31. Wu, W.J.; Chen, L.G.; Sun, F.R. On the “area to point” flow problem based on constructal theory. Energy Convers. Mgmt. 2007, 48, 101–105. [Google Scholar] [CrossRef]
  32. Wei, S.H.; Chen, L.G.; Sun, F.R. The volume-point constructal optimization for discrete variable cross-section conducting path. Appl. Energy 2009, 86, 1111–1118. [Google Scholar] [CrossRef]
  33. Lorenzini, G.; Biserni, C.; Dalpiaz, F.L.; Fagundes, T.M.; Rocha, L.A.O. Geometric optimization of T-shaped constructs coupled with a heat generating basement: A numerical approach motivated by Bejan’s constructal theory. J. Eng. Thermophys. 2017, 26, 485–497. [Google Scholar] [CrossRef]
  34. Lorenzini, G.; Biserni, C.; Rocha, L.A.O. Constructal design of X-shaped conductive pathways for cooling a heat-generating body. Int. J. Heat Mass Transfer 2013, 58, 513–520. [Google Scholar] [CrossRef]
  35. Lorenzini, G.; Biserni, C.; Rocha, L.A.O. Constructal design of non-uniform X-shaped conductive pathways for cooling. Int. J. Therm. Sci. 2013, 71, 140–147. [Google Scholar] [CrossRef]
  36. Hajmohammadi, M.R. Phi and Psi shaped conductive routes for improved cooling in a heat generating piece. Int. J. Therm. Sci. 2014, 77, 66–74. [Google Scholar] [CrossRef]
  37. Hajmohammadi, M.R.; Lorenzini, G.; Shariatzadeh, O.J.; Biserni, C. Evolution in the design of V-shaped highly conductive pathways embedded in a heat-generating piece. Trans. ASME J. Heat Transfer 2015, 137, 061001. [Google Scholar] [CrossRef]
  38. Feng, H.J.; Chen, L.G.; Xie, Z.H.; Sun, F.R. Constructal design for “+” shaped high conductive pathways over a square body. Int. J. Heat Mass Transfer 2015, 91, 162–169. [Google Scholar] [CrossRef]
  39. Lorenzini, G.; Barreto, E.X.; Beckel, C.C.; Schneider, P.S. Constructal design of I-shaped high conductive pathway for cooling a heat-generating medium considering the thermal contact resistance. Int. J. Heat Mass Transfer 2016, 93, 770–777. [Google Scholar] [CrossRef]
  40. Beckel, C.C.; Isoldi, L.A.; Santos, E.D.D.; Rocha, L.A.O. Constructal design of non-uniform T-shaped conductive pathways for cooling heat generating bodies. Vetor Rio Grande 2016, 26, 2–13. [Google Scholar]
  41. Lorenzini, G.; Barreto, E.X.; Beckel, C.C.; Schneider, P.S.; Isoldi, L.A.; Santos, E.D.; Rocha, L.A.O. Geometrical evaluation of T-shaped high conductive pathway with thermal contact resistance for cooling of heat-generating medium. Int. J. Heat Mass Transfer 2017, 108, 1884–1893. [Google Scholar] [CrossRef]
  42. Ghodoossi, L.; Egrican, N. Conductive cooling of triangular shaped electronics using constructal theory. Energy Convers. Manag. 2004, 45, 811–828. [Google Scholar] [CrossRef]
  43. Feng, H.J.; Chen, L.G.; Xiao, Q.H.; Sun, F.R. Global constructal optimization design for triangular assembly at micro and nanoscales. CIESC J. 2012, 63, 47–53. (In Chinese) [Google Scholar]
  44. Chen, L.G.; Wu, W.J.; Sun, F.R. Constructal re-optimization of heat conduction with the triangular elemental area. Int. J. Low-Carbon Technol. 2014, 9, 256–261. [Google Scholar] [CrossRef]
  45. Feng, H.J.; Chen, L.G.; Xie, Z.H.; Sun, F.R. “Volume-point” heat conduction constructal optimization based on minimization of maximum thermal resistance with triangular element at micro and nanoscales. J. Energy Inst. 2016, 89, 302–312. [Google Scholar] [CrossRef]
  46. Silva, A.K.D.; Vasile, C.; Bejan, A. Disc cooled with high-conductivity inserts that extend inward from the perimeter. Int. J. Heat Mass Transfer 2004, 47, 4257–4263. [Google Scholar] [CrossRef]
  47. Sharifi, F.; Ghaedmini, H.; Salimpour, M.R. Using incomplete variable cross-section highly conductive inserts for cooling a disc. Front. Heat Mass Transfer 2012, 3, 043005. [Google Scholar]
  48. Rocha, L.A.O.; Lorente, S.; Bejan, A. Constructal design for cooling a disc-shaped area by conduction. Int. J. Heat Mass Transfer 2002, 45, 1643–1652. [Google Scholar] [CrossRef]
  49. Rocha, L.A.O.; Lorente, S.; Bejan, A. Conduction tree networks with loops for cooling a heat generating volume. Int. J. Heat Mass Transfer 2006, 49, 2626–2635. [Google Scholar] [CrossRef]
  50. Xiao, Q.H.; Chen, L.G.; Sun, F.R. Constructal optimization for “disc-to-point” heat conduction without the premise of optimized last-order construct. Int. J. Therm. Sci. 2011, 50, 1031–1036. [Google Scholar] [CrossRef]
  51. Bahadormanesh, N.; Salimpour, M.R. Constructal design of high-emissivity radiation inserts embedded in a disc-shaped heat generation body. Appl. Therm. Eng. 2016, 112, 638–648. [Google Scholar] [CrossRef]
  52. Sharfi, F.; Salimpuor, R.M.; Campo, A. Cooling a solid disc with uniform heat generation using inserts of high thermal conductivity within the constructal design platform. Int. J. Therm. Environ. Eng. 2016, 12, 15–26. [Google Scholar]
  53. Chen, L.G.; Feng, H.J.; Xie, Z.H.; Sun, F.R. Constructal optimization for “disc-point” heat conduction at micro and nanoscales. Int. J. Heat Mass Transfer 2013, 67, 704–711. [Google Scholar] [CrossRef]
  54. Daneshi, M.; Zare, M.; Salimpour, M.R. Micro and nano-scale conductive tree-structures for cooling a disc-shaped electronic piece. Trans. ASME J. Heat Transfer 2013, 135, 031401. [Google Scholar] [CrossRef]
  55. Cetkin, E.; Oliani, A. The natural emergence of asymmetric tree-shaped pathways for cooling of a non-uniformly heated domain. J. Appl. Phys. 2015, 118, 024902. [Google Scholar] [CrossRef] [Green Version]
  56. Feng, H.J.; Chen, L.G.; Xie, Z.H.; Sun, F.R. Constructal design for a rectangular body with nonuniform heat generation. Eur. Phys. J. Plus 2016, 131, 274. [Google Scholar] [CrossRef]
  57. Feng, H.J.; Chen, L.G.; Xie, Z.H.; Sun, F.R. Constructal entransy dissipation rate minimization of a rectangular body with nonuniform heat generation. Sci. China Technol. Sci. 2016, 59, 1352–1359. [Google Scholar] [CrossRef]
  58. You, J.; Feng, H.J.; Chen, L.G.; Xie, Z.H. Heat conduction constructal optimization for nonuniform heat generating area based on triangular element. Int. J. Heat Mass Transfer 2018, 117, 896–902. [Google Scholar] [CrossRef]
Figure 1. Heat conductive model of a radial-pattern disc with nonuniform heat generating.
Figure 1. Heat conductive model of a radial-pattern disc with nonuniform heat generating.
Entropy 20 00685 g001
Figure 2. Simplified triangular element with nonuniform heat generation and constant cross-sectional high conductivity channel.
Figure 2. Simplified triangular element with nonuniform heat generation and constant cross-sectional high conductivity channel.
Entropy 20 00685 g002
Figure 3. Effect of p on T ˜ c , min versus k ˜ ϕ .
Figure 3. Effect of p on T ˜ c , min versus k ˜ ϕ .
Entropy 20 00685 g003
Figure 4. Effect of ϕ on T ˜ c , min versus k ˜ ϕ .
Figure 4. Effect of ϕ on T ˜ c , min versus k ˜ ϕ .
Entropy 20 00685 g004
Figure 5. Effect of p on L ˜ o p t versus k ˜ ϕ .
Figure 5. Effect of p on L ˜ o p t versus k ˜ ϕ .
Entropy 20 00685 g005
Figure 6. Effect of ϕ on L ˜ o p t versus k ˜ ϕ .
Figure 6. Effect of ϕ on L ˜ o p t versus k ˜ ϕ .
Entropy 20 00685 g006
Figure 7. Effect of p on Δ T ˜ 0 versus λ .
Figure 7. Effect of p on Δ T ˜ 0 versus λ .
Entropy 20 00685 g007
Figure 8. Effect of ϕ on Δ T ˜ 0 versus λ .
Figure 8. Effect of ϕ on Δ T ˜ 0 versus λ .
Entropy 20 00685 g008
Figure 9. Effect of k ˜ on Δ T ˜ 0 versus λ .
Figure 9. Effect of k ˜ on Δ T ˜ 0 versus λ .
Entropy 20 00685 g009
Figure 10. Effect of p on Δ T ˜ 0 , min versus N .
Figure 10. Effect of p on Δ T ˜ 0 , min versus N .
Entropy 20 00685 g010
Figure 11. Simplified triangular element with nonuniform heat generating and variable cross-sectional high conductivity channel.
Figure 11. Simplified triangular element with nonuniform heat generating and variable cross-sectional high conductivity channel.
Entropy 20 00685 g011
Figure 12. Effect of m on Δ T ˜ 0 , min versus N .
Figure 12. Effect of m on Δ T ˜ 0 , min versus N .
Entropy 20 00685 g012
Figure 13. Optimal construct comparison of triangular elements with different high conductivity channels: (a) constant cross-section, (b) variable cross-section.
Figure 13. Optimal construct comparison of triangular elements with different high conductivity channels: (a) constant cross-section, (b) variable cross-section.
Entropy 20 00685 g013

Share and Cite

MDPI and ACS Style

You, J.; Feng, H.; Chen, L.; Xie, Z. Constructal Optimization for Cooling a Non-Uniform Heat Generating Radial-Pattern Disc by Conduction. Entropy 2018, 20, 685. https://doi.org/10.3390/e20090685

AMA Style

You J, Feng H, Chen L, Xie Z. Constructal Optimization for Cooling a Non-Uniform Heat Generating Radial-Pattern Disc by Conduction. Entropy. 2018; 20(9):685. https://doi.org/10.3390/e20090685

Chicago/Turabian Style

You, Jiang, Huijun Feng, Lingen Chen, and Zhihui Xie. 2018. "Constructal Optimization for Cooling a Non-Uniform Heat Generating Radial-Pattern Disc by Conduction" Entropy 20, no. 9: 685. https://doi.org/10.3390/e20090685

APA Style

You, J., Feng, H., Chen, L., & Xie, Z. (2018). Constructal Optimization for Cooling a Non-Uniform Heat Generating Radial-Pattern Disc by Conduction. Entropy, 20(9), 685. https://doi.org/10.3390/e20090685

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