Linlin Geng

Periodic fluctuations in pump head are a common unsteady phenomenon in tubular pumps; however, their underlying energy dissipation mechanism remains insufficiently understood. This study conducted transient numerical simulations on a two-blade tubular-flow pump with highly twisted blades operating at the design flow condition. Using entropy production theory, the research quantitatively examined Rotor–Stator Interaction (RSI), vortex development, and hydraulic loss characteristics. The findings reveal that turbulent entropy production is the primary contributor to total energy dissipation, while entropy generated by wall friction is minimal. Although the impeller experiences the greatest absolute energy loss, the fluctuations in entropy production within the guide vane are significantly larger. This indicates that the energy loss represented by the entropy production in the guide vane primarily drives the periodic head fluctuations of the pump. High entropy production is concentrated near both the leading and trailing edges of the guide vane, exhibiting a trough-shaped radial distribution influenced by the leading-edge hub vortex and tip leakage vortex. Additionally, the transient changes in entropy production under RSI are governed by the periodic formation and strengthening of the guide vane passage vortex, along with the shedding and breakdown of the wake vortex. Velocity analysis shows that the circumferential movement of the impeller wake continuously modifies the instantaneous inflow conditions at the guide vane’s leading edge, causing periodic changes in the incidence angle that enhance the passage vortex while weakening the wake vortex. This study provides deep insights into the operational stability of pumps and pumping stations in water transfer projects.

The accuracy of the numerical simulation of unsteady cloud cavitation around a hydrofoil depends on the combination of the cavitation model, the Reynolds-Averaged Navier–Stokes (RANS) turbulence model and the exponential coefficient n of the Reboud’s correction. To assess the influence of such choices, three turbulence models, three cavitation models and two values of n have been combined to predict the shedding frequency and the morphology of cloud cavitation around the NACA65012 and NACA0009 hydrofoils. The comparison with analogous experimental results obtained in a cavitation tunnel indicates, firstly, that the same numerical set-up differs in accuracy depending on the hydrofoil geometry. Secondly, within the scope of the two tested hydrofoils and corresponding flow conditions, the Shear Stress Transport (SST) turbulence model and the value of n = 10 appear to be more accurate and more robust for all tested cases. And finally, the predicted shedding frequency is more sensitive to the selection of the turbulence model than to the cavitation model. If an erosion model is implemented, then it is found that the predicted potential energy distribution of the cavitating flow is sensitive to the selected cavitation model. In our case, the Sauer model gives a more accurate distribution and intensity of erosion power than the rest of the cavitation models.

A novel pumped storage system using centrifugal pumps to transfer water between reservoirs in coastal hydropower plants has significantly mitigated grid instability. However, frequent start–stop operations of large vertical centrifugal pumps, which serve as the core equipment, severely affect the operational stability of these systems. In this study, the intrinsic connection between the cavitating flow field and irreversible losses during the process was analyzed using the entropy production theory. The time–frequency characteristics of pressure pulsation in pump were analyzed by using the continuous wavelet transform. The results indicate that with the reduction in the flow rate and rotational speed, the sheet cavitation at the impeller inlet rapidly weakens until it vanishes. The cavity cavitation within the draft tube commences to emerge in the turbine mode. Separation vortices are formed due to the mismatch in the flow angle at the impeller outlet. These vortices induce local cavitation, causing both a rapid energy loss increase and high-amplitude, low-frequency pressure pulsations. During transient processes, flow instabilities induce high-amplitude, low-frequency pressure pulsations within the stay vane region, with maximum amplitude attained during runaway condition. The research results provide a theoretical foundation for the stable operation of centrifugal pumps.

