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.     

基于B样条的三维船体水动力数值计算

本文采用基于B样条的一种新的数值方法计算三维船体水动力，用B样条函数表达三维船体表面的几何形状以及流场中未知物理量的分布，为了验证该数值方法的可行性和精确度首先对处于无限流体域中的圆球体绕流问题进行了计算；其次计算了由ITTC所推荐的Wigley船型的兴波阻力以及以一攻角斜航时的操纵水动力；最后在一些假设下对两船作平行航行时的干扰水动力作了相应计算工作。数值计算结果与其它试验或理论结果在定量或定性上吻合良好。     

潜艇操纵性水动力数值计算中湍流模式的比较与运用

湍流模式的选取对潜艇操纵性水动力数值计算精度有着重要影响,采用六种湍流模式计算了SUBOFF主艇体及主艇体加指挥室围壳两种模型在一定漂角范围内的潜艇操纵性水动力,并与试验值进行了比较,结果表明SSTk-ω模型更为适合潜艇操纵性水动力计算。在此基础上对SUBOFF全附体模型在一定漂角和攻角范围内的艇体水动力进行预报计算,并对该计算方法应用于潜艇操纵性水动力预报计算的计算精度与适用范围进行了探讨。     

Numerical study on the hydrodynamic forces on a ship berthing to quay by taking free-surface effect into account

In the present study, numerical simulation of the berthing maneuver of a ship in the prescribed translational motion is performed. The transient viscous flow and hydrodynamic forces on the hull are calculated by solving the unsteady Reynolds-Averaged Navier–Stokes equations in overset grid system, and the free surface is captured using volume-of-fluid (VOF) approach. The present numerical results have been compared with previous computational results by Toda and available experimental data respectively. Since the effects of the quaywall and free surface are taken into consideration in the present study, it is found that the agreement is significantly better than that resulting from Toda’s 3D CFD based approach. Then the effects of various standoff distances between the ship and quaywall on the lateral forces are investigated. Meanwhile, the detailed transient flow features around the berthing ship are obtained, which are helpful to understand the interactional effects between the ship and quaywall. The present results may provide helpful guidance for vessels’ safe maneuvering in berthing motion and the design of fender system in the quay.     

Time domain simulations of the wave-induced motions of ships in maneuvering condition

The wave-induced motions of ships in maneuvering condition are numerically studied based on potential theory. The total disturbance potential is decomposed into a basic part and a perturbation part. The basic flow is evaluated based on the double-body model with a trailing vortex sheet. The perturbation flow is solved by using a time domain Rankine panel method to determine the hydrodynamic forces, and the wave-induced ship motions are then evaluated by an Adam–Moulton scheme. The solving process of the wave-induced motion is integrated with the maneuvering prediction by using a two–time scale model. Numerical tests are firstly carried out for a Series 60 ship, and the numerical results are compared with the experimental data to validate the numerical method for the basic flow. Then the wave-induced motions of the S-175 container ship in straight course and in turning condition are simulated; the numerical results are compared with the ITTC data and the experimental data, which show fairly good agreements.     

Locomotion by mechanical pectoral fins

The objectives of this study were the development of a new device for maneuvering an underwater vehicle using the mechanism of a fish swimming, an experimental and theoretical analysis of the hydrodynamic characteristics of the device, and its application to maneuvering a fish robot. Observations and experimental analysis of the pectoral fins of a black bass (Micropterus salmoides) revealed that the locomotion of the fish, such as swimming forward at low speed, swimming backward, and turning in a horizontal plane is generated by using a combination of a feathering motion and a lead-lag motion of the pectoral fins. A mechanical pectoral fin making a feathering motion and a lead-lag motion generates a thrust force in a range of phase differences between both motions. The unsteady vortex lattice method, including the effect of viscosity, can express fairly well the unsteady forces acting on the mechanical pectoral fin in the range of phase differences where it exerts the thrust force. The fish robot, consisting of a model fish body and a pair of mechanical pectoral fins, can not only swim forward and turn in almost the same position, but can also swim in a lateral direction without changing the yaw angle. Translation of an article that appeared in the Journal of The Society of Naval Architects of Japan, vol. 182 (1997): The original article won the SNAJ prize, which is awarded annually to the best papers selected from the SNAJ Journal, JMST, or other quality journals in the field of naval architecture and ocean engineering.     

