Next Article in Journal / Special Issue
Sensitivity Analysis of BLISK Airfoil Wear
Previous Article in Journal / Special Issue
Component-Specific Preliminary Engine Design Taking into Account Holistic Design Aspects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Characterization of Unsteady Forces Triggered by Cavitation on a Centrifugal Pump †

1
Chemical Propulsion, SITAEL S.p.A., 56121 Pisa, Italy
2
Department of Civil and Industrial Engineering, University of Pisa, 56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in Proceedings of the 17th International Symposium on Transport Phenomena and Dynamics of Rotating Machinery (ISROMAC 2017).
Int. J. Turbomach. Propuls. Power 2018, 3(2), 13; https://doi.org/10.3390/ijtpp3020013
Submission received: 6 February 2018 / Revised: 27 April 2018 / Accepted: 2 May 2018 / Published: 7 May 2018

Abstract

:
The paper presents an experimental campaign aimed at the characterization of the relationship between cavitation-induced instabilities and forces acting on the shaft, relevant to space application turbopumps. The experiments have been carried out on a six-bladed unshrouded centrifugal turbopump. Pressure fluctuations are analyzed in their frequency content for understanding the instability nature (axial, rotating) and their main characteristics (e.g., amplitude, rotating direction). The spectral analysis of the force components highlights a strong relationship of the z-component (along the rotating axis) with axial instabilities. On the other hand, rotating cavitation may involve force oscillations along all the three components leading to unwanted and dangerous fluctuating unbalances perpendicular to the rotating axis.

1. Introduction

Propellant feed turbopumps are a crucial component of all liquid propellant rocket engines due to the severe limitations associated with the design of dynamically stable, high power density machines capable of meeting the extremely demanding suction, pumping, and reliability requirements of modern STSs (space transportation systems).
The attainment of such high power/weight ratios is invariably obtained by running the impeller at the maximum allowable speed and low shaft torque. Current configurations in space applications are typically characterized by the presence of lighter, but also more flexible, shafts. Therefore, since the operation under cavitating conditions is tolerated, the turbopump is exposed to the onset of dangerous fluid dynamic and rotordynamic instabilities which may be triggered by cavitation phenomena. When designing a turbomachine, particularly if it has to operate at high rotational speeds, it is important to be able to predict the fluid-induced forces acting on the various components of the turbomachine. The study of radial and rotordynamic forces on turbomachines, by means of analytical/numerical and experimental approaches has been extensively carried out in the last 50 years by many researchers all around the world [1,2,3,4,5,6,7,8]. However, the experimental characterization of the cavitation influence on these phenomena is still very poor even if it is a common operational condition in space application. Moreover, it is extremely important since the occurrence of cavitation drastically modifies the inertia of the fluid surrounding the impellers and, in turn, the critical speeds of the machine.
Fluid instabilities developing in space inducers have been widely studied in the past [9,10,11,12,13,14,15,16,17]. However, poor evidences in open literature on experimental activities for centrifugal pumps do not allow a complete understanding of the cavitation induced flow phenomena. In fact, few detailed experimental results can be found in literature focused on the instabilities and their interaction with the system. However, rotating cavitation [18] and auto-oscillation [19] have been detected on centrifugal pumps and studied in previous works.
The experimental characterization of the unsteady fluid forces/moments acting on space turbopump impellers as a consequence of the onset of the most dangerous types of cavitation-induced instabilities is here described and analyzed.
This is done by employing a new methodology, consisting of measuring simultaneously the forces acting on the shaft, by means of a rotating dynamometer previously employed in past experimental activities [20,21,22], and the pressure field around the centrifugal pump at different locations, by means of piezoelectric pressure transducers. This technique allows for comparing and evaluating the influence of pressure fluctuations on the unsteady force component that the pump exerts on the shaft. The nature of the unsteady flow phenomena interacting with the system can be therefore assessed from the knowledge of mechanical stresses.

2. Materials and Methods

2.1. Experimental Apparatus

