The Action Mechanism of Rotor–Stator Interaction on Hydraulic and Hydroacoustic Characteristics of a Jet Centrifugal Pump Impeller and Performance Improvement

: The aim of this study was to investigate the action mechanism of the rotor–stator interaction (RSI) on the transient flow field and hydrodynamic noise field inside the impeller of jet centrifugal pumps (JCPs) and optimize effects of the guide vane on the hydraulic and hydroacoustic characteristics of the impeller. The numerical method of CFD (computational fluid dynamics) coupled with CFA (computational fluid acoustics) was used to analyze the correlation between the guide vane and the flow/sound performances of the impeller. The orthogonal test method, with the hydroacoustic performance of the impeller taken as the objective, was used to optimize the structural parameters of the guide vane for the stability of the hydraulic performance of the JCP. The results show that the RSI leads to a significant increase in the hydroacoustic level of the impeller, but it is indispensable for improving the hydraulic performance of the pump. The RSI effect on the fluctuation intensity of the transient flow field inside the impeller is much more sensitive than the time ‐ average, and the fluctuation intensity of the flow field is positively correlated with the vortex intensity inside the impeller. When the impeller geometry is constant, the evolution processes of the flow field inside the impeller are mainly related to the blade number of the guide vane; when the number of guide vanes is given, the RSI effect on the hydroacoustic characteristic of the impeller is characterized by a positive correlation between the total sound pressure level (SPL) and the fluctuation intensity of the flow field. The frequency spectrum characteristics of the hydroacoustic SPL of the impeller are not consistent with the pressure fluctuation characteristics inside the impeller. The pressure fluctuation characteristics are related not only to the blade number and speed of the impeller but also to its wake characteristics determined by the guide vane. The optimization scheme for the stable hydraulic performance of the JCP significantly reduced the total SPL of the impeller compared with the original scheme, which verifies the feasibility of using the weight matrix optimization method to obtain the global optimization scheme.


Introduction
The rotor-stator interaction (RSI) is a ubiquitous phenomenon in turbomachinery [1]. It spreads up-and downstream, resulting in the formation of strong flow field fluctuations in the vaneless area between the dynamic and static cascades, which have an important impact on the stability, safety, and robustness of machines. A certain degree of the RSI effect cannot be avoided, but excessive RSI effects can reduce the efficiency of the unit, cause fatigue damage of flow passage components, and increase vibration and noise, among other negative outcomes [2,3]. During the past several decades, researchers have used numerical or experimental methods to study the RSI effect on the transient flow field, hydraulic excitation force, vibration, noise, and other characteristics of vane pumps. Dong et al. [4] experimentally studied the effect of modifications to tongue and impeller geometries on the flow structure and the resulting noise in a centrifugal pump, and their results demonstrated that the primary sources of noise were associated with interactions between the nonuniform outflux from the impeller (jet/wake phenomenon) and the tongue. González et al. [5] numerically simulated the dynamic and unsteady flow effects due to the impeller-volute interaction inside a centrifugal pump and showed that the strong effects of the secondary flow were concentrated in radial positions close to the impeller exit. Additionally, pressure fluctuations at the blade passing frequency revealed that the blade tongue interacted with the flow at the impeller outlet plane, and the interaction had clear effects, such as increased fluctuation levels under off-design conditions, limited operation ranges, and increased losses. Guo et al. [6] measured pressure fluctuations and the radial fluid forces induced by the rotor-stator interaction and investigated their relationship, and they demonstrated that both the blade pressure fluctuations and the volute static pressures were circumferentially nonuniform (not axisymmetric) under off-design operating conditions. Rodriguez et al. [7] carried out a theoretical analysis to qualitatively predict and explain the moving blade passing frequency and the relationship between their amplitudes. They presented a way to characterize the vibration in the frequency domain resulting from the RSI. Using a one-dimensional approach, Nicolet et al. [8] modeled, simulated, and analyzed the rotor-stator interaction of a scale model of a Francis pumpturbine and a related test rig and found that the pressure patterns of the diametrical mode in the vaneless gap were dominated by the pressure patterns of the resonance. Zhu et al. [9,10] simulated the flow in the centrifugal pump with different volute casings by using a shear stress transport turbulent model, and their findings indicated that the blade passing frequency dominated the pressure fluctuation, radial hydraulic force, and torque in the impeller, whose fluctuations were periodic. Furthermore, they studied the rotor-stator interaction of three tongues of the centrifugal pump using dynamical subgrid-scale models. Rodriguez et al. [11] studied a 100 MW pump-turbine by making use of pressure measurements and vibrations measured in a guide bearing that rotated with the shaft. From their results, they recommended that the vibration caused by rotation with the shaft be measured to detect RSI characteristics in large hydraulic turbomachinery. Zhou et al. [12] simulated a volute-type centrifugal pump by using LES models to reveal the impeller-volute interaction under rotating stall conditions, and they found that the rotating stall phenomenon occurred when the flow rate was decreased to 0.70 of the nominal flow rate; as the flow rate was further decreased, the area of the stall cells increased. Under the rotating stall condition, the blade passing frequency (i.e., three times the rotation frequency) decreased by 50% because of alternate stalled and unstalled passages; when the rotating stall occurred, the amplitude of the pressure fluctuation was much higher than that at unstalled points. Zhang et al. [13] studied the pressure pulsation characteristics induced by the rotor-stator interaction in a nuclear reactor coolant pump. At the nominal flow, the flow rate distribution in the guide vanes was relatively uniform, and the amplitude of the pressure pulsation reached the minimum. Furthermore, the vortex shedding or backflow from the impeller blade exit had the same frequency as the pressure pulsation, and that the absolute value of phase differences reflected the vorticity intensity. Chalghoum et al. [14,15] analyzed the rotor-stator interaction (impeller and volute) and its effects on the complex internal flow field within a centrifugal pump for a wide range of operating flow rates, and they found that the pressure distributions around the volute tongue were inversely proportional to the evolution of the relative angle between the blade exit and volute tongue. As a result of the rotor-stator interaction, the amplitude and frequencies of the pressure fluctuations rose from the inlet to the outlet of the impeller for the whole flow range. At the volute tongue region, the pressure amplitude was much larger than that at other regions because of the presence of vortices.
Scholars have extensively researched the influence of the rotor-stator interaction on the hydraulic performance of hydraulic turbomachinery, but the effect on hydroacoustic performance has rarely been studied. In this study, a typical jet centrifugal pump (JCP) was taken as the study object to investigate the action mechanism of the rotor-stator interaction on its impeller hydraulic and hydroacoustic characteristics and optimize these parameters through multi-objective optimization.

