极地船舶冰区操纵性能的六自由度离散元分析

极地船舶操纵破冰性能是破冰船设计建造过程中需重点考虑的问题。为分析极地船舶操纵性能,本文发展了冰区船舶六自由度操纵破冰运动模型,采用扩展多面体离散元模拟海冰,舵桨操纵模型提供破冰力。开展平整冰区定常直航模拟计算,并与Lindqvist船体冰阻力经验公式展开对比验证;开展雪龙号敞水35°舵角操纵性仿真,并与实船试航结果进行对比。在此基础上,开展冰厚及舵角影响下的船体结构冰载荷及破冰轨迹计算,模拟得到操纵破冰航行中船体垂荡、横摇及纵摇运动时程;最后分析操纵破冰船体线载荷分布与直航破冰下的差异。     

CFD‐based multi‐objective optimization method for ship design

This paper concerns development and demonstration of a computational fluid dynamics (CFD)‐based multi‐objective optimization method for ship design. Three main components of the method, i.e. computer‐aided design (CAD), CFD, and optimizer modules are functionally independent and replaceable. The CAD used in the present study is NAPA system, which is one of the leading CAD systems in ship design. The CFD method is FLOWPACK version 2004d, a Reynolds‐averaged Navier–Stokes (RaNS) solver developed by the present authors. The CFD method is implemented into a self‐propulsion simulator, where the RaNS solver is coupled with a propeller‐performance program. In addition, a maneuvering simulation model is developed and applied to predict ship maneuverability performance. Two nonlinear optimization algorithms are used in the present study, i.e. the successive quadratic programming and the multi‐objective genetic algorithm, while the former is mainly used to verify the results from the latter. For demonstration of the present method, a multi‐objective optimization problem is formulated where ship propulsion and maneuverability performances are considered. That is, the aim is to simultaneously minimize opposite hydrodynamic performances in design tradeoff. In the following, an overview of the present method is given, and results are presented and discussed for tanker stern optimization problem including detailed verification work on the present numerical schemes. Copyright © 2006 John Wiley & Sons, Ltd.     

Unsteady RANS method for surface ship boundary layer and wake and wave field

Results are reported of an unsteady Reynolds‐averaged Navier–Stokes (RANS) method for simulation of the boundary layer and wake and wave field for a surface ship advancing in regular head waves, but restrained from body motions. Second‐order finite differences are used for both spatial and temporal discretization and a Poisson equation projection method is used for velocity–pressure coupling. The exact kinematic free‐surface boundary condition is solved for the free‐surface elevation using a body‐fitted/free‐surface conforming grid updated in each time step. The simulations are for the model problem of a Wigley hull advancing in calm water and in regular head waves. Verification and validation procedures are followed, which include careful consideration of both simulation and experimental uncertainties. The steady flow results are comparable to other steady RANS methods in predicting resistance, boundary layer and wake, and free‐surface effects. The unsteady flow results cover a wide range of Froude number, wavelength, and amplitude for which first harmonic amplitude and phase force and moment experimental data are available for validation along with frequency domain, linear potential flow results for comparisons. The present results, which include the effects of turbulent flow and non‐linear interactions, are in good agreement with the data and overall show better capability than the potential flow results. The physics of the unsteady boundary layer and wake and wave field response are explained with regard to frequency of encounter and seakeeping theory. The results of the present study suggest applicability for additional complexities such as practical ship geometry, ship motion, and maneuvering in arbitrary ambient waves. Copyright © 2001 John Wiley & Sons, Ltd.     

Transient response of a ship to an abrupt flooding accounting for the momentum flux

Numerical non-linear time domain simulation method for damaged ship motions is presented. Floodwater motion modelling is based on the lumped mass method with a moving free surface. The ship and floodwater motions are fully coupled. The variation of the floodwater mass is accounted for. A model to account for the flooding ingress transporting the momentum is presented. The experiments of abrupt flooding have shown that the ship may experience the first large roll towards the undamaged side, especially when a large undivided compartment is flooded. The presented time domain model is validated against the experimental data on the roll damping of the flooded ship and transient flooding. Two different initial stability conditions and two different compartment layouts are studied. Viscous dissipation of the floodwater motions is modelled with an equivalent friction coefficient. The impact of the viscous damping is studied. Transient flooding tests show that the inflow momentum has to be accounted for when the undivided compartment is flooded. The simulation model is capable of capturing the impact of the inflooding jet and the first roll on the opposite side of the damage is reproduced.     

