SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Effect of variable viscosity on vortex instability of non-Darcy free convection boundary layer flow adjacent to a non-isothermal horizontal surface in a porous medium

Ahmed M Elaiw14*, Ahmed A Bakr24 and Fuad S Ibrahim35

Author Affiliations

1 Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia

2 Department of Mathematics, Faculty of Science for Girls, King Khalid University, Abha, Saudi Arabia

3 Department of Mathematics, University Collage in Makkah, Umm-Alqura University, Makkah, Saudi Arabia

4 Department of Mathematics, Faculty of Science, Al-Azhar University (Assiut Branch), Assiut, Egypt

5 Department of Mathematics, Faculty of Science, Assiut University, Assiut, Egypt

For all author emails, please log on.

Boundary Value Problems 2012, 2012:26  doi:10.1186/1687-2770-2012-26

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2012/1/26


Received:1 October 2011
Accepted:27 February 2012
Published:27 February 2012

© 2012 Elaiw et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this article, we study the effect of variable viscosity on the flow and vortex instability of non-Darcian free convection boundary layer flow on a horizontal surface in a saturated porous medium. The wall temperature is a power function of the distance from the origin. The variation of viscosity is expressed as an exponential function of temperature. The transformed boundary layer equations, which are developed using a non similar solution approach, are solved by means of a finite difference method. The analysis of the disturbance flow is based on linear stability theory. The local Nusselt number, critical Rayleigh number and the associated wave number at the onset of vortex instability are presented over a wide range of wall to ambient viscosity ratio parameters μ*= μw/μ. The variable viscosity effect is found to enhance the heat transfer rate and destabilize the flow for liquid heating, while the opposite trend is true for gas heating.

Keywords:
vortex instability; porous media; variable viscosity; non-Darcy; free convection; finite difference method

1 Introduction

The problems of the vortex mode of instability in free and mixed convection flows over horizontal and inclined heated surfaces in saturated porous media have received considerable attention (see [1-11]). The instability mechanism is due to the presence of a buoyancy force component in the direction normal to the surface. The importance of the problems is due to a large number of technical applications, such as material transfer associated with the storage of radioactive nuclear waste, separation processes in the chemical industry, filtration, transpiration cooling, transport processes in aquifers, and ground water pollution. Some researchers considered Darcy model [1-5], others considered non-Darcy model [6-11]. A comprehensive literature survey on this subject can be found in the recent book by Nield and Bejan [12].

All of the studies mentioned previously dealt with constant viscosity. The fundamental analysis of convection through porous media with temperature dependent viscosity is driven by several contemporary engineering applications from the cooling of electronic devices to porous journal bearings and is important for studying the variations in constitutive properties. The effect of variable viscosity for convective heat transfer through porous media has also been widely studied (e.g., [13-19]). The effect of variations in viscosity on the instability of flow and temperature fields are discussed by Kassoy and Zebib [13] and Gray et al. [14]. Lai and Kulacki [15] considered the variable viscosity effect for mixed convection along a vertical plate embedded in a saturated porous medium. The effect of variable viscosity on non-Darcy, free or mixed convection flow on a horizontal surface in a saturated porous medium is studied by Kumari [16]. Jayanthi and Kumari [17] studied the same problem presented in Kumari [16] for vertical surface. Afify [18] studied the effects of non-Darcy, variable viscosity and Hartmann-Darcy number on free convective heat transfer from an impermeable vertical plate embedded in a thermally stratified, fluid saturated porous medium for the case of power-law variation in the wall temperature. In [15-18], the viscosity of the fluid is assumed to vary as an inverse linear function of temperature, whereas in [19], the variation of viscosity with temperature is represented by an exponential function.

The effects of variable viscosity on the vortex instability of horizontal free and mixed convection boundary layer flows in a saturated porous medium for isothermal and non-isothermal surfaces are studied by Jang and Leu [19] and Elaiw et al. [20], respectively. However, the effect of variable viscosity on the flow and vortex instability of non-Darcian free convection boundary layer flow over a non-isothermal horizontal plate does not seem to have been investigated. This motivated the present investigation.

