Research Article  Open Access
Analysis of Heat Transfer in Berman Flow of Nanofluids with Navier Slip, Viscous Dissipation, and Convective Cooling
Abstract
Heat transfer characteristics of a Berman flow of water based nanofluids containing copper (Cu) and alumina (Al_{2}O_{3}) as nanoparticles in a porous channel with Navier slip, viscous dissipation, and convective cooling are investigated. It is assumed that the exchange of heat with the ambient surrounding takes place at the channel walls following Newton’s law of cooling. The governing partial differential equations and boundary conditions are converted into a set of nonlinear ordinary differential equations using appropriate similarity transformations. These equations are solved analytically by regular perturbation methods with series improvement technique and numerically using an efficient RungeKutta Fehlberg integration technique coupled with shooting scheme. The effects of the governing parameters on the dimensionless velocity, temperature, skin friction, pressure drop, and Nusselt numbers are presented graphically and discussed quantitatively.
1. Introduction
The study of fluid flow and heat transfer between two porous boundaries has gained tremendous attention of researchers due to its wide applications in engineering and industrial processes. Some of the practical interests include problems dealing with transpiration cooling where the walls of a channel containing heated fluid are protected from overheating by passing cooler fluid over the exterior surface of the channel; fluid flow occurring during the separation of isotopes of Uranium235 and Uranium238 by gaseous diffusion in order to produce fuel for nuclear reactors; controlling boundary layer flow over aircraft wings by injection or suction of fluid out of or into the wing; lubrication of porous bearings; petroleum technology; ground water hydrology; seepage of water in river beds; purification and filtration processes; methods of decreasing rates of heat transfer in combustion chambers exhaust nozzles and porous walled flow reactors, and so forth. In a pioneering work, Berman [1] presented an exact solution of the NavierStokes equations that describes the steady twodimensional flow of an incompressible viscous fluid along a channel with parallel rigid porous walls, under the action of uniform suction or injection of fluid at the surface. Sellars [2] extended Berman’s work to high suction Reynolds number. Yuan [3] considered the flow in a channel with porous walls. He obtained solutions for small suction and injection values and asymptotic solution valid at large injection values. Terrill [4, 5] gave an exact series solution for the fully developed laminar flow in a pipe of circular cross section with porous wall driven by a spatially variable and time independent suction or injection at the surface. Studies on developing flow in porouswalled ducts with suction and injection effects were carried out by Sorour et al. [6] and Zaturska et al. [7]. In view of the above interests, several researchers have also investigated the heat transfer problems between two permeable parallel walls under different physical situations [8–10].
Moreover, the applications of conventional heat transfer fluids such as water and glycol mixture in engineering flow processes are limited due to their low thermal properties. A potential solution to improve these thermal properties is to add nanoparticles into the conventional fluids, hence forming the socalled nanofluids as coined by Choi [11]. Nanofluids contain thermally conducting submicron solid particles and have great potential as a highenergy carrier. The nanoparticles such as copper, alumina, titania, and copper oxide, unlike largersized particles, can be suspended stably within the conventional fluids without settling out of suspension. Nanofluids are free from numerous problems such as abrasion, clogging, and high pressure loss and are considered to be the nextgeneration working fluids in modern heat transfer technologies [12]. Experimental results [13–15] have shown that, even with small solid volume fraction of nanoparticles (usually less than 5%), the thermal conductivity of heat transfer fluids can be enhanced by 10–50%. Several authors [16–18] have also theoretically investigated the heat transfer enhancement of nanofluids under different physical conditions.
Meanwhile, advances in the manufacture of microdevices have enabled experimental investigation of fluid flow in nano and microscale, and many experimental results have provided evidences to support the slip condition [19, 20]. In order to describe the slip characteristics of fluid on the solid surface, Navier [21] introduced a more general boundary condition, namely, the fluid velocity component tangential to the solid surface, relative to the solid surface, which is proportional to the shear stress on the fluidsolid interface. The proportionality is called the slip length which describes the slipperiness of the surface. However, from the literature survey, it is found that no study has been conducted on the heat transfer characteristics of Berman flow of nanofluids with Navier slip, viscous dissipation, or convective cooling at the walls. Hence, the present study is an attempt in this direction. The flow of water base nanofluids containing copper (Cu) and alumina (Al_{2}O_{3}) as nanoparticles in a uniformly porous wall channel with Navier slip, viscous dissipation, and convective heat exchange with the ambient surrounding is investigated. Sections 2–4 give in more details the model nonlinear governing equations together with the analytical and numerical solution techniques employed to tackle the problem. In Section 5, we present graphically and discuss the main features of the flow and heat transfer characteristics in a range of governing parameters. Final conclusions are drawn in Section 6. To the best of our knowledge, the results of this paper are new and they have not been published before.
2. Problem Formulation
Consider a twodimensional steady flow of a viscous incompressible water base nanofluids containing copper (Cu) and alumina (Al_{2}O_{3}) as nanoparticles in a uniformly porous wall channel. The channel wall is subjected to Navier slip and convectively exchange heat with the ambient surrounding. We choose the Cartesian coordinates system in such a way that the axis is taken along the channel and the axis is normal to it as shown in Figure 1.
The governing equations which are those of conservation of mass, momentum, and energy are [1–5, 8–10] where are the velocity components of the nanofluid in the directions, respectively, uniform wall suction velocity, is the pressure, is the nanofluid temperature, is the effective dynamic viscosity of the nanofluid [22], is the effective thermal conductivity of the nanofluid [23], and and are the nanofluid density and the heat capacitance of the nanofluid, respectively, which are given by [13–18]. where is the nanoparticles solid volume fraction, is the reference density of the fluid fraction, is the reference density of the solid fraction, is the viscosity of the fluid fraction, is the thermal conductivity of the fluid fraction, is the specific heat at constant pressure, and is the thermal conductivity of the solid volume fraction.
Due to the symmetric nature of the flow, the boundary conditions at the channel centerline and at the porous wall may be written as [1–4] where is the ambient surrounding temperature, is the coefficient of heat transfer, and is the Navier slip coefficient. Introduce the stream function and vorticity into the governing equations (1)–(4) as follows:
After eliminating the pressure from (2) and (3), we obtain
The following dimensionless variables and parameters are introduced into (8), together with their corresponding boundary conditions: and we obtain with where is the flow Reynolds number such that represents wall suction and represents wall injection; Ec is the Eckert number; Bi is the Biot number; Pr is the base fluid Prandtl number; , , , and can be easily determined from the thermophysical properties of the base fluid and the nanoparticles; and is the Navier slip parameter such that corresponds to no slip, while full lubrication is described in the limit . We seek a similarity form of solution due to Berman [1]; that is,
Equations (10)(11) together with the boundary conditions in (12) then become where . The dimensionless fluid axial pressure gradient is given as follows: where
Other quantities of practical interest in this study are the local skin friction coefficient and the local Nusselt number Nu, which are defined as where is the skin friction and is the heat flux at the channel walls which are given by
Using (9) and (13), we substitute (22) into (21) and obtain
In the following section, the boundary value problem (14)–(18) was solved analytically using regular perturbation method and numerically by the RungeKuttaFehlberg method with shooting technique [24]. The results are utilised to compute the fluid pressure gradient, local skin friction, and local Nusselt number as highlighted in (20) and (23).
3. Perturbation Method
Due to the nonlinear nature of the model equations (14)–(18), it is convenient to form a power series expansion in the parameter ; that is,
Substituting the solution series in (24) into (14)–(18) and collecting the coefficients of like powers of , we obtain the following.
Zeroth Order. Consider with
Higher Order . Consider with
The equations are solved iteratively and the series solutions for the velocity and temperature fields are given as where . We remark here that if we set the parameters and in (29), we will recover the solution for the classical case of conventional fluid given in [1–4]. Using a computer symbolic algebra package (MAPLE), several terms of the above solution series in (29)–(31) are obtained. From (29)–(31) together with (20) and (23), we obtained the series solutions for the skin friction, Nusselt number, and axial pressure gradient as follows: where the expression for , is given in the appendix. We are aware that these power series solutions are valid for very small parameter values of . However, using HermitePadé approximation technique (see Makinde [25]) that is based on the series summation and improvement method, the usability of the extended solution series is improved beyond small parameter values of .
4. Numerical Procedure
An efficient finite difference approach based on RungeKuttaFehlberg method with shooting technique [24] has been employed to numerically solve the coupled nonlinear ordinary differential equations (14)–(16) subject to the boundary conditions (17)(18) for different values of governing parameters. The boundary value problem is first transformed into an initial value problem (IVP). Let
Substituting (33) into (14)–(18), we obtain a system of firstorder differential equations, respectively, as follows: subject to the initial conditions:
By applying the shooting method with the NewtonRaphson algorithm to guess the unspecified conditions , , , and in (35), the resulting initial value problem is then integrated numerically until the boundary conditions at are achieved. The step size is taken as and the convergence criteria were set to 10^{−7}.
5. Results and Discussion
The flow of water base nanofluids containing Cu and Al_{2}O_{3} as nanoparticles and their heat transfer characteristics in a symmetrically porous channel with Navier slip and convective cooling at the surface are investigated. The governing partial differential equations and the corresponding boundary conditions are converted into a set of nonlinear ordinary differential equations and tackled both analytically using the perturbation method coupled with series improvement technique and numerically using RungeKutta Fehlberg integration technique coupled with shooting scheme. Thermophysical properties of base fluid and nanoparticles are presented in Table 1. For pure water, the momentum diffusivity is dominant and convection is very effective in transferring within the fluid in comparison to pure conduction. Following [16–18], we take in the numerical computation. Note that when , no nanoparticle is present in the based fluid (water). The solid volume fraction in the base fluid is taken as to 0.3 (i.e., ranging from 0 to 30 percent). In order to get a clear insight into the entire flow structure and thermal development, we have assigned numerical values to other parameters encountered in the problem. Numerical solutions are displayed in Tables 2 and 3 together with Figures 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, and 23. The numerical values of suction driven normal velocity profile are displayed in Table 2. In the absence of nanoparticles and Navier slip at the channel walls, the results agreed well with the one already in the literature (see [1–3]) for the flow of conventional fluid in a symmetrical porous wall channel. However, Cuwater nanofluid is affected more by the combined effects of suction and Navier slip at wall in comparison to Al_{2}O_{3}water nanofluid. Table 3 shows the perfect agreement between the series solution and the numerical solution for the axial velocity profiles along the channel centerline with increasing concentration of nanoparticles for both Cuwater and Al_{2}O_{3}water nanofluids.