Figure 1 shows the general layout of the cavitating pump rotordynamic test facility (CPRTF, see [23] for further information) at SITAEL S.p.A. (Pisa premise, Italy) (formerly ALTA S.p.A.). The facility is specifically intended for the experimental analysis of the relevant phenomena on space application turbopumps with particular attention on cavitation-related phenomena. The facility operates with water under similarity conditions in order to simulate the behavior of cryogenic fluids [24].
In order to relate the fluid instabilities with the unsteady forces acting on the pump, the experiments have been carried out by jointly using a rotating dynamometer (extensively employed in past activities, [25]), and a suitable number of piezoelectric pressure transducers from PCB Piezotronics (Depew, NY, USA), mod. S112A22 3.45 bar range, flush-mounted on the machine casing. Table 1 reports the assessed dynamometer errors with a 95% confidence level.
The dynamometer has been placed between the pump and the rotor shaft in order to avoid the influence of seals and bearing forces. The positions of the pressure transducers have been chosen in order to evaluate and identify the nature of the fluid instability (rotational/axial), thus different axial and azimuthal positions have been exploited. For this purpose, three different axial positions have been made available for the transducers on the pump casing (Figure 2 and Figure 3): upstream, midstream, and downstream (at the diffuser). With reference to Figure 3, the relative positions of the pressure taps have been designed in such a way that, at the same instant, number-like taps should face the through-flow of the same blade channel. The taps exploited during the present campaign are reported in Table 2.
With reference to Figure 1 and Figure 2, the test setup includes:
  • two electromagnetic flowmeters (Fisher-Rosemount, Saint Louis, MO, US, model 8732E, range: 0–100 l/s, accuracy 0.5% FS (full scale)), for the measurement of the inlet and outlet flow rates;
  • an absolute pressure transducer ( p i n @IPT1) placed six diameters upstream, with respect to the blade leading edge (General Electric, Boston, MA, US, model Druck PMP1400, range 0–1 bar, accuracy 0.25% FS), for the assessment of the inlet cavitation number ( σ ) ;
  • an absolute pressure transducer ( p C @OPT1) placed in the test chamber/section in order to measure the influence of the pressure on the axial force measured by the dynamometer. In fact, there is a measured force related only to the pressure difference between the water pressure and the sealed dynamometer which is internally at atmospheric condition (General Electric, Boston, MA, US, model UNIK5000, range 0–6 bar, accuracy 0.1% FS);
  • one differential pressure transducer ( Δ p @IPT1 to OPT2) which measures the pump pressure rise between the inlet station placed at 6 diameters upstream of the blade leading edge, and the outlet station placed about two diameters downstream of the blade trailing edge (General Electric, Boston, MA, US, model UNIK5000, range 0–5 bar, accuracy 0.1% FS);
  • one temperature sensor PT100 (range 0–100 °C, accuracy 0.5 °C).
The test item employed in the present experimental campaign has been designed at SITAEL and is a six-bladed radial pump named VAMPUFF (Figure 4) whose main characteristics are summarized in Table 3. The design is based on the reduced order model developed at SITAEL S.p.A. (formerly ALTA S.p.A.) [26,27] which has been already applied to a similar radial impeller [28]. The clearance between the impeller and the casing has been set to 1 mm.

2.2. Pumping Performance

The tests are performed at room temperature (T ≅ 20 °C), while the inlet pressure is set well above the vapor pressure in order to avoid cavitation inception within the pump. Referring to the pump geometry shown in Table 3, the parameters used to characterize the performance are summarized in Table 4.

2.3. Cavitating Performance

Steady-state cavitating experiments are performed by maintaining a fixed inlet pressure condition during the tests. Continuous cavitating experiments are performed by linearly lowering the inlet pressure in order to reach different cavitating regimes; the pressure decrease rate is approximately 0.003 bar/s. Signals are acquired at 5000 samples per second and each continuous experiment is divided into sub-blocks of 5 s in order to allow a quasi-steady hypothesis as experimentally demonstrated in previous works [17,29] and highlighted in the following results. Per each sub-block the operating conditions and relevant phenomena do not change significantly. The correspondent frequency resolution ( Δ f ) is 0.2 Hz. Moreover, each sub-block is superposed with the following one by 50% in order to catch all the phenomena. The main parameters used to characterize the cavitating performance are in Table 5.

2.4. Flow Instabilities and Unsteady Fluid Forces

The pressure fluctuations of the flow through the pump triggered or not by cavitation have been investigated during the cavitating experiments. Cervone et al. [16], and Pace et al. [29] report the procedures already exploited by the authors and summarized here. As for the cavitating performance evaluation, the signals, acquired by means of the piezoelectric transducers located on the casing of the pump and the diffuser, are divided into sub-blocks. Per each sub-block, the signals are fast-Fourier transformed in order to obtain the frequency content of oscillating phenomena. At each frequency of interest, the analysis of the cross-spectrum phase of signals from transducers placed at the same axial location but at different azimuthal positions allows for understanding the physical nature of the phenomenon. In fact, oscillating axial phenomena involve signals characterized by zero-phase cross-spectrum. On the other hand, rotating phenomena lead to cross-spectrum phases ( φ ) proportional to the circumferential spacing Δ ϑ between transducers and proportional to the phenomenon number of lobes ( n ) . For instance, for a generic cross-spectrum from the transducer 1 to the transducer 2, it would be
n = | φ 12 / ( ϑ 2 ϑ 1 ) | ,
where negative values of φ 12 imply a phenomenon moving from transducer 1 to transducer 2. Aliasing in the azimuthal direction is eliminated by comparison of the cross-correlations of pressure signals from transducers with different angular separations. The analysis is validated based on the value of the coherence function between signals obtained from two generic transducers ( x , y ) γ x y = | S x y ( f ) | 2 / [ S x x ( f ) S y y ( f ) ] , where S is either the cross-spectrum or the auto-spectrum of the referred transducers ( x , y ) in the subscripts (see [30] for further information). In the present work, only phenomena with values of γ x y greater than 0.95 have been considered.
In order to understand the frequency content of the signals recorded through the dynamometer, it is useful to consider and compare the effect of different types of unsteady forces on the dynamometer rotating frame ( x y z ) rotating with speed Ω and on the fixed one ( XYZ ) , where z and Z are directed along the rotating axis. Therefore, the force evaluated in the two reference frames has the same z -component as well as the same total force acting on the plane X Y ( x y ) .
Considering Figure 5, the schematic proposed involves a fluctuating phenomenon generating a force with constant direction ( θ 0 ) with respect to the fixed reference frame, a mean value F 0 and an oscillating amplitude: F 1 = F 0 + F 1 ˜ cos ( ω 2 t ) . Furthermore, the schematic considers a constant amplitude force F 2 rotating with respect to the fixed reference frame with an angular velocity ω 1 . The force components in the rotating frame and in the fixed one are
F X = F 0 cos ( θ 0 ) + F ˜ 1 cos ( ω 2 t ) cos ( θ 0 ) + F 2 cos ( ω 1 t ) ,
F Y = F 0 sin ( θ 0 ) + F ˜ 1 cos ( ω 2 t ) sin ( θ 0 ) + F 2 sin ( ω 1 t ) ,
F x = F ˜ 1 2 cos [ ( ω 2 Ω ) t + θ 0 ] + F ˜ 1 2 cos [ ( ω 2 + Ω ) t θ 0 ] + F 2 cos [ ( ω 1 Ω ) t ] + F 0 cos ( Ω t θ 0 ) ,
F y = F ˜ 1 2 sin [ ( ω 2 Ω ) t + θ 0 ] F ˜ 1 2 sin [ ( ω 2 + Ω ) t θ 0 ] + F 2 sin [ ( ω 1 Ω ) t ] F 0 sin ( Ω t θ 0 ) .
Therefore, a frequency analysis of the force components would show two frequencies (ω1/2π, ω2/2π) in the fixed reference frame. On the other hand, the same analysis in the rotating frame would show multiple frequencies shifted from the above-mentioned ones by ±Ω/2π.

