船舶多自由度斜航运动水动力数值研究

计及多自由度运动的船舶斜航水动力预报对船舶航行安全具有重要意义。文章通过耦合求解船舶运动方程和雷诺平均N-S方程,并采用VOF方法和高精度自由面捕捉技术对作多自由度斜航运动船舶的粘性绕流场进行数值模拟。船舶动态平衡位置根据计算出的力和力矩来决定,得到包括升沉、纵倾和横倾在内的船舶浮态。文中采用的算例与爱荷华大学进行的模型试验相同,通过比较数值计算结果和试验值验证了该方法的有效性。对船模在受约束和自由运动两种状态下的船舶运动和流场进行模拟,通过比较分析船舶升沉、纵倾和横倾的影响。文中计算获得的详细流场细节特征,包括前体和舭部的涡以及船体表面上的压力,有助于理解船舶斜航运动浮态变化的机理。     

浅水中低速斜航船舶水动力预报及验证与确认分析

浅水中的斜航船舶受到浅水阻塞效应和不对称流的综合影响。为预报该运动中的船舶水动力,文章采用基于定常雷诺平均纳维—斯托克斯方程的计算流体动力学方法,对浅水中做斜航运动的船舶粘性绕流场进行数值模拟。考虑低航速运动的特点,忽略航速影响下的自由面兴波,由数值计算得到水动力系数在漂角影响下的变化规律。针对计算精度问题,在数值模拟中从验证和确认角度分析和评估计算结果：通过网格收敛性分析分析数值误差与不确定度;结合试验数据考察计算模型的误差。此外,从计算区域尺度、湍流模型、边界条件、船体下沉和纵倾作用方面对模型误差的影响因素进行探讨,可为改进计算模型、提高数值模拟精度提供参考依据。     

Time-domain TEBEM method for mean drift force and moment of ships with forward speed under the oblique seas

A numerical method for solving 3D unsteady potential flow problem of ship advancing in waves is put forward. The flow field is divided into an inner and an outer domain by introducing an artificial matching surface. The inner domain is surrounded by ship wetted surface and matching surface as well as part of the free surface. The free surface condition for the inner domain is formulated by perturbation about the double-body flow or uniform incoming flow assumption. The outer domain is surrounded by matching surface and the rest free surface as well as infinite far-field radiation boundary. The free surface condition for the outer domain is formulated by perturbation about uniform incoming flow. The simple Green function and transient free surface Green function are used to form the boundary integral equation (BIE) for the inner and outer domains, respectively. Taylor Expansion Boundary Element Method (TEBEM) is utilized to solve the double-body flow and inner domain and outer domain unsteady flow BIE. Matching conditions for the inner domain flow and outer domain flow are enforced by the continuity of velocity potential and normal velocity on the matching surface. Direct pressure integration on ship wetted surface is used to obtain the first-order and second-order wave forces (moments). The numerical predictions on the displacement, added resistance, sway mean drift force and yaw mean drift moment of the modified KVLCC2 ship at different forward speeds are investigated by the proposed TEBEM method. It is also compared with the other numerical results. The physical tank experiment results are also developed to validate the accuracy of numerical tank results. Compared with the experiment solutions, a good agreement can be obtained by TEBEM method.

     

LNG船舶靠泊对护舷的撞击力试验分析

为确保LNG船舶靠泊安全,对浙江LNG接收站工程进行了船舶靠泊物理模型试验。采用牵引船模及直流力矩电机调节速度的方法,模拟了船舶不同靠泊速度和不同偏心距的靠泊方式。试验表明,船舶靠泊角度为5°时,4个护舷受力不均,1#护舷先受力且受力最大;不同偏心距条件下的撞击力没有明显规律,且相差不大;靠泊速度为0.15 m/s条件下,静水靠泊和顶流靠泊时设计护舷型号可满足要求,横浪1.5m时靠泊则撞击力不满足要求。建议尽量减小靠泊角度,使得护舷受力均匀,以减少靠泊时的撞击力。     

Study of damaged ship motions taking into account floodwater dynamics

