• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Novel scaling law for estimating propeller tip vortex cavitation noise from model experiment *

    2017-03-14 07:06:26JisooParkWoojaeSeong

    Jisoo Park, Woojae Seong

    Department of Naval Architecture and Ocean Engineering, Seoul National University, Seoul, Korea, E-mail: online1616@snu.ac.kr

    Introduction

    Since tip vortex cavitation (TVC), usually the earliest type of cavitation to appear on a marine propeller, is related to cavitation inception speed (CIS)and underwater radiated noise (URN), the physical understanding of TVC is critically important in designing a propeller. In terms of URN, TVC noise is considered to be the most significant source of marine propeller noise because it occurs the earliest among other types of cavitation and significantly influences the noise spectrum over a broad frequency range[1].

    Because of the complex flow and noise field around a marine propeller, a theoretical approach to TVC noise estimation is practically unrealizable. Thus,estimation of prototype noise level is realized through extrapolation of the model-scale measurement made in a cavitation tunnel. For this purpose, some scaling laws on cavitation noise have been suggested, which have been applied in many studies.

    Levkovskii[2]initially proposed the cavitation noise scaling law, assuming cavitation as a single bubble. The derivation of the cavitation noise scaling law started from the fact that the ratio of radiated acoustic power to potential power of a single bubble is constant, which was based on experimental results.However, this scaling law could not be applied to the cavitation noise generated by a 3-D lifting surface since propeller parameters and boundary layer effect are not taken into account. For the cavitation attached to a propeller blade surface, such as sheet cavitation,the scaling law for propeller cavitation noise is well established by the ITTC noise estimation rule[3]. This scaling law has been applied to predict propeller cavitation noise from model tests in two studies[4,5]. In these studies, it was shown that the ITTC noise estimation rule had good agreement with the measured prototype cavitation noise level in the case of sheet cavitation developed on the blade surface.

    It is known that TVC development is significantly affected by the boundary layer and propeller loading[6]. However, since the ITTC noise estimation rule does not include these effects properly, it would be inadequate to apply the ITTC noise estimation rule to cavitation occurring continuously at a distance from the propeller such as TVC. Latorre[7]measured the fully developed TVC noise of scaled hydrofoils and applied the cavitation noise scaling method of Bojorheden and Astrom[8]and the Levkovskii method[2]. Oshima[9]modified the ITTC noise estimation rule to incorporate the boundary layer effect using the Reynolds number, which significantly affects TVC occurrence. Prototype power spectral density (PSD)estimated with the modified ITTC noise estimation rule was compared with measured prototype PSD for verification. In this study, the scaling exponent determined from model and prototype experiments was 0.150, quite different from other experimental results(summarized in Shen et al.[10]). More importantly, the estimated prototype PSD did not match well with the measured data in the high frequency range above 6 kHz-7 kHz. From these analyses, it is ambiguous as to whether the model or prototype experiments are incorrect or the modified ITTC noise estimation rule gives incorrect estimation result. Since model and prototype experiments were conducted based on the visual inspection of TVC using an observation window (the method considered accurate), it could be considered that the experiments were performed correctly. Thus, it would be reasonable to re-examine the modified ITTC noise estimation rule that reflects only the boundary layer effect and suggest that such a treatment may be inadequate and incomplete.

    To overcome this limitation of considering the boundary layer effect in applying the ITTC noise estimation rule for TVC noise, the scaling law for TVC noise is newly derived using the Rayleigh-Plesset equation, the Rankine vortex model, the lifting surface theory, and extra assumptions based on the physics of TVC, which will be detailed in Section 2.This will provide the broadband noise scaling law considering the hydro-acoustic basis of TVC. In Section 3, model-scale data are extrapolated with the modified ITTC noise estimation rule as well as the novel scaling law deduced in Section 2, and these two sets of extrapolated results are compared with the prototype measurements for verification. The causes of errors are also analyzed. Finally, the conclusions are provided in Section 4.

    1. Derivation of the scaling law for TVC noise

    TVC is a small bubble flume generated in the vortex core center and noise is emitted from the bubble growth, elongation, and collapse. Since it is impossible to consider the noise of each bubble when deriving the scaling law, derivation of the scaling law for TVC starts with consideration of the total power of the bubble flume as a whole. In addition, TVC results from nuclei that enter into a vortex trailing from the tip caused by crossover of fluid from the high pressure side to the low pressure side of the propeller. Therefore, to study this physical phenomenon it is necessary to bring all factors into play that underlying the physical basis of TVC, such as the bubble dynamics(Rayleigh-Plesset equation), vortex model (Rankine vortex model), the lifting surface theory, and the number of bubbles generated per unit time when deriving the scaling law for TVC noise. Since the novel scaling law derived in this study is generally based on Rayleigh-Plesset equation and the formulation for elongation of TVC is the same as that for initial radial growth of TVC (which can be identified in Arndt[11]), it can be normally applied to TVC noise estimation regardless of TVC development stage except for very long TVC.

    1.1 Rayleigh-Plesset equation

    The Rayleigh-Plesset equation, commonly used to describe the growth and collapse history of a bubble,is given for the bubble radius R in relation to time[12]as follows

    where ρ denotes the density of fresh water, R radial displacement of bubble, R˙ radial velocity of bubble, R˙ radial acceleration of bubble,vp vapor pressure, p∞ambient pressure at infinity,0gp partial pressure of gas at the initial state, γ the ratio of the specific heat at constant pressure to that at constant volume,0R equilibrium radius, A surface tension force of bubble wall, and μ dynamic viscosity. As secondary terms such as gas pressure,surface tension, and viscosity can be ignored, the Rayleigh-Plesset equation can be simplified as follows

    The above equation is integrated from initial time t = 0 to arbitrary time t and bubble wall velocity is expressed as

    In the growth stage of a bubble, R ( 0)/R( t) is assumed to be zero except at the initial instant (t = 0).Thus, the bubble growth velocity becomes constant after the initial stage, that is

    Denoting the time for the bubble to grow to its maximum size as τG, the maximum bubble radius,Rmaxbecomes

    1.2 Rankine vortex model

    TVC generally forms at the center of the vortex flow generated by the tip of a 3-D lifting surface. One of the simplest models that explains this vortex flow structure generated at the tip of a 3-D lifting surface is the Rankine vortex model[13]. In this model, the vortex flow is separated into two domains: the inside and outside of the core.

    1.2.1 Core inner region (r ≤ac)

    In the inner core region, vorticity is assumed to be constant and the viscous effect becomes dominant.The flow inside the core (r <ac) is considered as the rotational motion of a solid body. At distance r from the core’s center axis, tangential velocity, Vtis expressed as follows

    where Γ denotes circulation andca vortex core radius. Since the inner region of the vortex flow is rotational and viscous, Bernoulli’s equation cannot be applied. Euler’s equation for steady flow becomes

    Pressure co-relation between the core center (r = 0,subscript 1) and core radius (r = ac, subscript 2) is derived from Euler’s equation as

    Here, tangential velocity at the core center is zero(V1=0) since on the core axis the flow is considered irrotational.

    1.2.2 Core outer region (r >ac)

    In the outer core region, the flow is irrotational and inviscid. In this potential region (r >ac),tangential velocity is found from the circulation conservation (r Vt=constant) as follows

    The pressure relation between the core radius (r = ac,subscript 2) and infinity (subscript ∞) is derived from Bernoulli’s equation for irrotational flow as follows

    Here,tV at infinity is zero ( =0)V∞since r is very large compared withca. From Eq.(8) and Eq.(10), the minimum pressure at the vortex core center is

    From this relation, it is evident that the minimum pressure is formed at the vortex core center ( =0R ,subscript 1), which makes it plausible that the vortex cavitation is developed at the vortex core center.

    1.3 Rearrangement of variables to incorporate prope- ller parameters in the scaling law

    For the bubble collapse situation, (0)R in Eq.(3)is simply replaced withmaxR to produce

    The reciprocal of Eq.(12) gives the radius-time relationship, which is integrated over the interval Rmax≤R≤0 giving the time for complete collapse,τCas follows[12]

    where Δp denotes difference between vortex core center pressure and vapor pressure. The frequency of noise radiated from a single bubble collapse, f is inversely proportional to the bubble collapse time( f α1/τc). Thus, the frequency scaling law becomes Subscripts s and m denote prototype and model,respectively. For radiated sound pressure ( p ′), the assumption of a constant ratio of radiated sound power ( Prad=4 π r2( p ′2)/(ρ0c0))to bubble potential power ( Ppot=(4 π Rmax

    3( p ′2)/3)Δ p N0) is adopted[2,8],that is

    where0c denotes sound speed. Thus, the above equation can be rearranged into a scaling form for radiated sound pressure as

    The scaling laws, Eqs.(14) and (16), reveal that the ratio of maximum bubble radius ( Rmax,sRmax,m)and the ratio of the number of cavitation bubbles gene- rated per unit time ( N0,s/N0,m) need to be found in order to predict the radiated noise of the prototype from the model experiment.

    1.4 Ratio of maximum bubble radius in the model and prototype

    When applying the Rayleigh-Plesset equation(Eq.(1)) to real life situations, the driving pressure p∞is regarded as ambient pressure at the center of the bubble rather than the pressure at infinity[13]. In the case of TVC, the pressure at the bubble center is the minimum pressure in the vortex ( p∞= pmin= p1).Thus, the ratio of maximum bubble radius in the model and the prototype can be rearranged from Eq.(5)as follows

    Here, the term ( p1,s- pv) /( p1,m- pv) can be expressed in terms of propeller and flow parameters using the Rankine vortex model. The ratio of difference between vortex core center pressure and vapor pressure in the model and the prototype can be arranged as follows using Eq.(11).

    In the case of TVC, determination of noise measurement conditions is based on the relationship between cavitation numbers and the Reynolds numbers ( σ /σ =(R e/R e)k) to keep the geometri

    n,s n,m s m cal similarity of the TVC pattern in the model and the prototype. Thus, the ratio of difference between the static pressure in the cavitation tunnel and the vapor pressure in the model and the prototype can be arranged as follows:

    where σndenotes cavitation number defined in terms of propeller rotational speed and V0inflow speed. Since the circulation scales as Γs/Γm=(cs/cm) (V0,s/V0,m) = λ (V0,s/V0,m), where c denotes chord length, and the vortex core radius scales as a /a = δ / δ = (R e/R e)0.5k/λ[10], where δ

    c,mc,smssm

    denotes boundary layer thickness and λ scale ratio,the ratio of dynamic pressure of the tip vortex in the model and that in the prototype can be arranged as follows

    Substituting Eq.(19) and Eq.(20) into Eq.(18)yields

    Next, to obtain an expression for the ratio of growth time in the model and the prototype τG,s/τG,m,the position vector of TVC, p is approximated as follows

    where D denotes propeller diameter, n propeller rotational speed and U ship speed. The total length of TVC progression during growth time can be found from

    Fig.1 (Color online) Experimental setup for noise measurement and observation for the same TVC pattern of the model and the prototype. The position of a bubble reaching the maximum radius is at the same azimuth ()θ as from top dead center (TDC) in the model and the prototype

    From Fig.1, we know that the position of a bubble reaching the maximum radius is at the same azimuth ()θ as it is from top dead center (TDC) in the model and in the prototype. Hence, the ratio of length of TVC progression during growth time in the model and the prototype,GL is as follows

    According to Eq.(23) and Eq.(24), the ratio of growth time in the model and the prototype is as follows Substituting Eq.(21) and Eq.(25) into Eq.(17), the ratio of maximum bubble radius in the model and prototype is arranged as

    This is different from the existing assumption that maximum bubble radius only depends on the propeller diameter ( Rmax~D), whic h is simply dete rmine d by dimensionalanalysis.Thisassumptionhasbeencommonly used in applying cavitation noise scaling laws[2,3,9,14]. However, in the case of TVC, the maximum bubble radius depends not only on propeller diameter, but also on the Reynolds number, as indicated in the newly developed Eq.(26). An increase of the Reynolds number indicates a decrease of the boundary layer thickness, which reduces the vortex core radius. When the vortex core radius is reduced,the minimum pressure in the vortex core center is lowered, and the bubble reaches a larger maximum radius. Hence, the proportional dependence of the Reynolds number on the maximum bubble radius is reasonable in reflection of the physics of TVC.

    1.5 Ratio of number of cavitation bubbles generated per unit time in the model and prototype

    The number of cavitation bubbles generated per unit time0( )N is difficult to assess analytically or numerically. There are not many experimental results on0N because scaling effects such as the influence of Reynolds number and roughness of the blade surface have a large impact on the onset of cavities[8].

    Until now, it has been assumed in deriving the cavitation noise scaling law that0N is only proportional to the number of propeller revolutions. In the case of TVC, a cavity is generated around TDC where the angle of attack is high and this cavity progresses downward from the propeller tip, making a successive spiral cavity line. It can be assumed that the number of nuclei flowing into the vortex line is proportional to the number of nuclei dissolved in the liquid per unit volume. Furthermore, similitude of the spectra of nuclei[13]for inflow speed and propeller diameter should be applied to assume the same nuclei distribution in the tip vortex line of the two geometrically scaled experiments. Inflow speed is defined as Advance ratio, J is defined as

    From Eq.(27), the ratio of inflow velocity in the model and the prototype is rearranged as

    Here, propeller rotational speed is higher in the model than the prototype by approximately the magnitude of scale ratio (λ= Ds/Dm)to make the Reynolds number of the model as close to the Reynolds number of the prototype as possible. Furthermore, the advance ratio of the model and that of the prototype is assumed to be equal (Js= Jm) to keep dynamic similarity between the model and the prototype. From these aspects, it can be judged that the ratio of inflow velocity in the model to that of the prototype is almost 1 and thus adjustment of nuclei spectrum for propeller diameter solely needs to be applied. The cumulative spectra of nuclei are shifted in the vertical direction with a ratio (Ds/Dm)2-kin the case of TVC[15],where kdenotes scaling exponent. This adjustment is made to equalize the number of bubbles passing through the unit cross-sectional area in the vortex line of the model to that of prototype propeller. Here,nuclei sizes are not considered. Hence, the number of nuclei flowing into the vortex line ( n0( V0, D)) is proportional to the number of nuclei dissolved in the liquid per unit volume ( n0) and is a function of inflow velocity and propeller diameter (f( V0, D)) which corresponds to the adjustment of nuclei spectrum.

    If the prototype experiment were performed in

    lineshowninFig.2 (borrowed from),which corresponds to the cumulative nuclei spectrum of the prototype experiment performed in fresh water.The cumulative nuclei spectrum of prototype is the solid red line in Fig.2, which equates to the reference line multiplied by 10 since the prototype experiment is performed in sea water and the number of nuclei in

    sea water is approximately 10 timesthe number of

    nuclei in fresh water for most.There are some

    Fig.2 (Co lor onl ine) Adjustment of the cumulative spectrum ofnuclei with changi ng exp erimental c ond itions . If the pro totypeexperimentw ere pe rformed inf resh water,therefer enc e line wouldbe theci rcledblue lineco rrespondin gtothe cum ulativenu cle ispec trumof the prototypeex perime ntperformedinfreshwater.Thecumul ative nuclei sp ectru m of the proto type is the s quared redline,whichisthereferen ce linemu ltip lied by 10.The adjuste d cum ulative nucle isp ectrumof the mod el experimentisthedotted bluelineresultingfromthereference line multiplied by ( λ 2-k)

    techniques formeasuring nuclei distribution in aliquid,suchasan acousticbubble .However, nuclei spectrum measurement requires special equipment which was not available during the present model and prototype experiments. Therefore,the cumulative nuclei spectrum of the prototype is chosen as the reference line multiplied by 10,following Kamiirisa’s experimental result[16]. The adjusted cumulative nuclei spectrum of the model experiment is the dotted blue line in Fig.2, resulting from the reference line multiplied by ( Ds/Dm)2-kwhich is applied to assume that the number of bubbles passing through the tip vortex line of the model experiment is the same as prototype experiment. It is difficult to ascertain the vortex core center pressure in the model or the prototype precisely. However, the ratio of the vortex core center pressure between the model and the prototype can be calculated from Eq.(21) as 2.9. Therefore, if the vortex core center pressure of the model experiment is -0.4×105Pa, it can be calculated from Eq.(21) that the vortex core center pressure of the prototype experiment is-1.2×105Pa. Finally, the ratio of the number of nuclei flowing into the vortex line that is required to keep the same nuclei distribution in the model and theis determined to be 10.7 from Fig.2. This value is used in estimating the prototype PSD. The number of bubbles formed in the vortex line per unit time is inversely proportional to characteristic time (Eq.(25)).Growth time (τG) is the characteristic time in this system. Therefore, in TVC, the value of0( )N scales as follows

    1.6 Final rearrangement of the scaling law in TVC

    1.6.1 Frequency scaling

    Inserting Eq.(26) into Eq.(14), frequency scales as follows

    To factor in the effect of the cavitation number difference in the model and prototype of the same TVC shape, Eq.(31) is rearranged as follows

    Because of the dependence of the cavitation number on the Reynolds number (), the final form of frequency scaling is as follows

    This is consistent with the concept of blade passage frequency (BPF=)Bn× for the same number of blades in the model and the prototype, where B denotes number of propeller blades.

    1.6.2 PSD scaling

    Substituting Eq.(26) and Eq.(30) into Eq.(16),acoustic intensity scales as follows

    Here, it is assumed that water density ()ρ and sound speed0()c are same in both the model and the prototype experiment. Thus, water density and sound speed terms cancel out one another. As is done in the frequency scaling case, taking into consideration the effect of the cavitation number difference and the dependence of the cavitation number on the Reynolds number, Eq.(34) can be rearranged as follows

    Since PSD, G scales as

    where fΔ denotes frequency bandwidth. Equation(33) and Eq.(35) can be combined to finally arrive at the PSD scaling formula as follows

    2. Application and discussions

    In order to estimate the prototype PSD from the measured model data, a model scale experiment was performed in Samsung Ship Model Basin (SSMB)cavitation tunnel. To keep every variable but the cavitation number constant in the model experiment and to maintain the same thrust coefficient (kT) as the prototype in order to maintain dynamic similarity between the model and the prototype, the ambient pressure in the cavitation tunnel was adjusted. Sound pressures were measured using a flush-mounted hydrophone array on the hull above the propeller in the model and the prototype. The scaling exponent was determined from the cavitation numbers and the Reynoldsnumbersof the model and prototype experi-ments as

    The incipient cavitation number of the model and prototype was determined acoustically. Cavitation started to develop in the prototype at a higher cavitation number compared to that which was observed in the model experiment. Here, CIS in the prototype is 8.8 kts, determined with the acoustical signal processing technique. The TVC pattern at 8.8 kts of the prototype speed condition is most similar to the TVC pattern at the prototype 16 kts of the model speed condition. Thus, the model PSD at the prototype(16 kts speed condition) is used to estimate the prototype PSD at 8.8 kts for the validation of the scaling law derived in this study in comparison with the modified ITTC noise estimation rule. This observation is consistent with the fact that cavitation number for the same TVC shape in the model and the prototype is different. Therefore, to correct this disparity for the same TVC shape, a scaling exponent determined using this method is used to derive the estimate of the prototype PSD from the model PSD.The experimental details of the model and the prototype noise measurement are given in Ref.[18].

    The ITTC noise estimation rule can also be used to estimate prototype PSD if the cavitation number difference for The same TVC shape in the model and the prototype is reflected in the ITTC noise estimation rule through a scaling exponent[9]. For reference and clarity, the scaling laws for the modified ITTC noise estimation rule with viscous correction are given below.

    For frequency[9]

    For PSD[9]

    Prototype PSD estimated with the modified ITTC noise estimation rule with viscous correction is compared with the results estimated with the scaling law based on the physics of TVC.

    Table 1 Various ratios of propeller parameters

    Various ratios of propeller parameters are required to estimate the prototype PSD, as shown in Table 1. Model PSD is obtained through Welch’s method[19]and PSD levels are extracted using one third-octave band analysis.

    The results are shown in Fig.3. The solid line represents the measured prototype PSD, the squared line represents the estimated prototype PSD using the modified ITTC noise estimation rule with viscous correction (Eq.(39) and Eq.(40)), and the circled line represents the estimated prototype PSD with the scaling law derived in this paper (Eq.(33) and Eq.(37)).One interval in the vertical axis indicates 10 dB/Hz.Since the TVC noise is dominant in the frequency range around 1 kHz, it is necessary to compare the estimates in that range. In the frequency range from 0.5 kHz to 2.0 kHz, the novel scaling law-estimated PSD has a higher value than the prototype PSD. In the case of the model, the experiment is conducted in a cavitation tunnel. The geometry of the cavitation tunnel test sections typically leads to echoes and reverberation, and the tunnel’s prime mover, flowcontrol, and manipulation elements (turning vanes,flow straighteners, etc.) may unintentionally produce background noise at the same time[20]. On the other hand, the experiment for the prototype is conducted in the sea, where there is no background noise caused by echoes or reverberation. Therefore, it is plausible for PSD estimated from the model experiment conducted within the environment of a cavitation tunnel to be above the prototype PSD which is measured in an open environment. If the effect of the reverberant environment of a cavitation tunnel can be quantified more precisely, it should be possible to estimate the prototype PSD more precisely, reducing differences in the model and prototype experiments.

    Fig.3 (Color online) Estimated prototype PSD using both the modified ITTC noise estimation rule with viscous correction (squared line) and the derived novel scaling law law (circled line) along with the measured prototype PSD(solid line). One interval in the vertical axis indicates 10 dB/Hz

    Furthermore, since occurrence positions of the model and the prototype TVC do not perfectly coincide in a low axial velocity region of the wake field, the ratio of propagation distance in the model and the prototype will not be exactly proportional to the scale ratio. Thus, if there is a 15%± error in the ratio of the distance from the cavitation source to the receiver in the model and the prototype, the estimated prototype PSD might carry an error ranging from-1.4 dB to 1.2 dB.

    In this study, it is assumed that0N is proportional to inflow velocity and inversely proportional to propeller diameter, as shown in Eq.(30). However,this assumption may be incorrect due to the random process of cavitation inception. Therefore, the dependence of0N on inflow velocity and propeller diameter is to be verified through well-designed experiments in which control variables such as propeller parameters and environmental factors and nuclei distribution are properly maintained.

    3. Conclusions

    The modified ITTC noise estimation rule is used to extrapolate the prototype TVC noise from the model test data, employing viscous correction to account for the dependence of cavitation number on the Reynolds number. However, the viscous correction has no physical basis on facets of TVC such as the geometry of the inflow to the propeller tip and the vortex structure. In this study, the scaling law for propeller TVC noise is derived from much more physically sound hydrodynamics.

    The prototype PSD is extrapolated from the model data separately using the modified ITTC noise estimation rule with viscous correction and using the scaling law presented in this paper. The extrapolated results obtained using the novel scaling law show better agreement with the measured prototype PSD than the extrapolated results obtained using the modified ITTC noise estimation rule with viscous correction. In the vicinity of the TVC frequency region,extrapolated results tend to show a higher level of PSD than measured prototype PSD. This can be attributed to the fact that the model experiment is carried out in a more reverberant environment than the TVC conditions of the prototype. It is necessary to quantify the effects of the cavitation tunnel on noise levels, taking into account the detailed tunnel specificaions such as the geometry, material, and instruments in the tunnel and the source properties of the TVC.

    TVC is a complex physical phenomenon, so it is difficult to verify the structure of TVC through numerical or experimental methods. The novel scaling law derived in this paper suggests the possibility of verification for the physical phenomenon of TVC, since it provides improved accuracy in the estimation of prototype PSD.

    Acknowledgements

    This work was supported by Samsung Heavy Industries (SHI). The model experiment was carried out in Samsung Ship Model Basin and the prototype data was provided by SHI.

    [1] Strasberg M. Propeller cavitation noise after 35 years of study [C]. ASME Symposium on Noise and Fluids Engi- neering. Atlanta, USA, 1977.

    [2] Levkovskii V. L. Modeling of cavitation noise [J]. Soviet Physics Acoustics, 1968, 13(1): 337-339.

    [3] Alexandrov K., Bovis A. G., Huang T. T. ITTC Cavitation Committee Report [C]. 18th International Towing Tank Conference. Kobe, Japan, 1987, 183-190.

    [4] Bark G. Prediction of propeller cavitation noise from model tests and its comparison with full scale data [J].Journal of Fluids Engineering,1985, 107(1): 112-119.

    [5] Atlar M., Takinaci A. C., Korkut E. et al. Cavitation tunnel tests for propeller noise of a FRV and comparisons with full-scale measurements [C]. 4th International Symposium on Cavitation. Pasedena, USA, 2001.

    [6] Mccormick B. W. On cavitation produced by a vortex trailing from a lifting surface [J]. Journal of Basic Engineering, 1962, 84(3): 369-379.

    [7] Latorre R. TVC noise envelope - an approach to tip vortex cavitation noise scaling [J]. Journal of Ship Research,1982,26(1): 65-75.

    [8] Bjorheden O., Astrom L. Prediction of propeller noise spectra [C]. Detnorske Veritas Symposium on Hydrodyna- mics of Ship and Offshore Propulsion Systems. Oslo,Norway, 1977.

    [9] Oshima A. A. Study on correlation of vortex cavitation noise of propeller measured in model experiments and full scale [J]. J. the Society of Nav. Architects of Jpn.,1990,168: 89-96.

    [10] Shen Y. T., Gowing S., Jessup S. Tip vortex cavitation inception scaling for high Reynolds number application [J]. Journal of Fluids Engineering, 2009, 131(7): 1301-1306.

    [11] Arndt R. E. A. Nucleation and bubble dynamics in vortical flows [J]. Journal of Fluids Engineering, 2000,122(3): 488-493.

    [12] Blake W. K. Mechanics of flow-induced sound and vibration [M]. Orlando, USA: Academic Press, 1986, 370-396.

    [13] Lecoffre Y. Cavitation-bubble trackers [M]. Leiden, The Netherlands: Balkema, 1999, 76-114,154-158.

    [14] Lee S., Seo J. Cavitation tunnel test for propeller-induced noise and its comparison with full-scale measurements [C].15th International Congress on Sound and Vibration.Daejeon, Korea, 2008.

    [15] Park J., Park C., Choo Y. et al. Experimental study of scale effects on the scaling law for tip vortex cavitation noise[C]. 4th Joint US-European Fluids Engineering Division Summer Meeting. Chicago, USA, 2014.

    [16] Kamiirisa H. The effects of water quality characteristics on cavitation noise [C]. 4th International Symposium on Cavitation. Pasadena, USA, 2001.

    [17] Wu X., Wendel M., Chahine G. et al. Gas bubble size measurements in liquid mercury using an acoustic spectrometer [J]. Journal of Fluids Engineering, 2014,136(3): 031303.

    [18] Lee J. H., Jung J. K., Lee K. J. et al. Experimental estimation of a scaling exponent for tip vortex cavitation via its inception test in full- and model-ship [J]. Journal of Hydrodynamics, 2012, 24(5): 658-667.

    [19] Welch P. D. The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms [J]. IEEE Transactions on Audio Electroacoustics, 1967, 15(2):70-73.

    [20] Chang N. A., Dowling D. R. Ray-based acoustic localization of cavitation in a highly reverberant environment [J].Journal of Acoustical Society of America, 2009, 125(5):3088-3100.

    夜夜爽天天搞| 男女下面进入的视频免费午夜| 黄色视频,在线免费观看| 日本五十路高清| 国产成+人综合+亚洲专区| 尤物成人国产欧美一区二区三区| 最近在线观看免费完整版| 国产精品,欧美在线| 精品人妻1区二区| 网址你懂的国产日韩在线| 国产成人av激情在线播放| 成人av在线播放网站| 中文字幕熟女人妻在线| 男人舔奶头视频| 在线十欧美十亚洲十日本专区| 老司机福利观看| 俺也久久电影网| 舔av片在线| 精品乱码久久久久久99久播| 亚洲美女黄片视频| 国产精品一区二区三区四区免费观看 | 无限看片的www在线观看| 色av中文字幕| 国产视频内射| 好男人在线观看高清免费视频| 搡老岳熟女国产| 亚洲 欧美 日韩 在线 免费| 国产高清三级在线| 尤物成人国产欧美一区二区三区| 欧美一级毛片孕妇| 亚洲一区高清亚洲精品| 国产真实乱freesex| 成人午夜高清在线视频| 久久精品影院6| 欧美在线黄色| 欧美最新免费一区二区三区 | 久久久久九九精品影院| 国产伦人伦偷精品视频| 国产精品爽爽va在线观看网站| 香蕉av资源在线| 成人永久免费在线观看视频| 性欧美人与动物交配| www.www免费av| 国产三级黄色录像| 国产av不卡久久| www.色视频.com| 国产野战对白在线观看| 中文亚洲av片在线观看爽| 小蜜桃在线观看免费完整版高清| 免费看日本二区| 51国产日韩欧美| 精品无人区乱码1区二区| 十八禁人妻一区二区| 欧美一区二区国产精品久久精品| 热99在线观看视频| 男人舔奶头视频| 亚洲中文字幕日韩| 少妇熟女aⅴ在线视频| 国产欧美日韩精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 久久久久性生活片| 桃红色精品国产亚洲av| 日本a在线网址| 19禁男女啪啪无遮挡网站| 九九久久精品国产亚洲av麻豆| 精品国产三级普通话版| 国内揄拍国产精品人妻在线| 熟妇人妻久久中文字幕3abv| 国产v大片淫在线免费观看| 午夜福利视频1000在线观看| 我要搜黄色片| 久久香蕉国产精品| 久久精品国产亚洲av香蕉五月| 成人av一区二区三区在线看| 精品福利观看| 日韩欧美在线乱码| 午夜精品一区二区三区免费看| 国产精品久久久久久久久免 | 精品国产三级普通话版| 国产午夜精品久久久久久一区二区三区 | 五月玫瑰六月丁香| 亚洲片人在线观看| 伊人久久大香线蕉亚洲五| 国产探花极品一区二区| 国产精品1区2区在线观看.| 亚洲va日本ⅴa欧美va伊人久久| 看片在线看免费视频| 亚洲人成网站高清观看| 国产激情偷乱视频一区二区| 久久国产精品人妻蜜桃| 国产激情偷乱视频一区二区| 少妇的逼水好多| 看免费av毛片| 亚洲精品乱码久久久v下载方式 | 亚洲av成人精品一区久久| 亚洲精品456在线播放app | 国产高清视频在线播放一区| 国产乱人视频| 男人舔奶头视频| 欧美日韩黄片免| 国语自产精品视频在线第100页| av黄色大香蕉| 国产成人a区在线观看| 亚洲在线观看片| 精品一区二区三区av网在线观看| 国产在线精品亚洲第一网站| 亚洲av成人精品一区久久| 国产成人啪精品午夜网站| 一卡2卡三卡四卡精品乱码亚洲| x7x7x7水蜜桃| 亚洲国产色片| 麻豆国产av国片精品| 国产精华一区二区三区| 嫩草影院入口| 免费看十八禁软件| 嫩草影院入口| 欧美乱妇无乱码| 午夜亚洲福利在线播放| 99在线视频只有这里精品首页| 熟女人妻精品中文字幕| 亚洲男人的天堂狠狠| 窝窝影院91人妻| 少妇高潮的动态图| 日韩免费av在线播放| 精品久久久久久,| 久久精品夜夜夜夜夜久久蜜豆| 国产91精品成人一区二区三区| 国产高清视频在线观看网站| 亚洲,欧美精品.| 少妇的逼好多水| 亚洲七黄色美女视频| 免费人成视频x8x8入口观看| 俺也久久电影网| 国产精品美女特级片免费视频播放器| 成人av在线播放网站| 看免费av毛片| 免费一级毛片在线播放高清视频| 成熟少妇高潮喷水视频| 亚洲精品456在线播放app | ponron亚洲| 我的老师免费观看完整版| 精品人妻1区二区| 中文字幕人成人乱码亚洲影| 美女黄网站色视频| 九色成人免费人妻av| 极品教师在线免费播放| 亚洲av熟女| 久久精品国产99精品国产亚洲性色| 亚洲成a人片在线一区二区| svipshipincom国产片| 欧美精品啪啪一区二区三区| 精品一区二区三区人妻视频| 免费看a级黄色片| 丰满的人妻完整版| 国产一区二区三区在线臀色熟女| 亚洲精品影视一区二区三区av| 精品欧美国产一区二区三| 亚洲第一欧美日韩一区二区三区| 国产精品久久电影中文字幕| 99热这里只有是精品50| 一本综合久久免费| 一区二区三区激情视频| 国产三级在线视频| 日韩欧美 国产精品| 我的老师免费观看完整版| 国产精品综合久久久久久久免费| ponron亚洲| 一进一出好大好爽视频| 一进一出好大好爽视频| 免费电影在线观看免费观看| 床上黄色一级片| 亚洲18禁久久av| 69av精品久久久久久| 1024手机看黄色片| 制服丝袜大香蕉在线| 俄罗斯特黄特色一大片| 超碰av人人做人人爽久久 | 观看美女的网站| 亚洲精品久久国产高清桃花| 久久九九热精品免费| 国产亚洲av嫩草精品影院| 日韩欧美精品免费久久 | 露出奶头的视频| 国内毛片毛片毛片毛片毛片| 啪啪无遮挡十八禁网站| 亚洲精品久久国产高清桃花| 好男人在线观看高清免费视频| 亚洲人成伊人成综合网2020| 亚洲avbb在线观看| 中文字幕人妻熟人妻熟丝袜美 | 国产蜜桃级精品一区二区三区| 国产免费av片在线观看野外av| 97超视频在线观看视频| 波多野结衣高清作品| 久久精品综合一区二区三区| 黄色视频,在线免费观看| 欧美在线一区亚洲| 青草久久国产| 一进一出抽搐动态| 99久国产av精品| 精品国内亚洲2022精品成人| 精品久久久久久久人妻蜜臀av| 9191精品国产免费久久| 国产伦一二天堂av在线观看| av天堂中文字幕网| 人人妻,人人澡人人爽秒播| 日本 av在线| 19禁男女啪啪无遮挡网站| 国产精品久久久久久人妻精品电影| 久久天躁狠狠躁夜夜2o2o| 午夜精品在线福利| 久久草成人影院| av专区在线播放| 亚洲国产精品合色在线| 亚洲人成网站在线播放欧美日韩| ponron亚洲| 亚洲欧美日韩无卡精品| 欧美日本亚洲视频在线播放| 亚洲片人在线观看| 亚洲内射少妇av| 一个人免费在线观看电影| 最近最新免费中文字幕在线| 午夜激情福利司机影院| 亚洲不卡免费看| 老司机深夜福利视频在线观看| 亚洲 欧美 日韩 在线 免费| 99久久久亚洲精品蜜臀av| 国产午夜精品久久久久久一区二区三区 | 久久久久国内视频| 午夜精品久久久久久毛片777| 成熟少妇高潮喷水视频| 亚洲成av人片在线播放无| 精品国产亚洲在线| 亚洲片人在线观看| 免费观看人在逋| 非洲黑人性xxxx精品又粗又长| 色尼玛亚洲综合影院| 9191精品国产免费久久| 欧美精品啪啪一区二区三区| 久久香蕉精品热| 2021天堂中文幕一二区在线观| 精品国产三级普通话版| 久久久久久大精品| 51午夜福利影视在线观看| 国产欧美日韩一区二区精品| 欧美另类亚洲清纯唯美| 一级毛片女人18水好多| 手机成人av网站| 神马国产精品三级电影在线观看| 国产精品香港三级国产av潘金莲| 亚洲av美国av| 欧美在线一区亚洲| 99在线视频只有这里精品首页| 精品无人区乱码1区二区| 法律面前人人平等表现在哪些方面| 日日摸夜夜添夜夜添小说| 欧美激情久久久久久爽电影| 熟女电影av网| 国产av在哪里看| 亚洲内射少妇av| 久久久色成人| 级片在线观看| 三级男女做爰猛烈吃奶摸视频| 男人舔女人下体高潮全视频| 国产精品爽爽va在线观看网站| 亚洲国产色片| 国产成人av激情在线播放| 亚洲成人免费电影在线观看| avwww免费| 亚洲欧美日韩卡通动漫| 老司机午夜福利在线观看视频| 欧美大码av| 午夜免费男女啪啪视频观看 | 天天一区二区日本电影三级| 久久久久久人人人人人| 婷婷六月久久综合丁香| 欧美日韩中文字幕国产精品一区二区三区| 中亚洲国语对白在线视频| 日本撒尿小便嘘嘘汇集6| 国产色婷婷99| 黄色成人免费大全| 91av网一区二区| 亚洲无线在线观看| 最近最新免费中文字幕在线| 国产v大片淫在线免费观看| 日韩大尺度精品在线看网址| 成人国产综合亚洲| 精品国产亚洲在线| 夜夜看夜夜爽夜夜摸| 欧美+亚洲+日韩+国产| 香蕉av资源在线| а√天堂www在线а√下载| 欧美日韩福利视频一区二区| svipshipincom国产片| 欧美乱色亚洲激情| 欧美一区二区精品小视频在线| 欧美黄色片欧美黄色片| 夜夜爽天天搞| 又爽又黄无遮挡网站| 日韩欧美在线乱码| 免费在线观看影片大全网站| 精品国产美女av久久久久小说| 搡老妇女老女人老熟妇| 一进一出抽搐动态| 精品一区二区三区av网在线观看| 一个人看的www免费观看视频| 亚洲美女黄片视频| 麻豆国产97在线/欧美| 亚洲成av人片免费观看| 99久久成人亚洲精品观看| 一个人看视频在线观看www免费 | 老鸭窝网址在线观看| 免费看日本二区| 亚洲成av人片在线播放无| 91字幕亚洲| 久久精品影院6| 国产亚洲精品久久久com| 高清在线国产一区| 少妇高潮的动态图| 亚洲成人精品中文字幕电影| or卡值多少钱| 国产精品美女特级片免费视频播放器| 国产成人a区在线观看| 欧美日韩综合久久久久久 | 欧美一区二区国产精品久久精品| 青草久久国产| 欧美高清成人免费视频www| 亚洲精品成人久久久久久| 少妇裸体淫交视频免费看高清| 色视频www国产| 亚洲成av人片在线播放无| 国产av在哪里看| a级毛片a级免费在线| 色视频www国产| 757午夜福利合集在线观看| 亚洲欧美日韩无卡精品| 亚洲性夜色夜夜综合| 国产日本99.免费观看| 久久精品国产亚洲av涩爱 | 免费在线观看成人毛片| 91在线观看av| 日韩人妻高清精品专区| 淫秽高清视频在线观看| 欧美日韩一级在线毛片| 久久久久久久久中文| 精品人妻1区二区| 欧美性猛交黑人性爽| 51国产日韩欧美| 一级毛片女人18水好多| 日本黄大片高清| 综合色av麻豆| 国产精品 欧美亚洲| 国产乱人视频| 国产精品精品国产色婷婷| 久久久国产成人精品二区| 熟妇人妻久久中文字幕3abv| 黄色女人牲交| 国产精品久久久久久人妻精品电影| 国产成人av教育| 亚洲国产欧美网| 久久香蕉精品热| 国产一区在线观看成人免费| 最近最新中文字幕大全免费视频| 欧洲精品卡2卡3卡4卡5卡区| 国产免费一级a男人的天堂| 亚洲熟妇中文字幕五十中出| 午夜激情福利司机影院| 国产探花极品一区二区| 免费av毛片视频| 真人一进一出gif抽搐免费| 久久久久久久久中文| 亚洲 国产 在线| 超碰av人人做人人爽久久 | 日本在线视频免费播放| 麻豆久久精品国产亚洲av| 一二三四社区在线视频社区8| 操出白浆在线播放| 精品久久久久久久人妻蜜臀av| 日韩欧美国产一区二区入口| 色综合站精品国产| 老司机午夜福利在线观看视频| 久久香蕉精品热| 国产精品三级大全| 国产日本99.免费观看| 亚洲一区二区三区不卡视频| 日韩欧美精品v在线| 日韩av在线大香蕉| 99热这里只有精品一区| 国语自产精品视频在线第100页| 丰满人妻一区二区三区视频av | 首页视频小说图片口味搜索| 国产午夜精品论理片| 18禁黄网站禁片免费观看直播| 国产一区二区三区在线臀色熟女| 观看免费一级毛片| 亚洲乱码一区二区免费版| 亚洲精品国产精品久久久不卡| 在线观看日韩欧美| 国产精品98久久久久久宅男小说| 黄片小视频在线播放| 禁无遮挡网站| netflix在线观看网站| 国产精品香港三级国产av潘金莲| 久久久久国内视频| 亚洲中文日韩欧美视频| 精华霜和精华液先用哪个| 免费看十八禁软件| 免费在线观看日本一区| 亚洲精品一区av在线观看| 日本a在线网址| 日韩欧美国产在线观看| 久久久久亚洲av毛片大全| 91在线精品国自产拍蜜月 | 天堂动漫精品| 色视频www国产| 变态另类成人亚洲欧美熟女| 亚洲熟妇熟女久久| 色尼玛亚洲综合影院| 午夜免费成人在线视频| 一个人观看的视频www高清免费观看| 亚洲在线观看片| 最好的美女福利视频网| 色精品久久人妻99蜜桃| 精品久久久久久久人妻蜜臀av| 亚洲精品成人久久久久久| 美女cb高潮喷水在线观看| 久久久久久久久中文| 精品久久久久久,| 精品久久久久久久末码| 午夜日韩欧美国产| 天堂√8在线中文| 久久久久性生活片| 国产精品1区2区在线观看.| 长腿黑丝高跟| 国产成人欧美在线观看| 身体一侧抽搐| av视频在线观看入口| 黑人欧美特级aaaaaa片| 人人妻,人人澡人人爽秒播| 免费人成视频x8x8入口观看| 丁香六月欧美| 国产一区在线观看成人免费| 青草久久国产| 特级一级黄色大片| 国产av不卡久久| 欧美一级毛片孕妇| 99视频精品全部免费 在线| 极品教师在线免费播放| 国产色爽女视频免费观看| 亚洲成人精品中文字幕电影| 国产成+人综合+亚洲专区| 在线天堂最新版资源| 亚洲av日韩精品久久久久久密| 高清在线国产一区| 女人高潮潮喷娇喘18禁视频| 一本一本综合久久| 变态另类成人亚洲欧美熟女| 18禁在线播放成人免费| 亚洲在线自拍视频| 久久国产乱子伦精品免费另类| 麻豆国产av国片精品| 伊人久久精品亚洲午夜| 精品久久久久久久毛片微露脸| 精品久久久久久久久久免费视频| 草草在线视频免费看| x7x7x7水蜜桃| 亚洲电影在线观看av| 搡老妇女老女人老熟妇| 精品一区二区三区av网在线观看| 欧美av亚洲av综合av国产av| 亚洲欧美日韩高清专用| 一个人看的www免费观看视频| 乱人视频在线观看| 日本在线视频免费播放| 日本与韩国留学比较| 真人一进一出gif抽搐免费| 少妇熟女aⅴ在线视频| 欧美一区二区国产精品久久精品| 两个人看的免费小视频| 国产高清激情床上av| 成人鲁丝片一二三区免费| 悠悠久久av| 18禁裸乳无遮挡免费网站照片| 婷婷精品国产亚洲av| av在线天堂中文字幕| 成人精品一区二区免费| 高清日韩中文字幕在线| 美女高潮的动态| av在线蜜桃| 一区二区三区免费毛片| 两性午夜刺激爽爽歪歪视频在线观看| 色综合站精品国产| 真实男女啪啪啪动态图| 日韩成人在线观看一区二区三区| tocl精华| 国语自产精品视频在线第100页| 欧美性猛交╳xxx乱大交人| 久久国产精品人妻蜜桃| 中文字幕人成人乱码亚洲影| 亚洲成人久久爱视频| 不卡一级毛片| 99久久无色码亚洲精品果冻| 天堂av国产一区二区熟女人妻| 久久久成人免费电影| 亚洲欧美精品综合久久99| 久久久久久人人人人人| 夜夜爽天天搞| 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 两人在一起打扑克的视频| av在线蜜桃| 俺也久久电影网| 激情在线观看视频在线高清| 色综合站精品国产| 国产在线精品亚洲第一网站| avwww免费| 色视频www国产| 精品人妻偷拍中文字幕| 给我免费播放毛片高清在线观看| 国产老妇女一区| 亚洲黑人精品在线| 内射极品少妇av片p| 亚洲无线在线观看| 亚洲在线观看片| 国产精品久久久人人做人人爽| 亚洲成a人片在线一区二区| 精品电影一区二区在线| 禁无遮挡网站| 又紧又爽又黄一区二区| 久久6这里有精品| 老司机深夜福利视频在线观看| 午夜福利免费观看在线| 欧美一级a爱片免费观看看| 黄色女人牲交| 特大巨黑吊av在线直播| 亚洲精品在线观看二区| 国产淫片久久久久久久久 | 1024手机看黄色片| 九九在线视频观看精品| 69av精品久久久久久| 国产精品99久久久久久久久| 久久久国产精品麻豆| 脱女人内裤的视频| 欧美黄色片欧美黄色片| 久久伊人香网站| 黄色女人牲交| 校园春色视频在线观看| 国产真实伦视频高清在线观看 | 国产午夜福利久久久久久| 免费看十八禁软件| h日本视频在线播放| 亚洲欧美一区二区三区黑人| 在线播放无遮挡| 久久久久久大精品| 美女被艹到高潮喷水动态| 波多野结衣高清作品| 99热精品在线国产| 美女黄网站色视频| 国产高清videossex| 国产精品av视频在线免费观看| 中出人妻视频一区二区| 亚洲成人精品中文字幕电影| 老司机午夜十八禁免费视频| 亚洲专区中文字幕在线| 亚洲国产色片| ponron亚洲| 操出白浆在线播放| 国产高清视频在线播放一区| 最近在线观看免费完整版| 国内精品美女久久久久久| 久久午夜亚洲精品久久| 全区人妻精品视频| www.www免费av| 欧美最新免费一区二区三区 | 欧美日本亚洲视频在线播放| 在线免费观看不下载黄p国产 | 别揉我奶头~嗯~啊~动态视频| 亚洲 欧美 日韩 在线 免费| 久久性视频一级片| 国产一级毛片七仙女欲春2| 午夜福利视频1000在线观看| 亚洲成人免费电影在线观看| 亚洲精品久久国产高清桃花| 欧美最新免费一区二区三区 | 国产精品香港三级国产av潘金莲| 亚洲国产精品久久男人天堂| 国产91精品成人一区二区三区| 亚洲天堂国产精品一区在线| 久久精品国产清高在天天线| 亚洲在线自拍视频| 变态另类丝袜制服| 99久久精品一区二区三区| 久久久国产精品麻豆| 欧美高清成人免费视频www| 亚洲狠狠婷婷综合久久图片| 99久久99久久久精品蜜桃| 舔av片在线| 欧美日韩国产亚洲二区| 国内毛片毛片毛片毛片毛片| 日韩国内少妇激情av| 91久久精品国产一区二区成人 | 精华霜和精华液先用哪个| 午夜福利成人在线免费观看| 欧美bdsm另类| 久久伊人香网站| 欧美黄色片欧美黄色片| 亚洲av成人精品一区久久| 亚洲av成人av| 免费一级毛片在线播放高清视频| a级毛片a级免费在线| 香蕉久久夜色| 日本a在线网址| 老司机在亚洲福利影院| 久久久久九九精品影院| 精品人妻偷拍中文字幕|