3. Results and Discussion

3.1. Pumping Performance

Figure 6 shows the pumping performance obtained at three different rotating speeds. The data overlapping confirms the Reynolds independence of the tests (fully turbulent flow).

3.2. Cavitating Performance

Figure 7 reports the VAMPUFF pump cavitating performance obtained at different flow coefficients Φ . Per each flow coefficient, solid markers with black borders show results obtained from steady-state experiments, while blank markers report continuous experiments results per each sub-block. The overlapping of the data from the two different procedures confirms the goodness of the continuous data and the steady-state hypothesis which can reliably represent the cavitating performance. The rotational speed Ω 1750 rpm (183.3 rad/s, 29.2 Hz) prevent the exceeding of the maximum force measure capability of the dynamometer, allowing at the same time to obtain relevant cavitation phenomena. Other relevant parameters are reported in Table 6.

3.3. Flow Instabilities and Unsteady Fluid Forces

The systematic evaluation of the flow instabilities during the cavitating experiments highlighted the presence of rotating and axial phenomena for different flowrates, as briefly summarized in Table 7. In the table, at the same flowrate, oscillating phenomena showing analogous characteristics (rotating with 1 lobe, axial, etc.) are grouped together regardless of the operating regimes ( σ , frequency).
The presence of cavitation instabilities leads to unwanted forces on the shaft even at Φ = Φ D . However, at design condition the intensity of the detected oscillating forces is minor than at lower flowrates. For this reason, the present study only focuses on Φ = 0.108 ( Φ / Φ D = 0.9 ) in order to show some of the typical pressure-force spectra relationships found at all the tested regimes.
Figure 8 reports the cavitating performance at a nominal Φ N = 0.108 as well as the flow coefficient behavior during the inlet pressure decay. At the end of the experiment, the massive presence of cavitation leads to the breakdown which, in turn, leads to the flow coefficient drop.
The flow instabilities illustrated in this paper are basically connected with the presence of cavitation. According to [31], cavitation inception usually starts in the tip vortex generated at the blade inlet. Therefore, it is of no surprise that the major part of the found phenomena is clearly visible by the upstream pressure transducers, becoming less visible while moving downstream. On the other hand, axial phenomena effectively propagate from the upstream to the downstream and vice versa, becoming clearly visible in both stations. In the figures, Ω is intended as the rotating frequency.
Figure 9 reports the energy frequency content of three pressure transducers as representative of the three different stream locations: upstream, midstream, and downstream (according to Figure 2 and Figure 3) versus the cavitation number σ . The frequency energy content ( E f ) is directly related to the amplitude ( A f ) of the acting oscillating phenomenon as follows:
E f = A f 2 4 N s f s ,
where N s is the number of samples considered for the fast Fourier transform (FFT) while f s is the sampling frequency exploited during the experiment.
While the blade passage frequency is clearly visible at 6 Ω per each station, the other relevant phenomena are generally of major interest at a single station. In order to understand the physical nature of such oscillating phenomena, Figure 10, Figure 11 and Figure 12 report the phase of the cross-spectrum of the upstream, midstream, and downstream pressure transducers, respectively. The figures show only the phenomena of major interest characterized by an amplitude A f 100 Pa and a coherence γ x y 0.95 , focusing on the range of f / Ω = [ 0 , 4 ] where interesting phenomena have been found. The phase values allow for understanding the nature of the phenomenon [16,17]. The main outcomes of this analysis are summarized in Table 8.
Figure 13 shows the frequency energy content of the force components measured by the dynamometer with a force amplitude oscillation greater than 0.5 N. In particular, the figure reports the force components independent of the chosen reference frame (i.e., rotating or fixed), as described previously.
RC2 and RC3 do not generate relevant effects on the forces sensed by the dynamometer, therefore they won’t be further discussed. On the contrary, some of the instabilities reported in Table 8 clearly lead to unwanted oscillations of the force acting on the shaft. In particular, the two axial instabilities A1 and A2 generate oscillations of F z for corresponding values of σ f , while they are not visible at all on a plane perpendicular to the rotational axis ( XY ,   x y ) .
Like A1 and A2, the pressure distribution connected with the rotating cavitation-induced instability RC1 generates a fluctuating component on the rotational axis.
However, RC1 also leads to a fluctuating component on the plane XY   ( x y ) at the same operating regimes ( σ ) and frequencies, and with an intensity directly connected to the corresponding energy value.
In order to understand the RC1 effects on the plane XY   ( x y ) , it is useful to analyze the frequency content of the force component F X and F x defined in the fixed frame and in the rotating one, respectively (Figure 14). Let’s consider the schematic proposed in Figure 5 and in the following. The frequency energy content of F x y is influenced only by the amplitude and the acting frequency of F 1 ˜ . The pressure distribution due to the presence of cavitation in the form of RC1 generates a rotating force imbalance with a rotational velocity corresponding to the rotational velocity of the phenomenon itself, which is given by the pressure transducer analysis. Moreover, when this pressure distribution interacts with the (static) volute tongue at the impeller exit, it may lead to an oscillating unbalanced force at the same frequency as the phenomenon itself. The above considerations may be summarized as ω 2 = ω 1 = 2 π f R C 1 , therefore
F X = F 0 cos ( θ 0 ) + [ F ˜ 1 cos ( θ 0 ) + F 2 ] cos ( ω 1 t ) ,
which is coherent with the frequency content in Figure 14. Moreover, the corresponding force on the rotating frame is given by
F x = F ˜ 1 2 cos [ ( ω 1 Ω ) t + θ 0 ] + F 2 cos [ ( ω 1 Ω ) t ] + F ˜ 1 2 cos [ ( ω 1 + Ω ) t θ 0 ] + F 0 cos ( Ω t θ 0 ) ,
where there are only three acting frequencies (as shown in Figure 14, center):
  • ( ω 1 Ω ) / 2 π , whose intensity depends on the combination of F ˜ 1 , F 2 , θ 0 ;
  • ( ω 1 + Ω ) / 2 π ;
  • Ω / 2 π .

