自然空化流动数值模拟中参数取值影响的研究

探索了应用商业软件Fluent6.2进行自然空化计算的方法.通过设定不同气核质量分数对半球头回转体的自然局部空化作了模拟,研究了气核对自然空化计算结果的影响.通过设定边界湍流参数计算了平头回转体的自然超空泡,研究了湍流参数设定对自然空化计算结果的影响.通过对圆盘空化器的超空泡计算,得出结论:气核含量主要影响超空泡尾部形态,湍流参数主要影响空泡尺度.多个算例的计算结果与相应的实验结果以及经验公式计算结果进行了比较,符合良好.     

泡状流中空化与激波相互作用的数值模拟

文章针对定常可压缩泡状流中回转体的空化绕流进行了数值模拟,在给定来流速度和环境压力的条件下,对不同含气率下的空化与激波的相互作用进行了研究。首先,对含气率为0的情形进行计算并同试验数据进行对比,证明所采用的计算模型是可信的。其次,改变流场含气率从0至0.5进行计算,结果表明,随着含气率的增大,流场的可压缩性随之增强,体现为数值模拟得到了回转体绕流流场的3种波系,包括回转体头部处的脱体激波、分离面处的膨胀波和空泡尾端的斜激波。而空化与激波的相互作用则体现为:头激波削弱了空化效应,使得空化区域减小;由于空泡外形的影响,使得空泡尾端出现了斜激波;含气率超过一定值,空泡已经无法闭合在物面上,而是闭合在空泡尾端的激波面上,体现为空泡尾端壁面逆压梯度趋于平缓。计算结果揭示了泡状流中空化与激波相互作用的新的物理现象。     

回转体空化研究现状及展望

针对水下高速运动回转体周围必然产生空化现象的问题,阐述了空化现象的主要负面影响,详细介绍了回转体空化实验装置、回转体水下空泡流和出入水研究进展及成果,简要分析了回转体空化研究的发展趋势。     

水下航行器局部空泡流场的空泡尾流模型初探

基于二维翼型局部空泡绕流场的空泡尾流模型,提出了几种可用于轴对称回转体局部空泡流场计算的空泡尾流模型,并对这些模型进行了探索性的理沦分析,给出了应用其中的压力恢复闭合模型求解回转体局部空泡绕流场的算例,通过比较发现与实测值的一致性较好,说明本文提出的空泡尾流模型适合于求解轴对称回转体线型的局部空泡流场。     

超空泡回转体减阻特性研究

基于Fluent 6.0的气泡两相流模型对超空泡回转体的减阻特性进行了数值研究.内容包括:外形对阻力及超空泡形状的影响;超空泡控制;阻力系数随空泡数的变化规律;长细比对减阻率的影响.超空泡减阻机理的研究表明:超空泡不仅可以减小回转体的摩擦阻力,还可以减小回转体的压差阻力.在外形、长细比和空泡数以及工程可实现性等诸多因素中,存在着一个最佳组合,使减阻率最高.     

Level set方法在空化流动计算中的应用

文章将Level set方法应用于自然空化与通气空化流动计算中,通过耦合入相间的质量传输方程,建立了一个用于描述气、汽、液多相空化流动的数值模拟方法,分别对绕圆头回转体和Clark-Y型水翼的自然空化流动和绕圆盘的通气空化流动进行了数值模拟研究,并与实验结果进行了对比,研究结果表明:相较常用的多相流模型(mixture model),Lev-el set方法通过时间和空间上的压缩离散方案,减小了相间界面的扩散,准确地捕捉到了相间的界面,可以有效地运用于空化多相流动。     

通气空泡流中洞壁效应对空化数和压力场影响的研究

文章在均质平衡多相流模型的基础上耦合输运方程型空化模型,通过求解混合介质的RANS方程、RNG k-ε湍流输运方程以及各相的质量输运方程,采用通用CFD软件—FLUENT数值模拟了水洞中带圆盘空化器航行体模型的定常通气空泡流动,研究了圆形截面闭式空泡水洞中洞壁效应对通气空化数和压力场的影响。得到的阻塞空化数线性正比于圆盘水洞直径比,且与三维圆盘自然空泡流的势流近似解基本一致,分析了洞壁效应作用下空化流场内的压力分布特点,并根据计算结果拟合了一定适用条件下通气空泡长度、最大直径和模型阻力系数的近似公式。     

空化模型在非定常空化流动计算的应用评价与分析

