Publication Cover
Numerical Heat Transfer, Part A: Applications
An International Journal of Computation and Methodology
Volume 69, 2016 - Issue 9
193
Views
10
CrossRef citations to date
0
Altmetric
Original Articles

Buoyancy-induced flow and heat transfer in a porous annulus between concentric horizontal circular and square cylinders

, , , &
Pages 1029-1050 | Received 18 Mar 2015, Accepted 26 Jun 2015, Published online: 01 Feb 2016
 

ABSTRACT

This paper reports on natural convection heat transfer in a porous annulus between concentric horizontal circular and square cylinders. The heated inner circular cylinder is maintained at the uniform hot temperature Th, whereas the cooled outer square duct is held at the uniform cold temperature Tc. A pressure-based collocated finite-volume method is used to numerically investigate the effects on the total heat transfer of Rayleigh number (Ra), Prandtl number (Pr), Darcy number (Da), porosity (ϵ), and annulus aspect ratio (R/L). Results demonstrate that at low Ra values, conduction is the dominant heat transfer mode. Convection contribution to total heat transfer becomes more important beyond a critical Ra value, which decreases with an increase in Da and/or ϵ. Furthermore, an increase in the enclosure aspect ratio (R/L) leads to an increase in total heat transfer. A similar behavior is obtained with Prandtl number, where predictions indicate higher heat transfer rates at higher Pr values with its effect increasing as Ra increases. Streamlines and isotherms reveal flow separation for some of the reported cases. Limited computations are also performed for natural convection in a porous annulus between two horizontal concentric circular cylinders having the same inner and outer perimeters as the investigated enclosure. Comparison of the predicted average Nusselt number estimates with similar ones obtained in the original enclosure reveals a large percentage difference in values, demonstrating the strong influence of geometry on natural convection in enclosures.

Nomenclature

cp=

specific heat of fluid at constant pressure

dp=

pore diameter

Da=

Darcy number

F=

constant in Forchheimer’s extension

g=

gravitational acceleration

j=

unit vector along the y-axis

k=

thermal conductivity of fluid

K=

permeability of the porous media

L=

width of the square duct

n=

unit vector normal to the surface

Nu=

local Nusselt number

=

average Nusselt number

p=

dimensional pressure

P=

dimensionless pressure

Pr=

Prandtl number

R=

radius of the cylinder

Ra=

Rayleigh number

Si, So=

distance along the inner and outer enclosure surfaces

T=

dimensional temperature

u, U=

dimensional and dimensionless x-velocity component

v, V=

dimensional and dimensionless y-velocity component

v, V=

dimensional and dimensionless velocity vector

x, y=

dimensional coordinates

X, Y=

dimensionless coordinates

α=

thermal diffusivity

β=

thermal expansion coefficient

μ=

dynamic viscosity

ν=

kinematic viscosity

ρ=

density

θ=

dimensionless temperature

ϵ=

porosity

ψ=

stream function

Subscripts=
c=

cold wall

h=

hot wall

i=

condition at inner surface

o=

condition at outer surface

Nomenclature

cp=

specific heat of fluid at constant pressure

dp=

pore diameter

Da=

Darcy number

F=

constant in Forchheimer’s extension

g=

gravitational acceleration

j=

unit vector along the y-axis

k=

thermal conductivity of fluid

K=

permeability of the porous media

L=

width of the square duct

n=

unit vector normal to the surface

Nu=

local Nusselt number

=

average Nusselt number

p=

dimensional pressure

P=

dimensionless pressure

Pr=

Prandtl number

R=

radius of the cylinder

Ra=

Rayleigh number

Si, So=

distance along the inner and outer enclosure surfaces

T=

dimensional temperature

u, U=

dimensional and dimensionless x-velocity component

v, V=

dimensional and dimensionless y-velocity component

v, V=

dimensional and dimensionless velocity vector

x, y=

dimensional coordinates

X, Y=

dimensionless coordinates

α=

thermal diffusivity

β=

thermal expansion coefficient

μ=

dynamic viscosity

ν=

kinematic viscosity

ρ=

density

θ=

dimensionless temperature

ϵ=

porosity

ψ=

stream function

Subscripts=
c=

cold wall

h=

hot wall

i=

condition at inner surface

o=

condition at outer surface

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.