TY - JOUR

T1 - Linear instability of a horizontal thermal boundary layer formed by vertical throughflow in a porous medium

T2 - the effect of local thermal nonequilibrium

AU - Patil, P.M.

AU - Rees, D.A.S.

PY - 2013/9

Y1 - 2013/9

N2 - In this paper we investigate the onset of convection in a saturated porous medium where uniform suction into a horizontal and uniformly hot bounding surface induces a stationary thermal boundary layer. Particular attention is paid to how the well-known linear stability characteristics of this boundary layer are modified by the presence of local thermal nonequilibrium effects. The basic conduction state is determined and it is found that the boundary layer forms two distinct regions when the porosity is small or when the conductivity of the fluid is small compared with that of the solid. A linearised stability analysis is performed which results in an ordinary differential eigenvalue problem for the critical Darcy-Rayleigh number as a function of the wave number and the two nondimensional parameters, H and γ, which are associated with local thermal nonequilibrium. This eigenvalue problem is solved numerically by first approximating the equations by fourth order compact finite differences, and then the critical Rayleigh number is computed iteratively using the inverse power method and minimised over the wavenumber. The variation of the critical Rayleigh number and wavenumber with H and γ is presented. One of the unusual effects of local thermal nonequilibrium is that there exists a parameter regime within which the neutral curve is bimodal.

AB - In this paper we investigate the onset of convection in a saturated porous medium where uniform suction into a horizontal and uniformly hot bounding surface induces a stationary thermal boundary layer. Particular attention is paid to how the well-known linear stability characteristics of this boundary layer are modified by the presence of local thermal nonequilibrium effects. The basic conduction state is determined and it is found that the boundary layer forms two distinct regions when the porosity is small or when the conductivity of the fluid is small compared with that of the solid. A linearised stability analysis is performed which results in an ordinary differential eigenvalue problem for the critical Darcy-Rayleigh number as a function of the wave number and the two nondimensional parameters, H and γ, which are associated with local thermal nonequilibrium. This eigenvalue problem is solved numerically by first approximating the equations by fourth order compact finite differences, and then the critical Rayleigh number is computed iteratively using the inverse power method and minimised over the wavenumber. The variation of the critical Rayleigh number and wavenumber with H and γ is presented. One of the unusual effects of local thermal nonequilibrium is that there exists a parameter regime within which the neutral curve is bimodal.

UR - http://www.scopus.com/inward/record.url?scp=84881061205&partnerID=8YFLogxK

UR - http://dx.doi.org/10.1007/s11242-013-0181-7

U2 - 10.1007/s11242-013-0181-7

DO - 10.1007/s11242-013-0181-7

M3 - Article

AN - SCOPUS:84881061205

SN - 0169-3913

VL - 99

SP - 207

EP - 227

JO - Transport in Porous Media

JF - Transport in Porous Media

IS - 2

ER -