基于SHIPFLOW软件的方尾舰船阻力快速预报

SHIPFLOW软件被广泛应用于基于阻力的船型方案优选,但其对舰船阻力预报的适用性却少有研究。为确保舰船方案论证及优选结果的准确性,针对常规水面舰船方尾船型开展计算分析,提出适用于方尾舰船兴波阻力计算的网格划分方式以及获取兴波阻力系数的方法。为准确估算舰船阻力,可引入排挤厚度来修改船体型值,以在兴波阻力计算中等效计入粘性影响,并在此基础上形成舰船阻力快速预报方法。将某方尾舰船阻力数值计算结果与试验结果进行了对比,结果表明,预报误差基本控制在7%以内,验证了该方法的可行性。     

基于纵倾调整的集装箱船实海域航行阻力减少的研究

为应对"能效运营指标(EEOI)自愿应用导则",提出了在航行中采用调节纵倾的方法,降低营运能耗。选取了3100TEU集装箱船作为研究对象,应用Rankine源自由面势流理论计算不同纵倾下船舶兴波阻力系数以及总阻力系数,应用三维面元法计算了规则波中的波浪增阻,并针对实海域进行了増阻预报。研究表明,一定的艉倾能减小兴波阻力和总阻力。通过实海域増阻预报,服务航速下采用艉倾0.4°的纵倾航行会使波浪増阻增加4.21%,但总阻力仍能减少3.38%,较有效地降低了EEOI值。     

数值模拟现代水面舰船带自由面的湍流绕流场

为了研究带自由面的船舶湍流绕流场,选择了一艘ITTC推荐的公开船模作为模拟对象,数值求解RANS方程。此模型为带声呐导流罩和方形尾封板的复杂水面舰船。采用Ogrid块拓扑算法生成质量较高的纯六面体多块结构化网格,计算中采用RNG k-ε湍流模式和标准壁面函数,并使用VOF算法来捕捉自由面。空间离散采用QUICK格式,压力-速度解耦采用PISO算法。将阻力和波形的数值结果与实验数据相比较,总阻力系数误差约2%,船侧波形、船首自由面吻合良好,显示了CFD方法在船舶水动力学中预测带自由面湍流绕流场的有效性。     

基于CFD的船舶不同浪向下运动及增阻数值分析

为分析波浪对船舶快速性和耐波性的影响,必须对波浪中航行的船舶阻力增值进行准确预报。本文基于计算流体力学软件FINE/Marine建立了Wigley船模的数值模型,对不同规则波波长下的船体运动和波浪增阻进行了计算,并与试验结果进行对比,验证了数值模型的可行性与准确性。同时计算分析了船舶在规则波中航行时的波浪增阻与浪向之间的变化关系。通过研究发现：随着浪向角的增大船舶波浪增阻逐渐增加,在60°浪向角时波浪增阻达到最大值,浪向角对波浪增阻的影响较大。本文的研究方法可用于船舶有航速下的不同浪向波浪增阻的数值预报。     

粘性对方尾船波浪力的影响分析

航行船舶绕射问题的准确求解是预报船舶波浪中运动的基础.文章分别采用基于三维移动脉动源的边界元法和基于RANS方程的CFD方法求解了规则波中航行的方尾船模型的波浪力,并通过与试验结果比较验证了方法的可靠性.CFD计算时通过从波浪力中分离出粘性力,研究了粘性对方尾船波浪力的影响.结果表明,基于RANS方程的波浪力计算结果与试验结果和势流理论结果吻合良好,在长波中与试验值更接近.粘性对船模斜浪航行时的横摇力矩有较大影响.     

分段式尾压浪板对高速船阻力性能的影响

基于CFD技术,以排水型高速船Model 5b为模型,寻求改善高速船阻力性能的尾压浪板新形式。首先尾压浪板新形式的确定在静水条件下进行,然后在波浪条件下验证该压浪板的阻力性能。基于CFD软件建立三维数值波浪水池,静水条件下采用切割体网格技术预报船模的阻力性能。波浪条件下数值水池入口采用直接造波方法,尾部采用人工阻尼消波方法,自由面采用VOF方法处理,采用重叠网格技术预报船模的阻力性能以及运动响应。确定一种比常规压浪板阻力性能优良的分段式压浪板,为船舶节能附体的研究提供参考。     

Fundamental research on resistance reduction of surface combatants due to stern flaps

It is well known that a decrease in ship resistance may be achieved due to the installation of a stern flap. Therefore, so far, a considerable amount of research on stern flaps has been conducted. Previous research has demonstrated that the primary mechanism by which a stern appendage reduces resistance is a change in the pressure distribution over the aft body of the hull, and secondly through effects on the running attitude, near and far field wave generation, and local transom flow among other phenomena. However, the change in pressure distribution is influenced by the other components. Hence, there is still room for argument about the relative contribution of each component to the pressure distribution. Therefore, as the first step of the research, by conducting the model experiment in towing tank and CFD (Computational Fluid Dynamics) analysis, we examined the effect of running attitudes and wave making at the after portion of the hull on resistance reduction. As a result, it is concluded that a flap affects a change in the wave generated at the transom part and it could lead to a decrease in wave-making resistance.     

