Analysis of Transient Three-Dimensional Natural Convection in Porous Media

1981 ◽  
Vol 103 (2) ◽  
pp. 242-248 ◽  
Author(s):  
Y. T. Chan ◽  
S. Banerjee

Transient multidimensional natural convection in porous media has been studied using a numerical method based on the simplified marker and cell technique with local cancellation of low order, diffusional truncation errors. The conservation equations and boundary conditions were phrased in terms of the primitive variables, velocity and temperature. Differences in temperature between the fluid and the solid matrix are considered. Heat transfer between the solid and liquid phases was modelled by representing the porous medium as an assemblage of spherical particles, and solving the conduction problem within the spheres at every time step. Nusselt numbers at walls were calculated from the temperature and velocity profiles. Numerical results for heat transfer through fluid saturated porous media heated from below are in good agreement with published experiments. Consideration of heat transfer between the solid and fluid phase leads to Nusselt numbers that vary with the thermophysical properties of the solid material, even when the Rayleigh number and fluid thermophysical properties are kept constant. This is also observed in experiments. The calculations also show convective instabilities of the right period at high Rayleigh numbers.

1991 ◽  
Vol 113 (2) ◽  
pp. 391-399 ◽  
Author(s):  
E. David ◽  
G. Lauriat ◽  
P. Cheng

The problem of natural convection in differentially heated vertical cavities filled with spherical particles saturated with Newtonian fluids is investigated numerically. The Brinkman–Darcy–Ergun equation is used as the momentum equation, and the wall effect on porosity variation is approximated by an exponential function. The effect of variable stagnant thermal conductivities is taken into consideration in the energy equation. The formulation of the problem shows that the flow and heat transfer characteristics depend on six dimensionless parameters, namely, the Rayleigh and Prandtl numbers of the fluid phase, the dimensionless particle diameter, the conductivity ratio of the two phases, the bulk porosity, and the aspect ratio of the cavity. The influences of these parameters on the heat transfer rate are thoroughly investigated. The predicted Nusselt numbers are compared with existing experimental results. It is found that the computed Nusselt numbers based on the present model compare the best with experimental data.


1997 ◽  
Vol 08 (04) ◽  
pp. 879-888 ◽  
Author(s):  
R. G. M. Van Der Sman

A lattice-Boltzmann scheme for natural convection in porous media is developed and applied to the heat transfer problem of a 1000 kg potato packaging. The scheme has features new to the field of LB schemes. It is mapped on a orthorhombic lattice instead of the traditional cubic lattice. Furthermore the boundary conditions are formulated with a single paradigm based upon the particle fluxes. Our scheme is well able to reproduce (1) the analytical solutions of simple model problems and (2) the results from cooling down experiments with potato packagings.


1971 ◽  
Vol 11 (04) ◽  
pp. 399-405 ◽  
Author(s):  
Michel A. Combarnous ◽  
Paul Bia

Abstract An experimental study of thermal convection in a horizontal porous layer bounded by isothermal planes has been performed with and without a mean flow of the saturating fluid phase. The temperature distribution and heat transfer resulting from convection have been determined. The theoretical criterion for the onset of convection (Rayleigh number NRa > 4p2) has been verified. For low values of NRa (< 260) a regular pattern of convective cells has been observed which may be motionless or moving depending on the experimental conditions. For NRa values higher than 260, another convective state has been found that is mainly unstable. Numerical computations have been worked out that confirm the experimental results on the heat transfer and size of convective cells. INTRODUCTION Thermal convective currents may exist in a porous medium when the vertical component of the temperature gradient runs in the same direction as the gravity vector. Such a configuration is possible in some cases during a real in-situ combustion test4 or during hot fluid injection. A general understanding of convection is also of interest for analyzing the abnormal temperature gradients sometimes observed in oil and gas reservoirs.22 From the phenomenological standpoint a distinction must be made between the natural or free convection occurring in a closed volume and the mixed or combined free and forced convection arising when there is a mean flow of the fluid phase. Natural convection in porous media has been the subject of numerous experimental attempts to determine a convection criterion, e.g., by using a linear theory,1,2,9 and to measure the mean heat transfer.3,7 Numerical studies8,10,15,19,20 have also been made, and their results are not always in good agreement with experimental observations. However, no thorough study of mixed convection has been made except for a theoretical analysis of the convection onset criterion5 and some qualitative experimental observations.4 This paper presents some findings concerning both the natural and mixed convection of an incompressible fluid in a homogeneous horizontal porous layer bounded by two parallel impervious isothermal surfaces. The temperature for the cooler upper surface is T1 and that of the lower surface is T2=T1+?T. THEORETICAL ANALYSIS The effect of pressure on specific gravity is assumed to be negligible with regard to the effect of thermal expansion. So the relationship between density and temperature may be expressed as:Equation 1 In addition to complying with equations for motion and continuity, velocity and temperature distributions obey the heat transfer equation:Equation 2 By taking into account the additivity property of enthalpy, M may be expressed asEquation 3 On the contrary, khR is not just a simple function of khj and khs and must be measured or calculated with different kinds of models.6