The Model Pump
The model JCP had the following parameters: rated flow Qd = 2.5 m 3 /h, rated water head Hd = 23 m, rated efficiency ηd = 20%, rotational speed n = 2850 r/min, shaft passing frequency SPF = 47.5 Hz, blade passing frequency of the impeller BPFI = 285 Hz, and blade passing frequency of the guide vane BPFG = 237.5 Hz. The main geometric parameters of the impeller and guide vane are listed in Table 1. Two-and three-dimensional structure diagrams of the model pump are shown in Figure 1.

Research Routes
The research route followed in this study is shown in Figure 2. The impeller structure remained unchanged while analyzing the correlation between the main geometric parameters of the guide vane and the transient flow field inside the impeller and its induced noise. The hydroacoustic performance of the impeller was taken as the objective, and the guide vane was optimized to increase the stability of the hydraulic performance of the JCP.

Meshes
The whole computational domain includes the impeller, guide vane, jet, front pump cavity, back pump cavity, inlet pipe, outlet pipe, and pump body. ANSYS-ICEM (16.0 ANSYS, Inc., Canonsburg, PA, USA) was employed for the discrete computational domain, hexahedral structured meshes were employed for the computational domain of the impeller and guide vane, and tetrahedral unstructured meshes were employed for other flow passage components. The grid independence verification is shown in Table 2, divided into six mesh schemes according to the density of the grid. The grid convergence index (GCI) recommended by ASME is used to estimate the grid dispersion error [16]. GCI is an index with a 95% confidence interval, which represents the distance between the more dense grid and the progressive value in two comparative grids, and predicts the influence of further grid refinement on the solution. For Schemes 5 and 4 grids, the calculated head uncertainties are 0.19% and 6.43%, respectively, and efficiency uncertainties are 0.37% and 7.59%. In order to balance the relationship between computing accuracy and resources, Scheme 4 grid is selected for numerical calculation in this paper. The automatic wall treatment method is used to deal with the boundary layer flow, the method can generate a processing mode that is not sensitive to Y + and allow users to refine or coarsen the near-wall mesh arbitrarily [17]. The impeller and guide vane mesh of the model pump are shown in Figure 3.

