Experimental evaluation of numerical simulation of cavitating flow around hydrofoil

Cavitation in hydraulic machines causes different problems that can be related to its unsteady nature. An experimental and numerical study of developed cavitating flow was performed. Until now simulations of cavitating flow were limited to the self developed “in house” CFD codes. The goal of the work was to experimentally evaluate the capabilities of a commercial CFD code (Fluent) for simulation of a developed cavitating flow. Two simple hydrofoils that feature some 3D effects of cavitation were used for the experiments. A relatively new technique where PIV method combined with LIF technique was used to experimentally determine the instantaneous and average velocity and void ratio fields (cavity shapes) around the hydrofoils. Distribution of static pressure on the hydrofoil surface was determined. For the numerical simulation of cavitating flow a bubble dynamics cavitation model was used to describe the generation and evaporation of vapour phase. An unsteady RANS 3D simulation was performed. Comparison between numerical and experimental results shows good correlation. The distribution and size of vapour structures and the velocity fields agree well. The distribution of pressure on the hydrofoil surface is correctly predicted. The numerically predicted shedding frequencies are in fair agreement with the experimental data.     

冲蚀磨损与冲蚀、空蚀交互磨损的对比研究

基于计算流体力学(computational fluid dynamics,CFD)方法,数值对比研究了在模拟水轮机工况中不同转速时,冲蚀与冲蚀、空蚀交互作用时,转盘表面的流场(压力场、流体浓度分布),然后在转盘式磨损装置上,进行了汽液固三相冲蚀与空蚀交互磨损试验.结果表明:对于冲蚀磨损来说,随着转速增加,压力增大,在冲蚀作用下加入空蚀磨损,压力增高,磨损加剧;而对于交互磨损来说,随着转速增加,转盘表面的空蚀磨损区域从空化孔附近开始沿着转盘旋转的反方向偏移,并且磨损程度加剧;数值计算的气泡较多且压力梯度较高区域和试验转盘磨损区域基本一致,数值计算结果和试验结果吻合得较好.     

射流对绕水翼云空化流动抑制机理研究

为理解绕水翼云空化流动的发展机理和探究水翼吸力面开孔射流的影响,采用密度 修正的RNG $k$-$\varepsilon $湍流模型和Schnerr-Sauer空化模型对原始NACA66(mod) 水翼和采用射流后的 水翼的云空化非定常过程进行模拟和对比分析;采用在水翼吸力面近壁区设立监测线的方法对近壁区的流场进行监测,得到 近壁区汽相体积分数、回射流速度、压力及压力梯度的时空分布云图;开展了云空化流场特性的涡动力学分析,进而分析水 翼云空化的发生机理和射流抑制空化的抑制机理. 结果表明:游离型空泡在下游溃灭时产生强烈的局部高压,其向上游传播 导致前缘空穴的一次回缩,而空穴的二次回缩受回射流的影响. 回射流的发展区域受限于较高的压力梯度,高的压力梯度一 直存在,但回射流在一个周期内的首次出现需要时间的积累. 在水翼吸力面射流使得射流孔附近压力升高,弥补了由于空化 和绕流造成的压降,压力梯度增大,抗逆压能力增强,对回射流起到阻挡作用;另一方面,射流使得回射流区域面积和回射 流的强度也有所减小,从而对云空化的发展起到抑制的效果. $Q$准则的涡结构云图相比于汽相体积分数云图能显示复杂的 流动结构,前缘附着型空穴和尾缘游离型空穴内存在旋涡,回射流对空穴存在剪切作用造成空穴脱落. 而射流对空穴和回射 流的剪切和阻挡使云空化发展得到抑制.      

ANALYSIS OF UNSTEADY CAVITATION FLOW OVER HYDROFOIL BASED ON DYNAMIC MODE DECOMPOSITION$^{\bf 1)}$

空化是船舶和水下航行体推进器中经常发生的一种特殊流动现象,它具有强烈的非定常性,空化的发生往往会影响推进器的水动力性能和效率. 为探究绕水翼非定常空化流场结构,本文基于 Schnerr-Sauer 空化模型和 SST $k$-$\omega $ 湍流模型,开展绕二维水翼非定常空化流动数值预报与流场结构分析. 通过将数值预报的空泡形态演变和压力数据与试验结果对比,验证了建立的数值方法的有效性. 并基于动力学模态分解方法对空化流场的速度场进行模态分解,分析了各个模态的流场特征. 结果表明,第一阶模态对应频率为 0,代表平均流场;第二阶模态对应频率约为空泡脱落频率,揭示了空泡在水翼前缘周期性地生长与脱落,第三阶模态对应频率约为第二阶模态的 2 倍,揭示了两个大尺度旋涡在水翼后方存在融合行为. 第四阶模态对应频率约为第二阶模态的 3 倍,具有更高的频率,表征流场中存在一些小尺度旋涡的融合行为. 最后对不同空化数下的空化流场进行了模态分解分析,发现脱落空泡的旋涡结构随着空化数的减小而增大,第二阶模态频率随着空化数的减小而减小.     