This paper presents an application of shallow water theory to describe the motion of floodwater inside a rolling ship in damage condition. The time domain theoretical approach to the coupled problems of ship and water inside compartment motions is briefly described, including the method used to solve for the water motion characteristics and forces exerted on the ship. This approach is applied to the study of the behaviour of a passenger Ro–Ro ship in regular beam seas and numerical results are given for the intact and damaged conditions. Comparison is made with experimental results. For the damaged condition, the characteristics of the floodwater motion are studied in the time domain for a number of different wave frequencies. The shape of the free-surface and phase of water motion in relation to the ship roll motion are shown for several wave frequencies. The dynamic floodwater roll moment is also shown and compared with the static roll moment (flat horizontal free surface), allowing the conclusion that the dynamic roll moment is much larger than the static roll moment, for high wave frequencies, and is in phase opposition in relation to the roll motion.     

基于OpenFOAM的船舶与液舱流体晃荡在波浪中时域耦合运动的数值模拟

波浪中载液船舶运动激励舱内液体的晃荡,舱内液体晃荡产生的冲击力同时作用在舱壁上,进而影响船舶的运动姿态。波浪中船体水动力和时延函数是在势流理论范畴下采用切片法和脉冲响应函数方法计算获得的,液舱内液体非线性晃荡是基于粘性流理论实时计算模拟,两者耦合建立了波浪中载液船舶与液舱流体晃荡耦合的运动方程。论文基于开源CFD开发平台OpenFOAM,自主开发实现了船体运动与液舱晃荡的耦合计算程序,并进行了相应的数值模拟计算和验证工作。该方法完整地考虑了波浪、船体和液舱晃荡之间的耦合作用,并结合船体内外流场特点分别采用了势流和粘性流理论,具有较高的计算效率。通过数值模拟计算和模型实验研究表明,数值模拟计算能够清晰显现出液舱晃荡对船体全局运动影响,船体运动计算结果与模型实验结果吻合良好。     

波浪中母船干扰下潜水员梯的运动及水动力特性分析

潜水员梯是通过销轴连接在母船舷侧为潜水员执行潜水作业时上下船的平台,为保证在波浪中使用时的安全性,对其波浪作用力、运动响应及销轴处的动作用载荷设计有较高要求。为研究母船干扰下潜水员梯的动态响应特性,采用基于深水格林函数为积分内核的边界元方法,先计算母船及潜水梯这一整体的运动响应及其平均湿表面的辐射和绕射压力分布,依据刚体运动合成原理求解得到潜水梯重心处的运动。再以潜水员梯为研究对象,采用压力积分方法获得潜水员梯的波浪力作用,并将其运动、波浪力和销轴处的作用力联立构成动力学平衡方程,逆向求解出销轴处的动作用力。基于上述思路,系统研究了波长、航速及潜水员梯与母船相对位置变化时潜水员梯的运动、波浪作用力及销轴处作用力的变化规律,为潜水员梯的布置、销轴的选型及潜水作业时母船的操纵可提供指导意义。     

进水口形式对闸墙长廊道船闸输水系统影响的数值模拟

船闸进水口水流条件关系到输水系统的输水效率、船舶停泊安全,因此进水口的形式布置是船闸水力设计和模型试验的重要内容之一。针对两种常见的进水口布置形式,结合某船闸输水系统的具体条件,采用数值模拟的方法,从输水水力特性、水流流态方面对比分析了这两种进水口布置形式对于输水系统的影响,研究结果表明：对于闸墙长廊道船闸输水系统,相对于进水口外侧导墙拉直,进水口外侧导墙斜拉的工程措施在控制惯性超高和闸室水面最大上升速度、增强进水支孔流量过程稳定性、改善进水口水流流态方面更具优势。研究成果可为船闸输水系统进水口布置提供技术参考。     

A finite pointset method for the numerical simulation of free surface flow around a ship

In this paper, we present a Finite pointset method (FPM) for the numerical simulation of free surface flow around a ship in calm water. It is a Lagrangian and meshless particle scheme which is applied to the projection method for the incompressible governing equations. This requires the solution of Poisson problems in each time step, so a moving least squares (MLS) interpolants is used for the spatial derivatives in order to discretize the Poisson equation with pressure-Dirichlet condition of free surface flow in meshless structure. Meanwhile, an additional problem of the periodic particle locations redistribution in the present approach is still handled by MLS interpolants. With the proposed FPM technique, problems associated with the free surface flow around a ship are circumvented. A verification of numerical modeling is made using the Wigley hull and the validity of the proposed methodology is examined by comparing the detail of wave profile and wave-making resistance with Series 60 model. The results demonstrate that FPM is able to perform efficient and stable simulations of free surface flow around a ship.     

