小攻角下水下高速航行体超空泡流特性研究

为研究小攻角下水下高速航行体超空泡形态及水动力特性,利用商业CFD软件Fluent6.2,对小攻角下高速航行体超空泡流进行数值模拟,分析了空化数、攻角对水下高速航行体空泡形态以及水动力特性的影响规律.研究表明,攻角为能够明显改变空泡形态的轴对称性;攻角越大,航行体的阻力系数也增大,不利于超空泡的减阻,甚至会导致航行体的失稳.研究结果将为开展水下高速航行体超空泡实验研究提供理论参考.     

水下高速航行体超空泡减阻特性数值模拟研究

利用CFD商业软件Fluent6.2对水下高速航行体超空泡流进行了数值模拟研究,计算了带圆盘与圆锥两种头部空化器航行体的阻力特性,详细分析了空化器直径、锥角、航行体长细比对超空泡减阻特性的影响,并计算了高速水下航行体自然超空泡减阻率,结果表明,在超空泡形态下,带圆盘空化器头部的水下高速航行体,更加有利于超空泡的减阻,超空泡减阻率可达95%以上.研究结果为进一步研究水下高速航行体的结构设计和水动力布局提供了理论参考.     

三维空泡回射流的理论研究

空泡回射流以某个方向冲击空泡壁面是导致空泡发生不同脱落方式的重要原因,但是空泡回射流的方向目前认为是一个不确定量,需要人为指定.该文针对无界重力流场中的三维空泡,基于势流假设和积分方程推导了一个计算回射流角度和回射流截面面积的三角函数代数方程组,并且在小角度条件下给出了回射流角度和截面面积的理论解析式.研究表明,空泡回射流角度与空化器的攻角、阻力和空泡浮力有关.     

小攻角水下航行体定常空泡外形计算方法

水下航行体上的空泡在攻角、重力及航行体的外形影响下,空泡的迎水面外形与背水面外形并不对称,这种不对称空泡外形对空泡水动力或航行体弹道有极大影响.本文结合动量原理和空泡独立膨胀原理,在小攻角条件下建立了定常三维空泡外形的方法.     

攻角空化器的超空泡形态非对称特性研究

为了满足水下高速超空化航行体运动平衡性和稳定性要求,必须采用非对称超空泡流动模式,采用具有非零攻角的空化器是实现该要求的一种重要策略.基于冲量定理,理论分析了非零攻角空化器引起的超空泡轴向变形量随轴向长度的关系,获得了空化器攻角对超空泡非对称性及其影响因素之间的定性和定量关系,其一阶近似与同类文献的结论是一致的,数值模拟表明了所得结论的合理性和分析的正确性.研究表明,非零攻角空化器引起的超空泡轴线变形量与沿空泡长度的距离呈线性规律增加,并在空泡尾部达到最大,不可能单独利用空化器攻角沿空泡长度来中和重力影响.     

超空泡航行体操纵过程流体动力特性数值模拟研究

超空泡航行体操舵过程会引起空泡变形,导致航行体流体动力学特性发生变化。为了了解超空泡航行体操舵过程中航行体的动力学特性,文中采用基于欧拉两流体模型的CFD数值模拟方法及动网格技术对超空泡航行体空化器、尾舵操舵过程以及航行体攻角变化过程中的空泡形态及航行体瞬态流体动力特性变化规律进行了研究。研究结果表明空化器操舵过程中空化器升力随偏转角基本呈线性规律变化,对航行体尾部滑行力的影响相对于攻角变化对滑行力的影响为小量;尾舵操舵过程改变了空泡尾部流场,对于航行体尾部滑行力会产生重要影响。     

小攻角下航行体三维非定常空泡形态理论预示方法

小攻角下水下航行体高速运动时会形成不对称、非定常空泡,空泡形态是水下流体动力及弹道设计的主要依据。该文基于空泡独立膨胀原理,考虑横流对独立空泡发展的影响,建立了带攻角状态下空泡形态理论计算模型。针对典型工况开展计算获得了航行体迎背流面空泡长度、空泡压力变化过程,并与试验数据进行比对以验证模型的合理性。     

小攻角水下航行体定常空泡外形计算方法

水下航行体上的空泡在攻角、重力及航行体的外形影响下,空泡的迎水面外形与背水面外形并不对称,这种不对称空泡外形对空泡水动力或航行体弹道有极大影响。本文结合动量原理和空泡独立膨胀原理,在小攻角条件下建立了定常三维空泡外形的方法。     

基于均相流输运模型的二维水翼空化模拟

基于均相流输运模型对NACA661-012型水翼进行了空化数值模拟,计算了不同空化数K和不同攻角α下的升力系数与阻力系数.从流场、压力场和边界层的变化分析了空化的产生发展对水翼升力和阻力的影响.当攻角α大于8°以后,空化流动的不稳定性增强,空泡呈现出周期性的生长溃灭,并伴有升力系数和阻力系数的周期性波动.计算结果与实验进行了对比,结果吻合较好.     