To study the effects of fluid compressibility on the dynamics of a cavitating vortex street flow in a regime where the vortex shedding frequency increases as a result of the cavitation increase, the cavitating wake behind a wedge was simulated employing both incompressible and compressible solvers. To do this, a compressible cavitation model was implemented, modifying the Zwart-Gerber-Belamri (ZGB) incompressible solver and including a pressure limit and absorbing boundary conditions to prevent a non-physical pressure field. To validate the performance of the compressible model, preliminary simulations were carried out on a 1D Sod cavitating tube and the cavitating vortex shedding behind a circular body at laminar flow conditions. The results of the cavitating wake behind the wedge with the incompressible and the compressible solvers showed similar predictions in terms of pressure, vortex shedding frequency, and instantaneous and average vapor volume fraction profiles. In spite of this, differences were obtained in the energy content of the fluid force fluctuations on the body at higher frequencies, which appear to be better resolved and amplified when the compressibility model is considered. Overall, both solvers provided comparable results in terms of cavitation phenomena that are well aligned with experimental observations.

The cavitating wake behind a wedge has been simulated employing both incompressible and compressible solvers with the purpose of comparing their performance. To develop the compressible solver, the ZGB cavitation model for incompressible flows has been modified to consider the compressibility effects. First, the equations of state for the vapor and liquid phases have been implemented. Then, a pressure limit and absorbing boundary conditions have been set to prevent a non-physical pressure field. To test the developed compressible solver, a series of verification tests have been conducted comprising a 1D sod cavitating tube and the cavitating vortex shedding behind a circular body at laminar flow conditions. After that, the incompressible and the validated compressible solvers have been applied to simulate the cavitating wake behind a wedge. The obtained results show similar predictions with the two solvers in terms of pressure, vortex shedding frequency and instantaneous and average vapor volume fraction profiles. Moreover, it is interesting to note that the components of the fluid force fluctuations on the body at the lowest frequencies, as predicted by both solvers, are indistinguishable. However, the higher frequency fluctuations of the forces tend to be better resolved and amplified when the compressibility is considered. Overall, both flow solvers have provided comparable cavitation phenomena, yielding results well aligned with experimental observations.

The use of pumps as turbines has been gaining more and more attention in recent years. The present work mainly investigates the influence of blade wrap angle on the internal flow and pressure fluctuation characteristics of centrifugal pumps as turbines. Five different wrap angles (35°,45°, 55°, 65°, and 75°) for a forward-curved impeller were numerically analyzed under multiple operating conditions. The accuracy of numerical simulation was validated by experimental results. The results show that maximum efficiency is achieved with a blade wrap angle of 35°, and the highest efficiency flow point gradually decreases as the blade wrap angle increases. It is found by conducting entropy production theory analysis that the high-entropy production rate regions in PATs are concentrated in the volute tongue and impeller blade inlet regions, and that the entropy production rate at the impeller inlet region increases and then decreases as the blade wrap angle decreases. In addition, pressure pulsation was affected not only by dynamic and static interference but also by an irregular vortex around the impeller; its magnitude under Q_t is higher than 0.8Q_t for blade wrap angles of 55° and 75°. The primary frequency of pressure pulsation within the impeller is the axial frequency fn and its multiples, and the frequency with the largest amplitude is 3fn. The periodicity of vortices is closely related to the periodicity of pressure pulsation. And it is suggested that a PAT with a 35° blade wrap angle is advantageous for improving the stability of a turbine.

