1556-276X-6-240 1556-276X Nano Express <p>Constructal blade shape in nanofluids</p> BaiChaobaichao@hku.hk WangLiqiulqwang@hku.hk

Department of Mechanical Engineering, The University of Hong Kong, Pokfulam Road, Hong Kong

 Nanoscale Research Letters 1556-276X 2011 6 1 240 http://www.nanoscalereslett.com/content/6/1/240 2171175110.1186/1556-276X-6-240 6122010213201121320112011Bai and Wang; licensee Springer.This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Blade configuration of nanofluids has been proven to perform much better than dispersed configuration for some heat conduction systems. The analytical analysis and numerical calculation are made for the cylinder--shaped and regular-rectangular-prism--shaped building blocks of the blade-configured heat conduction systems (using nanofluids as the heat conduction media) to find the optimal cross-sectional shape for the nanoparticle blade under the same composing materials, composition ratio, volumetric heat generation rate, and total building block volume. The regular-triangular-prism--shaped blade has been proven to perform better than all the other three kinds of blades, namely, the regular-rectangular-prism--shaped blade, the regular-hexagonal-prism--shaped blade, and the cylinder--shaped blade. Thus, the regular-triangular-prism--shaped blade is selected as the optimally shaped blade for the two kinds of building blocks that are considered in this study. It is also proven that the constructal cylinder--regular-triangular-prism building block performs better than the constructal regular-rectangular-prism--regular-triangular-prism building block.

 NanofluidsYogesh Jaluria, Liqiu Wang, Stephen Choi and Oronzio Manca

Introduction

Nanofluids are mixtures of nanoparticles and base fluids, which have different thermal conductivities 12345. They were identified and proposed as a result of people's persistent pursuit for more and more efficient heat-transfer media. It should be noted that conventional heat transfer fluids have normally very low thermal conductivity, thus destroying much exergy during heat transport. At present, a great amount of attention is paid to studies on nanofluids, with the aim of addressing many unsolved issues 678910.

Constructal theory is a novel thought for nature and society 11121314151617. It tries to explain phenomena based on optimization, or natural selection in biological terms. One of its viewpoints is that two flow mechanisms are better than one 18. This is what one sees naturally in river basins, lung structure, and the percolation threshold effect 19. Our previous studies have proved that the blade configuration of nanofluids is much better than the dispersed configuration for the two kinds of disk-shaped heat conduction systems with different boundary conditions 2021. It is believed that the continuous nanoparticle blades with higher thermal conductivity serve as the second conduction mechanism, and its optimized cooperation with the base fluid of low thermal conductivity leads to the much better performance. In this study, the blade configuration of nanofluids is considered in detail by studying the influence of the shapes of high-conductivity blades in two kinds of building blocks of the total blade-configured heat conduction systems. This study is inspired by the need for the optimization of the cross section of duct for minimum flow resistance 11. By treating heat as a flow medium flowing in blades, one can also find the optimal shape for the blades, which offers minimum thermal resistance.

Optimal blade shape for two kinds of building blocks of blade-configured heat conduction systems

The cylinder--shaped and regular-rectangular-prism--shaped building blocks are studied in this article. The nanoparticle blade has four different shapes: regular triangular prism, regular rectangular prism, regular hexagonal prism, and cylinder. For conciseness, the word "regular" is omitted hereafter. A format of "building-block-shape--blade-shape" is used to indicate the eight kinds of building blocks. Figure 1 shows the two kinds of building blocks with cylinder blades. For all the eight kinds of building blocks, uniform heat generation rate occurs in the base fluid region. All the external surfaces, except the cross-sectional plane x = 0 of the blade, are adiabatic. The cross-sectional plane x = 0 of the blade serves as the heat sink with constant temperature. The composition of these two-material building blocks is fixed by volume fraction

<p>Figure 1</p>

Two kinds of heat conduction building blocks considered in this study: (a) cylinder--cylinder building block; (b) rectangular-prism--cylinder building block

