Journal of Nanotechnology

Volume 2018, Article ID 4349572, 9 pages

https://doi.org/10.1155/2018/4349572

## Numerical Research on Convective Heat Transfer and Resistance Characteristics of Turbulent Duct Flow Containing Nanorod-Based Nanofluids

^{1}School of Mechanical Engineering, Jiangnan University, Wuxi, China^{2}Jiangsu Key Laboratory of Green Process Equipment, Changzhou University, Changzhou, China^{3}School of Aeronautics and Astronautics, Zhejiang University, Hangzhou, China

Correspondence should be addressed to Jianzhong Lin; nc.ude.ujz.cilbup@nilzjcem

Received 25 October 2017; Accepted 21 January 2018; Published 6 March 2018

Academic Editor: Carlos R. Cabrera

Copyright © 2018 Fangyang Yuan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A coupled numerical model for nanorod-based suspension flow is constructed, and the convective heat transfer and resistance characteristics of the nanofluid duct flow are investigated. The numerical results are verified by experimental results and theoretical models. Most of nanorods are located randomly in the bulk fluid, while particles near the wall aligned with the flow direction. Friction factor of nanofluids with nanorods increases with higher particle volume concentration or aspect ratio, but the increment reduces when the Reynolds number gets larger. The relative Nusselt number is obtained to characterize the intensity of convective heat transfer. The results show that the Nusselt number of nanofluids increases when the particle volume concentration or aspect ratio becomes larger. Compared to increasing the aspect ratio of nanorods, increasing the particle volume concentration would be more effective on enhancing the convective heat transfer intensity in industrial applications although it will cause a slight increase of resistance.

#### 1. Introduction

Nanofluids, which are fluids with nanosized particles suspended, are considered to be next-generation mediums for advanced heat transfer or cooling technology [1]. Since Choi put forward the concept firstly [2], the study of nanofluids has been intensively performed and reported in worldwide. Abnormal heat transfer characteristics were reported by experiments including the effective heat conductivity, natural and forced convection heat transfer, and boiling heat transfer of nanofluids [3–5]. Factors that can influence thermophysical properties have been discussed in many literatures, such as nanoparticle size and shape, material of particle and base fluid, particle loading, temperature, and so on, but quantitative understanding is still a lack on the topic [1, 3]. Another hot topic is the flow resistance characteristic of nanofluids [6, 7]. The pumping power needed for the bulk fluid did not increase dramatically as the particle loading grows in experiments [5], whereas the viscosity of nanofluids shows a marked increase after adding nanoparticles into base fluid [1, 6]. The rheology of nanofluids also leads to the complexity on characterizing physical properties of bulk flow [8, 9].