1987 ◽  
Vol 109 (2) ◽  
pp. 356-362 ◽  
Author(s):  
J. T. Hong ◽  
Y. Yamada ◽  
C. L. Tien

This work examines analytically the effects of non-Darcian and nonuniform permeability conditions on the natural convection from a vertical plate in porous media. The non-Darcian effects, which include the no-slip and inertia effects, decrease the flow and heat transfer rate, while the nonhomogeneity effect enhances the heat transfer. For packed spheres, in particular, the nonhomogeneity in permeability due to the packing of spheres near the solid wall results in a strong flow-channeling effect that significantly increases the heat transfer. The effect of transverse thermal dispersion is also examined. This dispersion effect causes an increase in the heat transfer.


2010 ◽  
Vol 132 (6) ◽  
Author(s):  
Yuping Duan ◽  
S. F. Hosseinizadeh ◽  
J. M. Khodadadi

The effects of insulated and isothermal thin baffles on pseudosteady-state natural convection within spherical containers were studied computationally. The computations are based on an iterative, finite-volume numerical procedure using primitive dependent variables. Natural convection effect is modeled via the Boussinesq approximation. Parametric studies were performed for a Prandtl number of 0.7. For Rayleigh numbers of 104, 105, 106, and 107, baffles with three lengths positioned at five different locations were investigated (120 cases). The fluid that is heated adjacent to the sphere rises replacing the colder fluid, which sinks downward through the stratified stable thermal layer. For high Ra number cases, the hot fluid at the bottom of the sphere is also observed to rise along the symmetry axis and encounter the sinking colder fluid, thus causing oscillations in the temperature and flow fields. Due to flow obstruction (blockage or confinement) effect of baffles and also because of the extra heating afforded by the isothermal baffle, multi-cell recirculating vortices are observed. This additional heat is directly linked to creation of another recirculating vortex next to the baffle. In effect, hot fluid is directed into the center of the sphere disrupting thermal stratified layers. For the majority of the baffles investigated, the Nusselt numbers were generally lower than the reference cases with no baffle. The extent of heat transfer modification depends on Ra, length, and location of the extended surface. With an insulated baffle, the lowest amount of absorbed heat corresponds to a baffle positioned horizontally. Placing a baffle near the top of the sphere for high Ra number cases can lead to heat transfer enhancement that is linked to disturbance of the thermal boundary layer. With isothermal baffles, heat transfer enhancement is achieved for a baffle placed near the bottom of the sphere due to interaction of the counterclockwise rotating vortex and the stratified layer. For some high Ra cases, strong fluctuations of the flow and thermal fields indicating departure from the pseudosteady-state were observed.


1995 ◽  
Vol 117 (4) ◽  
pp. 910-917 ◽  
Author(s):  
T. J. Heindel ◽  
F. P. Incropera ◽  
S. Ramadhyani

Three-dimensional numerical predictions and experimental data have been obtained for natural convection from a 3 × 3 array of discrete heat sources flush-mounted on one vertical wall of a rectangular cavity and cooled by the opposing wall. Predictions performed in a companion paper (Heindel et al., 1995a) revealed that three-dimensional edge effects are significant and that, with increasing Rayleigh number, flow and heat transfer become more uniform across each heater face. The three-dimensional predictions are in excellent agreement with the data of this study, whereas a two-dimensional model of the experimental geometry underpredicts average heat transfer by as much as 20 percent. Experimental row-averaged Nusselt numbers are well correlated with a Rayleigh number exponent of 0.25 for RaLz ≲ 1.2 × 108.


10.2118/6813 ◽  
1968 ◽  
Author(s):  
K. Aziz ◽  
P.H. Holt ◽  
P.S. Karra

Sign in / Sign up

Export Citation Format

Share Document