A RANS-based technique for simulation of the flow near a rolling cylinder at the free surface

Accurate prediction of roll motion is very important in ship dynamics, particularly when the safety of the ship is in question. In this paper "COMET," a Reynolds averaged Navier–Stokes (RANS) solver, is used to simulate the flow in the vicinity of a rectangular cylinder rolling at the free surface and to demonstrate the potential of RANS-based techniques in complex flows with free surfaces. The method adopted is based on an unstructured collocated finite-volume technique, which uses a variation of the SIMPLE scheme for pressure correction. An advanced interface capturing technique known as high resolution interface capturing (HRIC) is used for tracking the free surface with reduced numerical diffusion. The results derived are compared with available numerical and experimental results, showing good agreement on added moment of inertia and damping coefficients, whilst accurately predicting other main features of the flow. However, further validation for other geometries and extension to 3-D hull forms are essential before applications to practical problems are realised. Received: January 26, 2000 / Accepted: October 2, 2000     

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.     

Self-propulsion computations using a speed controller and a discretized propeller with dynamic overset grids

A method that can be used to perform self-propulsion computations of surface ships is presented. The propeller is gridded as an overset object with a rotational velocity that is imposed by a speed controller, which finds the self-propulsion point when the ship reaches the target Froude number in a single transient computation. Dynamic overset grids are used to allow different dynamic groups to move independently, including the hull and appendages, the propeller, and the background (where the far-field boundary conditions are imposed). Predicted integral quantities include propeller rotational speed, propeller forces, and ship’s attitude, along with the complete flow field. The fluid flow is solved by employing a single-phase level set approach to model the free surface, along with a blended k−ω/k−ɛ based DES model for turbulence. Three ship hulls are evaluated: the single-propeller KVLCC1 tanker appended with a rudder, the twin propeller fully appended surface combatant model DTMB 5613, and the KCS container ship without a rudder, and the results are compared with experimental data obtained at the model scale. In the case of KCS, a more complete comparison with propulsion data is performed. It is shown that direct computation of self-propelled ships is feasible, and though very resource intensive, it provides a tool for obtaining vast flow detail.     

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.     

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