船舶在进出船厢运动中的附加质量的计算

采用三维有限元方法，计算船在进出船厢运动中的附加质量。文中对单圆柱体在无限域和半无限域中运动的附加质量及双圆柱体同心放置时内圆柱运动的附加质量分别进行了计算，计算值与理论值吻合很好；对船在进出船厢运动中的附加质量做了大量的计算，结果表明：随着船与船厢侧壁间距的减小，船的附加质量系数增大；随着船厢中水深的减小，船的附加质量系数增大，随着船离船厢封闭端的距离的减小，附加质量系数增大。     

ANALYSIS OF SHIP MANEUVERING PERFORMANCES AND ICE LOADS ON SHIP HULL WITH DISCRETE ELEMENT MODEL IN BROKEN-ICE FIELDS

采用离散元模型对碎冰区浮冰与船舶结构的相互作用进行了数值研究。碎冰由三维圆盘单元构成,并考虑其在海流作用下的浮力、拖曳力和附加质量。船体结构由一系列三角形单元组合构造。通过海冰与船体单元间的接触判断和接触力计算,确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰况（冰速、冰厚、冰块尺寸和密集度）以及航速条件下,海冰对船体的动力作用过程进行了数值分析,对比分析了以上因素对船体冰载荷的影响,可为冰区船舶的安全运行和结构设计提供一定的借鉴作用。     

Simulation of three-dimensional unsteady viscous free surface flow around a ship model

We present here a numerical method for solving the free surface flow around a ship at forward speed in calm water. The fluid is assumed to be Newtonian and the Reynolds-averaged Navier-Stokes equations are solved by a finite difference method. Modelization of turbulence is achieved by the algebraic model proposed by Baldwin and Lomax. Fully non-linear free surface conditions are satisfied in the model and a method to avoid the incompatibility between free surface conditions and no-slip conditions at the waterline is proposed. Numerical results obtained for a Wigley hull are compared with experimental results.     

A method to compute ship exhaust plumes with waves and wind

Thermal and concentration transport models are implemented in CFDShip‐Iowa version 4.5, a semi‐coupled solver for air/water free surface flow (Int. J. Numer. Meth. Fluids 2008; 58 (6):591–624), to investigate the exhaust plume around ship superstructures. An incompressible, variable density approximation is implemented where the density can change in all governing equations due to temperature variations only. The thermal and concentration models are tested for the cases of steady and unsteady flow with thermal and solution transport in a 2D square cavity, and for a 3D thermal plume in an open environment, showing good agreement between computational results and experimental data. To test the method in an extreme motions condition, the exhaust plume of the ONR Tumblehome model DTMB 5613 is studied, showing complicated vortical structures in air including a pair of counter‐rotating vortices downstream of the stack for cross‐flow, and bended bird‐plume shape in the symmetry plane and varying arc‐shape in axial sections both for temperature and NO_x concentration fields. Effects of smoke exhaust speed and wind speed on the temperature and concentration distributions are studied. Finally, a smoke downwash computation is performed for a ship free to move in 6 degrees of freedom in a sea state 8 condition. Copyright © 2010 John Wiley & Sons, Ltd.     

船舶碰撞数值仿真的附加质量模型

根据船舶碰撞的运动滞后和局部损伤特性 ,采用附加质量处理流体 结构耦合作用 ,用详细的有限元模型表达撞击船首和被撞船侧的直接涉撞区结构 ,而将非碰撞区的船体结构视为刚体 ,形成一个附加质量模型用于碰撞仿真计算。附加质量模型不仅建模工作量小 ,而且CPU时间也大大减少。算例表明 ,附加质量模型与流固耦合模型的分析结果具有良好的一致性 ,可以满足一般的工程精度要求。     

Development of a measuring system for segmented ship models

