Sixth International Symposium on Cavitation

6 downloads 8279 Views 635KB Size Report
cruising at service speed is presented in Fig. 6. ..... Brubakk, E., Smogeli, H., “QE2 from turbine to diesel - consequences for noise and ... Fluid Mechanics, Vol. ... Sponagle, N.C., “Noise from tip vortex and bubble cavitation”, DREA Tech.
International Shipbuilding Progress, Vol. 54, Number 2-3, 2007

Achievements and challenges in cavitation research on ship propellers Tom van Terwisga, Erik van Wijngaarden, Johan Bosschers and Gert Kuiper Maritime Research Institute Netherlands, Wageningen, The Netherlands E-mail: [email protected]

This paper reviews cavitation research on ship propellers performed in the recent past in MARIN’s Depressurized Towing Tank and cavitation tunnels. Cases are presented through which advances have been made possible in the design of propellers and ships, and unresolved issues are discussed that form the challenges of future research. Our focus is on inception, hull-pressure fluctuations, and erosion on propeller blades and rudders.

1.

Introduction

Cavitation occurs on nearly all ship propellers. It may lead to expensive problems if not acknowledged in an early design stage. The two most frequently occurring problems are noise and vibration in the afterbody due to cavitation-induced pressure fluctuations on the hull, and cavitation erosion on propeller blades and appendages. Early recognition of these adverse effects is important, not only to ensure compliance with contract requirements, but also because often cavitation has to be controlled at the cost of propeller efficiency. To ensure that the propulsor meets the requirements relating to comfort (vibration and noise) and safe and economic operation (erosion), model scale experiments or computations addressing cavitation are conducted prior to construction. As reliable computational techniques that incorporate all aspects of cavitation that relate to noise, vibration and erosion are lacking, scale model experiments are still the preferred choice. From experiments in the first cavitation tunnel, built by Parsons in 1895, insights were built up as to how one should model cavitation phenomena to allow for full scale predictions. This has led to the use of medium sized cavitation tunnels in the mid nineteen hundreds and eventually to the large cavitation tunnels of today which allow for complete ship models to be installed in front of the propeller. At the Maritime Research Institute Netherlands (MARIN) propeller cavitation behavior is investigated in the Depressurized Towing Tank (DTT) which allows for the testing of large scale ship models in the presence of a free surface. The DTT is in operation since 1972 [19] and has been modernized in 2001. A medium sized cavitation tunnel and a small high-

1

speed cavitation tunnel are also available. They are primarily used for background research on cavitating propellers and wing type geometries. Reliable predictions of full scale cavitation behavior have become more important with the ever increasing powers that must be absorbed by a single propeller. This necessitated the careful validation of predictions made by cavitation laboratories through full scale measurements and observations. This paper addresses a number of such validation studies and argues that despite the undeniable success of cavitation facilities around the world in qualitatively aiding in the design process, the reliability of predictions is not yet such that problems are avoided at all times. It must be understood that propeller designers, in order to meet the often contradicting requirements of high comfort and propulsive efficiency, are in need of a more quantitative improvement of prediction capabilities. Whilst showing areas where successes are scored, the paper also indicates topics where further research is necessary. In brief, cavitation research has allowed for today’s greatly improved design of propellers and afterbodies, resulting in much reduced levels of noise and vibration. Also our understanding of the phenomenological nature of cavitation erosion has greatly improved. Nevertheless, the state-of-the-art is such that several research issues are still outstanding and hamper our capability to predict cavitation behavior. In this respect one should think of hull-pressure fluctuations at orders of blade rate, which are observed to be sensitive with respect to environmental conditions, leading to varying accuracy in predictions. A similar statement holds for low frequency broadband pressure fluctuations. Also worth mentioning are the practical prediction of cavitation erosion due to a lack of detailed understanding of the erosion mechanism and the delay in inception of certain types of cavitation on model scale. Each of the main topics identified above is treated separately, starting with inception, then followed by hull-pressure fluctuations, and finally erosion.

2.

Cavitation inception

Cavitation inception is of direct importance to Navy vessels, because of the sudden increase in noise levels at the onset of cavitation. Indirectly, it is also important because a delay in cavitation inception is likely to affect predictions of erosion and radiated pressure fluctuations on the hull. Cavitation inception forces us to look into micro-scale phenomena, which are not easily observed or measured on propellers. Therefore, fundamental research on hydrofoils in dedicated scale model experiments is imperative on improving general understanding and predictions. The results of cavitation inception tests are often summarized in an inception diagram, i.e. a plot of the inception number against a parameter related to the operating condition of the propeller, e.g. a thrust or torque coefficient. The simplicity of such a representation can be misleading in the sense that the underlying phenomena are complex and sensitive with regard to environmental conditions. A review of phenomena and sensitivities associated with cavitation inception is given by Rood [31]. Performing repeatable cavitation inception tests with acceptable uncertainties is no matter of course. Experience and proper facilities with ample opportunities for visual observation are indispensable for obtaining meaningful results.

2

To better understand the inception behavior of sheet cavitation on propellers, a four year research project on the delay of sheet cavitation inception was started. The David Taylor Model Basin (DTMB) and MARIN were involved in the program which came to be known as ‘AFDEASR’. The program was a success in that sheet cavitation on the demonstrator propeller showed a significantly higher inception speed on full scale than tip vortex cavitation [20]. The next challenge was to delay the inception of tip vortex cavitation as well. The same model basins then started a new research program, called ‘TIPVOR’ (1996-2001), on the origin of tip vortex cavitation and the possibilities for its suppression [21]. With their focus on cavitation inception the two programs produced a wealth of knowledge and practical experience on the subject. The remainder of this section attempts to give a brief account of the state-of-the-art, largely based on the ‘AFDEASR’ and ‘TIPVOR’ programs, as well as some open literature. 2.1.

Repeatability