Numerical Settings and Boundary Condition
CFX (16.0, ANSYS Inc., Pittsburgh, PA, USA) was employed for the flow field calculation. The impeller flow fields were calculated in the rotating coordinate system in the multi-coordinate system, and the flow fields of the other components were calculated in the stationary coordinate system. The k-ω turbulence model was adopted, which is better than the k-ε model in the aspect of near-wall flow analysis [17], and GGI (general grid interface) technology was employed to exchange data between the rotor and stator. The automatic near-wall treatment method was applied to the flow near the wall. Pressure inlet and velocity outlet boundary conditions were set, and the inlet pressure value was set according to the value obtained from the experiment. All solid walls were characterized as non-slip walls, and roughness was set to 25 μm according to actual processing. The governing equations were solved by a second-order upwind formula, and time terms were discretized by a second-order implicit scheme. All residuals were less than 10 −5 . The time-step was set to 0.000117 s, which means that the impeller rotated about 2 degrees per time-step. First, steady calculations were carried out, and the results were taken as the initial value for unsteady calculations. Pressure fluctuation information files from 10 rotation cycles on the impeller surface were output for hydroacoustic calculations after the flow field showed a stable periodic change.

Time-Average and Fluctuation Intensity of the Flow Field
The unsteady flow characteristics of the model pump were analyzed by a statistical method. The unsteady physical quantities include two parts on each grid node (x, y, z) of the flow domain, the time-averaged component  , and the periodic component   , where the periodic component represents the variation in physical quantities in a rotation period of the impeller [18]. The two components are respectively expressed as follows: where n, N, and t0 are the grid nodes, the number of samples in an impeller rotation period, and the starting time of the rotation cycle, respectively.
The non-dimensional standard deviation of the periodic component   is defined as the flow field fluctuation intensity, which includes the pressure fluctuation intensity, velocity fluctuation intensity, turbulence energy fluctuation intensity, and so on. The pressure fluctuation intensity is expressed as The velocity fluctuation intensity is expressed as where U2 is the circumferential speed at the impeller outlet.

Method of Computational Fluid Acoustics (CFA)
In this paper, LMS Virtual.Lab 12.0 ( Siemens, Berlin, Germany ) and the acoustic finite element method (FEM) were used to solve the hydrodynamic noise of the impeller, and the wide-band calculation method of fan noise was adopted, which did not need a Fourier transform when importing pressure fluctuation data from the impeller surface. The element length of the acoustic FEM mesh was determined by Equation (5): where c and fmax are the speed of sound and the maximum computational frequency.
Because of the structural characteristics of the model JCP, the length of the mesh element was set to 6 mm. When the hydraulic noise was solved, the surface of other flow passage components was set as the total reflection wall. That is, there was no transmission of sound through the wall, and the sound only traveled upstream and downstream along the water. The boundary conditions of the inlet and outlet of the pump were defined as the property of total sound absorption. The monitoring point was set at the pump outlet pipe at 3 times the pipe-diameter.
For measuring noise, the 1/3 octave analysis method can reflect the spectrum characteristics of the noise source with relatively high accuracy so as to understand the generation mechanism and put forward reduction measures [19]. In this study, the A-weighted method, which is closest to reflecting the subjective feeling of the human ear, was employed to evaluate the hydroacoustic sound pressure level at each frequency. The total sound pressure level (SPL) was obtained by superposing the Aweighted SPL of each center frequency of 1/3 octave bands, which is the energy summation method shown in Equation (6): where LP is the total SPL after evaluation and summation, dBA; M is the summation number; and LPi is the SPL of each center frequency of 1/3 octave bands after evaluation.

Numerical Validation
Numerical calculations of the flow field and sound field were validated using the JCP test-bed of hydraulic performance and hydrodynamic noise in the Key Laboratory of Fluid Machinery and Systems, Gansu Province. The test system is shown in Figure 4. The position of the hydrophone at the pump outlet is the same as that of the noise monitoring point in the numerical calculation.   Figure 5a compares the curves of hydraulic performance between the experiment and the numerical simulation. The variation trend of the performance curve is consistent between numerical calculations and the experiment, and the values at different working points are similar. The maximum error of the water head is 4.6%, and efficiency is 4.1%: these results confirm the accuracy of the hydraulic calculation results.
In Figure 5b, the frequency response curves of the sound pressure level at the monitoring point are compared between numerical calculations and the experiment under rated conditions. As a whole, the level of noise from the experiment is higher than that from the calculation; both the hydroacoustic SPL of the guide vane and that of the impeller have a pronounced amplitude peak at BPFI. This is because the JCP contains mechanical noise and hydrodynamic noise that are from different flow passage components and have different characteristics. These different noises can be enhanced by overlapping each other, but they also can be weakened by subtracting from each other. In addition, the flow field calculation cannot completely and accurately capture the flow separation, secondary flow, vortex, or local slight cavitation.