Time resolved PIV and flow visualization of 3D sheet cavitation

Time-resolved PIV was applied to study fully developed sheet cavitation on a hydrofoil with a spanwise varying angle of attack. The hydrofoil was designed to have a three-dimensional cavitation pattern closely related to propeller cavitation, studied for its adverse effects as vibration, noise, and erosion production. For the PIV measurements, fluorescent tracer particles were applied in combination with an optical filter, in order to remove the reflections of the laser lightsheet by the cavitation. An adaptive mask was developed to find the interface between the vapor and liquid phase. The velocity at the interface of the cavity was found to be very close to the velocity predicted by a simple streamline model. For a visualization of the global flow dynamics, the laser beam was expanded and used to illuminate the entire hydrofoil and cavitation structure. The time-resolved recordings reveal the growth of the attached cavity and the cloud shedding. Our investigation proves the viability of accurate PIV measurements around developed sheet cavitation. The presented results will further be made available as a benchmark for the validation of numerical simulations of this complicated flow.     

基于动力学模态分解法的绕水翼非定常空化流场演化分析

空化是船舶和水下航行体推进器中经常发生的一种特殊流动现象,它具有强烈的非定常性,空化的发生往往会影响推进器的水动力性能和效率. 为探究绕水翼非定常空化流场结构,本文基于 Schnerr-Sauer 空化模型和 SST $k$-$\omega $ 湍流模型,开展绕二维水翼非定常空化流动数值预报与流场结构分析. 通过将数值预报的空泡形态演变和压力数据与试验结果对比,验证了建立的数值方法的有效性. 并基于动力学模态分解方法对空化流场的速度场进行模态分解,分析了各个模态的流场特征. 结果表明,第一阶模态对应频率为 0,代表平均流场;第二阶模态对应频率约为空泡脱落频率,揭示了空泡在水翼前缘周期性地生长与脱落,第三阶模态对应频率约为第二阶模态的 2 倍,揭示了两个大尺度旋涡在水翼后方存在融合行为. 第四阶模态对应频率约为第二阶模态的 3 倍,具有更高的频率,表征流场中存在一些小尺度旋涡的融合行为. 最后对不同空化数下的空化流场进行了模态分解分析,发现脱落空泡的旋涡结构随着空化数的减小而增大,第二阶模态频率随着空化数的减小而减小.      

非定常空化流场结构的实验研究

为深入研究非定常空化流场结构,本文用实验方法研究了绕Clark-Y型水翼的非定常空化流动现象.实验在空化水洞中进行,采用高速摄像技术观测了云状空化阶段的非定常空穴形态,并应用粒子成像测速系统(PIV)对绕水翼空化流场的速度场和涡量场等流动特性进行了同步的实验分析.研究表明:空化现象对流场结构有着重要的影响,在无空化和空...     

以超声振动原理模拟滑动轴承的气蚀磨损

采用超声振动原理建立了滑动轴承气蚀磨损模拟试验装置,并用该试验装置研究了滑动轴承气蚀磨损现象,通过分析比较大量实验研究结果,发现了一些新的气蚀现象,并对滑动轴承气蚀机理进行了分析。     

Experimental investigation of added mass effects on a hydrofoil under cavitation conditions

The influence of leading edge sheet cavitation and supercavitation on the added mass effects experienced by a 2-D NACA0009 truncated hydrofoil has been experimentally investigated in a hydrodynamic tunnel. A non-intrusive excitation and measuring system based on piezoelectric patches mounted on the hydrofoil surface was used to determine the natural frequencies of the fluid–structure system. The appropriate hydrodynamic conditions were selected to generate a range of stable partial cavities of various sizes and also to minimize the effects of other sources of flow noise and vibrations. The main tests were performed for different sigma values under a constant flow velocity of 14 m/s and for incident angles of both 1° and 2°. Additionally, a series of experiments in which the hydrofoil was submerged in air, partially and completely submerged in still water and without cavitation at 7 and 14 m/s were also performed. The maximum added mass effect occurs with still water. When cavitation appears, the added mass decreases because the cavity length is increased, and the added mass is minimum for supercavitation. A linear correlation is found between the added mass coefficients and the entrained mass that accounts for the mean density of the cavity, its dimensions and its location relative to the specific mode shape deformation.     

主动射流控制水翼空化的数值模拟与分析

