EFFECT OF A PROPELLER AND GAS DIFFUSION ON BUBBLE NUCLEI DISTRIBUTION IN A LIQUID

A multi-bubble dynamics code accounting for gas diffusion in the liquid and through the bubble wall was developed and used to study the modification of a bubble nuclei population dynamics by a propeller.The propeller flow field was obtained using a Reynolds-Averaged Navier-Stokes(RANS) solver and bubble nuclei populations were propagated in this field.The numerical pro-cedure enabled establishment of the possibility of production behind the propeller of relatively large visible bubbles starting from typical ocean nuclei size distributions.The resulting larger bubbles are seen to cluster in the blade wakes and tip vortices.Parametric investigations of the initial nuclei size distribution,the dissolved gas concentration,and the cavitation number were conducted to ide-ntify their effects on bubble entrainment and the resultant void fractions and bubble distribution modifications downstream from the propeller.Imposed synthetic turbulence-like fluctuations unto the average RANS flow field were also used to study the effect avera-ging in the RANS procedure has on the results.     

EXPERIMENTAL ESTIMATION OF A SCALING EXPONENT FOR TIP VORTEX CAVITATION VIA ITS INCEPTION TEST IN FULL- ANDMODEL-SHIP

Tip vortex cavitation noise of marine propeller became primary concerns to reduce hazardous environmental impacts from commercial ship or to keep the underwater surveillance of naval ships. The investigations of the tip vortex and its induced noise are normally conducted through the model test in a water cavitation tunnel. However the Reynolds number of model-test is much smaller than that of the full-scale, which subsequently results in the difference of tip vortex cavitation inception. Hence, the scaling law between model-and full-scales needs to be identified prior to the prediction and assessment of propeller noise in full scale. From previous researches, it is generally known that the incipient caivtation number of tip vortex can be represented as a power of the Reynolds number. However, the power exponent for scaling, which is the main focus of this research, has not been clearly studied yet. This paper deals with the estimation of scaling exponent based on tip vortex cavitation inception test in both full-and model-scale ships. Acoustical measurements as well as several kind of signal processing technique for an inception criterion suggest the scaling exponent as 0.30. The scaling value proposed in this study shows slight difference to the one of most recent research. Besides, extrapolation of model-ship noise measurement using the proposed one predicts the full-scale noise measurement with an acceptable discrepancy.     

NUMERICAL INVESTIGATION OF CAVITATION PERFORMANCE OF SHIP PROPELLERS

The cavitation performance of propellers is studied based on viscous multiphase flow theories. With a hybrid grid based on Navier-Stokes (N-S) and bubble dynamics equations, some recent validation results are presented in this paper in the predictions of the thrust, the torque and the vapor volume fraction on the back side of propeller blade for a uniform inflow. The numerical predictions of the hydrodynamic performance and the sheet cavitation under several operating conditions for two propellers agree with the corresponding measured data in general. The thrust and the torque are plotted with respect to the advance rate and the cavitation number. The cavitation performance breakdown is closely related to the strong sheet cavitation around propellers. The models with parameters modified are shown to predict the propeller cavitation well.     

Experimental measurement of tip vortex flow field with/without cavitation in an elliptic hydrofoil

In this paper, recent measurements of tip vortex flow with and without cavitation carried out in Cavitation Mechanism Tunnel of China Ship Scientific Research Center(CSSRC) are presented. The elliptic hydrofoil with section NACA 662-415 was adopted as test model. High-speed video(HSV) camera was used to visualize the trajectory of tip vortex core and the form of tip vortex cavitation(TVC) in different cavitation situations. Laser Doppler velocimetry(LDV) was employed to measure the tip vortex flow field in some typical sections along the vortex trajectory with the case of cavitation free. Stereo particle image velocimetry(SPIV) system was used to measure the velocity and vorticity distributions with and without cavitation. Series measurement results such as velocity and vorticity distributions, the trajectory of tip vortex core, the vortex core radius, cavity size and cavitation inception number were obtained. The results demonstrated that the minimum pressure coefficient in the vortex core obtained by flow field measurement was quite coincident with the tip vortex cavitation inception number obtained under the condition of high incoming velocity and low air content. And TVC would decrease the vortex strength comparing with the case without cavitation.     

LES NUMERICAL SIMULATION OF CAVITATION BUBBLE SHEDDING ON ALE 25 AND ALE 15 HYDROFOILS

