Simulation of Natural Convection Heat Transfer Enhancement by Nanoparticles in an Open-Enclosure Using Lattice Boltzmann Method

A numerical analysis is performed to investigate the laminar, free convection flow in an Open Enclosure Using Lattice Boltzmann Method (LBM) in the presence of Carbon nanotube and Cu nanoparticles. The problem is studied for different volume fractions of nanoparticles, and aspect ratio of the cavity for various Rayligh numbers. The volume fraction of added nanoparticles to water (as base fluid) is lower than 1% to make dilute suspensions. The study presents a numerical treatment based on LBM to model convection heat transfer of Carbon nanotube based nanofluids. Results show that adding a low value of Carbon nanotube to the base fluid led to significant enhancement of convection rate. Results show that adding nanoparticles to the base fluid enhances the rate of natural convection in a cavity. Make a comparison between Carbon nanotube and Cu-nanoparticles shows that the Carbon nanotube-nanoparticle has better performance to enhance convection rate at comparison with Cunanoparticles.


Introduction
The enhancement of fluids heat transfer is an interesting topic for different kinds of industrial and engineering applications.A well-known way to enhance the rate of convection heat transfer of conventional fluids such as water, oil and ethylene glycol which has low thermal conductivity [1] [2] is adding nano scale conductive particles.The added particles can be metals [3], non metals [4] or Carbon nanotubes [5] [6].Because of the high thermal conductivity of these particles, they can improve conductivity of the suspensions systematically.Nowadays, nanoparticle added fluids are known as nanofluid that first time Choi [7] named this kind of fluid suspensions.He presented enhancement of convection heat transfer by adding nanoparticles to the fluids.In recent years, many studies have been conducted to study heat transfer of nanofluids numerically and experimentally [7]- [12].Khanafer et al. [8] presented a heat transfer enhancement by adding nanoparticles to fluid in a two dimensional enclosures at natural convection regime for the different Grashof numbers.They presented an increase of ave Nu about 7.5%, 12%, 15.5% and 20%, respectively, by adding a volume fraction of added Cu nanoparticles to the Water equal to 4%, 6%, 8% and 10% when 4 10 Gr = .
Gherasim et al. [9] conducted to an experimental study and represented heat transfer enhancement possibility of coolants with suspended Al 2 O 3 nanoparticles dispersed in water inside a radial flow cooling device.Saleh et al. [10]

Gr =
and also for Cu nanoparticles at same condition the values are about 3.9%, 7.9%, 11.8%, 15.7%, 19.5%.In 1991, Iijima [13] discovered Carbon nanotubes as an allotrope of Carbon which is made of long-chained molecules of Carbon with Carbon atoms arranged in a hexagonal complex to form a tubular structure.In last decade Carbon nanotubes has been mentioned as attractive topic by many researches which is generally due to their special properties at physical view such mechanical and thermal properties.At this point of view, the Carbon nanotubes have extraordinary thermal properties such as thermal conductivity about twice as high as diamond [14] or thermal stability up to 2800 C in vacuum.The higher thermal conductivity of Carbon nanotubes relative to other nanoparticles led to that the nanofluids containing cylindrical Carbon nanotubes are expected to have better heat transfer properties compared with the other nanofluids with spherical nanoparticles [15] [16].Gavili et al. [17] simulated the mixed convection in a two-sided lid-driven differentially heated square cavity for nanofluids containing Carbon nanotubes for nanofluids.Natural convection in a square cavity and its fluid flow is a classical benchmark in heat transfer problems.Open cavities are a kind of 2-D cavity which has an open side.These kinds of cavities have special physics for both flow and temperature fields in open side because of outgoing of flow from this side.Many studies have been done on analysis of buoyant flows and their heat transfer in open cavities.Javam et al. [18] investigated stability of stratified natural convection flow in open cavities.Mohammad et al. [19] presented natural convection in an open ended cavity and slots they analyzed the effect of aspect ratio of cavity on heat transfer rate.They also presented a good procedure for simulating open boundaries in Lattice Boltzmann Method, their technique demonstrated the abilities of the LBM to simulate Natural convection in open cavities.The progress of using the Lattice Boltzmann Method (LBM) as a numerical technique to simulate the heat transfer and fluid flow has been obvious in the last decade [20]- [23].The Lattice Boltzmann Method has well-known advantages such easy implementation, possibility of parallel coding and simulating of complex geometries and fluid dynamic problems such as melting, fuel cell, porous media, nanofluids and etc.The convection heat transfer of different nanoparticle-based nanofluid in a cavity was mentioned in current years [24]- [27].Fattahi et al. [24] applied Lattice Boltzmann Method to study natural convection heat transfer of Al2O3 and Cu nanofluid in a cavity.Their results show that increasing the solid volume fraction led to a heat transfer enhancement at any Rayleigh number and also heat transfer increases with increase in Rayleigh number for a particular volume fraction.Nemati et al.  [29] is used.On other hand, the study applies the Patel model [30] to evaluate of the effective thermal conductivity of Cu-nanofluid.To the best of author knowledge, the effects of carbon nanotubes on flow and thermal fields of natural convection is an unknown perspective which is not understand heretofore.Therefore, in this study the used numerical method is LBM with coupled double population approach for flow and temperature fields.The effect of the volume fraction of the Cu and Carbon nanotube-particles on average Nusselt number, streamlines and temperature contours is investigated for various Rayligh numbers.Also the aspect ratio of the cavity is studied as an important geometric parameter of the 2D cavities in different condition to exhibit the role of this parameter on both heat transfer rate and fluid flow of the base fluid and the nanofluid.

