International Journal of Rotating Machinery

Volume 2012 (2012), Article ID 213907, 11 pages

doi:10.1155/2012/213907

Research Article

Analysis of Cavitation Instabilities in a Four-Blade Inducer

O. Coutier-Delgosha, A. Dazin, G. Caignaert, and G. Bois

Arts et Metiers ParisTech/LML Laboratory, 8 boulevard Louis XIV, 59046 Lille Cedex, France

Received 26 December 2011; Accepted 18 April 2012

Academic Editor: Meinhard Taher Schobeiri

Copyright © 2012 O. Coutier-Delgosha et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The cavitating behavior of a four-blade inducer tested in the LML laboratory large test facility is considered in the present paper. Experimental investigations based on unsteady pressure measurements and records from a six-component balance mounted on the inducer shaft are performed. Spectral analysis of the signals enables to detect several characteristic frequencies related to unbalanced two-phase flow patterns. The objective of the present paper is the understanding of the physical phenomena associated to these frequencies. Therefore, wavelet decomposition, flow visualizations, and direct analysis of the high-frequency force, moment, and pressure signals are applied. Results at nominal flow rate only are considered. Not only classical unbalanced cavitation patterns, but also unexpected flow organizations are discussed.

1. Introduction

Rocket engine turbopump inducers are designed to operate in cavitating conditions, because of the low pressure of LH2 and LOX in the tanks. However, development of cavitation must be controlled to avoid any reduction of the inducer performance. Attention is also focused on instabilities associated with cavitation, since they are responsible for unbalanced flow pattern that compromise the rotor equilibrium. Self-oscillation of the sheet cavities on the blades may also generate problematic pressure fluctuations at the inducer outlet.

The understanding of these different unsteady flow behaviors has not been completed, yet, although significant progress has been achieved in the last decades. Experimental investigations of the behavior of various inducers, characterized by a blade number varying usually between 2 and 5, have been conducted for example in Pisa University (Italy) in collaboration with Avio, in Arts et Metiers ParisTech and CREMHyG (France) in collaboration with SNECMA Moteurs, and in Osaka University (Japan) by Tsujimoto and colleagues in cooperation with the JAXA program. In these different studies, several recurrent unstable behaviors have been identified: (i) rotating cavitation is characterized by a nonsymmetrical two-phase flow pattern, which rotates faster or slower than the inducer, leading to so-called super-synchronous or subsynchronous regimes, respectively. Rotation speed in the inducer rotating frame is generally low, leading in the absolute frame to frequencies close to the inducer rotation one. Such phenomenon, which leads to significant radial forces, has been detected in nearly all tested inducers [1–3] (Tsujimoto et al. 1997). (ii) A surge mode oscillation characterized by a global pulsation of the vapor areas at low frequency: such oscillations have been investigated in simpler configurations of 2D Venturi profiles or 2D foil sections in cavitation tunnels and analyzed in details by many authors [4–7]. (iii) Rotating stall in backflow vortices: such behavior has been identified by Tsujimoto and colleagues [8] in the case of cavitation clouds extending upstream from the inducer inlet. It basically consists in a slow rotation of this detached cavitation, at frequency that depends on the number of vortices.

The occurrence of these different flow patterns strongly depends on the cavitation development in the machine. A typical sketch in the case of a four-blade inducer is given for example in Figure 1.

Figure 1: Sketches of cavitation patterns and performance evolution as the cavitation number decreases in a four-blade inducer [9].

At cavitation inception, only a steady and balanced flow pattern with one short attached cavity on each blade is observed from flow visualizations. When cavitation parameter is slightly decreased, a steady and alternate cavitating configuration appears with alternatively one short and one long cavity. For a lower cavitation parameter, just above breakdown, rotating cavitation can be identified: unbalanced attached cavities are observed in the different channels, their distribution rotating faster than the inducer. Finally near the breakdown of the inducer, a steady and balanced flow pattern with fully developed cavitation is observed.

More recently, supplementary unstable flow patterns have been identified in particular configurations: (i) higher-order rotating cavitation has been reported first by Fujii et al. [10], from Osaka University, in the case of a three-blade inducer. Frequency related to this phenomenon is 5 in the absolute frame. (ii) Higher-order cavitation surge: this axial instability has been found by the Osaka Team at frequency 5 , and also by researchers from Pisa University at frequencies 4.4 and 6.6 in the case of a two-blade inducer. In both cases, the phenomenon was detected at high flow coefficient.

