Table of Contents Author Guidelines Submit a Manuscript
Journal of Applied Mathematics
Volume 2012 (2012), Article ID 648675, 18 pages
http://dx.doi.org/10.1155/2012/648675
Research Article

Lie Group Analysis of Natural Convective Flow from a Convectively Heated Upward Facing Radiating Permeable Horizontal Plate in Porous Media Filled with Nanofluid

1School of Mathematical Sciences, Universiti Sains Malaysia, 11800 Penang, Malaysia
2Department of Engineering Sciences, Pakistan Navy Engineering College, National University of Sciences and Technology, Karachi 75350, Pakistan

Received 9 February 2012; Revised 19 March 2012; Accepted 20 March 2012

Academic Editor: Srinivasan Natesan

Copyright © 2012 Md. Jashim Uddin 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

Two-dimensional, steady, laminar and incompressible natural convective flow of a nanofluid over a connectively heated permeable upward facing radiating horizontal plate in porous medium is studied numerically. The present model incorporates Brownian motion and thermophoresis effects. The similarity transformations for the governing equations are developed by Lie group analysis. The transformed equations are solved numerically by Runge-Kutta-Fehlberg fourth-fifth order method with shooting technique. Effects of the governing parameters on the dimensionless velocity, temperature and nanoparticle volume fraction as well as on the dimensionless rate of heat and mass transfer are presented graphically and the results are compared with the published data for special cases. Good agreement is found between numerical results of the present paper and published results. It is found that Lewis number, Brownian motion and convective heat transfer parameters increase the heat and mass transfer rates whilst thermophoresis decreases both heat and mass transfer rates.

1. Introduction

Nanoparticles are made from various materials, such as oxide ceramics (Al2O3, CuO), nitride ceramics (AlN, SiN), carbide ceramics (SiC, TiC), metals (Cu, Ag, Au), semiconductors, (TiO2, SiC), carbon nanotubes, and composite materials such as alloyed nanoparticlesAl70Cu30 or nanoparticle core-polymer shell composites. Nanofluids aim to achieve the maximum possible thermal properties at the minimum possible concentrations (preferably < 1% by volume) by uniform dispersion and stable suspension of nanoparticles (preferably < 10 nm) in host fluids [1]. Present heat transfer industries require high performance heat transfer equipment. The idea of improving heat transfer performance of fluids with the inclusion of solid particles was first introduced by Maxwell [2]. But, suspensions involving milli or microsized particles create problems, such as fast sedimentation, clogging of channels, high pressure drop, and severe erosion of system boundaries. To overcome these difficulties Choi [3] used ultrafine nanoparticles with base fluid termed as nanofluid. Modern material technologies facilitated the manufacturing of nanometer-sizes particles. Nanofluids have superior thermophysical properties like high thermal conductivity, minimal clogging in flow passages, long-term stability and homogeneity. Nanofluids have several industrial applications such as in electronics, automotive, and nuclear applications where efficient heat dissipation is necessary. According to Schaefer [4], nanobiotechnology is a fast growing field of research and application in many domains such as in medicine, pharmacy, cosmetics, and agroindustry. Advances in nanoelectronics, nanophotonics, and nanomagnetics have seen the arrival of nanotechnology as a distinct discipline in its own right [5].

A good number of research papers have been published on nanofluids to understand their performance so that they can be used to enhance the heat transfer in various industrial applications. A review of convective transport in nanofluids was conducted by Buongiorno [6]. Khan and Aziz [7] studied natural convection flow of nanofluid past a vertical plate with uniform heat flux. The Cheng and Minkowycz [8] problem was investigated by Nield and Kuznetsov [9] for nanofluid where the model incorporates the effect of Brownian motion and thermophoresis. Kuznetsov and Nield [10] presented a similarity solution of natural convective boundary-layer flow of a nanofluid past a vertical plate. An analytical study on the onset of convection in a horizontal layer of a porous medium with the Brinkman model and the Darcymodel filled with a nanofluid was presented by Kuznetsov and Nield [11, 12]. Godson et al. [13] presented the recent experimental and theoretical studies on convective heat transfer in nanofluids, their thermophysical properties, and applications and clarifies the challenges and opportunities for future research. Vajravelu et al. [14] studied convective heat transfer in the flow of viscous Ag water and Cu water nanofluids over a stretching surface. Noghrehabadi et al. [15] studied effect of partial slip boundary condition on the flow and heat transfer of nanofluids past stretching sheet prescribed constant wall temperature. Very recently, Aziz and Khan [16] studied similarity analysis of natural convective flow of a nanofluid over a convectively heated vertical plate.

According to previous researchers, for example, Aboeldahab and Azzam [17] radiation must be considered in calculating thermal effects in many new engineering processes occurring at high temperatures, such as the nuclear reactor cooling system, gas turbines, the various propulsion devices for aircraft, missiles, satellites, and space vehicles and various devices for space technology, underground nuclear wastes disposal, and so forth. Due to diverse applications of radiation, many investigators investigate the effect of radiation on the hydrodynamic or hydromagnetic or hydroelectric boundary layer flow over different geometries under different boundary conditions. A few examples are the papers by Cortell [18], Bataller [19], and Ishak et al. [20]. Gbadeyan et al. [21] present a numerical analysis of boundary layer flow of a nanofluid due over a linearly stretching sheet in the presence of thermal radiation. Very recently, Chamkha et al. [22] investigated mixed convective boundary-layer flow over an isothermal radiating vertical wedge placed in a porous medium filled with a nanofluid numerically using Keller box method.