The use of segmented ship models to test and study various ship responses in model scale poses a challenge to instrumentation and test engineers. In recent years, the authors have developed a segmented ship model to study and trace ice loads acting on the ship hull. The model contains three segments at the bow. Each segment is supported at multiple points which enable the resolution of the location, magnitude, and direction of the external loads. The typical segment support system consists of four vertical, two transverse, and one longitudinal support points with uniaxial compression pression transducers. Stabilization of the segment is achieved by using three specially designed tension links acting in the vertical, transverse and longitudinal directions, respectively. The system is subjected to several levels of calibration which include individual transducer calibration, calibration using internal loading applied by the tension links, external calibration using a specially designed and built calibration rig capable of exerting normal and inclined loads, and calibration of the effects of buoyancy changes in a floating model. The results of calibration and ice-breaking tests indicate that position prediction of the external load can be made within 20 mm. The normal load of less than 100 Newtons can be determined within a few percentage points but the frictional load magnitude and direction are found to be subject to greater errors particularly for low friction factors of the order of 0.1. A.M. Nawwar (SEM Member), presently Associate Professor of Mechanical Engineering, Kuwait University, Kuwait 13060     

丰、枯水年磨刀门水道盐水上溯运动规律对比

通过一维纵向逐时咸界图的方法, 得到了磨刀门水道丰、枯水年的盐水运动规律图.对比分析了不同年份不同上游径流对磨刀门水道盐水上溯运动规律的影响. 分析表明, 丰、枯水年不同的上游径流量将改变磨刀门水道盐水入侵的上溯距离和上溯速度的快慢, 造成程度不同的咸潮灾害. 磨刀门水道中盐水运动的半月及日周期性变化规律不受上游径流的影响,均为小潮期间盐水逐步上溯, 大潮期间盐水逐步退出. 上游径流量的大小只影响磨刀门水道盐水上溯运动的速度和距离远近, 不改变磨刀门水道中盐水上溯运动的半月周期规律.      

The application of three-dimensional hydroelastic analysis of ship structures in Pekeris hydro-acoustic waveguide environment简

The hydroelastic analysis and sonoelastic analysis methods are incorporated with the Green＇s function of the Pekeris ocean hydro-acoustic waveguide model to produce a three-dimensional sonoelastic analysis method for ships in the ocean hydro-acoustic environment. The seabed condition is represented by a penetrable boundary of prescribed density and sound speed. This method is employed in this paper to predict the vibration and acoustic radiation of a 1 500 t Small Water Area Twin Hull （SWATH） ship in shallow sea acoustic environment. The wet resonant frequencies and radiation sound source levels are predicted and compared with the measured results of the ship in trial.     

URANS simulations for a high-speed transom stern ship with breaking waves

An exploratory study of high-speed surface ship flows is performed to identify modelling and numerical issues, to test the predictive capability of an unsteady RANS method for such flows, to explain flow features observed experimentally, and to document results obtained in conjunction with the 2005 ONR Wave Breaking Workshop. Simulations are performed for a high-speed transom stern ship (R/V Athena I) at three speeds Froude number (Fr) = 0.25, 0.43 and 0.62 with the URANS code CFDSHIP-IOWA, which utilizes a single-phase level set method for free surface modelling. The two largest Fr are considered to be high-speed cases and exhibit strong breaking plunging bow waves. Structured overset grids are used for local refinement of the unsteady transom flow at medium speed and for small scale breaking bow and transom waves at high-speeds. All simulations are performed in a time accurate manner and an examination of time histories of resistance and free surface contours is used to assess the degree to which the solutions reach a steady state. The medium speed simulation shows a classical steady Kelvin wave pattern without breaking and a wetted naturally unsteady transom flow with shedding of vortices from the transom corner. At higher speeds, the solutions reach an essentially steady state and display intense bow wave breaking with repeated reconnection of the plunging breaker with the free surface, resulting in multiple free surface scars. The high-speed simulations also show a dry transom and an inboard breaking wave, followed by outboard breaking waves downstream. In comparison to an earlier dataset, resistance is well predicted over the three speeds. The free surface predictions are compared with recent measurements at the two lowest speeds and show good agreement for both non-breaking and breaking waves.     

A numerical study of the turbulent flow around the stern of ship models

The present work is concerned with the numerical calculation of the turbulent flow field around the stern of ship models. The finite volume approximation is employed to solve the Reynolds equations in the physical domain using a body-fitted, locally orthogonal curvilinear co-ordinate system. The Reynolds stresses are modelled according to the standard k-ε turbulence model. Various numerical schemes (i.e. hybrid, skew upwind and central differencing) are examined and grid dependence tests have been performed to compare calculated with experimental results. Moreover, a direct solution of the momentum equations within the near-wall region is tried to avoid the disadvantages of the wall function approach. Comparisons between calculations and measurements are made for two ship models, i.e. the SSPA and HSVA model.     