These instabilities induce some strong radial forces that may perturb the rotor balance, and important pressure fluctuations in the lines. They must be quantified and controlled to avoid any major effect on the global pump behavior. Their detection is usually based, in all aforementioned references, on pressure transducers mounted on the inducer shroud at several angles in the same section. The analysis of the phase difference between the transducers enables to distinguish axial instabilities (characterized by zero-phase difference) from rotating cavitation (where a phase difference equal to the angle difference is obtained).

As can be seen, the physical characterization of the flow instabilities related to cavitation is still weak, which explains that only a few attempts have been performed to explain inception and mechanisms of rotating cavitation [11, 12]. Yoshida and colleagues state that inception of unstable behaviours is mainly governed by the sheet cavity length . Sheet cavities longer than the cascade throat remain equal because they are fully constrained by the adjacent blade. In opposition, sheet cavities that do not reach the cascade throat () are characterized by a degree of freedom, which may lead to unstable cavitation such as synchronous nonsymmetrical flow patterns. However, no consensus is presently obtained concerning too the mechanisms that are responsible for the inception of unbalanced flow patterns as the ones that control the successive synchronous and nonsynchronous regimes that are usually observed in inducers.

Conversely, several studies have been conducted to investigate the effects of the inducer design on its unsteady behavior: it has been found, for example, that the shape of the hub [13], of the leading edge [14], (Bakir et al. 2002), and even of the casing [15] have significant influence on cavitation instabilities. The effects of unequal spacing of the blades have been also investigated by numerical simulations [16]. All inducers have been usually manufactured with one to 4 blades in the inlet section.

In the present paper, the cavitating behaviour of a four-blade inducer (Figure 2) is investigated. Experiments have been performed in the large test facility of ENSAM Lille, with the collaboration of the CNES (French Space Agency) and SNECMA Moteurs. Although a large range of mass flow rates, inlet pressures, and rotation speeds have been tested, attention is focused presently on a few hydrodynamic conditions only, at nominal mass flow rate and reference rotation speed . The successive flow instabilities, which are encountered when the inlet pressure is progressively decreased, are analysed in details, in order to improve their characterization and physical understanding. Original force measurements performed with a balance mounted on the inducer shaft, associated with unsteady pressure measurements on the inducer shroud, enable to obtain data both in the absolute and inducer rotating frame, which enriches the analysis of the results.

Figure 2: Inducer geometry.

Note that all results presented in the paper are dimensionless for reasons of confidentiality. So, the “+” indicates that each variable is where is a reference value which is not given.

2. Experimental Setup

2.1. Test Facility

The LML laboratory large test facility devoted to the study of axial pumps in cavitating conditions has been used for the experiments. This two-stairs facility is equipped with a 200 kW motor that enables to reach an inducer rotation speed of 6000 rpm. The device upstairs (Figure 3) is mainly composed of the inducer to be tested and a stator downstream from the inducer that conducts the flow towards the discharge pipe. The inlet and discharge pipes are connected to two tanks. A free surface is maintained in the upstream one, so that the pressure can be controlled in the installation, from about 0.1 bars up to 16 bars. A part of the volume flow rate is taken at the inducer outlet for the axial equilibrium device, and then reinjected in the upstream tank. Downstairs are located two large resorption tanks used for the extraction of dissolved gas, a variable head loss device that enables to control the flow rate, and the vacuum pump used to decrease the pressure in the inlet tank. A heat exchanger is also connected to the facility to control the water temperature.

Figure 3: General view of the LML large test facility.

2.2. Acquisition Device

The inducer test section is equipped with a large variety of sensors devoted either to the characterization of steady flow properties, or to the analysis of the flow unsteady fluctuations.

Concerning steady flow properties, the torque and the rotation speed N are measured with a Torquemaster TM213 torque meter, while a Rosemount AP6 pressure transducer (range 0–6 bars) is used for the inlet pressure . A Rosemount DP7E22 differential pressure transducer (range 0–16 bars) connected between and is used to obtain directly the elevation . The main flow rate and the recirculating one are measured with Endress + Hauser electromagnetic flow meters. The motor shaft rotation speed is also controlled with a photoelectric cell.

