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

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

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

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

绕水翼空化流动及振动特性的实验研究

空化是发生在水力机械内部的一种水动力现象,其发展具有显著的非定常特性.空化流动中空穴的脱落以及溃灭会诱发结构振动,对水力机械的效率、噪声、安全性等造成影响. 研究空化流动中结构的振动特性具有重要的工程意义. 采用实验的方法研究了绕NACA66 水翼空化流动的空穴形态和水翼振动特性. 实验在一闭式空化水洞中进行. 采用高速摄像技术观测不同空化阶段的空穴形态,应用多普勒激光测振仪测量水翼的振动速度,并通过一套同步系统实现了高速相机和多普勒激光测振仪的同步触发和测量. 采用小波分析方法对不同空化阶段下的空穴形态和水翼振动数据在时域和频域中的特性进行了分析.对云状空化阶段的同步测量结果进行了研究,分析了振动与空穴发展过程的联系. 结果表明,随着空化数的降低,流场经历了无空化、初生空化、片状空化和云状空化4个阶段,水翼的振动强度呈逐渐增大趋势. 在片状空化和云状空化阶段,空穴脱落导致水翼振动,诱发的振动频率与空穴脱落频率相同. 对于云状空化,在附着型空穴生长阶段水翼发生高频小幅度振动,在空穴脉动和断裂脱落期间水翼表现为低频大幅振动.      

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

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

采用分段湍流模式研究绕水翼的空化流动

结合空化流动特点,建立了一个包含空间尺度信息的分段湍流模式。计算中,应用基于质量传输空化模型,分别采用三种湍流模型计算了绕Clark-y型水翼云状空化流动,得到了随时间变化的空泡形态以及升、阻力等流场和动力特性。通过与实验结果的对比,表明这三种湍流模型均能捕捉云状空化区域的空泡形态和空泡脱落的非定常细节。分段湍流模式能够更好地调整流场内的湍流黏性,更精确地预测空穴长度和空穴尾部水汽分布,与实验结果吻合较好。     

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

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

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

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

云状空化非定常脱落机理的数值与实验研究

结合数值计算和实验技术研究了云状空化的非定常脱落机理. 实验采用高速录像技术观察了绕Clark-y型水翼云状空化形态随时间的变化. 数值计算采用了汽-液两相的当地均相流模型,湍流封闭采用了一种修正的RNG k-\varepsilon方程,利用商业软件的二次开发技术,引入了一种与空化区域水汽相密度相关的系数,对湍流模型进行了修正. 计算和实验结果均表明:云状空化尾部存在着准周期性的涡状空化团的脱落;局部压强的增加是引起空穴断裂进而脱落的直接原因;压强升高是由于近壁处的反向射流引起的;空穴尾部的大规模的旋涡运动引起的近壁处的逆压梯度是引起反向射流的主要原因.      

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

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

非定常空化流动涡旋运动及其流体动力特性

采用大涡模拟方法对绕水翼云状空化的水动力特性和非定常流场结构进行研究. 基于实验结果对数值方法进行验证,分析空化与流场内部涡旋结构之间的相互作用以及对水翼动力特性的影响. 研究结果表明:大涡模拟方法可以准确模拟绕水翼流动的非定常过程. 在无空化条件下,升阻力系数存在斯特劳哈数St = 0.85 的主频波动,这是由水翼尾部涡旋结构的发展脱落引起的;在云状空化条件下,升阻力系数存在St = 0.34 的高能量密度低频波动,这是由大规模云状空泡团的发展和脱落引起的;云状空化阶段的升阻力系数在St = 0.5~1.5 的范围内都存在较高的波动,这是由于空化现象对水翼尾缘涡旋结构的发展和脱落产生影响,在不同发展阶段,空化现象不同程度地降低尾缘涡旋结构脱落频率.      

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.     

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.     

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

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

A numerical study of the influence of the Basset force on the statistics of LDV velocity data sampled in a flow region with a large spatial velocity gradient

 The motion of small particles, such as those typically used as seeding particles for tracer particle flow velocity measurement techniques, is studied numerically for a flow region with a large spatial velocity gradient. The influence of the Basset history integral on the statistics of results of particle motion calculations which are based on multi-disperse particle size distributions is investigated. The biasing of the measured velocity data, with regard to the actual flow velocity, which results as a consequence of such particle size distributions is discussed. It is found that the net effect of the Basset integral on the calculations is indeed to reduce the maximum RMS deviation associated with the multi-disperse distribution and that the relative reduction increases with a decreasing particle density. The main result of this study is, thus, that it is desirable to use light tracer particles not only because they more readily adjust to a changed flow velocity but in particular also because they tend to contribute less to the overall RMS deviation of velocity data sampled in a region with a large spatial gradient of the flow velocity. Received: 9 August 1996/Accepted: 25 January 1997     

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.     

