Towards the modeling of the ditching of a ground-effect wing ship within the framework of the SPH method

Ditching often takes place for a ground-effect wing (WIG) ship. During the ditching, the extreme load developed by water impacts may cause damages to structures, posing a great threat to the safety of crew and passengers. In the paper, a weakly compressible smoothed particle hydrodynamics (SPH) model combined with enhanced numerical techniques has been adopted to tackle the ditching problems. In order to handle the motion of a rigid body in the three-dimensional ditching problems, the six degrees-of-freedom (6-DOF) equations of motion are incorporated into the SPH scheme. The accuracy of the SPH model is validated through two benchmarks, respectively, the two-dimensional wedge water entry and the three-dimensional stone-skipping. The former is aimed to validate the prediction of pressures during the free surface impact while the latter is a good case for testing the coupling motions of the rigid body. Furthermore, the ditching of the real scale WIG ship under different conditions is simulated with the established SPH model, through which some useful conclusions are drawn.     

基于无网格MPS方法数值分析方舱破舱进水问题

采用移动粒子半隐式法(moving particle semi-implicit,简称MPS)对自由漂浮二维方舱的破舱瞬时进水过程进行数值模拟。首先,采用基于GPU平台自主开发的MPS软件模拟破舱进水问题,并与其他方法得到的数值模拟结果进行对比验证。然后,对该二维方舱的各种模型进行了数值模拟,其中开孔位于不同位置以表示舷侧不同高度下的损坏。此外,还研究了不同类型的挡板对破舱进水后方舱稳定性的影响。结果表明损坏的孔洞和内部挡板会影响损坏舱段的运动特性,开孔距静水面的距离越大引起舱段的横摇等运动幅度越大,垂直挡板比水平挡板对舱内洪水的影响更大。     

Numerical analysis on ship maneuvering coupled with ship motion in waves

This paper considers a numerical analysis of ship maneuvering performance in the presence of incident waves and resultant ship motion responses. To this end, a time-domain ship motion program is developed to solve the wave–body interaction problem with the ship slip speed and rotation, and it is coupled with a modular-type 4-DOF maneuvering problem. In this coupled problem, the second-order mean drift force, which can play an important role in the ship maneuvering trajectory, is estimated by using a direct pressure integration method. The developed method is validated by observing the second-order mean drift force, and planar trajectories in maneuvering tests with and without the presence of incident waves. The comparisons are made for two ship models, Series 60 with block coefficient 0.7 and the S-175 containership, with existing experimental data. The maneuvering tests observed in this study include a zig-zag test in calm water, and turning tests in calm water and in regular waves. The present results show a fair agreement of overall tendency in maneuvering trajectories.     

系泊船非线性波浪力时域计算:二维模型

为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算.     

Numerical simulations of sloshing flows with an elastic baffle using a SPH-SPIM coupled method

In this work, sloshing flows were successfully simulated by using a coupled numerical scheme between smoothed particle hydrodynamics (SPH) and smoothed point interpolation method (S-PIM) (SPH-SPIM coupled method). SPH is a Lagrangian particle method to solve flow fields while S-PIM is developed to deal with the structure dynamics. A coupling scheme is proposed, the key of which is that the fluid and solid fields are not necessary to be discretized by the same resolution. The stability, accuracy, convergence and conservation of the SPH-SPIM coupled method were validated by the case of hydrostatic water column on an elastic plate. Then, a wave impact problem was simulated to verify that the present SPH method worked well for sloshing flows. Finally, two sloshing problems with an elastic baffle were simulated, which validated the accuracy and stability of the method in predicting the fluid-structure interaction (FSI) features during the process of sloshing. It has been found that both the shape of the free surface and the large deformation of the elastic baffle can be well captured by the present method, which shows the potential of the present method to be a good candidate for simulating sloshing problems.     

An experimental investigation into the influence of the damage openings on ship response

Ship motions after damage are difficult to evaluate since they are affected by complex phenomena regarding fluid and structures interactions. The possibility to better understand how ship behavior in damage is influenced by these phenomena is important for improving ship safety, especially for passenger vessel.In this paper an experimental campaign is carried out on a passenger ferry hull, to show the effects of the water dynamics across damage openings on ship motions. Novel aspects of this research include the study of the effects of the damage position on the ship roll response. The study is carried out for still water and for beam regular waves at zero speed.Results from the experiments carried out underline that the roll behavior of a damaged ship is affected by the position of damage opening and not only by its size. Assuming the same final equilibrium conditions after flooding but characterized by different damage openings it is possible to observe how motions RAOs and roll decay characteristics modify according to the opening locations.     