5.1. Velocity Profiles with Parameter Variation
The effects of parameter variation on both the axial and normal velocity components are displayed in Figures 2–9. Generally, it is interesting to note that the effect of Navier slip is to cause flow reversal at the channel walls. In Figure 2, it is observed that Cuwater moves faster along the channel centerline region and subjected to higher flow reversal at the wall as compared to Al_{2}O_{3}water in the presence of suction. With Cuwater as the working nanofluid, it is observed that increasing nanoparticles volume fraction concentration from 0 to 30% increases both the axial velocity along the centerline region and the flow reversal at the channel walls as shown in Figure 3. Similar trend is observed in Figures 45 with a growing suction parameter and Navier slip parameter. Figures 6–9 show the nanofluids normal velocity profiles. With suction and Navier slip, the Cuwater moves faster to the wall as compared to Al_{2}O_{3}water as illustrated in Figure 6. A further increase in normal velocity towards the walls is observed with a growing in suction parameter, Navier slip parameter, and nanoparticles volume fraction as shown in Figures 7–9.
5.2. Temperature Profiles with Parameter Variation
Figures 10–16 illustrate the nanofluids temperature profiles across the channel with different parameter variation. Generally, a decrease in the fluid temperature near the channel walls is observed due to convective heat loss to ambient surrounding. It is noteworthy that the temperature of Cuwater nanofluid is generally higher than that of Al_{2}O_{3}water nanofluid under the same flow condition as shown in Figure 10. In Figure 11, it is observed that the nanofluid temperature decreases with growing in nanoparticles volume fraction. Similar effect of a decrease in nanofluid temperature is observed in Figure 12 with Cuwater as working nanofluid as the Biot number increases. This is expected, since an increase in Biot number indicates a rise in convective cooling due to heat loss to the ambient surrounding from the channel walls. Meanwhile, a combine increase in the suction, Navier slip, and viscous dissipation as shown in Figures 13–15 causes an increase in the nanofluid temperature. This may be attributed to the fact that , , and Ec increase the internal heat generation within the fluid due velocity gradient increases, leading to a rise in temperature. Figure 16 elucidates the temperature profiles with increasing axial distance along the channel. The nanofluid temperature decreases within the channel centerline region and increases near the wall region with increasing axial distance. Moreover, it is interesting to note that at point within the channel, the nanofluid temperature is independent of the axial distance.
5.3. Skin Friction, Pressure Gradient, and Nusselt Number
Figures 1718 depict the skin friction profiles for both Cuwater and Al_{2}O_{3}water nanofluids at the channel walls. The skin friction generally increases with an increase in nanoparticles volume fraction; however, it is noticed that the skin friction produced by Cuwater is more intense than the one produced by Al_{2}O_{3}water as shown in Figure 17. This is expected since the velocity gradient of Cuwater near the channel walls is higher than that of Al_{2}O_{3}water. In Figure 18, it is observed that the skin friction generally increases with an increase in Navier slip. Meanwhile, a growing in suction increases the skin friction, while a growing in injection decreases the skin friction. The pressure drop along the channel is illustrated in Figures 19 and 20. For both Cuwater and Al_{2}O_{3}water nanofluids, the pressure drop increases with increasing nanoparticles volume fraction. Interestingly, the pressure drop produced by Al_{2}O_{3}water is slightly higher than that of Cuwater as shown in Figure 19. Figure 20 shows a general increase in pressure drop with a rise in Navier slip. A growing in suction decreases the pressure drop, while a growing in injection increases the pressure drop along the channel. Figures 21–23 elucidate the heat transfer rate at the channel walls with different parameter variation. It is observed that the wall heat transfer rate (Nu) decreases with an increase in the nanoparticles volume fraction as shown in Figure 21. A slight increase in Al_{2}O_{3}water Nusselt number is noticed as compared to Cuwater Nusselt number. In Figure 22, it is seen that the increase in the Navier slip parameter results in increase in Nusselt number. This may be due to a rise in the nanofluid temperature gradient at the channel walls. Meanwhile, as the suction increases, the heat flux at the wall increases, while a decrease in wall heat is observed with a rise in injection. The strength of the wall heat flux is enhanced with increasing axial distance, viscous dissipation, and convective cooling as illustrated in Figure 23. This can be attributed to a rise in the temperature gradient due to convective heat exchange with the ambient along the channel walls.
6. Conclusion
The combined effects of viscous dissipation, Navier slip, and convective cooling on Berman flow and heat transfer of water base nanofluids containing Cu and Al_{2}O_{3} as nanoparticles are investigated. The nonlinear model problem is tackled both analytically using perturbation series method and numerically using RungeKutta Fehlberg integration technique coupled with shooting scheme. We summarize below some of the essential features of physical interest from the above analysis.(i)Cuwater nanofluid moves faster with enhanced flow reversal at the walls as compared to Al_{2}O_{3}water nanofluid.(ii)Nanofluids velocity and flow reversal at the walls increase with suction, , and .(iii)Cuwater produced higher temperature as compared to Al_{2}O_{3}water. The nanofluids temperature increases with suction, , and Ec, but it decreases with Bi and .(iv)The skin friction produced by Cuwater is more intense than that of Al_{2}O_{3}water. The skin friction increases with suction , , and but decreases with injection .(v)The pressure drop produced by Al_{2}O_{3}water in more than that of Cuwater. The pressure drop is enhanced by injection, , and but decreases by suction.(vi)The Nusselt number increases with suction, , Bi, Ec, and , but it decreases with injection and . A slight increase in Nu for Al_{2}O_{3}water is noticed as compared to Cuwater.
Appendix
Nomenclature
:  Velocity components 
:  Coordinates 
:  Nanofluid thermal conductivity 
:  Prandtl number 
:  local Biot number 
:  Ambient temperature 
:  Dimensionless stream function 
:  Wall suction velocity 
:  Temperature 
:  Dimensionless temperature 
:  Reynolds number 
:  Specific heat at constant pressure 
:  Solid fraction thermal conductivity 
:  Base fluid thermal conductivity 
:  Dimensionless axial velocity 
:  Eckert number 
:  Dimensionless axial coordinate 
:  Axial pressure gradient coefficient 
:  Channel half width 
:  Nusselt number 
:  Skin friction coefficient. 
:  Stream function 
:  Dimensionless temperature 
:  Nanofluid dynamic viscosity 
:  Nanofluid thermal diffusivity 
:  Dimensionless normal coordinate 
:  Nanofluid density 
:  Solid fraction density 
:  Base fluid kinematic viscosity 
:  Base fluid dynamic viscosity 
:  Solid volume fraction parameter 
:  Base fluid density 
:  Slip coefficient 
:  Navier slip parameter 
:  Vorticity 
:  Dimensionless stream function 
:  Dimensionless vorticity 
:  Dimensionless temperature. 
Conflict of Interests
The authors declare that there is no conflict of interests.
References
 A. S. Berman, “Laminar flow in channels with porous walls,” Journal of Applied Physics, vol. 24, pp. 1232–1235, 1953. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 J. R. Sellars, “Laminar flow in channels with porous walls at high suction Reynolds numbers,” Journal of Applied Physics, vol. 26, no. 4, pp. 489–490, 1955. View at: Google Scholar
 S. W. Yaun, “Further investigation of laminar flow in channels with porous walls,” Journal of Applied Physics, vol. 27, no. 3, pp. 267–269, 1956. View at: Google Scholar  MathSciNet
 R. M. Terrill, “An exact solution for flow in a porous pipe,” Journal of Applied Mathematics and Physics, vol. 33, no. 4, pp. 547–552, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. M. Terrill, “Laminar flow in a porous tube,” Journal of Fluids Engineering, vol. 105, pp. 303–306, 1983. View at: Google Scholar
 M. M. Sorour, M. A. Hassab, and S. Estafanous, “Developing laminar flow in a semiporous twodimensional channel with nonuniform transpiration,” International Journal of Heat and Fluid Flow, vol. 8, pp. 44–54, 1987. View at: Google Scholar
 M. B. Zaturska, P. G. Drazin, and W. H. H. Banks, “On the flow of a viscous fluid driven along a channel by suction at porous walls,” Fluid Dynamics Research, vol. 4, no. 3, pp. 151–178, 1988. View at: Google Scholar
 R. J. Pederson and R. B. Kinney, “Entranceregion heat transfer for laminar flow in porous tubes,” International Journal of Heat and Mass Transfer, vol. 14, pp. 159–161, 1971. View at: Google Scholar
 O. D. Makinde, “Extending the utility of perturbation series in problems of laminar flow in a porous pipe and a diverging channel,” Australian Mathematical Society B. Applied Mathematics, vol. 41, no. 1, pp. 118–128, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. D. Raithby, “Laminar heat transfer in the thermal entrance region of circular tubes and twodimensional rectangular ducts with wall suction and injection,” International Journal of Heat and Mass Transfer, vol. 14, pp. 224–243, 1971. View at: Google Scholar
 S. U. S. Choi, “Enhancing thermal conductivity of fluids with nanoparticles,” in Proceedings of the ASME International Mechanical Engineering Congress and Exposition, pp. 99–105, ASME, San Francisco, Calif, USA, 1995. View at: Google Scholar
 B. C. Pak and Y. I. Cho, “Hydrodynamic and heat transfer study of dispersed fluids with submicron metallic oxide particles,” Experimental Heat Transfer, vol. 11, no. 2, pp. 151–170, 1998. View at: Google Scholar
 D. Wen and Y. Ding, “Experimental investigation into the pool boiling heat transfer of aqueous based calumina nanofluids,” Journal of Nanoparticle Research, vol. 7, no. 2, pp. 265–274, 2005. View at: Google Scholar
 J. Wang, J. Zhu, X. Zhang, and Y. Chen, “Heat transfer and pressure drop of nanofluids containing carbon nanotubes in laminar flows.,” Experimental Thermal and Fluid Science, vol. 44, pp. 716–721, 2013. View at: Google Scholar
 Y. Xuan and Q. Li, “Investigation on convective heat transfer and flow features of nanofluids,” Journal of Heat Transfer, vol. 125, no. 1, pp. 151–155, 2003. View at: Google Scholar
 H. F. Oztop and E. AbuNada, “Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids,” International Journal of Heat and Fluid Flow, vol. 29, pp. 1326–1336, 2008. View at: Google Scholar
 O. D. Makinde, “Computational modelling of nanofluids flow over a convectively heated unsteady stretching sheet,” Current Nanoscience, vol. 9, pp. 673–678, 2013. View at: Google Scholar
 S. Kakac and A. Pramuanjaroenkij, “Review of convective heat transfer enhancement with nanofluids,” International Journal of Heat and Mass Transfer, vol. 52, no. 13/14, pp. 3187–3196, 2009. View at: Google Scholar
 P. Huang and K. S. Breuer, “Direct measurement of slip length in electrolyte solutions,” Physics of Fluids, vol. 19, Article ID 028104, 2007. View at: Publisher Site  Google Scholar
 M. J. Martin and D. I. Boyd, “Momentum and heat transfer in a laminar boundary layer with slip flow,” Journal of Thermophysics and Heat Transfer, vol. 20, no. 4, pp. 710–719, 2006. View at: Google Scholar
 C. L. M. H. Navier, “Mémoire sur les lois du mouvement des fluides,” Mémoires de l'Académie Royale des Sciences de l'Institut de France, vol. 6, pp. 389–440, 1823. View at: Google Scholar
 H. C. Brinkman, “The viscosity of concentrated suspensions and solution,” The Journal of Chemical Physics, vol. 20, pp. 571–581, 1952. View at: Google Scholar
 J. C. Maxwell, Electricity and Magnetism, Clarendon, Oxford, UK, 3rd edition, 1904.
 T. Y. Na, Computational Methods in Engineering Boundary Value Problems, Academic Press, New York, NY, USA, 1979. View at: MathSciNet
 O. D. Makinde, “HermitePadé approach to thermal radiation effect on inherent irreversibility in a variable viscosity channel flow,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2330–2338, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2014 O. D. Makinde 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.