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.     

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

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

局部空泡水翼性能研究

空泡对水翼的水动力特性具有重要的影响,使用数值模拟方法研究空泡特性及其机理有着广泛的工程应用价值。对二维NACA翼型的模拟结果表明,基于势流理论的边界元积分方法适用于局部空泡水翼的性能研究,翼型几何参数变化对空泡长度和升力性能均有显著影响,在水翼设计与使用过程中必须对空泡发生与变化给予充分关注。     

一种计算水翼水动力的三维面元法

本文用Rankine奇点面元法计算了深、浅水中三维水翼的定常升力绕流。水翼的厚度和升力效应分别以水翼表面分布的Rankine源和法向偶极子来模拟，在自由表面上也分布Rankine源，通过满足相应的边界条件和尾缘处的Kutta条件求出这些奇点强度。以在自由表面下作小攻角定常运动的水翼为例进行了计算，计算结果与试验结果和其他计算结果作了比较。     

三维水翼升力线理论(英文)

Prandtl’s lifting line theory was generalized to the lifting problem of a three-dimensional hydrofoil in the presence of a free surface. Similar to the classical lifting theory, the singularity distribution method was utilized to solve two-dimensional lifting problems for the hydrofoil beneath the free surface at the air-water interface, and a lifting line theory was developed to correct three-dimensional effects of the hydrofoil with a large aspect ratio. Differing from the classical lifting theory, the main focus was on finding the three-dimensional Green function of the free surface induced by the steady motion of a system of horseshoe vortices under the free surface. Finally, numerical examples were given to show the relationship between the lift coefficient and submergence Froude numbers for 2-D and 3-D hydrofoils. If the submergence Froude number is small free surface effect will be significant registered as the increase of lift coefficient. The validity of these approaches was examined in comparison with the results calculated by other methods.     

解三维水翼绕流的下潜涡环栅格法

在包含于机翼表面内的某次表面及尾涡面上分布法向偶极子,在翼面上满足物面边界条件,建立了解三维机翼绕流的下潜涡环栅格法.考虑到非线性的自由液面边界条件,用下潜涡环栅格法求解水翼绕流问题.通过算例验证了方法的正确性和程序的可靠性.本方法可用于水翼及水下安定翼和各种舵的水动力计算.     

基于翼剖面改型的空化抑制

为提高水翼抗空化的性能,对二维翼型的吸力面外形进行适当改造。首先通过数值计算对稳态无空化流场和稳态空化流场进行模拟,计算所得的吸力面压力系数与实验值吻合良好,验证了模型的可行性。在此基础上,采取基于阻碍回射流从而控制空化的思路,在翼型吸力面上设置微小方形凸起,并提出设置拱弧的新方案。通过对原翼型及两种改型的空化流场瞬态模拟,对比了不同时刻各模型气体体积分数云图所反映出的翼面空化程度差异。计算结果验证了阻流体对云状空化的抑制作用,同时表明设置拱弧阻流体的效果比方形阻流体好。     

厚翼剖面的空泡特性研究

本文用奇点法解厚翼剖面局部空泡绕流问题。求解时精确考虑了空泡表面是一条流线的运动学边界条件,空泡表面压力为常数的动力学边界条件,以及水翼湿表面是一条流线的运动学边界条件,用迭代法求出了空泡表面的形状、压力分布等水动力特性,并与实验结果作了比较,得较好结果。     

Hydrodynamic response of a passive shape-adaptive composite hydrofoil

The primary objective of this paper is to present cavitation tunnel tests performed on an optimised shape-adaptive composite hydrofoil and compare the results to other composite hydrofoils. The optimised composite hydrofoil was designed based on prior literature and was manufactured using an optimised ply orientation schedule and a pre-twist. In the same experiment schedule a composite hydrofoil that has a ply orientation that is opposite to the optimised hydrofoil was also tested. In addition to the cavitation tunnel experiments, the paper also presents results predicted from FEA and CFD simulations for the optimised hydrofoil and compares the results from numerical methods to experiments. The results show that the optimised hydrofoil has an improved L/D ratio and a delayed stall phenomenon compared to other hydrofoils. Furthermore, due to the pre-twisted optimised geometry, the hydrofoil does not suffer from loss of lift at low angles of attack. The experimental results demonstrated the importance of characterising the performance of flexible shape-adaptive hydrofoils based on the actual velocity of the flow in addition to the conventional characterisation based on Reynold's number. Additional numerical simulations were performed to investigate the hydrofoils observed load dependant deformation behaviour. These results clearly show that for the same Reynold's number, the hydrofoil can have an appreciably different response if the flow velocity is different.     

振动翼推进性能计算及试验研究

采用非定常边界元法基于格林定理构建振动翼计算模型,讨论相关的构建过程、程序结构;采用时变尾涡面和时间步进法在时域中处理振动翼的尾涡存储效应,在尾涡面上应用等压Kutta条件并通过New-ton-Raphson法进行迭代求解,对在无限域中的振动翼水动力特性进行计算和研究;通过数值计算结果与试验数据对比,证实本方法的有效性。     

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

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.     

用等压Kutta条件求解三维水翼和螺旋桨片空泡

