# CFD Study on Wall/Nanoparticle Interaction in Nanofluids Convective Heat Transfer

Authors: Mohammad Reza Tarybakhsh; Ali Akbar Lotfi Neyestanak; Hamed Tarybakhsh
Journal: Advances in Materials Science and Engineering
 (2013)

---

## Abstract

The Brownian motion of the nanoparticles in nanofluid is one of the potential contributors to enhance effective thermal conductivity and the mechanisms that might contribute to this enhancement are the subject of considerable discussion and debate. In this paper, the mixing effect of the base fluid in the immediate vicinity of the nanoparticles caused by the Brownian motion was analyzed, modeled, and compared with existing experimental data available in the literature. CFD was developed to study the effect of wall/nanoparticle interaction on forced convective heat transfer in a tube under constant wall temperature condition. The results showed that the motion of the particle near the wall which can decrease boundary layer and the hydrodynamics effects associated with the Brownian motion have a significant effect on the convection heat transfer of nanofluid.


---

## Body

## 1. Introduction

Fluid heating and cooling play significant roles in a lot of industrial processes such as refinery, petrochemical, power stations, and electronics [1]. Heat transfer in traditional fluids such as water, ethyleneglycol, and oil inherently has low thermal conductivity compared with that in the metals and metal oxides. Therefore, fluids with suspended solid particles are expected to have better heat transfer properties. But, for the particles with millimeter and micron size, there are some problems like sedimentation, clogging, and pressure drop. Choi employed the particles in nanometer dimensions as a suspended solution [2]. An increasing number of experimental investigations have demonstrated that nanofluids can dramatically increase the effective thermal properties of the base fluid [3–12]. Abu-Nada [13] investigated the effects of variable viscosity and thermal conductivity of a nanofluid (Al2O3-water) on the natural convective heat transfer. Sharma et al. [14] experimentally studied the convective heat transfer coefficient and pressure drop in the transient region for Al2O3-water nanofluid under a constant heat flux. They found that convective heat transfer increases by adding Al2O3 nanoparticles in water. Moraveji et al. [15] simulated water-Al2O3 nanofluid through a tube under a constant heat flux. They found that the heat transfer coefficient rises by increasing the nanoparticle concentration and Reynolds number. Furthermore, the heat transfer coefficient increases by particle diameter reduction.Because of these suitable properties and the excellent stability of these fluids, nanofluids present a promising alternative to traditional heat transfer fluids in a wide range of applications. However, the precise mechanisms that contribute to the observed enhancement are not currently well understood and no widely accepted explanations have been identified. As a result, in order to fully understand the mechanisms that govern the enhancement of these nanofluids and optimize the thermophysical properties, the mixing effect of the base fluid in the immediate vicinity of the nanoparticles caused by the Brownian motion of the nanoparticles was analyzed, modeled, and compared with existing experimental data available in the literature. Due to the lack of experimental tools required to isolate the effect of each mechanisms, many investigations have focused only on the Brownian motion effects of the nanoparticles, using either numerical simulation or experimental techniques [16–22]. Considering some of the molecular simulations, it was concluded that Brownian motion of nanoparticles has no impact on the effective thermal conductivity of the nanofluids [18], while the experimental study of nanofluids diffusion provides some evidence for the contribution of the Brownian motion effects, along with other factors [18–20]. In order to have better understanding of the governing phenomena and to determine the contribution of Brownian motion, a numerical simulation was conducted to determine the validity of the previously stated conclusions [23, 24].In this research, the convective heat transfer near the wall region of the tube flow containing water and Al2O3 nanofluid under a constant temperature was simulated using the Computational Fluid Dynamics (CFD) tools. Al2O3 nanoparticle with average diameters of 27 nm was used. The effects of the nanoparticle Brownian motion on the convective heat transfer coefficient were investigated near the wall.


## 2. The Simulation of Brownian Motion Effect on the Effective Thermal Convection

Visual observations of Brownian motion indicate that each nanoparticle can be modeled as having a local periodic motion within the suspension, as shown in Figure1. As illustrated, points A and C can be used to represent the farthest points of local periodic motion, and point B is the location at which the local periodic motion has the highest velocity. Based upon this diagram, the velocity and the range of influence can be determined from the following expressions:
(1)md2ydt2+by=0,x(t)=y0cosωt+v0ωsinωt,
where m is mass, x displacement, t time, b constant, x angular velocity, m0 viscosity, and x0 the initial location. UFD code in the model includes Brownian motion.Figure 1
Trajectory of nanoparticle.