Lattice Boltzmann Method for Flow and Thermal Fields
The LBM utilizes two distribution functions, for the flow and temperature fields.It uses modeling of movement of fluid particles to define macroscopic parameters of fluid flow.The distribution functions are obtained by solving the lattice Boltzmann Equation (LBE), which is a special form of the Kinetic Boltzmann Equation.
The basic form of the Lattice Boltzmann Equation with an external force by introducing BGK approximations can be written as follows for the both flow and the temperature fields [20]: , , , , eq t t g x e t t t g x t g x t g x t m τ and t τ are the dimensionless collision-relaxation times for the flow and tem- perature fields, respectively.They defined as follows [31]: and temperature field ( eq g α ) are calculated as follows in different α directions: ( ) ( ) For these model the values of e α and w α for various α directions will be, re- spectively: In order to incorporate external force in collision part of Lattice Boltzmann model (Equation ( 1)), radiation heat transfer and viscous dissipation are neglected at the numerical simulation.Therefore, to capture buoyancy force effects in the flow field, the Boussinesq approximation is applied.Thus, to model buoyancy force in Equation ( 1), the external force needs to be assumed as below in the needed direction: ( ) Finally, macroscopic variables can be calculated as follows: Flow density:

Nanofluids Modeling
At present study, nanofluid is assumed as a single phase fluid.Thermal diffusivity of nanofluid is as follows: ( ) The density, heat capacitance and thermal expansion of nanofluid can be defined as [32]: ( ) The effects of nanoparticles on the viscosity of nanofluids are introduced by the socalled apparent viscosity which is presented by app µ [28]:

Main Problem
In the present study, the effects of Carbon nanotube and Cu-nanoparticles on natural convection heat transfer in a 2D deep open cavity are investigated using the Lattice Boltzmann Method based on double population approach and bounce back method.
The horizontal walls of the cavity are assumed to be insulated while the left wall is maintained at a uniform temperature ( h T ) differentially higher than the open end side temperature ( c T ).The considered physical system is presented in Figure 2. The effect of Carbon nanotube/Cu-nanoparticles on both fluid flow and heat transfer is investi- temperature are tabulated in Table 1.
The dimensionless quantities are as follows: The boundary conditions on the solid walls are in the following forms: 0, 1, for 0, 0 1 0, 0, for 0 , 0 0, 0, for 0 , 1 The boundary conditions on the east open side ( , 0 1 x W y = ≤ ≤ ) are in the following forms: The boundary conditions in the LBM can be implemented through the distribution function.For both flow and thermal fields, the distribution functions out of the domain are known from the streaming process.The unknown distribution functions are those toward the domain.