Temperature measurements are performed on the upstream and downstream bearings of the inducer shaft, to control their increase at high rotation speed, and the water temperature is measured to regulate the heat exchanger operation so that a 25°C temperature is maintained in the flow. The torque and the rotation speed are acquired at frequency 100 Hz during 10 s, while all other signals are recorded at frequency 40 Hz during 140 s. From these records, mean and RMS values of the signals are calculated for each investigated flow condition. Note that the inlet and outlet pressures are measured in the inlet and delivery pipes more than one meter upstream and downstream from the inducer, so the calculated pressure head includes the head losses between the two pressure taps, including the one in the stator and in the downstream 90° bends that can be seen in Figure 4.

Figure 4: Scheme of the inducer test section including the acquisition equipment.

Absolute uncertainties on the mean flow characteristics are 0.5% of for , 0.05% of for , 0.25% of for , 2% of the noncavitating torque value for , and 0.2% of the measurement range for the upstream pressure and . Precision regarding the cavitation number and the head, torque, and flow rate coefficients depend on the flow conditions and will be indicated in the figures hereafter.

Unsteady flow properties are investigated with nine Kistler 701A piezo-electric pressure transducers (Figure 4): six of them are located in the upstream and delivery pipes, two other ones ( and ) are located at the inlet of the inducer test section, and the last one is between the inducer and the stator. All transducers are mounted flush to the internal pipe wall. Four accelerometers are also installed on the inducer casing to measure the vibrations in the axial and radial directions. A six-component balance is mounted on the shaft to obtain the axial force , the radial forces and , the torque , and also the bending moments and . These 19 parameters are recorded simultaneously at frequency 2048 Hz with the LMS CADA-X code. An antialiasing filtering (Butterworth, 800 Hz, −50 dB/octave) is applied during acquisition.

Relative uncertainty on the pressure and vibration measurements is found to be close to 1% in noncavitating conditions and 2% in cavitating conditions. Precision on the forces and moments measured by the balance is estimated to 1% of the measurement ranges.

2.3. Experimental Process

Unsteady flow properties are investigated by varying continuously the inlet pressure from 2 bars down to the 25% drop of the inducer head. The duration of the pressure decrease is about 360 s, which may result in quasisteady flow conditions (this point is checked in [17]). During this process, the variable head loss device position is not modified, so the mass flow rate progressively decreases together with the inducer elevation. The present study focuses on the results obtained at nominal mass flow rate and reference rotation speed . However, data recorded at other rotation speeds are also used, because (i) visualizations have been performed at 0.6 and (ii) intensity of some fluctuations decreases with , so they can be detected only at high values of .

3. Methodology

Three methods have been applied for the analysis of the instabilities.

(i) Wavelet decomposition: this technique enables to identify on a signal the time of inception of the successive characteristic frequencies. It is based on a decomposition of the initial signal into secondary signals . Each signal is obtained from a typical shape called “wavelet” and a discrete frequency corresponding to the period of the wavelet. For example, the secondary signal at frequency is composed of a succession of wavelets whose period is 1/, and magnitude varies according to the consistence of the initial signal with the shape of this wavelet. The residual, that is, - is then decomposed into a secondary signal related to frequency , and a new residual that will be decomposed again and so on. The initial signal is decomposed according to several discrete frequencies, until a residual of low magnitude is obtained. If high amplitude is observed locally on a secondary signal , then it means that a fluctuation at frequency close to is included in the initial signal. In the present study, decomposition into discrete wavelet (Daubechies no. 5) is applied to the pressure and force signals.

(ii) Study of the temporal signals: this method aims, through analysis of the force and pressure unsteady signals, to improve the understanding of the inception and development of the instabilities. The main objective consists in identifying transient unbalanced flow patterns, flow conditions at inception and/or vanishing of the instabilities, and phenomena related to the detected characteristic frequencies. The analysis focuses particularly on the phase difference between and , or and (mounted in the same cross-section with 90° difference between them). Attention is also paid to the mean values of and (the mean value of these two signals in noncavitating conditions has been set to zero). In this analysis, low-pass filters have been applied in order to clearly distinguish some characteristic frequencies. Cut-off frequencies will be indicated systematically hereafter.

4. Wavelet Decomposition

The present inducer is characterized by a large number of nonsymmetrical flow patterns: so-called “synchronous” (unbalanced flow pattern rotating at inducer speed), “subsynchronous,” “super-synchronous” (unbalanced flow pattern rotating slower or faster than the inducer, resp.), and alternate (unbalanced flow pattern with identical opposite sheet cavities) regimes have been identified. An unknown fluctuation at frequency 0.9 in the absolute frame has been also detected.

These flow instabilities have been identified from the waterfall plots of , , , and , as indicated in Figure 5 for N = .

Figure 5: Waterfall plot for and (, ).

