Mixed convection in porous media between vertical concentric cylinders

In this paper, the dimensionless parameters Gr^**/Re_p and H/d_p that govern mixed convection in a vertical annulus filled with porous media, are deduced through dimensional analysis, and the criteria for mixed convection is established thereby. The calculated results for detecting the non-Darcian effects are discussed. The predictions agree with the experiments. © 1999 Scripta Technica, Heat Trans Asian Res, 28(2): 95–101, 1999     

Free convection of composite porous medium in a vertical channel

A fully developed free convection flow of immiscible fluids in a vertical channel filled with a porous medium is analyzed in the presence of source/sink. The flow is modeled using the Darcy–Brinkman–Forchheimer equation model. The viscous and Darcy dissipation terms are included in the energy equation. The channel walls are maintained at two different constant temperatures. The transport properties of both fluids are assumed to be constant. Continuous conditions for velocity, temperature, shear stress, and heat flux of both fluids at the interface are employed. The resulting coupled nonlinear equations are solved analytically using regular perturbation method and numerically using finite difference method. The velocity and temperature profiles are obtained in terms of porous parameter, Grashof number, viscosity ratio, width ratio, conductivity ratio, and heat generation or heat absorption coefficient. It is found that the presence of porous matrix and heat absorption reduces the flow field. © 2011 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.20340     

Natural convection and its fractal for liquid freezing in a vertical cavity filled with porous medium

The numerical simulation for a freezing liquid-saturated porous media in a vertical cylindrical cavity under the third kind of thermal boundary condition is reported in this paper. It shows that the effect of natural convection in the liquid phase decreases the freezing layer thickness and the freezing front has a wave shape instead of a stable plane, with one or more pair of eddy cells. This indicates a fractal existence. © 1999 Scripta Technica, Heat Trans Asian Res, 28(3): 165–171, 1999     

Experimental study of natural convection heat transfer of air layers in vertical annuli under high rayleigh number conditions

In this paper experimental investigations of natural convection heat transfer of air layers in vertical annuli are presented. In these experiments, the surface of the inner cylinder is maintained at a constant heat flux condition and the outer cylinder is cooled in the atmosphere. In order to obtain the convective contribution, the overall heat transfer data are corrected for thermal radiation and axial conduction losses from the end plates in the annuli. Special emphasis, in these investigations, was placed on the high Rayleigh number region where no experimental data are available in the literature. Data were obtained for Rayleigh numbers greater than 10⁹. The radius ratios studied were 2.03 and 3.92, and the aspect ratios studied were 23.94 and 66.67. Heat transfer correlations for average Nusselt numbers were developed for different Rayleigh number regions. For the low Rayleigh number region the results of this paper agree with the correlations reported in the literature. Much needed data and correlations for the high Rayleigh number region are obtained for the first time. These results improve the predictive ability for the heat transfer characteristics in the high Rayleigh number region. ©1999 Scripta Technica, Heat Trans Asian Res, 28(1): 50–57, 1999     

Mathematical modeling of onset of convection in a porous layer with viscosity variation

IntroduthenIn a novel EOR (Enhanced Oil Recovery) method,diffusion-convechon may be employed as a drive toforce heavy, viscous oil to move. A ~ of gas isexpanded beneath oil, and gas is allowed tO diffuseupward. As oil dissolves gas, it becomes lighter.Therefore, an inverse density gradient is develOPed in oillayer gradually. When this density gradient exceeds acritical value, convechon is started. It may be concjudedthat develOPing a theory, which predictS the onset ofconvection, should…     

Natural convection along a vertical complex wavy surface

A previously proposed transformation has been applied to the natural-convection boundary layer along a complex vertical surface created from two sinusoidal functions, a fundamental wave and its first harmonic. The numerical results demonstrate that the additional harmonic substantially alters the flow field and temperature distribution near the surface. The conclusion that the averaged heat-transfer rate per unit-wetted wavy surface is less than that of a corresponding flat plate [L.S. Yao, Natural convection along a vertical wavy surface, ASME J. Heat Transfer 105 (1983) 465–468] has been confirmed for this more complex surface. On the other hand, the total heat-transfer rates for a complex surface are greater than that of a flat plate. The numerical results show that the enhanced total heat-transfer rate seems to depend on the ratio of amplitude and wavelength of a surface.     