Study on Ship Automatic Berthing System with Mooring Lines

This study describes an automatic berthing system with mooring lines. It is designed to be berthed by using mooring device on the upper deck of a ship. It is to berth once maintaining parallel with the quay by controlling both forward and aft breast lines. Berthing method is used through length adjustment of mooring lines connected between ship and quay by controlling the angular velocity and the torque of hydraulic motor in mooring device. The study is conducted under three changing conditions of draft, such as even-keel, rise of the gravity center and trim to stern. Variables affecting berthing stability are determined based on the control performance of each condition. Bond graphs method is used to model the system. Controller is designed as PID control method of reference-model algorithm. The control program is composed of synchronous control system based on the equations derived with the numerical analysis. The tank test is conducted to verify the usefulness of the control program.     

Numerical effects on organic-matter sedimentation and remineralization in biogeochemical ocean models

We investigate the effects of different vertical grid resolutions and algorithms for the calculation of particle sinking on the sedimentation and remineralization of particulate organic matter. Simulations carried out with an idealized 1D model of detritus sinking show that a coarse vertical resolution, such as used in many global biogeochemical models, tends to enhance the particle flux through numerical mixing within the vertical boxes, and thereby simulates deeper remineralization, compared to a model with a fine vertical resolution. This effect can be ameliorated by assuming a distribution of detritus within the individual grid boxes that corresponds to the prescribed sinking and remineralization parameters. Experiments of the different flux algorithms, carried out with 3D global biogeochemical models of different vertical grid resolution reveal impacts on simulated biogeochemical tracer distributions that are similar to those obtained by substantial variations in biogeochemical model parameters. Our results indicate that numerical schemes have to be considered when comparing biogeochemical parameter values of different models and also when porting biogeochemical models among different circulation models.     

Wash Waves Generated by High Speed Displacement Ships

It is difficult to compute far-field waves in a relative large area by using one wave generation model when a large calculation domain is needed because of large dimensions of the waterway and long distance of the required computing points. Variation of waterway bathymetry and nonlinearity in the far field cannot be included in a ship fixed process either. A coupled method combining a wave generation model and wave propagation model is then used in this paper to simulate the wash waves generated by the passing ship. A NURBS-based higher order panel method is adopted as the stationary wave generation model; a wave spectrum method and Boussinesq-type equation wave model are used as the wave propagation model for the constant water depth condition and variable water depth condition, respectively. The waves calculated by the NURBS-based higher order panel method in the near field are used as the input for the wave spectrum method and the Boussinesq-type equation wave model to obtain the far-field waves. With this approach it is possible to simulate the ship wash waves including the effects of water depth and waterway bathymetry. Parts of the calculated results are validated experimentally, and the agreement is demonstrated. The effects of ship wash waves on the moored ship are discussed by using a diffraction theory method. The results indicate that the prediction of the ship induced waves by coupling models is feasible.     

A Study of the Stability Properties in Simulation of Wave Propagation with SPH Method

When ordinary Smoothed Particle Hydrodynamics (SPH) method is used to simulate wave propagation in a wave tank, it is usually observed that the wave height decays and the wave length elongates along the direction of wave propagation. Accompanied with this phenomenon, the pressure under water decays either and shows a big oscillation simultaneously. The reason is the natural potential tensile instability of modeling water motion with ordinary SPH which is caused by particle negative stress in the computation. To deal with the problems, a new sextic kernel function is proposed to reduce this instability. An appropriate smooth length is given and its computation criterion is also suggested. At the same time, a new kind dynamic boundary condition is introduced. Based on these improvements, the new SPH method named stability improved SPH (SISPH) can simulate the wave propagation well. Both the water surface and pressure can be well expressed and the oscillation of pressure is nearly eliminated. Compared with other improved methods, SISPH can truly reveal the physical reality without bringing some new problems in a simple way.     

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.     

Experimental studies of a damaged ship section in forced heave motion