Method to Determine the Global Optimization Scheme
A weight matrix analysis method was employed to obtain the global optimization scheme of a multi-objective orthogonal experimental design. The comprehensive balance method or comprehensive scoring method involves extensive calculations and neglects the differences and connections between indexes, and the weight matrix analysis method overcomes these shortcomings. The influence degree of factors on the index and the primary and secondary order can be assessed by calculating the weights of different levels of factors on the experiment results to quickly reveal the global optimization scheme [20]. The method requires the construction of a three-layer data model with n factors and m levels for each factor according to the data structure, as shown in Table 3. Table 3. Data structure of the weight matrix analysis method.

Index Layer Indices of Experimental Investigation
In Table 2, the first layer is the index layer of the experimental investigation, and the matrix M (as shown in Equation (7)) is established. The average value of factor i at level j is defined as kij. If a larger index value corresponds to a better value, then kij = Kij; otherwise, kij = 1/Kij. The second layer is the factor layer, and the matrix T (as shown in Equation (8)) is established. The third layer is the level layer, and the matrix R (as shown in Equation (9)) is established, where ri is the range of factor i. 11 12 1 The weight matrix ω affecting the experimental index is defined as Because K11T1 is the ratio of level 1 to all levels of factor 1, the ω1 is product of K11T1 and R1, which reflects the influence degree of level 1 of factor 1 on the index.

