CFD simulation of 3-dimensional motion of a ship in waves: application to an advancing ship in regular heading waves

A new computational fluid dynamics simulation method has been developed for the unsteady motion of a ship advancing in waves. The objective is to evaluate the added resistance and predict the performance of a ship in waves. In this study, a finite volume method, in the framework of a boundary-fitted grid system, is employed. The motion of the ship is solved with six degrees of freedom by using the hydrodynamic forces and moments obtained from the solution of the simulation method. The marker–density–function method is employed to calculate the nonlinear free surface. This method is applied to the coupled motion problem of heaving and pitching. Received for publication on Nov. 15, 1999; accepted on Nov. 18, 1999     

基于时间步长自适应技术的海上浮基风电平台绕流数值模拟分析

圆柱绕流问题对于研究海上浮基风电平台在波浪和海流作用下的动力特性以及开发深海风能具有重要的理论和工程应用价值,很多不可压缩流体力学数值模拟方法都基于圆柱绕流的实验或计算结果进行验证。基于自适应时间步长理论及小雷诺数（Re=100）情况下,采用有限体积法,借助FLUENT软件中的用户自定义（UDF）功能,通过二次开发在FLUENT中实现圆柱绕流的数值模拟,并对计算结果与先前的研究结果作了比较,依此来验证时间步长自适应技术和精细边界层网格设计的合理性。计算结果表明了该方法能有效获得准确的流体动力学参数并提高数值模拟计算精度,为该领域的深入研究提供依据。     

考虑晃荡耦合船舶运动的晃荡压力算法（英文）

The effect of coupling between sloshing and ship motions in the evaluation of slosh-induced interior pressures is studied. The coupling between sloshing loads and ship motions is modelled through a hybrid algorithm which combines a potential flow solution based on transient Green function for the external ship hydrodynamics with a viscous flow solution based on a multiphase interface capturing volume of fluid(VOF) technique for the interior sloshing motion. The coupled algorithm accounts for full nonlinear slosh forces while the external forces on the hull are determined through a blended scheme of linear radiationdiffraction with nonlinear Froude-Krylov and restoring forces. Consideration of this level of nonlinearities in ship motions is found to have non-negligible effects on the slosh-coupled responses and slosh-induced loads. A scheme is devised to evaluate the statistical measure of the pressures through long-duration simulation studies in extreme irregular waves. It is found that coupling significantly influences the tank interior pressures, and the differences in the pressures between coupled and uncoupled cases can be as much as 100% or more. To determine the RAO over the frequency range needed for the simulation studies in irregular waves, two alternative schemes are proposed, both of which require far less computational time compared to the conventional method of finding RAO at each frequency, and the merits of these are discussed.     

Simulation of unsteady free-surface flow around a ship hull using a fully coupled multi-phase flow method

An accurate, efficient algorithm for solving free surface flows around ship hulls is described. Accuracy is achieved using a compressive advection discretization which maintains a sharp free surface interface representation without relying on a small time step. Efficiency is obtained using a solution algorithm which implicitly couples velocity, pressure, and volume fractions. The algorithm has been implemented in a computational fluid dynamics code called CFX and is validated by means of a comparisonwith experimental data of benchmark cases in both steady and transient conditions.     

Numerical study on active wave devouring propulsion

The possibility of extracting energy from gravity waves for marine propulsion was numerically studied by a two-dimensional oscillating hydrofoil in this study. The commercially available computational fluid dynamics software FLUENT was used for the unstructured grid based on the Reynolds-average Navier?CStokes equation. The free surface waves and motion of the flapping foil were implemented by customizing the FLUENT solver using a user-defined function technique. In addition, dynamic mesh technology and post processing capabilities were fully utilized. The validation of the model was carried out using experimental data for an oscillation hydrofoil under the waves. The results of the simulation were investigated in detail in order to explain the increase of propeller efficiency in gravity waves. Eight design parameters were identified and it was found that some of them greatly affected the performance of wave energy extraction by the active oscillating hydrofoil. Finally, the overall results suggested that when the design parameters are correctly maintained, the present approach can increase the performance of the oscillating hydrofoil by absorbing energy from sea waves.     