Although many researchers have reviewed lots of measured experimental data and demonstrated that particle shape plays an important role on convective heat transfer characteristic of nanofluids [10, 11], few models can be found in literatures proposed for nanofluid flow containing nonspherical nanoparticles. Nanorods, as well as nanotubes, nanofibers, and nanobelts are used more and more for nanofluid preparation because of larger aspect ratio and surface area compared with spherical nanoparticles. Such kind of one-dimensional nanomaterials are available for heat transfer application with the development of nanoscale synthesis technology [12]. Xie et al. [13] first found the thermal conductivity of nanofluids with elongated nanoparticles is superior to that of spherical nanoparticles. Murshed et al. [14] measured enhanced thermal conductivity of water-based nanofluid containing TiO_{2} nanorods with the aspect ratio of 4; they showed that the results give 12% higher than that predicted by classical Hamilton–Crosser model. Then, two models derived considering interfacial layer existence between particle and base fluid for the effective thermal conductivities of nanofluids with spherical and cylindrical nanoparticles [15]. Yang and Han [16] investigated the thermal conductivity of nanofluids dispersed by Bi_{2}Te_{3} nanorod in perfluoro-*n*-hexane and oil. They found the thermal conductivities of nanofluids decrease with increasing temperature due to nanorod aspect ratio. The results measured by Zhang et al. [17, 18] employing carbon nanotubes (CNTs) agreed well with the predictions of the unit-cell model proposed by Yamada and Ota [19] by considering shape factor for cylindrical particles. By using the excluded volume concept, Koo et al. [20] constructed a model which can reproduce the nonlinear increase of the thermal conductivity for carbon nanotube and nanofiber suspensions. Timofeeva et al. [21] studied the particle shape effect on thermophysical properties of alumina-EG/H_{2}O suspensions and concluded that elongated particles and agglomerates resulted in higher viscosity at the same volume fraction due to structural limitation of rotational enhancements in effective thermal conductivities diminished by interfacial effects. Molecular dynamics (MD) simulation based on the model developed by Ghosh et al. [22] have shown that cylindrical nanoparticles pick up thermal energy much faster than spherical nanoparticles during a collision with a block-shaped heat source [23]. Comparative experimental results reported by Yu et al. [10] indicated that convective heat transfer characteristic of nanofluids containing nonspherical nanoparticles is affected significantly by the shear-induced alignment and orientational motion of particles. By considering the thickness and thermal conductivity of the interfacial nanolayer, Jiang et al. [24] proposed a model which provides good predictions for the effective thermal conductivity of CNT-based nanofluids. Other existing models are summarized by Yang et al. [11], and three similar theoretical models are proposed for nanofluids with finite cylindrical particles by anisotropy analysis [25–27]. All of the three models give minimal relative errors to experimental results compared to other models. Lin et al. [7, 28] took numerical simulation on the phaolefins-based nanofluids fluid flow containing cylindrical Al_{2}O_{3} nanoparticles in laminar and turbulent pipes, respectively. The results show that the friction factor of flow decreases when the Reynolds number and particle aspect ratio become larger. The Nusselt number of convective heat transfer increases with the increase of Reynolds number, particle aspect ratio, and volume concentration. By comparing the defined performance evaluation criterion (PEC), nanofluids containing rod-like nanoparticles with large aspect ratio and a suitable particle volume concentration are more effective for convective heat transfer process at higher Reynolds number. A coupled numerical model is constructed by Yuan et al. [29] to simulate convective heat transfer and resistance characteristics of TiO_{2}/water nanofluids with cylindrical particles in laminar channel flow. The results show that the Nusselt number of nanofluid flow related to the Reynolds number, axial length, Prandtl number, and particle volume concentration and a fitted formula is proposed to predict the Nusselt number of nanofluid flow containing cylindrical nanoparticles.

As summarized above, high potential of rod-like nanoparticles in enhancing thermal conductivity of nanofluids has got more attention by researchers. Related reports focus mostly on viscosity and thermal conductivity properties, while progress made on convective heat transfer characteristics is still highly needed. In addition, lack of numerical modelling on nanofluid flow brings an obstacle to industrial applications. The objective of this work is to develop a coupled model for nanorod-based suspension flow and investigate the convective heat transfer and resistance characteristics of the nanofluid duct flow numerically. The orientation and movement of nanorod are considered, and key factors that can influence the friction factor and Nusselt number of the fluid flow are discussed in sequence.

#### 2. Model for Nanorod-Based Nanofluids

##### 2.1. Governing Equations

The duct flow for nanofluid is considered as impressible, and the governing equations include the continuity equation, modified Navier–Stokes equation with the additional term of nanorods, and the energy conservation equation [28, 30, 31]:where , *p*, *ρ*_{nf}, and *T* are the nanofluid velocity, pressure, density, and temperature, respectively; *μ* is the dynamic viscosity of the pure fluid; *μ*_{a} is the apparent viscosity; and are the mean second- and fourth-order tensors of particle orientation, respectively; is the rate-of-strain tensor; and *C*_{nf} is the thermal diffusivity coefficient of the nanofluid, (*k* is the turbulent kinetic energy, *ε* is the turbulent dissipation rate, *C*_{μ} = 0.09, and turbulent Prandtl number, Pr_{T} = 0.9) is the eddy thermal diffusivity coefficient.

Substituting instantaneous velocity, pressure, temperature, rate-of-strain tensor, and tensors of particle orientation which consist of mean and fluctuation part in (1–3) and averaging, we havewhere *U*_{i}, *P*, and are the mean nanofluid velocity, pressure, and temperature, respectively, is the fluctuation flow velocity, and is the mean rate-of-strain tensor.