This work documents a detailed series of experiments performed in a wave flume on a thin walled prismatic hull form. The model consists of a rectangular opening located on the side. The length of the model is slightly smaller than the flume breadth to achieve two-dimensional (2D) behavior in the experiments. Forced oscillatory heave tests in calm water have been carried out by varying the model-motion parameters and examining both intact and damaged conditions. Video recordings, measurements of the wave elevation inside the damaged compartment and of the force on the model were performed in all the experiments. The effect of damage opening in the model on hydrodynamic loads is examined by comparing with an intact section. A theoretical analysis is used to explain the behavior of added mass and damping coefficients in heave for a 2D damaged section. The presented results demonstrate occurrence of sloshing and piston mode resonances in the tests and their influence on the hydrodynamics loads of a damaged ship. Detailed physical investigations are presented at these resonance frequencies for the damaged section. Effect of filling level in the damage compartment, damage-opening length and air compressibility in the airtight compartment is examined. Nonlinear effects are documented and appear dominant, especially, for lowest filling level where we have shallow-water depth conditions in the damaged compartment. Resonance phenomena that can lead to significant local loads are identified for the shallow water condition. Air compressibility in the airtight compartment and floodwater act as a coupled system and influence inflow/outflow of floodwater in the compartment. It has a significant effect on local floodwater behavior in the damaged compartment.     

Study on the uncertainty of the available time under ship fire based on Monte Carlo sampling method

Available safety egress time under ship fire(SFAT) is critical to ship fire safety assessment,design and emergency rescue.Although it is available to determine SFAT by using fire models such as the two-zone fire model CFAST and the field model FDS,none of these models can address the uncertainties involved in the input parameters.To solve this problem,current study presents a framework of uncertainty analysis for SFAT.Firstly,a deterministic model estimating SFAT is built.The uncertainties of the input parameters are regarded as random variables with the given probability distribution functions.Subsequently,the deterministic SFAT model is employed to couple with a Monte Carlo sampling method to investigate the uncertainties of the SFAT.The Spearman’s rank-order correlation coefficient(SRCC) is used to examine the sensitivity of each input uncertainty parameter on SFAT.To illustrate the proposed approach in detail,a case study is performed.Based on the proposed approach,probability density function and cumulative density function of SFAT are obtained.Furthermore,sensitivity analysis with regard to SFAT is also conducted.The results give a high-negative correlation of SFAT and the fire growth coefficient whereas the effect of other parameters is so weak that they can be neglected.     

Flow computations for ship safety problems

A number of ship safety problems involve flow computations. Most of them do not require involved CFD calculations, but can be solved by relatively simple means which, nonetheless, are not state of the art. The problems discussed here are, in order of increasing sophistication: cross-flooding of damaged compartments; evacuation of persons on board a ship; sinking of a damaged ship in still water; accelerations and loads on free-fall lifeboats and damaged ship survival time in a seaway.     

A Corrected Solid Boundary Treatment Method for Smoothed Particle Hydrodynamics

Smoothed Particle Hydrodynamics method (SPH) has a good adaptability for simulating of free surface flow problems. However, there are some shortcomings of SPH which are still in open discussion. This paper presents a corrected solid boundary handling method for weakly compressible SPH. This improved method is very helpful for numerical stability and pressure distribution. Compared with other solid boundary handling methods, this corrected method is simpler for virtual ghost particle interpolation and the ghost particle evaluation relationship is clearer. Several numerical tests are given, like dam breaking, solitary wave impact and sloshing tank waves. The results show that the corrected solid boundary processing method can recover the spurious oscillations of pressure distribution when simulating the problems with complex geometry boundary.     

A time-domain method to evaluate the use of moving weights to reduce the roll motion of a ship