Conjugate natural convection in a fluid-saturated porous enclosure with two solid vertical partitions

Conjugate natural convection in a fluid-saturated square porous enclosure with two solid vertical partitions of finite and equal thickness equispaced from center of enclosure is investigated in this paper. The primary objective is to attenuate the Nusselt number (Nu) and hence the heat transfer rate across a differentially heated enclosure. Darcy's model is considered. Numerical computation is performed using successive accelerated replacement and explicit scheme. Partition ratio, partition length, thermal conductivity ratio, and modified Rayleigh number are the parameters under study. Fluid flow is analyzed by observing transient changes of streamlines and isotherms for partition length 0.3-1, thermal conductivity ratio 0.5-2, partition ratio 0.1-0.3 and modified Rayleigh number 100 and 1000 where partition ratio is the ratio of distance between center of enclosure and either of the partition center to the total length of the enclosure; while Nusselt number is calculated to estimate the heat transfer rate for each configuration. It is found that, employing a solid partition within the enclosure most definitely reduces the Nusselt number. The drop in Nusselt number is more for partition length 0-0.6 after which it does show a drop in Nu but only very subtle. Further, Nu is the least for partition ratio 0.2. Also, Nusselt number is proportional to thermal conductivity ratio which is the ratio of thermal conductivity of solid to porous medium.     

Nonsimilar solutions for mixed convection of water at 4°C over a vertical surface with a convection boundary condition in a porous medium

A boundary layer analysis has been presented for the mixed convection of water at 4°C over a vertical plate embedded in a porous medium. The Robin or convective boundary condition at the surface has been considered where the heat lost from the surface is the product of a heat transfer coefficient and the temperature difference between the surface and the free stream. The governing non‐similar boundary layer equations for both the forced and free convection dominated regimes were solved numerically by means of an implicit finite difference method. The friction factor and dimensionless heat transfer rate (Nusselt number) are presented for several values of the dimensionless heat transfer coefficient and buoyancy parameter. © 2012 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library (wileyonlinelibrary.com/journal/htj). DOI 10.1002/htj.21022     

Effects of buoyancy assisting and opposing flows on mixed convection boundary layer flow over a permeable vertical surface

The present study deals with new similarity solution of steady mixed convection boundary layer flow over a permeable surface for convective boundary condition. It has been shown that a self similar solution is possible when the mass transfer velocity at the surface of the plate varies like x^−1/2, where x is the distance from the leading edge of the solid surface. Two point boundary value problem governed by non-linear coupled ordinary differential equations have been solved numerically using implicit finite difference scheme in combination with the quasi-linearization technique. It is interesting to note that dual solutions exist for buoyancy assisting flow, besides that usually reported in literature for buoyancy opposing flow. Further, the buoyancy assisting force causes considerable overshoot in the velocity profile and the Prandtl number strongly affects the thermal boundary layer thickness including the surface heat transfer rate.     

Transient free convection in cold water past an infinite vertical porous plate

A finite-difference solution of the transient free convection flow of water near 4°C past an infinite vertical porous plate has been carried out. The most appropriate relation for the density in terms of the temperature, as derived by Gebhart and Mollendorf, is employed. It has been observed that, when the buoyancy force is in the upward direction, the velocity profiles are positive and opposite to those in the case where the buoyancy force is in the downward direction. Also, the skin friction is more when the buoyancy force is in the upward direction than when the buoyancy force is in the downward direction.     

Natural convection heat transfer of high temperature gas in an annulus between two vertical concentric cylinders

Water cooling panels have been adopted as the vessel cooling system of the High Temperature Engineering Test Reactor (HTTR) to cool the reactor core indirectly by natural convection and thermal radiation. In order to investigate the heat transfer characteristics of high temperature gas in a vertical annular space between the reactor pressure vessel and cooling panels of the HTTR, we carried out experiments and numerical analyses on natural convection heat transfer coupled with thermal radiation heat transfer in an annulus between two vertical concentric cylinders with the inner cylinder heated and the outer cylinder cooled. In the present experiments, Rayleigh number based on the height of the annulus ranged from 2.0 × 10⁷ to 5.4 × 10⁷ for helium gas and from 1.2 × 10⁹ to 3.5 × 10⁹ for nitrogen gas. The numerical results were in good agreement with the experimental ones regarding the surface temperatures of the heating and cooling walls. As a result of the experiments and the numerical analyses, the heat transfer coefficient of natural convection coupled with thermal radiation was obtained as functions of Rayleigh number, radius ratio, and the temperatures and emissivities of the heating and cooling wall surfaces. © 2005 Wiley Periodicals, Inc. Heat Trans Asian Res, 34(5): 293–308, 2005; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.20070     