Large-amplitude motion responses of a Ro–Ro ship to regular oblique waves in intact and damaged conditions

This article presents a nonlinear time-domain simulation method for the prediction of large-amplitude motions of a Ro–Ro ship in regular oblique waves in an intact and a damaged condition. Numerical computations and model tests have been carried out to investigate the dynamic motion responses of Ro–Ro ship Dextra to various wave amplitudes at three different wave headings. The results of numerical and experimental investigations for stern quartering waves are reviewed. Comparisons between predictions and measurements show good agreement except in the roll-resonant region. Nonlinear effects are significant in horizontal modes of motion, and resonant roll motion, and there is strong coupling between all modes of motion in the roll-resonant region for large-amplitude responses. On the other hand, the time-domain simulation technique suffers from numerical drift in horizontal modes of motion as wave amplitude increases. This is due to nonlinear equations of motion and the lack of a restoring force and moment in horizontal motion. Received: April 30, 2002 / Accepted: August 9, 2002 Acknowledgments. II Programme of the European Community Commission under contract No. BRPR-CT97-0513. Address correspondence to: H.S. Chan (hoi-sang.chan@ncl.ac.uk)     

规则波和不规则波中船舶艉砰击及其振动响应的试验研究

在拖曳水池中对某船舶进行了艉砰击及其振动响应的试验研究.在规则波以及不规则波中的零航速、艉随浪情况下观察到了严重的艉砰击现象.试验数据分析表明,合成弯矩可以分成由波浪载荷引起的低频成分以及由砰击载荷引起的高频成分.由于严重艉砰击载荷的作用,发现在某次规则波试验中合成弯矩比波浪弯矩要大出44%,在3.24m不规则波中合成弯矩增加了43%.不规则波中的试验数据统计表明合成弯矩分布范围服从Weibull分布.推导了服从Weibull分布随机变量的短期概率极值预报公式,针对试验数据进行了预报.还讨论了试验数据分析中的不确定性问题.试验研究表明,对于艉部平坦肥大的船舶,在设计和使用中需要引起对艉砰击及其振动响应问题的重视.     

方尾船绕射问题RANS求解研究

基于RANS方程,采用VOF方法捕捉自由液面,通过自定义函数UDF在Fluent平台上建立了三维数值波浪水池,模拟了规则波的生成、传播和消波,发现波陡对数值造波的精度有较大影响。求解了规则波中顶浪航行的方尾船模型的波浪力,并与试验结果和势流理论结果进行了比较。结果表明,基于RANS方程的波浪力计算结果与试验结果和势流理论结果吻合良好,在长波中与试验值更接近。     

Evaluation of added resistance in regular incident waves by computational fluid dynamics motion simulation using an overlapping grid system

A computational fluid dynamics simulation method called WISDAM-X was developed to evaluate the added resistance of ships in waves. The Reynolds-averaged Navier–Stokes (RANS) equation was solved by the finite-volume method and a MAC-type solution algorithm. An overlapping grid system was employed to implement rigorous wave generation, the interactions of ships with incident waves, and the resultant ship motions. The motion of the ship is simultaneously solved by combining the solution of the motion of the ship with the solution of the flow about the ship. The free surface is captured by treatment by the density-function method. The accuracy of WISDAM-X is examined by a comparison with experimental data from a container carrier hull form, and shows a fairly good agreement with respect to ship motion and added resistance. Simulations were also conducted for a bow-form series of a medium-speed tanker to examine the effectiveness of the WISDAM-X method as a design tool for a hull form with a smaller resistance in waves. It was confirmed that the WISDAM-X method can evaluate the added resistance with sufficient relative accuracy and can be used as a design tool for ships.     

基于RANS法和边界层理论预报三维船体阻力

基于船体阻力成因和分类,应用二因次RANS法并结合EFD技术为研究手段,研究了船体各阻力成分和兴波波形,考虑边界层第一层网格高度对数值计算结果的影响,通过数值模拟计算和兴波波形比较,分析表明第一层网格无量纲高度y+值充分影响数值结果精度,总阻力在高傅汝德数下其数值误差有增大趋势,同一航速下摩擦阻力值随y+增加略微增大,剩余阻力值却明显减少,且自由面兴波随傅汝德数增大首部波高愈发明显,船中及尾部区域兴波向外扩散,船首尾肩部开尔文波系明显且波峰后移,其趋势与EFD波形一致,具有实用性。     

压浪板与游艇航态及阻力性能关系的研究

安装压浪板会在游艇艉部形成局部高压区,产生附加水动升力从而改变游艇的航行姿态,并影响游艇的阻力性能。采用半理论半经验方法分析安装压浪板前后游艇的姿态及阻力变化,研究压浪板对游艇产生影响的作用机理以及二者之间的匹配关系;获得压浪板最佳尺寸及安装角度的计算方法。结合方艉尖舭型游艇实例,验证压浪板在工程应用中的价值。     