Cavitation inception is affected by a great number of parameters, such as the character of the wetted flow (laminar vs. turbulent), the occurrence of boundary layer separation and the density of nuclei. Inception often shows a strong sensitivity with regard to minor deviations in geometry and flow conditions (Fig. 1). The fundamental factors that affect cavitation inception are instantaneous pressure, availability of nuclei and time. Residence time, i.e. the time that a nucleus is subjected to a low pressure, plays a role if the bubbles in the flow do not have reached their critical radius beyond which they ‘explode’ into vaporous bubbles. Although the production of vapor is considered to be a near instantaneous process, the preceding growth of a nucleus filled with gas (e.g. air) is subject to inertia effects and occurs over a certain period of time. This period of growth may lead to scale effects, as the residence time on model scale is typically smaller than on full scale. With scale factor λ it scales as 1 λ for Froude scaling and as 1 λ 2 for Reynolds scaling. Considering the sensitivities of inception, the repeatability of the cavitation number at which inception on a propeller blade occurs is by no means guaranteed. Within the ‘TIPVOR’ program, DTMB and MARIN studied the repeatability and checked the effect of the respective laboratories on inception. With sand roughness applied to the tip and leading edge the repeatability was found to be within about 4% in inception number, based on tests performed at an air content of about 5.5 p.p.m. [34]. An unexplained bias error of some 20-25% was also reported for different test sessions. Differences between blades were found to account for 20%. This scatter was the same for a 320 mm propeller manufactured by MARIN and a 400 mm propeller made by another reputable institute. Regarding water treatment, a significant difference between ‘old’ and freshly supplied water in the cavitation tunnel could not be observed. An average change in inception number of 2-4% occurred after ‘old’ water was deaerated and subsequently aerated to the same standard air content of 5.5 p.p.m. Also no effect in inception number between two different observers could be found.

3

Leading Edge geometry Separation bubble with vortex cavitation

Leading Edge geometry

Bubble cavitation without separation

Katz [1984]

Cp Cp Leading Edge geometry Cavitation inception on a roughened surface Cp

Flat plate rectangular to mean flow

Vortex cavitation in a separated shear layer O’Hern [1986]

U

0

Fig. 1. Mechanisms of inception at the leading edge of a foil and in the shear layer of separated flow [16,28]. 7.00 suct. side 5168 (Cgm = 2.2 ppm) pres. side 5168 (Cgm = 2.2 ppm) suct. side 5168 (Cgm = 5.4 ppm) pres. side 5168 (Cgm = 5.4 ppm)

6.00

suct. side 6468 (Cgm = 2.2 ppm) pres. side 6468 (Cgm = 2.2 ppm) suct. side 6468 (Cgm = 5.6 ppm)

5.00

pres. side 6468 (Cgm = 5.6 ppm) suct. side 5168 (p0 = 151.68 kPa)

DTMB-tests

suct. side 5168 (p0 = 113.76 kPa) pres. side 5168 (p0 = 151.68 kPa)

σn

4.00

3.00

2.00

1.00 p0 = ambient pressure in flow Cgm = total gas content

0.00 0.00

0.05

0.10

0.15

0.20

0.25

0.30

0.35

0.40

0.45

KT

Fig. 2. Scatter in repeated experiments of cavitation inception in two facilities with two propeller models representing the same geometry (referred to as nos. 5168 and 6468).

4

As to the difference between MARIN’s cavitation tunnel and DTT, it should be noted that in the DTT cavitation the desinence condition is usually established, whereas in the cavitation tunnel cavitation inception is determined. The procedure followed for the determination of cavitation inception or desinence is also different [35]. These differences result in a cavitation desinence number from the DTT that is about 10% higher on average than the inception number found from the cavitation tunnel in a thrust coefficient, KT , range of 0.3 to 0.4. Concerning MARIN’s cavitation tunnel, with a cross-section of 900 by 900 mm, and DTMB’s open jet cavitation tunnel, measuring 1100 mm in diameter, the difference in inception number for suction side tip vortex cavitation did not exceed 11% in cavitation number, σ n , in the KT range of 0.27 to 0.38. Large differences occurred, however, in pressure side cavitation inception between the two facilities (Fig. 2). These differences may well be caused by the application of tip and leading edge roughness in MARIN’s cavitation tunnel and by the faintness of the cavitating vortex, when only a bubble trace is visible. 2.2.

Case study

After an indication of the uncertainty in the model test results had been obtained, a systematic optimization study was conducted, eventually leading to the redesign of a propeller for a naval supply vessel. This propeller was tested, built and fitted to the ship, after which extensive viewing trials were held in the Mediterranean. Fig. 3 shows a comparison between the predicted cavitation inception results for this newly designed propeller and full scale trial results. The figure compares visual inception data with scale model results from the DTT, extrapolated to full scale. One clear observation is that the predicted operating point with the lowest cavitation inception number for an off-design pitch of approximately -7 deg. corresponds well with the optimum operating point following from the full scale measurements. For a further comparison of the predicted and full scale cavitation inception curves, one should compare inception numbers at equal pitch setting. However, such a comparison is not straightforward as insufficient data are available for a proper interpolation between pitch angles, both for model and full scale. Since the ship is of the single screw type, a variation in operating points could only be obtained through a variation of propeller pitch. Error bars around the full scale observations indicate the estimated uncertainty in the derived torque coefficient, K Q . The difference in pressure side cavitation inception between the best and worst blade was determined to be about 2 degrees in pitch. A corresponding uncertainty of ±1 deg. in pitch has been assumed. The scatter between blades in cavitation inception point on suction side appeared to be smaller than on pressure side. The difference on suction side between the best and worst blade was approximately 1 deg. in pitch. All indicated inception points are based on inception of the ‘median’ blade.

5

10 9

dots : full scale lines : model scale

0 deg

-6.7

8 7

0 deg

σ

6

Full Scale Results

5-12.1 4 3

-11.0

-2.0

-10.0

-3.0

-8.9

2

-3.9 -5.0

-7.9

-5.8

1 0 0.0080

0.0130

Design Condition

0.0180

0.0230 KQ

-6.7 deg face sheet cav.

0.0280

Figures in the graph refer to deviation (in deg.) from design pitch angle

0.0330

0.0380

Fig. 3. Comparison of full scale visual cavitation inception with similar results from model tests in the DTT, the latter extrapolated to full scale values.

The correlation between predicted and observed cavitation inception on pressure side is very satisfactory for the lightest propeller loadings. This holds for both the character of the cavity and the inception numbers. It is noted that in the model tests the first type of cavitation to become visible on the pressure side was leading edge vortex cavitation. During full scale observations, the images were made at a lower resolution, which made it more difficult to discern leading edge vortex and sheet cavitation. Nevertheless, cavitation was observed in the same tip region as in the model tests, whether it was sheet or leading edge vortex cavitation. The correspondence between model scale prediction and full scale measurement becomes less favorable for operating points closer to the design operating point (i.e. the deepest point of the cavitation inception bucket, see Fig. 3). A possible explanation for this fact is that closer to the design condition the blade tip loading becomes lighter. From model scale observations one could clearly observe the low vorticity associated with the bubble trace shed from the leading edge. Therefore, it is likely that hardly any vorticity is shed in the bubble trace that trails the leading edge at pressure side. Still, all model scale inception points on pressure side were attributed to leading edge vortex cavitation, thus indicating that extrapolation from model to full scale should be done according to McCormick’s rule for cavitation inception of vortices [24]. Looking at the increasing deviation between model and full scale inception when getting closer to the design condition, it seems plausible that there is a transition from