A cavitation calculation scheme is developed and applied to ALE 15 and ALE 25 hydrofoils, based on the Bubble Two-phase Flow (BTF) cavity model with a Large Eddy Simulation (LES) methodology. The Navier-Stokes equations including cavitation bubble clusters are solved through the finite-volume approach with a time-marching scheme. Simulations are carried out in a 3-D field with a hydrofoil ALE 15 or ALE 25 at an angle of attack of 8° and cavitation number σ = 2.3 with a 2 × 10⁶, meshing system. With the time-marching, the cavitation bubble gradually grows to a steady lump shape and then produces an irregular small bubble behind the main cavitation bubble, finally shedding from the leading edge of the cloud cavitation structure. The calculated results including velocity field and pressure field are consistent with experiment data at the same Reynolds number and cavitation number. The vortex and reverse flow are observed on the hydrofoil surface.     

The effect of water quality on tip vortex cavitation inception

The inception of tip vortex cavitation is very sensitive to water quality. In order to quantify the effect of water quality on the inception of tip vortex cavitation, we develop a motion model to describe the migration and growth of nuclei in water. An analytical solution of migration of nuclei in a vortex flow is obtained so that the capture times of various nuclei can be given out directly. A criterion is built to determine the critical nucleus in a certain nuclei spectra distribution. Tensile strength of the critical nucleus is used to quantify the effect of water quality and correct the tip vortex cavitation inception number. Finally this change of cavitation inception number is compared with experimental results to validate our model.     

URANS COMPUTATION OF CAVITATING FLOWS AROUND SKEWED PROPELLERS

Cavitating flows around skewed propellers are investigated numerically by means of the unsteady Reynolds Averaged Navier-Stokes (RANS) Equation method. The standard k-ε turbulence and the modified Z-G-B cavitation models are employed. A measured nominal wake is used for the inlet velocity boundary condition. Predicted cavitating evolution processes and tip cavity patterns are compared with experimental observations. In addition, the influence of the skew angles on the cavitation and unsteadiness performances of propellers operating in a non-uniform wake is also studied. Results show that the modified Z-G-B cavitation model performs better to simulate the cavitating flow cases studied in this paper. Comparisons demonstrate that the skewed propeller with a skew angle of 20 o is the best choice for a given stern wake with a assigned thrust and the minimum force fluctuations.     

非均匀进流对螺旋桨空化水动力性能的影响

以全附体SUBOFF潜艇桨盘面处标称伴流作为七叶大侧斜桨的非均匀进流,分析了对空化引起的推力和力矩崩溃性能以及叶背梢涡空化初生的影响,描述了非均匀进流条件下空化时推力和力矩的脉动特征、桨叶空化面积和空化形态随周向位置的变化,并给出了伴流中螺旋桨是否出现叶面片空化的区间划分。空化模拟时采用改进Sauer空化模型和修正SST湍流模型。梢涡空化初生依据作者提出的当σ>σi时,叶截面压力系数分布不再改变的相对准则来判定。结果表明,伴流能够有效延迟该桨叶面片空化初生,但同时会使叶背片空化初生提前且空化范围明显增加;在对应敞水效率最高的工况点,伴流使得该桨因空化引起的推力下降起始点较均匀来流时提前量大于7.1%,且叶背梢涡空化初生提前量大于6.7%;即使在空化较为严重的情况下,该桨在伴流场中的推力和力矩脉动幅值也很小;由于桨叶旋向的影响,一周内桨叶空化面积对于艇尾左右舷空间并不完全对称。研究结果可用于分析大侧斜使艇尾桨空化初始航速提高一定程度的现象。     

PREDICTION OF UNSTEADY CAVITATION OF PROPELLER USING SURFACE-PANEL METHOD

This paper has predicted the range and volume of unsteady sheet cavitation of a propeller by using the surface panel method. The linearization in cavity thickness is adopted to reduce the computing time and storage space. The iteration scheme between chordwise strips has been used because the range and volume of cavitation are both unknown. The propeller cavitation range determined by the calculation method presented in this paper agrees with the observation results of cavity image at cavitation tunnel very well, and this proves the practicability of the method.     

NUMERICAL PREDICTION OF BLADE FREQUENCY NOISE OF CAVITA- TING PROPELLER

The blade frequency noise of a cavitating propeller in a uniform flow is analyzed in the time domain. The unsteady loading (of a dipole source) and the sheet cavity volume (of a monopole source) on the propeller surface are calculated by a potential-based surface panel method. Then the time-dependent pressure and the cavity volume data are used as the input for the Fowcs Williams-Hawkings formulation to predict the acoustics pressure. The integration of the noise source is performed over the true blade surface rather than the ideal blade surface without thickness. The noise characteristics of the cavitating propeller are discussed. With the sheet cavitation, the thickness (cavitation) noise is larger than the loading noise and is the dominant noise source.The noise directivity is not as clear as that of the noise under a non-cavitation condition. The cavitation noise is attenuated mores lowly than the non-cavitation noise.     