Compared to traditional centrifugal pumps, magnetic pumps are widely used in industries such as chemical, pharmaceutical, and petroleum due to their characteristics of leakage-free operation and the ability to transport toxic and corrosive fluids. However, the efficiency of magnetic pumps is relatively low. Improving the efficiency of pumps helps to reduce energy loss and lower industrial costs. In this study, a magnetic pump was chosen as the research subject. The study aims to improve the efficiency and stability of the magnetic pump by optimizing the impeller blades based on the load curve. A combined approach of a numerical simulation and experimental verification was used to investigate the impact of the anterior loading point (AL), posterior loading point (PL), and slope (SL) in the blade loading curve on the pump’s performance. The slope, which had the most significant impact on pump performance, was selected as the dependent variable to analyze the internal pressure pulsation and main shaft radial force of the magnetic pump. The research found that the hydraulic performance test results of the magnetic pump were in good agreement with the simulation results. When efficiency is used as the optimization objective, the anterior loading point should be moved as far back as possible, and the posterior loading point should be moved as far forward as possible. Through the study of internal pressure fluctuations and radial forces within the pump, the radial force distribution is sequentially as follows: the anterior loading method, posterior loading method, and middle loading method at a rated flow rate. The maximum pressure pulsation amplitude was near the volute casing diffuser area. Compared to the original pump, the optimized magnetic pump showed a 5.05% improvement in hydraulic efficiency under the rated conditions. This research contributes to enhancing the performance and efficiency of magnetic pumps, making them more suitable for various industrial applications.

The tip leakage vortex (TLV) has aroused great concern for turbomachine performance, stability and noise generation as well as cavitation erosion. To better understand structures and dynamics of the TLV, large-eddy simulation (LES) is coupled with a homogeneous cavitation model to simulate the cavitation flow around a NACA0009 hydrofoil with a given clearance. The numerical results are validated by comparisons with experimental measurements. The results demonstrate that the present LES can well predict the mean behavior of the TLV. By visualizing the mean streamlines and mean streamwise vorticity, it shows that the TLV dominates the end-wall vortex structures, and that the generation and evolution of the other vortices are found to be closely related to the development of the TLV. In addition, as the TLV moves downstream, it undergoes an interesting progression, i.e., the vortex core radius keeps increasing and the axial velocity of vortex center experiences a conversion from jet-like profile to wake-like profile.

Vortex cavitation can appear in the wake flow of hydrofoils, inducing unwanted consequences such as vibrations or unstable behaviors in hydraulic machinery and systems. To investigate the cavitation effects on hydrofoil vortex shedding, a numerical investigation of the flow around a 2D NACA0009 with a blunt trailing edge at free caviation conditions and at two degrees of cavitation developments has been carried out by means of the Zwart cavitation model and the LES WALE turbulence model which permits predicting the laminar to turbulent transition of the boundary layers. To analyze the dynamic behavior of the vortex shedding process and the coherent structures, two identification methods based on the Eulerian and Lagrangian reference frames have been applied to the simulated unsteady flow field. It is found that the cavitation occurrence in the wake significantly changes the main vortex shedding characteristics including the morphology of the vortices, the vortex formation length, the effective height of the near wake flow and the shedding frequency. The numerical results predict that the circular shape of the vortices changes to an elliptical one and that the vortex shedding frequency is significantly increased under cavitation conditions. The main reason for the frequency increase seems to be the reduction in the transverse separation between the upper and lower rows of vortices induced by the increase in the vortex formation length.

A numerical investigation of the erosion aggressiveness of leading edge unsteady cloud cavitation based on the energy balance approach has been carried out to ascertain the main damaging mechanisms and the influence of the free stream flow velocity. A systematic approach has permitted the determination of the influence of several parameters on the spatial and temporal distribution of the erosion results comprising the selection of the cavitation model and the collapse driving pressure. In particular, the Zwart, Sauer and Kunz cavitation models have been compared as well as the use of instantaneous versus average pressure values. The numerical results have been compared against a series of experimental results obtained from pitting tests on copper and stainless steel specimens. Several cavitation erosion indicators have been defined and their accuracy to predict the experimental observations has been assessed and confirmed when using a material-dependent damaging threshold level. In summary, the use of the average pressure levels during a sufficient number of simulated shedding cycles combined with the Sauer cavitation model are the recommended parameters to achieve reliable results that reproduce the main erosion mechanisms found in cloud cavitation. Moreover, the proposed erosion indicators follow a power law as a function of the free stream flow velocity with exponents ranging from 3 to 5 depending on their definition.