6

leading edge vortex cavitation to leading edge sheet cavitation with a corresponding transition in the correction for model scale effects in inception. This suggestion is confirmed by the good correspondence in predicted and measured inception, when it is assumed that there is no vorticity in the leading edge region for the pitch setting closest to the design condition (at a pitch reduction of -7.9 deg.), and the cavitation inception scales as if it were sheet cavitation, hence without any shift in the inception diagram. Observations on suction side show a somewhat different picture. Regarding the character of the incipient cavitation, it was observed that both at model and full scale the first type of cavitation to occur appears to be local vortex cavitation attached to the aft part of the tip [36]. At higher propeller loadings an isolated suction side sheet develops near the tip at the leading edge and occurs almost simultaneously with the local tip vortex at inception [21]. Despite this agreement in character of the cavity, inception of tip vortex cavitation on full scale occurs at significantly higher ship speeds than predicted. Based on the similarity in cavitation behavior between model and full scale, however, it is concluded that scale effects in the local propeller flow are not likely to cause the discrepancy. Scale effects in the wake are more likely to be responsible for the difference. These wake scale effects are expected to become apparent predominantly in suction side cavitation, as this type of cavitation occurs first in the wake peak, where one would expect the strongest scale effects. Pressure side cavitation inception occurs well outside the wake peak, where the local axial flow speed is at a maximum. Due to the aforementioned deviations around the optimum operating condition for cavitation free operation, the predicted inception speed of the vessel appeared to be lower than the actual inception speed determined during the trials. It is concluded that inception condition of pressure side cavitation is predicted very satisfactorily in this case. This may not be expected at first, as it was mentioned earlier that pressure side cavitation inception appears more prone to scale effects and a large scatter between facilities. The conclusion drawn here, however, only holds for pressure side cavitation in the propeller blade tip region, and for a model test procedure where roughness is applied to the leading edge and the tip of the propeller blade. It is also concluded that the art of making adequate full scale predictions is in distinguishing between vortex cavitation and vorticity free cavitation, be it in the form of bubbles or a sheet. This requires very careful observations. On the suction side the viscous scale effects in the wake peak should be corrected for to arrive at a better full scale prediction. This scaling can be most adequately achieved by computing both model and full scale wake, and then using the wake deficit at the observed location of inception to correct the inception speed prediction. This procedure is only valid if the scale effect in the wake is sufficiently small, otherwise, the locus and character of inception may change. 2.3.

Topics for further research

Having demonstrated the satisfactory performance of DTT predictions with full scale observations, a few remaining issues should be mentioned. These relate to the delayed cavitation inception for pressure side sheet cavitation at the inner radii of propeller blades and on rudders. In both cases, the delay may possibly be ascribed to

7

extended laminar flow in the leading edge region. A flow condition which fails to produce micro-vortices that instigate cavitation inception [31]. Delay of cavitation inception is a problem of all experimental facilities, including large cavitation tunnels. To solve this problem, experimental studies have been and are still being conducted for the DTT and cavitation tunnel. Initial studies have been made of the effect of propeller Reynolds number, Rnn , on cavitation inception in the Rnn range of 0.9 - 2.4 million in the cavitation tunnel. Within this range, no clear trend in change of inception condition with Rnn was found. Experimental values appeared to be within a band of 10% of KT . Also a similar study in the DTT with Reynolds numbers between 0.6 and 1.0 million showed no trend in change of inception point, but a larger scatter in KT was found. At the inner radii of the propeller blades at very low local Reynolds numbers, cavitation inception failed to occur. In the cavitation tunnel, however, at the same Reynolds number and positions on the blades, cavitation inception was indeed observed. Because the cavitation tunnel has a turbulence level which is at least ten times higher than that of the DTT, it is hypothesized that at these low local Reynolds numbers, further to leading edge roughness, a sufficient turbulence level is required to invoke cavitation. Other, yet unpublished, tests also suggest that turbulence levels, provided they contain the turbulence spectra that affect inception, indeed do have a significant effect. These findings correspond to results of investigations by Van Der Meulen and Yuan-Pei [25], who investigated the effect of an upstream wire on cavitation inception on a foil. Research in the coming years will be directed towards a better understanding of the scale effects mentioned and attempting to either account for them or find ways to alleviate or annul their effects. The advent of RANS codes for the computation of multi-phase flows is expected to enhance our understanding of the sensitivity of cavitation inception for e.g. local flow separation. However, it is noted that RANS codes are not expected to be capable of computing all micro-scale phenomena that appear to affect local separation and cavitation inception, such as transition from laminar to turbulent boundary layer flow or the onset of cavitation through microvortices in the boundary layer.

3.

Cavitation-induced hull excitation

Pulsating or collapsing cavities caused by varying propeller blade loading radiate pressure fluctuations that exert a vibratory excitation force on the ship’s afterbody. Thus, cavitation can become the primary cause of inboard noise and vibration. Propeller and ship designers have long been trying to improve their designs in order to keep excitation forces below certain limits that would cause noise and vibration problems or violate contract requirements. In the following an overview is given of important topics in this respect, followed by two items of special interest for the prediction of hullpressures, i.e. the scale effect on the effective wake field and broadband excitation.

8

3.1.

Overview

With the introduction of highly skewed propellers at the end of the seventies, a significant reduction in hull-pressure amplitudes at blade passage frequency (BPF) was obtained. For propellers with stringent noise requirements substantial tip unloading was applied, causing a further decrease in pressure amplitudes at BPF, be it at the expense of efficiency. The same can be said of the choice of propeller rotation direction. When the propeller rotates against the tangential flow, as a rule, an increase in efficiency must be weighted against a decrease in comfort level. Reduced hull-pressure amplitudes at BPF are sometimes accompanied by increased amplitudes at higher orders of BPF. These higher harmonic pressure fluctuations are caused, amongst others, by the interaction of the collapsing sheet cavity and the cavitating tip vortex which may lead to oscillations of the cavitating vortex and possibly break-up of the cavitating structure followed by rebounds. The break-up of the cavitating vortex structure is often called ‘bursting’, which may be a somewhat confusing term due to the analogy with non-cavitating vortex bursting which is a different phenomenon. It has been shown by various researchers that the air content of the testing facility has a significant effect on the higher harmonics and that they can be reduced by injecting air into the cavitating tip vortex. The pressure pulses that are radiated by the cavitating propeller are diffracted by the wetted hull. This transfer path may be interrupted by injecting air at the fluidstructure interface. This has proven to be an effective isolation measure [27]. Isolation is especially effective when the risk of vibration is in broadband excitation, rather than in excitation at orders of BPF, since then it becomes rather difficult to prevent structural resonances from occurring. The interest in broadband excitation due to propeller cavitation has increased much recently. For this reason a separate section is devoted to it. It is sometimes argued (e.g. in [39]) that contract specifications based on maximum pressures levels above the propeller at BPF are inadequate when they are influenced too much by non-cavitating contributions due to blade thickness and loading. This is because, in general, the latter effects are not contributing much to the integrated forces [13]. In this regard it is interesting to note that modern passenger ships show hullpressure amplitudes at BPF of only 1-2 kPa at maximum. A significant percentage thereof is caused by the effects of blade thickness and, although usually to a smaller extent, blade loading. An inverse boundary element method for acoustic wave diffraction may be applied to determine the acoustic source strength of the propeller [41]. This source strength may serve as a useful parameter in contract specifications, since it is not influenced by the non-cavitating contributions. Requirements based on the maximum pressure amplitude above the propeller are indeed biased in this respect and should be replaced by restrictions on the integrated vibratory force. However, model experiments and ship trials, performed with a few flush-mounted pressure transducers, only give limited information on the pressure distribution, and sometimes integration to an effective excitation force is not possible. Van Wijngaarden et al. [42] show how to use the source strength concept to obtain such an integral force, even when only sparse hull-pressure data are available.

