Fans Design and Operation of Centrifugal, Axial-Flow, and Cross-Flow Fans
Fans Design and Operation of Centrifugal, Axial-Flow, and Cross-Flow Fans
2014–9:47am] [1–11]
//blrnas3.glyph.com/cenpro/ApplicationFiles/Journals/SAGE/3B2/PIAJ/Vol00000/140060/APPFile/SG-PIAJ140060.3d (PIA) [PREPRINTER stage]
Original Article
Abstract
In this study, we investigated the effects of volute tongue geometry variation on the head, efficiency, and radial force of a
centrifugal pump. Numerical simulation modeling based on k ! turbulence with automatic near wall treatments was
used to simulate the turbulent flow. The effect of blade position with respect to the volute tongue on instantaneous
pump characteristics was investigated. The parametric studies were done for cutwater gap, tongue shape, and volute
tongue angle. Numerical results showed that the large cutwater gap caused lower radial force, especially at high flow
rates. Investigations using various volute tongue shapes indicated that the short tongue volute decreased the radial force
at design and low flow rates. Considering all aspects, the most satisfactory volute tongue angle was found to be 5 less
than the outlet velocity angle of the impeller; yielding about 40% lower radial force than others at the design point.
Keywords
Pump, volute, low specific speed, computational fluid dynamics, radial force, cutwater gap, tongue
somewhat widened. Furthermore, they found that this effects of cutwater gap, and tongue shape and angle on
replacement shifts the maximum efficiency point the performance and radial force of low-Ns centrifugal
towards a higher flow rate. Although Guo et al. con- pumps.
cluded that all variations are only affected by tongue
profile replacement; it seems that the cutwater gaps
Main pump geometry
have not been set the same for all tongue shapes in
their study. Uchida et al.5 made measurements to The experimental measurement of Kelder et al.9 was
monitor radial force and pump characteristic curves used for CFD validation throughout this study. The
for different cutwater gaps and tongue shapes. investigated centrifugal pumps have a low-Ns of 21 r/
According to their report, pump characteristic and min and their design parameters are summarized in
radial force are greatly affected by tongue area size Table 1.
and tongue shape. In addition, the measured max- Experiments were performed at flow rates of
imum instantaneous radial force was in agreement 82.5%, 100%, and 112% of the design flow rate. An
with the Stepanoff equation.6 Performing experiments LDV measurement system was used for velocity meas-
and numerical simulations, Dong et al.7 showed that urements. Velocity and static pressure data were col-
radial force depends on the geometry of the region lected at the locations shown in Figure 1, with the
where the gap between impeller and volute is small, pressure tabs located at the shroud side, just outside
including both tongue shapes and blade-tongue gap. the impeller yet in volute at mid-height of the outer
More recently, Torabi and Nourbakhsh8 studied the wall. Further details of the pump configuration and
effects of volute geometrical factor variation on radial experimental methods can be found in Kelder et al.9
force. They observed that the gap between impeller We designed eight volutes (details discussed in next
and tongue severely affects the radial force. More spe- sections) to study the effects of volute tongue geom-
cifically, they found that a 6% reduction in the cut- etry on head, efficiency and radial force. Three
water gap increases the amplitude of the static
pressure variation by 50% around the impeller exit.
Kelder et al.9 studied turbulent flow in the volute of a Table 1. Main parameters of pump.
low specific speed pump close to the design point.
Parameter Value
Both experimental and numerical approaches
showed that, under this condition, the core flow Metric specific speed, Ns (r/min) 21
behaves like a potential flow (i.e. no boundary layer Nominal flow coefficient, n 0.15
separation occurs). Thus, they employed inviscid Nominal head coefficient, n 0.124
equations to solve flow field through the low-Ns Impeller exit diameter, d2 (mm) 640
pump. Impeller exit width, b2 (mm) 25
Due to the established importance of tongue geom-
Number of blades, Z 7
etry for pump performance, here we investigated the
Figure 1. (a) kelder volute geometry and locations of velocity and static pressure measurement; (b) volute cross-sectional
dimensions.
Alemi et al. 3
geometrical parameters (shape, cutwater gap, and boundary condition effects on simulation of the
tongue angle) of volute tongue were considered. pump flow domain. Figure 2(b) shows the generated
grids for the pump impeller and volute. Due to com-
Numerical method and model plexity of geometry, the hexahedral (near walls) and
tetrahedral cells were used for inlet and outlet mod-
description
ules and tetrahedral cells were used for impeller and
Commercial software CFX-13.0 was used to solve the volute. Mesh clustering was applied close to some
full 3D Reynolds averaged Navier–Stokes (3D- important regions, such as leading and trailing edge
RANS) equations for various working conditions. of the blades and the volute tongue where flow separ-
The solver was a 3D CFD code in which the gov- ation may occur (Figure 2(b)). During grid generation
erning equations are discretized using the finite-ele- process, orthogonal quality, aspect ratio, and skew-
ment-based finite volume method. Such a hybrid ness were assured to be in desirable range. For steady
CFD method is advantageous because of its geometric simulation the grids between the impeller and volute
flexibility (due to the finite element method) and in one hand, and suction pipe and impeller on the
retention of the conservation properties of the finite other hand, were conjoined by means of a frozen
volume method, which results in low numerical error rotor interface. In other words, for steady simulation,
in the case of nonsmooth grids.10 For space discret- impeller rotation was considered and the frame of ref-
ization, the second-order upwind scheme was applied. erence was changed but relative orientation of the
For grid generation the 3D models of impellers and components across the interface was fixed. Both two
volute were first generated, and then exported to com- frames of reference were connected in such a way that
mercial software. Water at 25 C was selected as the each of them had a fixed relative position throughout
working fluid. Mass flow rate with flow direction and the calculation. This was equivalent to taking a snap-
constant pressure were implemented for inlet and out- shot from the flow field in an instance of time.
let boundary conditions respectively. Spence and Moreover, general grid interface was employed
Amaral-Texeira10 reported that these boundary con- between the volute and discharge duct.
ditions cause rapid convergence. As shown in Figure 3, grid independency studies
It was found that by increasing the turbulent inten- were performed using various mesh sizes. We found
sity at inlet from 5% to 10% and viscosity ratio from that increasing the number of cells from 7 105 to
10 to 100, the predicted head changes less than 0.01%. 1 106 results in a 0.8% change in the predicted
Thus, turbulent intensity and viscosity ratio were set head. Further mesh refinement from 1 106 cells to
to 5% and 10%, respectively. The convergence cri- 1:9 106 cells decreases the head by about 0.7%.
teria for numerical simulation were set at maximum Since the predicted head did not show a strong sensi-
residuals of 1 104 . tivity to grid refinement, the numerical results can be
For mesh generation, the pump was divided into considered as grid independent. All the numerical
four components: (i) inlet duct, (ii) impeller, (iii) results were obtained using the 1 106 cell mesh.
volute, and (iv) outlet duct (Figure 2(a)). Straight Unsteady simulation was done just in main geom-
ducts at inlet and outlet were used to prevent etry using k ! model at the design point and
Figure 2. (a) 3D models of total domain including: (i) inlet duct, (ii) impeller, (iii) volute, and (iv) outlet duct; (b) computational mesh
for the investigated pump.
Alemi et al. 5
The undulating pressure distribution seen in CFD models returned similar results, consistent with the
predictions (Figure 4) is a result of impeller blades experimental data of Kelder et al.9 It should be
effect. It is seen that all three turbulence models give noted that discrepancies between computations and
comparable results in both design and off-design con- experimental data could be partly due to the uncer-
ditions. Despite the presence of local discrepancies tainties in measurements.
between CFD predictions and measurements, both Kelder et al.9 employed Euler equations in their
the level and trend of CFD results are in overall agree- numerical predictions and found reasonable agree-
ment with the measurements. In Figure 5, the pre- ment between computational and experimental data.
dicted radial and tangential velocity components However, level of agreement between the numerical
were compared with the experimental data of Kelder and experimental data reported here is better than
et al.9 Comparisons were made for three working con- that of Kelder.
ditions (=n ¼ 0:825, 1.0, and 1.12) along four tra- In Figure, 6 the predicted coefficients using
verses (A, D, F, and H) of the volute (Figure 1). three turbulent models with steady simulation and
Although, some discrepancies between turbulence just k ! model with unsteady simulation were com-
models at =n ¼ 0:825 are seen here, all turbulence pared with experimental data for a wide range of flow
Figure 5. Comparison of predicted and measured nondimensional radial and circumferential velocities in the volute along different
traverses at =n ¼ 0.825, 1, and 1.12.
Figure 6. Comparison of predicted head and radial force using various turbulence models with the experimental data.
coefficients, including the design point and off-design respect to the tongue at design condition. The results
conditions. Firstly, in steady simulations, although indicate that when tongue locations vary between
the k " turbulence model returns more accurate two blades, the head magnitude, efficiency, and
results for the nondimensional head near the design radial forces will be altered by less than 2.5%,
point, the head predictions of the k ! model at off- 1.5%, and 8%, respectively. However, for steady
design conditions are closer to the experimental data, simulation of present work, in all cases the tongue
especially at high flow rate (Figure 6(a)). The k ! was similarly located in the mid-position of two
turbulence model returns radial force predictions in blades ( ¼ 25:7 Þ.
better agreement with experiments than other
models (Figure 6(b)). Secondly, as shown in
Figure 6, although unsteady simulation using k !
Volutes with various tongues profiles
turbulence model yields the closest head coefficients The influence of various volute tongue profiles on the
to the measurement data, the estimated radial forces radial force and pump characteristics was
by unsteady simulation are more or less similar to investigated.
steady simulation. As the principal focus of the cur-
rent study is radial force prediction at off-design con- Cutwater gap. The radial gap between the leading edge
ditions the k ! turbulence model with steady of the impeller blades and volute tongue is called the
simulation was employed for calculations. cutwater gap. This distance affects the interaction
between blade and tongue. Therefore, cutwater gap
could change the pump hydraulic behavior, especially
Impeller blade and volute tongue interaction
radial forces. It is expected that cutwater gap incre-
When the blade position relative to tongue ðÞ ment could decrease the radial force. To investigate
changes, pressure and velocity distribution in the this prediction, three volutes with different cutwater
volute vary due to interaction between the impeller gaps (equal to 3.2%, 5%, and 8.4% of the impeller
blade and volute tongue (i.e. pump characteristics radius) were designed. Volutes with maximum and
are instantaneous and unsteady). Figure 7 presents minimum cutwater gap are displayed in Figure 8
pump performance curves versus blade location with and numerical results for these volutes are shown in
Alemi et al. 7
Figure 9. It was observed that all three volutes have tongue angle of 1 , whereas the reverse condition
identical head and efficiency. As expected, large cut- occurs at a high flow rate. At design point, all
water gap generates lower radial force, especially at tongue angles produce more or less equal head and
design and high flow rate (Figure 9(c)). efficiency and the 1 volute tongue angle experiences
the least radial force. This observation is in good
Tongue angle. Volutes tongue angle (v ) (Figure 10) is agreement with Figure 12, which shows that the pres-
one of the principal geometrical parameters of a sure field is more uniform around the impeller with a
volute and can affect hydraulic performances of a 1 volute tongue angle at design condition, compared
pump. In order to avoid any shock and separation to other angles.
during volute design, the volute angle should be the-
oretically equal to the absolute angle of impeller out- Tongue shape. Three tongue shapes (short, middle, and
flow (2 ) at design conditions. In the present pump sharp) with different t (12 , 18 , and 22 ) were
test case, the absolute angle of impeller outlet velocity designed to study their influence on pump character-
at the design flow rate was obtained at about 6 . To istics (Figure 13). Numerical results of these volutes
investigate the effects of tongue angle variation on the are illustrated in Figure 14. Short tongue volutes pro-
pump performance, four volutes with various tongue duce lower radial force, especially at low capacity.
angles (1 , 6 , 11 , and 16 ) were designed (Figure 10). While all volute designs yield more or less the same
Numerical calculated results are demonstrated in head and efficiency, the volute with sharp tongue is
Figure 11. The volutes with tongue angle of 1 and slightly more efficient at off-design conditions.
16 give higher efficiency at low and high capacity
respectively. This observation might be expected as
the absolute angle of impeller outflow (2 ) at low cap-
Conclusion
acity decreases due to absolute velocity reduction. In this paper, the hydraulic characteristics of the cen-
This results in a better match between outflow and a trifugal pumps for various volute designs are
Alemi et al. 9
Figure 12. Pressure contours in various tongue angles at design flow rate ðn Þ.
Alemi et al. 11
11. Jones WP and Launder BE. The prediction of laminar- Q, Qn flow rate; normal or design flow rate
ization with a two equation model of turbulence. Int J (m3/s)
Heat Mass Transfer 1972; 15: 301–314. r4 volute inlet radius (mm)
12. Wilcox DC. Multiscale model for turbulent flows. In: ro radius of volute wall midline (mm)
AIAA 24th aerospace sciences meeting. Atlanta: AIAA,
r (r–r4)/(ro–r4) radial distance in volute,
1986.
nondimensional
13. Menter FR. Multiscale model for turbulent flows. In:
AIAA 24th fluid dynamics conference. Atlanta: AIAA, Re ð
d22 Þ= Reynolds number,
1993. nondimensional
14. Menter FR. Two-equation eddy-viscosity turbulence Tave: average pressure and viscous moment
models for engineering applications. AIAA J 1994; 32: on impeller (N.m)
1598–1606. u2 blade tip velocity (m/s)
Vr Vr =u2 , radial velocity, nondimensional
V V =u2 , circumferential velocity,
nondimensional
Appendix Z number of blades
Notation v tongue angle ( )
" dissipation rate (m2/s3)
b2 breadth of impeller at discharge (mm)
hyd: gQH=ðTave:
Þ, hydraulic efficiency
d2 outside diameter of impeller (mm)
(%)
Fr radial force (N)
kinematic viscosity (m2/s)
Fr Fr =ð0:5u22 d2 b2 Þ, nondimensional
! specific dissipation
rate (1/s)
g gravity (m/s2)
, n Q=
d22 b2 , flow coefficient; flow coef-
H, Hn head; normal or design head (m)
ficient at design flow rate Qn,
k turbulent kinetic energy (m2/s2)
nondimensional
N rotation speed (r/min)
, n gH=ð
2 d2 Þ, head coefficient; head
Ns NQ0:5 H0:75 , metric specific speed
coefficient at design head rate Hn,
(r/min)
nondimensional
P static pressure (Pa)
fluid density (kg/m3)
Pin inlet static pressure (Pa)
t angular position of tongue ( )
P ðP Pin Þ=u22 , nondimensional