可压缩空化流动空穴演化及压力脉动特性实验研究

空化流动具有高度的压缩性,空化流动非定常特性及其流体动力与压缩性密切相关.为研究可压缩空化流动空泡脱落的回射流和激波机制下周期性空穴结构演化及其诱导流体动力特性,本文采用多场同步测试方法对典型云状空化流动进行了实验研究,获得了文丘里管扩张段内部云状空化空穴形态演化及其诱导同步壁面压力脉动信号.并基于数字图像处理技术,对附着型片状空穴和脱落型云状空穴结构演化进行了精细化定量分析.结果表明:可压缩空化流动回射流机制下,空穴演化呈现附着型空穴生长-回射流产生及发展-附着型空穴失稳断裂及大尺度空泡云团产生脱落的非定常过程,激波机制下空穴演化具有附着型空穴生长-激波产生及传播-附着型空穴失稳断裂及大尺度空泡云团脱落的非定常特征,激波传播时间占空穴脱落周期小于回射流推进.激波与空穴相互作用导致空穴内部含气率瞬间大范围大幅度下降,诱导复杂流体动力.激波传播过程中,空泡内部压力脉动大幅增加,激波前缘诱导压力脉冲,而回射流推进过程中,壁面压力脉动相对平稳,回射流头部存在小幅增加.不同机制下空穴结构存在显著差异,具有不同的相间质量传输过程.     

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

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

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

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

绕水翼空化流动多尺度数值研究

空化的多尺度效应是一种涉及连续介质尺度、微尺度空化泡以及不同尺度间相互转化的复杂水动力学现象, 跨尺度模型的构建是解析该多尺度现象的关键. 本文基于欧拉-拉格朗日联合算法, 通过界面捕捉法求解欧拉体系下大尺度空穴演化, 通过拉格朗日体系下离散空泡模型求解亚网格尺度离散空泡的运动及生长溃灭. 同时, 通过判断空泡与网格尺度间的关系判定不同尺度空化泡的求解模型. 基于建立的多尺度算法对绕NACA66水翼空化流动进行模拟, 将数值结果与实验进行对比, 验证了数值计算方法的准确性. 研究结果表明, 离散空泡数量与空化发展阶段密切相关, 在附着型片状空穴生长阶段, 离散空泡数量波动较小, 离散空泡主要分布在气液交界面位置; 在回射流发展阶段, 离散空泡逐渐增加并分布在回射流扰动区; 在云状空穴溃灭阶段, 离散空泡数量增多且主要分布在气液掺混剧烈的空化云团溃灭区. 在各空化发展阶段, 离散空泡直径概率密度函数均符合伽玛分布. 空化湍流流场特性对拉格朗日空泡空间分布具有重要影响, 离散空泡主要分布在强湍脉动区、旋涡及回射流发展区域.      

收缩扩张管内液氮空化流动演化过程试验研究

本文基于低温空化试验平台研究了收缩扩张流道内液氮非定常空化流动的演化过程. 试验采用高时空分辨率的高速摄像机对77 K液氮在不同空化数σ下空穴结构的演变进行了精细化的分析和研究. 利用试验得到的空穴长度和面积等数据, 定量分析了液氮空化流动的非定常特性与时空演变规律. 研究结果表明: (1)在相似来流速度和温度条件下, 随着空化数的减小, 液氮空化流动呈现四种典型流型, 空穴长度在2.5 h以内为初生空化、空穴长度在2.5 h ~ 7.5 h之间为片状空化、空穴长度在7.5 h ~ 15 h之间为大尺度云状空化, 空穴长度超过15 h为双云状空化, 且在大尺度云状空化和双云状空化阶段均捕捉到了回射流现象; (2)液氮空化流动从初生空化到双云状空化, 脱落空穴的尺度逐渐增大, 空穴面积脉动的幅值和准周期均有所增加. 同时, 在大尺度云状空化与双云状空化阶段, 喉口处堵塞效应对空化流动的影响显著增强; (3)相比于初生空化, 片状空化、大尺度云状空化以及双云状空化中脱落空穴的移动距离依次增加了0.97倍、2.65倍与2.68倍, 溃灭时间依次增加了1.18倍、3.59倍与4.47倍, 但溃灭速度依次减小了0.10倍、0.20倍与0.30倍. 除此之外, 对于双云状空化阶段, 存在两种显著不同的脱落空穴演化过程.      

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

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

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

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

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

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

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

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

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

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

Passive control of cavitating flow around an axisymmetric projectile by using a trip bar

Quasi-periodical evolutions such as shedding and collapsing of unsteady cloud cavitating flow, induce strong pressure fluctuations, what may deteriorate maneuvering stability and corrode surfaces of underwater vehicles. This paper analyzed effects on cavitation stability of a trip bar arranged on high-speed underwater projectile. Small scale water tank experiment and large eddy simulation using the open source software Open FOAM were used, and the results agree well with each other. Results also indicate that trip bar can obstruct downstream re-entrant jet and pressure wave propagation caused by collapse, resulting in a relatively stable sheet cavity between trip bar and shoulder of projectiles.     

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.     