为了改善高速流动工况下水翼吸力面上流场的空化特性,提出了水翼表面主动射流对绕水翼周围流动加以控制的方法.基于密度分域滤波的FBDCM混合湍流模型联合Zwart-Gerber-Belamri空化模型,分析了来流空化数为0.83,来流攻角为8°,射流位置距水翼前缘为x=0.19c时,主动射流对于水翼吸力面上流动的空化特性和水动力特性影响.对回射流的强度进行了量化分析,以探究回射流与流场空化特性的关系.数值分析结果表明,在射流水翼吸力面上的时均空泡体积为原始水翼的1/15,使得流场内空化流动由云空化状态转变为较为稳定的片空化状态,显著地削弱了云空化的发展.此外,射流极大地改善了水翼的水动力性能,使得水翼的升阻比较原始水翼提高了22.9%,空泡的脱落频率减少了26.2%,空泡脱落所引起的振幅减小了9.1%.射流大幅降低了水翼吸力面上低压区面积,水翼吸力面上流体的逆向压力减小,回射流强度降低;同时,射流使水翼吸力面上的边界层减薄,增强了流动的抗逆压梯度能力,一定程度上阻挡了回射流向水翼前缘的流动,这也从机理上分析了主动射流抑制空化的原因.      

The influence of surface roughness on cloud cavitation flow around hydrofoils

The aim of this study is to investigate experimentally the effect of surface roughness on cloud cavitation around Clark-Y hydrofoils. High-speed video and particle image velocimetry(PIV) were used to obtain cavitation patterns images(Prog. Aerosp. Sci. 37: 551–581, 2001), as well as velocity and vorticity fields. Results are presented for cloud cavitating conditions around a Clark-Y hydrofoil fixed at angle of attack of α = 8?for moderate Reynolds number of Re = 5.6 × 10~5. The results show that roughness had a great influence on the pattern, velocity and vorticity distribution of cloud cavitation. For cavitating flow around a smooth hydrofoil(A) and a rough hydrofoil(B), cloud cavitation occurred in the form of finger-like cavities and attached subulate cavities, respectively. The period of cloud cavitation around hydrofoil A was shorter than for hydrofoil B.Surface roughness had a great influence on the process of cloud cavitation. The development of cloud cavitation around hydrofoil A consisted of two stages:(1) Attached cavities developed along the surface to the trailing edge;(2) A reentrant jet developed, resulting in shedding and collapse of cluster bubbles or vortex structure. Meanwhile, its development for hydrofoil B included three stages:(1) Attached cavities developed along the surface to the trailing edge, with accumulation and rotation of bubbles at the trailing edge of the hydrofoil affecting the flow field;(2) Development of a reentrant jet resulted in the first shedding of cavities. Interaction and movement of flows from the pressure side and suction side brought liquid water from the pressure side to the suction side of the hydrofoil, finally forming a reentrant jet. The jet kept moving along the surface to the leading edge of the hydrofoil, resulting in large-scale shedding of cloud bubbles. Several vortices appeared and dissipated during the process;(3) Cavities grew and shed again.     

考虑空化影响的水润滑橡胶合金轴承冲蚀磨损研究

用计算流体力学方法,数值模拟水润滑轴承空化-冲蚀交互作用时气液固三相流场的动力学特性(压力场、速度场、气含率分布),然后在水润滑轴承摩擦磨损装置上,进行试件磨损试验,并观察试件表面形貌.结果表明:考虑空化影响后,水润滑轴承整个流场压力分布更接近实际;数值模拟所得的流场压力、速度、气含率最大值,均出现在发生空化的位置附近,其余位置基本不变,说明交互磨损比单一磨损严重.观察试件表面磨痕,存在短程犁沟、空蚀针孔、麻点状气蚀坑和蚀坑,磨痕呈现规律性,磨痕与轴转速的方向基本一致.试验结果和数值计算吻合较好,证明了理论分析的正确与合理.上述仿真与试验初步探讨了水润滑轴承空化与冲蚀交互作用的磨损机理与影响因素.     

Multiphase fluid dynamics and transport processes of low capillary number cavitating flows

To better understand the multiphase fluid dynamics and associated transport processes of cavitating flows at the capillary number of 0.74 and 0.54, and to validate the numerical results, a combined computational and experimental investigation of flows around a hydrofoil is studied based on flow visualizations and time-resolved interface movement. The computational model is based on a modified RNG k-ε model as turbulence closure, along with a vapor-liquid mass transfer model for treating the cavitation process. Overall, favorable agreement between the numerical and experimental results is observed. It is shown that the cavi- tation structure depends on the interaction of the water-vapor mixture and the vapor among the whole cavitation stage, the interface between the vapor and the two-phase mixture exhibits substantial unsteadiness. And, the adverse motion of the interface relates to pressure and velocity fluctuations inside the cavity. In particular, the velocity in the vapor region is lower than that in the two-phase region.     

Une technique de suivi d'interfaces pour la modélisation de la cavitation par poche