4. Conclusions and Future Work

The experimental characterization of pressure fluctuations on a centrifugal pump for various operating regimes has been here presented and discussed. The paper focuses on the unsteady forces acting on the driving shaft which are triggered by such fluctuations. The contemporary measurement of the unsteady forces and the pressure field on centrifugal pumps constitutes an innovative feature of the present study. Although the paper presents data coming from a specific regime in terms of flow coefficient ( Φ = 0.9 Φ D ), similar behaviors have been detected for analogous phenomena observed during other regimes. In particular, the presence of rotating phenomena may generate fluctuations of the force in the plane perpendicular to the rotational axis as well as on the rotational axis itself. On the perpendicular plane, the frequency content associated with the generated force suggests that the force can be approximated by the sum of a component with a fixed direction and fluctuating intensity together with a purely rotating component. The frequency of both components is the same as for the source phenomenon. A key role in such behavior is most likely played by the single-tongue volute, whose intrinsic asymmetry interacts with the rotating phenomena, leading to force amplitude fluctuations. At Φ = Φ D the outlet flow has a nominal zero incidence angle with the tongue, which is most likely the reason why the force triggered by cavitation instabilities is also reduced. This consideration confirms the goodness of the pump design approach described in [26], also from a rotordynamic point of view when cavitation instabilities arise. However, experimental campaigns are needed to better understand the interaction between the presence of cavitation instabilities and forces in presence of different shaped volutes and diffusers. Furthermore, the paper highlights that special attention should be paid when the pump design includes any stator (e.g., a vaned diffuser) especially when the operating conditions may include cavitation.
On the other hand, axial phenomena lead to fluctuations of the force component directed along the rotational axis only.
Interestingly, not all the phenomena measured by the pressure transducers lead to relevant fluctuations of the forces sensed by the dynamometer. Thus, different sensors may underline specific characteristics and effects of the fluid fluctuations. For instance, a strong fluctuation involving a small portion of the blade channels does not necessary lead to a strong fluctuation of the forces on the shaft. On the other hand, phenomena with smaller fluctuation intensities but characterized by large area of interest may involve stronger oscillations of the forces on the shaft which are more dangerous from a structural point of view. Therefore, it is advisable to fully characterize the fluid instabilities and their effects by means of different sensors under the foreseen operating regimes in order to correctly design the system.
The systematic analysis of the signals from the pressure and force sensors with different methods, such as wavelet decomposition, may give some useful insight into the inception and end of each flow pattern although it does not improve the understanding of the fully-developed flow instabilities itself [32]. Therefore, future activities may include such analysis especially when the pump configuration allows for optical access (e.g., an inducer setup at the CPRTF) in order to exactly match the sensors measurements to the cavitation structures.
To date, it is not clear which eventual unstable forces would arise with an inducer ahead of a centrifugal stage (typical in space application) under the presence of cavitation instabilities. However, the present results suggest that the massive presence of such instabilities for an inducer (if compared to a centrifugal stage) may lead to possible dangerous conditions when the inducer is coupled with static elements typical of centrifugal stages (e.g., the volute tongue). Further investigations are needed in order to better clarify these aspects.

