FMFP10 - AM - 04 Fluid Flow Simulations Within Rotating Annulus
FMFP10 - AM - 04 Fluid Flow Simulations Within Rotating Annulus
FMFP10 - AM - 04 Fluid Flow Simulations Within Rotating Annulus
4th International
Conference
Fluid Mechanics
and Power
Fluid Power
ProceedingsProceedings
of the 37th International
& 4th &
National
Conference
on Fluidon
Mechanics
and Fluid
December 16-18, 2010, IIT Madras, Chennai, India.
FMFP10 - AM - 04
FLUID FLOW SIMULATIONS WITHIN ROTATING ANNULUS
Arun K Sukumaran
College of Engineering
Thiruvananthapuram, Kerala
Reji R.V
College of Engineering
Thiruvananthapuram, Kerala
[email protected]
K.S.Santhosh
College of Engineering
Thiruvananthapuram, Kerala
ABSTRACT
The fluid flow associated with heat transfer in
rotating annuli exists in many engineering
applications. The present study presents the
effect of the variation of rotational Reynolds
number (Re) and axial Reynolds number (Rea)
on heat transfer and fluid flow in rotating
annulus for a radius ratio (j= ri/ro) 0.5. The
computational work using a CFD software
FLUENT is validated with experimental work.
The wall of inner rotating cylinder is kept at
constant heat flux and stationary outer wall is at
operating temperature. Results for various Re
and Rea are presented. Axial velocity, swirl
velocity, inner wall temperature profiles and
Nusselt number variations are also presented.
Literature Review
Convective heat transfer from a horizontal
rotating cylinder has been studied by many
authors. Most of them concluded that rotation
has no significant effects on the heat transfer
coefficient for low rotational speeds and that the
heat transfer is governed by free convection.
However, for very high rotational speeds, the
forced convection is the predominant heat
transfer regime and the average Nusselt number
on the cylinder surface is often obtained with
correlations such as Nu=aReb where a and b
being constants.
and
Nu
=8.854Pr Re
respectively.
NUMERICAL ANALYSIS
In the present work, numerical analysis of fluid
flow and heat transfer in a rotating annulus with
various radius ratios and various axial flow is
considered. The inner cylinder rotates with a
uniform angular velocity while outer cylinder is
kept stationary. The surface of the inner
cylinder is subjected to uniform heat flux and
surface of outer cylinder is in isothermal
(operating) condition. A validation study was
conducted to check the compatibility of CFD
software FLUENT 6.3 for predicting heat
transfer in rotating annulus. The experimental
study mentioned in Sheng-Chung Tzeng, 2006
was taken and the dimension and boundary
conditions are being exactly same as that in
experimental setup. The present analysis
presents how the variation of rotational
Reynolds number (Re) and axial Reynolds
number (Rea) influence the heat transfer and
fluid flow in rotating annulus.
Physical Model
The present study is the fluid flow and heat
transfer through rotating annulus with axial
flow. Figure 1 shows the physical model of the
present study, the inner cylinder is rotating with
an angular velocity rad/sec and outer cylinder
is kept stationary. Assuming the flow is
incompressible and steady, the outer cylinder is
kept in isothermal condition and inner cylinder
wall subjected to constant heat flux. The inlet
velocity of the fluid is uniform ui and entering at
Governing equations
The incompressible Navier-Stokes equations,
energy equation coupled with k-e turbulence
model in single rotating frame were considered
for this case. Additional terms originated from
the relative frame formulation were treated a
source terms in the FLUENT formulations.
Axisymmetric swirl
Steady
Standard k-epsilon
turbulence model
Simple (semi implicit
pressure linked equation)
Air
Copper
Copper
Validation
The experimental work of Sheng Chung Tzeng,
2006 was taken for validation of the
methodology.
Property
Density
Air
Specific
heat
Thermal
conductivit
y
Viscosity
Thermal
expansion
coefficient
Density
Specific
Copp heat
er
Thermal
conductivit
y
Density
Specific
Bake heat
lite
Thermal
conductivit
y
Unit
Method
Value
Kg/m3
Boussin
esq
1.225
J/kg-K
Constant
1006.
43
W/mK
Constant
0.024
2
Kg/ms
Constant
1.789
x10-05
1/K
Constant 0.003
Kg/m3
Constant 8978
J/kg-K
Constant 381
W/mK
Constant 387.6
Kg/m3
Constant 1280
J/kg-K
Constant 1590
W/mK
Constant 0.23
Stationary
Stationary
Rotating
Stationary
Type
Isothermal
Isothermal
Constant heat flux
850 W/m2
Isothermal
Axisymmetric swirl
Steady
Standard k-epsilon
turbulence model
Simple (semi implicit
pressure linked equation)
Air
Bakelite
Copper
Copper
Copper
, 2400 Re
45,000.) and obtained from the numerical
analysis is shown in fig. 4. The numerical results
are slightly over predicted for all cases of
Reynolds number. The error of numerical value
with correlated value in between 4 to 15%.
Reason of this error may be due to discretization
and numerical errors.
CONCLUSIONS
A numerical study and a detailed parametric
analysis on the problem have been conducted on
the rotating heated annulus.
The current
problem is validated using an experimental work
and the analysis is in good agreement with
correlation. The parameters like axial and swirl
velocity components, temperature distribution
and Nusselt number variations are studied and
discussed.
As rotational Reynolds number increases
the magnitude of axial velocity will
decreases only at low axial Reynolds
number but at high axial Reynolds
numbers the effect of rotation will
diminish.
The axial Reynolds number has no
significant effect on swirl velocity. As
the rotational Reynolds number increases
the hydrodynamic boundary layer
become thinner and is nearer to inner
rotating wall.
As
rotational
Reynolds
number
increases, wall temperature decreases
which may due to high rate of transfer of
thermal energy from the inner wall to
fluid.
As axial Reynolds number increases
more heat is transferred from inner
rotating wall to the fluid which is
flowing through annulus.
As
rotational
Reynolds
number
increases, Nusselt number also increases
for low axial Reynolds number and at
high axial Reynolds number the effect
diminishes.
As rotational speed increases more heat
is transferred from shaft to fluid, which
NOMENCLATURE
Ming-I Char, Yuan-Hsiung Hsu, 1998
Numerical prediction of turbulent mixed
convection in a concentric horizontal rotating
annulus with low-re two-equation models,
International Journal of Heat and Mass Transfer,
Vol. 41, pp. 1633-1643.
Greek symbols
- Density of fluid, kg/m3
w - Wall shear stress, N/m2
- Angular velocity of inner cylinder, rad/s
Subscripts/superscripts
a- Axial
- Rotational
i Inner, Inlet
o - Outer, Outlet
w Wall
s Swirl
ot Operating temperature
REFERENCES
Escudier, I W Gouldson, 1995. Concentric
annular flow with center body rotation of a
newtonian and a shear- thinning liquid,
International Journal of Heat and Fluid Flow
Vol. 16, 156-162.
Joo-Sik Yoo, 1998. Mixed convection of air
between two horizontal concentric cylinders