For example, super-synchronous rotating cavitation is clearly detected for . It is identified by the low frequency that is obtained in the rotating frame from the signal, which is correlated with the peaks that can be seen from the analysis of , at frequency higher than the one of the inducer rotation. This pattern basically consists of large-scale fluctuations of the radial force components and (Figure 6), at frequency that progressively decreases from a few tens of Hertz down to zero when is decreased. A 90° phase difference between and is observed, which suggests that these fluctuations may be associated to a rotating instability. Note that the values of the forces are not given in ordinate for confidentiality reasons: only the zero is indicated to facilitate the discussion.

Figure 6: Time signals of and during the super-synchronous regime (, , cut-off frequency 400 Hz).

For higher cavitation number (), a stable nonsymmetrical flow pattern is detected. It is mainly characterized by a clear frequency on the spectrum analysis of and (see Figure 5), which is usually due to different sheet cavity sizes on successive blades, but identical sheet cavities on opposite blades. Since this arrangement is stable in the rotating frame, no frequency is detected in the waterfall plot of or .

In this section we focus on results at speed 1.2 , in order to have the clearest available fluctuations. Indeed, it has been shown by Coutier-Delgosha et al. (2009) that the rotation speed has some effects on the intensity of the radial loads, but does not modify the instability types.

The decomposition into discrete wavelet (Daubechies n°5) is applied to the signal (Figure 7). Frequency is 682 Hz in the present case, and = 1.33 Hz. The analysis is performed during the pressure decrease between  s and  s, that is from the inception of the first instability regime (alternate cavitation) to the performance breakdown. For high frequencies, high amplitude is nearly systematically obtained, whereas below 100 Hz, several points can be noticed.

(i)Close to 40 Hz, two sequences are characterized by a significant increase of the intensity: between  s and 230 s (), and between  s and 320 s ().(ii)The inception of the super-synchronous instability can be observed clearly in the decomposition at frequency  Hz, and also in the adjacent decompositions, with a lower intensity. It is observed that this characteristic frequency progressively decreases, until a short synchronous regime (which can be detected here at the lowest frequency) is obtained, for 253 s < < 255 s.

Figure 7: Wavelet decomposition of signal .

The same process is applied to the unsteady pressure . The main resulting information is a very clear unsteadiness characterized by a frequency close to 20/40 Hz for 210 s < t < 230 s. Such frequency was already observed during the same time on the signal, but not so clearly. Figure 7 shows the waterfall plot for at rotation speed 1.2 . A high-intensity signal at frequency close to 30 Hz is obtained, between the regime of alternate cavitation and the inception of the super-synchronous instability, that is, for . This corresponds to 210 s < t < 230 s. It can be noticed that this frequency was much weaker at rotation speed (see Figure 5). Similarly, the second instability regime detected on the wavelet decomposition (for 300 s < < 320 s) clearly appears on the waterfall plot at 1.2 , while it can hardly be detected at speed (cf. Figure 5).

It can be concluded that the wavelet decomposition does not enable any better understanding of the mechanisms related to the different instabilities. Conversely, it can be used to detect precisely their inception, as can be seen for example in Figure 9 in the case of the super-synchronous regime.

5. Analysis of the Temporal Signals

Several types of instabilities are encountered during the process of pressure decrease at 1.2 and . In the present paper, attention is focused on the regime of alternate cavitation.

This regime is obtained for 130 s  s, that is, for . It is characterized by sheet cavities of different sizes on two consecutive blades but identical on opposite ones. This nonsymmetrical arrangement results in a clear frequency on the spectrum analysis of and (see Figure 8). Conversely, no evidence of this unbalanced flow pattern can be found in the waterfall plot of or , since it is stable in the inducer rotating frame.

(a)

(b)

Figure 8: Waterfall plot of (a) and (b) (, 1.2 ).

Figure 9: Wavelet decomposition for . Zoom on the time interval 225 s  s for  Hz (, ).

Supplementary information can be obtained in the present case from the analysis of the temporal signals: if a low-pass filter with cut-off frequency 100 Hz is applied to the signals of and , in order to get rid of the high-frequency fluctuations, it becomes clear that mean values of and are significantly different during all the regime (Figure 10). Note that the values of the forces are not given in ordinate for confidentiality reasons: only the zero is indicated to facilitate the discussion.

Figure 10: Temporal signals of and with low pass filter (cut-off frequency 100 Hz), , 1.2 .