CFD Analysis of Ship-to-Ship Hydrodynamic Interaction

A numerical study of ship-to-ship interaction forces is performed using a commercial CFD code, and the results are compared with experimental data and with the results of a panel method analysis. Two ship models have been used in the interaction forces analysis: a tug and a tanker, advancing parallel to each other with different lateral distances and two different values of the fluid depth. Computations are carried out with four different flow models: inviscid and viscous flow with the free surface modeled as a rigid wall and inviscid and viscous flow with the deformable free surface. A fair agreement was obtained with available experimental data and results obtained by panel method. The influence of viscosity in the computations is found to be comparatively weak, while the wavemaking effects may be important, at small magnitude of the horizontal clearance.     

船船水动力相互作用的CFD分析研究（英文）

A numerical study of ship-to-ship interaction forces is performed using a commercial CFD code, and the results are compared with experimental data and with the results of a panel method analysis. Two ship models have been used in the interaction forces analysis: a tug and a tanker, advancing parallel to each other with different lateral distances and two different values of the fluid depth. Computations are carried out with four different flow models: inviscid and viscous flow with the free surface modeled as a rigid wall and inviscid and viscous flow with the deformable free surface. A fair agreement was obtained with available experimental data and results obtained by panel method. The influence of viscosity in the computations is found to be comparatively weak, while the wavemaking effects may be important, at small magnitude of the horizontal clearance.     

面向自主靠泊的船舶航迹估计方法

针对无人船在靠泊过程中自身运动估计和周围环境重建的需求,提出一种基于水面区域去除的船舶航迹估计方法。分析船舶靠泊场景中水域图像的特征,提出水岸线检测方法,将该水岸线检测方法与视觉里程计(Visual Odometry,VO)方法相结合,提出一种泊位水域环境中的船舶航迹估计方法,采用实船试验进行验证。现场试验表明:该方法在自主靠泊场景中对船舶航迹进行估计是可行的,试验结果与真实数据(Real Time Kinematic,RTK)值对比表明:在忽略尺度的情况下,算法估计轨迹的精度可满足靠泊要求。     

Finite-volume simulation of flows about a ship in maneuvering motion

A finite-volume method of computing the viscous flow field about a ship in maneuvering motion was developed. The time-dependent Navier-Stokes equation discretized in the generalized boundary-fitted curvilinear coordinate system is solved numerically. A third-order upwind differencing scheme, a marker and cell (MAC)-type explicit time marching solution algorithm and a simplified subgrid scale (SGS) turbulence model are adopted. The simulation method is formulated, including the movement of a computational grid fitted to the body boundary that allows computation of the flow field around a body under unsteady motion. To estimate the maneuvering ability of a ship, the accurate prediction of the hydrodynamic forces and moments of the hull is important. Therefore, experimental methods of finding the hydrodynamic forces of a ship in maneuvering motion, such as the oblique towing test, the circular motion test (CMT) and planar motion mechanism (PMM) test, were established. Numerical simulation methods for those captive model experiments were developed introducing computational fluid dynamics (CFD). First, numerical methods for steady oblique tow and steady turn simulation were developed and then extended to unsteady forced motion. Simulations were conducted about several realistic hulls, and the results were verified by comparisons with measured results obtained in model experiments. Hydrodynamic forces and the moment, the longitudinal distribution of the hydrodynamic lateral force, and the pressure distribution on the hull surface showed good agreement.     

40万吨级超大型油轮船舶系泊动态响应的试验研究

为探讨超大型油轮在现有码头停靠、作业时的系泊情况,采用船舶系泊物理模型试验方法,分析研究了40万吨级油轮系泊作业时在长周期波及波浪、风、流联合作用下,系泊船舶的动态响应及对系缆力、护舷撞击力的影响。得出了40万吨级油轮在系泊作业过程中在不同波高及周期的波浪作用下运动量、缆力、撞击力的变化规律,提出在试验条件下40万吨级油轮停靠码头系泊作业的系缆方式及应注意问题。可为40万吨级油轮系泊作业提供依据。     

