Resistance and propulsion characteristics of various commercial ships based on CFD results

This paper uses computational tools to examine the speed performance of various types of commercial ships including resistance and propulsion characteristics. Eight commercial ships built in the last decade were selected for the study. They include four large-sized container carriers, one bulk carrier, one VLCC, and two LNG carriers. The Reynolds averaged Navier-Stokes equation has been utilized, and the computations were executed under the same conditions of the model tests to predict the speed performance, i.e., resistance and self-propulsion. The self-propulsion point was obtained from load-varying tests. The speed performance was predicted based on the model-ship performance analysis method of the revised ITTC’78 method. The limiting streamlines on the hull, wave characteristics around the model ship, and the wake characteristics on the propeller plane were also investigated. After completing the computations, a series of model tests were conducted to evaluate the accuracy of the computational predictions. The predictions clearly reveal the differences in the resistance and propulsion characteristics regarding the various types of commercial ships, and may be applicable to hull-form design.     

碎冰条件下单点系泊船舶动态响应特性研究

为了研究不同单点系泊位置对极地系泊船舶在来冰角度改变时动态响应的影响,运用离散元方法(DEM),计算了船舶在60%密集度下,冰速为1.0 m/s、冰厚为1 m时不同角度的船冰阻力,结合系泊分析软件二次开发功能,计算分析系泊船舶在不同来冰角度和不同系泊位置下受上述冰情影响下的响应。通过对模拟结果与试验结果的对比,认为离散元方法(DEM)的数值模拟结果与试验结果吻合较好;单点系泊位置为距船头L/4时,系泊船舶能以最快速、最稳定的方式对来冰方向的变化做出响应。     

Influence of Ice Size Parameter Variation on Hydrodynamic Performance of Podded Propulsor

During ice-breaking navigation, a massive amount of crushed ice blocks with different sizes is accumulated under the hull of an ice-going ship. This ice slides into the flow field in the forward side of the podded propulsor, affecting the surrounding flow field and aggravating the non-uniformity of the propeller wake. A pulsating load is formed on the propeller, which affects the hydrodynamic performance of the podded propulsor. To study the changes in the propeller hydrodynamic performance during the ice podded propulsor interaction, the overlapping grid technique is used to simulate the unsteady hydrodynamic performance of the podded propulsor at different propeller rotation angles and different ice block sizes. Hence, the hydrodynamic blade behavior during propeller rotation under the interaction between the ice and podded propulsor is discussed. The unsteady propeller loads and surrounding flow fields obtained for ice blocks with different sizes interacting with the podded propulsor are analyzed in detail. The variation in the hydrodynamic performance during the circular motion of a propeller and the influence of ice size variation on the propeller thrust and torque are determined. The calculation results have certain reference significance for experiment-based research, theoretical calculations and numerical simulation concerning ice podded propulsor interaction.     

A comparative study of numerical simulations for fluid–structure interaction of liquid-filled tank during ship collision

Although International Maritime Organization (IMO) has taken many measures to minimize ship collisions, ships carrying liquid cargo sometimes do get struck by other vessels. The outflow of crude oil causes very serious consequences to the environment. In such cases it is necessary to analyze the response of structure of struck liquid cargo-filled tank to account for fluid–structure interaction accurately. In this paper, numerical simulation of collision between a container ship with double hull very large crude carrier (VLCC) is presented. Three different numerical simulation mothods were adopted to model fluid–structure interaction in liquid-filled cargo tank, namely arbitrary Lagrangian–Eulerian finite element method, Lagrangian finite element method and linear sloshing model. The numerical simulation results reveal that the fluid–structure interaction of liquid cargo-filled tank has a significant effect on the motion and structural response of the struck cargo tank. Compared with the calculation results of ALE FE method, the linear sloshing model underestimates the influence of fluid–structure interaction of liquid cargo tank while the Lagrangian–Eulerian finite element method may be considered as the practical method for engineering applications as it provided more reasonable results with a relatively low central processing unit (CPU) time.     

Calculation Methods of Added Resistance and Ship Motion Response Based on Potential Flow and Viscous Flow Theory

