螺旋桨非定常轴承力计算

采用扰动速度势面元法计算螺旋桨非定常轴承力,桨叶、桨毂和尾涡面由双曲四边形面元进行离散,对时域内非定常问题的求解采用时间步进迭代方法,建立了满足桨叶随边非定常等压库塔条件的非线性迭代结构,使迭代求解更加有效、快速和稳定。     

面元法预报螺旋桨表面非定常压力分布

本文建立了扰动速度势面元法预报螺旋桨表面非定常压力分布的理论和数值方法。该方法把桨叶和桨毂表面离散为若干四边形双曲面元,每个面元上布置等强度源汇和偶极子分布,螺旋桨尾涡面也离散为布置等强度偶极子的四边形双曲面元。所有的时域参数通过傅里叶级数展开转化为频域参数,使得在时间域内的求解转化为每一谐调阶上的求解。桨叶随边处通过迭代和采用广义逆矩阵方法在每一谐调阶上实现非线性等压库塔条件。桨叶表面非定常压力     

全方向推进器非定常水动力性能的面元预报方法

研究了全方向推进器非定常水动力性能的面元预报方法,基于螺旋桨面元法建立了全方向推进器的非定常水动力性能计算的数学模型,对全方向推进器的非定常水动力性能进行了数值预报。采用了关于扰动速度势的基本积分微分方程,并采用双曲面元以消除面元间的缝隙。用Newton-Raphson迭代过程在桨叶随边满足压力Kutta条件。在计算面元的影响系数时,应用Morino导出的解析计算公式加快了数值计算的速度。为避免数值求导中的奇异性,用Yanagizawa方法求得物体表面上的速度分布。本文计算结果与日本水池模型试验结果、升力线方法计算结果及升力面方法计算结果进行了对比。     

带定子导管螺旋桨定常和非定常性能预估的基于速度势的面元法

开发了一个带定子导管螺旋桨定常和非定常水动力性能预报的数值计算方法,该方法采用基于速度势的面元法分别处理绕对带桨毂螺旋桨、导管和定子周围的流动,并通过迭代计算处理它们之间的相互干扰影响.为简化计算,采用了诱导速度的周向平均值来考虑螺旋桨、导管和定子之间定常水动力相互干扰影响.它们之间的非定常水动力相互干扰通过时间域内迭代方法求解,然而螺旋桨泄出涡、导管泄出涡和定子泄出涡之间的非定常干扰采用了时域内简化处理,这样的简化处理既保证了计算精度,又大幅度地减少了计算时间.对于螺旋桨、导管和定子的计算程序分别依次逐个运行直至每个部分的水动力收敛.算例考核表明该数值方法的计算结果与试验数据吻合较好.     

非定常螺旋桨表面压力面元法计算

本文建立了非定常螺旋桨势函数面元法的理论方法和数值方法，发展了相应的计算机程序，采用通过迭代实现的非线性压力库塔条件，计算了一常规桨和一大侧斜桨的桨叶表面非定压力，并与测试结果及人的计算结果进行了比较。     

用RANS方法预报螺旋桨非定常空泡

基于空泡混合模型用RANS方法对螺旋桨非定常空泡进行了预报.采用尺度函数以渐进扩张的方式对桨叶表面进行网格划分,在导边、随边、叶根、叶梢等部位加密,而在桨叶表面中部单元逐渐增大.为了更好地考虑湍流边界层的流动特性,构建了棱柱层和四面体混合非结构化网格.为节约划分网格的时间,建立了一种几乎不需人为干预的快速划分螺旋桨计算域网格的方法.对非均匀进口来流用自定义函数进行设置,空泡的预报结果与试验结果吻合较好,证明该方法在螺旋桨非定常空泡计算方面的可行性和良好的应用前景.     

螺旋桨水动力性能的数值预报方法

基于速度势的低阶面元法预报螺旋桨的水动力性能。选用四边形双曲面元对桨叶进行离散以消除面元间的缝隙,基本积分方程由格林公式导出。在面元上布置等强度源汇和偶极子。采用线性尾涡并在每个尾涡面元上布置等强度的偶极子。利用Newton-Raphson迭代过程满足桨叶随边非线性等压kutta条件,使桨叶上下表面的压力在随边处一致。利用Morino计算影响函数的解析公式,采用Yanagizawa方法求得物体表面上的速度分布,并对普通桨和大侧斜桨进行了数值预报。     