Although the values of the forces are not given in ordinate for confidentiality reasons, it is clear that fluctuates around a negative average value , while is characterized by a positive mean value close to . It implies that the present nonsymmetrical flow pattern causes in addition a constant mean radial force, which does not rotate in the inducer rotating frame. It means that cavitation on opposite blades may not be fully identical. However, the frequency on the pressure waterfall plots shows that two identical cavitation patterns separated by 180° are also present inside the inducer. The term “cavitation pattern” is voluntary used, since it may concern too sheet cavitation on the blade suction side as tip cavitation close to the inducer shroud.

The most plausible cavitation pattern is thus composed of two identical cavitation areas on two opposite blades (1 and 3), while the two other ones (2 and 4) are both of different size. In this configuration, a constant radial load is obtained in the inducer rotating frame because of the dissymmetry between blades 2 and 4. This load, represented schematically in Figure 11, may result (excepted in a particular orientation of axes y and z) in different mean values of the components and .

Figure 11: Plausible configuration for . Rectangles represent schematically the respective lengths of the sheet cavities on the four blades (1.2 , ).

Another feature of the alternate cavitation regime can be remarked on the extended waterfall plots of or (see Figure 12). Frequency exactly equal to is obtained for .

Figure 12: Extended waterfall plot of (, ).

The following points can be checked from the analysis of the waterfall plots obtained at various rotation speeds and mass flow rates.

(i)It is a frequency, which changes with the inducer rotation speed.(ii)It can be detected mainly above the nominal flow rate, that is, within the range of mass flow rate that leads to most of the regimes of radial loads. It is particularly strong at high mass flow rate (see Figure 13).(iii)It is not specific to the alternate cavitation flow pattern: all regimes of developed and stable (i.e. nonrotating) cavitation result in the occurrence of the frequency. It is observed for example in the case presented in Figure 13 (, 1.18 ) that this frequency is obtained during three successive periods: (1) alternate cavitation, (2) short synchronous regime that occurs just after the super-synchronous one, and (3) regime of highly developed cavitation during the final performance breakdown. (iv)The frequency can be detected even in noncavitating conditions, but strongly attenuated.

Figure 13: Extended waterfall plot of (1.18 , ).

From these observations, no obvious explanation can be found regarding the occurrence of the frequency. Conversely, some points enable to improve the understanding of this phenomenon. For example, no peak at frequency is observed on the waterfall plots of the radial load intensity Fr, whatever the flow configuration is. It means that the radial load does not fluctuate; only its components do. Besides, the frequency is clearly detected on the waterfall plots of the bending moments and , which suggests that this instability is related to a shaft bend stress. On the other hand, the frequency peak is very small on the waterfall plots of the axial load , which eliminates any direct connection with an axial effort on the shaft.

Three different hypotheses can be postulated.

(i)A possible explanation is a rotating perturbation that would rotate at frequency in the reference steady frame, which would generate a frequency in the inducer rotating frame. Such phenomenon has been observed in experiments conducted in the Osaka University with a three-blade inducer, resulting in a frequency in the inducer frame [10]. It is called high-order rotating cavitation by Tsujimoto and colleagues. Although the present detected frequency is not the same, this example shows that the hypothesis of a perturbation with high rotation speed is plausible. However, it is not the most probable one, because of the exact value that is obtained.(ii)It can be the result of a rotor/stator interaction at the inducer outlet. A method proposed by Brennen [18] enables to obtain the characteristic frequencies related to such interaction. In the present case, the stator in the discharge pipe is equipped with 13 blades, the inducer has 4 blades, so frequencies , , , and should appear on the plots with similar intensities. So, frequency should not be predominant. Still it must be noticed that a frequency range 250 Hz–300 Hz has been previously identified as resonant frequencies of the installation. The peak may thus be of higher intensity than the other ones if it belongs to this range. However this hypothesis is also rejected, because the peak is of high intensity for all rotation speeds, from 0.6 up to 1.2 , while the value is not systematically included into the range of resonant frequencies.(iii)It can be related to a fluctuation of the radial load angle during the inducer rotation. More precisely, it can be anticipated that several fluctuations occurring during each inducer revolution may result in the measurement of the frequency . To investigate this point, a more in-depth analysis of the temporal signals of and has been conducted during the alternate cavitation regime for flow conditions and 1.18 , where frequency is remarkably clear on the waterfall plots. Figure 14 presents the general aspect of the two signals ( in grey in black), which confirms the difference of mean level that was mentioned previously.