文章基于试验观测数据评价了Kubota与Singhal两种空化模型在非定常云状空化流动数值模拟中的应用。采用商业软件的二次开发技术将两种空化模型引入了计算软件,针对绕Clark-Y水翼的云状空化流动进行数值计算,并与水洞试验结果进行了对比。结果表明,两种空化模型计算的云状空化阶段的流场时均速度分布以及涡量分布具有明显的不同。采用Singhal空化模型可以得到和试验观测更加相近的空化旋涡区与空泡云的旋涡分离的脱落形式。     

大攻角水下航行体自然空泡流的数值模拟研究

基于相分数输运方程型的均质平衡流空化模型,采用有限体积法研制了大型空泡流计算程序,对大攻角下运行的水下航行体三维空泡流进行了数值模拟,并与实验结果进行了对比.首次将非线性涡粘湍流模式与基于Rayleigh-Pies-set方程的TEM型空化模型相结合,建立了自然空泡流的数学模型.采用基于SIMPLE的压力-速度-密度耦合修正算法、二阶精度三时间层格式以及基于延迟修正的高阶对流TVD格式.计算模拟了0.2～0.6空化数、4°～20°攻角的不同工况,得到的三维空泡形状及压力分布与实验结果相符.研究了大攻角下航行体周向上的空泡形态分布特征,给出了多种空泡尺度和升阻系数与空化数和攻角之间的关系.通过定量分析发现,空泡的不对称性导致航行体某些部位受力集中,表明高速带空泡运动的航行体在大攻角运动中其结构将受到巨大的水动力载荷.计算还发现,大攻角下的阻力系数与空化数之间的关系和零攻角条件下刚好相反,并根据空泡的不对称性从形状阻力与粘性阻力的关系上对这种现象作出了解释.     

螺旋桨梢涡空泡初生及尺度效应研究

应用气泡动力学方程和气泡运动微分方程研究了螺旋桨梢涡空泡的初生问题,并根据简化Reyleigh-Plesset方程推导了不同尺度模型空化初生空泡数换算公式,建立了空化初生尺度换算模型,研究了螺旋桨梢涡空泡初生尺度效应问题。螺旋桨梢涡流场应用RANS方法求解,湍流模型为经过旋转和曲率修正的代数雷诺应力模型(EARSM-CC)。计算结果表明,文中建立的数值方法能够准确预报出螺旋桨梢涡流场分布;螺旋桨梢涡空泡初生空泡数计算结果高于试验观察值;应用文中建立的空化初生尺度换算模型得到的结果与数值模拟结果基本吻合,而且其结果偏于安全。     

二维水翼型空化流的数值计算(英文)

In order to predict the effects of cavitation on a hydrofoil, the state equations of the cavitation model were combined with a linear viscous turbulent method for mixed fluids in the computational fluid dynamics (CFD) software FLUENT to simulate steady cavitating flow. At a fixed attack angle, pressure distributions and volume fractions of vapor at different cavitation numbers were simulated, and the results on foil sections agreed well with experimental data. In addition, at the various cavitation numbers, the vapor fractions at different attack angles were also predicted. The vapor region moved towards the front of the airfoil and the length of the cavity grew with increased attack angle. The results show that this method of applying FLUENT to simulate cavitation is reliable.     

水下爆炸冲击波产生的区域空化效应二次加载的并行数值研究

船舶等水面结构受到水下爆炸作用时,其附近水域会产生空化现象,包括局部空化和区域空化。以某水面船为对象,数值模拟了该船受水下爆炸作用时的动态响应。比较了考虑空化效应和不考虑空化效应下船体结构的不同响应,分析了区域空化效应对此船的二次加载作用,二次加载引起的加速度峰值与冲击波引起的加速度峰值在同一量级上。使用高性能计算机并行计算比较了不同流体网格密度模型的计算结果。结果对比显示,为了更精确地模拟空化效应及其对结构的影响,加密水域网格是必要的。     

带前置定子的导管桨片空泡数值预报及梢涡空泡起始判别研究

对水下航行体后带前置定子的导管桨片空泡和梢涡空泡起始进行研究。首先,采用基于混合网格的RANS求解器结合Singal空泡模型,数值模拟带前置定子导管桨片空泡形态,并与试验结果进行对比,预报的片空泡形态与试验结果吻合良好,表明用该方法预报水下航行体后带前置定子导管桨片空泡可行;其次,对水下航行体后带前置定子的导管桨多个方案全湿流动进行了数值模拟,通过梢部端面的最小压力值判断空泡起始,与试验结果进行比对,其吻合良好,表明可用该方法判别梢涡空泡起始。     