螺旋桨非定常性能预估的面元法

本文建立了一个预估螺旋桨非定常性能的基于速度势的非定常面元法,该方法中采用了双曲面元,避免了面元之间的缝隙,并采用迭代方法实现非线性的等压库塔条件。文中以第22届国际水池会议(ITTC)推进技术委员会推荐的HSP桨与D4679桨和挪威船舶技术研究所(MARINTEK)提供的PF-W-B桨为算例进行了考核计算,与试验结果和其他研究工作者的计算结果比较表明,螺旋桨的非定常性能预报的结果是令人满意的。     

拖式吊舱推进器的非定常水动力性能

应用升力面理论涡格法和面元法,建立了拖式吊舱推进器非定常水动力性能的数值计算方法。螺旋桨桨叶采用升力面理论涡格法计算,吊舱舱体及支架采用HESS-SMITH面元法计算,螺旋桨与吊舱及支架之间的相互影响通过迭代计算来处理。针对拖式吊舱推进器,通过系统的计算和分析,研究了螺旋桨负荷、吊舱和支架诱导速度各分量以及标称与实效诱导速度对其水动力性能的影响。研究表明,就吊舱及支架的实效诱导速度而言,轴向及周向诱导速度主要由支架引起,径向诱导速度主要由吊舱舱体引起。当考察吊舱推进器的定常水动力性能时,可略去吊舱诱导速度的径向及周向分量;考察非定常性能时,可略去径向分量,但应考虑周向分量的影响。以吊舱及支架的标称诱导速度作为进流,将导致非定常推力、扭矩的平均值降低,脉动量幅值减小,因此,虽然标称诱导速度容易得到,但据此进行吊舱推进器的性能预报或设计都会引起一定的误差。非定常水动力的脉动幅值取决于船尾伴流与吊舱诱导速度的相对比例,略去吊舱诱导速度会导致桨叶非定常力的变化特征发生较大变化。     

螺旋桨面元法数值库塔条件的改进

本文采用面元法和改进的为加塔条件预报螺旋桨定常性能。螺旋桨表面和它的尾涡面离散为四边形双曲面元法和每个面无上布置等强度衣子和源汇分布,通过对桨叶随边处非线性压库塔条件迭代方法的改进,建立了更为有铲的数值库塔条件的迭代方法,本方法与其它数值方法的计算结果和试验非常吻合,而且具有稳定的收敛特性     

导管调距桨的定常性能预估

本文建立了一个预估导管可调螺距螺旋桨水动力性能的数值计算方法，即螺旋浆用升力面理论、导管采用面元法，通过迭代计算考虑浆和导管的相互影响。引入了一个修正的螺旋浆尾涡模型，来模拟尾涡片的扭曲变形及分离现象，对导管桨性能预估的各种影响因素分析了系统研究，并考虑了桨毂对性能的影响。对JD简易导管桨和导管调距桨（JD7704导管＋JDC三叶可调螺距螺旋浆）分别进行了计算，并与实验结果进行了比较。结果表明，本文所建立的方法可较好地预估导管桨的水性能，精度比以往有较大的提高。     

基于面元法的吊舱推进器定常性能研究

基于势流理论面元法建立了吊舱推进器定常性能的计算方法.分别建立螺旋桨和吊舱的积分方程,通过在表面上布置双曲面元将方程离散为以面元上偶极强度为未知量的矩阵.螺旋桨和吊舱之间的相互影响通过迭代计算来处理.Newton-Raphson迭代过程被用来在桨叶随边满足压力Kutta条件.为避免数值求导中的奇异性,用柳泽(Yanagizawa)方法求得物体表面的速度分布.支架作为升力体处理,并通过迭代计算更新支架的尾涡形状.计算了拖式吊舱推进器的定常水动力性能,与实验结果的比较表明,计算误差在5%以内.分析了舱体对螺旋桨的影响,舱体的伴流会引起螺旋桨的载荷增大.     

