基于Rayleigh-Plesset方程的空化模型改进与应用

为了提高Schnerr-Sauer模型预测空化流场的能力,提出了一种修正的球形空泡动力学模型,结合均相流假设建立了一种基于Rayleigh-Plesset方程的改进空化模型。联立修正后的滤波器湍流模型,分别采用改进后空化模型和Schnerr-Sauer模型对绕二维Clark-Y型水翼非定常空化流动进行计算,得到了不同空化数条件下翼型吸力面上空穴形态随时间变化规律及水动力系数特征。与已有文献的实验结果进行对比分析,结果表明:在准静态片状空化阶段,改进的空化模型计算得到的片状空穴长度基本上维持不变,在空穴尾部捕捉到明显的相间界面,与实验描述较为一致;在云状空化阶段,改进的空化模型模拟的空化区域范围更广,预测的瞬时升力系数与实验值的变化趋势吻合程度更高,且其对云空化周期性演变过程中一些流动细节的预测能力更好。     

基于Rayleigh-Plesset方程的改进空化模型及应用

为了提高Schnerr-Sauer模型预测空化流场的能力,提出了一种修正的球形空泡动力学模型,结合均相流假设建立了一种基于Rayleigh-Plesset方程的改进空化模型。联立修正后的滤波器湍流模型,分别采用改进后空化模型和Schnerr-Sauer模型对绕二维Clark-Y型水翼非定常空化流动进行计算,得到了不同空化数条件下翼型吸力面上空穴形态随时间变化规律及水动力系数特征。通过与已有文献的实验结果进行对比分析,结果表明:在准静态片状空化阶段,改进的空化模型计算得到的片状空穴长度基本上维持不变,在空穴尾部捕捉到明显的相间界面,与实验描述较一致;在云状空化阶段,改进的空化模型模拟的空化区域范围更广,预测的瞬时升力系数与实验值的变化趋势吻合程度更高,且其对云空化周期性演变过程中一些流动细节的预测能力更好。     

基于Rayleigh-Plesset方程的改进空化模型及应用

为了提高Schnerr-Sauer模型预测空化流场的能力,提出了一种修正的球形空泡动力学模型,结合均相流假设建立了一种基于Rayleigh-Plesset方程的改进空化模型。联立修正后的滤波器湍流模型,分别采用改进后空化模型和Schnerr-Sauer模型对绕二维Clark-Y型水翼非定常空化流动进行计算,得到了不同空化数条件下翼型吸力面上空穴形态随时间变化规律及水动力系数特征。通过与已有文献的实验结果进行对比分析,结果表明:在准静态片状空化阶段,改进的空化模型计算得到的片状空穴长度基本上维持不变,在空穴尾部捕捉到明显的相间界面,与实验描述较一致;在云状空化阶段,改进的空化模型模拟的空化区域范围更广,预测的瞬时升力系数与实验值的变化趋势吻合程度更高,且其对云空化周期性演变过程中一些流动细节的预测能力更好。     

基于Rayleigh-Plesset方程的改进空化模型及应用

为了提高Schnerr-Sauer模型预测空化流场的能力,提出了一种修正的球形空泡动力学模型,结合均相流假设建立了一种基于Rayleigh-Plesset方程的改进空化模型。联立修正后的滤波器湍流模型,分别采用改进后空化模型和Schnerr-Sauer模型对绕二维Clark-Y型水翼非定常空化流动进行计算,得到了不同空化数条件下翼型吸力面上空穴形态随时间变化规律及水动力系数特征。通过与已有文献的实验结果进行对比分析,结果表明:在准静态片状空化阶段,改进的空化模型计算得到的片状空穴长度基本上维持不变,在空穴尾部捕捉到明显的相间界面,与实验描述较一致;在云状空化阶段,改进的空化模型模拟的空化区域范围更广,预测的瞬时升力系数与实验值的变化趋势吻合程度更高,且其对云空化周期性演变过程中一些流动细节的预测能力更好。     

基于Rayleigh-Plesset方程的改进空化模型及应用

为了提高Schnerr-Sauer模型预测空化流场的能力,提出了一种修正的球形空泡动力学模型,结合均相流假设建立了一种基于Rayleigh-Plesset方程的改进空化模型。联立修正后的滤波器湍流模型,分别采用改进后空化模型和Schnerr-Sauer模型对绕二维Clark-Y型水翼非定常空化流动进行计算,得到了不同空化数条件下翼型吸力面上空穴形态随时间变化规律及水动力系数特征。通过与已有文献的实验结果进行对比分析,结果表明:在准静态片状空化阶段,改进的空化模型计算得到的片状空穴长度基本上维持不变,在空穴尾部捕捉到明显的相间界面,与实验描述较一致;在云状空化阶段,改进的空化模型模拟的空化区域范围更广,预测的瞬时升力系数与实验值的变化趋势吻合程度更高,且其对云空化周期性演变过程中一些流动细节的预测能力更好。     

水翼表面布置射流水孔抑制空化

