An accurate numerical analysis for the onset of thermal convection in a two
-layer system is presented. The system comprises a saturated porous layer o
ver which lies a layer of the same fluid. The layered system is heated from
below, the upper (fluid) surface is free to the atmosphere, and convection
driven by surface tension is allowed for. The eigenvalues and eigenfunctio
ns for the instability problem are derived by utilizing a D-2 Chebyshev tau
method (J. J. Dongarra, B. Straughan, and D. W. Walker, 1996, Appl. Numer.
Math. 22, 399-435). This allows us to obtain highly accurate eigenvalues a
nd eigenfunctions in a very efficient manner. The onset of convection is se
en to have a bimodal nature in which convection may be dominated by the por
ous medium or by the fluid, depending on the depths of the relative layers
and the strength of the tension in the fluid surface. The effect of surface
tension is investigated in detail and it is found that for the parameter (
d) over cap (=depth of fluid layer/depth of porous layer) very small, the s
urface tension has a strong effect on convection dominated by the porous me
dium, whereas for (d) over cap larger the surface tension effect is observe
d only with the fluid made. (C) 2001 Academic Press.