Testing and calculations of the motion of a model ship with a direct-flow wave propulsor

A new design of the wave propulsor is presented. In this design the thrust mechanism is due to the interaction between the waves and the ship structure elements rather than to ship’s motions. To verify the possibility of using a rigidly fixed inclined plate as the ship wave propulsor a model catamaran was constructed in the Institute of Mechanics of Moscow State University. The effect of the upwave motion of the ship, whose mean velocity is a nonmonotonic function of the wavelength, is studied. As the plate edge pierces the water surface, the ship starts to move in the opposite direction, that is, downwave. The experimentally observable effects are also revealed in the numerical simulation using the XFlow software package which involves the meshless lattice Boltzmann method. On the basis of the calculated results it is shown that the upwave motion effect is due to a variation in the hydrostatic force component in the case of wave breakdown.     

Time series prediction of nonlinear ship structural responses in irregular seaways using a third-order Volterra model

To predict the nonlinear structural responses of a ship traveling through irregular waves, a third-order Volterra model was applied based on the given irregular data. A nonlinear wave–body interaction system was identified using the nonlinear autoregressive with exogenous input (NARX) technique, which is one of the most commonly used nonlinear system identification schemes. The harmonic probing method was applied to extract the first-, second- and third-order frequency response functions of the system. To achieve this, a given set of time history data of both the irregular wave excitation and the corresponding midship vertical bending moment for a certain sea state was fed into the three-layer perceptron neural network. The network parameters are determined based on the supervised training. Next, the harmonic probing method was applied to the identified system to extract the frequency response function of each order. While applying the harmonic probing method, the nonlinear activation function (i.e., the hyperbolic tangent function) was expanded into a Taylor series for harmonic component matching. After the frequency response functions were obtained, the structural responses of the ship under an arbitrary random wave excitation were easily calculated with rapidity using a third-order Volterra series. Additionally, the methodology was validated through the in-depth analysis of a nonlinear oscillator model for a weak quadratic and cubic stiffness term, whose analytic solutions are known. It was confirmed that the current method effectively predicts the nonlinear structural response of a large container carrier under arbitrary random wave excitation.     

Mean velocity and moments of turbulent velocity fluctuations in the wake of a model ship propulsor

Pod drives are modern outboard ship propulsion systems with a motor encapsulated in a watertight pod, whose shaft is connected directly to one or two propellers. The whole unit hangs from the stern of the ship and rotates azimuthally, thus providing thrust and steering without the need of a rudder. Force/momentum and phase-resolved laser Doppler anemometry (LDA) measurements were performed for in line co-rotating and contra-rotating propellers pod drive models. The measurements permitted to characterize these ship propulsion systems in terms of their hydrodynamic characteristics. The torque delivered to the propellers and the thrust of the system were measured for different operation conditions of the propellers. These measurements lead to the hydrodynamic optimization of the ship propulsion system. The parameters under focus revealed the influence of distance between propeller planes, propeller frequency of rotation ratio and type of propellers (co- or contra-rotating) on the overall efficiency of the system. Two of the ship propulsion systems under consideration were chosen, based on their hydrodynamic characteristics, for a detailed study of the swirling wake flow by means of laser Doppler anemometry. A two-component laser Doppler system was employed for the velocity measurements. A light barrier mounted on the axle of the rear propeller motor supplied a TTL signal to mark the beginning of each period, thus providing angle information for the LDA measurements. Measurements were conducted for four axial positions in the slipstream of the pod drive models. The results show that the wake of contra-rotating propeller is more homogeneous than when they co-rotate. In agreement with the results of the force/momentum measurements and with hypotheses put forward in the literature (see e.g. Poehls in Entwurfsgrundlagen für Schraubenpropeller, 1984; Schneekluth in Hydromechanik zum Schiffsentwurf, 1988; Breslin and Andersen in Hydrodynamics of ship propellers, 1996; Schneekluth and Bertram in Ship design for efficiency and economy, 1998), the co-rotating propellers model showed a much stronger swirl in the wake of the propulsor. The anisotropy of turbulence was analyzed using the anisotropy tensor introduced by Lumley and Newman (J Fluid Mech 82(1):161–178, 1977). The invariants of the anisotropy tensor of the wake flow were computed and were plotted in the Lumley–Newman-diagram. These measurements revealed that the anisotropy tensor in the wake of ship propellers is located near to the borders of the invariant map, showing a large degree of anisotropy. They will be presented and will be discussed with respect to applications of turbulence models to predict swirling flows.     

