Authors: Linmin Li, Yabiao Niu, Guolai Wei, Sivakumar Manickam, Xun Sun, Zuchao Zhu
Categories: Original Research Article, Cavitation, Eulerian-Lagrangian, Large eddy simulation, Multiscale modeling, Noise characteristics
Source: Ultrasonics Sonochemistry
Keywords: Cavitation, Noise characteristics, Eulerian-Lagrangian, Multiscale modeling, Large eddy simulation
Hydrodynamic cavitation is a common phenomenon in fluid machinery [1], [2], [3], [4] which primarily occurs at locations where local pressures are lower than saturated vapour pressures, resulting in several adverse effects, like pressure loss, erosion, vibration, and noise [5], [6], [55]. For example, cavitation noise caused by propellers has emerged as a major source of underwater radiated noise emanating from ships, posing a serious threat to the marine environment [7], [8], [56], [63]. On the contrary, hydrodynamic cavitation can be utilized to enhance the reaction process of various chemical and environmental applications [58], [59], [60], [62]. Full-scale measurements are the most accurate method of measuring noise; however, these measurements are limited by experimental facilities and scaling assumptions. It is, therefore, imperative to develop practical computational methods for predicting noise signatures that address these concerns.
There has been a continuous development of numerical models over the last few years, and the noise characteristics of cavitation have attracted increasing attention [9], [10], [11], [12], [13], [61]. A method of inverse tracing was developed by Kao [14] to address radiated noise resulting from unsteady cavitation of propeller blades. A significant difference was observed between the predicted and measured radiated noise due to an overestimation of the size of the blade cavity. A computational fluid dynamics (CFD) method has been demonstrated by McIntyre et al. [15] to accurately predict the qualitative trend of sound spectra produced by propeller cavitation. Nevertheless, the Schnerr-Sauer cavitation model overestimated the amount of radiated noise, indicating that the cavitation models need to be improved. Dong et al. [16] investigated composite hydrofoils' dynamic response and acoustic properties using cavitation-induced vibrations. In their studies, the composite materials were found to suppress cavitation, decrease pressure pulsations in the flow field, and increase lift-to-drag ratios. Kalikatzarakis et al. [17] developed a hybrid model that integrated physical and data-driven simulations to model propeller cavitation and noise. Using the semi-empirical formula, the near-field propeller cavitation noise prediction results were also relatively accurate [18]. Additionally, several studies have demonstrated that combining CFD and acoustic analogies based on the Williams and Hawkings (FW-H) equation is effective [19]. To investigate unsteady cavitation around a NACA0015 hydrofoil and its associated noise, Yu et al. [20] combined a filter-based turbulence model with the Zwart-Gerbera-Belamri (ZGB) cavitation model. To predict hydrodynamic noise, caused by tip-vortex cavitation (TVC) of submarine propellers, Ku et al. [21] employed delayed detached eddy simulation (DDES) with adaptive mesh refinement and the FW-H equation, and concluded that a higher pitch angle resulted in more noise. It has been reported by Testa [22] that an improved hydroacoustic formulation based on the standard FW-H equation can accurately predict noise magnitude and direction based on the emitted noise and sheet cavitation dynamics. An investigation by Sezen et al. [23] examined the influence of biofouling roughness on propeller cavitation and underwater radiation noise using the Reynolds-averaged Navier-Stokes (RANS) solver coupled with the porous FW-H equation.
The combination of CFD and FW-H equations can effectively predict hydrodynamic noise caused by cavitation, but accurate prediction of the noise spectrum still has some limitations. In an attempt to predict the acoustic behaviour of TVC, Sezan et al. [24] combined the RANS-DES approach with the porous FW-H equation. RANS solver prediction was found to be insufficient compared to DES solver prediction. Using the multiphase DDES model, Li et al. [25] investigated the underwater radiated noise of a ship operating at its design speed, and they found that the method underestimated the levels of broadband noise in the frequency range of 50–112 Hz, where the contribution of TVC was significant. The maximum under-prediction in this range was about 28 dB at 72 Hz, and at frequencies above 200 Hz, the broadband noise became more and more under-predicted with increasing frequency. Lidtke et al. [26] performed scale-resolving CFD simulations to predict cavitating propeller noise in off-design conditions. They found that accurately modeling a wide range of frequencies, including tonal and broadband components, was still challenging for scale-resolving simulations. This is due to the need for high spatial resolution when resolving the small length scales involved in high-frequency noise generation mechanisms.
Klapwijk et al. [27] investigated the effect of synthetic inflow turbulence on error estimates under steady inflow conditions using a wing in wet and cavitating conditions. An Eulerian–Lagrangian one-way coupled approach was used by Ku et al. [28] to study the inception of TVCs and the noise produced by underwater submarine propellers. The model's validity was verified by comparing the predicted acoustic pressure spectrum with the measured one. Lidtke et al*.* [29] simulated sheet cavitation and its noise on a NACA0009 hydrofoil using the large eddy simulation (LES) method with the Schnerr-Sauer cavitation model. As a result, they found that the radiated noise correlated with the second temporal derivative of the total cavitation volume and fluctuations in load.
Moreover, it has been demonstrated that it is impossible to predict high-frequency noise without considering small-scale bubble structures. Additionally, Zhang et al. [30] found that the noise predictions made by the LES and FW-H hybrid methods and the LES and Powell vortex sound theory were similar. The numerical research regarding cavitation noise has demonstrated the possibility of accurately simulating the noise induced by resolvable structures. However, high-frequency noise generated by small-scale bubble structures cannot be accurately predicted.
Research has demonstrated that a multiscale approach incorporating Lagrangian particle tracking and interface-resolving methods [31], [32], [33], [34] can effectively resolve large- and small-scale phase structures. Ling et al. [34] proposed a multiscale model in which the large-scale interfaces were resolved using the volume of fluid (VOF) method to investigate droplet formation in atomization. The small droplets were tracked using a Lagrangian point-particle model. Li et al. [31] achieved a discrete-continuum transition algorithm for simulating the large and small bubbles and their transformation in a gas-stirring process, where small bubbles were tracked using the Lagrangian method and large bubbles were resolved by VOF. The multiscale bubble distribution was well characterized.
To improve the prediction of hydrodynamic cavitation, Xu et al. [35] and Cheng et al. [36] modified the Schnerr-Sauer cavitation model to enhance the TVC prediction by taking the effect of non-condensable gas into account, the gas bubbles were added in the flow field and tracked using the Lagrangian approach to calculate the local total bubble pressure for calculating the local mass transportation. Li et al. [33], [37] proposed a multiscale simulation algorithm to reproduce the complex phase structure of tip-leakage cavitating flows and analyzed the impact of model parameters using the coupled Eulerian-Lagrangian method with considering the transformation between large cavities and small bubbles. Furthermore, Yakubov et al. [38] proposed and validated a hybrid MPI/OpenMP algorithm for parallelising the Eulerian-Lagrangian approach to simulate cavitating flows and improve the computation efficiency. Based on the collapse dynamics of Lagrangian bubbles, Peter and Moctar [39] examined cavitation-induced erosion by applying a multiscale Eulerian-Lagrangian model. An investigation of the flow around a sharp-edged bluff body was conducted by Ghahramani et al. [40] using a hybrid cavitation model, which combined a mixture model with a Lagrangian bubble model. They have pointed out that the correct estimation of small-scale cavities can be as important as that of large-scale structures, because the cavitation noise, erosion and strong vibrations occur at small time and length scales.
Specifically, this work was focused on developing a multiscale model, to replicate the cavitation microstructure, and further to overcome the macroscale model's limitations on predicting cavitation noise. Our work explores the multiscale cavitation phenomenon and its associated noise around a hydrofoil that is part of the NACA66 (National Advisory Committee for Aeronautics). The Eulerian-Lagrangian multiscale cavitation model and the LES approach were employed for the transient cavitating flow, and the noise was predicted using the FW-H acoustic analogy equation. By comparing the experiments with the simulations, we observed that the multiscale approach performs better than the VOF method in predicting multiscale features. Lastly, the relationship between the evolution of the multiscale cavitation structure and the noise signal was examined.
A VOF method can be used for simulating the macroscale morphology of cavitation. As a result, the governing equations can be expressed as follows (Eqs. (1), (2)):
where F→s is the surface tension [41] and F→i is the resultant force interacting with the microbubbles for the momentum coupling. ρm and μm represent the density and viscosity of the mixture, respectively, which are defined as follows (Eqs. (3), (4)):
The mass transfer rate m˙ is given by the Schnerr–Sauer cavitation model as follows (Eq. (5)) [42]:
where m˙e and m˙c are the evaporation and condensation terms, respectively. Ce and Cc are the coefficients of evaporation and condensation fixed at 1.0 and 0.2, respectively [43], and pv is the saturated vapour pressure set to 3169 Pa. RB is the bubble radius represented as follows (Eq. (6)):
where nB is the bubble number density, which is 1 × 10^11^ [43], [57].
Under the current computing capability, the microscale bubbles in cavitating flows cannot be directly resolved by the VOF method due to its dependence on mesh resolution. Microscale bubbles can be tracked using a Lagrangian tracking method called the discrete bubble model (DBM), represented as follows (Eq. (7)):
where F→Drag, F→VM, and F→PG are the drag force, virtual-mass force, and pressure gradient force, respectively, represented as follows (Eqs. 8–10):
where dB represents the discrete bubble diameter, CD is the drag coefficient [44], and CVM is the coefficient of virtual-mass force (equals 0.5).
To model the nucleation process, cavitation cores are randomly assigned in a region where the pressure is below the nucleation pressure to simulate the growth and collapse of discrete bubbles. Zwart-Gerber-Belamri's cavitation model [45] assumes an initial bubble radius of 1 × 10^-6^ m. Furthermore, using the simplified Rayleigh-Plesset equation, we can simulate the size evolution of discrete bubbles as follows (Eq. (11)):
where Csp, Cevp and Ccond are the coefficients of saturated pressure, evaporation and condensation in the discrete bubble model, and have values of 1.25, 0.5 and 0.01, respectively [37].
The volume fraction of discrete bubbles is evaluated by integrating the volume of all bubbles in each cell and the total volume fraction including both macroscale and microscale parts is calculated as (Eq. (12)):
Then the total volume of cavity V is calculated by volume integrating in the flow field as (Eq. (13)):
where Vcell represents the cell volume.
LES can simulate large-scale vortices as well as model small-scale vortices directly. An LES model can reproduce the detailed information of unsteady flows compared to a RANS model. The mass and momentum conservation equations for vapour/two-phase liquid flows are shown below (Eqs. (14), (15)).
where u¯i and u¯j are the filtered velocities in the i and j directions, respectively, F is the force acting on the fluid including the interacting force from discrete bubbles, and τij is the sub-grid scale (SGS) stress term, which is defined as follows (Eq. (16)):
The eddy-viscosity model is expressed as follows (Eq. (17)):
In previous studies [46], the wall-adapting local eddy-viscosity (WALE) model was applied, and the SGS stress was expressed as follows (Eq. (18)):
where Cw is the empirical coefficient of the WALE model equal to 0.325, and the Sijd and S¯ij are defined as follows (Eqs. (19), (20)):
An acoustic analogy is widely used for engineering simulations and applications based on the FW-H equation, which is extended from Lighthill's acoustic analogy [47] and developed by Ffowcs-Williams and Hawkings [19] to incorporate the influence of surfaces in arbitrary motion. Using integrals of source field quantities, the sound at any arbitrary location can be determined. Di Francescantonio et al. [48] presented a new boundary integral formulation to evaluate the radiated noise, which permits the selection of permeable surfaces within the fluid domain to integrate. By integrating the permeable surfaces, noise sources in the volume enclosed by the surfaces are accounted for without performing volume integration [49]. This method has been proven to be effective in hydroacoustic analysis for cavitating flows [9], [50], and is also employed in the present work.
The porous FW-H equation can be derived from the Navier-Stokes equations and continuity equation as follows (Eq. (21)):
where p′ is defined as p′=p-p0, which is the sound pressure in the far field. δf represents the *Dirac-*δ function. Pij and Tij represent the compressive stress tensor and Lighthill’s stress tensor, respectively. c represents the sound speed and ρ0 is the ambient density.
A solution to Eq. (21) can be expressed as the sum of the three terms at any arbitrary receiver x as follows (Eq. (22)):
The three terms on the right-hand side of Eq. (22) are represented as pT′x,t is the thickness of the noise source (monopole). Cavitation development and collapse are responsible for the thickness acoustic source. Alternatively, it is referred to as cavitation noise. pL′x,t is the loading (dipole) noise source. In addition to low-frequency noise, it is an unsteady pulse force generated by turbulent and inhomogeneous flow fields on the hydrofoil surface. pQ′x,t describes the quadrupole noise resulting from the hydrofoil's boundary layer and vortex shedding, which is a high-frequency noise.
Following Eqs. (23), (24), (25), (26) consider the intermediate acoustic variables.
For a stationary source surface S, pT′x,t and pL′x,t are expressed as follows (Eqs. (27), (28)):
where ...ret denotes terms evaluated at the retarded time. Since the receiver locations are far from the source surface S, pT′x,t can be expressed as follows (Eq. (29)):
where ret¯ and r¯ represent the average retarded time and the average distance of all noise sources on the surface S, respectively, and V represents the cavity volume in the source surface S [51]. Thus, pT′x,t is proportional to the second derivative of the cavity volume.
Fig. 1 illustrates the computational domain and boundary conditions. Dimensions of the domain are 850 mm × 75 mm × 250 mm. This hydrofoil (NACA66) has a chord length (C) of 100 mm, an angle of attack of 8°, and a span length of 0.75C. Hydrofoils are positioned at the center in the direction of height. From the leading edge, the inlet and outlet of the calculation domain are situated at 3.5C and 5C, respectively. In this case, the inlet velocity is 10 m/s, and the outlet pressure is based on the cavitation number. Validation is provided by the experiment conducted by Wu et al. [52], whereas the simulation is conducted with the condition of σ = 1.5. To ignore the wall effect, the hydrofoil surfaces are designed as no-slip walls, and the side walls are designed as free-sliding walls.
Fig. 1 Computational domain and boundary conditions.
In the current model, the computational domain is discretised into structured hexahedron cells, as shown in Fig. 2. During the refinement procedure, the mesh at the cavitation areas is ensured to have acceptable orthogonal quality and improve the simulation accuracy. Approximately 100 × 120 meshes are located adjacent to the hydrofoil walls. The first layer's cell height is 0.012 mm, and its growth ratio is 1. A refined mesh is used near the hydrofoil's leading edge and tail to capture the cavitation process's intricate details.
Fig. 2 The distribution of meshes in the computational domain for (a) side walls and (b) hydrofoils.
In the present work, the momentum equation is solved using the bounded central differencing spatial discretisation algorithm, and the bounded second-order implicit time formulation is utilised for the LES method. To enhance the simulation stability, the implicit scheme is used for solving the equation of the volume fraction in the VOF method, which is discretised as (Eq. (30)):
where n is the index for the previous time step, U is the volume flux through the face and the subscript f represents the face value.
For pressure–velocity coupling, the pressure-implicit with splitting of operators (PISO) scheme, which is suitable for transient simulations, is utilised. In this scheme, the momentum equation using the operator H for the finite-difference representation is first solved using the pressure field Pn from the old time step as (Eq. (31)):
where ρui∗ is the flux calculated by the intermediate velocity field ui∗ which may not satisfy the zero-divergence condition. So a corrector step is performed using a new velocity field ui∗∗ and a corresponding new pressure field P∗ obtained by the pressure equation seeking the zero-divergence condition, which is taken as (Eqs. (32), (33)):
A further corrector step can be introduced to obtain ui∗∗∗ and P∗∗. By conducting corrector steps until a pre-defined tolerance is met, the new velocity can be consistent with the new pressure.
Moreover, two values of time step (1 × 10^-5^ s and 5 × 10^-6^ s) are set to test the time step independency, the convergence criteria for all variables are set as 1 × 10^-4^. The computational code for VOF-DBM coupled simulations is developed using user defined functions (UDF) embed into the FLUENT software. The VOF and DBM equations are solved for continuous vapour phase and discrete bubbles, respectively. In the region with vapour volume fraction larger than 0.6 (as the maximum volume fraction of discrete phase is about 0.6), the discrete bubbles are removed and the mass transfer rate is calculated using the source term in Eq. (5). Otherwise, the growing and collapsing terms of discrete bubbles in the dispersed region are modelled using Eq. (11). The algorithm for transformation between VOF and DBM [37] is also implemented for bubbles to grow into continuous ones and large cavities to collapse into discrete bubbles. Thus, the coupling between VOF and DBM is achieved.
For the FW-H approach to calculate the far-field sound pressure signal, noise integration surfaces must be established for tracking velocity and pressure variations. The setup location of the noise integration surface is shown in Fig. 3. To ensure that the cavitation area does not extend beyond the integration surface when the cavity drops off, it is necessary to set the noise integration surface as large as possible to cover the whole cavitation area. Further, the locations of the noise receivers are indicated in Fig. 3 as R1-5. Several noise receiving points are set up around the hydrofoil at varying distances to measure sound pressure distribution in the same direction. For compact sound sources, the receivers and integration surfaces must satisfy the Minimum Frequency Parameter (MFP) criteria [53].
Fig. 3 A diagrammatic representation of the FW-H integration surface and the locations of the noise receivers.
An investigation of the effect of mesh resolution on simulation results is conducted by comparing three coarse, medium, and fine. Fig. 4 shows the distribution of y^+^ values on the hydrofoil surface using different meshes. Based on the y^+^ values (0 ∼ 3) at the surface of the hydrofoil, it is evident that the LES requirements for simulations are generally met. An overview of the parameters used for mesh generation and the comparison of average drag and lift coefficients are presented in Table 1. Observations show that the coarse mesh simulations of Cl and Cd have large errors. In contrast, those of medium and fine meshes are closer to each other. Further simulations use the medium mesh to balance computing demands and simulation accuracy.
Fig. 4 Y^+^ values for (a) coarse, (b) medium, and (c) fine meshes.
Fig. 5 illustrates the temporal evolution of the lift coefficient Cl and drag coefficient Cd for a hydrofoil simulated under the fine and medium meshes. Lift and drag coefficients fluctuate periodically as a result of cavitation shedding. Fig. 6 shows the power spectrum density (PSD) of Cd under the medium and fine meshes. In this case, we find that the PSD of Cd follows the −7/3 power law. As shown in Fig. 6**,** the main frequency of the fluctuation is 45.58 Hz, while the corresponding period is 21.93 ms**.** According to the experimental data [52], the period of cavitation shedding is 22.29 ms. The discrepancy between the simulation and the experiment may result from the 1 mm gap between the hydrofoil and the side wall, as well as from the roughness of the side wall (not considered in the simulation). The difference between the experimental and simulation cavitation shedding periods is approximately 1.6%, indicating that the simulation is accurate.
Fig. 5 Temporal evolution of (a) lift coefficient C
land (b) drag coefficient Cdfor hydrofoils with fine and medium meshes.
Fig. 6 Power spectral density (PSD) of C
d.
To investigate the sensitivity of the time step, Fig. 7 shows the results of sensitivity analysis of the time step performed in the present work. From both the curves of lift coefficients and the cavitation features in a typical cycle, it can be found that the two results using 0.01 ms and 0.005 ms are close to each other. The results indicate that the time step used in the present simulation is adequate as the implicit scheme for VOF allows a relatively large time step.
Fig. 7 Evolution of (a) C
land (b) cavitation features using different time steps.
Fig. 8 compares cavitation characteristics observed experimentally and predicted by numerical simulation using the macroscale (VOF) and multiscale (VOF-DBM) approaches, with a cavitation number of 1.5. The surface of hydrofoil experiences periodic cycles of cavitation development, shedding, and collapse. It is evident from the evolution of the cavity volume that cavitation is shed periodically. During each cycle of cavity evolution, four stages can be the development of sheet cavitation, the cutoff and collapse of sheet cavitation, the shedding of cloud cavitation, and the collapse of cloud cavitation. As shown in Fig. 8**,** the sheet cavitation develops from the leading edge to the trailing edge and grows to a location where X/C = 0.8 from 0/10 T to 6/10 T. Cavitation generates re-entrant jets at the tail of the cavity, which travels from the sides of the trailing edge toward the center of the leading edge. As illustrated at 9/10 T in Fig. 8**,** when the re-entrant jets converge at the hydrofoil's midsection, the cavity's tail cuts off and falls off on the downstream side. A re-entrant jet caused by an adverse pressure gradient in the downstream flow field is responsible for causing the cavity shedding process. VOF-DBM results in comparison to VOF results reveal numerous discrete bubbles around the trailing edge of the hydrofoil. Experimental observations support the location of discrete bubbles. Accordingly, the VOF-DBM method can accurately predict the transition between large cavities and microscale bubbles in the cavitating flows. Compared to VOF results, these results are more consistent with the experiment.
Fig. 8 Comparison of cavitation characteristics between (a) experimental results [52] and simulations predicted by (b) macroscale (VOF) and (c) multiscale (VOF-DBM) (σ = 1.5, α = 8°, U~∞~ = 10 m/s).
Fig. 9 illustrates the evolution of cavity volume Vvap over several typical cycles calculated using VOF and VOF-DBM methods. The results show that the cavity volume predicted by the VOF-DBM method is larger than that predicted by the VOF method; however, the variation period of the cavity volume remains unchanged. In particular, the cavity volume calculated by the VOF-DBM method is approximately 20% larger than that calculated by the VOF method when discrete bubbles are considered. Accordingly, it can be concluded that the macroscale VOF method cannot capture the distribution and transport of discrete microbubbles at the current mesh resolution. Multiscale models, on the other hand, are capable of acquiring these microscale structures.
Fig. 9 Typical cycles of vapor volume evolution.
Using the VOF and VOF-DBM methods, Fig. 10 illustrates the pressure distribution on the surface of the hydrofoil as well as the cavity shapes at two representative instants. In Fig. 10**(b)**, the black dotted circle indicates that VOF-DBM predicts greater low-pressure areas near the trailing edge than VOF. The multiscale model also captures microscale bubbles; negative pressure should accompany bubble transport. Furthermore, the multiscale model can capture the pressure surge caused by microbubble collapse. Based on this study, the VOF-DBM method can capture more detailed characteristics of the transient cavitating flow, including the microscale structures and their effects on the flow field.
Fig. 10 Comparison of the pressure distribution on the hydrofoil surface and cavitation patterns in (a) macroscale (VOF) and (b) multiscale (VOF-DBM) simulations.
Acoustic signals generated by the flow around a hydrofoil, including cavitation noise, can be analysed using an acoustic integration surface around the hydrofoil. Fig. 11 illustrates the noise pressure signal p′ obtained from the monitoring point at R3 over several cavitation periods, where p′ represents the total noise according to Eq. (21). The total noise is composed of cavitation noise, dipole noise and quadrupole noise, as shown in Eq. (22). The maximum values of total noise calculated by both the VOF and VOF-DBM methods are approximately 200–300 Pa. Due to the inclusion of the sound pressure signal induced by microbubbles in the VOF-DBM method, the acoustic signal calculated by this method has a greater frequency and amplitude.
Fig. 11 Total sound pressure signal at R3 calculated using two (a) VOF and (b) VOF-DBM.
Fig. 12 illustrates the Sound Pressure Level (SPL) of p′ R3 based on the frequency spectrum obtained by both methods. In this study, 45.58 Hz was the dominant frequency of the total noise, consistent with the cavitation shedding frequency. This study illustrates that cavitation is the main noise source, while the attached cavity shedding dominates the maximum SPL. It is evident that the SPL calculated by the VOF method and that calculated by the VOF-DBM method are similar in the low-frequency band. Nevertheless, the SPL calculated by the VOF-DBM method is higher than that obtained by the VOF method in the high-frequency band, especially in the frequency band around 10 kHz, which illustrates that the frequency of cavitation noise caused by microbubbles is about 10 kHz.
Fig. 12 Calculation of SPL for R3 using the (a) VOF and (b) VOF-DBM methods.
The SPL measurements of the monitoring points at R1-5 located above the hydrofoil were plotted, as shown in Fig. 13**(a), to examine the differences in noise between different receiving points. In general, the SPL decreases as the distance from the hydrofoil increases. The dominant frequencies are the same in all monitoring points, and the dominant Strouhal number (fC/U∞) is 0.456. Moreover, we observed that the decaying slope of the sound signal spectrum, monitored at all points, falls between f
^-1^ ∼ f
^-2^, which is consistent with that described in the experimental investigation [54]. Fig. 13(b)** shows the evolution of overall SPL along with the logarithmic distance in the y direction. It illustrates that the SPL decreases linearly with the logarithmic distance. Furthermore, it is found that the noise is a continuous spectrum with signals ranging from low frequencies to high frequencies and the SPL at low frequencies is relatively high.
Fig. 13 (a) Vertical SPL for receivers at R1-5, (b) evolution of overall SPL at R1-5, along with logarithmic distance.
As shown in Fig. 14**,** the cavitation (thickness) noise signal at receiver R3 is represented by the VOF and the VOF-DBM methods, where pT′ represents the thickness noise or cavitation noise, as shown in Eq. (22). As compared to the total noise p′ in Fig. 11, the cavitation noise signal is approximately a third of the total noise pressure signal. A cavitation noise signal exhibits strong pulses in every cycle when its pT′ value is around the negative peak. In contrast, the strong pulses in the total noise signal when the p′ value is around the positive peak are not present in the cavitation noise signal, indicating that the strong pulses are not the result of thickness noise caused by the growth and collapse of a cavity. The results from VOF and VOF-DBM indicate that the VOF-DBM method can predict a greater proportion of high-frequency signals, indicating that microbubbles play a significant role in producing high-frequency noise. Unfortunately, the traditional VOF method does not consider microbubbles and other macroscale methods. Additionally, the fundamental frequency of the cavitation noise is the same as the cavitation shedding frequency. Investigating the relationship between the cavitation noise and the unsteady cavity features would be useful.
Fig. 14 Cavitation noise pressure signals at R3 calculated using two methods (a) VOF and (b) VOF-DBM.
The volume evolution of discrete bubbles is shown in Fig. 15**(a), it is found that the evolution period is also 21.93 ms, consistent with the cavitation's primary shedding period. As shown in Fig. 15(b)**, the Power Spectral Density (PSD) of the cavitation noise induced by discrete microbubbles is primarily characterised by high PSD values at high frequencies. In the process of collapse, the cavity transforms into microscale bubbles that are transported downstream, and the change in pressure around the discrete bubbles results in their collapse. Due to the high frequency of collapse caused by many microbubbles, the cavitation noise caused by microbubbles is also in the high-frequency band and cannot be ignored. The highest PSD of the cavitation noise generated by discrete microbubbles is approximately 10 kHz.
Fig. 15 (a) Discrete bubble volume evolution over four cycles and (b) PSD of noise by discrete bubbles.
As shown in Fig. 16**(a), the cavitation noise signal and the cavity volume evolution simulated using the VOF-DBM method are plotted to examine the relationship between the cavitation noise and the multiscale cavitation features.** As shown in Fig. 16**(b)**, the multiscale simulation and the experiment were compared at three different instants in the evolution curve plotting. Comparing the evolutions of the cavitation noise signal and the cavity volume, the maximum and minimum cavities correspond to the negative and positive peaks of the cavitation noise pressure signal, respectively. Cavitation noise consists of strong pulses when the cavity volume is large and begins to decrease. From this, the strong pulse of noise may have arisen due to the collapse of a large-scale cavity during the transport of a detached cavity.
Fig. 16 (a) Cavitation noise signal at R3 and volume evolution of the cavity over four cycles, and (b) the three marked instants from both simulation and experiment.
Radiated noise from cavitation is closely related to cavity structures, particularly for strong pulses resulting from cavity collapse. In order to study the details of cavitation features when the strong pulse occurs, the cavitation noise signal and the cavity volume evolution in one cycle are plotted together, and the cavitation features around the two corresponding instants are shown in Fig. 17. There are two signature cavitation noise pulses which are designated as I and II. There is a direct correlation between the intensity of cavitation noise and the acceleration of cavity volume. It can be observed from the curves that the positive cavitation noise pressure is associated with an increase in cavity growth rate. In contrast, the negative cavitation noise pressure is primarily associated with a decrease in cavity growth rate or an increase in cavity collapse rate, except for pulses. According to Fig. 16, the maximum cavity volume occurs when the attached cavity reaches its maximum size and the detached cavity transforms into a cloud of cavities. As seen in Fig. 17, the noise pulse arises after the moment of maximum cavity volume, which is because of the first collapse of the large cloud cavity, as seen in Fig. 17**(b)** where the cavitation structure evolution around instant I is shown for revealing the noise pulse mechanism. According to the cavitation structure evolution at instant II, some cavity clouds are found to grow and collapse again, resulting in a second strong noise pulse during the cloud cavity transportation. Compared to the first strong noise pulse, the second one is weaker. Furthermore, in addition to the two main noise pulses, other pulses are also found during other cavitation cycles, as the cloud cavity generates and collapses throughout the entire cavity evolution process.
Fig. 17 (a) Cavitation noise at R3 during a typical cycle and (b) simulated cavity structures at the same time.
This study uses the VOF-DBM coupled multiscale cavitation model to investigate the cavitation structures over a wide range of time and space scales and the corresponding cavitation noise characteristics on a NACA66 hydrofoil. The LES approach and FW-H equations are also employed to predict the noise pressure signals in the cavitating flow. VOF and VOF-DBM methods are compared in their ability to predict cavitation structures and noise characteristics. The noise characteristics induced by discrete microbubbles are taken into account, and the relationship between cavitation noise characteristics and multiscale cavity features is investigated. Several key conclusions can be summarised as
VOF-DBM can provide a more accurate indication of the transition between large cavities and microscale bubbles in cavitating flows compared to VOF-VOF. VOF-DBM simulation results are more consistent with the experiment than VOF results, as VOF-DBM can more accurately reproduce the bubble cloud transport in the wake region.
Main cavity shedding dominates the maximum SPL of the noise. In the low-frequency band, the VOF method provides similar results to the VOF-DBM method. In the high-frequency band, the SPL calculated by the VOF-DBM method is greater than that calculated by the VOF method. Induced by microbubbles, cavitation noise is concentrated around 10 kHz in frequency. Based on the simulation results in the high-frequency band, the multiscale model is more accurate at predicting the noise caused by cavitating flows.
According to the cavitation noise pressure signals, the maximum cavity volume corresponds to the negative peak, while the minimum cavity volume corresponds to the positive peak. During the collapse of the cloud cavity, a typical noise characteristic is observed, including two intense pulses. In most cases, the first pulse occurs when a large detached cavity collapses, which happens shortly after the attached cavity grows to its largest proportions. As a result of some cavity clouds growing and then collapsing again, the second pulse occurs.
Linmin Li: Investigation, Resources, Methodology, Writing – original draft. Yabiao Niu: Investigation, Visualization, Writing – original draft. Guolai Wei: Writing – review & editing. Sivakumar Manickam: Writing – review & editing. Xun Sun: Conceptualization, Writing – review & editing. Zuchao Zhu: Conceptualization, Supervision.
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
The authors are grateful for the support of the National Natural Science Foundation of China (Grant Nos. U21A20126, 52006197, 52276032), the Natural Science Foundation of Zhejiang Province (Grant No. LQ21E060012), and the Key Research and Development Program of Zhejiang Province (Grant No. 2022C01146). Thanks to Dr. Qin Wu from Beijing Institute of Technology for providing the experimental observations.
Xun Sun, Email: xunsun@sdu.edu.cn.
Zuchao Zhu, Email: zhuzuchao@zstu.edu.cn.