Two kinds of heat conduction building blocks considered in this study: (a) cylinder--cylinder building block; (b) rectangular-prism--cylinder building block.

It is assumed that ϕ ≪ 1, and the thermal conductivity ratio of nanoparticle material and base fluid material is fixed and large. The thermal contact resistance is not considered.

In order to study the influence of the blade shape, the materials of the base fluid and nanoparticle, volumetric heat generation rate, and volumes of the eight kinds of building blocks are also fixed, besides the volume fraction and thermal conductivity ratio; however, the slenderness is free to vary to achieve the constructal system (building block) overall temperature difference. Here, the slenderness refers to the ratio of the radius to length for the cylinder building block, and the ratio of the circumscribing cylinder radius to length for the rectangular prism building block. For the simplest cylinder--cylinder building block, analytical analysis can be made; the system overall temperature difference, the constructal system overall temperature difference, and the constructal slenderness can be obtained analytically. Based on a slenderness range predicted by the analytic result, the numerical calculation is then conducted for all the eight kinds of systems to obtain, as accurately as possible, the results for comparison among different blade shapes and different building block shapes.

Analytical analysis for cylinder--cylinder building block

Owing to the much higher thermal conductivity of nanoparticle material, heat conduction inside this kind of building block can be considered to consist of two one-dimensional routes-radial conduction inside cross-sectional planes of the base fluid region and axial conduction along the blade.

For the radial heat conduction inside base fluid, the governing equation and boundary conditions are

and