Double diffusive natural convection along an inclined wavy surface in a porous medium

This paper studies the double diffusive natural convection near an inclined wavy surface in a fluid saturated porous medium with constant wall temperature and concentration. A coordinate transformation is employed to transform the complex wavy surface to a smooth surface, and the obtained boundary layer equations are solved by the cubic spline collocation method. Effects of angle of inclination, Lewis number, buoyancy ratio, and wavy geometry on the heat and mass transfer characteristics are studied. Results show that a decrease in the angle of inclination leads to a greater fluctuation of the local Nusselt and Sherwood numbers. Moreover, increasing the angle of inclination tends to increase the total heat and mass transfer rates.     

The onset of transient convection in bottom heated porous media

The theory of transient convection in bottom heated porous media under constant heat flux (CHF) condition or fixed surface temperature (FST) condition is advanced and verified by computational fluid dynamics (CFD) simulations. The use of κ∗, instead of κ_m tends to artificially inflate the value of Rayleigh number by about 30%. A new transient Rayleigh number for unsteady-state heat conduction was defined to predict the onset of transient convection in porous media, which were successfully simulated. The critical transient Rayleigh number from the simulation for CHF was about 29.60, which is close to the theoretical value of 27.1 calculated by Ribando and Torrance in 1976. In the case of FST, the critical transient Ra_c was found to be 30.9, which is close to the theoretical value of 32.3. The critical times of onset for simulations were predicted with good accuracy. The prediction of the critical wavelengths of the emerging plumes were fair for the 2D simulations. Any experiment to verify the linear stability analysis for thermal instability must simultaneously concur in the three eigenvalue parameters, namely the Biot number, the critical wavenumber and the corresponding critical Rayleigh number, apart from the physical boundaries. The average maximum transient Nusselt number was found to be 3.41 for CHF and 3.5 for FST respectively.     

Transient free convective flow through a vertical cylinder filled with a porous material

This article aims to explore the impressive impact of emerging parameters on transient fully evolved free convective flow inside a vertical cylinder containing a porous material. The mathematical formulation of the model related to the considered physical circumstance is presented under compatible boundary conditions. Closed‐form solutions are received for the velocity field, the temperature distribution, mass flux, skin friction, and the Nusselt number in terms of Bessel functions and modified Bessel functions of the first kind. Impressive effects of parameters such as the Darcy number $Da$, Prandtl number $\text{Pr}$, viscosity ratio $M$, and also time $t$ on both the velocity and temperature distribution have been explored employing graphs and tables. It is irradiated by analysis that flow erection, heat transfer rate, skin friction, and mass flux are admirably impacted by the Prandtl number, the Darcy number, viscosity ratio parameter, and time. It is found that both the velocity and temperature field profiles rise with the rising value of time and ultimately attain their steady state. Moreover, the Prandtl number and the viscosity ratio parameter reduce the velocity profiles, while the reverse phenomenon occurs with the Darcy number.     

Natural convection heat and mass transfer along a vertical wavy surface

A numerical study of natural convection heat and mass transfer along a vertical wavy surface has been performed. The wavy surface is maintained at uniform wall temperature and constant wall concentration. A simple coordinate transformation is employed to transform the complex wavy surface to a flat plate. A marching finite-difference scheme is used for the analysis. The buoyancy ratio N, amplitude-wavelength ratio α, and Schmidt number Sc are important parameters for this problem. The numerical results, including the developments of skin-friction coefficient, velocity, temperature, concentration, Nusselt number as well as Sherwood number along the wavy surfaces are presented. The effects of the buoyancy ratio N and the dimensionless amplitude of wavy surface on the local Nusselt number and the local Sherwood number have been examined in detail.     

Effect of variable viscosity on mixed convection boundary layer flow over a vertical surface embedded in a porous medium

