ABSTRACT
Composite cavities formed by a clear space, a layer of porous material, and a solid plate can be engineered for controlling the overall heat transfer across the enclosure. Using different layer dimensions, as well as distinct porous and solid materials, the value of the cavity Nusselt number can be modified with regard to traditional Nu ∝ Ran behavior, which is encountered either in completely empty cavities or in cavities fully fitted with porous materials. Motivated by such novel application, this work presents a study about turbulent natural convection in a composite concentric annulus. The annulus is assumed to be two-dimensional and positioned horizontally, being isothermally heated at the inner cylinder and cooled from the outer surface. Laminar flow is considered in addition to the turbulent regime, which is handled via the standard k–ε model. The wall treatment applied is the High Reynolds approach. The Two-Energy Equation Model (2EEM) is utilized in the porous section. The transport equations are discretized using the control-volume method. The system of algebraic equations is relaxed via the Semi Implicit Pressure-Linked Equations (SIMPLE) algorithm. A new numerical methodology is applied to resolve all three layers in a single computational domain by establishing two temperature sets, defined according to the location inside the composite structure. Nusselt number behavior shows that for Rayleigh number up to 104 there is no significant variation between the laminar and turbulence models, although the differences increase when the flow gets more intense and/or the porous material becomes more permeable. When comparing the effects of Rayleigh number, Darcy number, porosity, and thermal conductivity ratio between the solid and the fluid on Nu, the results indicate that the solid-phase properties have a greater influence in enhancing the overall heat transferred through the cavity.
Nomenclature
cF | = | Forchheimer coefficient |
c1,2,3,k,μ | = | model constants |
cp | = | specific heat |
D | = | deformation rate tensor, |
Da | = | Darcy number, |
D | = | particle diameter |
g | = | gravity acceleration vector |
= | generation rate of | |
= | generation rate of | |
h | = | heat transfer coefficient |
I | = | unit tensor |
K | = | permeability, |
k | = | turbulent kinetic energy |
= | fluid thermal conductivity | |
= | solid thermal conductivity | |
= | conductivity tensor due to dispersion | |
= | conductivity tensor due to tortuosity | |
= | characteristic length, | |
ni | = | unit vector normal to Ai |
Nu | = |
|
Pr | = | Prandtl number |
R | = | internal to external radius ratio, |
r* | = | nondimensional radius, |
ri | = | internal radius |
ro | = | external radius |
Ra | = |
|
Ram | = |
|
T | = | temperature |
= | nondimensional temperature of the fluid phase or the solid block | |
= | nondimensional temperature of the solid material of the porous matrix | |
u | = | microscopic velocity |
= | Darcy or superficial velocity (volume average of u) | |
x,y | = | Cartesian coordinates |
α | = | thermal diffusivity |
β | = | thermal expansion coefficient |
= | representative elementary volume | |
= | fluid volume inside | |
ε | = | dissipation rate of turbulent kinetic energy, |
μ | = | dynamic viscosity |
ν | = | kinematic viscosity |
ρ | = | density |
= | nondimensional constants | |
ϕ | = | porosity, |
Special characters | = | |
φ | = | general scalar variable |
= | time average | |
= | time fluctuation | |
= | intrinsic average | |
= | volume average | |
= | spatial deviation | |
= | absolute value (Abs) | |
= | general vector variable | |
= | effective value, | |
= | solid/fluid | |
= | hot/cold | |
= | macroscopic value | |
= | transpose |
Nomenclature
cF | = | Forchheimer coefficient |
c1,2,3,k,μ | = | model constants |
cp | = | specific heat |
D | = | deformation rate tensor, |
Da | = | Darcy number, |
D | = | particle diameter |
g | = | gravity acceleration vector |
= | generation rate of | |
= | generation rate of | |
h | = | heat transfer coefficient |
I | = | unit tensor |
K | = | permeability, |
k | = | turbulent kinetic energy |
= | fluid thermal conductivity | |
= | solid thermal conductivity | |
= | conductivity tensor due to dispersion | |
= | conductivity tensor due to tortuosity | |
= | characteristic length, | |
ni | = | unit vector normal to Ai |
Nu | = |
|
Pr | = | Prandtl number |
R | = | internal to external radius ratio, |
r* | = | nondimensional radius, |
ri | = | internal radius |
ro | = | external radius |
Ra | = |
|
Ram | = |
|
T | = | temperature |
= | nondimensional temperature of the fluid phase or the solid block | |
= | nondimensional temperature of the solid material of the porous matrix | |
u | = | microscopic velocity |
= | Darcy or superficial velocity (volume average of u) | |
x,y | = | Cartesian coordinates |
α | = | thermal diffusivity |
β | = | thermal expansion coefficient |
= | representative elementary volume | |
= | fluid volume inside | |
ε | = | dissipation rate of turbulent kinetic energy, |
μ | = | dynamic viscosity |
ν | = | kinematic viscosity |
ρ | = | density |
= | nondimensional constants | |
ϕ | = | porosity, |
Special characters | = | |
φ | = | general scalar variable |
= | time average | |
= | time fluctuation | |
= | intrinsic average | |
= | volume average | |
= | spatial deviation | |
= | absolute value (Abs) | |
= | general vector variable | |
= | effective value, | |
= | solid/fluid | |
= | hot/cold | |
= | macroscopic value | |
= | transpose |