Fluid flow and heat transfer in porous media have many engineering applications such as postaccidental heat removal in nuclear reactors, solar collectors, drying processes, storage of radioactive nuclear waste, heat exchangers, geothermal energy recovery and crude oil extraction, ground water pollution, thermal energy storage, building construction and flow through filtering media, separation processes in chemical industries [23]. Reviews of the fundamental theories and experiments of thermal convection in porous media with practical applications are presented in the books by Nield and Bejan [24], Vadasz [25], Vafai [26]. The classical problem of free convective flow in a porous medium near a horizontal flat plate was first investigated by Cheng and Chang [27]. After his pioneering works several authors such as Chang and Cheng [28], Shiunlin and Gebhart [29], Merkin and Zhang [30] and Chaudhary et al. [31] have extended the problem in various aspects. Gorla and Chamkha [32] presented a boundary layer analysis for the free convection flow of nanofluid over a horizontal upward facing plate in a porous medium numerically. Khan and Pop [23] extended this problem for nanofluid. Above investigators considered isothermal or isoflux thermal boundary conditions. However, the idea of using the thermal convective heating boundary condition was introduced by Aziz [33] to analyze Blasius flow. Following him, several authors, for example, Yao et al. [34], Uddin et al. [35], Magyari [36], and Yacob et al. [37] among others, used this boundary condition to study convective phenomena.

Above investigators found similarity solutions via dimensional analysis which can find only one particular type of similarity independent variable of the form 𝜂=𝑐𝑦𝑥𝑟, where 𝑟 is a numerical constant and 𝑐 is a dimensional constant [38]. However, if one deals with the governing partial differential equations by Lie group analysis, then one can obtain former similarity transformation as well as some new forms [39, 40]. Sometime it is extremely difficult to transform the PDEs to ODEs by using dimensional analysis. On the other hand, reduction of PDEs with boundary conditions to ODEs is much easier by use of Lie group analysis. The number of independent variables of PDEs can be reduced by one if the PDEs remain invariant under Lie group of transformations and the new system contains one less independent variable than the original one. This methodology can be applied (𝑛1) times to reduce a boundary value problem of PDEs having 𝑛 number of independent variables to a boundary value problem of ODEs. The solution of reduced equations is much easier than the solutions of the original system of PDEs [41]. Hence, Lie group of transformations may be considered as the generalization of dimensional analysis. It is successfully applied in many areas such as in mathematical physics, applied and theoretical mechanics and applied mathematics and in the transport phenomena [42, 43]. Avramenko et al. [44] presented that the symmetrical properties of the turbulent boundary-layer flows and other turbulent flows are studied utilizing the Lie group theory technique.

The aim of the present study is to investigate the effect of thermophoresis, the Brownian motion, radiation and the thermal convective boundary condition on the boundary layer flow of a nanofluid over an upward facing radiating permeable horizontal plate numerically. A possible application of this problem is in the design of furnace where the transfer of heat from surfaces occurs simultaneously by radiation and convection. Also, the interaction of solar radiation with the earth’s surface fabricates complex free convection patterns and hence complicates the studies associated with the weather forecasting and marine environment for predicting free convection patterns in oceans and lakes. Using similarity transformations developed by Lie group analysis, the governing partial differential equations are reduced to a set of coupled nonlinear ordinary differential equations with the corresponding boundary conditions. The effect of emerging flow controlling parameters on the dimensionless axial velocity, the temperature, the nanoparticle volume fraction, the rate of heat transfer, and the rate of nanoparticle volume fraction is investigated and shown graphically and discussed.

2. Formulation of the Problem

We consider a two-dimensional (𝑥,𝑦) laminar free convective boundary layer flow past a permeable upward facing horizontal plate with radiation effects in a porous media filled with a nanofluid (Figure 1). The temperature 𝑇 and the nanoparticle volume fraction 𝐶 take constant values 𝑇𝑤 and 𝐶𝑤 at the boundary whilst 𝑇 and 𝐶 at free stream. Bottom of the plate is heated by convection from a hot fluid at temperature 𝑇𝑓 which gives a variable heat transfer coefficient 𝑓(𝑥). It is assumed that 𝑇𝑓>𝑇𝑤>𝑇. The Oberbeck-Boussinesq approximation is employed. The following four field equations represent the conservation of mass, momentum, thermal energy, and nanoparticles, respectively. The field variables are 𝑉: Darcy velocity vector, 𝑇: the temperature, and 𝐶: the nanoparticle volume fraction [23]:𝑉=0,(2.1)𝜌𝑓𝜀𝜕𝑉𝜇𝜕𝑡=𝑃𝐾𝑉+𝐶𝜌𝑃𝜌+(1𝐶)𝑓1𝛽𝑇𝑇𝑔,(2.2)𝜌𝐶𝑃𝑓𝜕𝑇+𝜕𝑡𝑉𝑇=𝑘𝑚2𝑇+𝜀𝜌𝐶𝑃𝑃𝐷𝐵𝐷𝐶𝑇+𝑇𝑇+𝑇𝑇16𝜎1𝑇33𝜅1𝜕2𝑇𝜕𝑦2,(2.3)𝜕𝐶+𝜕𝑡𝑉𝐶=𝐷𝐵2𝐷𝐶+𝑇𝑇2𝑇.(2.4) We write 𝑉=(𝑢,𝑣).

648675.fig.001
Figure 1: Coordinate system and flow model.