NUMERICAL ANALYSIS OF THE INFLUENCE OF THE TIP CLEARANCE FLOWS ON THE UNSTEADY CAVITATING FLOWS IN A THREE- D IMENSIONAL INDUCER

To clarify the influences of the tip clearance flows on the unsteady cavitating flow,the three-dimensional unsteady cavitating flows through both the two-dimensional cascades and the three-dimensional inducer with and without tip clearance are performed numerically.The governing equations for the compressible fluid flow with the DES turbulence model are employed with the assumption of the isentropic process of liquid phase.The evolution of cavities is represented as the source/sink of vapor phase.The basic equations in the curve linear coordinate are solved by the finite difference method.As the results of the three-dimensional cavitating flows through the two-dimensional cascades,the tip clearance flows from the pressure side to the suction side of the blade produces the tip vortex cavitation,which affects the sheet cavitation on the leading edge of the next blade and enhances the blockage effect near the casing than the flows without tip clearance.On the other hand,in the case of the three-dimensional inducer,the large backflow cavitation is observed around the inlet of the inducer,where the cavities are developed on the casing by the tip clearance flows.The large pressure gradient between the non-cavitating pressure side and the cavitating suction side enhances the tip clearance flows.The calculation considering the tip clearance reproduces the developed cavitation region similar to that of experimental visualizations.Additionally,the backflow cavitation rotates with the speed slower than the rotation speed of the inducer.Then,the rotation of backflow cavitation causes the periodic fluctuation of the outlet pressure greater than that of the inlet pressure.     

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

The tip vortex cavitation(TVC) noise of marine propellers is of interest due to the environmental impacts from commercial ships as well as for the survivability of naval ships. Due to complicated flow and noise field around a marine propeller, a theoretical approach to the estimation of TVC noise is practically unrealizable. Thus, estimation of prototype TVC noise level is realized through extrapolation of the model TVC noise level measured in a cavitation tunnel. In this study, for the prediction of prototype TVC noise level from a model test, a novel scaling law reflecting the physical basis of TVC is derived from the Rayleigh-Plesset equation, the Rankine vortex model, the lifting surface theory, and other physical assumptions. Model and prototype noise data were provided by Samsung Heavy Industries(SHI) for verification. In applying the novel scaling law, similitude of the spectra of nuclei is applied to assume the same nuclei distribution in the tip vortex line of the model and the prototype. It was found that the prototype TVC noise level predicted by the novel scaling law has better agreement with the prototype TVC noise measurement than the prototype TVC noise level predicted by the modified ITTC noise estimation rule.     

高负荷螺旋桨水动力性能试验研究

以某大型集装箱船高负荷螺旋桨为研究对象,分别采用大侧斜、反弯扭叶片和毂帽鳍等技术进行螺旋桨优化设计,以改善螺旋桨的空泡和水动力性能.设计了5型螺旋桨,分别在空泡水筒进行了模拟伴流场中的空泡和激振力试验,并在拖曳水池进行了敞水试验.试验结果表明,采用大侧斜可有效减小螺旋桨诱导船艉激振力;反弯扭叶片螺旋桨空泡性能较好,有较低的螺旋桨诱导船艉激振力水平;桨毂帽鳍可改善毂涡空泡,模型试验状态下可提高推进效率2%左右.     

VISCOUS SCALE EFFECTS ON PROPELLER TIP VORTEX CAVITATION

ＶＩＳＣＯＵＳ　ＳＣＡＬＥ　ＥＦＦＥＣＴＳ　ＯＮ　ＰＲＯＰＥＬＬＥＲ　ＴＩＰ　ＶＯＲＴＥＸ　ＣＡＶＩＴＡＴＩＯＮＶＩＳＣＯＵＳＳＣＡＬＥＥＦＦＥＣＴＳＯＮＰＲＯＰＥＬＬＥＲＴＩＰＶＯＲＴＥＸＣＡＶＩＴＡＴＩＯＮ￥ＹｅＹｕａｎ－ｐｅｉ；ＬｕＦａｎｇ；Ｓ...     

QUADRANT ANALYSIS OF BUBBLE INDUCED VELOCITY FLUCTUATIONS IN A TRANSITIONAL BOUNDARY LAYER

Our previous study showed that the frictional drag decreases with increasing void fraction at Re > 1300, while it increases at Re < 1000. Decomposition of the Reynolds shear stress also implied that bubbles induce isotropy of turbulence. In order to confirm our previous analysis and to further investigate flow fields in the vicinity of bubbles, we analyze velocity fluctuations on the quadrant space in the streamwise and transverse directions ( u′ – v′ plane). Here, we focus on two specific Reynolds numbers ...     

Experimental study of the interaction between the spark-induced cavitation bubble and the air bubble