Numerical analysis of unsteady cavitating turbulent flow and shedding horse-shoe vortex structure around a twisted hydrofoil

Cavitating turbulent flow around hydrofoils was simulated using the Partially-Averaged Navier–Stokes (PANS) method and a mass transfer cavitation model with the maximum density ratio (ρ_l/ρ_v,clip) effect between the liquid and the vapor. The predicted cavity length and thickness of stable cavities as well as the pressure distribution along the suction surface of a NACA66(MOD) hydrofoil compare well with experimental data when using the actual maximum density ratio (ρ_l/ρ_v,clip = 43391) at room temperature. The unsteady cavitation patterns and their evolution around a Delft twisted hydrofoil were then simulated. The numerical results indicate that the cavity volume fluctuates dramatically as the cavitating flow develops with cavity growth, destabilization, and collapse. The predicted three dimensional cavity structures due to the variation of attack angle in the span-wise direction and the shedding cycle as well as its frequency agree fairly well with experimental observations. The distinct side-lobes of the attached cavity and the shedding U-shaped horse-shoe vortex are well captured. Furthermore, it is shown that the shedding horse-shoe vortex includes a primary U-shaped vapor cloud and two secondary U-shaped vapor clouds originating from the primary shedding at the cavity center and the secondary shedding at both cavity sides. The primary shedding is related to the collision of a radially-diverging re-entrant jet and the attached cavity surface, while the secondary shedding is due to the collision of side-entrant jets and the radially-diverging re-entrant jet. The local flow fields show that the interaction between the circulating flow and the shedding vapor cloud may be the main mechanism producing the cavitating horse-shoe vortex. Two side views described by iso-surfaces of the vapor volume fraction for a 10% vapor volume, and a non-dimensional Q-criterion equal to 200 are used to illustrate the formation, roll-up and transport of the shedding horse-shoe vortex. The predicted height of the shedding horse-shoe vortex increases as the vortex moves downstream. It is shown that the shape of the horse-shoe vortex for the non-dimensional Q-criterion is more complicated than that of the 10% vapor fraction iso-surface and is more consistent with the experiments. Further, though the time-averaged lift coefficient predicted by the PANS calculation is about 12% lower than the experimental value, it is better than other predictions based on RANS solvers.     

Scale effect on unsteady cloud cavitation

No experiment was conducted, yet, to investigate the scale effects on the dynamics of developed cavitating flow with periodical cloud shedding. The present study was motivated by the unclear results obtained from the experiments in a Venturi-type section that was scaled down 10 times for the purpose of measurements by ultra-fast X-ray imaging (Coutier-Delgosha et al. 2009). Cavitation in the original size scale section (Stutz and Reboud in Exp Fluids 23:191–198, 1997, Exp Fluids 29:545–552 2000) always displays unsteady cloud separation. However, when the geometry was scaled down, the cavitation became quasi steady although some oscillations still existed. To investigate this phenomenon more in detail, experiments were conducted in six geometrically similar Venturi test sections where either width or height or both were scaled. Various types of instabilities are obtained, from simple oscillations of the sheet cavity length to large vapor cloud shedding when the size of the test section is increased. It confirms that small scale has a significant influence on cavitation. Especially the height of the test section plays a major role in the dynamics of the re-entrant jet that drives the periodical shedding observed at large scale. Results suggest that the sheet cavity becomes stabile when the section is scaled down to a certain point because re-entrant jet cannot fully develop.     

Numerical investigation of periodic cavitation shedding in a Venturi

Unsteady partial cavitation is mainly formed by an attached cavity which presents periodic oscillations. Under certain conditions, instabilities are characterized by the formation of vapour clouds convected downstream the cavity and collapsing in higher pressure region. Two main mechanisms have been identified for the break-off cycles. The development of a liquid re-entrant jet is the most common type of instabilities, but more recently, the role of pressure waves created by the cloud collapses has been highlighted. This paper presents one-fluid compressible simulations of a self-sustained oscillating cavitation pocket developing along a Venturi geometry. The mass transfer between phases is driven by a void ratio transport equation model. The importance of traveling pressure waves in the physical mechanism is put in evidence. Moreover, the importance of considering a non-equilibrium state for the vapour phase is exhibited.     

Sheet Cavitation and Cloud Cavitation, Re-Entrant Jet and Three-Dimensionality

The unsteady behaviour of a sheet cavity on a hydrofoil under stationary conditions is important with respect to both noise generation and erosion. Observations of two-dimensional as well as three-dimensional cavities are performed. From the observation it is concluded that the formation of a re-entrant jet is of major importance to explain the unsteadiness. Numerical simulation of the two-dimensional cavitating flow is performed, showing good agreement with the observations.