Here 𝜌𝑓 is the density of the base fluid, 𝜇 is the dynamic viscosity of the base fluid, 𝛽 is the volumetric expansion coefficient of nanofluid, 𝜌𝑝is the density of the nanoparticles, (𝜌𝐶𝑃)𝑓 is the heat effective heat capacity of the fluid, (𝜌𝐶𝑃)𝑃 is the effective heat capacity of the nanoparticle material, 𝑘𝑚 is effective thermal conductivity of the porous medium, 𝜀 is the porosity, 𝐾 is permeability of the porous media, 𝑔 is the gravitational acceleration, 𝜎1 is the Sefan-Boltzman constant, and 𝑘1 is the Rosseland mean absorption coefficient. Here 𝐷𝐵 stands for the Brownian diffusion coefficient and 𝐷𝑇 stands for the thermophoretic diffusion coefficient. To ignore an advective term and a Forchheimer quadratic drag term in the momentum equation, we assumed that the flow is slow.

Consider a steady state flow. In keeping with the Oberbeck-Boussinesq approximation and an assumption that the nanoparticle concentration is dilute, and with a suitable choice for the reference pressure, we can linearize the momentum equation and write (2.2) as𝜇0=𝑃𝐾𝜌𝑉+𝑃𝜌𝑓𝐶𝐶+1𝐶𝜌𝑓𝛽𝑇𝑇𝑔.(2.5) Making the standard boundary layer approximation based on an order of magnitude analysis to neglect the small order terms, we have the governing equations𝜕𝑢𝜕𝑥+𝜕𝑣𝜕𝑦=0,(2.6)𝜕𝑃𝜕𝑥𝜇=𝐾𝑢,(2.7)𝜕𝑃𝜕𝑦𝜇=𝐾𝑣+1𝐶𝜌𝑓𝑔𝛽𝑇𝑇𝜌𝑃𝜌𝑓𝑔𝐶𝐶,(2.8)𝑢𝜕𝑇𝜕𝑥+𝑣𝜕𝑇𝜕𝑦=𝛼𝑚𝜕2𝑇𝜕𝑦2𝐷+𝜏𝐵𝜕𝐶𝜕𝑦𝜕𝑇𝜕𝑦+𝐷𝑇𝑇𝜕𝑇𝜕𝑦2+16𝜎1𝑇33𝜌𝑐𝑝𝑓𝜅1𝜕2𝑇𝜕𝑦2,(2.9)𝑢𝜕𝐶𝜕𝑥+𝜈𝜕𝐶𝜕𝑦=𝐷𝐵𝜕2𝐶𝜕𝑦2+𝐷𝑇𝑇𝜕2𝑇𝜕𝑦2,(2.10) where 𝛼𝑚=𝑘𝑚/(𝜌𝑐𝑃)𝑓is the thermal diffusivity of the fluid and 𝜏=𝜀(𝜌𝐶𝑃)𝑝/(𝜌𝐶𝑃)𝑓 is a parameter.

The boundary conditions are taken to be [35]𝑣=𝑣𝑤𝑥,𝑘𝜕𝑇𝜕𝑦=𝑓𝑥𝑇𝑓𝑇𝑤,𝐶=𝐶𝑤,at𝑦=0,𝑢0,𝑇𝑇,𝐶𝐶as𝑦.(2.11) Here 𝑣𝑤(𝑥): mass transfer velocity. The following nondimensional variables are introduced to make (2.6)–(2.11) dimensionless𝑥=𝑥𝐿𝑅𝑎,𝑦=𝑦𝐿,𝑢=𝑢𝐿𝛼𝑚𝑅𝑎,𝑣=𝑣𝐿𝛼𝑚,𝜃=𝑇𝑇Δ𝑇,𝜙=𝐶𝐶,Δ𝐶Δ𝑇=𝑇𝑓𝑇,Δ𝐶=𝐶𝑤𝐶,(2.12) where 𝐿is the plate characteristic length and 𝑅𝑎=𝑔𝐾𝛽(1𝐶)Δ𝑇𝐿/(𝛼𝑚𝜈)is the Rayleigh number. A stream function 𝜓 defined by𝑢=𝜕𝜓𝜕𝑦,𝑣=𝜕𝜓𝜕𝑥,(2.13) is introduced into (2.6)–(2.11) to reduce number of dependent variables and equations. Note that (2.6) is satisfied identically. We are then left with the following three dimensionless equations:𝜕2𝜓𝜕𝑦2+𝜕𝜃𝜕𝑥𝑁𝑟𝜕𝜙𝜕𝑥=0,𝜕𝜓𝜕𝑦𝜕𝜃𝜕𝑥𝜕𝜓𝜕𝑥𝜕𝜃𝜕𝜕𝑦2𝜃𝜕𝑦2𝑁𝑏𝜕𝜃𝜕𝑦𝜕𝜙𝜕𝑦𝑁𝑡𝜕𝜃𝜕𝑦2𝜕𝑅2𝜃𝜕𝑦2=0,𝐿𝑒𝜕𝜓𝜕𝑦𝜕𝜙𝜕𝑥𝜕𝜓𝜕𝑥𝜕𝜙𝜕𝜕𝑦2𝜙𝜕𝑦2𝑁𝑡𝜕𝑁𝑏2𝜃𝜕𝑦2=0.(2.14)