Figure 14: Temporal signals of (in blue) and (in green) (, 1.18 ).

A zoom on a shorter time period is performed in Figure 15(a), to obtain only 25 rotations of the inducer. Oscillations at frequency can be clearly observed on both signals. Two secondary oscillations of lower amplitude are also detected on the signal, which may explain the presence of frequency in the spectrum. In the case of , these secondary oscillations do not appear so clearly. A second zoom is performed to focus on six rotations of the inducer only (see Figure 15(b)). Regular oscillations at frequencies and are observed: the two secondary peaks are regularly positioned between the main oscillations. On this figure, secondary peaks can be also seen in the signal of component (Figure 15(c)). It can be also noticed that the phase difference between and varies.

(i)The phase difference at low frequency , which corresponds to the gap between the vertical full lines in Figures 15(b) and 15(c), is small. It equals 90°, since it is related to the inducer weight.(ii)The phase difference between the secondary peaks is not the same. It can be estimated that it is close to 180°. It means that the phenomenon that generates the secondary oscillations results in fluctuations of and characterized by phase opposition. It must be mentioned that the analysis of various periods has enabled to draw this conclusion, while interpretation of Figures 15(b) and 15(c) is not so clear. On the basis of these observations, no definitive interpretation can be found for frequency .

(a)

(b)

(c)

Figure 15: Zoom on the signals of (in blue) and (in green), 1.18 , , and low pass filter with cut-off frequency 400 Hz.

However, a possible explanation, which could be validated only from high-speed visualizations, is given hereafter. It is based on a probable oscillation of the two sheet cavities of different sizes no. 2 and no. 4 denoted in Figure 11. These two areas may fluctuate around a mean size, and even reverse (i) because of the vertical pressure gradient due to the gravity within the inducer. The Froude number based on the inducer diameter and the axial velocity is about 0.2, which confirms that the gravity may have some influence on the flow during its passage in the machinery; (ii) or by effect of the two-phase structures which may be different at the top and at the bottom of the inducer, which would be also an indirect effect of the gravity.

In this case, the oscillation of components and at frequency , due to the inducer weight, may be coupled to this inversion of sheet cavities 2 and 4, which would occur twice during each rotation of the inducer, leading to a more complex signal. A scheme of the evolutions of and during a single rotation of the inducer is proposed in Figure 17. It is supposed here that the radial load related to the nonsymmetry between blades 2 and 4 is oriented according to the scheme given in Figure 16.

Figure 16: Expected orientation of and , and the blades (1.18 , ).

Figure 17: Scheme of the qualitative evolutions of (blue curve) and (red curve) during a single rotation of the inducer.

This particular orientation of and enables to explain the 180° phase difference between the two components in the secondary oscillations. Indeed, in such configuration, the inversion of the sheet cavities 2 and 4, which makes an inversion of the associated radial load, results in opposite variations of and . It must be noticed that the blade orientations drawn in Figure 17 represent schematically the direction of the radial load generated by nonidentical sheet cavities on opposite blades. It can be anticipated that these directions are close to the lines that would join the leading edge of the opposite blades.

It is observed in Figure 15 that the oscillations of are of larger amplitude than those of , which suggests that the angle between each component and the blades is not 45°. As indicated in Figure 16, the direction of may be close to the direction of the radial load related to the nonsymmetrical flow pattern on blades 2 and 4. The expected evolutions of and during a single inducer rotation are drawn in Figure 17 on the basis of this hypothesis.

The eight schemes at the top of Figure 16 show the eight successive positions of sheet cavities on blades 2 and 4, during a revolution of the machinery, by step of 45°. It can be observed that two inversions of the small and big sheet cavities occur during one period, each time they cross the horizontal direction. Between two inversions, the difference of size between sheet cavities 2 and 4 first increases reaching a maximal value, and then decreases. It is maximal when the radial load is vertical, since in such configuration the pressure gradient, which controls the pulsation of the cavities according to our assumption, has the strongest effect on the sheet cavities located on blades 2 and 4.

On the central part of Figure 17, the curves in full line correspond to the evolutions of and without the radial load generated by the unbalanced flow pattern on blades 2 and 4: a sinusoidal evolution; due to the respective orientations of the gravity, , and , with a phase difference between the load components, is obtained in that case. This behaviour is modified by the fluctuations of the sheet cavities on blades 2 and 4, leading to the curves drawn in doted lines.

The evolution of is composed of the following steps.