To improve the energy efficiency of ships and to predict ship motion response under actual sea conditions, the far-field theory, strip theory, and Fujii and Takahashi’s modified semi-empirical method are based and studied to calculate the wave-induced added resistance. Firstly, a new modified formula based on the Maruo method is presented to calculate the radiation added resistance for the ship with a complex surface. Meanwhile, some calculation details such as the Green function, the shape of the sections (shape below the still water level or shape below the wave level) in the strip theory, and so on are discussed. Finally, the CFD method is used to simulate the motions of the hull and the added resistance, and the results of the CFD method and those of other numerical methods are analyzed and compared with the experiment results. The modified method in the paper can predict the added resistance in waves for the complex-hull-surface ships well and quickly.

     

Some methods to obtain the added resistance of a ship advancing in waves

Added resistance in waves is an important part of ship dynamic due to its economical effect on ship exploitation. Can the ship sustain speed in a rough sea state? If not, this will produce delays and economical losses if added resistance is not taken into account in the propulsion design. There are not many simple methods to obtain the added resistance in waves of a ship, and the validity of the results is not always good enough for different ships. In this paper, some theories that can be used to predict the added resistance of a ship are studied and validated against seakeeping tests of some monohull models. Tests and results focus in head seas, which are the most severe for the added resistance. Experimental results are compared with numerical calculations and conclusions about the range of application of the presented theories are obtained.     

冰脊形成过程的离散元模拟及影响因素分析

冰脊对极地船舶及海洋工程结构的冰载荷设计及冰区安全运行具有重要的影响,其几何形态是重要的考虑因素。为研究冰脊的几何形态,本文基于Voronoi切割算法构造扩展多面体海冰单元,并使其相互冻结形成平整冰;采用离散单元法模拟海冰在相对运动时挤压形成冰脊的动力过程,统计分析脊帆高度、龙骨宽度、龙骨深度、龙骨水平倾角等主要几何参数,并确定这些几何参数间的对应关系。将离散元模拟结果同冰脊现场观测资料对比分析以验证扩展多面体离散元方法模拟冰脊形成过程的可行性。在此基础上讨论冰厚、冰速、海冰强度等因素对冰脊形态的影响。本文采用离散元方法对冰脊形成过程及几何特性的数值分析可为深入开展冰脊的形成机理及其对船舶、海洋工程结构物的冰载荷分析提供参考依据。     

GPU-Based DEM Simulations of Global Ice Resistance on Ship Hull During Navigation in Level Ice

The ice resistance on a ship hull affects the safety of the hull structure and the ship maneuvering performance in icecovered regions.In this paper,the discrete element method(DEM)is adopted to simulate the interaction between level ice and ship hull.The level ice is modeled with 3D bonded spherical elements considering the buoyancy and drag force of the water.The parallel bonding approach and the de-bonding criterion are adopted to model the freezing and breakage of level ice.The ship hull is constructed with rigid triangle elements.To improve computational efficiency,the GPU-based parallel computational algorithm was developed for the DEM simulations.During the interaction between the ship hull and level ice,the ice cover is broken into small blocks when the interparticle stress approaches the bonding strength.The global ice resistance on the hull is calculated through the contacts between ice elements and hull elements during the navigation process.The influences of the ice thickness and navigation speed on the dynamic ice force are analyzed considering the breakage mechanism of ice cover.The Lindqvist and Riska formulas for the determination of ice resistance on ship hull are employed to validate the DEM simulation.The comparison of results of DEM,Lindqvist,and Riska formula show that the DEM result is between those the Lindqvist formula and Riska formula.Therefore the proposed DEM is an effective approach to determine the ice resistance on the ship hull.This work can be aided in the hull structure design and the navigation operation in ice-covered fields.     

Unsteady ship waves in shallow water of varying depth based on Green–Naghdi equation