URANS simulations of static and dynamic maneuvering for surface combatant: part 1. Verification and validation for forces, moment, and hydrodynamic derivatives

Part 1 of this two-part paper presents the verification and validation results of forces and moment coefficients, hydrodynamic derivatives, and reconstructions of forces and moment coefficients from resultant hydrodynamic derivatives for a surface combatant Model 5415 bare hull under static and dynamic planar motion mechanism simulations. Unsteady Reynolds averaged Navier–Stokes (URANS) computations are carried out by a general purpose URANS/detached eddy simulation research code CFDShip-Iowa Ver. 4. The objective of this research is to investigate the capability of the code in regards to the computational fluid dynamics based maneuvering prediction method. In the current study, the ship is subjected to static drift, steady turn, pure sway, pure yaw, and combined yaw and drift motions at Froude number 0.28. The results are analyzed in view of: (1) the verification for iterative, grid, and time-step convergence along with assessment of overall numerical uncertainty; and (2) validations for forces and moment coefficients, hydrodynamic derivatives, and reconstruction of forces and moment coefficients from resultant hydrodynamic derivatives together with the available experimental data. Part 2 provides the validation for flow features with the experimental data as well as investigations for flow physics, e.g., flow separation, three dimensional vortical structure, and reconstructed local flows.     

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

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

Circular motion tests and uncertainty analysis for ship maneuverability

Circular motion test data and uncertainty analysis results of investigations of the hydrodynamic characteristics of ship maneuvering are presented. The model ships used were a container ship and two tankers, and the measured items were the surge and sway forces, yaw moment, propeller thrust, rudder normal and tangential forces, pitch and roll angles, and heave. The test parameters were the oblique angle and yaw rate for the conditions of a hull with a rudder and propeller in which the rudder angle was set to zero and the propeller speed was set to the model self-propulsion conditions. Carriage data showing the accuracy of the towing conditions in the circular motion test are also presented. It was confirmed that the uncertainties in the hydrodynamic forces such as the surge and sway forces, yaw moment, rudder tangential and normal forces, and propeller thrust were fairly small. The reported uncertainty analysis results of the circular motion test data may be beneficial in validating data quality and in discussing reliability for simulation of ship maneuvering performance.     

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

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

Prediction of ship steering capabilities with a fully nonlinear ship motion model. Part 1: Maneuvering in calm water

This paper introduces a new method for the prediction of ship maneuvering capabilities. The new method is added to a nonlinear six-degrees-of-freedom ship motion model named the digital, self-consistent ship experimental laboratory (DiSSEL). Based on the first principles of physics, when the ship is steered, the additional surge and sway forces and the yaw moment from the deflected rudder are computed. The rudder forces and moments are computed using rudder parameters such as the rudder area and the local flow velocity at the rudder, which includes contributions from the ship velocity and the propeller slipstream. The rudder forces and moments are added to the forces and moments on the hull, which are used to predict the motion of the ship in DiSSEL. The resulting motions of the ship influence the inflow into the rudder and thereby influence the force and moment on the rudder at each time step. The roll moment and resulting heel angle on the ship as it maneuvers are also predicted. Calm water turning circle predictions are presented and correlated with model test data for NSWCCD model 5514, a pre-contract DDG-51 hull form. Good correlations are shown for both the turning circle track and the heel angle of the model during the turn. The prediction for a ship maneuvering in incident waves will be presented in Part 2. DiSSEL can be applied for any arbitrary hull geometry. No empirical parameterization is used, except for the influence of the propeller slipstream on the rudder, which is included using a flow acceleration factor.     

不对称流动对航道中航行船舶的影响