Nomenclature

Latin
A Axial
A f Fluctuation amplitude
E f Fluctuation energy
F ˜ Force intensity fluctuation
F Force
F S Full scale
N Blade number
Q Volume flow rate
R C Rotating cavitation
R e Reynolds
S Power spectrum density
T Temperature
f Frequency
n Flow instability lobes
p Pressure
r Radius
z Axial length, rotating axis
Greek
Δ Variation
Φ Flow coefficient
Ψ Head coefficient
Ω Rotating speed, rotating frequency
γ x y Coherence between sensors x and y
ν Fluid kinematic viscosity
ρ Fluid density
σ Cavitation number
ω Force angular frequency
φ Cross-spectrum phase
θ Force initial phase
ϑ Azimuthal direction
Subscripts
D Design
N Nominal
T Total
c Chamber
i n Inlet
v Vapor
x y z Rotating frame
X Y Z Fixed frame

Author Contributions

D.V., A.P. and R.H. designed the experiment. D.V., G.P. and R.H. carried out the experimental activities. D.V., G.P., A.P. and R.H. analyzed the data. L.d. supported and inspired the entire project. D.V. wrote the paper. All the authors contributed to writing-review and editing.

Acknowledgments

The present work has been possible thanks to the supports of the European Space Agency along several years. The authors would like to express their great gratitude to Giorgio Saccoccia and to Gianni Pellegrini. Special gratitude goes to Lucio Torre that has been a great mainstay of the whole Chemical Propulsion Team for many years and to Kai Su.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bhattacharyya, A.; Acosta, A.J.; Brennen, C.E.; Caughey, T.K. Rotordynamic Forces in Cavitating Inducers. J. Fluids Eng. 1997, 119, 768. [Google Scholar] [CrossRef]
  2. Bhattacharyya, A. Internal Flows and Force Matrices in Axial Flow Inducers. Ph.D. Thesis, California Institute of Technology, Pasadena, CA, USA, 1994. [Google Scholar]
  3. Franz, R.J. Experimental Investigation of the Effect of Cavitation on the Rotordynamic Forces on a Whirling Centrifugal Pump Impeller. Ph.D. Thesis, California Institute of Technology, Pasadena, CA, USA, 1989. [Google Scholar]
  4. Rosenmann, W. Experimental investigations of hydrodynamically induced shaft forces with a three bladed inducer. In Proceedings of the ASME Symposium on Cavitation in Fluid Machinery, Chicago, IL, USA, 7–11 November 1965; pp. 172–195. [Google Scholar]
  5. Ehrich, F.; Childs, D. Self-Excited Vibrations in High Performance Turbomachinery. Mech. Eng. 1984, 106, 66–79. [Google Scholar]
  6. Hergt, P.; Krieger, P. Radial Forces in Centrifugal Pumps with Guide Vanes. Proc. Inst. Mech. Eng. 1969, 184, 101–107. [Google Scholar] [CrossRef]
  7. Suzuki, T.; Prunières, R.; Horiguchi, H.; Tsukiya, T.; Taenaka, Y.; Tsujimoto, Y. Measurements of Rotordynamic Forces on an Artificial Heart Pump Impeller. J. Fluids Eng. 2007, 129, 1422–1427. [Google Scholar] [CrossRef]
  8. Brennen, C.E.; Franz, R.; Arndt, N. Effects of Cavitation on Rotordynamic Force Matrices. In Proceedings of the 3rd Conference on Advanced Earth-to-Orbit Propulsion Technology; Richmond, R.J., Wu, S.T., Eds.; NASA: Washington, DC, USA, 1988; pp. 227–239. [Google Scholar]
  9. Hashimoto, T.; Yoshida, M.; Watanabe, M. Experimental Study on Rotating Cavitation of Rocket Propellant Pump Inducers. J. Propuls. Power 1997, 13, 488–494. [Google Scholar] [CrossRef]
  10. Tsujimoto, Y.; Yoshida, Y.; Maekawa, Y.; Watanabe, S.; Hashimoto, T. Observations of Oscillating Cavitation of an Inducer. J. Fluids Eng. 1997, 119, 775–781. [Google Scholar] [CrossRef]
  11. Zoladz, T. Observations on rotating cavitation and cavitation surge from the development of the Fastrac engine turbopump. In Proceedings of the 36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Las Vegas, NV, USA, 24–28 July 2000. [Google Scholar]
  12. Tsujimoto, Y.; Semenov, Y. New Types of Cavitation Instabilities in Inducers. In Proceedings of the 4th International Conference on Launcher Technology “Space Launcher Liquid Propulsion”, Liege, Belgium, 3–6 December 2002. [Google Scholar]
  13. Subbaraman, M.; Patton, M. Suppressing Higher-Order Cavitation Phenomena in Axial Inducers. In Proceedings of the 6th International Symposium on Cavitation, Wageningen, The Netherlands, 11–15 September 2006. [Google Scholar]
  14. Kamijo, K.; Yoshida, M.; Tsujimoto, Y. Hydraulic and mechanical performance of LE-7 LOX pump inducer. J. Propuls. Power 1993, 9, 819–826. [Google Scholar] [CrossRef]
  15. Cervone, A.; Bramanti, C.; Rapposelli, E.; Torre, L.; d’Agostino, L. Experimental Characterization of Cavitation Instabilities in a Two-Bladed Axial Inducer. J. Propuls. Power 2006, 22, 1389–1395. [Google Scholar] [CrossRef]
  16. Cervone, A.; Torre, L.; Pasini, A.; d’Agostino, L. Cavitation and Flow Instabilities in a 3- Bladed Axial Inducer Designed by Means of a Reduced Order Analytical Model. In Proceedings of the 7th International Symposium on Cavitation, Ann Arbor, MI, USA, 17–22 August 2009. [Google Scholar]
  17. Torre, L.; Pace, G.; Miloro, P.; Pasini, A.; Cervone, A.; d’Agostino, L. Flow Instabilities on a Three Bladed Axial Inducer at Variable Tip Clearance. In Proceedings of the 13th International Symposium on Transport Phenomena and Dynamics of Rotating Machinery, Honolulu, HI, USA, 4–9 April 2010. [Google Scholar]
  18. Friedrichs, J.; Kosyna, G. Rotating Cavitation in a Centrifugal Pump Impeller of Low Specific Speed. J. Fluids Eng. 2002, 124, 356–362. [Google Scholar] [CrossRef]
  19. Yamamoto, K. Instability in a Cavitating Centrifugal Pump. JSME Int. J. Ser 2 Fluids Eng. Heat Transf. Power Combust. Thermophys. Prop. 1991, 34, 9–17. [Google Scholar] [CrossRef]
  20. Valentini, D.; Pace, G.; Torre, L.; Pasini, A.; d’Agostino, L. Influences of the Operating Conditions on the Rotordynamic Forces Acting on a Three-Bladed Inducer Under Forced Whirl Motion. J. Fluids Eng. 2015, 137, 071304–071304-10. [Google Scholar] [CrossRef]
  21. Valentini, D.; Pace, G.; Pasini, A.; Torre, L.; Hadavandi, R.; d’Agostino, L. Fluid-induced rotordynamic forces on a whirling centrifugal pump. Eur. J. Mech. B/Fluids 2017, 61, 336–345. [Google Scholar] [CrossRef]
  22. Valentini, D.; Pace, G.; Torre, L.; Pasini, A.; d’Agostino, L. Pumping and Suction Performance of a Whirling Inducer. In Proceedings of the Space Propulsion 2014, Cologne, Germany, 19–22 May 2014. [Google Scholar]
  23. Pace, G.; Pasini, A.; Torre, L.; Valentini, D.; d’Agostino, L. Cavitating Pump Rotordynamic Test Facility at ALTA S.p.A.: Upgraded Capabilities of a Unique Test Rig. In Space Propulsion Conference 2012; ESA: Bordeaux, France, May 2012. [Google Scholar]
  24. Pasini, A. Pumping Performance Similarity, Cavitation-Induced Instabilities and Fluid-Induced Rotordynamic Forces in Tapered Inducers. Ph.D. Thesis, University of Pisa, Pisa, Italy, 2010. [Google Scholar]
  25. d’Agostino, L. Comparison of rotordynamic fluid forces in axial inducers and centrifugal turbopump impellers. IOP Conf. Ser. Mater. Sci. Eng. 2016, 129, 012001. [Google Scholar] [CrossRef]
  26. d’Agostino, L.; Torre, L.; Pasini, A.; Baccarella, D.; Cervone, A.; Milani, A. A Reduced Order Model for Preliminary Design and Performance Prediction of Tapered Inducers: Comparison with Numerical Simulations. In Proceedings of the 44th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Hartford, CT, USA, 21–23 July 2008. [Google Scholar]
  27. d’Agostino, L.; Torre, L.; Pasini, A.; Cervone, A. On the Preliminary Design and Noncavitating Performance Prediction of Tapered Axial Inducers. J. Fluids Eng. 2008, 130, 111303:1–111303:8. [Google Scholar] [CrossRef]
  28. Valentini, D.; Pasini, A.; Pace, G.; Torre, L.; d’Agostino, L. Experimental Validation of a Reduced Order for Radial Turbopump Design. In Proceedings of the 49th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, San Jose, CA, USA, 14–17 July 2013. [Google Scholar]
  29. Pace, G.; Valentini, D.; Pasini, A.; Torre, L.; Fu, Y.; d’Agostino, L. Effects of Geometry on Flow Instabilities of Different Three-Bladed Axial Inducers. J. Fluids Eng. 2015, 137, 041304. [Google Scholar] [CrossRef]
  30. Bendat, J.S.; Piersol, A.G. Random Data: Analysis and Measurement Procedures; Wiley Series in Probability and Statistics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2010; ISBN 978-1-118-03242-8. [Google Scholar]
  31. Brennen, C.E. Hydrodynamics of Pumps; Concepts ETI, Inc. and Oxford University Press: Oxford, UK, 1994. [Google Scholar]
  32. Coutier-Delgosha, O.; Dazin, A.; Caignaert, G.; Bois, G. Analysis of Cavitation Instabilities in a Four-Blade Inducer. Int. J. Rotating Mach. 2012, 2012. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The cavitating pump rotordynamic test facility (CPRTF) layout in SITAEL S.p.A. Abbreviations: OPT1, outlet pressure tap 1; OPT2, outlet pressure tap 2; IPT1, inlet pressure tap 1; IPT2, inlet pressure tap 12.