SYNCHRONOUS EFFECT OF SLIPPING HEAVY LOADS ON RO-RO SHIP ROLLING IN WAVES

Common effect of wave and slip of internal vehicles will make rolling of the roll-on ship serious. This is one of the important reasons for overturn of ro-ro ships. The multibody system with a floating base is composed of ro-ro ship and slipping vehicles. Takes the rolling angle of the ship and the transverse displacements of the slipping vehicles on desk as freedoms. Making use of the analysis of apparent gravitation and apparent buoyancy, the wave rolling moment is derived. By means of dynamic method of multibody system, dynamic equations of the system are established. Taking a certain channel ferry as an example, a set of numerical calculation have been carried out for rolling response of the multibody system with a floating base of a ro-ro ship and displacements response of the slipping vehicles under common effect of free slipping vehicles and wave, and a conclusion has been drawn that the motion of the numerous free slipping heavy loads will trend to be synchronous under restraining of the side-wall bulkhead with time because of repeated collision.     

The influence of forward speed on ship motions in abnormal waves: Experimental measurements and numerical predictions

Ship encounters with abnormal waves are increasingly well documented and it is therefore important to be able to model such encounters in order to assess the risks involved and whether there is a requirement for more stringent design rules.This paper presents the results of an experimental investigation into the influence of abnormal waves on a vessel travelling with forward speed in irregular seas. The vessel studied in this case is a naval frigate travelling at a range of speeds. To put the motions measured in abnormal waves into context comparisons are made to those in random seas with a variety of significant wave heights, both non-severe and severe. A further objective is to compare experimental measurements with motion predictions from both a two-dimensional linear strip theory and a three-dimensional partly nonlinear seakeeping model.Results demonstrate that abnormal waves may not necessarily be the worst conditions that a ship can encounter. However, accelerations derived from the rigid body motions appear to be in excess of rules values. This has implications for design due to the unexpected nature of abnormal wave occurrence and the consequent likelihood of encountering such a wave at a higher speed (hence in a more severe operating condition) than a random sea of an equivalent height.The three-dimensional partly nonlinear model demonstrates improved agreement with experimental measurements of rigid body motions, compared to the two-dimensional strip theory. It is therefore considered to have greater potential as a design tool for abnormal wave encounters. Further validation with a wide range of sea states and vessel types is required.     

绕椭球的低速流动研究

理解和预测绕椭球的流动对指导飞行器和潜艇等交通工具的设计具有很强的工程意义. 近年来, 针对椭球绕流开展了大量的实验和数值模拟研究. 对有攻角下椭球绕流分离的定性描述和定量研究, 促进了对三维分离的辨识和拓扑研究. 文章对流场特性进行了分析, 介绍了分离对气动力、噪声、尾迹的影响, 以及实验条件对流动的影响. 上述定常流动与非定常机动过程之间存在明显差异, 非定常机动过程不能作为定常或准定常问题处理, 在机动过程中, 分离出现明显延迟, 气动力出现明显变化. 随后介绍了数值模拟在求解绕椭球流动中的进展, 当前求解雷诺平均的N-S方程湍流模式仍然是解决绕椭球大范围分离流动的主要工程方法, 大涡模拟和分离涡模拟等也逐渐得到了广泛应用. 受限于计算能力, 直接数据模拟只能用于较低雷诺数, 在高雷诺数流动中还不适用. 非定常机动过程的数值模拟较定常状态, 与实验结果的差距要大一些. 最后, 介绍了对椭球绕流场转捩的研究进展, 对T-S转捩与横流转捩的机理和辨识已经较为准确, 数值模拟结果与实验结果基本相符, 但对再附转捩的认识还不够清晰, 尤其是迎风面, 因此椭球绕流转捩的研究还需要依靠实验.