A tracking method is presented for the modeling of partial and supercavitation. The velocity and pressure fields in the cavitating flow are computed by a Navier–Stokes solver using a pseudo-compressibility method. The cavity flow is computed from the velocity field by a tracking method based on a volume of fluid technique (VOF). The method is illustrated by several computations, two cases of partial cavitation on a hydrofoil and a case of a cavitating body emerging at a free surface.     

Modelling the erosion efficiency of cavitation cleaning jets

A model is presented of the erosion damage and cleaning efficiency of cavitation and conventional cleaning jets. The effect of supply pressure, jet velocity, standoff distance, nozzle diameter and type of fluid are considered and assessed in relation to experimental findings.     

翼型空化绕流数值研究

空化是发生在流体机械上的复杂过程，理论研究遇到很大困难。本文引入合适的空化数值模型，将空腔界面近似为自由面，用界面构造精度较高的流体体积方法求解空腔位置，通过直接求解原始变量的NavuerStokes方程，数值模拟了无界域中空化在翼型上发生、发展和脱落的周期过程；并分析了空化产生对翼型表面的压力分布、翼型收到的阻力和升力的影响。结果表明，空化出现在翼型上表面；由于空化的产生，翼型表面压力分布不稳定，导致升力、阻力和流场压力出现波动，这是实际中产生噪声和损失的主要原因。     

Investigation of a re-entrant jet reflection at an inclined cavity closure line

Cavitation on two-dimensional hydrofoils with swept leading edges always displays some 3-dimensional effects. It is well known that the cavity closure on such hydrofoil is not perpendicular to the channel walls, but is curved in a distinctive pattern. The cavitation pocket is longer in the region where the hydrofoil is the shortest. Also the dynamics of cavitation is very distinctive. In the region where the hydrofoil is the longest attached and steady cavitation with no cloud separation exists. On the other side, where the hydrofoil is the shortest, cavitation cloud separations occur.Different explanations for this pattern were proposed in the past but they have not jet been clearly confirmed neither experimentally nor by numerical simulation.In the present paper a clear explanation supported by the numerical simulation and also by experimental measurements, is given.     

绕水翼超空化流动形态与速度分布

为揭示超空化流场结构特性,利用高速全流场显示技术,观察了绕hydronautics水翼的超空化流动形态,并利用数字粒子图像测速仪(DPIV)测量了其速度分布. 在测量空穴内部流速分布时,采用空化流场中的空化泡作为示踪粒子来显示流动结构. 结果表明:随着空化数的降低,超空化流动显现出了明显的阶段特征,其中水汽混合相和汽相的分布决定了空化区域的形态与流速分布;空化区和主流区的汽液交界面处存在着较大的速度梯度;低速分布区域随着空化数的降低由水翼吸力面中后部向水翼下游移动;在空化区域内部,水汽混合区的速度相对较低,而汽相区则与主流区有着相近的速度分布.关键词超空化水翼、DPIV、高速摄像、空化形态、流速分布      

Interfacial dynamics‐based modelling of turbulent cavitating flows,Part‐2: Time‐dependent computations

The interfacial dynamics‐based cavitation model, developed in Part‐1, is further employed for unsteady flow computations. The pressure‐based operator‐splitting algorithm (PISO) is extended to handle the time‐dependent cavitating flows with particular focus on the coupling of the cavitation and turbulence models, and the large density ratio associated with cavitation. Furthermore, the compressibility effect is important for unsteady cavitating flows because in a water–vapour mixture, depending on the composition, the speed of sound inside the cavity can vary by an order of magnitude. The implications of the issue of the speed of the sound are assessed with alternative modelling approaches. Depending on the geometric confinement of the nozzle, compressibility model and cavitation numbers, either auto‐oscillation or quasi‐steady behaviour is observed. The adverse pressure gradient in the closure region is stronger at the maximum cavity size. One can also observe that the mass transfer process contributes to the cavitation dynamics. Compared to the steady flow computations, the velocity and vapour volume fraction distributions within the cavity are noticeably improved with time‐dependent computations. Copyright © 2004 John Wiley & Sons, Ltd.     

绕栅中水翼空化流动的数值和实验研究

采用数值计算和实验研究的方法研究了绕水翼和栅中水翼的非定常空化流动. 实验采用高速录像技术分别观察了绕水翼和栅中水翼云状空化形态随时间的变化, 测量了升阻力, 并对测量数据进行了频率分析. 计算时空化模型选用了能比较准确描述旋涡空化非定常特性的Kubota模型, 湍流模型采用能准确捕捉流场非定常特性的FBM模型. 计算模型的可靠性用实验结果进行验证. 结果表明, 计算与实验的结果基本一致, 相比绕单个水翼的空化流动, 绕栅中水翼的空穴厚度比较薄, 翼型近壁处的逆压梯度较小, 反向射流的速度较小, 且水汽混合区速度梯度较小, 空穴的脱落周期变长, 平均升阻力系数较小