In (5), the nanofluid density iswhere the subscripts “f” and “p” stand for pure fluid and nanoparticles, respectively, and *ϕ*_{e} is the effective particle volume concentration, which is used for replacing the nominal particle volume fraction *ϕ*, in consideration of the particle aggregation:where *s* and *s*_{e} are the effective size of aggregates and primary particle size, respectively, and *f*_{i} is the fractal index and varies from 1.5 to 2.45 for those of rod-like nanoparticles. Yu et al. [10] gives the value as: *s*/*s*_{e} = 1.48 and *f*_{i} = 1.95 according to experimental data.

Batchelor’s theory was extended by Mackaplow and Shaqfeh [32] to account for two-body interactions and gave the apparent viscosity in (5):where *n* and *l* are the number density and half length of nanorods, respectively, *μ* is the dynamic viscosity of the pure fluid, and *r* is the particle aspect ratio, which is the ratio of length of diameter.

The thermal diffusivity coefficient in (6) is [33, 34]where the heat capacitance of the nanofluid is given by

In (10), the classical Hamilton–Crosser model is employed for predicting thermal conductivity of nanorod-based nanofluid as follows:in which *k*_{f} and *k*_{p} are the thermal conductivity of pure fluid and nanoparticles, respectively, and *K* is the shape factor given by *K* = 3/*γ*, where *γ* is the particle sphericity, defined as the ratio of the surface area of a sphere with the same volume as that of the particle and the surface area of the particle. The model predicts well the thermal conductivity of nanofluid with ZnO nanoparticles with experimental data by Ferrouillat et al. [35] in rod-like shape.

The Reynolds stress tensor in (5) iswhere eddy viscosity . For solving (5), the -equation and -equation for turbulent flow are given bywhere *C*_{1} = 1.44, *C*_{2} = 1.92, *σ*_{k} = 1.0, and *σ*_{ε} = 1.3.

##### 2.2. Nanoparticle Orientation Distribution

The mean second- and fourth-order tensors of particle orientation in (5) describe the orientation distribution of nanoparticle population in flow field; they can be calculated by Advani and Tucker [36] aswhere *p*_{i} is a unit vector parallel to the particle’s axis and is the probability density function for particle orientation at any position with **p** being the orientation vector. The governing equation of is the modified Fokker–Planck equation considering the rotary effect caused by Brownian movement [37, 38]:where *D*_{rB} is the Brownian rotary diffusion coefficient, is the gradient operator projected onto the surface of the unit sphere, and is the particle angular velocity. Cintra and Tucker [39] expressed aswhere is the vorticity tensor, . The last term in (18) is introduced to model the behavior at higher concentrations, in which *D*_{rI} is a rotary diffusion coefficient resulting from particle-particle interactions. Folgar and Tucker [40] suggested that *D*_{rI} is isotropic and given by , where the interaction coefficient *C*_{I} is fixed to 0.01.

For cylindrical particles, the Brownian rotary diffusion coefficient in (17) is

With aspect ratio *r*, the rotational friction coefficients around long axis and short axis are [41, 42].with

Substituting instantaneous and which consist of mean and fluctuation part in (17) and averaging, we have the mean equation of probability density functions for the particle orientation:where and are the dispersion coefficients of linear and angular displacement in which *ν* is the kinetic viscosity of the pure fluid [43], which appear based on the dispersion mechanism:

Integral on the angle after multiplying component of unit vector *p*_{i}, (22) is transformed to the equation about :

The mean fourth-order particle orientation tensor needs to be closured by model. Considering the shearing-stretch feature of duct flow, the orthotropic closure model [44] is fit for (24):in which det(*a*_{ij}) is the determinant of second-order tensor **a**_{2}.

##### 2.3. General Dynamics Equation for Rod-Like Nanoparticles

The particle volume concentration *ϕ* in (7), (11), and (12) is considered distribution inhomogeneity in the flow. In order to obtain the distribution of *ϕ* in the flow, it is necessary to solve the general dynamics equation for nanorods:

Based on the Reynolds average, the above equation is transformed towhere is the mean particle volume distribution function based on the volume and the last term on the left-hand side represents the change in resulting from turbulent diffusion and is usually assumed to be [45]:where *ε*_{i} is the eddy diffusion coefficient which is a function of position and , where *v*_{t} and Sc_{t} are the turbulent viscosity and the Schmidt number, respectively.

*D*_{tB} in (26) is the Brownian translational diffusion coefficient for rod-like nanoparticles:and the translational friction coefficients parallel and perpendicular to the long axis are [41, 42]with

Substituting (29) in (27), we have averaged general dynamics equation for nanoparticles. With the moment transformation which involves multiplying the equation by *v*^{k} and then integrating over the entire size regime, the equation for rod-like nanoparticles finally becomeswhere the zero moment and the first-order momentdescribe the total particle number density and volume of nanoparticles at a given point.

#### 3. Numerical Method

Schematic of a turbulent duct flow containing nanorod-based nanofluid is shown in Figure 1. The length and diameter of the channel are *L* and *D*, respectively.