### 2.1. CFD Model

The influence range for this local periodic motion model can be calculated by solving the governing equations of the convection caused by the motion of Al2O3 nanoparticles in three dimensions. Figure 2 shows the geometry of CFD model. The Brownian motion of nanoparticle in nanofluid flow with a velocity of 1 m/s and a temperature of 24°C in a tube was modeled by CFD method and the result was discussed. Water flow in a small tube (diameter = 0.2 micron, length = 1 micro) with nanoparticle motion near the wall was studied. The particle has a constant velocity of 1 m/s in x-direction and a periodic motion in y-direction.Figure 2
Geometry and boundary condition of system.The surface of nanoparticles was assumed to be adiabatic, and wall 2 was supposed to be a heated surface, with a constant temperature of 24.1°C, as shown in Figure2. For the purpose of parametric study, initial temperature of water flow, Tbulk,water, was set at 24°C with a velocity of 1 m/s. To study the convective heat transfer enhancement effect, temperature gradient was assumed between fluid and wall temperatures; in addition, particle motion effects were modeled near the wall.To model fluids and heat transfer, the two-dimensional computational grids consisting unstructured cells were used separately in some cases. It has been proven that the unstructured cell technology is a significant improvement in terms of meshing flexibility and simulation time to perform a complete simulation. The mesh size surrounding the particle is too fine to predict the variation of pressure, velocity, and temperature with high accuracy, but it has normal size near the walls to increase the calculation speed (Figure3).Figure 3
Mesh model of base case.In Figure3 the mesh is generated with 1.25 million mesh elements. The nanoparticle has a diameter of 26 nm including the possible adsorption layer [25] which is the actual regime for mutual interaction between liquid molecules and nanoparticle surface. The flow domain volume is 200 × 1000 nm. For this model, the transient state was modeled. The actual time needed to compute is approximately 24 h in order to achieve a convergent steady-state solution. A refined grid analysis may have resulted in a greater mesh refinement in the region of nanoparticles, but the current grid appeared to have enough accuracy for this single nanoparticle case study.Transient CFD code was prepared based on the SIMPLE algorithm and second order upwind method allows calculating with Cartesian and cylindrical coordinates. Modeling of the Brownian motion via CFD requires geometry specification through identifying the computational grid, numerical solution strategy, and specification conditions.The base fluid phase considered here was assumed to be Newtonian with laminar flow and to have constant physical properties except heat capacity and conductivity coefficient. Nanoparticle periodic motion model and the governing equations (continuity, momentum, and energy) for the convection caused by nanoparticles motion in two-dimensional form can be determined as shown below:(2)∂ρ∂τ+∂ρu∂x+∂ρυ∂y=0,ρ(∂u∂τ+u∂u∂x+υ∂u∂y)=Fx-∂p∂x+μ(∂2u∂x2+∂2u∂y2),ρ(∂u∂τ+u∂u∂x+υ∂u∂y)=Fx-∂p∂x+μ(∂2u∂x2+∂2u∂y2),ρCpdTdτ=∂∂x(k∂T∂x)+∂∂y(k∂T∂y)+αυTdPdτ+μΦ,
where ρ is the density, u, υ, w, Fx, and Fy are the velocity components and external forces in the x-, y-directions, respectively, P is the static pressure, Cp is the specific heat capacity, T is the temperature, k is the thermal conductivity, αυ is the expansion coefficient, μ is the viscosity, and Φ is the diffusion matrix. The grid model was generated by Gambit and solved by Ansys Fluent 6.3 and Brownian motion was added to the fluent model by UFD file. The transient state was modeled using a time step of dt=1e-7 s and a total time of 0.001 s. To study the Brownian motion effect, initially the water fluid flow in a tube containing no nanoparticle and that in the case of single nanoparticle were studied besides the isothermal, pressure, and velocity fields. Finally, the effect of multiple nanoparticles was studied, too.


## 2.1. CFD Model