respectively, where, kf is the thermal conductivity of the base fluid; q''' is the volumetric heat generation rate; Tc is the temperature at the interface of the blade and base fluid at x = L0 (L0 is the length of the building block); r0 is the radius of the inner blade; and R0 is the outer radius of the building block. By solving Equations (2) and (3), the radial temperature distribution of the cross-sectional plane x = L0 of the base fluid region is obtained:

For the axial conduction inside the blade, the governing equation and boundary conditions, respectively, are

and

where kp is the thermal conductivity of the nanoparticle material, and T0 is the heat-sink temperature at the cross-sectional plane x = 0 of the blade. By solving Equations (5) and (6), the axial temperature distribution along the blade is obtained as

Thus, the overall temperature difference of the building block is

Nondimensionalizing this overall temperature difference with (constant), one has

Where is the nondimensional system (building block) overall temperature difference, and is the ratio of the thermal conductivities of the nanoparticle material and the base fluid material. By substituting Equation (1) in Equation (9), becomes

which indicates that the building block's overall temperature difference (or thermal resistance) depends on its slenderness under the same composing materials, composition ratio, volumetric heat generation rate, and total volume. Figure 2 typifies this dependence at ϕ = 0.05 and = 641.6667 (thermal conductivity ratio of copper and water). By minimizing this nondimensional system overall temperature difference with respect to , the nondimensional constructal system overall temperature difference can be obtained:

<p>Figure 2</p>

Variation of the nondimensional system overall temperature difference with slenderness for cylinder--cylinder building block (analytical result, ϕ = 0.05 and = 641.6667)

Variation of the nondimensional system overall temperature difference with slenderness for cylinder--cylinder building block (analytical result, ϕ = 0.05 and = 641.6667).

and

At , the best-performing cylinder--cylinder building block can be obtained. If one specifies ϕ = 0.05 and = 641.6667, then the optimal slenderness will be 0.240618, and the nondimensional constructal system overall temperature difference will be 0.296778, which is the lowest point in Figure 2. Similarly, for the other kinds of building blocks considered here, there also exists such a best-performing slenderness .

Numerical calculation for all the eight kinds of building blocks

The actual heat conduction in the cylinder--cylinder building block is of course not a simple combination of two one-dimensional conductions. For the other kinds of building blocks, the flow of heat is even more complex. In order to have as accurately as possible results for the comparison, a finite volume computational fluid dynamics (CFD) code 22 is used for obtaining numerical results for all the eight kinds of building blocks.

The conservation of energy equations are

and

for the base fluid and the blade, respectively. Here,

For rectangular-prism--series building blocks, R0 stands for the radius of the circumscribing cylinder, for which one has (where a is the side length of the rectangular cross section). Constant, , is introduced to use fewer grids to achieve accurate enough results. The "+1" in the nondimensional temperature expression is introduced to ensure that the computing process does not touch the limit of 0 K. At the interface between the kf region and kp blade, the continuity of heat flux requires that

where ñ is the nondimensional normal vector; . As all the external surfaces, except the plane x = 0 of the kp blade, are adiabatic,

and for the plane x = 0 of the kp blade, one has

Note that is exactly the nondimensional system overall temperature difference, , as shown in Equation (9), where is the maximal nondimensional temperature in the heat conduction building blocks.

It is specified that ϕ = 0.05 and = 641.6667 for the numerical calculation, the value used in the analytical analysis. The finite volume CFD code is chosen because of its efficiency and flexibility to generate a large number of results for various geometries which differ slightly from each other. Hexahedron grids are used to mesh the cylinder--hexagonal-prism and cylinder--cylinder building blocks, whereas all the other building blocks are meshed with wedge grids. Appropriate grid number is determined by doubling the interval number in , , and directions each time, until the change of temperature becomes less than 0.05% (the maximal temperature, , is used specifically for this criterion). Table 1 shows an example of how this grid independence is reached.

<p>Table 1</p>

Grid-independence check (cylinder--triangular-prism building block, R0/L0 = 0.25)

Number of grids

Changing of

2500

0.271978

-0.010181 × 100

20120

0.269209

-8.803569 × 10-4

167360

0.268972

-3.346073 × 10-5

1248160

0.268963

The nondimensional system overall temperature difference is shown in Figures 3 and 4 for the eight kinds of building blocks. The optimal slenderness and nondimensional constructal system overall temperature difference for the cylinder--cylinder building block are 0.25 and 0.289577, respectively. Comparing with the approximate analytical results, the differences are only 3.75% and 2.49%, which confirms the accuracy of the finite volume CFD code.

<p>Figure 3</p>

Numerical results for the cylinder--series building blocks and the analytical result for the cylinder--cylinder building block (ϕ = 0.05 and = 641.6667)

Numerical results for the cylinder--series building blocks and the analytical result for the cylinder--cylinder building block (ϕ = 0.05 and = 641.6667).

 <p>Figure 4</p>

Numerical results for the rectangular-prism--series building blocks (ϕ = 0.05 and = 641.6667)

Numerical results for the rectangular-prism--series building blocks (ϕ = 0.05 and = 641.6667).

All the eight kinds of heat conduction building blocks show strong performance dependence on the slenderness of the building blocks. Under the specific values for ϕ and , the optimal slendernesses for the cylinder--series and rectangular-prism--series building blocks are 0.25 and 0.3, respectively, within a resolution of 0.05. Below the optimal slenderness, the performance dependence on increases in the order from cylinder--cylinder building block to cylinder--hexagonal-prism building block, cylinder--rectangular-prism building block, and then cylinder--triangular-prism building block. Above the optimal slenderness, this trend changes its direction. A similar situation happens for the rectangular-prism--series building blocks. When the value is small enough ( for both the cylinder--series and rectangular-prism--series building blocks), the blade shape begins to have very weak effect on the system's performance, due to the diminishing role played by the building block's cross-sectional area when compared to its length. It can be seen from both Figures 3 and 4 that the four nondimensional system overall temperature difference curves almost collapse into one curve. For the cylinder--series building blocks, the collapsing curves approach the analytical result of the cylinder--cylinder building block as expected since the analytical result becomes more and more accurate as the value decreases. The analytical result shown in Figure 3 has the same ϕ and values as those for the numerical calculation (the same curve as Figure 2). Thus, the accuracy of the finite volume CFD code is verified again. It should be noted that both the cylinder--series and rectangular-prism--series building blocks have fixed optimal slenderness. Therefore, the triangular-prism--shaped blade always performs the best among the four kinds of blades considered. Besides, the performance difference between successively shaped blades decreases from triangular-prism blade to cylinder blade, which is consistent with their surface area changing under the same volume, as shown in Table 2. Here, the surface area of the blade is exactly the interfacial area between the base fluid and the blade. Thus, it is very likely that the performance difference can be attributed to the interfacial area difference: larger interfacial area will facilitate more low-temperature surfaces for the heat generation region, thereby lowering the maximal temperature in that region.

<p>Table 2</p>

Perimeters of four blade cross sections having the same area of π

Cross section shape

Perimeter

Triangle

= 8.080642

Rectangle

= 7.089816

Hexagon

= 6.597817

Circle

2π = 6.283186

Therefore, under the same composing materials ( = 641.6667), composition ratio (ϕ = 0.05), volumetric heat generation rate and total volume, the cylinder--triangular-prism and rectangular-prism--triangular-prism building blocks with slenderness values of 0.25 and 0.3, respectively, should be used for achieving the lowest system overall temperature difference (or, system thermal resistance) in practical applications. Both the pursuits of energy and material savings make this aim very significant. Furthermore, if one can also set the outer shape of the heat conduction building block free (often constrained by efficient packing and manufacturing), a comparison can be made between the constructal (both the blade shape and slenderness having been optimized) cylinder--triangular-prism building block and constructal rectangular-prism--triangular-prism building block. Since the total volumes for the cylinder--series and rectangular-prism--series building blocks are and , respectively, to ensure that the comparison is based on the same total building-block volume, the nondimensional constructal system overall temperature difference of the cylinder--triangular-prism building block is divided by π2/3

and the nondimensional constructal system overall temperature difference of the rectangular-prism--triangular-prism building block is divided by 22/3

Thus, the constructal cylinder--shaped heat conduction building block performs better than the rectangular-prism--shaped building block.

Conclusions

Inspired by the duct cross section optimization for minimum flow resistance, the shape of the nanoparticle blade is optimized for the cylinder--shaped and rectangular-prism--shaped building blocks of the blade-configured heat conduction systems (blade configuration of nanofluids) based on the same composing materials, composition ratio, volumetric heat generation rate, and total building block volume. The four kinds of blade shapes are triangular prism, rectangular prism, hexagonal prism, and cylinder. For the cylinder--cylinder building block, analytical analysis can be conducted. Explicit expressions for the system overall temperature difference, constructal system overall temperature difference, and constructal slenderness can be obtained. Then, based on the slenderness range predicted by the analytical result, numerical calculations are performed for the eight kinds of building blocks to obtain as accurately as possible results for comparison. One specifies that ϕ = 0.05 and = 641.6667 for the numerical calculation.

The performances of the eight kinds of building blocks depend strongly on the building-block slenderness. The constructal slendernesses leading to minimum system overall temperature differences (system thermal resistances) are 0.25 and 0.3, respectively, for the cylinder--series and rectangular-prism--series building blocks. For both the cylinder--series and rectangular-prism--series building blocks, the triangular-prism--shaped blade performs the best among all the four kinds of blades considered. This is explained by the size of interfacial area sustained by the four kinds of blades with a fixed volume. Also, the constructal cylinder--triangular-prism building block is proved to perform better than the constructal rectangular-prism--triangular-prism building block at the same composing materials, composition ratio, volumetric heat generation rate and total building-block volume.

Abbreviations

CFD: computational fluid dynamics.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

Both authors contributed equally.

Acknowledgements

The financial support from the Research Grants Council of Hong Kong (GRF718009 and GRF717508) and the CRCG of the University of Hong Kong (Project 10400920) is gratefully acknowledged.

 <p>Nanofluids: from vision to reality through research</p>ChoiSUSJ Heat Transfer200913103310610.1115/1.3056479<p>Nanofluids: synthesis, heat conduction, and extension</p>WangLQWeiXHJ Heat Transfer200913103310210.1115/1.3056597<p>The thermal conductivity of aqueous nanofluids containing ceria nanoparticles</p>BeckMPYuanYHWarrierPTejaASJ Appl Phys201010706610110.1063/1.3330506<p>Nanofluid convection in microtubes</p>LeeJGharagozlooPEKoladeBEatonJKGoodsonKEJ Heat Transfer201013209240110.1115/1.4001637<p>A comparative study of thermal behavior of iron and copper nanofluids</p>SinhaKKavlicogluBLiuYMGordaninejadFGraeveAJ Appl Phys200910606430710.1063/1.3225574<p>Effect of aggregation and interfacial thermal resistance on thermal conductivity of nanocomposites and colloidal nanofluids</p>EvansWPrasherRFishJMeakinPPhelanPKeblinskiPInt J Heat Mass Transfer2008511431143810.1016/j.ijheatmasstransfer.2007.10.017<p>Role of Brownian motion in the enhanced thermal conductivity of nanofluids</p>JangSPChoiSUSAppl Phys Lett2004844316431810.1063/1.1756684<p>Enhanced thermal conductivity of nanofluids: a state-of-the-art review</p>ÖzerinçSKakaçSYazıcıoğluAGMicrofluid Nanofluid20108145170<p>A benchmark study on the thermal conductivity of nanofluids</p>BuongiornoJVenerusDCPrabhatNMcKrellTTownsendJChristiansonRTolmachevYVKeblinskiPHuLWAlvaradoJLBangICBishnoiSWBonettiMBotzFCecereAChangYChenGChenHSChungSJChyuMKDasSKDi PaolaRDingYLDuboisFDzidoGEapenJEscherWFunfschillingDGalandQGaoJWGharagozlooPEGoodsonKEGutierrezJGHongHPHortonMHwangKSIorioCSJangSPJarzebskiABJiangYRJinLWKabelacSKamathAKedzierskiMAKiengLGKimCKimJHKimSLeeSHLeongKCMannaIMichelBNiRPatelHEPhilipJPoulikakosDReynaudCSavinoRSinghPKSongPXSundararajanTTimofeevaETritcakTTuranovANVan VaerenberghSWenDSWitharanaSYangCYehWHZhaoWHZhouSQJ Appl Phys200910609431210.1063/1.3245330<p>The limiting behavior of the thermal conductivity of nanoparticles and nanofluids</p>TejaASBeckMPYuanYHWarrierPJ Appl Phys201010711431910.1063/1.3354094BejanAShape and Structure, From Engineering to NatureCambridge: Cambridge University Press2000BejanALorenteSDesign with Constructal TheoryHoboken: Wiley2008<p>Unifying constructal theory of tree roots, canopies and forests</p>BejanALorenteSLeeJJ Theor Biol200825452954010.1016/j.jtbi.2008.06.02618647610<p>The constructal law and the thermodynamics of flow systems with configuration</p>BejanALorenteSInt J Heat Mass Transfer2004473203321410.1016/j.ijheatmasstransfer.2004.02.007<p>Constructal theory of generation of configuration in nature and engineering</p>BejanALorenteSJ Appl Phys200610004130110.1063/1.2221896<p>Constructal theory: from engineering to physics, and how flow systems develop shape and structure</p>ReisAHAppl Mech Rev20065926928210.1115/1.2204075<p>Constructal design for cooling a disc-shaped area by conduction</p>RochaLAOLorenteSBejanAInt J Heat Mass Transfer2002451643165210.1016/S0017-9310(01)00269-1<p>The constructal law of design and evolution in nature</p>BejanALorenteSPhilos Trans R Soc B20103651335134710.1098/rstb.2009.0302<p>"Percolation effect" and mechanical behavior of carbon black filled polyethylene</p>ChodakIKrupaIJ Mater Sci Lett1999181457145910.1023/A:1006665527806<p>Constructal allocation of nanoparticles in nanofluids</p>BaiCWangLQJ Heat Transfer201013205240410.1115/1.4000473<p>Constructal structure of nanofluids</p>BaiCWangLQJ Appl Phys201010807431710.1063/1.3486501<p>FLUENT 6.3 Documentation, Fluent Inc</p>