浅水条件下船舶通过船闸时的水动力性能数值研究

文章通过应用CFD方法数值模拟在浅水条件下通过船闸的船舶粘性绕流,对船舶通过船闸时的水动力性能进行了数值预报研究。通过UDF编程定义船舶的运动,使用动网格方法和滑移交界面技术进行船舶运动过程中的网格更新,计算作用在船体上的水动力,并由计算得到的水动力求得船体下沉和纵倾。为了验证所采用的数值方法,以一艘通过比利时泽布吕赫Pierre Vandamme船闸的船舶为例,在模型尺度下进行了计算,并将计算结果和佛兰德水利研究所的模型试验基准数据进行了比较。通过分析不同船速、偏心距和水深条件下的数值结果,给出了这些因素对船舶通过船闸时的水动力性能的影响。该文研究结果可为浅水条件下船舶通过船闸时的安全操纵和控制提供一定的指导。     

浅水中船舶水动力特性数值计算

本文对浅水中船舶水动力特性进行数值计算研究.采用RANS方程结合RNG K-ε两方程湍流模型,对一方形系数0.6的系列60船模在浅水中的阻力、升沉、纵倾和兴波进行数值计算,其中自由面采用VOF方法处理;计算中,水深Froude数范围0.6～1.8,包含了临界和超临界水深Froude数.数值计算得到的阻力、升沉和纵倾与模型试验结果及采用三维扩展Boussinesq方程的计算结果进行了比较分析,吻合较好,部分计算结果得到改进.     

船舶与浮式码头碰撞过程仿真

基于任意拉格朗日-欧拉法,描述了靠泊船与浮式码头碰撞过程中流固耦合关系。采用耦合法建立了浮式码头与靠泊船碰撞的仿真模型,真实模拟了船舶与浮式码头碰撞的过程,并与试验结果进行对比。由此可知,靠泊船碰撞对于浮式码头的运动响应、连接器所承受载荷、锚缆张力影响较为显著,流体区域也将产生较强的扰动,在实际工程应用中应加以注意。同时,也证明任意拉格朗日-欧拉法在流固耦合的研究领域优势明显。     

纵摇容器中液体晃荡的非线性数值模拟

在海洋工程领域,液体晃荡是一种普遍存在的物理现象。对于船舶而言,转动比平动有着更重要的影响。该文针对纵摇容器中的液体晃荡问题,采用高阶边界元方法建立自由水面满足完全非线性边界条件的时域数学模型。通过大地坐标系和随体坐标系之间的坐标变换,使得计算域仅控制在随体坐标系内。求解中采用半混合欧拉—拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过与已发表试验和数值结果的对比,验证了建立模型的准确性。进而开展大量数值试验研究容器纵摇运动频率、纵摇转动中心和容器中布置一竖直隔板对晃动波面与荷载的影响。     

月池共振对超深水钻井船运动性能的影响

为了研究月池内流体共振对船体运动的影响,以一艘超深水钻井船为研究对象,采用三维势流理论计算钻井船月池内自由液面和船体的运动,通过月池内自由液面的运动信息,得到月池各种模态下的共振频率,并以此来检验相关理论和经验公式对月池共振频率预报的准确性。月池底部打开和关闭2种状态下钻井船的运动特性显示,月池的活塞共振运动会显著增加船体的垂荡运动,而月池的晃荡共振运动对船体运动的影响较小。     

Numerical analysis on free surface waves and stern viscous flow of a ship model

This paper deals with numerical techniques for computing the viscous flow past a ship hull with and without a free surface using a Reynolds-averaged Navier-Stokes solver with global conservation. In the first technique, a coarse grid is used to find an approximate solution to the free surface problem. Interpolation of a fine grid is subsequently carried out, and a more exact solution, particularly in the boundary layer and wake, is obtained. In the second technique, a modified Baldwin-Lomax model is introduced to compute the viscous flow with and without a free surface. These numerical techniques are applied to simulations of the flow around a Series 60 and an SR196C ship model. The results are compared with measurement data, and the usefulness of the numerical techniques is demonstrated.