The steady mixed convection boundary layer flow over a vertical impermeable surface embedded in a porous medium when the viscosity of the fluid varies inversely as a linear function of the temperature is studied. Both cases of assisting and opposing flows are considered. The transformed boundary layer equations are solved numerically by a finite difference method. Numerical results for the flow and heat transfer characteristics are obtained for various values of the mixed convection parameter ε and the variable viscosity parameter θ_e. It has been found that in the opposing flow case, dual solutions exist and boundary separation occurs.     

Developing buoyant convection in vertical porous annuli with unheated entry and exit

The open‐ended vertical double‐passage annular space between three vertical concentric coaxial cylinders is an important physical configuration portraying many practical applications. Hence, in the present analysis, the developing buoyant convection in vertical double‐passage annuli filled with fluid‐saturated porous media is studied numerically by imposing unheated entry and unheated exit thermal boundary conditions. The numerical solutions of the mathematical model equations are found through finite difference technique. The velocity profiles in radial as well as axial directions and temperature profiles have been depicted for vast range of nondimensional numbers, baffle position, and heating and un‐heating ratio. The velocity and thermal gradients decreases as heating section length decreases. Maximum velocity and heat transport occurs in a narrow annular passage rather than equal or wider passages. The presence of porosity causes a reduction in flow velocities and thermal gradients.     

MHD free convection flow along a vertical wavy surface with heat generation or absorption effect

We formulate the problem of free convection from a vertical wavy surface embedded in a uniform porous medium in the presence of an external magnetic field and internal heat generation or absorption effects. Using the appropriate transformations, the boundary layer equations are reduced to non-linear partial differential equations. The transformed boundary layer equations are solved numerically using Runge–Kutta integration scheme with the shooting technique. We have focused our attention on the evaluation of the local Nusselt number Nu_x, dimensionless velocity, f′, and temperature, θ. The governing parameters are the amplitude of the waviness of the surface, a, ranging from 0.0 (flat plate) to 0.3, and the heat generation absorption parameter Q ranging from − 0.25 to 0.25, and magnetic parameter Mn, ranging from 0.0 to 2.0. The effect of all these parameters are discussed and plotted.     

Hall and ion-slip effects on MHD free convection flow of rotating Jeffrey fluid over an infinite vertical porous surface

An attempt has been made to explore Hall and ion-slip effects on an unsteady magnetohydrodynamic rotating flow of an electrically conducting, viscous, incompressible, and optically thick radiating Jeffrey fluid past an impulsively vertical moving porous plate. Analytical solutions of the governing equations are obtained by Laplace transform technique. The analytical expressions for skin friction, Nusselt number, and Sherwood number are also evaluated. The velocity, temperature, and concentration distributions are displayed graphically in detail. From engineering point of view, the changes in skin friction, Nusselt number, and Sherwood number are observed with the computational results presented in a tabular manner. It is observed that the effects of rotation and Hall current tend to accelerate secondary velocity and decelerate primary velocity throughout the boundary layer region. Thermal and concentration buoyancy forces tend to accelerate both velocity components. Thermal radiation and thermal diffusion tend to enhance fluid temperature throughout the boundary layer region. Rotation and Jeffrey fluid parameters tend to enhance both stress components.     

Natural convection and radiation heat transfer in a vertical porous layer with a hexagonal honeycomb core (Part 1: numerical analysis)

Combined heat transfer characteristics were obtained numerically for three-dimensional natural convection and thermal radiation in a long and wide vertical porous layer with a hexagonal honeycomb core. We assumed that the porous layer was both homogeneous and isotropic. The pure Darcy law for fluid flow and Rosseland's approximation for radiation were used. The numerical methodology was based on an algebraic coordinate transformation technique and the transformed governing equations were solved using the SIMPLE algorithm. The effect of radiation on the heat transfer characteristics was investigated over a wide range of radiation numbers and temperature ratios for two Darcy-Rayleigh number values (Ra^* = 100 and 1000) and for a fixed aspect ratio of H/L = 1. The results are presented in the form of combined radiation and convection heat transfer coefficients and are compared with the corresponding values for pure natural convection. © 1999 Scripta Technica, Heat Trans Asian Res, 28(4): 278–294, 1999