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

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

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

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

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

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

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

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

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

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

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成;(2) U型空泡团头部溃灭;(3) U型空泡团腿部溃灭.在U型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象,抑制了下一周期的空化发展.     

绕三维水翼云状空化现象的实验研究

为了研究云状空化阶段空穴发展和脱落的机理,采用实验的方法对绕三维水翼云状空化流动进行了研究.实验在高速水洞中进行,采用高速摄像技术研究了不同空化阶段的空穴形态,并测量了翼型所受的升阻力,并对上述数据进行了频谱分析.结果发现:在云状空化阶段,观测到空穴的产生-发展-脱落-溃灭的准周期性变化;并捕捉到空泡脱落时附着在翼型前...     

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成; (2) U型空泡团头部溃灭; (3) U型空泡团腿部溃灭.在U 型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象, 抑制了下一周期的空化发展.      

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

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

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

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

Velocity field measurements of cavitating flows

A particle Image Velocimetry (PIV) system has been developed to study the microfluid mechanics of cavitating flows. Planar PIV was used to examine the non-cavitating flow in the thin boundary layer near a hydrofoil surface for the cases of a naturally developing boundary layer and a boundary layer stimulated to turbulence by roughness near the foil leading edge. PIV was also used to examine the flow near the surface of individual cavitation bubbles and incipient attached cavitation. A system was devised to create a single nucleus in the flow upstream of a hydrofoil, and planar PIV was used to study the flow around the resulting traveling cavitation bubble. Velocity vectors were determined close to the solid surfaces and the gas/liquid interfaces of the bubbles. Seeding of the flow with particles did not result in the addition of active cavitation nuclei.     

LDA在减压箱空化特性试验研究中的应用

本文用新型的光纤激光测速仪(Fibrc-Optic LDA)对减压箱内闸墩附近紊动流场的空化特性进行了测试研究.结果表明:只有采用断面上各点流速的统计平均值作为参考断面的平均流速才能保证足够的精度;速度效应对空化初生的影响可不考虑,但紊流强度的变化对其影响却不容忽视。文中将不同尺寸闸墩的流速系数用统一经验公式表达,发现流速系数的变化与紊流强度密切相关,由此得出的紊流强度随流速而变化的曲线与实测值符合较好。本文还对 LDA 在类似减压箱的流场应用提出了几点体会和建议。     

Flow regime transitions due to cavitation in the flow through an orifice

This paper presents both experimental and theoretical aspects of the flow regime transitions caused by cavitation when water is passing through an orifice. Cavitation inception marks the transition from single-phase to two-phase bubbly flow; choked cavitation marks the transition from two-phase bubbly flow to two-phase annular jet flow.

It has been found that the inception of cavitation does not necessarily require that the minimum static pressure at the vena contracta downstream of the orifice, be equal to the vapour pressure liquid. In fact, it is well above the vapour pressure at the point of inception. The cavitation number [σ = (P₃ − P_v)/(0.5 pV²); here P₃ is the downstream pressure, P_v is the vapour pressure of the liquid, ρ is the density of the liquid and V is the average liquid velocity at the orifice] at inception is independent of the liquid velocity but strongly dependent on the size of the geometry. Choked cavitation occurs when this minimum pressure approaches the vapour pressure. The cavitation number at the choked condition is a function of the ratio of the orifice diameter (d) to the pipe diameter (D) only. When super cavitation occurs, the dimensionless jet length [L/(D - d); where L is the dimensional length of the jet] can be correlated by using the cavitation number. The vaporization rate of the surface of the liquid jet in super cavitation has been evaluated based on the experiments.

Experiments have also been conducted in which air was deliberately introduced at the vena contracta to simulate the flow regime transition at choked cavitation. Correlations have been obtained to calculate the critical air flow rate required to cause the flow regime transition. By drawing an analogy with choked cavitation, where the air flow rate required to cause the transition is zero, the vapour and released gas flow rate can be predicted.     

Coherent and turbulent process analysis of the effects of a longitudinal vortex on boundary layer detachment on a NACA0015 foil

In this paper, the influence of a single tip vortex on boundary layer detachment is studied. This study offers a preliminary approach in order to better understand the interaction between a propeller hub vortex and the rudder installed in its wake. This configuration belongs to the field of marine propulsion and encompasses such specific problem as cavitation inception, modification of propulsive performances and induced vibrations. To better understand the complex mechanisms due to propeller–rudder interactions it was decided to emphasize configurations where the hub vortex is generated by an elliptical 3-D foil and is located upstream of a 2-D NACA0015 foil at high incidences for a Reynolds number of 5×10⁵. The physical mechanisms were studied using Time Resolved Stereoscopic Particle Image Velocimetry (TR-SPIV) techniques. Particular attention was paid to the detachment at 25° incidence and a detailed cartography of the mean and turbulent properties of the wake is presented. Proper Orthogonal Decomposition (POD) analysis was applied in order to highlight the unsteady nature of the flow using phase averaging based on the first POD coefficients to characterize the turbulent and coherent process in the near wake of the rudder.     

瞬态空化泡演变过程的数值模拟