The purpose of this article is to examine in details the effect of temperature-dependent viscosity on the flow and vortex instability of a horizontal non-Darcian free convection boundary layer flow in a saturated porous medium with variable wall temperature. This is accomplished by considering the Ergun model equation of motion. The variation of viscosity with temperature is represented by an exponential function, which is more accurate than linear function especially for large temperature differences. The transformed boundary layer equations, which are given using a non similar solution approach, are solved by means of a finite difference method. The stability analysis is based on the linear stability theory. The resulting eigenvalue problem is solved using a finite difference scheme.

2 Analysis

2.1 The main flow

We consider a semi-infinite non-isothermal horizontal surface (Tw) embedded in a porous medium (T) as shown in Figure 1, where x represents the distance along the plate from its leading edge, and y the distance normal to the surface. The wall temperature is assumed to be a power function of x, i.e., Tw(x) = T+ Axλ, where A is a constant and λ is the parameter representing the variation of the wall temperature. In order to study transport through high porosity media, the original Darcy model is improved by including inertia. For the mathematical analysis of the problem we assume that, (i) local thermal equilibrium exists between the fluid and the solid phase, (ii) the physical properties are constant, except for the viscosity μ and the density ρ in the bouncy force, (iii) Ergun non-Darcy model [21] is considered, and (iv) the Boussinesq approximation is valid. With these assumptions, the governing equations are given by:

thumbnailFigure 1. The physical model and coordinate system.

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M1">View MathML</a>

(1)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M2">View MathML</a>

(2)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M3">View MathML</a>

(3)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M4">View MathML</a>

(4)

where ρ = ρ[1 - β(T - T)] is the fluid density. Note that the second term on the left-hand side of Equations (2) and (3) represents the inertia force, where K* is the inertia coefficient in Ergun model. As K* = 0, Equations (2) and (3) reduce to Darcy model.

The viscosity of the fluid μ is assumed to vary with temperature according to the following exponential form

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M5">View MathML</a>

(5)

where μis the absolute viscosity at ambient temperature and A1 is a constant adopted from the least square fitting for a particular fluid. Formula (5) is a generalization of that used in [19], where the wall temperature is taken to be constant.

The pressure terms appearing in Equations (2) and (3) can be eliminated through cross-differentiation. The boundary layer assumption yields ∂/∂x ≪ ∂/∂y and v u. With ψ being a stream function such that u = ∂ψ/y and v = -∂ψ/∂x, the Equations (1)-(4) become

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M6">View MathML</a>

(6)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M7">View MathML</a>

(7)

The boundary conditions are defined as follows

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M8">View MathML</a>

(8)

By introducing the following similarity variables:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M9">View MathML</a>

(9)

Equations (5)-(7) become

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M10">View MathML</a>

(10)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M11">View MathML</a>

(11)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M12">View MathML</a>

(12)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M13">View MathML</a>is the wall to ambient viscosity ratio parameter, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M14">View MathML</a> is the local Rayleigh number, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M15">View MathML</a> is the Rayleigh number based on the pore diameter and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M16">View MathML</a> is the Ergun number.

The boundary conditions are transformed as follows:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M17">View MathML</a>

(13)

where the primes denote the derivatives with respect to η. The velocity components and the heat transfer coefficient in term of the local Nusselt number are given by:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M18">View MathML</a>

(14)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M19">View MathML</a>

(15)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M20">View MathML</a>

(16)

2.2 The disturbance flow

The standard linear stability method yields the following:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M21">View MathML</a>

(17)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M22">View MathML</a>

(18)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M23">View MathML</a>

(19)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M24">View MathML</a>

(20)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M25">View MathML</a>

(21)

where the subscripts 0 and 1 signify the mean flow and disturbance components, respectively.

Following the order of magnitude analysis method described in detail by Hsu et al. [1], the terms ∂u1/∂x and ∂2T1/∂x2 in Equations (17) and (21) can be neglected. The omission of ∂u1/∂x in Equation (17) implies the existence of a disturbance stream function Ψ1, such as

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M26">View MathML</a>

(22)

Eliminating P1 from Equations (18)-(20) with the aid of (22) leads to

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M27">View MathML</a>