本文利用能够处理不对称自由面流动的改进的Dawson's方法来研究当船偏离航道航行时航道岸壁及水深对船舶所受的水动力的影响.计算结果与模型试验进行了比较,结果表明当航道水深不是很浅时,本文提出的方法和计算程序可以给出比较满意的结果,同时还可以发现当船舶偏离中心航道航行时,船两侧湿表面压力以及水面变化也是不对称的,水深以及岸边距离都能影响船体所受水动力及船波.     

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

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

基于NURBS的船体与桥墩间水动干扰力及力矩的非定常计算

本文在一些假设下完成了船体与桥墩之间水动干扰力及力矩的数值计算工作，所采用的方法是基于非均匀有理B样条(NURBS)的三维Rankine源边界元法，即将船体、桥墩的几何形状以及流场中待求的物理量分布完全表达成NURBS的形式。给出了船体通过圆柱形桥墩时的水动干扰力及力矩的数值计算结果，并进行了一些分析。     

舰船呆木设计及对操纵性影响研究

针对排水型舰船尾部线型特征,提出船体尾部呆木线型设计的数学函数,定义了呆木线型的特征量,并设计了系列呆木线型;运用有效展舷比概念扩充了船体水动力经验公式,得到呆木线型对船体水动力影响的估算方法;采用"分离型"操纵性数学模型预报了系列呆木对舰船操纵性的影响,得到呆木面积与航向稳定性提高成正比等结论。     

船舶操纵水动力导数的数值求解及敏感度分析

[目的]为了兼顾船舶操纵运动预报的成本与精度,基于数值计算方法,结合水动力导数敏感度分析,提出一种船舶操纵运动预报方法。[方法]首先,求解RANS方程,应用流体体积（VOF）法捕捉自由液面,采用动态网格方法对DTMB 5 415船型进行约束运动的数值计算,并将回归得到的线性水动力导数与试验值进行对比,验证数值方案的有效性;然后,基于MMG分离建模方法建立DTMB 5 415船模的操纵数学模型,并利用龙格-库塔算法进行求解,对船舶回转和Z形操纵运动进行仿真;最后,分析水动力导数对这2种操纵运动的敏感度。[结果]结果显示,采用所提方法得到的操纵轨迹和衡准参数仿真结果与试验结果一致,回转运动参数的平均误差为5.1%,Z形操纵运动参数的平均误差为11.7%,较文献使用CFD进行自航船模模拟得到运动参数的精度与计算成本均有所改善;水动力敏感度分析结果也验证了部分非线性水动力导数对操纵性衡准的影响较小,可以使用经验公式进行估算。[结论]研究表明采用所提方法进行船舶操纵性预报方法可行,可在满足工程应用精度的同时大大减少计算成本,尤其适用于船舶设计阶段的操纵性预报及优化。     

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.     

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

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

浸没在水中的海洋运载器数值模拟研究（英文）

This article presents a numerical study of the forces induced by hydrodynamic impact, that is, the impact of a part of the bottom of the hull on the water surface. The prediction of these efforts is often based on numerical simulations to determine the shock intensity of a structure on the surface of a weakly compressible fluid(for example, water). The short duration of the impact is also investigated in this work. This phenomenon occurs especially when a ship encounters a harsh and difficult sea conditions. Under such conditions, it is important to know how to predict the hydrodynamic forces applied to the structure to correctly optimize the ship elements during its design stage or to prevent possible damage. Indeed, various factors such as speed of the ship and height of the swell can cause the hull to partially emerge and then fall violently onto the water surface, which is referred to by naval personnel as tossing or slamming causing vibrations, stresses, and fatigue to the structural elements of the ship. In this work, we present an example of phenomenon modeling and then a numerical study of the different geometries(dihedron) that play a role in different sections of the bow. Then, we compare our present results with the theoretical and experimental results of other researchers in the field. The average interval impact time for a dihedral model corresponding to the section of the chosen ship and other experimental and theoretical data is in good agreement with the experimental and theoretical measurements.     

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

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