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.

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.

To accelerate the integration of fluctuating renewable energy technologies in the power systems, it is necessary to increase the flexibility of hydropower by operating turbines at off-design conditions. Unfortunately, this strategy causes deleterious flow phenomena such as von Kármán’s vortices at the wake of the vanes and/or blades. When their shedding frequency lies in the vicinity of a structure’s natural frequency, lock-in occurs and vibration amplitudes increase significantly. Moreover, if cavitation occurs at the centers of these vortices, the structure’s dynamic response will be modified. In order to understand this interaction and to avoid its negative consequences, the vibration behavior of a NACA 0009 hydrofoil under a torsional lock-in condition was numerically simulated for cavitation-free and cavitating-flow conditions. The results showed that the presence of vortex cavitation modified the formation and growth process of shed von Kármán vortices in the near-wake region which, in turn, caused an increase of the work performed by the hydrofoil deformation on the surrounding flow and a sharp decrease of the maximum vibration amplitude under resonance conditions.

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.