9

3.2.

Scale effect on effective wake field

The quality of the wake field is essential for a good propeller design. From a vibration point of view improvements are made by paying close attention to the afterbody shape and the orientation of appendages in order to make the propeller inflow more uniform. For the accurate prediction of vibratory hull-pressure forces the correct representation of the effective wake field on model scale is a prerequisite. However, differences between model and full scale wakes are inherent to model scale experiments performed at Reynolds numbers much lower than on full scale. An extensive correlation study, called ‘CoCa’, of hull-pressure fluctuations predicted at model scale and measured at full scale, confirmed the importance of correctly modeling the effective wake [23]. When the wake scale effect is small, e.g. in passenger ships and other twin screw vessels, the prediction of the maximum pressure amplitude at BPF is generally good (see Fig. 4), whereas for some classes of single-screw ships the correlation is unsatisfactory.

1.60

Costa Atlantica, 75% & 50% MCR, full scale vs. model scale normalized amplitudes of hull pressure pulses @ 1st blade rate frequency (fs: full scale; ms: model scale)

1.40 1.20 1.00 0.80

Full Scale Trials @ small pod angle i ti

Model Scale tests 1xOld 2xNew DTT Full Scale Trials

0.79 0.72

0.69

0.69

0.71

0.70

0.72

Model Scale

0.66

0.60

0.53

0.54

0.54

0.59

0.40 0.20 0.00 75% 75% 75% 75% 75% fs fs fs fs fs

74% 75% 75% ms ms ms

50% 50% fs fs

50% 50% ms ms

Fig. 4. Correlation of normalized max. hull-pressures for ‘Costa Atlantica’ @ 75% MCR and 50% MCR [23].

By default, the model test conditions are chosen such that the correct average propeller loading is attained and measured pressure fluctuations are scaled without any correction. As the sector of the propeller disc from where the cavitation originates is around the centerline wake peak, a correct level of the pressure fluctuations can be predicted only if in this sector the wake distribution, and hence the blade loading, are properly represented. When this is not the case, e.g. when the average wake scale effect is small, whilst the scale effect on the wake peak depth is large (as e.g. in container ships), there is a clear tendency to overpredict the maximum pressure amplitude at BPF. On the other hand, if the average wake scale effect is of the same magnitude as the scale effect on the wake peak depth, the correlation between model and full-scale hullpressure amplitudes is much better.

10

As an improvement on the choice of model test conditions, the blade loading in the upper wake peak sector may be adjusted to that in full scale by requiring equal advance ratios in that area. In this manner the propeller becomes less heavily loaded on average, but the loading of the blades passing the wake peak is much better simulated. The above procedure improves the accuracy of predictions of hull-pressure levels. However, the influence of the scale effect on the circumferential gradient of the wake field may still not have been taken into account properly. For instance, the widening of the wake peak may influence the time history of the cavitation volume. Falcão de Campos et al. [9] studied the effect of the wake peak gradient on the cavity volume history in 2D by applying the boundary element method to the cavitating flow over a blade section placed in a gust. To perform a similar study in 3D, using realistic wake fields, is the next challenge. For the time being, another computational approach may be adopted. Holtrop and Kuiper [13] suggest to determine the wake scale effect on the hull-pressure fluctuations by computing the hull-pressure amplitudes caused by the cavity dynamics for both model and full scale wake distributions, and to use the ratio found as a correction factor on the pressure levels measured in the model scale experiment. This method requires that the results of the hull-pressure fluctuation measurements can be divided into a noncavitating and a cavitating part, since only the latter needs to be corrected. Furthermore, the computation of the pressure amplitudes needs to be reliable and should take into account the effects of the various forms of cavitation. At present, the method can only be applied to sheet cavitation. Both procedures, which may be used in conjunction, require knowledge of the scale effect on the wake distribution, in particular the scale effect on the centerline wake peak depth. Some typical data have been collected from a number of RANS computations in a wide range of Reynolds numbers, which may serve to predict the scale effect on the depth of the wake peak. More computational data are needed in this respect. As the number of RANS calculations is growing rapidly improved insight will become available in due course. For now, the following simple prediction rule, proposed by Holtrop, may be used,

λ ⎡Vmin ⎤ ⎡Vmin ⎤ = ⎢ V ⎥ −⎢ V ⎥ ⎣ ⎦ ship ⎣ ⎦ model 50 with Vmin denoting the minimum velocity in the wake peak. 3.3.

Broadband excitation, its appearance

A phenomenon which has been observed in numerous cases is that of so-called broadband pressure fluctuations. It has become of increasing importance, for which reason it is treated here in some detail. Unlike BPF orders which are of a tonal nature, broadband pressure fluctuations show a hump like increase of the pressure levels in the spectrum. As the hump extents over a limited frequency range, typically between 20 and 100 Hz, it should actually be interpreted as a narrowband stochastic process. The

11