Prediction of Hydrodynamic Performance of Marine Propellers by Surface Panel Method

The potential based low order surface panel method is used to predict the hydrodynamic performance of marine propellers. In present method the hyperboloidal quadrilateral panels are employed to avoid the gap between the panels. The influencecoefficients of panels are calculated by Morino's analytical formulations for increasing numerically calculating speed.The pres-sure Kutta condition is satisfied on the trailing edge of propeller blade by Newton-Raphson iterative procedure.Therefore the     

基于速度势迭代的面元法预报对转桨性能

采用基于诱导速度势迭代的面元法预报对转桨的定常水动力性能,在计算中,前桨和后桨的相互干扰通过影响系数实现;前桨尾涡面距后桨表面太近时的异变影响系数通过Lagrange插值求出。同时采用基于诱导速度迭代的面元法预报对转桨水动力性能,并与基于诱导速度势法进行对比,计算结果表明,采用诱导速度势迭代的方法计算的结果与试验值吻合良好。与基于诱导速度迭代的方法相比,计算结果精确,且可以节省计算时间。     

A method for numerical calculation of propeller hydrodynamics in unsteady inflow

The hydrodynamic performance of a propeller in unsteady inflow was calculated using the surface panel method. The surfaces of blades and hub were discreted by a number of hyperboloidal quadrilateral panels with constant source and doublet distribution. Each panel's corner coordinates were calculated by spline interpolation between the main parameter and the blade geometry of the propeller. The integral equation was derived using the Green Formula. The influence coefficient of the matrix was calculated by the Morino analytic formula. The tangential velocity distribution was calculated with the Yanagizawa method, and the pressure coefficient was calculated using the Bonuli equation. The pressure Kutta condition was satisfied at the trailing edge of the propeller blade using the Newton-Raphson iterative procedure, so as to make the pressure coefficients of the suction and pressure faces of the blade equal at the trailing edge. Calculated results for the propeller in steady inflow were taken as initialization values for the unsteady inflow calculation process. Calculations were carried out from the moment the propeller achieved steady rotation. At each time interval, a linear algebraic equation combined with Kutta condition was established on a key blade and solved numerically. Comparison between calculated results and experimental results indicates that this method is correct and effective.     

可调螺距螺旋桨水动力性能分析

利用面元法分析可调距螺旋桨的水动力性能.计算过程中,采用较为简捷的关于扰动速度势的基本积分微分方程,并采用双曲面形状的面元以消除面元间的的间隙,Newton-laphson迭代过程被用来在桨叶随边满足压力Kutta条件,使桨叶上下表面的的压力在随边有良好的一致性,同时用模拟物体真实行状的面元法来解决调距桨在螺距变化时的叶剖面畸变的问题.用Morino导出的解析计算公式来计算面元的影响系数,加快了数值计算的速度.以无厚度线性尾涡模拟桨叶泄出涡.调距螺旋桨最佳转轴位置由理论方法求出,使得桨叶的转叶矩为零.计算过程中计入了桨毂的影响,并分析了桨毂对桨叶表面压力分布的影响.最后给出了调矩螺旋桨水动力性能随随螺距的变化规律,并和试验结果作了比较分析.     

粘流方法的螺旋桨尾流场数值分析(英文)

The computational fluid dynamics(CFD) method is used to numerically simulate a propeller wake flow field in open water.A sub-domain hybrid mesh method was adopted in this paper.The computation domain was separated into two sub-domains,in which tetrahedral elements were used in the inner domain to match the complicated geometry of the propeller,while hexahedral elements were used in the outer domain.The mesh was locally refined on the propeller surface and near the wake flow field,and a size function was used to control the growth rate of the grid.Sections at different axial location were used to study the spatial evolution of the propeller wake in the region ranging from the disc to one propeller diameter(D) downstream.The numerical results show that the axial velocity fluctuates along the wake flow;radial velocity,which is closely related to vortices,attenuates strongly.The trailing vortices interact with the tip vortex at the blades’ trailing edge and then separate.The strength of the vortex shrinks rapidly,and the radius decreases 20% at one diameter downstream.