本文采用有限体积法通过求解不可压缩的雷诺平均(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.     

Level set方法在空化流动计算中的应用

文章将Level set方法应用于自然空化与通气空化流动计算中,通过耦合入相间的质量传输方程,建立了一个用于描述气、汽、液多相空化流动的数值模拟方法,分别对绕圆头回转体和Clark-Y型水翼的自然空化流动和绕圆盘的通气空化流动进行了数值模拟研究,并与实验结果进行了对比,研究结果表明:相较常用的多相流模型(mixture model),Lev-el set方法通过时间和空间上的压缩离散方案,减小了相间界面的扩散,准确地捕捉到了相间的界面,可以有效地运用于空化多相流动。     

一种计算水翼水动力的三维面元法

本文用Rankine奇点面元法计算了深、浅水中三维水翼的定常升力绕流。水翼的厚度和升力效应分别以水翼表面分布的Rankine源和法向偶极子来模拟，在自由表面上也分布Rankine源，通过满足相应的边界条件和尾缘处的Kutta条件求出这些奇点强度。以在自由表面下作小攻角定常运动的水翼为例进行了计算，计算结果与试验结果和其他计算结果作了比较。     

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

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

Fully three-dimensional ship seakeeping computations with a surge-corrected Rankine panel method

A 3-d seakeeping code uses first-order Rankine panels with special numerical integration on the ship's hull and Rankine point source clusters above the free surface. The code computes the motions of the ship in regular waves of small height (linearized). The steady flow is captured without simplification by solving the fully nonlinear wave-resistance problem first. A special treatment of the surge motion considers the influence of periodic quantities on thrust and resistance, and improves surge motion predictions. Radiation and open-boundary conditions are enforced by staggered grids. Results for the ITTC standard test case S-175 containership agree well with experiments except for very long waves. The importance of capturing the three-dimensional steady flow contributions is also demonstrated.     

URANS and DES analysis for a Wigley hull at extreme drift angles

Vortical structures and associated instabilities for flows around the Wigley hull for a wide range of drift angles (10° ≤ α ≤ 60°) with free surface are identified and analyzed. Quantitative verification and validation are conducted on three systematically refined grids with comparison to the experimental data for α = 10°. Analysis of the flow pattern shows a strong correlation between the vortical structures and free-surface wave elevation. For α = 10° and 30°, the flows remain steady and vortices are generated at the keel and fore and aft perpendiculars of the hull. The strength and complexity of these vortices increase with increasing α. At α = 45°, flow becomes unsteady without any significant change in the main flow pattern. At α = 60°, a complex and unsteady flow field on the leeward side of the hull is formed with a large recirculation region from the aft to the fore end, which prevents the flow coming from below the keel from moving up and generating the keel vortices observed at lower drift angles. Karman-like and helical instabilities are analyzed. The effect of Froude number is more apparent for large than for small drift angles.     

A local model for the simulation of two–dimensional spilling breaking waves

A numerical model for the simulation of two dimensional spilling breaking waves is described. It is derived from a previous model which, in turn, takes its underlying ideas from the Cointe and Tulin theory of steady breakers. With respect to the former model, the present one is local, i.e., the inception, extension, and geometry of the breakers are determined through the local shape of the water surface. The model has been implemented in a RANSE code, which was developed for the simulation of ship flows, through a modification in the boundary conditions. This yields a simple and effective way to reproduce the breakers influence on the underlying flow. The resulting code has been used to simulate the flow past a submerged hydrofoil. The numerical results are compared with those of the previous model and with the experimental data obtained by Duncan.     

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.

     

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

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

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

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.     

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

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

Assessment of the volume of fluid method for free-surface wave flow

This study was concerned with the free-surface wave flow around a surface-piercing foil. The volume of fluid method implemented in a Navier–Stokes computational fluid dynamics code was employed. Three widely used discretization schemes for the volume of fluid method were assessed for a test case that involved general ship waves, spilling breaking waves in front of the leading edge, and bubbly free surfaces in separated regions. A single computational approach was selected for the comparison, and a grid-dependence study was carried out. The computational results were validated against existing experimental data, showing good agreement. The validation results suggest that all three discretization schemes perform well, but the best and most efficient results were obtained using the high-resolution interface capturing scheme.     

破损船舶进水模拟

本文应用RANS方程进行了破损舰艇的进水模拟。使用了流体体积法、两方程湍流模型及壁函数处理自由面附近湍流流动，给出了后台阶流动、 系列60以及破损Burgundy船的CFD模拟计算结果，捕捉了船舶破损相关的基本流动特征。     

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.     

用naoe-FOAM-SJTU求解器数值模拟船舶在波浪上的操纵运动问题（英文）

Ship maneuvering in waves includes the performance of ship resistance, seakeeping, propulsion, and maneuverability. It is a complex hydrodynamic problem with the interaction of many factors. With the purpose of directly predicting the behavior of ship maneuvering in waves, a CFD solver named naoe-FOAM-SJTU is developed by the Computational Marine Hydrodynamics Lab(CMHL) in Shanghai Jiao Tong University. The solver is based on open source platform OpenFOAM and has introduced dynamic overset grid technology to handle complex ship hull-propeller-rudder motion system. Maneuvering control module based on feedback control mechanism is also developed to accurately simulate corresponding motion behavior of free running ship maneuver. Inlet boundary wavemaker and relaxation zone technique is used to generate desired waves. Based on the developed modules, unsteady Reynolds-averaged Navier-Stokes(RANS) computations are carried out for several validation cases of free running ship maneuver in waves including zigzag, turning circle, and course keeping maneuvers. The simulation results are compared with available benchmark data. Ship motions, trajectories, and other maneuvering parameters are consistent with available experimental data, which indicate that the present solver can be suitable and reliable in predicting the performance of ship maneuvering in waves. Flow visualizations, such as free surface elevation, wake flow, vortical structures, are presented to explain the hydrodynamic performance of ship maneuvering in waves. Large flow separation can be observed around propellers and rudders. It is concluded that RANS approach is not accurate enough for predicting ship maneuvering in waves with large flow separations and detached eddy simulation(DES) or large eddy simulation(LES) computations are required to improve the prediction accuracy.