The boundary conditions in (2.11) become𝜕𝜓𝜕𝑥=𝐿𝑣𝑤(𝑥)𝛼𝑚,𝜕𝜃𝜕𝑦=𝑓𝐿𝑘(1𝜃),𝜙=1at𝑦=0,𝜕𝜓𝜕𝑦0,𝜃0,𝜙0as𝑦.(2.15) Five parameters in (2.14) are 𝑁𝑡,𝑁𝑏,𝑁𝑟,𝑅, and 𝐿𝑒 and they stand for the thermophoresis parameter, the Brownian motion parameter, the buoyancy ratio parameter, radiation parameter, and the Lewis number, respectively, which are defined by𝑁𝑡=𝜏𝐷𝑇Δ𝑇/𝛼𝑚𝑇,𝑁𝑏=𝜏𝐷𝐵Δ𝐶/𝛼𝑚,𝐿𝑒=𝛼𝑚/𝐷𝐵,𝜌𝑁𝑟=𝑃𝜌𝑓Δ𝐶/𝜌𝑓𝛽1𝐶Δ𝑇,𝑅=16𝜎1𝑇33𝜌𝑐𝑝𝑓𝑘1𝛼𝑚.(2.16)

3. Symmetries of the Problem

By applying Lie group method to (2.14), the infinitesimal generator for the problem can be written as𝑋=𝜉1𝜕𝜕𝑥+𝜉2𝜕𝜕𝑦+𝜏1𝜕𝜕𝜓+𝜏2𝜕𝜕𝜃+𝜏3𝜕𝜕𝜙(3.1) where the coordinates (𝑥,𝑦,𝜓,𝜃,𝜙)transformed into the coordinates (𝑥,𝑦,𝜓,𝜃,𝜙). The infinitesimals 𝜉1,𝜉2,𝜏1,𝜏2, and 𝜏3 satisfies the following first order linear differential equations𝑑𝑥𝑑𝜀=𝜉1𝑥,𝑦,𝜓,𝜃,𝜙,𝑑𝑦𝑑𝜀=𝜉2𝑥,𝑦,𝜓,𝜃,𝜙,𝑑𝜓𝑑𝜀=𝜏1𝑥,𝑦,𝜓,𝜃,𝜙,𝑑𝜃𝑑𝜀=𝜏2𝑥,𝑦,𝜓,𝜃,𝜙,𝑑𝜃𝑑𝜀=𝜏3𝑥,𝑦,𝜓,𝜃,𝜙.(3.2) Using commercial software Maple 13, it was found that the forms of the infinitesimals are𝜉1=𝑐1𝑥+𝑐2,𝜉2=23𝑐1𝑦+𝑐3,𝜏1=13𝑐1𝜓+𝑐6,𝜏2=𝑐4𝜃,𝜏3=𝑐5𝜙,(3.3) where 𝑐𝑖(𝑖=1,2,3,4,5,6) are arbitrary constants. Hence, the equations admit six finite parameters Lie group transformations. It is apparent that the parameters 𝑐2,𝑐3, and c6correspond to the translation in the variables 𝑥,𝑦, and 𝜓, respectively. It is also observed that the parameters 𝑐1,𝑐4, and 𝑐5 correspond to the scaling in the variables 𝑥,𝑦,𝜓,𝜃, and 𝜙, respectively. The generators corresponding to the infinitesimal given by (3.3) are as follows:𝑋1𝜕=𝑥+2𝜕𝑥3𝑦𝜕+1𝜕𝑦3𝜓𝜕𝜕𝜓,𝑋2=𝜕𝜕𝑥,𝑋3𝜕=𝑦𝜕𝑦,𝑋4𝜕=𝜃𝜕𝜃,𝑋5=𝜕𝜕𝜙.(3.4) We consider scaling transformations and hence set 𝑐2=𝑐3=𝑐6=0.

Thus the infinitesimals become𝜉1=𝑐1𝑥2,𝜉2=23𝑐1𝑦,𝜏1=13𝑐1𝜓,𝜏2=𝑐4𝜃,𝜏3=𝑐5𝜙.(3.5) In terms of differentials, we have𝑑𝑥𝑐1𝑥=𝑑𝑦(2/3)𝑐1𝑦=𝑑𝜓(1/3)𝑐1𝜓=𝑑𝜃𝑐4𝜃=𝑑𝜙𝑐5𝜙.(3.6) Here 𝑐10.

3.1. Similarity Transformations

From (3.6), 𝑑𝑥/𝑐1𝑥=𝑑𝑦/(2/3)𝑐1𝑦, which on integration 𝑦𝑥2/3=constant=𝜂(say).(3.7a) Similarly, 𝑑𝑥/𝑐1𝑥=𝑑𝜓/(1/3)𝑐1𝜓leads to𝜓𝑥1/3=constant=𝑓(𝜂)(say),thatis,𝜓=𝑥1/3𝑓(𝜂),(3.7b) where 𝑓 is arbitrary function of 𝜂.

Equations 𝑑𝑥/𝑐1𝑥=𝑑𝜃/𝑐4𝜃and 𝑑𝑥/𝑐1𝑥=𝑑𝜙/𝑐5𝜙 lead to𝜃=𝑥𝑐4/𝑐1𝜃(𝜂),𝜙=𝑥𝑐5/𝑐1𝜙(𝜂).(3.7c)Thus from (3.7a)–(3.7c) we obtain the following similarity transformations:𝑦𝜂=𝑥2/3,𝜓=𝑥1/3𝑓(𝜂),𝜃=𝑥𝑐4/𝑐1𝜃(𝜂),𝜙=𝑥𝑐5/𝑐1𝜙(𝜂).(3.8) Now, to make sure that 𝜃0,𝜙0 as 𝜂, set 𝑐4=𝑐5=0.

