In this paper, we investigate the effect of vertical throughflow on the onset of convection in a horizontal layer of a non-Darcy porous medium saturated with a nanofluid. A normal mode analysis is used to find solutions for the fluid layer confined between parallel plates with free-rigid boundaries. The criterion for the onset of stationary and oscillatory convection is derived. The analysis incorporates the effects of Brownian motion, thermophoresis and a convective boundary condition. The effects of the concentration Rayleigh number, Lewis number, Darcy number and modified diffusivity ratio on the stability of the system are investigated.
Keywords:throughflow; convective boundary condition; oscillatory convection; linear stability analysis
In the last several years, an innovative technique for improving the heat transfer characteristics by adding ultra fine metallic particles in common fluids such as water and oil has been investigated. The term nanofluid refers to these kinds of fluids which have applications in automotive industries, energy saving devices, nuclear rectors, etc., Choi . Nanoparticles also have medical applications including cancer therapy and nano-drug delivery (Shaw and Murthy ). Buongiorno  noted that the absolute velocity of nanoparticles can be viewed as the sum of the base fluid velocity and a relative (slip) velocity. He concluded that in the absence of turbulence, Brownian diffusion and thermophoresis would be important. A lot of work has been done on nanofluids; see, for instance, Nadeem et al., Malik et al., Nadeem et al. and the review by Das et al..
The effect of a magnetic field on flow and heat transfer problems is important in industrial applications, such as in the buoyant upward gas-liquid flow in packed bed electrodes (Takahashi and Alkire ), sodium oxide-silicon dioxide glass melt flows (Guloyan ), reactive polymer flows in heterogeneous porous media , electrochemical generation of elemental bromine in porous electrode systems (Qi and Savinell ).
The convective boundary condition is more general and realistic in engineering and industrial processes such as transportation cooling processes, material drying, etc. It therefore seems appropriate to use the convective boundary condition to study other boundary layer flow situations. Aziz  studied the Blasius flow over a flat plate with a convection thermal boundary condition. Ishak  investigated the effects of suction and injection on a flat surface with convective boundary condition. The study of the boundary layer flows over flat surfaces under convective surface boundary condition has attracted the attention of many researchers, such as Makinde and Aziz , Yao et al..
Nanofluids have great potential as coolants due to their enhanced thermal conductivities. The enhancement of effective thermal conductivity was confirmed by experiments conducted by many researchers (Masuda et al.). Instability of nanofluids in natural convection was studied by Tzou . Tzou  studied thermal instability of nanofluids in natural convection. Thermal instability in nanofluids in a porous medium has been a topic of interest due to potential applications of such flows in food and chemical processes, petroleum industry, bio-mechanics and geophysical problems. Buongiorno’s model was applied to the problem of the onset of instability in a porous medium layer saturated with a nanofluid by Nield and Kuznetsov [19,20]. Kuznetsov and Nield  used the Brinkman model to study thermal instability in a horizontal porous layer saturated with a nanofluid. Other related studies of thermal instability in a porous medium saturated with a nanofluid include those by Kuznetsov and Nield  and Nield and Kuznetsov [23-26].
In this study, we extend the work by Nield and Kuznetsov  to a non-Darcy porous medium saturated with a nanofluid. We have considered a convective boundary condition in place of an isothermal condition. The effect of the Biot number, magnetic parameter, Brinkman-Darcy parameter on thermal instability has been studied.
2 Mathematical formulation
We consider a horizontal layer of porous medium confined between the planes and where the z-axis is vertically upwards. Each boundary wall is assumed to be permeable to the throughflow and perfectly thermally conducting. Radiation heat transfer between the sides of walls is negligible when compared with other modes of heat transfer. The size of nanoparticles is small as compared to the pore size of the matrix. The nanoparticles are spherical, and the nanofluid is incompressible and laminar. It is assumed that nanoparticles are suspended in the nanofluid using either a surfactant or surface charge technology, preventing the agglomeration and deposition of these on the porous medium. The porosity of the medium is denoted by ϵ and the permeability by K. The temperatures at the lower and upper wall are and with . The nanoparticle volume fractions are at the lower wall and at the upper wall, and it is assumed that the difference is small in comparison with . A uniform magnetic field of strength B is imposed normal to the plate. Using the modified Brinkman model and the Oberbeck-Boussinesq approximation, the conservation equations for mass, momentum, energy and nanoparticles are as follows:
where is the velocity vector, ρ is the density of the fluid, t is the time, p is the hydraulic pressure, ϕ is the volume fraction of nanoparticles, is the density of nanoparticles, β is the coefficient of thermal expansion, is effective viscosity, μ is viscosity and is electric conductivity. In the energy equation, is the heat capacity of the fluid in the porous medium, is the heat capacity of nanoparticles and is thermal conductivity. In the equation of continuity for nanoparticles, is the Brownian diffusion coefficient, given by the Einstein-Stokes equation and is the thermophoretic diffusion coefficient of nanoparticles.
We note that in equation (2),
is the convective derivative. We introduce non-dimensional variables as
where is thermal diffusivity of the fluid. It is assumed that the wall is heated by convection from a hot fluid with temperature and heat transfer coefficient . Under these conditions, the thermal field is written as
while the boundary conditions for velocity and nanoparticle concentration are
After substituting equations (5) into (1)-(4), the resulting nondimensional equations are written as follows (ignoring primes):
The boundary conditions are written as
The dimensionless parameters in equations (8)-(12) are the Prandtl number Pr, the Darcy number Da, the Vadasz number Va, the density Rayleigh number Rm, the Rayleigh-Darcy number Ra, the concentration Rayleigh number Rn, the Brinkman-Darcy number , the magnetic parameter M, the Lewis number Le, the modified diffusivity ratio , the modified particle-density increment , the Peclet number Q and the Biot number Bi. These parameters are defined, respectively, by
We note here that the parameter Rm is a measure of the basic static pressure gradient.
2.1 Basic solution
A time-independent quiescent solution is obtained in the z direction only and has the form
Assumptions that Le is very large (of order 102 to 103, see Buongiorno ), equations (9)-(11) now reduce to
Solving equations (15) and (16) with boundary condition (12) gives
where and . In the limit , we obtain
As , the thermal boundary condition at the lower plate reduces to (isothermal condition), and the base solution becomes
the same results were obtained by Nield and Kuznetsov .
2.2 Perturbation solution
To study the stability of the system, we superimpose infinitesimal perturbations on the basic state solution,
Substituting (21) in equations (8)-(11), and linearizing by neglecting products of primed quantities, we get the following equations:
with the boundary conditions
The derivatives of and are
For regular fluids, the parameters Rn, and are zero and the third term in equation (25) is absent since . For , and in the absence of a magnetic field, the equations reduce to the familiar Horton-Roger-Lapwood problem with throughflow. Taking the curl of equation (23) and simplifying, we obtain
3 Normal modes and stability analysis
Differential equations (24), (25), (29) and boundary conditions (26) constitute a linear boundary value problem that can be solved using the method of normal modes. The perturbation quantities are of the form
where l and m are wave numbers in the x and y directions and n is the growth rate of the disturbances. Substituting into the differential equations, we obtain
subject to the boundary conditions
Here and is a dimensionless resultant wavenumber. We use a Galerkin-type weighted residuals method to obtain approximate solutions. We choose trial functions , , ; , that satisfy the boundary conditions and consider
Substituting equation (35) into equations (31)-(33) gives a system of 3N algebraic equations in 3N unknowns. The vanishing of the determinant of coefficients produces an eigenvalue equation for the system. Regarding Ra as the eigenvalue, we find Ra in terms of the other parameters. It is interesting to note that the convective boundary condition is applicable for rigid-free and rigid-rigid boundary conditions. In this study we mainly focus on rigid-free boundary and discuss the case of stationary and oscillatory convection. The vanishing of the shear-stresses tangent to the surface and continuity equation gives the boundary conditions
We assume that the solutions to W, Θ and Φ are of the form
which satisfies the boundary conditions. Using the above transformation in equations (31)-(33) and integrating with respect to z from to , the system of equations may be written as
where and , , , and are
For the non-trivial solution, the determinant of the augmented matrix is equal to zero, s is a dimensionless growth factor. We put , where ω is real and is the dimensionless frequency.
3.1 Stationary convection
In the case of non-oscillatory convection, . We considered the first approximation, i.e., , which gives the non-oscillatory stability boundary as
From equation (39), we get the stationary Rayleigh number as
where and are given by
The critical cell size at the onset of instability is obtained when
Solving the above equation, we get a polynomial in of order 2 and the critical is written as
where the values of , and are given by
We calculate the corresponding critical Rayleigh number using the above critical value of α for stationary convection.
3.2 Oscillatory convection
In case of the oscillatory convection, . The stability boundaries are
In order for ω to be real, it is necessary that
From equations (43) and (44), Ra (independent of ω) and ω (independent of Ra) are written as
The critical cell size at the onset of instability is obtained by solving . The critical value of α for oscillatory convection is calculated from
The critical Rayleigh number is obtained using the above critical value of α for oscillatory convection.
4 Results and discussion
The critical Rayleigh numbers for stationary and oscillatory convection are calculated from equations (40) and (46), respectively. The influence of the Biot number, magnetic parameter, Darcy number, porosity of the medium on the stationary Rayleigh number and oscillatory Rayleigh number are shown in Figures 1-6, where we have also shown the critical stationary and oscillatory Rayleigh numbers for different Darcy and Biot numbers. In the present problem, we mainly focused on the influence of Biot number on stationary and oscillatory critical Rayleigh numbers.
Figure 1. Variation of the stationary Rayleigh number with wave number for different value of (a)Bi,M( , ) and (b)Da,ε( , ) with , , , , , .
The effects of the Biot number, magnetic parameter, Darcy number and porosity on the stationary Rayleigh number are shown in Figure 1, where it is evident that increases with the magnetic parameter but decreases with the Biot number. Hence the magnetic parameter exerts a stabilizing influence on the stationary convection regime, but Biot number destabilizes stationary convection. The Darcy number has a stabilizing effect on stationary convection, while the porosity parameter has a destabilizing influence on the stationary convection for fixed Biot numbers. This finding is in line with the results reported by Chand and Rana .
Figure 2 shows the influence of the modified diffusivity ratio λ, the modified particle-density increment and the Lewis number on the stationary Rayleigh number. The stationary Rayleigh number increases with , and the Lewis number, and this helps to stabilize the stationary convection regime. On the other hand, λ destabilizes stationary convection. Thus, in the absence of the Darcy-Brinkman and magnetic parameters, the stability boundary depends on and as earlier suggested by Nield and Kuznetsov . The throughflow (due to Peclet number Q) also assists in increasing the critical stationary Rayleigh number.
Figure 2. Variation of the stationary Rayleigh number with wave number for different value of (a)λ, ( , ) and (b) ,Le( , ) with , , , , , , , , , .
The influence of the Biot number, magnetic parameter, Peclet number and Vadasz number on the oscillatory Rayleigh number is shown in Figure 3. The critical oscillatory Rayleigh number increases with the magnetic parameter while decreasing with parameters Bi, Va and Q. Hence the magnetic parameter stabilizes the oscillatory convection, while the other parameters are destabilizing to the oscillatory regime.
Figure 3. Variation of the oscillatory Rayleigh number with wave number for different value of (a)Bi,M( , ) and (b)Q,Va( , ) with , , , , , .
The critical stationary and oscillatory Rayleigh number as a function of the Darcy number is shown in Figure 4. The critical stationary Rayleigh number decreases with the Biot number. For , the critical Rayleigh number has a minimum value that depends on the Biot number as shown in Figure 5. For larger Biot numbers (isothermal condition), the critical Rayleigh number is small. The critical value of is 64.9325, 44.4513 and 42.3685 for , respectively. The critical values, on the other hand, are 562.3439, 286.1277 and 258.3124 for , respectively.
Figure 4. Influence of the Biot number on (a) and (b) as a function ofDawith , , , , , , , .
Figure 5. Effect of the Biot number on critical as a function ofDawith , , , , , , , .
We have defined Rn so that it is positive when the particle density increases upwards (the destabilizing situation). From Figure 6, we note that Ra takes a negative value when Rn is sufficiently large. In this case, the destabilizing effect of nanoparticle concentration is so large that the bottom of the fluid layer must be cooled relative to the top to produce a state of neutral stability as earlier found by Kuznetsov and Nield  in the absence of a magnetic field and higher Biot numbers. In the present problem, a state of neutral stability appeared when .
Figure 6. Effect of the Biot number on critical as a function ofDawith , , , , , , , .
In this study we used linear stability to investigate the onset of thermal instability in a non-Darcy porous medium saturated with a nanofluid, and with a convective boundary condition. We have determined the effects of various embedded parameters such as the Biot number, the magnetic parameter and the Darcy number on the critical Rayleigh number for the onset of both oscillatory and stationary thermal instabilities. We have shown that increasing the Darcy number and the magnetic parameter has the effect of increasing the critical Rayleigh number for the onset of thermal instabilities in the case of stationary convection, while increasing the Biot number and the porosity is destabilizing to the stationary regime. The modified diffusivity ratio, particle density increment and the Lewis number help to stabilize stationary convection. Oscillatory convection was found not to be as sensitive to the fluid and physical parameters as stationary convection.
The authors declare that they have no competing interests.
All the authors were involved in carrying out this study.
The authors wish to thank the University of KwaZulu-Natal for financial support.
Buongiorno, J: Convective transport in nanofluids. J. Heat Transf.. 128, 240–250 (2006). Publisher Full Text
Takahashi, K, Alkire, R: Mass transfer in gas-sparged porous electrodes. Chem. Eng. Commun.. 38, 209–227 (1985). Publisher Full Text
Guloyan, YA: Chemical reactions between components in the production of glass-forming melt. Glass Ceram.. 60, 233–235 (2003). Publisher Full Text
Liu, H, Thompson, KE: Numerical modelling of reactive polymer flow in porous media. Comput. Chem. Eng.. 26, 1595–1610 (2002). Publisher Full Text
Qi, J, Savinell, RF: Analysis of flow-through porous electrode cell with homogeneous chemical reactions: application to bromide oxidation in brine solutions. J. Appl. Electrochem.. 23, 873–886 (1993). Publisher Full Text
Aziz, A: A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition. Commun. Nonlinear Sci. Numer. Simul.. 14, 1064–1068 (2009). Publisher Full Text
Ishak, A: Similarity solutions for flow and heat transfer over a permeable surface with convective boundary condition. Appl. Math. Comput.. 217, 837–842 (2010). Publisher Full Text
Makinde, OD, Aziz, A: MHD mixed convection from a vertical plate embedded in a porous medium with a convective boundary condition. Int. J. Therm. Sci.. 49, 1813–1820 (2010). Publisher Full Text
Yao, S, Fang, T, Zhong, Y: Heat transfer of a generalized stretching/shrinking wall problem with convective boundary conditions. Commun. Nonlinear Sci. Numer. Simul.. 16, 752–760 (2011). Publisher Full Text
Masuda, H, Ebata, A, Teramae, K, Hishinuma, N: Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles. Netsu Bussei. 7, 227–233 (1993). Publisher Full Text
Tzou, DY: Thermal instability of nanofluids in natural convection. Int. J. Heat Mass Transf.. 51, 2967–2979 (2008). Publisher Full Text
Nield, DA, Kuznetsov, AV: Thermal instability in a porous medium layer saturated by a nanofluid. Int. J. Heat Mass Transf.. 52, 5796–5801 (2009). Publisher Full Text
Kuznetsov, AV, Nield, DA: Effect of local thermal non-equilibrium on the onset of convection in a porous medium layer saturated by a nanofluid. Transp. Porous Media. 83, 425–436 (2010). Publisher Full Text
Kuznetsov, AV, Nield, DA: Thermal instability in a porous medium layer saturated by a nanofluid: Brinkman model. Transp. Porous Media. 81, 409–422 (2010). Publisher Full Text
Kuznetsov, AV, Nield, DA: The onset of double-diffusive nanofluid convection in a layer of a saturated porous medium. Transp. Porous Media. 85, 941–951 (2010). Publisher Full Text
Nield, DA, Kuznetsov, AV: The onset of convection in a horizontal nanofluid layer of finite depth. Eur. J. Mech. B, Fluids. 29, 217–223 (2010). Publisher Full Text
Nield, DA, Kuznetsov, AV: The onset of convection in a layer of cellular porous material: effect of temperature-dependent conductivity arising from radiative transfer. J. Heat Transf.. 132, (2010) Article ID 074503
Nield, DA, Kuznetsov, AV: The onset of double-diffusive convection in a nanofluid layer. Int. J. Heat Fluid Flow. 32, 771–776 (2011). Publisher Full Text
Nield, DA, Kuznetsov, AV: The effect of vertical through flow on thermal instability in a porous medium layer saturated by a nanofluid. Transp. Porous Media. 87, 765–775 (2011). Publisher Full Text
Chand, R, Rana, GC: On the onset of thermal convection in rotating nanofluid layer saturating a Darcy-Brinkman porous medium. Int. J. Heat Mass Transf.. 55, 5417–5424 (2012). Publisher Full Text