The two-dimensional study of the interaction between liquid loshing and elastic structures

Sloshing phenomenon in the liquid cargo carriers has caught the attention of researchers as the interaction between the sloshing waves and structure is one of the key point and difficulty in the study of sloshing. In this paper, we captured the free surface with a volume of fluid (VOF) method and then calculated the motions and responses of the structure by adopting the Reynolds-averaged Navier-Stokes (RANS) equations for the whole fluid domain. With the use of user defined functions (UDF) in Fluent, the interaction between fluid and structure was then simulated. As a reasonable simplification, the authors studied the response of a single cantilever in a tank under sloshing loads; Further study should pay more attention to the mechanisms of interaction between sloshing waves and elastic structures.     

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.     

Improving longitudinal motion prediction of hybrid monohulls with the viscous effect

A new method improves prediction of the motion of a hybrid monohull in regular waves. Stem section hydrodynamic coefficients of a hybrid monohull with harmonic oscillation were computed using the Reynolds Averaged Navier-Stokes Equations (RANSE). The governing equations were solved using the finite volume method. The VOF method was used for free surface treatment, and RNGК-ε turbulence model was employed in viscous flow calculation. The whole computational domain was divided into many blocks each with structured grids, and the dynamic process was treated with moving grids. Using a 2-D strip method and 2.5D theory with the correction hydrodynamic coefficients allows consideration of the viscous effect when predicting longitudinal motion of a hybrid monohull in regular waves. The method is effective at predicting motion of a hybrid monohull, showing that the viscous effect on a semi-submerged body cannot be ignored.     

舰船气泡尾流数值模拟(英文)

Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture multiphase model to simulate the bubbly wake around the KCS hull. The realizable k-ε two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.     

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

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

有限元法在非线性波浪与结构物相互作用中的应用(英文)

We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.     

海上浮式结构物的水动力仿真研究

基于粘性计算流体力学的方法建立三维数值波浪水池,模拟有限振幅波的传播,并计算三维球体在规则波环境下所受的波浪力。采用两相不可压的RANS方程求解非定常不可压缩粘性流体,并采用流体体积函数(VOF)法对自由面进行动态模拟。通过编写用户自定义函数(UDF)设置边界入口速度和波高,实现在波浪水池尾部1~2倍波长区域消波,最终求出有限振幅波的模拟结果以及规则波中三维球体所受的波浪力,该结果与势流理论边界元法得到的结果在趋势上吻合良好。该研究方法为模拟分析其他海洋结构物在波浪中的水动力奠定基础,丰富与扩展了数值波浪水池的应用。     

水锤冲击时管路系统流固耦合响应的特征线分析方法研究

本文以Wiggert 和Hatfield[3]的特征线分析方法为基础,研究管路在水锤冲击下考虑泊松耦合时流体和结构的瞬态响应.推导了分别对应于管中流体压缩波和管壁中纵波的特征关系式和相容方程,并联合采用空间插值和显式时间插值进行数值求解.针对经典的水锤压力冲击的算例,将边界和初始条件离散化,编制了MATLAB程序进行计算,获得管道中流体压力和流速及轴向应力和振动速度的时程曲线,计算结果与理论分析相当吻合.根据计算结果,提出了对于实际管路设计和水锤防护有益的结论.鉴于特征线法本身的优势,有望在管路系统抗冲击设计分析和防护研究中得到进一步的应用.     

高速三体船波浪中运动与增阻CFD计算研究

基于数值波浪水池技术,对波浪中高速三体船运动及增阻进行CFD计算研究。控制方程—RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波。首先对规则波顶浪中单个主船体的运动和增阻进行了计算,并与高速细长体理论计算结果进行了比较;随后进行了三体船运动和增阻计算,分析了侧片体对主船体阻力增加的影响。为高速三体船耐波性研究提供了数值工具。     

高速滑行艇在规则波中运动的RANSE数值模拟(英文)