(23)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M28">View MathML</a>

(24)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M29">View MathML</a>

(25)

As in Hsu et al. [1], we assume that the three-dimensionless disturbances for neutral stability are of the form

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M30">View MathML</a>

(26)

where a is the spanwise periodic wave number. Substituting Equation (26) into Equations (23)-(25) yields

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M31">View MathML</a>

(27)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M32">View MathML</a>

(28)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M33">View MathML</a>

(29)

Equations (27)-(29) are solved based on the local similarity approximation [22], wherein the disturbances are assumed to have weak dependence in the streamwise direction (i.e., ∂/∂x ≪ ∂/∂η).

On applying the following transformations

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M34">View MathML</a>

(30)

into Equations (27)-(29) leads to the following:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M35">View MathML</a>

(31)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M36">View MathML</a>

(32)

with the boundary conditions

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M37">View MathML</a>

(33)

where G1 - G4 are given by

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M38">View MathML</a>

For fixed values of ξ, λ, k, Er, Rad, and μ*, the solution F and Θ is an eigenfunction for the eigenvalue Rax. When μ* = 1, Er = Rad = ξ = 0, Equations (31) and (32) reduce to those presented in [1] and when μ* = 1, ξ = 1, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M39">View MathML</a>, λ = 0.5 they reduce to those given in [8]. It may also remarked that for λ = Er = Rad = ξ = 0, Equations (31) and (32) reduce to those presented in Jang and Leu [19].

3 Numerical scheme

In this section, we compute the approximate values of Rax for Equations (31) and (32) with boundary conditions (33). An implicit finite difference method is used to solve first the base flow Equations (11) and (12) with boundary conditions (13). The results are stored for a fixed step size h, which is small enough to make an accurate linear interpolation between mesh point. The domain is 0 ≤ η ≤ η, where ηis the edge of the boundary layer of the basic flow. For a positive integer N, let h = η/N and ηi = ih, i = 0, 1, . . . , N. The problem is discretized with standard centered finite differences of order two, following Usmani [23]. The solution of the eigenvalue problem is achieved by using the subroutine GVLRG of the IMSL library, see [24].

4 Results and discussion

Numerical results for the local Nusselt number, neutral stability curves, critical Rayleigh numbers and associate wave numbers at the onset of vortex instability for wall to ambient viscosity ratio parameter μ* ranging from 0.1 to 10 and Ergun number Er ranging from 0.1 to 0.6 are presented. The effect of the non uniform temperature profile on the wall is also studied, and corresponds to variation of the parameter λ from 0.5 to 2. It is known that as the temperature is increased, the gas viscosity increases and the liquid viscosity decreases [19]. Therefore, for a heated wall, values of μ* > 1 corresponds to the case of gas heating, values of μ* < 1 corresponds to the case of liquid heating.

The local Nusselt number as a function ξ for various values of μ*, Er, and λ is shown in Figures 2 and 3. The heat transfer rate is observed to be more than the constant-viscosity case for liquids μ* < 1 and it has the opposite trend for gases μ* > 1. Further, a lower local Nusselt number is seen to occur with greater values of Er and ξ. Furthermore, a higher Nusselt number is seen to occur with greater values of λ.

thumbnailFigure 2. The local Nusselt number as a function of ξ for selected values of μ*.

thumbnailFigure 3. The local Nusselt number as a function of ξ for selected values of Er and λ.

Figures 4 and 5, respectively, show the neutral stability curves, in terms of Rayleigh numbers Rax and the dimensionless wave number k for selected values of μ* and ξ. Observe that as μ* and ξ increase, the neutral stability curves shift to a higher Rayleigh number Rax. On the other hand, for liquid heating μ* < 1, the neutral stability curves shift to lower Rayleigh numbers and higher wave numbers, indicating a destabilization of the flow, whereas for gas heating μ* > 1, the opposite trend is true.

thumbnailFigure 4. Neutral stability curves for selected values of μ*.

thumbnailFigure 5. Neutral stability curves for selected values of ξ.

