基于滤波器湍流模型的水翼空化数值模拟

采用界面捕获法,在均相流模型假设下,引入滤波函数修正了RNG k-ε湍流模型,基于正压流体假设的状态方程空化模型.采用SMPLEC算法,数值求解雷诺平均的Navier-Stokes方程,模拟了二维NACA66(Mod)翼型的定常、非定常空化流动.计算得到的翼型表面压力分布与试验结果基本一致,较好地模拟了非定常空化云的初生、发展、断裂和脱落的周期性过程,非定常流动的斯特劳哈数与试验结果相吻合.给出了2个时刻的翼型表面速度矢量场,分析了空泡内部流动结构,结果发现翼型空泡尾部的反向射流是引起翼型表面局部压力升高,进而导致空泡发生断裂及空化云脱落的重要原因.     

滤波器湍流模型在超空化流动计算中的应用评价

对滤波器湍流模型(FBM)在超空化流动计算中的应用进行分析和评价.在CFX软件中,分别采用FBM模型和标准κ-ε模型对Hydronautics翼形的超空化流动进行模拟,结合水洞实验结果从空泡形态和流场结构两方面对两种模型的计算结果作对比分析.结果表明,FBM模型对空泡形态和空泡尺度的模拟更符合实际.FBM模型可以捕捉超空化流动中较大尺度的涡团结构,并能更清晰地模拟出空化区尾部涡团交替脱落的非定常细节,同时.FBM模型对超空化流动中反向射流现象的模拟更敏感.     

绕水翼空化流动的数值模拟

基于正压流体假设的状态方程空化模型,考虑空泡流可压缩性的影响修正了RNGk～ε湍流模型,采用SMPLEC算法,数值求解Reynolds平均的Navier-Stokes方程,模拟了二维NACA66(Mod)翼型的定常、非定常空化流动。计算得到的翼型表面压力分布与试验结果相吻合,较好地模拟了非定常空化云的长大、断裂、涡状空化团脱落和破裂的周期性过程,分析讨论了空化云的脱落机理,指出反向射流是引起空化云脱落的重要原因。     

空泡半径非线性变化的空化模型及应用

为了提高Schnerr-Sauer空化模型模拟空化流动的能力,基于Rayleigh-Plesset方程和均相流假设建立了一种考虑空泡半径随时间非线性变化的空化模型.联立滤波器湍流模型,分别采用该非线性空化模型和原始Schnerr-Sauer模型对绕二维Clark-Y翼型空化流进行了非定常数值计算,得到了时均化升阻力系数随空化数的变化曲线、云空化时空泡周期性演变及其沿弦长方向速度分布规律.与其他的实验结果进行对比分析,结果表明:采用该非线性模型计算得到的时均化升阻力系数与实验值符合程度更好,并能较准确地捕捉升阻力系数在不同空化阶段的变化规律;该非线性模型能准确地模拟翼型尾缘空化云的成长、脱落和溃灭的全过程,且能更准确地模拟翼型尾流的速度场.     

非均匀流场螺旋桨空泡数值模拟

以某油船螺旋桨为研究对象,基于Rayleigh-Plesset方程,采用Schnerr-Sauer空化模型和可实现的k-ε两层湍流模型,利用计算流体力学(CFD)软件STAR-CCM+模拟了非均匀来流条件下螺旋桨的空泡形态.通过对螺旋桨叶梢区域进行有效合理的网格加密,以较少的网格数量成功捕获梢涡空泡.数值计算与试验结果对比表明:计算结果准确再现了桨叶进出伴流区域空泡初生、发展和溃灭的整个历程;每个相位角叶背片空泡形态与试验观察吻合,片空泡面积相差在5%以内;该数值方法虽然能够捕捉到梢涡空泡,但是还不能对梢涡空泡的非定常特性和空间结构进行准确预测.基于上述结果,该方法适用于非均匀流螺旋桨空泡流动模拟.     

绕三维扭曲水翼非定常空泡脱落特性的数值研究