ship structure may easily be excited at its natural frequencies within this frequency region. In general, vibration problems due to broadband pressure fluctuations are observed in two different situations, described below. The first situation involves ships with controllable pitch propellers sailing at low speed with significantly reduced pitch [29]. Due to the reduction in pitch, the loading at the leading edge of the propeller will become negative near the leading edge and a strong leading edge vortex may arise on the face of the propeller in combination with a vortex of opposite rotation direction trailing from the tip region. Either one or both vortices may cavitate. Due to their mutual interaction, short waves arise on the cores which ultimately cause the vortices to spiral around each other. If one vortex is significantly weaker than the other, the weakest vortex will break-up, causing the formation of ring vortices around the strongest vortex. If the strength of both vortices is more equal, both will break-up or collapse, thereby causing large hull-pressure pulses. The second situation involves ships with stringent noise and vibration requirements that are running at service speed. In order to meet comfort requirements, the propeller tip is unloaded and the wake field is smoothly varying in tangential direction. One of the first examples of this situation is found on cruise liner ‘QE2’ [8]. For this particular ship the highest broadband levels were occurring aft of the propeller suggesting that the pressure amplitudes are related with a collapse of the cavitating tip vortices. Fig. 5 shows the hull-pressure spectra of ‘QE2’. Even though the broadband character is not exactly reproduced by model tests, it was sufficiently present to make a sound judgment for new propeller designs, which eventually solved the problem for the ship. Such a pragmatic approach is still used for the propeller design of new cruise vessels of the highest comfort class, such as ‘M/S Color Fantasy’ [14].

Fig. 5. Comparison of full scale (left) and model scale (DTT) (right) measured broadband hull-pressures [8].

12

The broadband character of the pressure pulses in the frequency domain is a result of stochastic variations of the cavity dynamics. Due to the periodicity of the blade passing below the hull, the low frequency hull-pressure spectrum is typically dominated by values at harmonics of BPF. If the acoustic source would be perfectly repetitive one would only see values at orders of BPF. In reality, however, significant variability occurs due to geometric variations between the propeller blades, temporal variations in the ship’s wake field and variability of the cavitation dynamics including inception, desinence and rebounds. Due to such stochastic variations of the cavity behavior, ‘leakage’ of power from the tonals to the sidebands will occur. A detailed discussion of the influence of stochastic variations on cavitation noise is presented by Baiter [1]. The influence of random amplitude and phase modulation on a periodic signal is also discussed by Bark [2]. The amplitude modulation causes a continuous spectrum of which the envelope is identical to that of the isolated acoustic phenomenon. The amplitudes of the tonals at orders of BPF are not affected. However, random phase modulation will decrease the tonals and distribute or ‘smear’ the power over a range of frequencies; the smearing becoming most effective at higher frequencies. Depending on the spectral shape of the isolated acoustic phenomenon, the tonals at higher harmonics may even ‘drown’ in the continuous spectrum and this will appear as a hump with a maximum located at a higher harmonic. An example of model scale measured pressure spectra for two different ships cruising at service speed is presented in Fig. 6. Ship ‘A’ has a propeller with a high tip loading and a pronounced wake peak, which results in a strong interaction between the sheet cavity, its re-entrant jet, and the cavitating tip vortex. Due to the interaction a collapse of the cavitating tip vortex is observed. As the collapse is highly repetitive, strong tonals are observed at the higher harmonics with only some smearing at the root. For the fourth BPF and higher orders, however, the BPF values are only slightly above the broadband levels. Ship ‘B’ has a propeller with an unloaded tip and a more smoothly varying wake field. The cavity pattern is then characterized by a cavitating tip vortex, which, when the blade is in the top position, grows and turns into an attached cavity of elliptical shape. Due to the growth and shrinking of this cavity, the cavitating tip vortex becomes distorted and may start to radiate noise. However, as no collapse occurs, the radiated noise is not as high as for ship ‘A’. Furthermore, it appears that there is a strong variability of the cavity dynamics. Analysis of the hull pressure variation from revolution to revolution for frequencies between the second and sixth order of BPF shows that the standard deviation of the pressure at each blade position is an order of magnitude greater than the average value. The variation is well described by a Gaussian probability distribution. In Fig. 7, an example is presented of the ensemble average of the pressure taken over all revolutions and the pressure variation for one distinct revolution. In the individual revolution, a ‘ringing’ phenomenon composed of higher order frequency components can be clearly identified. These events vary significantly in amplitude and phase from blade to blade and from revolution to revolution. As a result, the ensemble average of the pressure signal is mainly determined by the first and second order of BPF and the higher harmonics are hardly visible.

13

160.0

ship A

150.0

ship B

Lp dB Pa^2/Hz

140.0

130.0

120.0

10 dB 110.0

100.0

0

1

2

3

4

5

6

7

BPF harmonic

Pressure [Pa]

Fig. 6. Examples of hull pressure spectra for two different ships, obtained from DTT experiments.

0

90

180 blade position [deg]

270

360

Fig. 7. Typical hull-pressure variation with blade position for ship ‘B’ for a single revolution (gray) and averaged over a number of revolutions (black).

3.4.

Broadband excitation, its possible causes

In some papers describing experimental data for the low frequency radiated noise of cavitating tip vortices it is suggested that a strong local increase in the spectrum is related to distortions of the cavitating vortex core, propagating as so-called Kelvin waves on the core surface [7,32]. The presence of these waves for different modes on cavitating vortex cores and their dispersion relation was first found by Lord Kelvin [33], while acoustic aspects were derived by Morozov [26]. The lowest order mode is the breathing mode which can be represented by a line array of monopoles. The following mode is a bending mode causing a ‘serpentine’ motion of the cavitating core which can be represented by a line array of dipoles. Higher order modes are represented by higher order poles which have very low radiated noise efficiency and are therefore less interesting from a radiated noise point of view. So far, however, in literature only

14

the dispersion relation has been analyzed to derive the frequency components of the various wave modes. The relation for the frequency, f n , of mode n scales as

fn ∝

Vt 1 2 ( p0 − pv ) = r r ρ

where r is the radius of the cavitating core, p0 − pv denotes the difference between free stream and vapor pressure, ρ is the density, and Vt is the tangential velocity component at the core which is related to the free stream pressure through Bernoulli’s equation. It is observed that, with an increasing amount of cavitation due to higher propeller loading, the hump in the spectrum increases in height and its centre moves to lower orders of BPF (Fig. 8), which is consistent with the formulation given above. At some point, however, the broadband levels no longer change, and only an increase in level of the tonals at BPF harmonics is observed.

145.0

propeller thrust

Lp dB Pa^2/Hz

135.0

125.0

115.0

10 dB

105.0

0

1

2

3

4

5

6

7

BPF harmonic

Fig. 8. Example of the influence of propeller thrust on hull-pressure spectrum for ship ‘B’.

Empirical relations between the diameter of the cavitating vortex core, the cavitation number, and the propeller loading for a propeller in the open water condition appear hard to develop due to e.g. the diffusion of non-condensable gas into the vortex core [7,36]. Instead, a simple Rankine vortex model can be used, and the relation for the mode frequency is then written as,

fn ∝

σn KT

15

fb

