Thermoconvective instabilities in a porous medium bounded by two concentric horizontal cylinders

Abstract

The study of natural convection in a saturated porous medium bounded by two concentric, horizontal, isothermal cylinders reveals different types of evolution according to the experimental conditions and the geometrical configuration of the model. At small Rayleigh numbers the state of the system corresponds to a regime of pseudo-conduction. The isotherms are coaxial with the cylinders. At larger Rayleigh numbers a regime of steady two-dimensional convection sets in between the two cylinders. Finally, for Rayleigh numbers above the critical Rayleigh number Ra*c the phenomena become three-dimensional and fluctuating. The appearance of these different regimes depends, moreover, on the geometry considered and, in particular, on two numbers: R, the ratio of the radii of the cylinders, and A, the ratio of the length of the cylinders to the radius of the inner one. In order to approach these experimental observations and to obtain realistic theoretical models, several methods of solving the equations have been used.The perturbation method yields information about the thermal field and the heat transfer between the cylinders under conditions close to the equilibrium state.A numerical two-dimensional model enables us to extend the range of investigation and to represent properly the phenomena when steady convection appreciably modifies the temperature distribution and the velocities within the porous layer.Neither of these models allows account to be taken of the instabilities observed experimentally above a critical Rayleigh number Ra*c. For this reason, a study of stability has been carried out using a Galerkin method based on equations corresponding to an initial state of steady convection. The results obtained show the importance of three-dimensional effects for the onset of fluctuating convection. The critical transition Rayleigh number Ra*c is thus determined in terms of the ratio of the radii R by solving an eigenvalue problem.A numerical three-dimensional model based on the method of finite elements has thus been developed in order to point out the different types of evolution with time. Steady two-dimensional convection and fluctuating three-dimensional convection have been actually found by calculation. The solution of the system of equations by the method of finite elements is briefly described.The experimental and theoretical results are then compared and a general physical interpretation is given.