The critical Rayleigh number and wave number are plotted as a function of Rad for selected values of μ* and λ in Figures 6 and 7. Observe that, as the variable viscosity parameter μ* or λ increases, the critical Rayleigh <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M40">View MathML</a> increases. Also, as Rad is increased the critical Rayleigh number <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M40">View MathML</a> is decreased. Further, the critical wave number k* is increased as Rad and λ are increased and μ* is decreased. From these figures it can be seen that, the inertia effects reduce the heat transfer rate and destabilize the flow. Finally, we conclude that, the variable viscosity effect enhances the heat transfer rate and destabilizes the flow for liquid heating, while the opposite tend is true for gas heating.

thumbnailFigure 6. Critical Rayleigh number as a function of Rad for selected values of μ* and λ.

thumbnailFigure 7. Critical wave number as a function of Rad for selected values of μ* and λ.

5 Conclusions

The non-Darcy free convection flow on a semi-infinite non-isothermal horizontal plate embedded in a porous medium with variable viscosity is investigated. The effects of variable viscosity characterized by the parameter μ* on the flow and vortex instability are examined. The governing partial differential equations are transformed to a non similar form by introducing appropriate transformations and are solved numerically using an implicit finite difference scheme. The resulting eigenvalue problem is solved by using a finite difference scheme. The effects of all involved parameters on the local Nusselt number, critical Rayleigh and associated wave number are presented. The results show that, for liquid heating, the variable viscosity effect enhances the heat transfer rate and destabilizes flow; for gas heating, the opposite is true.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All authors participated in the design of the study and performed the numerical simulation. All authors conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Nomenclature

a spanwise wave number

d mean particle diameter or pore diameter

f dimensionless base state stream function

F dimensionless disturbance stream function

g gravitational acceleration

i complex number

k dimensionless wave number

K permeability of porous medium

K* inertial coefficient in Ergun Equation

Nux local Nusselt number

P pressure

Rax local Rayleigh number

Rad Rayleigh number based on the pore diameter

T fluid temperature

u, v, w volume averaged velocity in the x, y, z directions

x, y, z axial, normal and spanwise coordinates

Greek symbols

α thermal diffusivity

β volumetric coefficient of thermal expansion

η pseudo-similarity variable

θ dimensionless base state temperature

Θ dimensionless disturbance temperature

μ dynamic viscosity of the fluid

ξ non similarity parameter

λ exponent in the wall temperature variation

ν kinematic viscosity

μ* wall to ambient viscosity ratio, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2012/1/26/mathml/M41">View MathML</a>

ψ stream function

Subscripts

w conditions at the wall

conditions at the free stream

0 basic undisturbed quantities

1 disturbed quantities

Superscripts

* critical value

' differentiation w.r.t η

Acknowledgements

The authors were grateful to the anonymous reviewers for constructive suggestions and valuable comments, which had improved the quality of the article.