The influence range for this local periodic motion model can be calculated by solving the governing equations of the convection caused by the motion of Al2O3 nanoparticles in three dimensions. Figure 2 shows the geometry of CFD model. The Brownian motion of nanoparticle in nanofluid flow with a velocity of 1 m/s and a temperature of 24°C in a tube was modeled by CFD method and the result was discussed. Water flow in a small tube (diameter = 0.2 micron, length = 1 micro) with nanoparticle motion near the wall was studied. The particle has a constant velocity of 1 m/s in x-direction and a periodic motion in y-direction.Figure 2
Geometry and boundary condition of system.The surface of nanoparticles was assumed to be adiabatic, and wall 2 was supposed to be a heated surface, with a constant temperature of 24.1°C, as shown in Figure2. For the purpose of parametric study, initial temperature of water flow, Tbulk,water, was set at 24°C with a velocity of 1 m/s. To study the convective heat transfer enhancement effect, temperature gradient was assumed between fluid and wall temperatures; in addition, particle motion effects were modeled near the wall.To model fluids and heat transfer, the two-dimensional computational grids consisting unstructured cells were used separately in some cases. It has been proven that the unstructured cell technology is a significant improvement in terms of meshing flexibility and simulation time to perform a complete simulation. The mesh size surrounding the particle is too fine to predict the variation of pressure, velocity, and temperature with high accuracy, but it has normal size near the walls to increase the calculation speed (Figure3).Figure 3
Mesh model of base case.In Figure3 the mesh is generated with 1.25 million mesh elements. The nanoparticle has a diameter of 26 nm including the possible adsorption layer [25] which is the actual regime for mutual interaction between liquid molecules and nanoparticle surface. The flow domain volume is 200 × 1000 nm. For this model, the transient state was modeled. The actual time needed to compute is approximately 24 h in order to achieve a convergent steady-state solution. A refined grid analysis may have resulted in a greater mesh refinement in the region of nanoparticles, but the current grid appeared to have enough accuracy for this single nanoparticle case study.Transient CFD code was prepared based on the SIMPLE algorithm and second order upwind method allows calculating with Cartesian and cylindrical coordinates. Modeling of the Brownian motion via CFD requires geometry specification through identifying the computational grid, numerical solution strategy, and specification conditions.The base fluid phase considered here was assumed to be Newtonian with laminar flow and to have constant physical properties except heat capacity and conductivity coefficient. Nanoparticle periodic motion model and the governing equations (continuity, momentum, and energy) for the convection caused by nanoparticles motion in two-dimensional form can be determined as shown below:(2)∂ρ∂τ+∂ρu∂x+∂ρυ∂y=0,ρ(∂u∂τ+u∂u∂x+υ∂u∂y)=Fx-∂p∂x+μ(∂2u∂x2+∂2u∂y2),ρ(∂u∂τ+u∂u∂x+υ∂u∂y)=Fx-∂p∂x+μ(∂2u∂x2+∂2u∂y2),ρCpdTdτ=∂∂x(k∂T∂x)+∂∂y(k∂T∂y)+αυTdPdτ+μΦ,
where ρ is the density, u, υ, w, Fx, and Fy are the velocity components and external forces in the x-, y-directions, respectively, P is the static pressure, Cp is the specific heat capacity, T is the temperature, k is the thermal conductivity, αυ is the expansion coefficient, μ is the viscosity, and Φ is the diffusion matrix. The grid model was generated by Gambit and solved by Ansys Fluent 6.3 and Brownian motion was added to the fluent model by UFD file. The transient state was modeled using a time step of dt=1e-7 s and a total time of 0.001 s. To study the Brownian motion effect, initially the water fluid flow in a tube containing no nanoparticle and that in the case of single nanoparticle were studied besides the isothermal, pressure, and velocity fields. Finally, the effect of multiple nanoparticles was studied, too.


## 3. Result and Discussion

### 3.1. Case 1: Laminar Fluid Flow without Nanoparticle

In the first case, To show the effect of nanoparticle interaction with wall in fluid flow, first the model was performing without nanoparticle and then the nanopartical interaction effect was studied. The boundary condition was defined as mentioned above. Figure4 illustrates the temperature gradient expected in the flow tube. According to the temperature counter the parameters gradient is very high in the entrance zone and not fully developed.Figure 4
Case 1, temperature gradient (time = 0.001 S).