Hence the similarity transformations are𝑦𝜂=𝑥2/3,𝜓=𝑥1/3𝑓(𝜂),𝜃=𝜃(𝜂),𝜙=𝜙(𝜂).(3.9) Thus the velocity component 𝑢,𝑣 can be expressed as𝑓𝑢=𝑥1/31,𝑣=3𝑥2/3𝑓2𝜂𝑓,(3.10) where primes indicate differentiation with respect to similarity independent variable 𝜂. It is worth citing that the similarity transformations in (3.9) are consistent with the well-known similarity transformations reported in the paper of Cheng and Chang [27] for 𝜆=0 in their paper, which support the validity of our analysis.

3.2. Governing Similarity Equations

Substituting the transformations in (3.9) into the governing (2.14) leads to the following nonlinear system of ordinary differential equations:𝑓23𝜂𝜃𝑁𝑟𝜙=0,(1+𝑅)𝜃+13𝑓𝜃+𝑁𝑏𝜃𝜙+𝑁𝑡𝜃2𝜙=0,+𝐿𝑒3𝑓𝜙+𝑁𝑡𝜃𝑁𝑏=0(3.11) subject to the boundary conditions𝑓(0)=𝑓𝑤,𝜃[](0)=𝐵𝑖1𝜃(0),𝜙(0)=1,𝑓()=𝜃()=𝜙()=0.(3.12) Here 𝑁𝑏=0 means there is no thermal transport due to buoyancy effects created as a result of nanoparticle concentration gradients and 𝑓𝑤=𝐿𝑣𝑤/3𝛼𝑚,𝑓𝑤>0 corresponds to suction and 𝑓𝑤<0 corresponds to injection, 𝐵𝑖=𝐿𝑓/𝑘 is the Biot number. It is mentioned that, for a true similarity solution, we must have𝑓=𝑓0𝑥2/3,𝑣𝑤=𝑣𝑤0𝑥2/3,(3.13) where (𝑓)0 and (𝑣𝑤)0 are constants.

4. Comparisons with the Literature

It is worth citing that in case of impermeable nonradiating plate (𝑓𝑤=𝑅=0) and for isothermal plate (𝐵𝑖), the problem under consideration reduces to the problem which has been recently investigated by Khan and Pop [23] and Gorla and Chamkha [32]. It is also worth mentioning that in case of impermeable non-radiating plate (𝑓𝑤=𝑅=0) and for constant wall temperature (𝐵𝑖), in the absence of buoyancy force (𝑁𝑟=0), thermophoresis (𝑁𝑡=0) and in the absence of Brownian motion (𝑁𝑏=0), the problem under consideration reduces to the problem which was investigated by Cheng and Chang [27] for 𝜆=0 in their paper. It is further noted that in case of non-radiating plate (𝑅=0), the problem under consideration reduces to the problem which was recently investigated by Uddin et al. [35].

5. Physical Quantities

The parameters of engineering interest are the local skin friction factor 𝐶𝑓𝑥, the local Nusselt number Nu𝑥, and the local Sherwood number Sh𝑥, respectively. Physically, 𝐶𝑓𝑥 indicates wall shear stress, Nu𝑥 indicates the rate of heat transfer whilst Sh𝑥 indicates the rate of mass transfer. These quantities can be calculated from following relations:𝐶𝑓𝑥=2𝜇𝜌𝑈2𝑟𝜕𝑢𝜕𝑦𝑦=0,Nu𝑥=𝑥𝑇𝑓𝑇𝜕𝑇𝜕𝑦𝑦=0,Sh𝑥=𝑥𝐶𝑤𝐶𝜕𝐶𝜕𝑦𝑦=0.(5.1) By substituting from (2.12) and (3.9) into (5.1), it can be shown that physical quantities can be put in the following dimensionless form:Ra𝑥Pr𝐶𝑓𝑥=2𝑓(0),Ra1/3𝑥Nu𝑥=𝜃(0),Ra1/3𝑥Sh𝑥=𝜙(0).(5.2) Here Ra𝑥=𝑔𝐾𝛽(1𝐶)Δ𝑇𝑥/(𝛼𝑚𝑣)is the local Rayleigh number, Pr=𝑣/𝛼𝑚is the Prandt number for porous media, and 𝑈𝑟=(1𝐶)𝑔𝐾𝛽Δ𝑇/𝛼𝑚 is reference velocity in porous media. Note that the local skin friction factor, the local Nusselt number, and the local Sherwood number are directly proportional to the numerical values of 𝑓(0),𝜃(0) and 𝜙(0), respectively.

6. Results and Discussion

The set of coupled nonlinear similarity Equations (3.11) with boundary conditions in (3.12) forms a two-point boundary value problem and has been solved numerically using an efficient Runge-Kutta-Fehlberg fourth-fifth order numerical method under Maple 14. To highlight the important features of the flow velocity, temperature, nanoparticle volume fraction, the heat transfer rate, and the nanoparticle volume fraction transfer rate, the obtained numerical results are displayed graphically. Numerical computations are done for0𝑅5,1𝑓𝑤1,0.1𝑁𝑏0.5,0.1𝑁𝑡0.5, 0𝑁𝑟0.5, 0𝐵𝑖5.0and1𝐿𝑒5. The results of the dimensionless heat transfer rates𝜃(0) and the dimensionless nanoparticle volume fraction rate -𝜙(0) are compared with the most recent results reported by Gorla and Chamkha [32] for special case in Table 1 and found to be in excellent agreement with each of values of 𝑁𝑟,𝑁𝑏, and 𝑁𝑡. This supports the validity of our other graphical results for dimensionless velocity, temperature, nanoparticle volume fraction, heat transfer, and nanoparticle volume fraction transfer rates.