References

  1. Hsu, CT, Cheng, P, Homsy, GM: Instability of free convection flow over a horizontal impermeable surface in a porous medium. Int J Heat Mass Transfer. 21, 1221–1228 (1978). Publisher Full Text OpenURL

  2. Jang, JY, Chang, WJ: Vortex instability in buoyancy-induced inclined boundary layer flow in a saturated porous medium. Int J Heat Mass Transfer. 31, 759–767 (1988). Publisher Full Text OpenURL

  3. Elaiw, AM: Vortex instability of mixed convection boundary layer flow adjacent to a non-isothermal inclined surface in a porous medium with variable permeability. ZAMM Z Angew Math Mech. 88, 121–128 (2008). Publisher Full Text OpenURL

  4. Hassanien, IA, Salama, AA, Elaiw, AM: Variable permeability effect on vortex instability of mixed convection flow in a semi-infinite porous medium bounded by a horizontal surface. Appl Math Comput. 146, 829–847 (2003). Publisher Full Text OpenURL

  5. Hassanien, IA, Salama, AA, Elaiw, AM: Variable permeability effect on vortex instability of a horizontal natural convection flow in a saturated porous medium with variable wall temperature. ZAMM Z Angew Math Mech. 84, 39–47 (2004). Publisher Full Text OpenURL

  6. Elaiw, AM, Ibrahim, FS, Bakr, AA: The influence of variable permeability on vortex instability of a horizontal combined free and mixed convection flow in a saturated porous medium. ZAMM Z Angew Math Mech. 87, 528–536 (2007). Publisher Full Text OpenURL

  7. Chang, WJ, Jang, JY: Non-Darcian effects on vortex instability of a horizontal natural convection flow in a saturated porous medium. Int J Heat Mass Transfer. 32, 529–539 (1989). Publisher Full Text OpenURL

  8. Chang, WJ, Jang, JY: Inertia effects on vortex instability of a horizontal natural convection flow in a saturated porous medium. Int J Heat Mass Transfer. 32, 541–550 (1989). Publisher Full Text OpenURL

  9. Hassanien, IA, Salama, AA, Moursy, NM: Non-Darcian effects on vortex instability of mixed convection over horizontal plates in porous medium with surface mass flux. Int Commun Heat Mass Transfer. 31, 231–240 (2004). Publisher Full Text OpenURL

  10. Hassanien, IA, Salama, AA, Moursy, NM: Inertia effect on vortex instability of horizontal natural convection flow in a saturated porous medium with surface mass flux. Int Comm Heat Mass Transfer. 31, 741–750 (2004). Publisher Full Text OpenURL

  11. Lee, DH, Yoon, DY, Choi, CK: The onset of vortex instability in laminar convection flow over an inclined plate embedded in a porous medium. Int J Heat Mass Transfer. 43, 2895–2908 (2000). Publisher Full Text OpenURL

  12. Nield, DA, Bejan, A: Convection in Porous Media. Springer Science + Business Media, New York, NY, USA (2006)

  13. Kassoy, DR, Zebib, A: Variable viscosity effects on the onset of convection in porous media. Phys Fluids. 18, 1649–1651 (1975). Publisher Full Text OpenURL

  14. Gray, J, Kassoy, DR, Tadjeran, H, Zebib, A: The effect of significant viscosity variation on convective heat transport in water saturated porous media. J Fluid Mech. 117, 233–249 (1982). Publisher Full Text OpenURL

  15. Lai, FC, Kulacki, FA: The effect of variable viscosity on convective heat transfer along a vertical surface in a saturated porous medium. Int J Heat Mass Transfer. 33, 1028–1031 (1990). Publisher Full Text OpenURL

  16. Kumari, M: Effect of variable viscosity on non-Darcy free or mixed convection flow on a horizontal surface in a saturated porous medium. Int Commun Heat Mass Transfer. 28, 723–732 (2001). Publisher Full Text OpenURL

  17. Jayanthi, S, Kumari, M: Effect of variable viscosity on non-Darcy free or mixed convection flow on a vertical surface in a fluid saturated porous medium. Mech Res Commun. 33, 148–156 (2006). Publisher Full Text OpenURL

  18. Afify, AA: Effects of variable viscosity on non-Darcy MHD free convection along a non-isothermal vertical surface in a thermally stratified porous medium. Appl Math Model. 31, 1621–1634 (2007). Publisher Full Text OpenURL

  19. Jang, JY, Leu, JS: Variable viscosity effect of on the vortex instability of free convection boundary layer flow over a horizontal surface in porous medium. Int J Heat Mass Transfer. 36, 1287–1294 (1993). Publisher Full Text OpenURL

  20. Elaiw, AM, Ibrahim, FS, Bakr, AA, Salama, AA: Effect of variable viscosity on vortex instability of mixed convection boundary layer flow adjacent to a non-isothermal horizontal surface in a porous medium. Arabian J Sci Eng. 36, 1517–1528 (2011). Publisher Full Text OpenURL

  21. Ergun, S: Fluid flow through packed columns. Chem Eng Prog. 48, 89–94 (1952)

  22. Sparrow, EM, Quack, H, Boerner, CT: Local nonsimilarity boundary layer solutions. AIAA J. 8, 1936–1942 (1970). Publisher Full Text OpenURL

  23. Usmani, RA: Some new finite difference methods for computing eigenvalues of two-point boundary value problems. Comput Math Appl. 9, 903–909 (1985)

  24. IMSL: References Manual. IMSL, Inc., Houston, TX (1990)