### 3.2. Case 2: Laminar Fluid Flow with One Nanoparticle

In this case, a nanoparticle was located near wall 2 at the same distance from the inlet and outlet. The particle moves through fluid inx-direction (V=1 m/s) and has very small periodic motion in y-direction. The boundary and initial conditions were temperature (walls 1 and 2 :  24.1°C, inlet nanofluid temp.: 24°C), velocity (inlet nanofluid: 1 m/s, initial nanoparticle: 001 m/s), and dt=1e-7. Figure 5 shows the typical simulated velocity, temperature, and pressure field around a nanoparticle after the first stage of modeling. As indicated, mixing flow has the enhancement effect of Brownian motion in nanofluid on convective heat transfer in comparison with the results of Case 1. As the nanoparticle surface was treated as adiabatic, there is an isothermal curvature because of the microconvection induced by Brownian motion of nanoparticle.Case 2, velocity, temperature and pressure gradient, and mesh model (time = 0.001 S).
(a)
(b)
(c)
(d)The pressure field at the top of Figure5 shows that the pressure in front of the nanoparticle increases very quickly and it is larger than that behind the nanoparticle.Figure6 illustrates the heat transfer coefficient in walls 1 and 2 versus tube length. Figure 6(a) shows that heat transfer coefficient changes through the whole length of the tube and Figure 6(b) indicates changes in the range of nanoparticles.Case 2, convective heat transfer coefficient on walls 1 and 2.
(a)
(b)Increase in convective heat transfer coefficient near the wall 2 shows the effect of nanoparticle Brownian motion. This motion causes reduction in heat and fluid boundary layer near the wall. Furthermore, the microconvection of vibration motion leads to increasing the convective heat transfer coefficient, as obviously seen in the figure.


### 3.3. Case 3: Laminar Fluid Flow with Two Adjacent Nanoparticles and the Effect of Brownian Motion

In this case, the effect of two adjacent nanoparticles, instead of one nanoparticle was studied. The distance between two particles is 30 nm. The boundary and initial conditions are the same as Case 1.Figure7 illustrates the velocity, temperature, and pressure gradient.Case 3, velocity, temperature, pressure gradient, and mesh model (time = 0.001 S).
(a)
(b)
(c)
(d)Compared with the results and heat transfer rate related to the model of Figures5 and 7, it is clear that the second nanoparticle causes increase in heat transfer from the wall to the fluid which is shown by the counters of these two figures. This improvement is due to the mutual effects of the nanoparticles. Therefore, increase in nanoparticle concentration because of microconvection enhancement leads to higher heat transfer coefficient in the fluid.The model of two adjacent nanoparticles clarifies how the induced microconvection influences the convective heat transfer coefficient of base fluid. Also, the modeling of two adjacent nanoparticles clarifies how the induced microconvection effects the heat transfer capability.Figure8 illustrates the heat transfer coefficient of walls 1 and 2. As seen in the figure, blue parts which indicate convective heat transfer coefficient near wall 2 show larger growth rate compared with that in Figure 6 in the case of one nanoparticle. So, it can imply the effect of more particles in the base fluid that researchers have mentioned [23, 24].Case 3, convective heat transfer coefficient on walls 1 and 2.
(a)
(b)When two nanoparticles are close to each other, the influence on the area will be more than twice due to the hydrodynamic interaction. This, in turn, increases the heat transfer capacity of the nanofluid at macroscale.


## 3.1. Case 1: Laminar Fluid Flow without Nanoparticle

In the first case, To show the effect of nanoparticle interaction with wall in fluid flow, first the model was performing without nanoparticle and then the nanopartical interaction effect was studied. The boundary condition was defined as mentioned above. Figure4 illustrates the temperature gradient expected in the flow tube. According to the temperature counter the parameters gradient is very high in the entrance zone and not fully developed.Figure 4
Case 1, temperature gradient (time = 0.001 S).


## 3.2. Case 2: Laminar Fluid Flow with One Nanoparticle