Correlations between the Guide Vane and Hydraulic and Hydroacoustic Characteristics of the Impeller
The values of the main geometric parameters for the original positive guide vane are shown in Table 4. When analyzing the influence of a certain parameter on the hydraulic and hydroacoustic characteristics of the impeller, other parameters should remain unchanged.  Table 5 shows the correlations between the guide vane and hydraulic and hydroacoustic characteristics of the impeller. As shown in the table, the guide vanes play a very important role in the hydraulic performance of the pump and the hydroacoustic characteristics of the impeller. Without positive and return guide vanes, the water head, hydraulic efficiency, the shaft power of the pump, and the SPL of the impeller are greatly reduced; this also holds if only the positive guide vane removed.         Figure 9 shows the evolution processes (the processes of the impeller rotating 60°; the same applies below) of the flow field inside the impeller and guide vane for three models. For the model without positive and return guide vanes, the evolution process of vortices inside the impeller is much different from that for the other two models. When the impeller is rotating, the vortex shapes change little in a single flow channel and distribute symmetrically in six channels. For the model without positive guide vanes, the variation law of the transient flow field is similar to that for the prototype, but the strength and size are smaller. For the prototype, the presence of the positive guide vane causes substantial channel blockage in the interface area of the dynamic and static cascades. By combining these results with Figure 8, we find that the frequency spectrum characteristics of the hydroacoustic SPL of the impeller are related not only to the blade number and speed of the impeller but also to its wake characteristics determined by the guide vane.  The analysis results in Figures 6-9 show that the RSI effect leads to a significant increase in the hydroacoustic SPL of the impeller, but it is indispensable for improving the flow field inside the impeller and the hydraulic performance of the pump. Table 6 shows the correlations between the number of guide vanes and the hydraulic and hydroacoustic characteristics of the impeller. The data show that as the blade number of the guide vane increases, the energy head of the impeller outlet and the shaft power increase monotonously, while the total SPL of the impeller and the hydraulic performance of the pump first increase and then decrease.  Figure 10 shows the time-average and the fluctuation intensity curves of pressure and velocity on the blade surface of the impeller with different numbers of guide vanes. As shown in the figures, as Z2 changes, the time-average curves of both the pressure and velocity rise linearly by a similar law along the increasing radius. The time-average curve of pressure is lower when Z2 is 3, which is consistent with the results in Table 5. The fluctuation intensity of both pressure and velocity are affected very significantly by Z2 and drastically decline as it increases.      Figure 13 shows the evolution processes of the flow field inside the impeller and guide vane with different numbers of guide vanes. The evolution rules of the flow field are observed to have obvious differences when Z2 is different. The scale of vortices inside the impeller and the degree of flow blockage at the inlet of the guide vane are the most pronounced when Z2 is 3, which results in extremely serious hydraulic losses. Thus, the hydraulic performance of the pump is much lower in this scenario than in the other two schemes in Table 5. With the increase in Z2, the strength of the vortices inside the impeller weaken gradually. A comparison with Figure 10 shows a positive correlation between the fluctuation intensity of the transient flow field on the blade surface and the intensity of the vortices inside the impeller: that is, the greater the intensity of the vortices inside the impeller, the greater the fluctuation intensity of the pressure and velocity on the blade surface.  Table 7 shows the correlations between the wrap angle of the positive guide vane and the hydraulic and hydroacoustic characteristics of the impeller. The data show that as the wrap angle increases, the energy increment and the shaft power of the impeller show a monotonous increasing trend. The total SPL decreases gradually, and the hydraulic performance of the pump first increases and then decreases. That is because the bigger wrap angle and the longer blade can guide the flow more continuously than shorter ones. Means, the less turbulences loss will occur and the performances should be higher; however, if the blades are too long, the blade roughness will cause more friction loss, and the efficiency will decrease again.  Figure 14 shows the time-average and the fluctuation intensity curves of pressure and velocity on the blade surface of the impeller with different wrap angles of the guide vane. The time-average curves of the pressure and velocity are highly similar along the radius with different φ2, but the curves of fluctuation intensity are affected very significantly and decline drastically as the wrap angle increases.     Figure 17 shows the evolution processes of the flow field inside the impeller and guide vane with different wrap angles of the guide vane; the evolution rules of the flow field are basically the same in the dynamic and static cascades. With the increase in wrap angles, the size and intensity of vortices inside the impeller and the degree of flow blockage at the inlet of the guide vane decrease gradually. This is because when the wrap angle is larger, the change gradient of the blade angle is smaller, the ability of the blades to control the flow is stronger, and the flow inside the guide vane is smoother, resulting in a better fluid transition between the dynamic and static cascades.  Table 8 shows the correlations between the inlet angle of the positive guide vane and the hydraulic and hydroacoustic characteristics of the impeller. With the increase in the inlet angle, the energy increment of the impeller shows a monotonous increasing trend, while the hydraulic performance of the pump and the total SPL of the impeller decline gradually, and the shaft power changes little.  Figure 18 shows the time-average and the fluctuation intensity curves of pressure and velocity on the blade surface of the impeller with different inlet angles of the positive guide vane. Although the inlet angles of the guide vane do not have a very significant effect on the time-average and fluctuation intensity of pressure and velocity, the fluctuation intensity is still more sensitive than the time-average. With the increase in βG1, the curves of the fluctuation intensity of the velocity show an overall decline, and the fluctuation intensity of the pressure first decreases and then increases near the impeller outlet.      Table 9 shows the correlations between the inlet diameter of the positive guide vane and the hydraulic and hydroacoustic characteristics of the impeller. As the inlet diameter increases, the energy head of the inlet and outlet and the energy increment increase monotonously, whereas the total SPL declines gradually, the shaft power changes are small, and the hydraulic performance of the pump first decreases and then increases.  Figure 22 shows the time-average and the fluctuation intensity curves of pressure and velocity on the blade surface of the impeller with different inlet diameters of the positive guide vane. With the increase in DG1, the time-average of pressure and velocity on the blade surface of the impeller change little, the inlet diameter of the positive guide vanes have low impact to the flow of the impeller runner, therefore the divergence of the time-average lines of pressure is small. While their fluctuation intensity decreases very dramatically. This is because the increase in the inlet diameter leads to an increase in the clearance between the dynamic and static cascades, and the RSI effect between the guide vane and the impeller is weakened.    The analyses in Sections 3.1-3.5 show that a good guide vane can improve the transient flow field structure inside the impeller, reduce the flow field distortion phenomena, and improve the hydroacoustic characteristic of the impeller and the hydraulic performance of the JCP. When the impeller geometry is constant, the evolution process of the flow field inside the impeller is mainly related to the blade number of the guide vane, and other parameters only affect the strength and size of the vortices without changing their basic law. When the number of guide vanes is given, the RSI effect on the hydroacoustic characteristic of the impeller is characterized by a positive correlation between the total SPL and the fluctuation intensity of the flow field.

Optimization Design of the Guide Vane
The total SPL of the impeller was taken as the objective to optimize the effects of the rotor-stator interaction on the hydraulic/hydroacoustic characteristics of the JCP impeller, and an orthogonal experimental method was used to optimize the structural parameters of the guide vane for the stability of the hydraulic performance of the JCP.