A method to evaluate the use of actively controlled moving weights on board ships to reduce roll motion is developed. The weights can simulate in principle anti-roll-tank systems, or they can be considered a possible anti-roll device in their own right. The ship, the moving weight, and the control device are considered components of a single dynamic system. The full eight-degree-of-freedom set of coupled governing equations for the complete dynamic system is derived. And a three-degree-of-freedom non-linear approximation for the roll motion only (MOTSIM) is derived from these eight equations. The reduced set of equations is used to determine the influence of various parameters and to evaluate control strategies. A PID controller is developed to command the position of the weight and a servomechanism model is used to predict its actual position. Then, the moving-weight system is incorporated into LAMP (Large–Amplitude–Motion Program), a computer code that integrates the governing equations of the sea and the motion of the ship interactively and simultaneously and predicts the motion of the ship in the time domain. A comparison of the results from the two simulations shows that there is fairly good correlation between the simple and complex models, but the simple model is a little optimistic in predicting the effectiveness of the moving-weight system. The results predict that the moving-weight system can be an excellent roll-suppressing device when the moving weight is as small as 1% of the displacement of the ship and the maximum distance the weight moves is as small as 15% of the half-beam.     

撞击载荷作用下水下沉箱防护设施损伤分析

面对通航区很可能发生的船舶抛锚、落物以及沉船等事故造成的危险,使用新型沉箱防护设施对水下生产系统进行安全防护,但目前撞击载荷下沉箱防护设施的安全评估标准并无明确界定。采用有限元方法对水下生产系统沉箱安全防护设施的撞击损伤进行了深入研究,考虑不同落物形状、下落速度、落物质量、碰撞位置对沉箱防护设施撞击损伤的影响。结果表明,锚、长方体落物对沉箱撞击最严重;撞击速度和落物质量与撞击应力成正比;对于非对称开孔的沉箱顶盖,距离开孔较近区域为应力敏感区,受到撞击时撞击应力较大。     

A Time-Domain Coupled Model for Nonlinear Wave Forces on A Fixed Box-Shaped Ship in A Harbor

A 2-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear wave forces on a fixed box-shaped ship in a harbor.The domain is divided into an inner domain and an outer domain.The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations.The other area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions.Along the interface boundaries between the inner domain and the outer domain,the volume flux is assumed to be continuous and the wave pressures are equal.Relevant physical experiment is conducted to validate the present model.It is shown that the numerical results agree with the experimental data.Compared with the coupled model with the flow in the inner domain governed by the Laplace equation,the present coupled model is more efficient and its solution procedure is more simple,which is particularly useful for the study on the effect of the nonlinear wave forces on a fixed box-shaped ship in a large harbor.     

砾石单元法与SPH耦合模型在数值波浪水槽中的应用研究

为研究以流体粒子描述波浪运动,以固体单元描述砾石运动的两相介质大变形运动,在港口、海岸工程科学研究中具有重要意义。本文提出砾石单元法（GEM）,介绍了光滑粒子动力学方法（SPH）和GEM的基本原理,阐述了GEM与离散单元法（DEM）的异同之处,说明了采用SPH方法与GEM构建波浪砾石耦合运动数学模型的方法和过程。应用SPH方法建立数值波浪水槽,用GEM模拟波浪作用下堆积砾石的滚落、坍塌变形,构建了SPH方法与GEM耦合数学模型。模拟了水槽造波和波浪生成过程和波浪作用下砾石的滚落、坍塌变形,并与物理模型试验成果进行了比较,结果基本吻合。本文提出的GEM法具有模拟单相堆积砾石运动和堆积砾石与流体粒子耦合多相介质运动的功能,是对DEM法的补充和改善。本文提出的堆积力学球概念和拟序排列求解方法是砾石单元法的重要组成部分。     

南海北部海洋生态模型的参数分析及遗传算法优化

海洋生态系统动力学模型是研究海洋生态环境的重要手段。随着模型的发展, 生态参数取值不确定性增加, 对模型结果的影响逐渐增大, 因此模型参数优化显得尤为重要。本研究在南海北部应用一维物理-生态耦合模型, 通过对模型生态参数进行敏感性分析, 获取关键生态参数, 利用遗传算法对参数进行优化。结果表明, 模型中的敏感参数主要集中于浮游植物生长和浮游动物生长、摄食和死亡以及碎屑沉降等过程。针对以上参数利用遗传算法优化, 发现仅加入表层卫星数据, 模型表层和垂向模拟误差分别降低27.80%和21.40%; 加入垂向观测数据, 表层和垂向模拟误差分别降低14.90%和32.70%。遗传算法应用于海洋生态模型的关键参数优化研究, 所获取的参数对模型有明显的改善效果, 提高了耦合模型对生态系统的模拟精度, 为参数优化在三维模型中的应用提供了依据。