采用边界积分方程方法，对无粘流体中三个空化泡以及自由面附近二个空化泡相互作用的演变过程进行了数值模拟。计算中边界用二阶有限元离散，影响系数矩阵非对角线元素用六点高斯数值积分方法计算，对第一类、第二类完全椭圆积分用高次多项式近似，对计算系数矩阵对角线元素中遇到的奇异积分进行了特殊处理。结果表明，在不同的给定参数下，空化泡的溃灭形态各异，柱状射流和环形射流都有可能发生，使空化泡演变成双泡或环形泡。     

A numerical scheme for Euler–Lagrange simulation of bubbly flows in complex systems

An Eulerian–Lagrangian approach is developed for the simulation of turbulent bubbly flows in complex systems. The liquid phase is treated as a continuum and the Navier–Stokes equations are solved in an unstructured grid, finite volume framework for turbulent flows. The dynamics of the disperse phase is modeled in a Lagrangian frame and includes models for the motion of each individual bubble, bubble size variations due to the local pressure changes, and interactions among the bubbles and with boundaries. The bubble growth/collapse is modeled by the Rayleigh–Plesset (RP) equation. Three modeling approaches are considered: (a) one‐way coupling, where the influence of the bubble on the fluid flow is neglected, (b) two‐way coupling, where the momentum‐exchange between the fluid and the bubbles is modeled, and (c) volumetric coupling, where the volumetric displacement of the fluid by the bubble motion and the momentum‐exchange are modeled. A novel adaptive time‐stepping scheme based on stability‐analysis of the non‐linear bubble dynamics equations is developed. The numerical approach is verified for various single bubble test cases to show second‐order accuracy. Interactions of multiple bubbles with vortical flows are simulated to study the effectiveness of the volumetric coupling approach in predicting the flow features observed experimentally. Finally, the numerical approach is used to perform a large‐eddy simulation in two configurations: (i) flow over a cavity to predict small‐scale cavitation and inception and (ii) a rising dense bubble plume in a stationary water column. The results show good predictive capability of the numerical algorithm in capturing complex flow features. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Statistical structure of the velocity field in cavitating flow around a 2D hydrofoil

Transformation of flow turbulence structure with cavitation occurrence, determination of the flow conditions favorable for nucleation of cavitation bubbles, influence of the statistical structure of turbulence on this process and the inverse effect of cavitation on the flow dynamics are challenging problems in modern fluid mechanics. The paper reports on the results of statistical processing of the velocity fields measured by a PIV technique in cavitating flow over a 2D symmetric hydrofoil for four flow conditions, starting from a cavitation-free regime and finishing by unsteady cloud cavitation. We analyze basic information on the statistical structure of velocity fluctuations in the form of histograms and Q-Q diagrams along with profiles of the mean velocity and turbulent kinetic energy. The research reveals that the flow turbulence pattern and distributions of turbulent fluctuations change significantly with the cavitation development. Under unsteady cloud cavitation conditions, the probability density function of the fluctuating velocity has a two-mode distribution, which indicates switching of two alternating flow conditions in a region above the hydrofoil aft part due to periodic passing of cavitation clouds. Behaviors of the mean and most probable velocities unexpectedly appear to be different with a monotonous increase of the incoming flow velocity. This finding must be caused by modification of the skewness coefficient of the fluctuating velocity.     

Influence of bubble size on micro-bubble drag reduction

Micro-bubble drag reduction experiments were conducted in a turbulent water channel flow. Compressed nitrogen was used to force flow through a slot injector located in the plate beneath the boundary layer of the tunnel test section. Gas and bubbly mixtures were injected into a turbulent boundary layer (TBL), and the resulting friction drag was measured downstream of the injector. Injection into tap water, a surfactant solution (Triton X-100, 20 ppm), and a salt-water solution (35 ppt) yielded bubbles of average diameter 476, 322 and 254 μm, respectively. In addition, lipid stabilized gas bubbles (44 μm) were injected into the boundary layer. Thus, bubbles with d ⁺ values of 200 to 18 were injected. The results indicate that the measured drag reduction by micro-bubbles in a TBL is related strongly to the injected gas volumetric flow rate and the static pressure in the boundary layer, but is essentially independent of the size of the micro-bubbles over the size range tested.     

湍流边界层拟序结构的实验研究

20世纪60年代后, 先后从流动显示发现了快慢斑、猝发、上升流、下扫流和多种涡结构等湍流边界层的拟序结构. 它们对湍流边界层的摩阻、传热传质和湍动能的产生等特性有重要影响. 涡结构是上述拟序结构的核心, 它影响其它拟序结构的发展和演变. 发卡涡通常被认为是基本涡结构. 发卡涡等涡结构的再生, 是湍流边界层拟序结构能够自持续的必要的因素.壁面低速流上升产生猝发, 是湍流边界层湍能的主要来源; 条件采样是测量猝发频率和其它拟序结构出现频率的重要手段. 流动显示对湍流边界层拟序结构作了大量定性观察, 有许多减阻和增加传热率等应用性研究在此基础上发展起来. 80年代后, 出现了测量湍流边界层的瞬时流速矢量场的多热线法和PIV技术, 三维PIV技术可望将来为湍流边界层的实验研究带来重大进展. 本文评述了流动显示法、多热线法和PIV技术的优点和不足之处, 以及它们在对湍流边界层拟序结构的研究中的贡献.