Flow around a cylinder at the beginning of the critical region

Results of experiments with a turbulent flow around a transversely aligned circular cylinder located at identical distances from the walls of a rectangular channel are reported. Data on averaged velocity fields around the cylinder are obtained by means of particle image velocimetry (PIV). Based on these fields, the near wake behind the cylinder is studied, and the kinematic characteristics for flow regimes with and without cavitation are compared. Based on the vector fields of averaged velocity, the angles of separation of the boundary layer from the cylinder surface in the considered flow regimes are determined. The drag coefficients of the cylinder for different flow regimes are calculated. It is demonstrated that the vortex region behind the cylinder and the drag coefficient of the cylinder increase in the case with cavitation. It is also shown that vortex shedding from the cylinder may be irregular, despite the fact that this process is quasi-periodic for most of the time.     

Flow around a sphere with material cross flow at low reynolds numbers

An analytical solution is carried out for the problem of the flow around a sphere with material cross flow at Reynolds numbers less than 1 and a blowing velocity less than the free stream velocity. The method of asymptotic expansions of Pearson and Proudman is used for the solution. Expressions are obtained for the distribution of the current and velocity component functions as well as for the aerodynamic drag coefficient of the sphere. It is shown that blowing diminishes the sphere drag, where its influence will increase as the Reynolds number grows.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 103–109, May–June, 1972.     

Combined experimental observation and numerical simulation of the cloud cavitation with U-type flow structures on hydrofoils

Sheet/cloud cavitation is an important topic that is a very common type of cavitation in turbo-machinery and marine propeller. Up to now we still have limited understanding of the cavitation shedding dynamics and cloud cavity formation and development. The present study used experimental and numerical studies to gain a better understanding of the complex physics involved in this problem. A series of experimental observations around hydrofoils are carried out in the cavitation tunnel of the China Ship Scientific Research Center (CSSRC) to illustrate the spatial–temporal evolution of the cloud cavity in detail. The results demonstrate that U-type flow structures are common in cloud cavities and can be divided into three stages and the closure line in a sheet cavity often has a convex–concave profile. Reentrant flows occur in the convex region with the jet direction normal to the contour edge so the shedding is mainly caused by the converging reentrant flows. Further analysis demonstrated that there was a striking difference with the cavity growth suppressed substantially in the twisted hydrofoil case if compared with straight hydrofoil and the effect of side entrant jets might make the cavity more uniform across the span. Numerical simulations were used to simulate the formation and development of the cloud cavity. The results show that the strong adverse pressure gradient in the stagnation region at the downstream end of the attached cavity forces the re-entrant flows into the vapor structure with a radially-diverging re-entrant jet and a pair of side-entrant jets, which causes the cavity shedding. Further analyzes of the local flow fields show that the interactions between the circulating flow and the shedding vapor cloud may be the main reason for the formation of the U-type cloud cavity structures.     

The boundary layer on a rapidly rotating cylinder

The study considers plane steady flow of an incompressible fluid around a circular cylinder rotating in a homogeneous free stream. On the basis of an asymptotic analysis of the Navier-Stokes equations for high Reynolds numbers, it is shown that at a certain value of the angular velocity of the cylinder an interaction arises between the flow in the boundary layer and the external potential flow. A solution is obtained numerically which describes the flow in the region of interaction.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 36–45, September–November, 1987.     

Flow past a plate lattice with developed cavitation

Problems of streamline cavitation flow past a lattice were studied in [1–8] using the Kirchhoff scheme. In this scheme the magnitude of the velocity at the free surface is equal to the stream velocity behind the lattice, and the cavitation number is zero (for a lattice the relative velocity and the cavitation number are defined from the stream velocity behind the lattice). In [4, 7] a solution is given of the problem of flow past a lattice using a scheme with an Efros-Gilbargreturn streamline, which permits considering arbitrary cavitation numbers; however, a unique solution is not given. Some other streamline schemes are mentioned in [8].In the following we consider the cavitational flow of an ideal incompressible inviscid and weightless fluid past an infinite lattice of flat plates, using the streamline wake model previously utilized by Wu [9] in studying cavitational flow past an isolated obstacle. In accordance with this model, the streamlines which separate from the body and bound the cavity behind it pass into two curvilinear infinitely long walls, along which the pressure increases and approaches the pressure in the undisturbed stream.It is further assumed that in the hodograph plane there corresponds to the curvilinear walls a cut along some line and that the complex potential takes the same values at points lying on opposite sides of the cut. In particular, at the points of contact of the streamlines with the curvilinear walls the complex potential is the same. In the Wu scheme the latter condition leads to vanishing of the velocity circulation along the contour CABC₁ (Fig. 1).In conclusion the author wishes to thank N. V. Yurtaeva for the accurately performed numerical work.