为研究水翼表面布置射流水孔对空化的抑制效果,采用RNG k-ε湍流模型结合全空化模型对8°攻角NACA0066(MOD)水翼进行二维非定常空化流场计算,分析当射流比分别为0.1,0.2和0.3时的水翼表面空化流场结构及其演化过程的流动特性,得到了不同射流条件下水翼升阻力系数的变化趋势.计算结果表明:布置射流水孔可以有效阻挡回射流从翼型尾部向翼型头部的运动,使得空泡脱落现象明显减弱,从而抑制了片空化向云空化的发展,且随射流比增加,射流对回射流和云空化的抑制作用更加显著;射流的注入同时使得升力系数和阻力系数的波动幅度明显降低,有效减弱了升阻力的周期性振荡;虽然时均升力和时均阻力受射流的影响有下降趋势,但升阻比呈上升趋势,水翼升阻比由无射流孔时的6.11增加到射流比为0.3时的7.06,布置射流孔还起到了减阻的作用.     

绕水翼空泡脱落及其空化脉动特性

为改善水力机械抗空化性能,采用数值模拟与试验相结合的方法对NACA0015水翼非定常云空泡脱落机理进行研究,分析不同空化阶段下非定常空泡结构对翼型表面产生的空化脉动规律,探讨空化非定常过程中压力脉动产生的主要原因.结果表明:基于密度分域滤波器的湍流模型(FBDCM)能较好地模拟水翼表面空泡周期性脱落的非定常过程;在攻角...     

非定常空化流动研究现状与进展

针对非定常空化流动数值计算模型,总结了不同汽液相间传输模型对模拟空穴发展和空泡脱落过程的影响,发现不同空化模型对于非定常空化流场的不同区域的预测有着明显的差异.建立了一种基于混合密度分域的空化流动计算模型,合理预测了空泡脱落时刻翼型尾部的不稳定脉动区.通过考虑空化流动的湍流多尺度与多相局部可压缩特性,提出了一种基于混合密度分域的混合湍流模型,预测结果与试验结果吻合更好.综合回顾试验与数值计算方面的研究工作,围绕空化流动特性与动力学问题,主要从非定常流场空穴形态、流动结构、空化流体动力等几个方面对非定常空化流动机理及相应研究成果进行了较为详细地评述和讨论.最后,基于非定常空化研究的发展趋势和进一步研究可能面临的突出问题,对非定常空化流动研究的发展方向进行了探讨与展望.     

绕水翼加速流空化特性数值模拟

为了解翼型加速流动过程中空化特性,采用分离涡湍流模型(Detached eddy simulation,DES)和均相流空化模型对攻角为5.8°的NACA66(mod)水翼进行数值模拟,分析了空化数为0.99、对应雷诺数为8×105时绕二维水翼的非定常流动。通过模拟2种不同加速度(5 m/s2和2.5 m/s2)情况下空化演变过程和流场结构变化特征,得出了加速过程特有的变化规律:不同加速流下空泡都先在翼型的前缘产生,经过一段时间发展壮大,在翼型尾翼处分离,前缘处空泡减少,尾翼处空泡增加并向后衍生,直至破裂。空化产生于总加速时间的0.6倍左右,在总加速时间的1.12倍左右结束第一个周期。加速度越小,升力系数振荡范围越小,阻力系数振荡范围越大,空化发展过程越缓慢。     

考虑不可凝结气体的空化流模型及数值模拟

以均相流假设为基础建立了一种基于输运方程的空化流模型,该模型在质量传输方程中不仅考虑了蒸发和凝结的机理,而且考虑了不可凝结气体的影响.采用RNG k-ε湍流模型,并引入与混合密度相关的修正函数对湍流涡黏性系数进行修正.应用文中的空化流模型,对NA-CA66翼型进行了定常空化流动数值模拟,翼型吸力面的压力系数分布曲线与试验结果吻合很好.在此基础上,进一步研究了模型中不可凝结气体质量分数以及不同进口湍动能和湍流耗散率对空泡形成和发展的影响,确定了模型中不可凝结气体质量分数、进口湍动能和湍流耗散率的合理取值.应用空化流模型对非定常空化流动进行了数值模拟,数值模拟结果清晰地反映了翼型表面空化云的初生、成长、脱落和溃灭的全过程,并指出反向射流是引起空化云脱落的重要原因.非定常计算得到的斯特劳哈数与试验相吻合,进一步验证了该模型在空化流数值计算中的可靠性.     

三维ALE15翼型空化流动数值模拟

在均相流假设下,考虑流体压力和速度湍流脉动、不可凝结性气体的影响,采用完全空化模型计算空化流场的相变,引入密度函数对RNG k-ε湍流模型的湍流粘性系数进行修正,提出了一种空化流动的数值模型和计算方法。根据试验条件给定的参数,采用提出的数值模型和计算方法,数值模拟了空化数为2.3时ALE15翼型定常空化流动。计算得到的不同剖面速度分布与试验数据吻合较好,验证了该数值模型和计算方法的一致性。不同剖面上,远离翼型表面的速度与主流区速度接近,沿着流动方向,远离翼型表面的速度逐渐减小,这与空泡形成的阻碍有关。空泡尾部出现较大的漩涡区,靠近翼型表面的速度为负值,这与反向射流的作用有关。     

