1. Introduction
Kaplan turbines are primarily utilized in hydropower stations with low head and large flow [
1]. They are distinguished by the capability to adjust the turbine blades in response to varying water flow conditions, thereby ensuring optimal operational efficiency. However, Kaplan turbines can cause vibration and pressure pulsation during operation, stemming from hydraulic, mechanical, or electromagnetic factors. Excessive vibration can result in mechanical fatigue, structural damage, and equipment failure. Similarly, pressure pulsation may lead to resonance within the turbine’s internal structure. Both phenomena are directly linked to the operational safety of the power station. This will not only impact its operational efficiency but also has the potential to result in equipment damage and even pose a risk of safety accidents. During turbine operation, periodic pressure fluctuations occur due to water flow, de-flow, secondary water impact, and the interaction between runner blade and water flow. These pressure pulsations exert forces on the turbine’s structural components, causing periodic changes that may induce mechanical vibration. Particularly when the frequency of these pressure pulsations aligns closely with or matches the natural frequency of the turbine structure, resonance can occur and exacerbate vibration issues.
The root causes of turbine vibration include hydrodynamic effects, asymmetry or damage to blades or components, instability in water flow, and structural resonance. The main measures for controlling vibration are adjusting blade angles, optimizing structural design, and using damping materials. In terms of the analysis of vibration causes and methods for controlling vibration, a Kaplan turbine unit exhibits significant vibration and generates high levels of noise during no-load operation at high speeds. Moraga et al. [
2] conducted modal analysis testing on the runner to address this issue, as well as high-speed no-load testing on the unit under accelerated conditions. The results indicated that the source of the high amplitude vibration and excessive noise was located near the head cover. By modifying the diameter and distribution of the air inlet, both vibration and noise levels were significantly reduced. Roig et al. [
3] conducted a no-load speed test of the guide vane and runner blade combination at various angles on a reduced-size Kaplan turbine model, under both non-cavitating and cavitating conditions. The study revealed that the presence of multiple simultaneous vortex structures in the vaneless region led to torsional vibration of the rotor. Furthermore, it was observed that the induced structural response on the low-pressure side of the turbine increased when operating under cavitation conditions. Dunca et al. [
4] conducted a field test to determine the optimal three-dimensional combination of runner blade opening and guide blade opening for different water head values. The optimal cam for vibration levels when operating in different states was determined.
In the realm of numerical simulation, advanced technologies such as computational fluid dynamics (CFD) are commonly utilized to forecast and analyze the flow field characteristics of hydraulic turbines [
5]. These technologies are also combined with vibration characteristics to comprehend the mechanism behind vibration generation. The study of fluid-structure coupling vibration characteristics of Kaplan turbine runner is relatively limited compared to other types. This is attributed to the fact that the blades rotate with changes in load, and the interaction within the blades may impact the additional mass and damping of the runner. Zhang et al. [
6] conducted a comparison between natural frequencies predicted by acoustic FSI and those from one-way FSI analysis in order to analyze how internal blade interaction influences additional mass and damping mechanisms. Puolakka et al. [
7] proposed a simplified method for calculating the unsteady vibration of turbine shafting based on unsteady airfoil theory. The method assumes that the runner, as a rigid body, oscillates in rotation and axial heave, and decomposes the reaction force into additional mass and damping. It was observed that over the frequency range of interest, there were moderate changes in increased mass and little changes in increased damping.
In terms of experimental research, vibration data are collected through model tests or prototype tests to verify the accuracy of numerical simulations and the effectiveness of vibration control strategies. For experimental monitoring methods, due to the limited number of studies on the feasibility of applying vibration measurement-based diagnostic techniques to turbine units and the absence of appropriate sensors installed on the units, Pennacchi et al. [
8] utilized this diagnostic tool to address the issue of hydraulic instability and proposed a highly automated condition monitoring system. Feng et al. [
9] conducted a study to measure the vibration signals of the rotor and draft pipe of a Kaplan turbine under various cavitation states. They utilized the multifractal detrended fluctuation analysis (MF-DFA) method to analyze the vibration signals and found that both the runner and draft pipe exhibited multifractal characteristics.
Pressure pulsation can lead to turbine vibration and the vibration of its attached structure. Prolonged and intense vibration poses a significant threat to the safe operation of the unit, with the potential to cause fatigue damage to the structure. Due to its unique structure and working principle, Kaplan turbines are particularly sensitive to pressure pulsation. Hence, the examination of pressure pulsation in these turbines is imperative to ensure the stability of unit operations.
Several scholars have conducted in-depth investigations into pressure pulsation by integrating experimental measurements with numerical simulations. Experimental studies have yielded valuable data under actual operating conditions, while numerical simulations have offered a more comprehensive understanding of flow dynamics and the ability to forecast pressure pulsations across various operational scenarios [
10]. Minakov et al. [
11] employed the Reynolds Mean Navier–Stokes (RANS) model, separate eddy simulation model (DES), and large eddy simulation model (LES) in their numerical simulations and laboratory experiments on a Kaplan turbine runner model. In contrast to the Limited Eddy Viscosity Model (LEVM), the RSM, DES, and LES could accurately replicate the intensity of pressure pulsations. Liu et al. [
12] utilized a three-dimensional full-channel unsteady turbulence method incorporating dynamic and static interaction to predict the pressure pulsation of a Kaplan turbine runner. They conducted tests on a test platform to compare the calculation results and found that the largest low-frequency pressure pulsation was generated by the draft tube. Rivetti et al. [
13] employed a combination of computational simulation and empirical measurements to investigate the pressure pulsation in a large Kaplan runner, encompassing the entire prototype. Pressure sensors were positioned at both the guide vane outlet and the draft pipe inlet. The findings revealed that the mean pressure exhibited fluctuations corresponding to variations in the guide vane channel within the rotating reference frame.
In the realm of numerical simulation, the characteristics of pressure pulsation have been analyzed, and the mechanism of pressure pulsation has been discussed. Wu et al. [
14] conducted a simulation of the unsteady flow in a prototype Kaplan turbine full channel by integrating the Reynolds-averaged Navier–Stokes equation with the RNG k-epsilon turbulence model. The simulation results revealed a sudden increase in pressure pulsation frequency in the draft tube, which may have been attributed to the resonance characteristics of the prototype. Yan et al. [
15] conducted an analysis of the internal flow characteristics of Kaplan turbine runner blades and draft tubes. The results revealed that when the guide vane opening was small, low-frequency and high-amplitude pressure pulsations occurred on the high- and low-pressure edges of blades. Additionally, the discrimination number was proposed as a quantitative measure for evaluating the flow uniformity of draft tubes. Altimemy et al. [
16] conducted a simulation of the internal flow structure of a Kaplan turbine operating at the peak of the designed flow rate in OpenFOAM. The simulation results indicated that the strong swirling flow away from the runner had a detrimental effect on the pressure pulsation. Luo et al. [
17] carried out a three-dimensional unsteady numerical analysis of a prototype Kaplan turbine operating under high head conditions. The findings revealed that the presence of eddy currents in the vaneless zone led to significant pressure pulsations when the guide vane opening was small under high head conditions. Iovanel et al. [
18] conducted a study on the limitations of unsteady Reynold-averaged Navier–Stokes (RANS) simulation in modeling internal flow within axial turbines. The study evaluated the classical turbulence model, the curvature correction model, and scaled the adaptive simulation-Shear Stress transport (SAS-SST) turbulence model. The findings revealed that while the pressure fluctuation frequency on the flow channel blade could be accurately captured, there was a significant underestimation of its amplitude.
For experimental measurement, Amiri et al. [
19] measured the unsteady pressure of blades and fixed components in the Kaplan turbine model during load change. It was discovered that the change of partial load caused the draft tube surge, forming a rotating vortex rope with two forms of fluctuation components, including rotation and plunging. This phenomenon caused flow oscillation throughout the entire flow passage, leading to local pressure fluctuations caused by the rotating vortex rope. Amiri et al. [
20] also conducted measurements on the unsteady pressure of the runner blades of the Kaplan turbine model. The findings revealed that the formation of a rotating vortex rope by the turbine during partial load operation was identified as one of the primary sources contributing to the oscillating force exert on the runner blades. This phenomenon subsequently led to pressure pulsations at two sub-synchronous frequencies within the runner. Dehkharqani et al. [
21] conducted experiments on the prototype Kaplan turbine during start-up, focusing on measuring the unsteady pressure and strain of the runner blade and shaft, as well as the bending and torsional strain of the shaft. The results indicated that low-frequency pressure fluctuations of the runner blade occurred after the guide vane was opened from its completely closed position.
In an effort to enhance the mitigation of pressure pulsation and propose methods for controlling fluctuations, Starecek et al. [
22] developed an approach to enhance operational stability in Kaplan turbines. The design incorporated six non-uniform blade flow channels, each featuring varying axial distances between blades and guide vanes. This configuration effectively mitigated pressure pulsation resulting from the interaction between the rib and the runner. Angulo et al. [
23] conducted air injection experiments on both a prototype and a model type of Kaplan turbine in order to investigate the impact of varying jet flow velocities on mitigating pressure pulsation in the draft tube and exhaust ring. The findings revealed that air injection effectively reduced pressure pulsation in the prototype experiment, while the extent of efficiency decrease in the model experiment was found to be overestimated. Rivetti et al. [
24] conducted air injection experiments on a Kaplan turbine runner at various loads and net heads, and they observed that within a specific range, the level of pulsation decreased in direct correlation with the rate of air flow.
Based on previous studies, the vibration and internal flow characteristics of Kaplan turbines during operation have a significant impact on the stable operation and safety of the turbines. However, there is a lack of comprehensive research combining experimental measurement and numerical simulation in this area. This paper deeply investigates the vibration and internal flow characteristics of Kaplan turbines using a combination of field tests, sensor monitoring, data analysis, and the computational fluid dynamics (CFD) method. The goal is to reveal the law and influencing factors of vibration and pressure fluctuation, providing theoretical basis and technical support for optimizing turbine design and improving operation stability.
4. Conclusions
(1) In the vibration test of the head cover and the upper and lower frame, it was observed that at a blade angle of 35°, the vibration displacement of the head cover in both the X and Y directions initially showed a slight decrease, followed by an increase as the guide vane angle increased. However, there was only a small increase in vibration when the guide vane angle was between 44.3° and 49.8°. On the other hand, at a blade angle of 35°, there was an overall increase in vibration of the head cover in both X and Y directions with an increase in the guide vane angle. Additionally, it was found that regardless of whether the blade angle was 35° or 55°, the vibration displacement of the head cover in X direction exceeded that in Y direction.
(2) When the blade angle was 35° and 55°, the vibration displacement of the upper frame in the X and Y directions surpassed that of the lower frame. Furthermore, it was observed that the vibration displacement of the upper frame in the X direction exceeded that of the lower frame in the Y direction. Interestingly, there was little disparity between the vibration displacement of the upper and lower frames in the Y direction. When the blade angle was set at 35°, the vibration displacement of the upper frame initially decreased and then increased as the guide blade angle widened. In contrast, the vibration displacement of the lower frame first increased and then stabilized as the guide blade angle expanded. However, when the blade angle was adjusted to 55°, both the upper and lower frame experienced an increase in vibration displacement as the guide blade angle widened, albeit at a decreasing rate.
(3) When the blade angle was 35° and 55°, five representative working conditions covering the experimental guide blade angle range were selected for numerical simulation. A total of 10 working conditions were carried out to cover this range. Due to the regulating effect of the guide vane on the direction of water flow, the water flow into the vaneless zone became more uniform and curved. As the blade angle increased, both the velocity in the vaneless zone and runner also increased. The velocity distribution tended to become more uniform with an increase in the guide vane angle. However, under conditions of a small guide vane angle, vortices appeared near the high-pressure edge of the runner in the vaneless zone, leading to relatively disorganized flow lines in that area. In terms of frequency components in pressure pulsation frequency domain at these vortices within the vaneless zone, multiple frequencies appeared due to their complex behavior. For a blade angle of 35°, dominant frequencies in the vaneless zone were either 1 time the blade passing frequency or runner rotating frequency. When BA = 55°, regardless of different guide vane angles, the main frequency remained at 1 time the blade passing frequency.
(4) When the guide vane angle was small (BA = 35°, GVA = 37.6°; BA = 55°, GVA = 42.9°, 46.9°), the frequency components of the suction surface and pressure surface of the runner blade became more complex, with the most intense pressure pulsation intensity. In other working conditions, the main frequency of the suction surface and pressure surface was either 1 or close to 1 time the runner rotating frequency, and the pressure pulsation intensity was relatively weak when there was a large guide vane angle in the working condition.
(5) The velocity flow line of the guide vane exhibited a relatively linear and concentrated trajectory, with a noticeable velocity gradient from the inlet to the outlet. As the blade angle increased, there was a gradual rise in speed. When the guide vane angle was small (BA = 35°, GVA = 37.6°; BA = 55°, GVA = 42.9°, 46.9°), a portion of the fluid’s pressure energy at the inlet was transformed into kinetic energy, resulting in an accelerated flow out of the guide vane and consequently leading to higher outlet velocities for the guide vane. These phenomena were accompanied by multiple low-frequency components characterized by varying amplitudes. In other operational conditions, the main frequency of the guide vane was 1 blade frequency, and the secondary frequency was mostly 4 times the transfer frequency.
(6) One of the most critical limitations of the current study is the lack of direct comparisons between experiments and numerical simulations. The CFD calculations only provide pressure results. Future studies could aim to better correlate vibration measurements with CFD pressure through the fluid–structure coupling (FSI) method and investigate the consequences of different combinations of blade angle and guide blade angle (e.g., structural deviation, stress, fatigue, etc.).