Result and Discussion
The presence of nanotube particles with high thermal conductivity in suspensions has different effects on temperature and flow fields.The effect of added nanotubes and Cu nanoparticles on flow and temperature fields of the Water (the base fluid) is presented in Figure 3 for the open cavity with 1 A = at different Rayligh numbers.At the present study to be sure that nanofluid behavior is completely same as single phase fluid, the volume fraction value of added nanotubes is less than 1% to make a dilute suspension, therefore it can be acceptable that the overall shape of the contours are same for both pure fluid and nanofluid with added nanotube particles.Also, the maximum value of stream-function in the flow field is presented in this figure for Water, Water-Cu and Water-Carbon nanotube nanofluid at different case studies.
The presented values for maximum stream-function shows that adding the 1% volume fraction of Carbon nanotube-nanoparticles increase the maximum value of streamfunction about 57%, 59% and 75% respectively for   The overall view of this figure shows that the increase of cavity's aspect ratio leads to reduce of flow strength.The distance between hot and cold side in an enclosure is one of most important parameter of natural convection regime, the increase of this distance change the natural convection regime to a conduction manner in the problem domain.

Conclusions
The

[ 25 ]
applied LBM to investigate the effect of Cu, Al 2 O 3 and CuO nanoparticles on mixed convection in a lid-driven cavity.Their results show that adding nanoparticles increase the rate of mixed convection heat transfer of the base fluid for all tested Reynolds numbers.Abu-Nada and Chamkha[26] conducted to study of mixed convection heat transfer of Water-Al 2 O 3 nanofluid in an inclined cavity.Heat transfer enhancement due to increase of nanoparticles volume fraction at different Richardson and Grashof numbers was presented in their results.Their results show that ave Nu on hot wall enhances about 3.3%, 8.4%, 11.9% and 17.2% respectively for 2%, 5%, 7% angle of the cavity is equal to zero and also this values are 3%, 8.5%, 12% and 17% for 2 Ri =.The investigation of the heat transfer enhancement by adding the Carbon nanotube to the base fluid in an open-ended cavity is the main aim of the present study.Also a comparison between one type of Carbon nanotubes i.e. single wall Carbon nanotubes (SWCNT) and Cu-nanoparticles is one of the major tasks of the study.In this simulation the effective conductivity and viscosity were calculated based on the new theoretical models.The presented model by Masoumi et al. [28] is applied for effective viscosity of the Cu and Carbon nanotube-nanofluid.To simulate thermal conductivity of Carbon nanotube-nanofluid a new theoretical model presented by Sabbaghzadeh and Ebrahimi By considering D2Q9 model as shown in Figure1, for applied lattice scheme for both flow and temperature fields, equilibrium distribution functions for flow field ( eq f α )

Figure 1 .
Figure 1.Discrete velocity vectors for the D2Q9 model of LBM.

Figure 2 .
Figure 2. Schematic geometry of the problem.

Figure 5 .
Figure 5. Streamlines and temperature contours of the base fluid, Carbon nanotube and Cu-nanofluid ( 1% ϕ = ) at effect of Carbon nanotube and Cu nanoparticles on natural convection heat transfer in an open-end enclosure was studied numerically.The problem was investigated at different aspect ratios of the cavity (1 4 A ≤ ≤ ) and the volume fractions of Carbon na- notube and Cu-nanoparticles ( 0 1% ϕ ≤ ≤ ) when Rayligh number vary from 10 3 to 10 5 .Some of most important results that have been achieved in this study are as follows:  Results show that adding a low value of Carbon nanotube to the base fluid led to significant enhancement of convection heat transfer. The heat transfer rate is at closed relation with thermal conductivity of the suspensions therefore the use of nanoparticles with better thermal conductivity leads to better heat transfer enhancement at base fluid. Make a comparison between Carbon nanotube and Cu-nanoparticles shows that the Carbon nanotube-nanoparticle has better performance to enhance convection rate. The aspect ratio of the cavity plays an important role on natural convection heat transfer.An increase of this parameter leads to heat transfer reduction in a 2D deep open cavity. Rayligh number loses its importance on natural convection in the high aspect ratios.
studied convection heat transfer of a nanofluid-filled trapezoidal enclosure.They investigated the effect of different volume fraction of Al 2 O 3 and Cu nanoparticles on heat transfer enhancement of Water as base fluid.They reported an enhancement of natural convection heat transfer about 5 10