(i)It can be anticipated that the radial load will systematically be oriented from the small cavity towards the big one (i.e., from high-pressure areas towards low pressures). So, configuration 1 in Figure 17 is characterized by a component related to the flow dissymmetry that is positive, and even almost maximal, since blades 2 and 4 are close to the vertical direction. Between steps 1 and 3, the radial load reaches its maximal value, and then it decreases, so the curves in full and dotted lines become closer. It must be reminded that all charts here focus on a qualitative behaviour, while the real evolutions of and may depend strongly on the relative values of the radial load related to the inducer weight and the one related to nonsymmetrical cavitation on blades 2 and 4. (ii)Just after configuration 3, blades 2 and 4 cross the horizontal direction, so an inversion of the sheet cavity lengths is obtained. Consequently, the radial load due to the dissymmetry, and thus the sign of the component related to this load, is also inverted. It results in a sudden decrease of , under the curve in full line. Then, reincreases progressively until step 5 is reached, since its direction becomes closer to the one of the gravity. In the same time, the radial load due to the dissymmetry also increases, but it is negative, so it slightly counterbalances the effect of the inducer weight. (iii)Steps 6 to 8 lead to the same variations of , but in the opposite direction: between steps 5 and 6, components related to the inducer weight and to the dissymmetry are both decreased. When the leading edges of blades 2 and 4 cross the horizontal direction (just after configuration 7), sheet cavities invert again, so the radial load that they generate also does, and increases over the full line curve. Then, the effect of the weight decreases, leading to the decrease of .

The curve for is obtained by a similar process: from steps 1 to 3, the radial load due to sheet cavities on blades 2 and 4 is negative and decreases. After configuration 3, it inverts, so it becomes positive and increases more and more when the sheet cavities get closer to the vertical direction. After configuration 5, the radial load decreases again, and also does and comes closer to the curve in full line. Just after step 7, sheet cavities invert again, and the component related to the dissymmetry becomes negative and decreases down to its value at step 1.

This scenario generates for both components and three successive fluctuations during each rotation of the inducer, which is consistent with frequency that is observed on the waterfall plots. In the configuration presented in Figure 17, the three pulsations of occur at equal distance to each other, which explains the intense peak at frequency . Oscillations of component are not so regular during one rotation, which is consistent with observations from Figure 15, and which also explains that the spectral line at is weaker for than for .

6. Conclusion

The instabilities detected in a four-blade inducer have been discussed in the present paper by wavelet decomposition and analysis of the temporal signals of the two components and of the radial load on the shaft. It has been found that wavelet decomposition does not improve the understanding of the flow instabilities, while it enables a precise detection of inception and stop of each unbalanced flow pattern. An unexpected spectral line at , which occurs during alternate cavitation, has been discussed by direct analysis of the temporal signals of and . A possible scenario, based on the assumption that the gravity may have—directly or indirectly—some effect on the sheet cavities on the blades, is proposed. Only high-speed video would enable to confirm this hypothesis.

Nomenclature

| Frequency (Hz)
: | Frequency of the inducer rotation (Hz)
: | Axial component of the force on the shaft (N)
: | Radial components of the force on the shaft (N)
: | Rotation speed ()
: | Reference rotation speed ()
: | Pressure elevation (Pa)
: | Absolute pressure in the inlet pipe (Pa)
: | Absolute pressure in the delivery pipe (Pa)
: | Total absolute pressure in the inlet pipe (Pa)
: | Total absolute pressure in the delivery pipe (Pa)
: | Vapor pressure (Pa)
: | Fluctuating pressures at the inducer inlet (Pa)
: | Fluctuating pressure at the inducer outlet (Pa)
: | Torque on the inducer shaft (Nm)
: | Bending moments on the shaft (Nm)
: | Volume flow rate (/s)
: | Nominal flow rate (/s)
: | Main flow rate (in the delivery pipe) (/s)
: | Secondary flow rate (for axial equilibrium) (/s)
: | Inducer tip radius (m)
: | Time (s)
| Outlet ratio hub radius/tip radius (—)
| Head coefficient ()/() (—)
| Torque coefficient /() (—)
| Flow rate coefficient Q/( ()) (—)
| Cavitation number ()/() (—).

Acknowledgments

The present study was performed in the frame of contractual activity with SNECMA Moteurs and the CNES (French Space Agency). The technical staff of the LML laboratory was much involved in the development and the operation of the test facility. The authors wish to thank especially J. Choquet and P. Olivier for their collaboration.

References

