利用直接边界元法求解兴波问题

本文描述了利用直接边界元法求解定常兴波问题的理论和数值方法，本方法利用线性化自由面条件，辐射条件基于Dawson的四点差分格式进行处理。本文给出了Wigley船的兴波阻力和波高的数值结果，并与Dawson方法和试验结果进行了比较，结果显示本文的结果合理可靠。     

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.     

船体强度三维有限元分析的载荷直接计算和验证

采用载荷直接计算法对鞭船船体强度进行三维有限元分析，并用该船的实航测试数据和“船体梁理论”进行了验证和讨论。     

A time-domain higher-order boundary element method for 3D forward-speed radiation and diffraction problems

A time-domain higher-order boundary element method for seakeeping analyses in the framework of linear potential theory is newly developed. Ship waves generated by two modified Wigley models advancing at a constant forward speed in calm water or incident waves and the resultant radiation and diffraction forces are computed to validate this code. A rectangular computational domain moving with the same forward speed as the ship is introduced, in which an artificial damping beach is installed at an outer portion of the free surface except the downstream side for satisfying the radiation condition. The velocity potential on the ship hull and the normal velocity on the free surface are calculated directly by solving the boundary integral equation. An explicit time-marching scheme is employed for updating both kinematic and dynamic free-surface boundary conditions, with an embedding of a second-order upwind difference scheme for the derivative in the x-direction to stabilize the calculation. Extensive results including the exciting forces, added mass and damping coefficients, wave profiles, and wave patterns for blunt Wigley and slender Wigley hulls with forward speed are presented to validate the efficiency of the proposed 3D time-domain approach. The corresponding physical tests of the radiation and diffraction problems in a towing tank are also carried out. Computed numerical results show good agreement with the corresponding experimental data and other numerical solutions.     

Underwater acoustic field and pressure fluctuation on ship hull due to unsteady propeller sheet cavitation

The main objective of this paper is to develop an efficient numerical method which can predict the underwater acoustic field and pressure fluctuation on a ship hull due to unsteady propeller sheet cavitation by linear acoustic theory. In addition, the noise scattered from the ship hull and reflected from the free surface are included. Concerning the computation of the acoustic field induced by unsteady sheet cavitation and forces of a marine propeller, a method is derived without making any approximation about the distance function between the noise source and field point. Thus, this method can be used to predict acoustic pressure at both far and near fields, and this is very important for the scattering problem because the ship hull is located very close to the propeller. For the computation of the scattering problem, a more efficient and robust method is derived in time domain, which can treat multi-frequency waves scattered from underwater obstacles. The acoustic fields of a container ship radiated by the propeller and scattered from the ship hull with free surface is investigated in this paper. The pressure fluctuations of low blade rate on the ship hull induced by the propeller are also computed by the present method and are found to be similar to the results obtained by a panel method satisfying the Laplace equation for the points near the propeller due to the small retarding time. However, for the points on the ship hull away from the propeller, the differences of the results between two methods will increase.     

港口非线性波浪的实验研究

进行了港口非线性波浪对方箱形船体作用的三维实验研究,给出了在不同入射波条件下无因次点压力在船体表面上的分布,分析了压力随周期和波高的变化规律。结果表明:压力幅值在波浪周期较大时,出现一阶压力减小,高阶压力增大的强非线性现象,甚至还出现了二阶压力大于一阶压力的情况。在港口工程和大型浮体设计中对长周期波浪应注意考虑波浪的非线性特征。     

绕船体自由面周围三维粘性流场的数值模拟

本文采用有限体积法通过求解不可压缩的雷诺平均(RANS)方程数值模拟了包括兴波的三维船体周围的粘性流场,湍流模式使用了子网格尺度模式(SGS)和Baldwin-Lomax模式相结合的混合模式.对于自由表面的处理,采用了任意拉格朗日-欧拉方法,网格为不仅与物体表面贴体,而且与自由表面贴体的动网格,即随着自由表面的变化要不断地重新划分网格.虽然此方法需要很长的计算时间,但能较好地描述船体的兴波情况.本文计算了系列60船模在Fn=0.316,Re=1.9×106时带自由面的粘性流场,计算结果与试验结果吻合较好.，An finite - volume method solving of incompressible RANS Equations is developed for the numeri-cal simulation of three - dimensional viscous flow with free - surface about a ship. A hybrid turbulence mod-el by combining the sub - grid scale (SGS) model and the Baldwin- Lomax model is used in this pa-per. The arbitrary - Lagrange - Euler formulation is devised for the treatment of free surface boundary condi-tion. The boundary - fitted coordinate system is fitted not only to the ship hull surface but also to the freesurface, so the computational grid is regenerated to follow the free surface deformation. Although this methodneeds spend much longer computational time, but it can fairly display wave pattern made by a ship. The nu-merical simulation has been carried out for viscous flow with free surface past Series 60 ship hull at Fn =0. 316, Re = 1.9 × 106. The computed results are in reasonable agreement with experimental measure-ments.     