水翼间隙非定常涡空化流动的试验和数值研究

采用高速全流场显示技术对水翼叶顶间隙区域非定常旋涡空化流动进行观测,并提出一种旋涡空化模型来对空化流动进行计算.研究结果表明,提出的空化模型可以较好地捕捉非定常泄漏涡空化.此外,基于计算结果对水翼前缘云状附着空化与间隙涡空化的相互作用进行重点分析.结果表明,附着空化的发展和脱落会引起水翼叶顶附近迎背流面压力差值的变化,从而导致叶顶区域射流剪切层空化具有显著的非定常特性;附着空化末端闭合处的压力波动会造成泄漏空化涡强度的波动;附着空化的断裂所引起的吸力面局部高压会导致泄漏空化涡强度的显著降低.     

翼型间隙空化数值模拟的湍流模型比较分析

为了比较不同湍流模型对间隙泄漏空化计算的影响及间隙泄漏空化的流动特性,利用CFD方法对一侧有间隙的翼型绕流流动进行分析.选取具有不同攻角和间隙的翼型作为研究对象,并利用3种常用的湍流模型,结合修正后的空化模型,对翼型间隙泄漏涡引起的空化进行计算.计算结果表明:修正密度的RNG k-ε模型对间隙泄漏空化的形态预测及泄漏涡附近主流方向速度分布的预测更为准确;受壁面条件和翼型尾迹的影响,泄漏涡所在位置的主流方向速度在间隙为10 mm时具有较大值.小间隙条件下,翼型间隙内侧的剪切空化更为强烈,翼型上方的泄漏空化距翼型上表面更远.间隙增大至10 mm时,较高的涡中心速度使泄漏涡空化发展距离更远、溃灭更慢.翼型叶顶的圆角对间隙流动有一定的引导作用,使间隙内侧及其附近的速度分布更为均匀,并对空化的发生位置产生影响.     

非定常空化流动涡旋特性分析

基于试验和数值结果对三维非定常空化流动的涡旋特性进行分析.试验采用高速摄像系统对三维空穴发展过程进行观测.数值计算中采用混合密度分域的湍流模型和基于相间质量传输的空化模型.与试验结果对比,表明该计算方法能够准确捕捉云状空化的非定常过程：附着空穴的初生和发展,附着空穴的断裂,以及大尺度脱落空泡团的形成.基于计算结果,采用涡量分析法并引入复杂流场中的涡旋结构判据研究了空化现象的发生和发展对涡量输运过程的影响.结果表明：反向射流的形成和发展导致速度梯度的变化,从而引起涡量产生项的非定常变化;气液两相的相互转换会导致流场内体积变化率以及密度梯度的变化,同样引起体积变化率项和斜压矩项的非定常变化.采用速度梯度张量不变量Q和R分析空化现象的涡旋特性可以发现：在水翼前端附着空穴区域,旋转效应大于形变效应;脱落空泡位置附近,旋转效应和形变效应共同支配云状空泡运动.     

热敏感流体空化数值模拟及流动特性分析

为描述和捕捉热敏感流体空化的流动特性,文中对Zwart-Gerber-Belamri(ZGB)空化模型进行了热力学修正.通过将ZGB模型与修正的ZGB空化模型对液氮绕二维水翼空化流动的模拟结果与Hord的试验数据进行对比. 结果表明:修正的ZGB空化模型对于热敏空化流动的模拟结果与试验数据更加吻合,特别是能够较好的预测水翼表面温度和压力的分布. 采用修正的ZGB空化模型对氟化酮绕NACA 0015三维水翼的空化流动进行研究,并将空腔脱落的演化过程与Kelly的试验数据进行对比,模拟结果可以合理预测到试验过程中三维水翼附近空腔的脱落及其演化过程,进一步证明了修正的ZGB空化模型对于不同热敏流体的适用性.最后,对热敏流体空化周期性脱落的动态特征进行了识别和分析.结果表明:Ω的等值面与脱落的空腔形状相似,这表明脱落的空腔区域中呈现大规模的旋涡运动.空腔的生长脱落使水翼壁面产生明显的温降,B因子能够有效的预测腔内温降.     

绕水翼非定常空化流场数值模拟

基于Reynolds平均法,采用RNG k-ε模型及Mixture两相流模型,对绕攻角10°的ys930水翼非定常空化流动进行了数值模拟,分析了空化数分别为10,08,05时的空化流场结构、流动特性及空泡演化过程等．结果表明:云状空泡由两部分组成,前一部分为空泡主体,稳定附着于水翼吸力面上,其内部充满水蒸气,压力为汽化压力;后一部分为空泡附体,为周期性气液两相运动区域．空泡前端位置基本稳定,空泡主体长度随时间变化先增大后减小;空泡主体长度由空化数和回射流强度共同决定．空泡厚度随时间变化先增大后减小,空化数越小,空泡能达到的最大厚度越大,同时最大厚度出现位置越靠近水翼尾缘．回射流的强度与空化数成反比,空化数越小,回射流强度越大,来流与回射流的相互作用决定了空泡附体脱落的位置．