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𝑇3∞3𝜅1𝜕2𝑇𝜕𝑦2,(2.3)𝜕𝐶+⃗𝜕𝑡𝑉⋅∇𝐶=𝐷𝐵∇2𝐷𝐶+ğ‘‡ğ‘‡âˆžî‚¶âˆ‡2𝑇.(2.4) We write ⃗𝑉=(𝑢,𝑣).

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𝑇3∞3𝜌𝑐𝑝𝑓𝜅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𝑇3∞3𝜌𝑐𝑝𝑓𝑘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 𝑐1≠0.

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),Ra−1/3𝑥Nu𝑥=âˆ’ğœƒî…ž(0),Ra−1/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.

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).

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.

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.

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)).

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.

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.