Experiments are carried out by using high-speed photography to investigate the interaction between the spark-generated cavitation bubble and the air bubble in its surrounding fluid. Three problems are discussed in detail： the impact of the air bubble upon the development of the cavitation bubble, the evolution of the air bubble under the influence of the cavitation bubble, and the change of the fluid pressure during the development of a micro jet of the cavitation bubble. Based on the experimental results, under the condition of no air bubble present, the lifetime of the cavitation bubble from expansion to contraction increases with the increase of the maximum radius. On the other hand, when there is an air bubble present, different sized cavitation bubbles have similarity with one another generally in terms of the lifetime from expansion to contraction, which does not depend on the maximum radius. Also, with the presence of an air bubble, the lifetime of the smaller cavitation bubble is extended while that of the bigger ones reduced. Furthermore, it is shown in the experiment that the low pressure formed in the opposite direction to the cavitation bubble micro jet makes the air bubble in the low pressure area being stretched into a steplike shape.     

On the numerical simulations of vortical cavitating flows around various hydrofoils

This paper reviews the numerical models of various cavitating flows around hydrofoils. Numerical models relating to cavitation flows, including mass transfer models and turbulence models, are summarized at first. Then numerical results and analysis of flow characteristics for the cavitating flows around twisted hydrofoils, truncated hydrofoil and tip leakage are discussed respectively. For mean flow fields, Reynolds averaged Navier-Stokes(RANS) simulation associated with a kind of nonlinear turbulence model is found to be an economic and robust numerical approach for different kinds of cavitating flows including cloud cavitation, tip cavitation and tip leakage cavitation. To predict the fluctuations of pressure and velocity, large eddy simulation(LES) and detached eddy simulation(DES) are two effective approaches. Finally, a few open questions are proposed for future research.     

NUMERICAL SIMULATION OF BUBBLE FLOW INTERACTIONS

Bubble flow interaction can be important in many practical engineering applications. For instance, cavitation is a problem of interaction between nuclei and local pressure field variations including turbulent oscillations and large scale pressure variations. Various types of behaviours fundamentally depend on the relative sizes of the nuclei and the length scales of the pressure variations as well as the relative importance of bubble natural periods of oscillation and the characteristic time of the field pressure variations. Similarly, bubbles can significantly affect the performance of lifting devices or propulsors. We present here some fundamental numerical studies of bubble dynamics and deformation, then a practical method using a multi-bubble Surface Averaged Pressure (DF-Multi-SAP©) to simulate cavitation inception and scaling, and connect this with more precise 3-D simulations. This same method is then extended to the study of two-way coupling between a viscous compressible flow and a bubble population in the flow field.     

PREDICTION OF PODDED PROPELLER CAVITATION USING AN UNSTEADY SURFACE PANEL METHOD

The unsteady sheet cavitation of podded propeller was predicted by using a surface panel method. The interaction between propeller and pod was treated with the iterative calculation of induced velocity potential, and the method of induced velocity potential can save a great deal of storage and computation time compared to the method of induced velocity. The induced velocity potential of unit singularity on every pod panel to every key blade panel and of unit singularity on every key blade panel and its wake panel to every pod panel were calculated when the key blade is at every angle position. Based on the wake model of the conventional single propeller, a new wake model of podded propeller was constructed. The propeller is analyzed only on the key blade in order to save computation time and memory space. The method can be used to calculate the hydrodynamics performance and cavitation of propeller in uniform and non-uniform inflows. It can give the unsteady force and cavitation shape of propeller. The propeller cavitation range determined by the present method agrees with the observation results of cavity image given in cavitation tunnel well, and this proves the practicability of the method.     

INTERACTION OF A CAVITATION BUBBLE AND AN AIR BUBBLE WITH A RIGID BOUNDARY

The motion of a spark-induced cavitation bubble and an air bubble near a rigid boundary is experimentally studied by using high-speed photography. Several dimensionless parameters are used to describe the geometrical configuration of the bubble-bubble-boundary interaction. The bubble-bubble interaction can be considered in two different conditions. The cavitation bubble will collapse towards the air bubble if the air bubble is relatively small, and away from the air bubble if the air bubble is relatively large. The two zones are identified in the bubble-boundary interaction, and they are the danger zone and the safety zone. The relative position, the bubble-boundary distance and the bubble-bubble distance play important roles in the bubble-bubble-boundary interaction, which can be considered in several conditions according to the responses of the bubbles. Air jets are found to penetrate into the cavitation bubbles. The cavitation bubble and the air bubble (air jet) move in their own way without mixing. The motion of a cavitation bubble may be influenced by an air bubble and/or a rigid boundary. The influence of the air bubble and the influence of the boundary may be combined, like some thing of a vector.