K. Kamijo, T. Shimura, and M. Watanabe, “An experimental investigation of cavitating inducer Instability,” ASME Paper 77-WA/FW-14, 1977.

B. Goirand, A. Mertz, F. Joussellin, and C. Rebattet, “Experimental investigation of radial loads induced by partial cavitation with liquid hydrogen inducer,” in Proceedings of the 3rd International Conference on Cavitation (ImechE '92), vol. C453/056, pp. 263–269, Cambridge, UK, 1992.

J. de Bernardi, F. Rossellini, and A. Von Kaenel, “Experimental analysis of instabilities related to cavitation in turbopump inducer,” in Proceedings of the 1st International Symposium on Pump Noise and Vibrations, pp. 91–99, Paris, France, 1993.

R. A. Furness and S. P. Hutton, “Experimental and theoretical studies of two-dimensional fixed-type cavities,” vol. 97, pp. 515–522, 1975.

Q. Le, J. P. Franc, and J. M. Michel, “Partial cavities: global behavior and mean pressure distribution,” ASME Transactions Journal of Fluids Engineering, vol. 115, no. 2, pp. 243–248, 1993.

J. B. Leroux, O. Coutier-Delgosha, and J. A. Astolfi, “A joint experimental and numerical study of mechanisms associated to instability of partial cavitation on two-dimensional hydrofoil,” Physics of Fluids, vol. 17, no. 5, article 052101, pp. 1–20, 2005.

O. Coutier-Delgosha, B. Stutz, A. Vabre, and S. Legoupil, “Analysis of cavitating flow structure by experimental and numerical investigations,” Journal of Fluid Mechanics, vol. 578, pp. 171–222, 2007.

Y. Tsujimoto, “Cavitation instabilities in inducers,” Tech. Rep. AVT-143 RTO AVT/VKI Lecture Series, von Karman Institute, Rhode-Saint-Genèse, Belgium, 2006.

F. Joussellin, J. De Bernardi, B. Goirand, and Y. Delannoy, “Analyse par films ultrarapides de poches de cavitation sur l'Inducteur de la Turbopompe à Hydrogène d'un Moteur Fusée,” in 5ème colloque de visualisation et de traitement d'images en Mécanique des Fluides, Poitiers, France, 1992.

A. Fujii, S. Azuma, Y. Yoshida, Y. Tsujimoto, H. Horiguchi, and S. Watanabe, “Higher order rotating cavitation in an inducer,” International Journal of Rotating Machinery, vol. 10, no. 4, pp. 241–251, 2004.

O. Coutier-Delgosha, Y. Courtot, F. Joussellin, and J. L. Reboud, “Numerical simulation of the unsteady cavitation behavior of an inducer blade cascade,” AIAA Journal, vol. 42, no. 3, pp. 560–569, 2004.

Y. Yoshida, Y. Sasao, K. Okita, S. Hasegawa, M. Shimagaki, and T. Ikohagi, “Influence of thermodynamic effect on synchronous rotating cavitation,” ASME Transactions Journal of Fluids Engineering, vol. 129, no. 7, pp. 871–876, 2007.

F. Bakir, S. Kouidri, R. Noguera, and R. Rey, “Experimental analysis of an axial inducer influence of the shape of the blade leading edge on the performances in cavitating regime,” ASME Transactions Journal of Fluids Engineering, vol. 125, no. 2, pp. 293–301, 2003.

Y. Yoshida, Y. Tsujimoto, D. Kataoka, H. Horiguchi, and F. Wahl, “Effects of alternate leading edge cutback on unsteady cavitation in 4-bladed inducers,” ASME Transactions Journal of Fluids Engineering, vol. 123, no. 4, pp. 762–770, 2001.

Y. D. Choi, J. Kurokawa, and H. Imamura, “Suppression of cavitation in inducers by J-Grooves,” ASME Transactions Journal of Fluids Engineering, vol. 129, no. 1, pp. 15–22, 2007.

H. Horiguchi, T. Takashina, and Y. Tsujimoto, “Theoretical analysis of cavitation in inducers with unequally spaced blades,” JSME International Journal, Series B, vol. 49, no. 2, pp. 473–481, 2006.

O. Coutier-Delgosha, G. Caignaert, G. Bois, and J. B. Leroux, “Influence of the blade number on inducer cavitating behavior,” Journal of Fluids Engineering. In press

C. E. Brennen, “Pump vibration, concepts,” in Hydrodynamics of Pumps, chapter 8, NREC and Oxford University Press, 1994.