Figure 1. The cavitating pump rotordynamic test facility (CPRTF) layout in SITAEL S.p.A. Abbreviations: OPT1, outlet pressure tap 1; OPT2, outlet pressure tap 2; IPT1, inlet pressure tap 1; IPT2, inlet pressure tap 12.
Ijtpp 03 00013 g001
Figure 2. VAMPUFF pump test setup.
Figure 2. VAMPUFF pump test setup.
Ijtpp 03 00013 g002
Figure 3. Piezoelectric pressure transducers positions.
Figure 3. Piezoelectric pressure transducers positions.
Ijtpp 03 00013 g003
Figure 4. VAMPUFF radial impeller.
Figure 4. VAMPUFF radial impeller.
Ijtpp 03 00013 g004
Figure 5. Schematic of a generic unsteady force acting on the shaft with respect to the rotating frame ( x y ) and the fixed one ( X Y ). F 1 , force with a fixed direction θ 0 and oscillating amplitude; F 2 , rotating force with velocity ω 1 ; Ω , rotating frame speed.
Figure 5. Schematic of a generic unsteady force acting on the shaft with respect to the rotating frame ( x y ) and the fixed one ( X Y ). F 1 , force with a fixed direction θ 0 and oscillating amplitude; F 2 , rotating force with velocity ω 1 ; Ω , rotating frame speed.
Ijtpp 03 00013 g005
Figure 6. VAMPUFF total head ( Ψ T ) versus flow coefficient ( Φ ) at different rotating speed ( Ω ), T ≅ 20 °C.
Figure 6. VAMPUFF total head ( Ψ T ) versus flow coefficient ( Φ ) at different rotating speed ( Ω ), T ≅ 20 °C.
Ijtpp 03 00013 g006
Figure 7. VAMPUFF cavitating performance as head coefficient ( Ψ ) versus cavitation number ( σ ) at different flow coefficients ( Φ ) (T = 20 °C; Ω = 1750 rpm).
Figure 7. VAMPUFF cavitating performance as head coefficient ( Ψ ) versus cavitation number ( σ ) at different flow coefficients ( Φ ) (T = 20 °C; Ω = 1750 rpm).
Ijtpp 03 00013 g007
Figure 8. Cavitating performance (head coefficient Ψ versus flow coefficient Φ ) and flow coefficient evolution, nominal Φ N = 0.108 (90% Φ D ).
Figure 8. Cavitating performance (head coefficient Ψ versus flow coefficient Φ ) and flow coefficient evolution, nominal Φ N = 0.108 (90% Φ D ).
Ijtpp 03 00013 g008
Figure 9. Frequency energy content [Pa2∙s] of pressure transducers placed upstream (left), midstream (center), and downstream (right). The red circled identified phenomena are reported in Table 8. Abbreviations: UP, upstream; MID, midstream; DWN, downstream; RC, rotating cavitation; A, axial.
Figure 9. Frequency energy content [Pa2∙s] of pressure transducers placed upstream (left), midstream (center), and downstream (right). The red circled identified phenomena are reported in Table 8. Abbreviations: UP, upstream; MID, midstream; DWN, downstream; RC, rotating cavitation; A, axial.
Ijtpp 03 00013 g009
Figure 10. Cross-spectrum phase ( φ ) of the upstream pressure transducers,   Φ = 0.108 (90% Φ D ). Red circled phenomena are summarized in Table 8.
Figure 10. Cross-spectrum phase ( φ ) of the upstream pressure transducers,   Φ = 0.108 (90% Φ D ). Red circled phenomena are summarized in Table 8.
Ijtpp 03 00013 g010
Figure 11. Cross-spectrum phase ( φ ) of the midstream pressure transducers, Φ = 0.108 (90% Φ D ). circled phenomena are summarized in Table 8.
Figure 11. Cross-spectrum phase ( φ ) of the midstream pressure transducers, Φ = 0.108 (90% Φ D ). circled phenomena are summarized in Table 8.
Ijtpp 03 00013 g011
Figure 12. Cross-spectrum phase ( φ ) of the downstream pressure transducers, Φ = 0.108 (90% Φ D ). Red circled phenomena are summarized in Table 8.
Figure 12. Cross-spectrum phase ( φ ) of the downstream pressure transducers, Φ = 0.108 (90% Φ D ). Red circled phenomena are summarized in Table 8.
Ijtpp 03 00013 g012
Figure 13. Frequency energy content [N2∙s] of the force along the axis ( F z ) , its perpendicular component ( F x y ) , and the total force ( F x y z ) acting on the dynamometer. Red circled phenomena are summarized in Table 8.
Figure 13. Frequency energy content [N2∙s] of the force along the axis ( F z ) , its perpendicular component ( F x y ) , and the total force ( F x y z ) acting on the dynamometer. Red circled phenomena are summarized in Table 8.
Ijtpp 03 00013 g013
Figure 14. Frequency energy content [N2∙s] of the force component along the X-axis in the absolute fixed reference frame ( F X , left), along the x-axis in the rotating reference frame ( F x ) and its perpendicular component ( F x y ) acting on the dynamometer. Red circled phenomena are summarized in Table 8.
Figure 14. Frequency energy content [N2∙s] of the force component along the X-axis in the absolute fixed reference frame ( F X , left), along the x-axis in the rotating reference frame ( F x ) and its perpendicular component ( F x y ) acting on the dynamometer. Red circled phenomena are summarized in Table 8.
Ijtpp 03 00013 g014
Table 1. Force measurement errors with a confidence level of 95%.
Table 1. Force measurement errors with a confidence level of 95%.
Force ComponentError
(95% Confidence)
F x ±1.19 N
F y ±1.49 N
F z ±1.80 N
Table 2. Pressure taps exploited for the present campaign.
Table 2. Pressure taps exploited for the present campaign.
PositionTap
Upstream1, 2, 3, 4
Midstream1, 7
Downstream1, 2, 3, 7
Table 3. Most relevant geometrical and operational parameters of the centrifugal pump.
Table 3. Most relevant geometrical and operational parameters of the centrifugal pump.
ParameterUnitSymbolValue
Design flow coefficient @r2[--] Φ D   0.120
Number of blades[--] N 6
Outlet radius[mm] r 2 105.00
Inlet tip radius[mm] r T 1 55.50
Inlet hub radius[mm] r H 1 31.42
Axial length[mm] z H 2 67.00
Inlet tip blade angle[deg] γ T 1 47.70
Inlet backsweep angle[deg] χ 1 0
Diffuser outlet radius[mm] r 3 126
Design Rotating Speed[rad/s]Ω157.1
[rpm]1500
Design volume flowrate[m3/s] Q d e s 0.022
Mean blade height[mm] h m 18.84
Tip solidity[--] σ T 2.64
Incidence tip angle @ design[deg] α 18.90
Outlet tip blade mean angle[deg] γ T 2 89.94
Outlet tip backsweep angle[deg] χ 2 64.00
Table 4. Pumping performance parameters.
Table 4. Pumping performance parameters.
Total Head Coefficient Ψ T = Δ p T / ρ Ω 2 r 2 2
Flow Coefficient Φ = Q / Ω r 2 3
Reynolds Number R e = 2 Ω r 2 2 / ν
Table 5. Cavitating performance parameters.
Table 5. Cavitating performance parameters.
Head Coefficient Ψ = Δ p / ρ Ω 2 r 2 2
Cavitation Number σ = ( p i n p v ) / 0.5 ρ Ω 2 r 2 2
Table 6. Relevant parameters for the cavitating experiments.
Table 6. Relevant parameters for the cavitating experiments.
Sampling Frequency ( f s ) 5000 Hz
Temperature20 °C
Experiment Duration (cont.)240 s
Sub-Block Duration5 s
Table 7. Identified phenomena at different flowrates.
Table 7. Identified phenomena at different flowrates.
Φ(Φ/ΦD)Identified Phenomena
0.096 (0.8)Rotating (1 lobe and 2 lobes); axial
0.108 (0.9)Rotating (1 lobe and 2 lobes); axial
0.120 (1.0)Rotating (1 lobe); axial
0.132 (1.1)Rotating (1 lobe); axial
0.144 (1.2)Axial
Table 8. Summary of the found instabilities at Φ N = 0.108. Abbreviations: RC, rotating cavitation; A, axial.
Table 8. Summary of the found instabilities at Φ N = 0.108. Abbreviations: RC, rotating cavitation; A, axial.
IDFrequency Range
f
(f [Hz])
σ RangeCharacteristicsMajor Station
A10.1–0.38
(3–11)
0.16–0.12AxialUp
RC12.57–2.74
(75–80)
0.2–0.16Rotating 1 lobeUp
RC22.33–2.85
(68–83)
0.16–0.11Rotating 2 lobesMid
A23.22–3.81
(94–111)
0.16–0.10AxialDown
RC33.26–3.46
(95–101)
0.19–0.14Rotating 1 lobeMid