Based on Green–Naghdi equation this work studies unsteady ship waves in shallow water of varying depth. A moving ship is regarded as a moving pressure disturbance on free surface. The moving pressure is incorporated into the Green–Naghdi equation to formulate forcing of ship waves in shallow water. The frequency dispersion term of the Green–Naghdi equation accounts for the effects of finite water depth on ship waves. A wave equation model and the finite element method (WE/FEM) are adopted to solve the Green–Naghdi equation. The numerical examples of a Series 60 (C_B=0.6) ship moving in shallow water are presented. Three-dimensional ship wave profiles and wave resistance are given when the ship moves in shallow water with a bed bump (or a trench). The numerical results indicate that the wave resistance increases first, then decreases, and finally returns to normal value as the ship passes a bed bump. A comparison between the numerical results predicted by the Green–Naghdi equation and the shallow water equations is made. It is found that the wave resistance predicted by the Green–Naghdi equation is larger than that predicted by the shallow water equations in subcritical flow , and the Green–Naghdi equation and the shallow water equations predict almost the same wave resistance when , the frequency dispersion can be neglected in supercritical flows.     

Optimized support vector regression algorithm-based modeling of ship dynamics

This study contributes to solving the problem of how to derive a simplistic model feasible for describing dynamics of different types of ships for maneuvering simulation employed to study maritime traffic and furthermore to provide ship models for simulation-based engineering test-beds. The problem is first addressed with the modification and simplification of a complicated and nonlinearly coupling vectorial representation in 6 degrees of freedom (DOF) to a 3 DOF model in a simple form for simultaneously capturing surge motions and steering motions based on several pieces of reasonable assumptions. The created simple dynamic model is aiming to be useful for different types of ships only with minor modifications on the experiment setup. Another issue concerning the proposed problem is the estimation of parameters in the model through a suitable technique, which is investigated by using the system identification in combination with full-scale ship trail tests, e.g., standard zigzag maneuvers. To improve the global optimization ability of support vector regression algorithm (SVR) based identification method, the artificial bee colony algorithm (ABC) presenting superior optimization performance with the advantage of few control parameters is used to optimize and assign the particular settings for structural parameters of SVR. Afterward, the simulation study on identifying a simplified dynamic model for a large container ship verifies the effectiveness of the optimized identification method at the same time inspires special considerations on further simplification of the initially simplified dynamic model. Finally, the further simplified dynamic model is validated through not only the simulation study on a container ship but also the experimental study on an unmanned surface vessel so-called I-Nav-II vessel. Either simulation study results or experimental study results demonstrate a valid model in a simple form for describing the dynamics of different types’ ships and also validate the performance of the proposed parameter estimation method.     

深远海养殖工船最小推进功率研究

针对十万吨级深远海养殖工船在波浪中的阻力和运动响应开展了船模试验与模拟计算,并对其最小推进功率进行了校核。以试验流体力学（EFD）模型试验与模型尺度计算流体力学（CFD）模拟计算的结果进行对比作为方法验证,将实尺度CFD模拟计算的结果直接用于最小推进功率的校核。结果显示,EFD模型试验结果与模型尺度CFD计算结果误差在10%左右,计算精度满足工程要求。在规定的恶劣海况下该船型波浪增阻占总阻力的比例最高可达56.3%,螺旋桨转矩可达最大转矩的63.9%。等级2简化评估法得到的最小推进功率为等级1线评估法给出的最小功率线值的58.9%。研究表明实尺度CFD模拟计算可直接用于船舶最小推进功率的校核,此深远海养殖工船使用等级2简化评估方法进行校核更容易满足规范要求,其最小推进功率应不小于6 833 kW。     

预测细粒酒精模型冰物理和力学参数的实验关系

依据模型冰物理性质控制其力学性质的原则 ,利用能量平衡观点建立了冰内未冻结液体含量同制冰过程中的冰池温度、降温及升温时间、模型冰厚度 (喷雾时间×流量 )的理论关系 ,以其量纲表达式作为预测模型冰物理和力学参数的综合控制指标 ;通过系统测量酒精溶液冰点温度、冰池大厅气温、模型冰温度、冰内未冻液含量、冰密度、冰压缩强度、冰弯曲强度和冰变形模量 ,获得了不同冰点温度模型冰的物理和力学参数同该综合控制指标的实验关系 ;实现了用人为可控制因子定量控制模型冰物理和力学参数的目的     

极地吊舱推进船舶回转性能的离散元分析