The Scheme and Results of the Orthogonal Experiment
The L9(3 4 ) standard orthogonal table [21] was used to construct the experimental design of the structural parameters of the guide vane: with the three factor levels discussed in Section 3 of the paper, nine schemes of the guide vane structure were obtained, and then their flow field and sound field were calculated. Table 10 specifies the levels of the orthogonal experimental factors, and A, B, C, and D are factor codes that correspond to Z2, φ2, βG1, and DG1, and Table 11 is the experimental schemes and results.

The Intuitionistic Analysis
The range of values reflects the influence degree of the change in the factor level on the index. The greater the range, the greater the influence of the selected factor on the index. Therefore, the most significant factor corresponds to the largest range. Because of the positive correlation between the head and efficiency, the hydraulic performance of the pump only reflects efficiency. A range analysis was carried out on the numerical calculation results in Table 10. Range analysis was carried out by Equations (12) and (13), where ki is the average value of j levels of factor i; s is range.
The results are shown in Table 12. As shown in the table, the sequence of the influence degree of the guide vane parameters on the hydraulic efficiency of the pump is A, D, C, B, and the sequence of the influence degree of the same parameters on the hydroacoustic performance of the impeller is D, C, B, A.  Figure 26 presents the effective curves of the optimization indexes. For the hydroacoustic performance of the impeller, the smaller the index, the better the performance; therefore, the relative optimal parameter combination is A1B3C3D3. For hydraulic efficiency, the higher the index, the better the performance, so the relative optimal combination is A2B2C2D3.

The Global Optimization Scheme
The matrix analysis method described in Section 2.6 was employed to obtain the global optimization scheme for both the hydraulic efficiency of the pump and hydroacoustic performance of the impeller. The weight matrices corresponding to hydraulic efficiency and the total SPL were calculated by Equations (7)- (11). For hydraulic efficiency, Kij = kij; for the hydroacoustic SPL, Kij = 1/kij. The results are shown in Equations (14) and (15), and Equation (16) is the global matrix, which is the average of Equations (14) and (15 These results show that the primary and secondary order of influence of different factors on the global performance is D, A, C, B. Factors A and B are optimal at level 2, factor C is optimal at level 1, and factor D is optimal at level 3. Therefore, the global optimal scheme is A2B2C1D3: that is, the number of blades Z2 = 5, wrap angle φ2 = 68°, inlet angle βG1 = 5°, and inlet diameter DG1= 127 mm. Figure 27 compares the performance of the prototype and that of the optimization scheme. After optimization, the hydraulic performance of the pump is stable under all conditions, and the hydroacoustic performance of the impeller is significantly reduced, so the desired objectives are realized. (a) Hydraulic performance of the JCP (b) Hydroacoustic performance of the impeller Figure 27. Performance comparison between the prototype and optimization scheme.

Conclusions
In this paper, a numerical method that couples CFD/CFA was used to analyze the correlation between the guide vanes and the hydraulic/hydroacoustic characteristics of an impeller. In addition, the hydroacoustic performance of the impeller was taken as the objective to optimize the structural parameters of the guide vane for the stability of the hydraulic performance of the JCP using the orthogonal test method. Several conclusions are drawn from the results: (1) The RSI leads to a significant increase in the hydroacoustic level of the impeller, but it is indispensable for improving the hydraulic performance of the pump; a good guide vane can improve the transient flow field structure inside the impeller, reduce the flow field distortion phenomena, and improve the hydraulic and hydroacoustic performance of the impeller. (2) The fluctuation intensity of the transient flow field inside the impeller is much more sensitive than the time-average to the RSI effect; inside the impeller, it has a positive correlation between the fluctuation intensity of the flow field and the vortex intensity. (3) When the impeller geometry is constant, the evolution processes of the flow field inside the impeller are mainly related to the blade number of the guide vane, and other structural parameters only affect the strength and size of vortices without changing their basic law. When the number of guide vanes is given, the RSI effect on the hydroacoustic characteristic of the impeller is characterized by a positive correlation between the total SPL and the fluctuation intensity of the flow field. (4) The frequency spectrum characteristics of the hydroacoustic SPL of the impeller are related not only to the blade number and speed of the impeller but also to its wake characteristics determined by guide vane. Moreover, the frequency spectrum characteristics of the hydroacoustic SPL are not consistent with the pressure fluctuation characteristics inside the impeller. (5) Under the design conditions, the optimization scheme for the stable hydraulic performance of the JCP significantly reduces the total SPL of the impeller compared with the original scheme, which verifies the feasibility of using the weight matrix optimization method to obtain the global optimization scheme.