in which f b denotes the BPF. The trend is confirmed through the analysis of a limited number of model tests. It is noted that a similar relation for the centre frequency of the hump is derived by Raestad [30] in his empirical model for tip vortex cavitation noise using full scale data, be it for the fact that his relation is based on a bubble collapse frequency. Kelvin waves have been shown to occur on the viscous core of non-cavitating vortices and may grow in the presence of strain as shown e.g. by Leweke and Williamson [22]. For cavitating vortices, the relation between standing Kelvin waves and observed distortions of the cavitating core has been investigated by Keller and Escudier [17] e.g. The general trends were well captured by the theory, but quantitatively only a moderate correlation with experiment was obtained. Even though the above suggests the presence of Kelvin waves on cavitating tip vortices of ship propellers, their importance for hull-pressure pulses should be further investigated. In this respect note the well known fact that subsonic waves on an infinite cylinder do not radiate noise at low frequencies [15]. The cavitating tip vortex distortions may have several causes, of which some, such as the interaction between two counter-rotating vortices, have already been mentioned. More examples have been presented by Kuiper [18] and Van Wijngaarden et al. [40]. Here, only the cavity dynamics causing the spectra in Fig. 6 are briefly discussed. The attached tip vortex cavity may grow on the blade forming a cavity, which has an elliptical cross-section. While shrinking, the elliptical shape is shed downstream, causing a ribbon shape appearance of the tip vortex cavitation. An example is given in Fig. 9. This type of cavitation may be present on ships with weak tonals and low frequency broadband hull-pressures such as ship ‘B’ of Fig. 6.

Fig. 9. Example of tip vortex cavitation in a smoothly varying wake field.

Fig. 10. Tip vortex cavity which is about to burst due to the interaction with the sheet re-entrant jet.

16

The most common observed distortion of the tip vortex, however, is caused by the interaction with a sheet cavity on the blade, which can take many different shapes. A cavitating vortex structure may emerge from the sheet or the re-entrant jet may form a vortex structure, which interacts with the tip vortex. Another possibility is that the reentrant jet destroys the tip vortex when large velocity components are directed towards the tip. An example where the sheet collapses towards the tip, thereby destroying the tip vortex is shown in Fig. 10. This is typically the behavior that is observed on ships with high tonal amplitudes at the higher harmonics, such as ship ‘A’ of Fig. 6. Current research at MARIN is directed towards further understanding of the stochastic variations of hull-pressure fluctuations on model and full scale in order to improve the interpretation of pressure spectra using Fourier analysis. Additionally, further theoretical and experimental study of low frequency radiated noise from cavitating vortices is needed. This should involve cavitating core deformations, the relation to Kelvin waves as well as the collapse and rebounds of cavitating structures.

4.

Cavitation erosion

Of all adverse effects of cavitation, erosion is perhaps the most complicated to assess. The difficulty is caused by the energy cascade that eventually leads to erosion. This cascade of energy conversions, from macro-scale cavity through micro-scale bubble collapse, and pressure wave generation to material damage, is described in a generic model by Fortes-Patella et al. [11,12]. The present understanding of the hydrodynamic part of this cascade of energy is concisely described by Bark et al. [3]. The challenge in research on cavitation erosion is first to improve the quality of erosion risk assessments. A second objective is to predict the erosion rate, thereby including the last step of the energy cascade. However, the first aim is not yet met to full satisfaction. This is illustrated by a statement of Bark et al. [4]: “Some full scale results from EROCAV indicate however that quite badly looking pressure side cavitation is in fact not always very erosive”, a statement which is also confirmed by Ligtelijn et al. [23]. This statement reveals the weakness of current experimental practice, a practice based on careful visual observations, linking the observer’s interpretation of the dynamics of the cavity on model scale to a prediction of the risk of erosion. The contribution of the EU-sponsored ‘EROCAV’ project consists of a set of basic guidelines and nomenclature to assist the observer in analyzing the dynamics of the cavity pictures, and then translating them into an assessment of the risk of erosion. These guidelines are presented in the ‘EROCAV’ handbook [3], which is an important contribution to the enhancement of erosion risk assessment. The reliability of this method, however, still depends strongly on the experience of the observer and of the quality of the pictures. Bark et al. [4] present a brief review of erosion prediction techniques that were evaluated within the ‘EROCAV’ project. We mention the ‘High Speed Video Method’, the ‘Paint Test Method’, the ‘Impact Method’, the ‘Soft Metal Method’, and the ‘Acoustic Emission Method’. The ‘High Speed Video Method’ and the ‘Paint Test Method’ were reported to be the most reliable methods at present. However, it should

17

be recognized that these methods are still of a qualitative nature and the quality of prediction still depends on the experience of the observer. As the ‘High Speed Video Method’ is concerned, Tukker and Kuiper [37] concluded that high speed video is the preferred way of making observations. It was argued that a conventional video camera, applying a time lapse principle to collect images, may yield misleading results. This is caused by the fact that in erosion risk assessment cavity dynamics is important, something that is not captured by time lapse methods. Until recently, the latter methods were used at low frame rates (25 or 30 Hz), sometimes supplemented by high resolution stills from photo cameras. It is concluded that high speed video significantly improves the visual judgment of the risk of cavitation erosion. The ‘Impact Method’ and the ‘Acoustic Emission Method’ are of a more quantitative nature. In the ‘Impact Method’, high frequency pressure transducers are used to determine both the fluid and the structure borne noise up to about 200 kHz. Exploratory experiments with this technique are reported by Foeth and Kuiper [10]. The authors applied the technique to single bubble collapses near a surface as well as to sheet cavitation. The experiments with a single bubble collapse were made in a dedicated set up, whilst the experiments with sheet cavitation were made on a NACA 0018-45 foil in MARIN’s high speed cavitation tunnel (40 by 80 mm measuring section). Observations on the foil with the sheet cavity indicate that the collapse of the main cavity structure migrates from the foil surface at decreasing cavitation numbers, which is expressed in a decrease in the ratio of structure versus fluid borne noise. Further studies are necessary to show whether this ratio can indeed be used as a signal to indicate the risk of erosion. The ‘Acoustic Emission Method’ is a technique largely similar to the ‘Impact Method’. The technique is further developed for full scale application at Lloyd’s Register [6]. The acoustic signal associated with impacting cavitation is measured by means of high frequency pressure transducers. The technique not only aims at indicating erosion, but also in localizing it. However, further development of this technique is needed for it to be used on ships in service. The ‘Soft Metal Method’ uses a soft metal insert in the body that is to be tested, with the aim to increase cavitation erosion on model scale. As such, this method has much in common with the ‘Paint Test Method’. A difference is that the ‘Soft Metal Method’ might have the capability of predicting damage rates. However, more research is necessary for the method to be used in a quantitative way. For the latter purpose this method is not yet recommended in practice. As it must be realized that a reliable prediction of the erosion rate is not yet within reach [5] our aim must be limited to the prediction of the risk of cavitation erosion. Although an improvement in the assessment of the risk of erosion has been made through the use of high speed video cameras and their improved resolution, a more quantitative assessment is still to be aimed for. Therefore, MARIN has chosen to further develop the ‘Impact Method’, which promises to produce quantitative data that can be used in conjunction with case independent threshold levels. Thus an independent judgment can be made in addition to the one made on the basis of high speed video pictures.

