Publication Cover
Numerical Heat Transfer, Part A: Applications
An International Journal of Computation and Methodology
Volume 72, 2017 - Issue 2
222
Views
2
CrossRef citations to date
0
Altmetric
Original Articles

Buoyancy-induced flow, heat, and mass transfer in a porous annulus

, , , &
Pages 107-125 | Received 27 Apr 2017, Accepted 07 Jul 2017, Published online: 22 Aug 2017
 

ABSTRACT

This paper reports on double-diffusive natural convection heat transfer in a porous annulus between concentric horizontal circular and square cylinders. A pressure-based segregated finite volume method is used to solve the problem numerically. The diffusion fluxes are discretized using the MIND fully implicit scheme. Furthermore, a modified pressure correction equation is derived that implicitly accounts for the nonorthogonal diffusion terms, which are usually neglected in the standard SIMPLE algorithm. Results indicate that convection effects increase with an increase in Rayleigh number, Darcy number, porosity, and enclosure aspect ratio. Further, at low Darcy values, porosity has no effect on the flow, temperature, and concentration fields.

Nomenclature

C=

main grid point at an element centroid

cp=

specific heat of fluid at constant pressure

dCF=

vector joining the two points C and F

dCF=

magnitude of dCF

dp=

pores diameter

Da=

Darcy number

eCF=

unit vector in the direction of dCF

E=

distance vector in the direction of dCF

E=

magnitude of E

F=

neighbor of element C; also constant in Forchheimer’s extension

g=

gravitational acceleration

h=

local convection heat transfer coefficient

=

average convection heat transfer coefficient

j=

unit vector along the y-axis

k=

thermal conductivity of fluid

K=

permeability of the porous media

L=

width of square duct

Le=

Lewis number

n=

unit vector normal to surface

N=

buoyancy ratio: Ras/Rat

NC, NF=

locations used in the calculation of the nonorthogonal part of the diffusive flux

Nu, =

local and average Nusselt number

P, P=

dimensional and dimensionless pressure

Pr=

Prandtl number

R=

radius of cylinder

Ra=

Rayleigh number

R/L=

enclosure aspect ratio

S=

dimensional solute concentration; also magnitude of S

S=

surface vector

Si, So=

distance along the inner and outer enclosure surfaces

Sh, =

local and average Sherwood number

T=

dimensional temperature

T=

vector equal to SE

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

ρ=

density

θ=

dimensionless temperature

ε=

porosity

ψ=

stream function

σ=

dimensionless solute concentration

Γ=

diffusion coefficient

Subscripts=
c=

cold wall or convection heat transfer

C=

refers to main grid point

f=

refers to element face

F=

refers to the F grid point

h=

hot wall

i=

condition at inner surface

m=

mass transfer

o=

condition at outer surface

s=

refers to solute concentration

t=

refers to temperature

Nomenclature

C=

main grid point at an element centroid

cp=

specific heat of fluid at constant pressure

dCF=

vector joining the two points C and F

dCF=

magnitude of dCF

dp=

pores diameter

Da=

Darcy number

eCF=

unit vector in the direction of dCF

E=

distance vector in the direction of dCF

E=

magnitude of E

F=

neighbor of element C; also constant in Forchheimer’s extension

g=

gravitational acceleration

h=

local convection heat transfer coefficient

=

average convection heat transfer coefficient

j=

unit vector along the y-axis

k=

thermal conductivity of fluid

K=

permeability of the porous media

L=

width of square duct

Le=

Lewis number

n=

unit vector normal to surface

N=

buoyancy ratio: Ras/Rat

NC, NF=

locations used in the calculation of the nonorthogonal part of the diffusive flux

Nu, =

local and average Nusselt number

P, P=

dimensional and dimensionless pressure

Pr=

Prandtl number

R=

radius of cylinder

Ra=

Rayleigh number

R/L=

enclosure aspect ratio

S=

dimensional solute concentration; also magnitude of S

S=

surface vector

Si, So=

distance along the inner and outer enclosure surfaces

Sh, =

local and average Sherwood number

T=

dimensional temperature

T=

vector equal to SE

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

ρ=

density

θ=

dimensionless temperature

ε=

porosity

ψ=

stream function

σ=

dimensionless solute concentration

Γ=

diffusion coefficient

Subscripts=
c=

cold wall or convection heat transfer

C=

refers to main grid point

f=

refers to element face

F=

refers to the F grid point

h=

hot wall

i=

condition at inner surface

m=

mass transfer

o=

condition at outer surface

s=

refers to solute concentration

t=

refers to temperature

Additional information

Funding

The financial support provided by the University Research Board of the American University of Beirut URB 2017-18 is gratefully acknowledged.

Log in via your institution

Log in to Taylor & Francis Online

PDF download + Online access

  • 48 hours access to article PDF & online version
  • Article PDF can be downloaded
  • Article PDF can be printed
USD 61.00 Add to cart

Issue Purchase

  • 30 days online access to complete issue
  • Article PDFs can be downloaded
  • Article PDFs can be printed
USD 716.00 Add to cart

* Local tax will be added as applicable

Related Research

People also read lists articles that other readers of this article have read.

Recommended articles lists articles that we recommend and is powered by our AI driven recommendation engine.

Cited by lists all citing articles based on Crossref citations.
Articles with the Crossref icon will open in a new tab.