Share and Cite

MDPI and ACS Style

Valentini, D.; Pace, G.; Pasini, A.; Hadavandi, R.; D’Agostino, L. Experimental Characterization of Unsteady Forces Triggered by Cavitation on a Centrifugal Pump. Int. J. Turbomach. Propuls. Power 2018, 3, 13. https://doi.org/10.3390/ijtpp3020013

AMA Style

Valentini D, Pace G, Pasini A, Hadavandi R, D’Agostino L. Experimental Characterization of Unsteady Forces Triggered by Cavitation on a Centrifugal Pump. International Journal of Turbomachinery, Propulsion and Power. 2018; 3(2):13. https://doi.org/10.3390/ijtpp3020013

Chicago/Turabian Style

Valentini, Dario, Giovanni Pace, Angelo Pasini, Ruzbeh Hadavandi, and Luca D’Agostino. 2018. "Experimental Characterization of Unsteady Forces Triggered by Cavitation on a Centrifugal Pump" International Journal of Turbomachinery, Propulsion and Power 3, no. 2: 13. https://doi.org/10.3390/ijtpp3020013

APA Style

Valentini, D., Pace, G., Pasini, A., Hadavandi, R., & D’Agostino, L. (2018). Experimental Characterization of Unsteady Forces Triggered by Cavitation on a Centrifugal Pump. International Journal of Turbomachinery, Propulsion and Power, 3(2), 13. https://doi.org/10.3390/ijtpp3020013

Article Metrics

Back to TopTop