吊舱推进装置在极地船舶中被广泛应用,其转舵模块可以带动螺旋桨摆动从而产生任意方向的推进力,使船舶操纵更为灵活。提出了一种吊舱推进船舶冰区操纵的离散元方法,对具有吊舱推进装置的冰区船舶破冰过程进行了数值模拟。以“雪龙2”号破冰船为研究对象,计算分析了船舶定速直航时的冰阻力,并通过与Lindqvist经验公式的对比验证了冰阻力计算的合理性。在此基础上进一步对船舶在不同冰厚、不同吊舱转向角下的回转破冰运动进行了离散元模拟,分析了回转半径与船长比值随冰厚的变化规律。计算结果表明：船舶的回转性能随冰厚的增大而降低,并随吊舱转向角的增大而显著提高。     

On the parametric rolling of ships using a numerical simulation method

This paper has shown a numerical motion simulation method which can be employed to study on parametric rolling of ships in a seaway. The method takes account of the main nonlinear terms in the rolling equation which stabilize parametric rolling, including the nonlinear shape of the righting arm curve, nonlinear damping and cross coupling among all 6 degrees of freedom. For the heave, pitch, sway and yaw motions, the method uses response amplitude operators determined by means of the strip method, whereas the roll and surge motions of the ship are simulated, using nonlinear motion equations coupled with the other 4 degrees of freedom. For computing righting arms in seaways, Grim's effective wave concept is used. Using these transfer functions of effective wave together with the heave and pitch transfer functions, the mean ship immersion, its trim and the effective regular wave height are computed for every time step during the simulation. The righting arm is interpolated from tables, computed before starting the simulation, depending on these three quantities and the heel angle. The nonlinear damping moment and the effect of bilge keels are also taken into account. The numerical simulation tool has shown to be able to model the basic mechanism of parametric rolling motions. Some main characteristics of parametric rolling of ships in a seaway can be good reproduced by means of the method. Comprehensive parametric analyses on parametric rolling amplitude in regular waves have been carried out, with that the complicated parametric rolling phenomena can be understood better.     

螺旋桨对海冰切削作用的DEM-FEM耦合分析

当船舶在冰区航行时,螺旋桨会与海冰相互碰撞并导致桨叶的变形和损坏,进而影响船舶的航行安全。为研究海冰与螺旋桨的相互作用过程,采用离散元(DEM)—有限元(FEM)耦合方法构建海冰—螺旋桨切削模型。海冰和螺旋桨模型分别采用具有黏结—破碎特性的球体离散单元和8节点六面体有限单元构造。基于该DEM-FEM耦合模型讨论了不同切削深度下,螺旋桨所承受冰载荷的特点和规律;最后,研究了螺旋桨切削海冰过程中进速系数、推力系数、扭矩系数之间的对应关系,并讨论了海冰—螺旋桨相互作用过程中冰压力、Mises应力和变形的分布特点。以上研究可为寒区船舶安全航行和螺旋桨设计提供有益的参考。     

艏艉线型优化对船舶阻力性能的影响

为提高母型船阻力性能,以船体阻力性能为优化对象,基于改造母型船法,研究船舶球鼻艏以及船尾线型的改变对船舶阻力性能的影响.采用高度集成化的Tribon系统、可视化绘图软件Auto CAD及CFD(Computational Fluid Dynamics)通用前处理软件ICEM联合建模的方法来建立船体模型.通过模拟计算结果与实验值的对比分析,验证CFD技术在船舶阻力性能预报中的合理性和有效性.通过对比3种不同球鼻艏时的船体阻力得知:从阻力性能方面考虑,对于低速丰满型船舶选用普通型球鼻艏以及中高速船舶采用上翘型球鼻艏均可以获取较好地减阻效果.同时比较不同航速下尾部线型对船体总阻力的影响表明,选优后的方形尾在相同的航速下阻力低、消耗的功率小、形状效应小、黏压阻力和摩擦阻力也相对较小.     

Ice resistance measurements in ridges with IB APU in the Baltic sea

