Search tips
Search criteria 


Logo of nanoreslettspringer open web sitethis articlemanuscript submissionregistrationjournal front pagespringer open web site
Nanoscale Res Lett. 2013; 8(1): 178.
Published online Apr 17, 2013. doi:  10.1186/1556-276X-8-178
PMCID: PMC3660296
Simulation of forced convection in a channel with nanofluid by the lattice Boltzmann method
Nor Azwadi Che Sidik,corresponding author1 Maysam Khakbaz,1 Leila Jahanshaloo,1 Syahrullail Samion,1 and Amer Nordin Darus1
1Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, UTM Skudai, Johor, 81310, Malaysia
corresponding authorCorresponding author.
Nor Azwadi Che Sidik: azwadi/at/; Maysam Khakbaz: maysam_khakbaz/at/; Leila Jahanshaloo: leilajahanshaloo/at/; Syahrullail Samion: syahruls/at/; Amer Nordin Darus: amer/at/
Received February 11, 2013; Accepted March 29, 2013.
This paper presents a numerical study of the thermal performance of fins mounted on the bottom wall of a horizontal channel and cooled with either pure water or an Al2O3-water nanofluid. The bottom wall of the channel is heated at a constant temperature and cooled by mixed convection of laminar flow at a relatively low temperature. The results of the numerical simulation indicate that the heat transfer rate of fins is significantly affected by the Reynolds number (Re) and the thermal conductivity of the fins. The influence of the solid volume fraction on the increase of heat transfer is more noticeable at higher values of the Re.
Keywords: Nanofluid, Volume fraction, Nusselt number, Reynolds number, Lattice Boltzmann method
Recent years have eyewitnessed a blossom flourishing in the evolvement of electronics, communications, and auto-computing industries, and this bearing is irrefutably continuing in this century. The cooling of electrical, mechanical, and electronic components has become troublesome in today's fast-growing technologies. Inasmuch as the significance of heat exchangers in tremendous engineering applications, the subject of potential heat transfer enhancement in these devices has received sizeable attention in practice and research. On account of the fact that the consistency of the electronic components commodiously increases, conspicuous lack of heat transfer enhancement both in macro- and microscales of channels is realized. Encountering a fluid flow by utilizing transverse surfaces in a channel is a prevalent method that is used to intensify the rate of heat transfer from heated surfaces.
Alamyane and Mohamad [1] studied the forced convection heat transfer in a channel with extended surfaces. The effects of the Reynolds number (Re) and the fin height and spacing on the fluid flow and the heat transfer were examined. Yang et al. [2] simulated the forced convection in a parallel plate channel. Constant temperature was considered in both upper and lower walls, and a transverse object was located at the lower channel wall. The effects of the Reynolds number, the thermal conductivity ratio of the fluid, and the fin profile area on the fluid flow and the heat transfer rate were analyzed. The study results showed that the heat transfer enhancement with an increment of the Reynolds number and the thermal conductivity ratio of the fluid at various fin profiles. Yang et al. [3] numerically investigated the effect of mix convection heat transfer in an inclined parallel plate channel with a transverse object at the bottom wall. In this research, the effects of thermal conductivity, Reynolds number, the fin profile, and the channel inclination on the heat transfer rate at various Richardson numbers were examined. They discovered that the ace aspect ratio of the fin was related to the fin with utmost heat transfer at various Reynolds and Richardson numbers.
Young and Vafai [4] observed the impact of controlling parameters on the cooling of heated channels with mounted objects. Concentrating on the effect of altering the dimensions of the object, the thermal conductivity, the heating method, and the Re was embraced. They deduced that the fluid flow and heat transfer are affected by the geometry and material of the object, and a correlation for the average Nusselt number was proposed as a function of the controlling parameters.
Meinders and Hanjalic [5] experimentally investigated the effect of the cubes' arrangement on the turbulent fluid flow. They comprehended that the flow stream was affected by the distance between the objects owing to the fact of augmenting the flow velocity. Moreover, amelioration in velocity distribution and heat transfer than the staggered distribution case was found for flow over inline cubes. Yan et al. [6] experimentally investigated the influence of short surface-mounted objects at the top of a flat plate on the heat transfer enhancement. Scrutinizing was done on the effect of varies cross sections, spacing and numbers of objects, and the Reynolds number. They perceived that the heat transfer was incremented when the height of the object is comparatively equal to half of the channel height.
In an experimental investigation by Yuan et al. [7], the heat transfer and friction characteristics of a channel which were attached by winglets were examined. Heat transfer from the channel was achieved to be noticeably augmented by using winglets in comparison with conventional channels with rectangular transverse objects. For a high Reynolds number, the heat transfer was enhanced by a factor of 2.7 to 6 times of the smooth channel.
Utilizing nanofluids for the purpose of enhancing the heat transfer in thermal systems is another alternative technique [8]. The thermal performance of different types of nanofluids has been the subject of many recent studies on forced, natural, and mixed convection problems. Several explorations have studied natural convection of nanofluids in cavities [9,10]. They argued that the addition of nanoparticles in the fluid indisputably increase the natural convection heat transfer.
Chein and Huang [11] analyzed the cooling of two silicon microchannel heat sinks with a water-Cu nanofluid. The heat transfer and fraction coefficients were based on the theoretical models and the experimental correlations. They realized that the heat transfer performance of microchannels was greatly improved when nanofluids were added into base fluid as coolants without any extra pressure drop.
Recently, Santra et al. [12] numerically investigated the effect of water-Cu nanofluid through parallel plate channel in laminar forced convection. A cold nanofluid was sent through the channel, and the walls of the channel were isothermally heated. The effects of the Reynolds number and the solid volume fraction on the heat transfer were studied by considering the fluid to be Newtonian and non-Newtonian. They observed that the rate of heat transfer increased with an increase of the Reynolds number and the solid volume fraction. The increase in the heat transfer was approximately the same for both scenarios.
The lattice Boltzmann method (LBM) is another numerical method that is often used to simulate flow problems. LBM have been used for more than 2 decades as an alternative numerical technique. In LBM, it is intended to model fluids as a collection of particles, which successively undergo collision and propagation over a discrete lattice mesh. Several lattice Boltzmann models have been proposed for the incompressible Navier–Stokes equations. A collision model was proposed by Bhatnagar et al. [13] to simplify the analysis of the lattice Boltzmann equation, which leads to the so-called lattice BGK model. Remarkable efforts have been conducted by many researchers that made this numerical method more attractive for fluid dynamics modeling, e.g., [14,15]. For more details about LBM and its application, kindly refer to the aforementioned publications.
Most of the researches cited above considered the heat transfer enhancement by adding either the fin or using nanofluids. The main objective of this study is to examine both of these effects on the heat transfer performance. In general, previous works were performed to investigate different cases of nanofluid flow and heat transfer in channels with mounted objects by focusing on changing geometries, arrangement, and dimensions of the objects. However, more efforts are needed in order to optimize the controlling parameters for best heat transfer enhancement.
Problem definition
The geometry of the problem is shown in Figure Figure1.1. A cold mixture of base fluid (water) and the nanoparticles (alumina) is forced to flow into a channel that is heated from its bottom and kept at a constant high temperature, while the top wall is insulated. The channel aspect ratio is fixed at L/H = 15. The Prandtl number is taken as 7.02, and the Reynolds numbers are 10, 50, and 100, whereas the extended surfaces' height to space ratio l/S is 0.2, and the ratio between the objects' height to the channel's height l/H is 0.2.
Figure 1
Figure 1
A schematic plot of flow in a channel.
The flow is assumed as Newtonian, laminar, two-dimensional, and incompressible. In addition, it is assumed that the cold mixture of base fluid (water) and the solid spherical nanoparticles (alumina) is in thermal equilibrium, and it flows at the same velocity as a homogenous mixture.
Numerical simulation
The D2Q9 LBM model is used to simulate fluid flow in two-dimensional channel with uniform grid size of δx × δy. The lattice Boltzmann equation (known as LBGK equation) with single relaxation time can be expressed as [13]
equation M1
which can be reformulated as
equation M2
where equation M3 and τf as the single relaxation time of the fluid, fi represents the particle distribution function, ei is the particle streaming velocity, and equation M4 is the local equilibrium distribution function. For D2Q9 model equation M5 is given by [8]
equation M6
where ρ is the density of the fluid and ωi is the weight function, which has the values of equation M7, equation M8 for i = 1 to 4, and equation M9 for i = 5 to 8. The macroscopic fluid flow velocity in lattice units is represented by u. In the LBM, the fluid macroscopic quantities such as density, ρ, and flow momentum, ρu are calculated using the distribution function fi, and given by equation M10 and equation M11, respectively. The streaming speed for particles in coordinate (x and y) directions (i.e., 1 to 4, see Figure Figure2)2) can be expressed as ei = cos(π/2 (i − 1)), sin(π/2 (i − 1)), whereas particles in diagonal directions (i.e., 5 to 8 in Figure Figure2)2) have velocities of equation M12; however, the particle in the lattice center is at rest and has no streaming speed, i.e., e0 = 0.
Figure 2
Figure 2
A schematic plot showing the thermal boundary conditions of the problem.
The thermal part is simulated using another distribution function for the temperature. For instance, g is used to simulate the distribution function of the dependent variable (temperature) in the lattice Boltzmann equation, and an approach similar to that used to simulate the fluid flow is utilized to simulate the temperature distribution. In addition, the algorithm suggested by Succi [15] is adopted throughout this work. The kinetic equation for the temperature distribution function with single relaxation time is given by:
equation M13
which can be written in the form
equation M14
Where gi represents the temperature distribution function of the particles, equation M15 is the local equilibrium distribution function of the temperature, and equation M16, where τt is the single relaxation time of the temperature distribution. Thus, the equilibrium distribution function of the thermal part is given by [15]:
equation M17
where, ϕ is the macroscopic temperature and equation M18 is the speed of sound. The diffusion coefficient can be obtained as a function of the relaxation time and given by equation M19. The macroscopic temperature is then computed from:
equation M20
A uniform lattice of 100 × 1,500 is used to perform all of the simulations. However, the number of lattices was doubled to test the grid dependency results.
Since the inlet velocity of the flow is specified, the inward distribution functions should be computed at the boundary. In the D2Q9 model, the values of the distribution functions pointing out of the domain at the inlet boundary (i.e., f3, f6, f7 in Figure Figure2)2) are known from the streaming step, and the only unknowns are (f1, f5, f8) as well as the fluid density ρ. Following the work of Zou and He [16], the inlet density and the distribution functions can be obtained from:
equation M21
The unknown distribution functions are calculated using
equation M22
An extrapolation scheme is used to simulate the outlet flow condition, which can be represented as fi(Nx, t) = fi(Nx − 1, t), i = 3, 6, 7. The bounce-back scheme is used to specify the boundary conditions on solid surfaces (no-slip boundary), in which the distribution functions pointing to the fluid are equal to those pointing out of the domain. The thermal boundary conditions for this case are given in Figure Figure2.2. For constant wall temperature (the lower wall temperature is constant), the unknown functions are obtained using the following equation [15]:
equation M23
The left-hand boundary (channel inlet) is kept at a constant temperature (Dirichlet boundary condition) and set to a dimensionless value of zero. Therefore, the resulted equations of the unknown distribution functions on the left boundary are given by gi(0, t) = − gi + 2(0, t), for i = 1, 5, and 8. For the adiabatic boundary condition, the gradient of the dependent variable normal to the boundary should be zero, i.e., [partial differential] [var phi]/[partial differential] y = 0. The distribution functions are found to be in the following form [15]:
equation M24
A second-order extrapolation similar to the one given in [17] is used to obtain the values of the unknown distribution functions for the right-hand side boundary (channel outlet) as follows:
equation M25
The local Nusselt number (Nux) is computed using the following equation:
equation M26
where Lc is the characteristic length and ϕwall is the wall constant temperature. The mean temperature ϕm is given by:
equation M27
The effective density of the nanofluid is
equation M28
where ϕ is the solid volume fraction. The effective dynamic viscosity of the nanofluid given by Brinkman [18] is
equation M29
The thermal diffusivity of the nanofluid is
equation M30
The heat capacitance of the nanofluid is
equation M31
keff is the effective thermal conductivity of the nanofluid and is determined using the model proposed by Patel et al. [19]. For the two-component entity of spherical particle suspension, the model gives:
equation M32
where ks and kf are the thermal conductivities of dispersed Al2O3 nanoparticles and pure water.
equation M33
where us is the Brownian motion velocity of the nanoparticles given by:
equation M34
where kb = 1.3087×10−23JK−1 is the Boltzmann constant.
Code validation and computational results
For the purpose to ensure that the obtained results are proper and that the code is free of errors, a flow of cold air in a two-dimensional heated channel was taken as a benchmark test. Both upper and lower walls were heated. The comparisons were carried up between the dimensionless velocity and temperature fields at different locations in the channel as shown in Figures  3 and and4.4. The obtained results were found to be identical to the results of [20].
Figure 3
Figure 3
Velocity and profiles at different cross sections.
Figure 4
Figure 4
Temperature profiles at different cross sections.
Figure Figure55 shows the effect of Reynolds on the temperature profiles at the same cross sections for Re = 10, 50, and 100. The figures depicted that the temperature profiles are less sensitive to the change in Reynolds compared to the velocity profiles.
Figure 5
Figure 5
Velocity and temperature profiles at different Re.
The effects of the Reynolds number and the solid volume fraction on the heat transfer, isotherms, and streamlines are studied. Figure Figure66 presents the streamlines and the isotherms for the Al2O3-water nanofluid (ϕ = 0.05) and pure water at different Reynolds number (Re = 10, 50, and 100).
Figure 6
Figure 6
Streamlines and isotherms for the Al2O3-water nanofluid and pure water at different Reynolds number. (A) Streamline plots at (a) Re = 10, (b) Re = 50, and (c) Re = 100. (B) Isotherm plots at Re = 10 and (a) [var phi] = 0.0 and (b) [var phi] = 0.05. ( (more ...)
The streamlines show that as the Reynolds number increases, the vortices that are formed behind the fins become larger and stronger. This can be more clearly illustrated in Figure Figure55 where the horizontal velocity in the middle section between fins is presented. At Re = 10, the velocity is consistently positive. However, as the Reynolds number increases, the flow velocity becomes negative. This is an indication of flow reversal. The strong vortex at high numbers enhances the heat transfer from left face objects to right face objects and the wall between the two fins. This difference, however, becomes noticeable at higher Re.
At low Reynolds numbers, the conduction is the dominating mechanism of heat transfer. Therefore, the isotherms stretch above the fins and take a large area in the channel. As Re increases, the convection becomes the dominating mechanism, and the strong cold inlet flow pushes the isotherms near the bottom wall. The comparison between the isotherms of the nanofluid and pure water shows that in each point of the channel, the nanofluid temperature is higher than the pure water. It is due to the nanofluid's higher thermal conductivity.
The current investigation is wrapped with the analysis of the effect of the Reynolds number and percentage of nanoparticle volume fraction on the heat transfer enhancement in the channel. Figure Figure77 and Table  1 display values of average Nusselt number at various Reynolds numbers and solid volume fraction from 0% to 5%. These figures demonstrate that the Nusselt number increases with the Reynolds number for values of volume fraction tested in the present study. For example, at Re = 100, in the addition of volume fraction of 5%, the average Nusselt number increases about 17%. High Reynolds number results in high energy transport through the fluid and cause irregular motion of nanoparticle. The higher solid volume fraction further stimulates the flow and contributes to higher Nusselt number as shown in the figure. The presence of nanoparticles also increases the rate of heat transfer by conduction mode through the flow.
Figure 7
Figure 7
Average Nusselt number for various Re.
Table 1
Table 1
Average Nusselt number for various Reynolds number and solid volume fraction
LBM was applied to simulate forced convection heat transfer in two-dimensional channel including extended surfaces to investigate the effect of changing different parameters such as Reynolds number (10, 50, and 100) and nanofluid (Al2O3) volume fractions (0.0, 0.03, and 0.05). The results showed that as the Reynolds number increases, the rate of heat transfer also increases. The formation of vortices both in front and behind the objects enhances the heat transfer process. As the solid volume fraction increases, the heat transfer is enhanced for all values of the Reynolds numbers. This enhancement is more significant at high Reynolds numbers. The heat transfer rate of the fins increases with the thermal conductivity ratio of the fin to pure water. This enhancement has a finite limit. At this limit, the temperature at all surfaces of the fins approach the wall temperature. In this condition, the fins behave like constant temperature of heat sources.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
MK, LJ, and SS conceived the study and checked the grammar of the manuscript. NACS and AND drafted the manuscript. All authors read and approved the final manuscript.
This research is financially supported by the Ministry of Higher Education of Malaysia through Fundamental Research Grant Scheme, FRGS Vot no. 4L074.
  • Alamyane AA, Mohamad AA. Simulation of forced convection in a channel with extended surfaces by the lattice Boltzmann method. Comput Math Appl. 2010;8:2421–2430. doi: 10.1016/j.camwa.2009.08.070. [Cross Ref]
  • Yang MH, Yeh RH, Hwang JJ. Forced convective cooling of a fin in a channel. Energy Convers Manage. 2010;8:1277–1286. doi: 10.1016/j.enconman.2010.01.003. [Cross Ref]
  • Yang MH, Yeh RH, Hwang JJ. Mixed convective cooling of a fin in a channel. Int J Heat Mass Transfer. 2010;8:760–771. doi: 10.1016/j.ijheatmasstransfer.2009.10.012. [Cross Ref]
  • Young TJ, Vafai K. Convective cooling of a heated obstacle in a channel. Int J Heat Mass Transfer. 1998;8:3131–3148. doi: 10.1016/S0017-9310(97)00323-2. [Cross Ref]
  • Meinders ER, Hanjalic K. Experimental study of the convective heat transfer from inline and staggered configuration of two wall-mounted cubes. Int J Heat Mass Transfer. 2002;8:465–482. doi: 10.1016/S0017-9310(01)00180-6. [Cross Ref]
  • Yan WM, Hsieh RC, Soong CY. Experimental study of surface-mounted obstacle effects on heat transfer enhancement by using transient liquid crystal thermograph. J Heat Transfer. 2002;8:762–769. doi: 10.1115/1.1459729. [Cross Ref]
  • Yuan ZX, Tao WQ, Yan XT. Experimental study on heat transfer in ducts with winglet disturbances. Heat Transfer Eng. 2003;8:76–84.
  • Fakhreddine SO, Rachid B. Heterogeneous nanofluids: natural convection heat transfer enhancement. Nanoscale Res Lett. 2011;8:222–232. doi: 10.1186/1556-276X-6-222. [PMC free article] [PubMed] [Cross Ref]
  • Veeranna S, Lakshmi NS. AL2O3-based nanofluids: a review. Nanoscale Res Lett. 2011;8:456–471. doi: 10.1186/1556-276X-6-456. [PMC free article] [PubMed] [Cross Ref]
  • Saeed ZH, Seyyed HN, Elham T, Javad S. Numerical investigation of Al2O3/water nanofluid laminar convective heat transfer through triangular ducts. Nanoscale Res Lett. 2011;8:179–188. doi: 10.1186/1556-276X-6-179. [PMC free article] [PubMed] [Cross Ref]
  • Chein R, Huang G. Analysis of microchannel heat sink performance using nanofluids. Appl Therm Eng. 2005;8:3104–3114. doi: 10.1016/j.applthermaleng.2005.03.008. [Cross Ref]
  • Santra AK, Sen S, Chakraborty N. Study of heat transfer due to laminar flow of copperewater nanofluid through two isothermally heated parallel plates. Int J Therm Sci. 2009;8:391–400. doi: 10.1016/j.ijthermalsci.2008.10.004. [Cross Ref]
  • Bhatnagar PL, Gross EP, Krook M. A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Phys Rev. 1954;8:511–525. doi: 10.1103/PhysRev.94.511. [Cross Ref]
  • Peter V. Heat transfer augmentation in nanofluids via nanofins. Nanoscale Res Lett. 2011;8:154–166. doi: 10.1186/1556-276X-6-154. [PMC free article] [PubMed] [Cross Ref]
  • Succi S. Applied lattice Boltzmann method for transport phenomena, momentum, heat and mass transfer. Can J Chem Eng. 2007;8:946–947.
  • Zou Q, He X. On pressure and velocity boundary conditions for the lattice Boltzmann BGK model. Phys Fluids. 1997;8:1591–1598. doi: 10.1063/1.869307. [Cross Ref]
  • He Y, Qi C, Hu Y, Qin B, Li F, Ding Y. Lattice Boltzmann simulation of alumina-water nanofluid in a square cavity. Nanoscale Res Lett. 2011;8:184–191. doi: 10.1186/1556-276X-6-184. [PMC free article] [PubMed] [Cross Ref]
  • Brinkman HC. The viscosity of concentrated suspensions and solution. J Chem Phys. 1952;8:571–581. doi: 10.1063/1.1700493. [Cross Ref]
  • Patel HE, Sundararajan T, Pradeep T, Dasgupta A, Dasgupta N, Das SK. A micro-convection model for thermal conductivity of nanofluids. Pramana J Phys. 2005;8:863–869. doi: 10.1007/BF02704086. [Cross Ref]
  • Kays WM, Crawford ME, Weigand B. Convective Heat and Mass Transfer. 4. Boston: McGraw Hill; 2005.
Articles from Nanoscale Research Letters are provided here courtesy of