Experimental study on the nonlinear pressure acting on a high-speed vessel in irregular waves

The present study focuses on the nonlinear behavior of pressure on the hull surface of a high-speed vessel in irregular waves, particularly the pressure responses of alternately wet and dry areas near the waterline and on the bow zone. The vessel has high deadrise angles that may be subject to slight impact and water pile-up effects. A series of experiments in regular and irregular head waves were conducted, and the validity of applying Volterra modeling was investigated. In a previous article using experimental data in regular waves, it was confirmed that the approximate third-order Volterra model adequately simulated the variation of pressure responses in regular waves of different steepness up to a wave amplitude with a wavelength ratio of 0.01, even for the highly nonlinear pressures acting on the abovementioned areas of the hull surface. In this article, further validation for the second part of the study was obtained using experimental data in irregular waves. The frequency response functions obtained from the previous study’s experimental data in regular waves were applied to the third-order Volterra model by combining the input of irregular waves to simulate the responses in irregular waves of sea state five. Then, the spectra and statistics were analyzed. For the motions, accelerations, and pressure responses in irregular waves (as well as for the simulated time histories) variance spectra and statistics such as cumulative distributions of peak values and probability density functions were compared with the experimental results. It was confirmed that even for highly nonlinear and non-Gaussian pressures on the abovementioned areas of the hull surface, the approximate third-order Volterra model simulates the pressure responses in irregular head waves up to a sea state of five with adequate accuracy on deterministic and statistical bases.     

不同海况下艏部砰击及鞭状效应的试验与数值分析

为了更深入地研究船舶的鞭状效应,在拖曳水池中对某船进行了艏部砰击及鞭状效应分段模型试验研究。提出了一种可以考虑砰击力的非线性水弹性计算方法。并改进了传统的分段模型,采用变截面梁对船体刚度进行模拟以更好地接近实船。在规则波迎浪下观察到了严重的艏部砰击现象。试验数据表明,当波高从5.6m增大到21m时,由于鞭状效应原因,总弯矩相比低频波浪弯矩的增大值从24.64%增长到92.02%。最后,将不同海况下的测量结果与基于线性与非线性水弹性理论的计算结果进行了比较分析,初步验证了文中方法和程序在预报船体波浪载荷中的适用性。     

集装箱船整船有限元结构分析

文章以一艘１７００箱集装箱船为例,阐述了整船有限元结构分析方法。先建立全船有限元模型和质量模型,再用三维流体动力计算程序进行波浪随机载荷的长期预报,并在此基础上导出设计波参数组,最后,在全船有限元模型上 计算得出船体结构在各个设计波上的应力分布和变形结果,所得到的船体结构有限元分析结果对同类型集装箱船的设计和强度分析有一定的参考价值,对其它类型的船舶结构强度分析也有一定的借鉴意义。     

Comparison of numerical and experimental results of nonlinear wave-induced vertical ship motions and loads

 A nonlinear time-domain procedure is presented which is used to calculate the vertical responses of a container ship advancing in head waves. The method assumes linear radiation forces represented by time convolution of memory functions, infinite frequency added masses, and radiation restoring coefficients. The nonlinear hydrostatic restoring and Froude–Krilov forces are computed exactly over the instantaneous wetted surface of the ship's hull. Forces due to green water on deck are calculated using the momentum method. Nonlinear effects are identified on different vertical ship responses, namely on the heave and pitch motions, the vertical accelerations, and the vertical bending moment. These non-linear effects are expressed by the variation of the transfer function with the wave amplitude, the higher-order harmonics of the time signals, the offset of the time series, and the asymmetry of the peaks. The numerical results and the quantified nonlinear effects are compared with experimental results showing an ability to reproduce the main nonlinear effects. Received: December 17, 2001 / Accepted: January 31, 2002     

船体在波浪中的非线性水动压力

本文提出了一个预报船体在波浪中大幅运动时非线性水动压力场的二维时域理论。船体扰动势用时域自由面格林函数和在入射波下的瞬时湿表面上的分布源求解；与非线性水动压力场相匹配的船体运动用差分法求得。为提高计算效率和避免数值过程发散，采用了改进的数值模型和方案。通过线性理论计算与模型试验结果的比较，指出了线性切片理论在预报水动压力场时的不足，水动压力与波高的非线性关系及正负水动压力沿船体表面的分布在Wigley船的计算比较中得到了说明。初步计算表明，该理论的实用化发展前景是令有鼓舞的。相应的计算机程序可在PC机上运行。     

