Thermal and fingering convection in superposed fluid and porous layers.

Abstract

Thermal and fingering convection in a horizontal porous layer underlying a fluid layer was studied using linear stability analysis, experiment (for the thermal convection case only), and nonlinear simulation. For the thermal convection case, the linear analysis shows that when the fluid layer is thin, convection is largely confined to the porous layer. When the fluid layer thickness exceeds 15% of the porous layer thickness, convection is localized in the fluid layer and the critical wavelength is dramatically reduced. Experimental investigations were then conducted in a test box 24 cm x 12 cm x 4 cm high to substantiate the predictions. The ratio of the thickness of the fluid layer to that of the porous layer, d, varied from 0 to 1. The results were in good agreement with predictions. To investigate supercritical convection, a nonlinear computational study was carried out. It was found that for d ≤ 0.13, the Nusselt number increases sharply with the thermal Rayleigh number, whereas at larger values of d, the increase is more moderate. Heat transfer rates predicted for d = 0.1 and 0.2 are in good agreement with the experimental results. For salt-finger convection at R(m) ≤ 1, the critical value of the solute Rayleigh number R(sm) decreases as d increases; the convection is unicellular. For 5 ≤ R(m) ≤ 10, the critical R(sm) initially decreases with d, and then remains almost constant for larger values of d; multicellular convection prevails at high d. For 20 ≤ R(m) ≤ 50, the critical R(sm) first decreases and then increases as d increases from 0 to 0.1. When d > 0.1, the critical R(sm) decreases slowly with d and remains almost constant for d ≥ 0.4. In the nonlinear computations for R(m) = 1, periodic convection sets in at a value of R(sm) between ten and eleven times the critical value. For the case of R(m) = 50, an aperiodic oscillation occurs when R(sm) is between four and five times the critical value. For the superposed layer cases d = 1 and 0.5, the convection characteristics are similar to those of thermal convection when R(m) = 0.01. For R(m) = 1, it was found that the onset of salt-finger convection is oscillatory. For R(m) = 50, the nonlinear code failed to obtain satisfactory results.