18

The above-mentioned techniques only work adequately if the character of the cavitation observed and measured on model scale is similar to that on full scale. Ideally, this similarity constraint should hold for both the macroscopic and microscopic structure of the cavity. More studies on the similarity of these structures for model and full scale are needed to get a better understanding of those cases where cavitation erosion occurred, but was not foreseen, not even from high quality experiments. The ‘EROCAV’ project has already greatly contributed to this need, but it is widely felt among the members of this project that additional high quality correlation data are needed.

5.

Conclusions

The present paper gives a review of cavitation research on ship propellers performed by MARIN in the recent past. Investigations regarding cavitation inception, cavitation-induced vibratory hull-forces and cavitation erosion are treated. Below, a summary is given of achievements made in these areas and outstanding issues that challenge future research. Regarding cavitation inception, it is concluded that a rise in the inception speeds for naval vessels has been obtained, as well as a better knowledge regarding the fundamental mechanisms controlling the inception process. Thus, it becomes possible to find the right balance between the inception of e.g. sheet and tip vortex cavitation in propeller design. There are, however, still problems remaining with reproducing cavitation inception at model scale, predominantly at the inner radii on the pressure side of propeller blades, on shaft struts and rudders. In the cases mentioned, cavitation inception in model experiments is delayed compared to full scale observations. Hypotheses for these phenomena have been suggested and work is ongoing to verify them. Cavitation-induced hull-pressure fluctuations at blade passing frequency have been significantly reduced through improved propeller and afterbody design. Passenger vessels are now reaching amplitudes as low as 1 kPa, and rarely show values much higher than 2 kPa. The large cavitation facilities are capable of predicting these levels quite well. However, when the wake scale effect becomes strong in the sector of the propeller disc where cavitation occurs, as with the large container vessels of today, the accuracy deteriorates and there is a tendency to overpredict. Methods to counteract this tendency are proposed. For propellers showing tip vortex cavitation, the hull-pressure spectrum is often characterized by a high broadband region. In a qualitative sense, this broadband region can be captured by model tests, which enables comparative propeller analysis. Quantitatively, however, differences occur which emphasize the need for further research on the variability of the pressure signals due to tip vortex cavitation and for mathematical models for the radiated noise of cavitating vortices. More reliable predictions of the risk of cavitation erosion have recently been obtained using observations made by means of high speed video. Furthermore, a better insight in cavitation erosion has been gained, but there is still an undesirable uncertainty in the prediction of the risk of erosion. A quantitative, less expert dependent procedure

19

is needed in the assessment of the risk of cavitation erosion on propellers and rudders. To this end, MARIN is further developing the so-called ‘Impact Method’ for use in their Depressurized Towing Tank. Notwithstanding the improved understanding of the fundamentals of cavitation on propellers and rudders, it must be realized that some scale and laboratory effects are difficult, if not impossible, to overcome. Reliable methods to treat the scale effects are therefore needed, which may not come from simple empirical rules, but may need dedicated solvers for the solution of the complicated flows over propulsors and rudders. MARIN is involved in the development of an analysis tool based on potential flow for cavitating propellers within the Cooperative Research Ships (CRS) [38]. Furthermore, MARIN is developing a propeller dedicated RANS solver for cavitating flows in cooperation with HSVA (Hamburg Model Basin). This code aims to accurately predict flow phenomena, such as local separation and vorticity generation. A better insight in and prediction of these flow phenomena is expected to raise the reliability of cavitation related predictions further.

Acknowledgements This review could only be made through the participation of MARIN in a number of sponsored Joint Research Projects. In particular, the following projects have contributed to this paper. The four year TIPVOR project, supported by the Office of Naval Research, the US Navy and the Royal Netherlands Navy. For both projects, model tests were carried out at NSWCCD (a.k.a. the David Taylor Model Basin), Carderock Division, W.-Bethesda, MD, USA and at MARIN, Wageningen, The Netherlands. The three year project CoCa (Correlation of Cavitation) was financially sponsored by the Dutch government and Wärtsilä Propulsion Netherlands. The four year EROCAV project was funded by the European Community under the ‘Competitive and Sustainable Growth’ Program (1998-2002), Project number GRD1-2000-25089 (see www.erocav.de). The authors are much indebted to all those people who have contributed to our knowledge in the areas treated in this paper. Especially the contributions by Mr. Jan Holtrop, head of the propeller design group at MARIN, are much appreciated.

References [1] [2] [3] [4] [5]

Baiter, H-J., “On different notions of cavitation noise and what they imply”, 2nd int. symp. on cavitation and multi-phase flow, ASME, Anaheim, CA, USA, 1986. Bark, G., “On the mechanisms of propeller cavitation noise”, Ph.D. thesis, Chalmers University of Technology, Göteborg, Sweden, 1988. Bark, G., Berchiche, N., Grekula, M., “Application of principles of observation and analysis of eroding cavitation – The EROCAV observation handbook”, EROCAV Report, Dept. of Naval Architecture, Chalmers University of Technology, Göteborg, Sweden, 2004. Bark, G., Friesch, J., Kuiper, G., Ligtelijn, J.T., “Cavitation Erosion on Ship propellers and Rudders”, PRADS 2004, Lübeck, Germany, Aug. 2004. Billet, M.L. et al., Report of the Specialist Committee on Cavitation Erosion on Propellers and Appendages on High Powered/High Speed Ships, Proc. of the 24th ITTC, Vol. 2, Edinburgh, Scotland, Sept. 2005.

20

[6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30]