由于空泡计算的复杂性,国内外以前计算螺旋桨空泡时大多采用的是Morino's Kutta条件,但严格来说,Morino's Kutta只适用二维问题,对三维水翼和螺旋桨空泡的计算会带来一定误差,而等压Kutta条件用Newton-Raphson迭代过程,确保在螺旋桨随边上下表面的压力相等,可以消除Mofino's Kuaa条件所带来的误差,提高空泡在梢部的预报精度.文中用等压Kutta条件预报了三维水翼和螺旋桨的片空泡,并将两种Kutta条件得到的计算结果同试验结果进行了比较.     

高速轻型穿浪双体船纵向运动改善措施研究

针对在改善高速轻型穿浪双体船(WPC)迎浪中波长与船长接近时纵向运动幅度较大的缺点,采用了理论计算与模型试验相结合的方法,对250 t级穿浪双体船开展了水翼改善纵向运动的理论和试验研究,分析了水翼形式、尺度和安装位置等对纵向运动的影响规律。数值计算和试验结果的比较表明,计及水翼—船体水动力干扰影响的切片理论可满足WPC加水翼后波浪中纵向运动计算的需要,但在纵向运动响应峰值处数值计算结果偏高。模型试验表明,250 t级WPC加装水翼后,迎浪纵摇和垂荡有义幅值可减少20%～30%。     

浅窄航道非线性兴波阻力数值计算

应用Rankine面元法计算船舶在浅窄航道中作匀速直航运动时的非线性兴波阻力．基于势流理论，在船体表面、自由面以及岸壁分布Rankine源来表达流场速度势，利用迭代方法满足非线性自由面边界条件，采用映像法满足水底边界条件，采用网格错位法满足辐射条件．以Wigley数学船型为例进行计算，通过改变航速和航道宽度参数，得到相应的兴波阻力．通过分析不同航速下航道宽度变化对兴波阻力的影响，得到了一些有意义的结论．     

Numerical analysis of cavitating propeller and pressure fluctuation on ship stern using a simple surface panel method “SQCM”

This paper presents a calculation method for the pressure fluctuation induced by a cavitating propeller. This method consists of two steps: the first step is the calculation of propeller sheet cavitation, and the second step is the calculation of pressure fluctuation on the ship stern. It is for practicality that we divide the method into two steps but do not calculate these steps simultaneously. This method is based on a simple surface panel method “SQCM” which satisfies the Kutta condition easily. The SQCM consists of Hess and Smith type source panels on the propeller or cavity surface and discrete vortices on the camber surface according to Lan’s QCM (quasi-continuous vortex lattice method). In the first step, the cavity shape is solved by the boundary condition based on the free streamline theory. In order to get the accurate cavity shape near the tip of the propeller blade, the cross flow component is taken into consideration on the boundary condition. In the second step, we calculate the cavitating propeller and the hull surface flow simultaneously so as to calculate the pressure fluctuation including the interaction between the propeller and the hull. At that time, the cavity shape is changed at each time step using the calculated cavity shape gotten by the first step. Qualitative agreements are obtained between the calculated results and the experimental data regarding cavity shape, cavity volume and low order frequency components of the pressure fluctuation induced by the cavitating propeller.     

Automated manufacture of optimised shape-adaptive composite hydrofoils with curvilinear fibre paths for improved bend-twist performance

Composite marine propellers improve hydrodynamic efficiency by inducing bend-twist coupling and allowing for passive pitch changes. One critical limitation, however, is the extent to which a composite propeller blade can deform and cause a pitch change without incurring structural failure. Recent numerical studies showed that curvilinear tows could improve the structural response of a composite blade by lowering its deflection or stress and strain required to induce a pitch change, but no experimental validation has been carried out before. The current study, thus, presents the manufacture of composite sandwich hydrofoils made with steered tows using automated fibre placement and validates the curvilinear tow benefits. Two hydrofoils were optimised with straight and curved fibre path layups, respectively and were manufactured for mechanical testing. The manufacturing complications arising from steering curvilinear tows in a three-dimensional convex mould are also discussed in the paper. The study found that significant tow buckling occurred near the tool cavity edge due to excessive steering radius during manufacture. The follow-up structural cantilevered tests showed that the experimental results were consistent with the FE predictions despite the presence of some manufacturing defects. The experiment agreed that the hydrofoil manufactured with curved tows achieved a similar tip twist but a significant reduction in deflection and critical principal strains compared to the hydrofoil made with straight tows. The use of a foam core reduced the overall weight of the sandwich hydrofoils by about 25% compared to that of a fully-carbon composite hydrofoil, and the numerical analysis showed that the core shear failure induced by transverse shear stresses was unlikely to occur.     

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

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

基于MSC.Patran的水翼升力求解方法研究

根据已有的水翼升力面理论计算方法及成果,对水翼模型进行部分简化,并借助MSC.Patran实现了水翼模型和流场计算域的参数化建模,设置边界和工况后,对水翼航行中的升力使用有限元法进行迭代求解。与已有规律和试验结果对比表明,使用参数化建模和有限元法可以计算出水翼的升力以及自由面的形变数据。不仅减少了手工劳动,而且预报了水翼的升力及分布。     

近自由表面大展弦比超空泡水翼的研究

本文采用非线性超空泡水翼二维解,考虑了重力效应及水翼兴波速度势,应用升力线理论对近自由表面三维超空泡水翼进行了研究,并以平板水翼为例作了计算,所得的结果比Furuya1975年的工作有所推进。