tab1
Table 1: Comparison of present results with Gorla and Chamkha [32] for different values of buoyancy and nanofluid parameters.
6.1. Velocity Profiles

Figures 2 and 3 exhibit the dimensionless axial velocity profiles for various values of the emerging flow controlling parameters. Dimensionless velocity and corresponding velocity boundary layer thickness are decreased with increasing values of the mass transfer velocity both for radiating (𝑅=5) and nonradiating (𝑅=0) plate. In Figure 2(a), it is found that dimensionless velocity increases with the increasing of the radiation parameter. It is apparent from Figure 2(b) that the dimensionless velocity decreases with rising of the buoyancy ratio parameter. The velocity is reduced with the suction; reverse phenomena are observed in case of the injection, as expected. The dimensionless velocity is elevated with rising of the Biot number and the Lewis number (Figure 3).

fig2
Figure 2: Effect of (a) radiation parameter, (b) buoyancy ratio parameter on the dimensionless velocity for different values of suction/injection parameter.
fig3
Figure 3: Effect of (a) Biot number, (b) Lewis number on the dimensionless velocity for different values of suction/injection parameter.
6.2. Temperature Profiles

Variation of the dimensionless temperature and corresponding thermal boundary layer thickness with radiation parameter, suction/injection parameter, the Biot number, thermophoresis, and Brownian motion parameters is shown in Figures 4 and 5, respectively.Temperature is increased with the increasing of radiation and Boit number (Figure 4). Physically, higher Biot number increases nanoparticle volume fraction as nanoparticle volume fraction distribution is driven by temperature distribution. The fluid on the right surface of the plate is heated up by the hot fluid on the left surface of the plate, making it lighter and flowing faster.

fig4
Figure 4: Effect of (a) radiation parameter, (b) Biot number on the dimensionless temperature for different values of suction/injection parameter.
fig5
Figure 5: Effect of (a) Brownian motion parameter, (b) thermophoresis parameter on the dimensionless temperature for different values of suction/injection parameter.

Note that the temperature increases with the increasing of the Brownian motion and thermophoresis parameters when the plate is permeable or not (Figure 5). From Figures 4 and 5, it is apparent that like regular fluid suction/injection parameter reduces the dimensionless temperature as expected.

6.3. Nanoparticle Volume Fraction Profiles

Figure 6 illustrates the impact of the controlling parameters on the dimensionless nanoparticle volume fraction inside the corresponding boundary layer. Dimensionless nanoparticle volume fraction is reduced due the enhance of both the radiation and the Lewis number when the plate is permeable or not (Figures 6(a) and 6(b)). Finally, from Figure 6, we found that the suction/injection parameter reduces the dimensionless nanoparticle volume fraction as in the case of regular fluid.

fig6
Figure 6: Effect of (a) radiation parameter, (b) Lewis number on the dimensionless nanoparticle volume fraction for different values of suction/injection parameter.
6.4. Heat Transfer Rate

The effect of various controlling parameters on the dimensionless heat transfer rate from a permeable horizontal upward facing plate with the thermal convective boundary condition in porous media is shown in Figure 7. From Figure 7(a), it is noticed that the dimensionless heat transfer rate decreases with an increase in thermophoresis and radiation parameter whilst it increases with the increasing of the suction velocity. It is further found from Figure 7(b) that the dimensionless heat transfer rate decreases with an increase in thermophoresis and buoyancy ratio parameter for permeable plate. We also noticed that heat transfer rate is a decreasing function of the radiation parameter (Figure 7(a)) whilst it is increasing function of the Boit number (Figure 7(b)).

fig7
Figure 7: Variation of local heat transfer rate with (a) Brownian motion, suction/injection, and radiation parameters, (b) thermophoresis, buoyancy ratio, and Biot number.
6.5. Nanoparticle Volume Fraction Rate

Figure 8 shows the effect of the radiation, the suction, thermophoresis, buoyancy ratio, and Lewis number parameters on the dimensionless nanoparticle volume fraction transfer rate from a permeable horizontal upward facing radiating plate in porous media. From Figure 8(a), we observed that the dimensionless nanoparticle volume fraction rate increases with an increase in Brownian motion, suction, and the radiation parameter. It is also found from Figure 8(b) that the dimensionless nanoparticle volume fraction rate decreases with an increase in both the thermophoresis and buoyancy ratio parameter for a permeable plate. It is further observed form Figure 8(b) that the Lewis number increases the dimensionless nanoparticle volume fraction rate, as in regular fluid.

fig8
Figure 8: Variation of local mass transfer rate with (a) Brownian motion, radiation, and suction/injection parameters, (b) Lewis number, buoyancy ratio, and thermophoresis parameters.

7. Conclusions