为了进一步探究绕三维扭曲水翼空化流动结构的非定常演变特性,基于均质流模型,采用大涡模拟(LES)方法,并耦合Zwart空化模型对绕扭曲水翼(NACA0009)的非定常空化流动进行了数值计算。结果表明:与试验结果对比发现,该数值计算方法可以很好地捕捉到绕NACA0009水翼的非定常空泡形态的演变过程,且计算得到的水翼升力系数的时均值及其特征频率均与实验值吻合较好。非定常空泡的脱落过程分为主脱落和二次脱落两个阶段,主脱落过程表现为大尺度空泡团的断裂,呈U型结构卷起并脱落,此时水翼表面的附着空泡呈凹形结构。随着时间的推移,附着空泡的两侧发生小尺度空泡团的二次脱落,产生局部高压,抑制了当地空泡的生长速度,从而使附着空泡由凹形结构逐渐发展为凸形结构。     

基于实验观测的分域湍流模型在通气超空化中的评价

为建立一种可以准确高效预测通气超空化流动的湍流模型,结合数值计算和实验结果,对绕锥头回转体通气超空化流动特性进行研究.实验采用高速录像观察了通气空化随时间的流场变化;数值计算中,分别应用标准k-ε湍流模型和密度分域的湍流模型计算了绕锥头回转体通气空化流动.其中,密度分域的湍流模型是在实验观测的基础上建立,即在空泡的前端含气量较大的区域应用DCM模型,以体现附着型空穴的可压缩性;在空泡尾部含气量较大的雾状空泡区域应用FBM湍流模型,以捕捉多尺度的空泡涡团结构.研究结果表明:与标准k-ε湍流模型相比,基于密度分域的湍流模型计算的结果与实验观测的现象基本吻合,有效减小了通气空化空泡区域内的湍流黏性,可以捕捉空泡区域内多尺度旋涡结构的演化过程,进而可以准确地预测通气超空化空泡断裂脱落的非定常流动细节.      

考虑当地流动特征的Zwart-Gerbera-Belamri空化模型改进及应用

为减小Zwart-Gerbera-Belamri空化模型中经验参数取值对数值模拟结果的影响,基于当地流动特征修正了Zwart-Gerbera-Belamri空化模型.联立Realizable k-ε湍流模型,分别采用修正前后的Zwart-Gerbera-Belamri空化模型对二维NACA0009 MOD水翼空化流场进行了数值模拟,空化数范围为0.75~0.90.并采用修正前后的Zwart-Gerbera-Belamri空化模型对二维Clark-Y水翼的初生空化、片状空化、云空化进行非定常数值模拟,空化数范围为3.00~0.55.将数值模拟结果与已有实验结果进行对比,结果表明:引入当地流动特征可以有效减小Zwart-Gerbera-Belamri空化模型中的气泡半径及气核体积分数的取值对数值模拟结果的影响,且修正后的空化模型能够捕捉不同空化数下绕Clark-Y翼型的时均升阻力系数的变化趋势.      

基于PANS模型的水翼非定常空化特性研究

采用一种局部时均化模型(PANS)对绕二维Clark-y型水翼的空化流动进行数值计算.分别研究了空化模型中最大密度比对空化计算的影响,PANS模型中关键控制参数fk对云空化非定常特性的影响和云空化阶段空泡形态及升力系数随时间的演变,并与实验结果进行对比.结果表明:最大密度比会影响液相与气相之间的质量传递速率,将最大密度比提高,可大大改善数值预测的精度;随着fk值的减小,PANS模型可释放更多的湍流尺度,较好地预测云空化中强烈的非定常特性;由于云空化呈现出产生-发展-脱落-溃灭的准周期性变化,使得水翼的升力系数变化剧烈,计算出时均升力系数为0.708,与实验值0.760相差约7%;回射流与主流相互作用,在水翼尾部形成的一对旋转方向相反的旋涡是云状空泡形成与发展的主要原因.     

基于密度分域的混合湍流模型的应用与评价

结合数值计算和实验结果评价了一种基于密度分域的混合湍流模型.实验采用高速录像和LDV测量技术观察了绕Clark-y型水翼云状空化随时间的流场结构变化和速度分布.数值模拟中,基于非定常空化的流动现象,建立了一种基于密度分域的混合湍流模型,并分别应用不同的湍流模型计算了绕水翼非定常云状空化流动.通过与实验结果的对比表明,分域湍流模型能够更好地调整流场内的湍流黏性,精确捕捉到空穴形态和空泡脱落的非定常细节,时均流场速度分布与实验值基本吻合.