Ice ridges form a difficult obstacle for ice navigation. Despite this fact, the resistance of ships in ridges has been investigated very little.A trial was performed with the Finnish icebreaker APU (propulsion machinery power of 8·8 MW) in April 1974 in the Baltic, in order to measure its resistance in ridges.Most of the 11 ridges were so massive that the ship could not penetrate through them by a continuous mode, the ramming mode had to be used.When the ship was penetrating a ridge, its speed and the propeller revolutions were registered for determining the resistance. A doppler radar was used for measuring ship spped.In case the ship was stopped in the ridge it was extracted, reversed and accelerated for the next ram. The time history was registered for determining the ship's speed of advance.Before starting a test, the ridge profile above the water level was measured. This was simply done by using a levelling instrument.The main object of the test was to determine the ship's speed of advance. The results, i.e. ice resistance, ship speed and ridge characteristics were analysed on three levels: momentary values average values for rams and average values for a series of rams in one ridge.The test series presented in the article was the first one in full scale in which the speed of the advance of a ship moving by ramming was determined and the ridge profiles were mapped. The measuring system developed worked well. As results the tests gave data of the ship's resistance and of the ship's speed of advance in ridges.     

锥体导管架海洋平台冰激振动的DEM-FEM耦合分析及高性能算法

在寒区海域,冰荷载是影响海洋平台安全运行的主要环境荷载之一,由其引起的冰激振动给平台结构及其上部设备带来了严重危害。为分析不同冰况下平台的振动响应,本文建立了导管架海洋平台冰激振动的离散元（DEM）-有限元（FEM）耦合模型。采用具有黏接-破碎性能的等粒径球体离散单元对海冰的破碎特性进行模拟,通过由梁单元构建的海洋平台有限元模型获得结构的振动响应。在离散元与有限元的接触区域中实现了两个模型间计算参数的传递。为提高该耦合模型的计算效率和规模,发展了基于动力子结构方法的DEM-FEM耦合模型。为验证该耦合模型的有效性和可行性,将不同冰况下得到的冰荷载与ISO19906和JTS 144-1-2010标准进行了对比。结果表明,计算得到的冰荷载与标准相近,且冰厚与冰荷载呈二次非线性关系。同时,从冰速和冰厚两方面对比了渤海四桩腿JZ20-2 MUQ锥体导管架平台冰激振动加速度的数值结果和现场实测数据,发现冰速与振动加速度呈线性关系,冰厚与振动加速度呈二次非线性关系,并且振动加速度与冰速和冰厚平方的乘积呈线性关系。     

Numerical simulation of PMM tests for a ship in close proximity to sidewall and maneuvering stability analysis

As the maneuverability of a ship navigating close to a bank is influenced by the sidewall, the assessment of ship maneuvering stability is important. The hydrodynamic derivatives measured by the planar motion mechanism (PMM) test provide a way to predict the change of ship maneuverability. This paper presents a numerical simulation of PMM model tests with variant distances to a vertical bank by using unsteady RANS equations. A hybrid dynamic mesh technique is developed to realize the mesh configuration and remeshing of dynamic PMM tests when the ship is close to the bank. The proposed method is validated by comparing numerical results with results of PMM tests in a circulating water channel. The first-order hydrodynamic derivatives of the ship are analyzed from the time history of lateral force and yaw moment according to the multiple-run simulating procedure and the variations of hydrodynamic derivatives with the ship-sidewall distance are given. The straight line stability and directional stability are also discussed and stable or unstable zone of proportional-derivative (PD) controller parameters for directional stability is shown, which can be a reference for course keeping operation when sailing near a bank.     

Full scale CFD seakeeping simulations for case study ship redesigned from V-shaped bulbous bow to X-bow hull form

Increasing propulsion efficiency, safety, comfort and operability are of the great importance, especially for small ships operating on windy sites like the North Sea and the Baltic Sea. Seakeeping performance of ships and offshore structures can be analysed by different methods and the one that is becoming increasingly important is CFD RANS. The recent development of simulation techniques together with rising HPC accessibility enables performance of advanced seakeeping simulations for ships in a full scale. The paper presents CFD seakeeping analysis for a case study vessel in two variants: V-shaped bulbous bow hull form (as built) and innovative hull form (X-bow type). The study presents the influence of redesigning the ship on selected seakeeping aspects. The advanced CFD model, with the application of overset mesh technique, was described in detail. Selected numerical results were validated on the basis of experimental testing in a towing tank and showed good agreement. The approach demonstrated here of performing the CFD seakeeping simulations for the analysis of ship performance in a full scale and close to real loading conditions has direct application to the design process as well as in determination of optimal operational parameters of any ship.