We studied numerically a 2-D steady laminar viscous incompressible boundary layer flow of a nanofluid over an upward facing horizontal radiating permeable plate placed in the porous media considering the thermal convective boundary condition. The governing boundary layer equations are transformed into highly nonlinear coupled ordinary differential equations using similarity transformations developed by Lie group analysis, before being solved numerically. Following conclusions are drawn:(i)the dimensionless velocity, the temperature, and the concentration decrease in case of the suction and increase in case of the injection; the phenomenon is reversed,(ii)the Brownian motion, radiation, thermophoresis, and buoyancy ratio parameters decrease the heat transfer rate whilst the suction parameter and the Biot number enhance the heat transfer rate,(iii)the radiation, Lewis number, Brownian motion, and the suction parameters cause to enhance nanoparticle volume fraction rate whilst thermophoresis and buoyancy ratio parameters lead to decreasing nanoparticle volume fraction rate.

References

  1. S. Kakaç, S. Özerinç, and A. G. Yazıcıoğlu, “Enhanced thermal conductivity of nanofluids: a state-of-the-art review,” Microfluid Nanofluid, vol. 8, pp. 145–170, 2010. View at Google Scholar
  2. J. C. Maxwell, A Treatise on Electricity and Magnetism, Clarendon Press, Oxford, UK, 2nd edition, 1981.
  3. S. U. S. Choi, “Enhancing thermal conductivity of fluids with nanoparticles,” in Developments and Applications of Non-Newtonian Flows, D. A. Siginer and H. P. Wang, Eds., pp. 99–105, American Society of Mechanical Engineers, Baltimore, Md, USA, 1995. View at Google Scholar
  4. H. Schaefer, Nanoscience: The Science of the Small in Physics, Engineering. Chemistry, Biology and Medicine (Nanoscience and Technology), Springer, 2010.
  5. R. P. E. Sharma, Nanostructuring Operations in Nanoscale Science and Engineering, Kal.Mcgraw Hill, 2010.
  6. J. Buongiorno, “Convective transport in nanofluids,” ASME Journal of Heat Transfer, vol. 128, pp. 240–250, 2006. View at Google Scholar
  7. W. A. Khan and A. Aziz, “Natural convection flow of a nanofluid over a vertical plate with uniform surface heat flux,” International Journal of Thermal Science, vol. 50, pp. 1207–1214, 2011. View at Google Scholar
  8. P. Cheng and W. J. Minkowycz, “Free convection about a vertical flat plate embedded in a porous medium with application to heat transfer from a dike,” Journal of Geophysics Research, vol. 82, no. 14, pp. 2040–2044, 1977. View at Publisher · View at Google Scholar
  9. D. A. Nield and A. V. Kuznetsov, “The Cheng-Minkowycz problem for natural convective boundary-layer flow in a porous medium saturated by a nanofluid,” International Journal of Heat Mass Transfer, vol. 52, pp. 5792–5795, 2009. View at Google Scholar
  10. A. V. Kuznetsov and D. A. Nield, “Natural convective boundary-layer flow of a nanofluid past a vertical plate,” International Journal of Thermal Science, vol. 49, pp. 243–247, 2010. View at Google Scholar
  11. A. V. Kuznetsov and D. A. Nield, “Thermal instability in a porous medium layer saturated by a nanofluid: brinkman model,” Transport in Porous Media, vol. 81, no. 3, pp. 409–422, 2010. View at Publisher · View at Google Scholar
  12. A. V. Kuznetsov and D. A. Nield, “Effect of local thermal non-equilibrium on the onset of convection in a porous medium layer saturated by a nanofluid,” Transport in Porous Media, vol. 83, pp. 425–436, 2010. View at Google Scholar
  13. L. B. Godson, D. Raja, D. L. Mohan, and S. Wongwisesc, “Enhancement of heat transfer using nanofluids—an overview,” Renewable and Sustainable Energy Reviews, vol. 14, pp. 629–641, 2010. View at Google Scholar
  14. K. Vajravelu, K. V. Prasad, J. Lee, C. Lee, I. Pop, and R. A. van Gorder, “Convective heat transfer in the flow of viscous Ag water and Cu water nanofluids over a stretching surface,” International Journal of Thermal Sciences, vol. 50, pp. 843–851, 2011. View at Google Scholar
  15. A. Noghrehabadi, P. Rashid, and G. Mohammad, “Effect of partial slip boundary condition on the flow and heat transfer of nanofluids past stretching sheet prescribed constant wall temperature,” International Journal of Thermal Sciences, vol. 54, pp. 253–261, 2012. View at Google Scholar
  16. A. Aziz and W. A. Khan, “Natural convective boundary layer flow of a nanofluid past a convectively heated vertical plate,” International Journal of Thermal Science, vol. 52, pp. 83–90, 2012. View at Google Scholar
  17. E. M. Aboeldahab and G. E. A. Azzam, “Thermal radiation effects on MHD flow past a semi-infinite inclined plate in the presence of mass diffusion,” Heat and Mass Transfer, vol. 41, no. 12, pp. 1056–1065, 2005. View at Google Scholar
  18. R. Cortell, “Suction, viscous dissipation and thermal radiation effects on the flow and heat transfer of a power-law fluid past an infinite porous plate,” Chemical Engineering Research and Design, vol. 89, pp. 85–93, 2011. View at Google Scholar
  19. R. C. Bataller, “Radiation effects in the Blasius flow,” Applied Mathematics and Computation, vol. 198, no. 1, pp. 333–338, 2008. View at Publisher · View at Google Scholar
  20. A. Ishak, N. A. Yacob, and N. Bachok, “Radiation effects on the thermal boundary layer flow over a moving plate with convective boundary condition,” Meccanica, vol. 46, no. 4, pp. 795–801, 2011. View at Publisher · View at Google Scholar
  21. J. A. Gbadeyan, M. A. Olanrewaju, and P. O. Olanrewaju, “Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition in the presence of magnetic field and thermal radiation,” Australian Journal of Basic and Applied Sciences, vol. 5, no. 9, pp. 1323–1334, 2011. View at Google Scholar
  22. A. J. Chamkha, S. Abbasbandy, A. M. Rashad, and K. Vajravelu, “Radiation effects on mixed convection over a wedge embedded in a porous medium filled with a nanofluid,” Transport in Porous Media, vol. 91, pp. 261–279, 2012. View at Google Scholar
  23. W. A. Khan and I. Pop, “Free convection boundary layer flow past a horizontal flat plate embedded in a porous medium filled with a nanofluid,” ASME Journal of Heat Transfer, vol. 133, no. 9, Article ID 094501, 4 pages, 2011. View at Publisher · View at Google Scholar
  24. D. A. Nield and A. Bejan, Convection in Porous Media, Springer, New York, NY, USA, 3rd edition, 2006.
  25. P. Vadasz, Emerging Topics in Heat and Mass Transfer in Porous Media, Springer, New York, NY, USA, 2008.
  26. K. Vafai, Porous Media: Applications in Biological Systems and Biotechnology, CRC Press, New York, NY, USA, 2010.
  27. P. Cheng and I. D. Chang, “Buoyancy induced flows in a saturated porous medium adjacent to impermeable horizontal surfaces,” International Journal of Heat Mass Transfer, vol. 19, pp. 1267–1272, 1976. View at Google Scholar
  28. I. D. Chang and P. Cheng, “Matched asymptotic expansions for free convection about an impermeable horizontal surface in a porous medium,” International Journal of Heat Mass Transfer, vol. 26, pp. 163–173, 1983. View at Google Scholar
  29. D. S. Shiunlin and B. Gebhart, “Buoyancy-induced flow adjacent to a horizontal surface submerged in porous medium saturated with cold water,” International Journal of Heat Mass Transfer, vol. 29, pp. 611–623, 1986. View at Google Scholar
  30. J. H. Merkin and G. Zhang, “On the similarity solutions for free convection in a saturated porous medium adjacent to impermeable horizontal surfaces,” Wärme-und Stoffübertragung, vol. 25, pp. 179–184, 1990. View at Google Scholar
  31. M. A. Chaudhary, J. H. Merkin, and I. Pop, “Natural convection from a horizontal permeable surface in porous medium-numerical and asymptotic solutions,” Transport in Porous Media, vol. 22, no. 3, pp. 327–344, 1996. View at Publisher · View at Google Scholar
  32. R. S. R. Gorla and A. Chamkha, “Natural convective boundary layer flow over a horizontal plate embedded in a porous medium saturated with a nanofluid,” Journal of Modern Physics, vol. 2, pp. 62–71, 2011. View at Google Scholar
  33. A. Aziz, “A similarity solution for laminar thermal boundary over a flat plate with a convective boundary condition,” Communication in Nonlinear Science and Numerical Simulation, vol. 14, pp. 1064–1068, 2009. View at Google Scholar
  34. S. Yao, T. Fang, and Y. Zhong, “Heat transfer of a generalized stretching/shrinking wall problem with convective boundary conditions,” Communication in Nonliear Science and Numerical Simulation, vol. 16, pp. 752–760, 2011. View at Google Scholar
  35. M. J. Uddin, W. A. Khan, and A. I. M. Ismail, “Free convection boundary layer flow from a heated upward facing horizontal flat plate embedded in a porous medium filled by a nanofluid and with a convective boundary condition,” Transport in Porous Media, vol. 92, no. 3, pp. 867–881, 2012. View at Google Scholar
  36. E. Magyari, “Comment on “a similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition” by A. Aziz, Comm.Non.Sci. Numer. Simul.' (2009).14, 1064'68,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, pp. 599–601, 2011. View at Google Scholar
  37. N. A. Yacob, A. Ishak, and I. Pop, “Falkner-Skan problem for a static or moving wedge in nanofluids,” International Journal of Thermal Sciences, vol. 50, pp. 133–139, 2011. View at Google Scholar
  38. G. W. Bluman and S. Kumei, Symmetries and Differential Equations, vol. 81 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1989.
  39. M. Jalil, S. Asghar, and M. Mushtaq, “Lie group analysis of mixed convection flow with mass transfer over a stretching surface with suction or injection,” Mathematical Problems in Engineering, vol. 2010, Article ID 264901, 14 pages, 2010. View at Publisher · View at Google Scholar
  40. M. Pakdemirli and M. Yurusoy, “Similarity transformations for partial differential equations,” SIAM Review, vol. 40, no. 1, pp. 96–101, 1998. View at Publisher · View at Google Scholar
  41. G. W. Bluman and S. C. Anco, Symmetry and Integration Methods for Differential Equations, Springer, New York, NY, USA, 2009.
  42. P. J. Olver, Applications of Lie Groups to Differential Equations, vol. 107 of Graduate Texts in Mathematics, Springer, New York, NY, USA, 1986.
  43. H. Ibragimov and V. F. Kovalev, Approximate and Re-Norm Group symmetry, Springer, New York, NY, USA, 2009.
  44. A. A. Avramenko, S. G. Kobzar, I. V. Shevchuk, A. V. Kuznetsov, and L. T. Iwanisov, “Symmetry of turbulent boundary-layer flows: investigation of different eddy viscosity models,” Acta Mechanica, vol. 151, pp. 1–14, 2001. View at Google Scholar