水面舰船迎浪航行时大幅运动预报的切片算法

将预报船舶运动和波浪载荷的切片理论加以扩展，应用于有限深水中船舶迎浪航行时大幅纵向运动和波浪载荷的时域求解。预报结果体现出考虑了湿表面及化后引起的船舶运动和受力的非线性性质。本计算方法简便、实用，适用于船舶纵向大幅运动预报。     

采用Green-Naghdi方程和有限元法计算浅水船波

采用波动方程／有限元法求解Green-Naghdi(G-N)方程计算船舶在有限水深区域的兴波和波浪阻力。把行驶船舶对水面的扰动作为移动压力直接加在Green-Naghdi方程里，以描述运动船体和水面的相互作用，并经此来计算不面波动、船底水动压力和波浪阻力。G－N方程比浅水方程增加一个非线性的频散项，以补充有限水深对浅水船波的影响。采用随船运动网格的有限方法，以Series 60 CB=0.6船作为算例给出浅水船波的计算结果，并与浅水方程的结果进行了比较。计算结果表明，当船速小于临界速度时，由于频散的影响，G－N方程级出的船后尾波波高比浅水方程的结果大，同时波浪阻力也比浅水方程的结果有所提高。当船速大于临界速度时，G－N方程的计算结果与浅水方程基本相同，频率散射无明显影响。     

船桨干扰定常空化性能数值模拟

文章采用CFD技术计算了船桨干扰定常空化性能。首先计算了不同空化数下的NACA66(MOD)水翼和DTRC4381螺旋桨的定常空化性能,同时预报了KCS船的伴流场、阻力性能和船体表面压力分布。为了验证计算船后螺旋桨空化流的可行性,文中对KCS船和KP505螺旋桨进行了整体考虑,并计算了不同空化数下的船后螺旋桨的空化覆盖面积。计算表明,文中采用的计算方法合理,计算结果与试验结果吻合。     

数值波浪水池及顶浪中船舶水动力计算

基于粘流理论建立了三维数值波浪水池,模拟了非线性波浪,并对规则波顶浪中前进的拘束船模的水动力进行了计算.数值模拟中,控制方程-RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波.给出了非线性规则波的模拟结果以及规则波顶浪中前进的拘束船模的水动力计算结果,并与理论解及DUT(Delfi University of Technology)的试验数据进行了比较,二者吻合良好.     

大开口船波浪载荷长期预报和弯扭强度整船有限元分析

大开口船全船弯扭联合变形与应力的精确计算，必须在整船结构模型上完成。本文以一艘5万吨级大开口船为例，用三维流体动力计算程序进行了波浪随机载荷的长期预报，并在此基础上导出设计波参数组，进而在全船整体结构有限元模型上计算了船体结构在各设计波上的应力分布，并采用嵌入精细舱口角区有限元网络的方法，在整船分析的同时计算出舱口角的应力集中值，获得了船体结构强度的详尽信息。文中阐述了波浪载荷的特点，设计波的确定，浮体完整结构计算的惯性平衡及大开口船的全船计算方法。     

Computational fluid dynamics (CFD) prediction of bank effects including verification and validation

Restricted waters impose significant effects on ship navigation. In particular, with the presence of a side bank in the vicinity of the hull, the flow is greatly complicated. Additional hydrodynamic forces and moments act on the hull, thus changing the ship's maneuverability. In this paper, computational fluid dynamics methods are utilized for investigating the bank effects on a tanker hull. The tanker moves straight ahead at a low speed in two canals, characterized by surface piercing and sloping banks. For varying water depth and ship-to-bank distance, the sinkage and trim, as well as the viscous hydrodynamic forces on the hull, are predicted by a steady state Reynolds averaged Navier–Stokes solver with the double model approximation to simulate the flat free surface. A potential flow method is also applied to evaluate the effect of waves and viscosity on the solutions. The focus is placed on verification and validation based on a grid convergence study and comparisons with experimental data. There is also an exploration of the modeling errors in the numerical method.     

绕船体粘性自由面流动的数值计算

考虑自由面的船舶粘性流动计算由于其复杂性和对计算机硬伯的较高要求，只是到了最近几年才得到较大发展，而在国内尚属起步阶段。本文使用VOF方法计算了绕Wigley船型的粘性自由面流动，流模式选择了SSTκ-ω模型。计算获得了船体上波形以及船波等值线图，同时考察了三个不同横截面上的流动状况。计算得到的阻力还和发表的试验结果作了对比，结果令人鼓舞。     

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

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