This paper presents a study on the numerical simulation of planing crafts sailing in regular waves. This allows an accurate estimate of the seas keeping performance of the high speed craft. The simulation set in six-degree of freedom motions is based on the Reynolds averaged Navier Stokes equations volume of fluid (RANSE VOF) solver. The trimming mesh technique and integral dynamic mesh method are used to guarantee the good accuracy of the hydrodynamic force and high efficiency of the numerical simulation. Incident head waves, oblique waves and beam waves are generated in the simulation with three different velocities (Fn =1.0, 1.5, 2.0). The motions and sea keeping performance of the planing craft with waves coming from different directions are indicated in the flow solver. The ship designer placed an emphasis on the effects of waves on sailing amplitude and pressure distribution of planing craft in the configuration of building high speed crafts.     

粘性流体晃荡与弹性结构的相互耦合作用

提出了处理粘性流体晃荡与弹性结构相互耦合作用的理论及相应数值计算方法.流体运动采用N-S方程描述,控制方程采用有限差分法离散,并由超松弛迭代法求解.液体自由表面通过流体体积法进行重构.为了考虑液舱结构变形对液体晃荡的影响,应用了FAVOR技术.由此建立了描述液体晃荡与结构相互耦合作用的水弹性理论,并由相应的计算方法进行了系列运算,以考察二维液舱不同结构刚度对液体晃荡的影响.计算结果有助于进一步理解此类问题的物理现象.     

滑行艇高速航行时的数值模拟(英文)

Planing vessels are applied widely in civil and military situations.Due to their high speed,the motion of planning vessels is complex.In order to predict the motion of planning vessels,it is important to analyze the hydrodynamic performance of planning vessels at high speeds.The computational fluid dynamic method(CFD) has been proposed to calculate hydrodynamic performance of planning vessels.However,in most traditional CFD approaches,model tests or empirical formulas are needed to obtain the running attitude of the planing vessels before calculation.This paper presents a new CFD method to calculate hydrodynamic forces of planing vessels.The numerical method was based on Reynolds-Averaged Navier-Stokes(RANS) equations.The volume of fluid(VOF) method and the six-degrees-of-freedom equation were applied.An effective process was introduced to solve the numerical divergence problem in numerical simulation.Compared with experimental results,numerical simulation results indicate that both the running attitude and hydrodynamic performance can be predicted well at high speeds.     

SWT方法在内部充水弹性球壳散射的应用

采用Sommerfeld-Watson变换(SWT)方法研究了内部充水弹性球壳的散射,应用绕数积分求根方法在频率域复平面和波数域复平面求解环绕波的频散方程。比较了弹性球壳在内部真空和内部充水两种情况下环绕波相速度曲线和衰减曲线,结果表明内部流体负载使得相对于内部真空情况有大量流体环绕波产生,相邻流体环绕波相速度曲线的"排斥"现象是明显的,流体波向周围媒质的再辐射作用在其相速度接近内部真空球壳弹性环绕波相速度的时候达到局部极值。用镜反射波和环绕波的叠加合成反向散射形态函数,与简正级数解符合良好。     

船行波对直立岸壁的作用研究

以Wigley数学船型为对象,采用传统的一阶Rankine源面元方法,对在浅水情况下的Wigley单体船、双体船以及三体船的艉浪进行数值计算,并用得到的浅水船行波对岸壁的作用力进行计算,结果表明该方法有一定的可行性。     

基于OpenFOAM的螺旋桨敞水性能预报方法

为深入研究螺旋桨周围的粘流问题,基于面向对象的开源CFD计算平台OpenFOAM,选择DTMBP4119桨作为对象,利用RANS方程计算了桨的敞水性能,并分别考察了网格依赖性和离散格式的影响。同时分析了不同半径处叶剖面的压力分布及桨前后的流场速度分布。为保证对流项离散的稳定性和高精度性,采用了二阶NVD格式Gamma混合差分格式,选取k-wSST模型作为湍流模型,压力速度耦合采用SIMPLE松弛算法。计算结果与试验数据吻合较好。研究建立了基于OpenFOAM的螺旋桨敞水性能预报方法。