Boorsma, A., Presentation to the EROCAV workshop on the Acoustic Emission Method for the detection and location of cavitation erosion, EROCAV workshop, Val de Reuil, France, May 2004. Briançon-Marjollet, L., Merle, L., “Inception, development and noise of a tip vortex cavitation”, 21st Symp. on Naval Hydrodynamics, Trondheim, Norway, 1996. Brubakk, E., Smogeli, H., “QE2 from turbine to diesel - consequences for noise and vibration”, IMAS 88, The design and development of passenger ships, pp. 169-176, 1988. Falcão de Campos, J.A.C., Vaz, G., Bosschers, J., “Influence of the reduced frequency on the 2-D partial sheet cavitation behavior for a foil in a gust”, 6th int. symp. on cavitation CAV2006, Wageningen, NL, Sept. 2006. Foeth, E.J., Kuiper, G., “Exploratory experiments to determine flow and structure borne noise of erosive cavity implosions”, ASME Heat Transfer/Fluid Eng. Summer Conference, Charlotte, NC, USA, July 2004. Fortes-Patella, R., Reboud, J.L., Briançon-Marjollet, L., “A phenomenological and numerical model for scaling the flow aggressiveness in cavitation erosion”, EROCAV Workshop, Val de Reuil, France, May 2004. Fortes-Patella, R., Challier, G., Reboud, J.L., “Cavitation erosion mechanism: Numerical simulations of the interaction between pressure waves and solid boundaries”, CAV2001, Pasadena, CA, USA, June 2001. Holtrop, J., Kuiper, G., “Cavitation and vibration research in the new depressurized towing tank”, Sprechtag Kavitation, Schiffbautechnische Gesellschaft, Hamburg, Germany, 2003. Hämäläinen, R., Lönnberg, B., Arén, P., Pettersson, G., “Highest Comfort Class design for ‘M/S Color Fantasy’, the world’s largest ever cruise-liner with car-deck”, 1st int. ship noise & vibration conf., London, UK, June 2005. Junger, M.C., Feit, D., “Sound, structures and their interaction”, 2nd edition, MIT Press, Cambridge, MA, USA, 1986. Katz, J., “Cavitation phenomena within regions of flow separation”, Journal of Fluid Mechanics, Vol. 140, 1984. Keller, J.J., Escudier, M.P., “Theory and observations of waves on hollow-core vortices”, Journal of Fluid Mechanics, Vol. 99, part 3, pp. 495-511, 1980. Kuiper, G., “New developments around sheet and tip vortex cavitation on ship’s propellers”, Fourth Int. Symp. on Cavitation, Pasadena, CA, USA, 2001. Kuiper, G., “Cavitation testing of marine propellers in the NSMB Depressurized Towing Tank”, Inst. Mech. Eng., Conf. on Cavitation, Edinburgh, Scotland, 1974. Kuiper, G., Jessup, S.D., “A propeller design method for unsteady conditions”, Trans. SNAME, 1993. Kuiper, G., Terwisga, T.J.C. van, Zondervan, G.J., Jessup, S.D., Krikke, E.M., “Cavitation inception tests on a systematic series of two-bladed propellers”, 26th Symp. on Naval Hydrodynamics, Rome, Italy, Sept. 2006. Leweke, T., Williamson, C.H.K., “Cooperative elliptic instability of a vortex pair”, Journal of Fluid Mechanics, Vol. 360, pp. 85-119. Ligtelijn, J.T., Wijngaarden, H.C.J. van, Moulijn, J.C., Verkuyl, J.B., “Correlation of cavitation: comparison of full scale data with results of model tests and computations”, SNAME Annual Meeting, Washington D.C., USA, 2004. McCormick, J., “On cavitation produced by a vortex trailing from a lifting surface”, ASME Journal of Basic Engineering, 1962. Meulen, J.H.J. van der, Yuan-Pei, Y., “Cavitation inception scaling by roughness and nuclei generation”, 14th Symp. on Naval Hydrodynamics, Ann Arbor, MI, USA, 1982. Morozov, V.P., “Theoretical analysis of the acoustic emission from cavitating line vortices”, Sov. Phys. Acoust., Vol. 19, No. 5, March-April 1974, pp. 468-471. Ødegaard, J., “Practical experience with reduction of propeller-induced noise and vibration”, Proc. of the 2nd int. ship noise & vibration conf. 2006, London, UK, June 2006. O’Hern, T., Acosta, A.J., “Cavitation observations in a turbulent shear flow”, 21st meeting of the American Towing Tank Conf., Washington D.C., USA, 1986. Okamura, N., Fujino, R., Tanaka, T., “An experimental investigation of the mechanism and the pressure of counter-rotating vortices on a CPP at the off-design condition” 20th Symp. on Naval Hydrodynamics, Santa Barbara, CA, USA, pp. 341- 354, 1994. Raestad, A-E., “Tip vortex index – an engineering approach to propeller noise prediction”, The Naval Architect, July/August 1996.

21

[31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42]

Rood, E.P., “Review – Mechanisms of Cavitation Inception”, J. of Fluids Eng., Vol. 113, June 1991. Sponagle, N.C., “Noise from tip vortex and bubble cavitation”, DREA Tech. Memorandum 90/202. Thomson, W. (Lord Kelvin), “Vibrations of a columnar vortex”, Proceedings of the Royal Society of Edinburgh, Scotland, March 1, 1880; Phil. Mag. X. 1880, pp. 152-165. Terwisga, T.J.C. van, Borg, E., “Tip vortex cavitation - experimental study on DTMB propeller model 5168 and MARIN propeller model 6468”, MARIN Report no. 13509-4-CT, January 1997. Terwisga, T.J.C. van, Borg, E., “Further cavitation tests on MARIN propeller 6468 - test series 2”, MARIN Report no. 13509-10-CT (figure 11), January 1997. Terwisga, T.J.C. van, Kuiper, G., Rijsbergen, M.X. van, “On experimental techniques for the determination of tip vortex cavitation on ship propellers”, Proc. of 1999 ASME/JSME Fluids Eng. Division Summer Meeting, San Francisco, USA, July 1999. Tukker, J., Kuiper, G., “High speed video observations and erosive cavitation”, PRADS 2004, Lübeck, Germany, Aug. 2004. Vaz, G., Bosschers, J., ”Modelling of three dimensional sheet cavitation on marine propellers using a boundary element method”, 6th int. symp. on cavitation CAV2006, Wageningen, NL, Sept. 2006. Wijngaarden, H.C.J. van, “Recent developments in predicting propeller-induced hull pressure pulses”, Proc. of the 1st int. ship noise & vibration conf. 2005, London, UK, June 2005. Wijngaarden, E. van, Bosschers, J., Kuiper, G., “Aspects of the cavitating propeller tip vortex as a source of inboard noise and vibration”, FEDSM2005-77271, ASME Fluids Engineering Division Summer Meeting and Exhibition, Houston, TX, USA, 2005. Wijngaarden, H.C.J. van, “Determination of propeller source strength from hull-pressure measurements”, Proc. of the 2nd int. ship noise & vibration conf. 2006, London, UK, June 2006. Wijngaarden, H.C.J. van, Brouwer, J., Zijlstra, M., “A spatial extrapolation technique for sparse measurements of propeller-induced hull-pressures using an inverse acoustic boundary element method”, 6th int. symp. on cavitation CAV2006, Wageningen, NL, Sept. 2006.

22