超空泡航行体楔形舵片流体动力学特性数值模拟

舵片是保证超空泡航行体运动稳定性和控制航行弹道的重要部分。文章基于均质平衡流模型和SST(Shear Stress Transport)湍流模型,计算了单独舵片的流体动力特性,并与试验数据进行了对比,结果符合较好,验证了计算模型的有效性。基于此方法,计算了单独舵片发生空化后在不同操舵状态下的非定常流体动力变化。结果表明,在攻角相同时,操舵状态下舵片的非定常升力系数和定常结果差别不大,而非定常阻力系数大于定常结果,并且操舵速度越快,阻力系数越大。另外计算了舵片发生空化后的流体动力系数,结果显示在攻角相同时,舵片的阻力系数和升力系数均小于其在全湿状态下的结果;在空化状态下,舵片升力系数的斜率小于全湿状态,并且舵片升力系数的斜率是变化的,存在某临界攻角,攻角大于此临界值时,升力系数的斜率减小,而此临界攻角恰好为舵片的吸力面刚刚出现空化时的攻角;操舵状态下舵片的阻力系数和升力系数的变化规律与定常结果一致,但是数值偏小。     

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

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.     

基于混合迭代法的二维水翼空泡性能分析(英文)

In order to study cavitation characteristics of a 2-D hydrofoil, the method that combines nonlinear cavitation model and mixed-iteration is used to predict and analyze the cavitation performance of hydrofoils. The cavitation elements are nonlinearly disposed based on the Green formula and perturbation potential panel method. At the same time, the method that combines cavity shape for fixed cavity length (CSCL) iteration and cavity shape for fixed cavitation number (CSCN) iteration is used to work out the thickness and length of hydrofoil cavitations. Through analysis of calculation results, it can be concluded that the jump of pressure and velocity potentially exist between cavitation end area and non-cavitations area on suction surface when cavitation occurs on hydrofoil. In certain angles of attack, the cavitation number has a negative impact on the length of cavitations. And under the same angle of attack and cavitation number, the bigger the thickness of the hydrofoil, the shorter the cavitations length.     

螺旋桨设计参数对桨叶片空泡性能的影响分析

文章基于扰动速度势面元法建立了在均流条件下螺旋桨桨叶片空泡数值预报方法,空泡模型采用压力恢复闭合模型。通过对5600TEU集装箱船螺旋桨空泡的数值预报,以及与试验结果的比较,验证了该方法的可行性。该方法能够较为快速准确地预报螺旋桨桨叶片空泡,可用于分析参数对螺旋桨空泡性能的影响,为抑制螺旋桨空化设计提供基础。在此基础上重点分析了桨叶侧斜、纵倾以及桨叶剖面型式对螺旋桨空泡性能的影响,计算表明加大侧斜能够减少空泡面积,空泡向外半径偏移;桨叶剖面的设计对空泡性能影响较大,优化设计桨叶剖面可以有效减少空泡面积,提高螺旋桨抗空化能力;纵倾向压力面弯曲的分布形式可以改善梢部的压力分布,减少叶梢附近空泡长度,从而可望减少由空泡引起的脉动压力。     

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

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

涂层对空化流动特性影响的实验研究

利用实验的方法研究了涂层对绕水翼空化流动特性的影响。分别针对喷涂环氧涂层和氟碳涂层的 Clark-Y 型水翼,采用高速摄像装置观察了不同空化阶段的空化流动形态。研究结果表明：（1）在初生空化阶段,当σ=1.82时,沿环氧涂层水翼表面展向排列着初生空泡,而氟碳涂层水翼还处于无空化状态,说明相对于环氧涂层,氟碳涂层对空化现象的产生有一定的抑制作用,氟碳涂层水翼初生空化数为1.50;（2）在片状附着型阶段,当σ小于1.63时,绕环氧涂层水翼的空化先于氟碳涂层水翼发展至片状空化,绕水翼空化流动产生大量分散空泡,沿水翼表面向后运动过程中逐渐长大,在高压区溃灭后形成小空泡并以马蹄涡形式继续运动。同一空化数下,绕环氧涂层水翼空化流动的空泡长度大于氟碳涂层水翼。但随空化数降低,两者空泡长度逐渐接近,说明环氧涂层在片状空化阶段对空化的抑制作用逐渐增强;（3）σ=0.87时空化发展至云状空化阶段,空化流动伴随周期性的云状空泡的脱落,绕环氧涂层水翼的空化流动周期及无量纲空化面积均小于氟碳涂层水翼,说明涂层对空化的非定常变化也有一定的抑制作用,且环氧涂层强于氟碳涂层。     

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

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

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.     

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

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.