On nonlinear ship waves and wave resistance calculations

A finite difference method based on the Euler equations is developed for computing ship waves and wave resistances. Time marching is carried out using a time-splitting fractional-step method. The second-order central difference is used to discretize the spatial differentials, while the convection terms are discretized by the QUICK scheme. A body- and free-surface-fitted grid system with a cell-centered stencil is used. A Poisson equation for the pressure increment at each time step is solved to enforce mass conservation. The method is validated by comparing the numerical results with the experimental data for a Wigley parabolic hull. The characteristics of ship waves, such as the wave profile along the hull, the wave pattern on the free surface, the pressure distribution on the hull surface, and the wave-making resistance are reasonably predicted. The calculated results are in good agreement with the experimental data. Received for publication on Oct. 13, 1998; accepted on Jan. 20, 1999     

沉管管节浮运过程中波浪附加阻力的水动力学分析

波浪引起的附加阻力是管节海上浮运过程中总阻力的重要组成。采用基于势流理论的水动力学软件AQWA计算了管节浮运过程中的波浪附加阻力。通过建立管节及附属结构的三维有限元模型,确定管节的几何、物理性质,建立湿表 面模型进行水动力学分析。对影响附加阻力的各因素,包括水深、航速、波浪谱及波浪要素等,进行了参数敏感性分析。计算结果表明,采用P-M谱的计算值较采用JONSWAP谱有所不同,附加阻力随水深增加而减小,阻力值与波浪周期及波浪入射角均有一定的关系,与波浪高度的平方成正比,且与拖航速度成正比。     

Numerical simulation of three-dimensional breaking waves

Three-dimensional (3D) wave breaking around bodies of complex geometry has been numerically investigated by use of two types of Navier-Stokes solvers, namely the finite-difference and the finite-volume methods employing rectangular and curvilinear coordinate systems, respectively. Both methods employ the density-function technique to capture the free surface location and can cope with complicated free surface configurations such as breaking waves. The accuracy of the density-function method is examined through the comparison with experimental results, and it is confirmed to be satisfactory when the grid spacing and the time increment are sufficiently small. New computational methods are applied to several problems including 3D breaking waves around ships and wave diffraction around offshore structures. The computed results show good agreement with experimental results indicating that wave breaking phenomena are successfully simulated. The qualitative accuracy, however, could be improved by including the dissipating effect of breaking waves.     

三体船兴波问题的数值计算

以基于NURBS(non-uniformra tionalB-splines,非均匀有理B样条)的广义高阶面元法讨论了三体船兴波波形与兴波阻力的数值计算问题。对两种不同船型的三体船进行了计算,计算结果与试验结果的比较表明了该方法的有效性。文中还讨论了高速方尾三体船侧体的横向偏距和纵向偏距对兴波阻力的影响。     

船舶在波浪中阻力增加预报研究进展

从数值计算和模型试验两个方面,介绍船舶在波浪中阻力增加预报的研究进展,分析各种方法的基本原理,提出阻力增加预报的发展方向.     

短波顶浪中KVLCC2船模阻力增加CFD计算研究

波浪中航行船舶阻力增加,特别是短波中的阻力增加,是船舶界关注的焦点之一,也是船舶水动力学界研究的热点之一。论文采用基于RANSE的数值波浪水池技术,针对KVLCC2船型,开展了短波顶浪中船舶阻力增加的数值计算研究。与模型试验结果的比较表明,文中的CFD方法能够相当准确地计算短波顶浪中航行船舶的阻力增加;对船体各部分波浪增阻的分析表明,船体艏段产生的波浪增阻占主导地位,艉段的波浪增阻很小,而平行中体段对波浪增阻几乎无贡献。     

Global wave loads on intact and damaged Ro-Ro ships in regular oblique waves

A nonlinear time-domain simulation method is presented for the prediction of dynamic global wave loads on a Ro-Ro ship at zero speed in regular oblique waves in an intact and a damaged condition. Numerical computations and model tests have been carried to investigate the structural responses of Ro-Ro ship Dextra to various wave amplitudes at three different wave headings (DTR-4.1-NEW-12.98, DEXTREMEL project BE97-4375, 1998; DTR-4.2-NEW-11.99, DEXTREMEL project BE97-4375, 1999). The results of numerical and experimental investigations for stern quartering waves are reviewed. Comparisons between predictions and measurements for global wave loads at the midship section of the intact and the damaged Ro-Ro ship show that the agreement between the theory and experiment for dynamic horizontal and vertical bending moments is excellent. On the other hand, correlation between the predictions and measurements for dynamic vertical shear force is better than that for dynamic horizontal shear force. Nevertheless, the calculated torsion moment values are higher than the measured values. As the wave amplitude is not small, the positive and negative peaks of global wave loads are no longer equal to each other as found in both the calculations and experiments. The dynamic vertical global wave loads in the damaged condition are larger than that in the intact condition.