In this case, a nanoparticle was located near wall 2 at the same distance from the inlet and outlet. The particle moves through fluid inx-direction (V=1 m/s) and has very small periodic motion in y-direction. The boundary and initial conditions were temperature (walls 1 and 2 :  24.1°C, inlet nanofluid temp.: 24°C), velocity (inlet nanofluid: 1 m/s, initial nanoparticle: 001 m/s), and dt=1e-7. Figure 5 shows the typical simulated velocity, temperature, and pressure field around a nanoparticle after the first stage of modeling. As indicated, mixing flow has the enhancement effect of Brownian motion in nanofluid on convective heat transfer in comparison with the results of Case 1. As the nanoparticle surface was treated as adiabatic, there is an isothermal curvature because of the microconvection induced by Brownian motion of nanoparticle.Case 2, velocity, temperature and pressure gradient, and mesh model (time = 0.001 S).
(a)
(b)
(c)
(d)The pressure field at the top of Figure5 shows that the pressure in front of the nanoparticle increases very quickly and it is larger than that behind the nanoparticle.Figure6 illustrates the heat transfer coefficient in walls 1 and 2 versus tube length. Figure 6(a) shows that heat transfer coefficient changes through the whole length of the tube and Figure 6(b) indicates changes in the range of nanoparticles.Case 2, convective heat transfer coefficient on walls 1 and 2.
(a)
(b)Increase in convective heat transfer coefficient near the wall 2 shows the effect of nanoparticle Brownian motion. This motion causes reduction in heat and fluid boundary layer near the wall. Furthermore, the microconvection of vibration motion leads to increasing the convective heat transfer coefficient, as obviously seen in the figure.


## 3.3. Case 3: Laminar Fluid Flow with Two Adjacent Nanoparticles and the Effect of Brownian Motion

In this case, the effect of two adjacent nanoparticles, instead of one nanoparticle was studied. The distance between two particles is 30 nm. The boundary and initial conditions are the same as Case 1.Figure7 illustrates the velocity, temperature, and pressure gradient.Case 3, velocity, temperature, pressure gradient, and mesh model (time = 0.001 S).
(a)
(b)
(c)
(d)Compared with the results and heat transfer rate related to the model of Figures5 and 7, it is clear that the second nanoparticle causes increase in heat transfer from the wall to the fluid which is shown by the counters of these two figures. This improvement is due to the mutual effects of the nanoparticles. Therefore, increase in nanoparticle concentration because of microconvection enhancement leads to higher heat transfer coefficient in the fluid.The model of two adjacent nanoparticles clarifies how the induced microconvection influences the convective heat transfer coefficient of base fluid. Also, the modeling of two adjacent nanoparticles clarifies how the induced microconvection effects the heat transfer capability.Figure8 illustrates the heat transfer coefficient of walls 1 and 2. As seen in the figure, blue parts which indicate convective heat transfer coefficient near wall 2 show larger growth rate compared with that in Figure 6 in the case of one nanoparticle. So, it can imply the effect of more particles in the base fluid that researchers have mentioned [23, 24].Case 3, convective heat transfer coefficient on walls 1 and 2.
(a)
(b)When two nanoparticles are close to each other, the influence on the area will be more than twice due to the hydrodynamic interaction. This, in turn, increases the heat transfer capacity of the nanofluid at macroscale.


## 4. Conclusions

This paper discussed the Brownian motion effect as one of the main factors for enhancing forced convective heat transfer coefficient of nanofluids. The corresponding temperature, pressure, and velocity fields were simulated using CFD model and a finite-difference algorithm. The simulations for single and adjacent nanoparticles were discussed in detail. The results clearly indicated that microconvection/mixing induced by the Brownian motion of nanoparticles could significantly affect the macroconvective heat transfer capability of the nanofluids. This information is especially interesting when accompanied with other works on the variation of the viscosity of nanofluids, due to the Brownian motion and hydrodynamic interaction between nanoparticles [26–30]. The results indicated the effect of the microconvection in these types of suspensions and demonstrate that Brownian motion is one of the key factors in relation to the observed high effective thermal convection of nanofluids. Also, results showed that the temperature of two adjacent nanoparticles is higher than that of one nanoparticle, which means that the hydrodynamic interaction between the nanoparticles dramatically enhances the thermal transport capability and the nanoparticles, not only in front of, but also behind, the nanoparticles. It is also apparent that two adjacent nanoparticles in close proximity will have a greater influence on the temperature field than two single nanoparticles, far apart from each other, due to the interaction.



---

*Source: 651365-2013-11-07.xml*