ABSTRACT
In this study, mixed convection in a cavity that has a fluid and superposed porous medium with an adiabatic rotating cylinder is numerically investigated. The bottom horizontal wall is heated and the top horizontal wall is cooled while the remaining walls are assumed to be adiabatic. An adiabatic rotating cylinder is inserted inside the cavity. The governing equations are solved by the Galerkin weighted residual finite element method. The effects of Rayleigh number (between 103 and 106), angular rotational speed of the cylinder (between 0 and 6,000), Darcy number (between 10−5 and 10−2), cylinder sizes (between R = 0.1 and R = 0.3) and three different vertical locations of the cylinder on the fluid flow and heat transfers characteristics are numerically investigated. It is observed that the cylinder size has a profound effect on the local and averaged heat transfer. The local and averaged heat transfers generally increase and the convection is more effective in the upper half of the cavity as the Rayleigh number and Darcy number enhance. The averaged heat transfers increases with the cylinder size until Ra = 105. The averaged heat transfer increases almost linearly with the angular rotational velocity of the cylinder and the increase rate becomes higher as the cylinder size increases. The local and averaged heat transfers enhances/deteriorate as the cylinder approaches the upper/lower wall of the cavity.
Nomenclature
Cp | = | specific heat at constant pressure (J kg−1 K−1) |
Da | = | Darcy number, K/L2 |
g | = | gravitational acceleration (m s−2) |
Gr | = | Grashof number, |
L | = | side length of the enclosure (m) |
h | = | Local heat transfer coefficient(W m−2K−1) |
k | = | Thermal conductivity (W m−1 K−1) |
K | = | Permeability of the porous medium (m2) |
Nu | = | Nusselt number, Nu = hL/kf |
p, P | = | pressure (N/m2), dimensionless pressure |
Pr | = | Prandtl number, Pr = νf/αf |
r | = | radius of the inner cylinder (m) |
R | = | Dimensionless radius of the inner cylinder |
R = r/L | = | |
Ra | = | Rayleigh number, gβf(Th − Tc)L3/νfαf |
Re | = | Reynolds number, Re=ωr L/νf |
Ri | = | Richardson number, Ri = Gr/Re2 |
T | = | temperature (K) |
u, v | = | Velocity components (ms−1) |
U, V | = | Dimensionless velocities components, V = vL/αf, U = uL/αf |
x, y | = | Dimensional coordinates (m) |
X, Y | = | Dimensionless coordinates, X = x/H, Y = y/H |
Greek symbols | = | |
α | = | thermal diffusivity (m2 s−1) |
β | = | rhermal expansion coefficient (K−1) |
ε | = | porosity of the porous layer |
ν | = | Kinematic viscosity (m2s−1) |
θ | = | Dimensionless temperature, θ = (T–Tc) / (Th–Tc) |
Ψ | = | dimensionless stream function |
ω | = | rotational speed of the inner cylinder (rad s−1) |
Ω | = | dimensionless rotational speed of the inner cylinder, Ω = ωL2/αf |
ρ | = | Density (kg m−3) |
μ | = | Dynamic viscosity (N.s m−2) |
Subscripts | = | |
c | = | cold |
f | = | fluid |
h | = | hot |
o | = | cavity center, standard conditions |
p | = | porous |
Nomenclature
Cp | = | specific heat at constant pressure (J kg−1 K−1) |
Da | = | Darcy number, K/L2 |
g | = | gravitational acceleration (m s−2) |
Gr | = | Grashof number, |
L | = | side length of the enclosure (m) |
h | = | Local heat transfer coefficient(W m−2K−1) |
k | = | Thermal conductivity (W m−1 K−1) |
K | = | Permeability of the porous medium (m2) |
Nu | = | Nusselt number, Nu = hL/kf |
p, P | = | pressure (N/m2), dimensionless pressure |
Pr | = | Prandtl number, Pr = νf/αf |
r | = | radius of the inner cylinder (m) |
R | = | Dimensionless radius of the inner cylinder |
R = r/L | = | |
Ra | = | Rayleigh number, gβf(Th − Tc)L3/νfαf |
Re | = | Reynolds number, Re=ωr L/νf |
Ri | = | Richardson number, Ri = Gr/Re2 |
T | = | temperature (K) |
u, v | = | Velocity components (ms−1) |
U, V | = | Dimensionless velocities components, V = vL/αf, U = uL/αf |
x, y | = | Dimensional coordinates (m) |
X, Y | = | Dimensionless coordinates, X = x/H, Y = y/H |
Greek symbols | = | |
α | = | thermal diffusivity (m2 s−1) |
β | = | rhermal expansion coefficient (K−1) |
ε | = | porosity of the porous layer |
ν | = | Kinematic viscosity (m2s−1) |
θ | = | Dimensionless temperature, θ = (T–Tc) / (Th–Tc) |
Ψ | = | dimensionless stream function |
ω | = | rotational speed of the inner cylinder (rad s−1) |
Ω | = | dimensionless rotational speed of the inner cylinder, Ω = ωL2/αf |
ρ | = | Density (kg m−3) |
μ | = | Dynamic viscosity (N.s m−2) |
Subscripts | = | |
c | = | cold |
f | = | fluid |
h | = | hot |
o | = | cavity center, standard conditions |
p | = | porous |