船舶柴油机喷油泵柱塞穴蚀机理的数值仿真研究

针对船舶柴油机喷油泵柱塞螺旋槽上方壁面的穴蚀现象,建立了喷油泵柱塞与套筒模型,利用FLUENT软件对喷油泵内流场进行数值仿真计算,重点研究局部低压区、湍流涡团以及气泡析出、成长和溃灭过程。研究结果表明,柱塞螺旋槽与套筒回油孔构成的重叠开口的大小和形状对燃油流场具有决定性影响,产生低压区和湍流,生成的气泡溃灭时对柱塞螺旋槽上方壁面和回油孔内壁产生机械冲击,造成穴蚀。     

三维水翼片空泡尺度效应研究

针对三维水翼的片空泡尺度效应问题,建立LES湍流模型对三维水翼的片空泡生成情况进行计算研究。模拟了片空泡在一个周期内的生成和脱落过程,空泡计算结果与试验做了详细对比,发现LES湍流模型的计算结果与实验结果吻合较好,能较为完整的显示片空泡的脱落过程。研究了来流速度和模型尺度的变化对片空泡产生的影响,发现速度对初始片空泡数基本无影响,而模型尺度对片空泡的形成和发展有影响,但该影响并不大。在此基础上推导了新的初始片空泡数的尺度效应公式。     

G4135柴油机燃烧过程数值模拟的研究

采用三维CFD软件,以G4135直喷式柴油机为对象,运用多维燃烧模型,对缸内燃烧过程进行数值模拟。计算了不同曲轴转角下的气体速度、压力、温度、氧气浓度和碳烟浓度等,计算所得数据与相关文献相吻合。结果表明FLUENT所提供的燃烧模型可以作为预测柴油机缸内燃烧的一种有效手段,并为进一步研究同类型柴油机的燃烧性能提供了依据。     

水翼部分空化流动二维数值仿真研究（英文）

In the present study, a new approach is applied to the cavity prediction for two-dimensional(2D) hydrofoils by the potential based boundary element method(BEM). The boundary element method is treated with the source and doublet distributions on the panel surface and cavity surface by the use of the Dirichlet type boundary conditions. An iterative solution approach is used to determine the cavity shape on partially cavitating hydrofoils. In the case of a specified cavitation number and cavity length, the iterative solution method proceeds by addition or subtraction of a displacement thickness on the cavity surface of the hydrofoil. The appropriate cavity shape is obtained by the dynamic boundary condition of the cavity surface and the kinematic boundary condition of the whole foil surface including the cavity. For a given cavitation number the cavity length of the 2D hydrofoil is determined according to the minimum error criterion among different cavity lengths, which satisfies the dynamic boundary condition on the cavity surface. The NACA 16006, NACA 16012 and NACA 16015 hydrofoil sections are investigated for two angles of attack. The results are compared with other potential based boundary element codes, the PCPAN and a commercial CFD code(FLUENT). Consequently, it has been shown that the results obtained from the two dimensional approach are consistent with those obtained from the others.     

Numerical simulation of the flow around two-dimensional partially cavitating hydrofoils

In the present study, a new approach is applied to the cavity prediction for two-dimensional （2D） hydrofoils by the potential based boundary element method （BEM）. The boundary element method is treated with the source and doublet distributions on the panel surface and cavity surface by usethe of the Dirichlet type boundary conditions. An iterative solution approach is used to determine the cavity shape on partially cavitating hydrofoils. In the case of a specified cavitation number and cavity length, the iterative solution method proceeds by addition or subtraction of a displacement thickness on the cavity surface of the hydrofoil. The appropriate cavity shape is obtained by the dynamic boundary condition of the cavity surface and the kinematic boundary condition of the whole foil surface including the cavity. For a given cavitation number the cavity length of the 2D hydrofoil is determined according to the minimum error criterion among different cavity lengths, which satisfies the dynamic boundary condition on the cavity surface. The NACA 16006, NACA 16012 and NACA 16015 hydrofoil sections are investigated for two angles of attack. The results are compared with other potential based boundary element codes, the PCPAN and a commercial CFD code （FLUENT）. Consequently, it has been shown that the results obtained from the two dimensional approach are consistent with those obtained from the others.