A 3D fully coupled analysis of nonlinear sloshing and ship motion

This study investigates the coupling effects of six degrees of freedom in ship motion with fluid oscillation inside a three-dimensional rectangular container using a novel time domain simulation scheme. During the time marching, the tank-sloshing algorithm is coupled with the vessel-motion algorithm so that the influence of tank sloshing on vessel motions and vice versa can be assessed. Several factors influencing the dynamic behavior of tank–liquid system due to moving ship are also investigated. These factors include container parameters, environmental settings such as the significant wave height, current velocity as well as the direction of wind, wave and flow current acting on the ship. The nonlinear sloshing is studied using a finite element model whereas nonlinear ship motion is simulated using a hybrid marine control system. Computed roll response is compared with the existing results, showing fair agreement. Although the two hull forms and the sea states are not identical, the numerical result shows the same trend of the roll motion when the anti-rolling tanks are considered. Thus, the numerical approach presented in this paper is expected to be very useful and realistic in evaluating the coupling effects of nonlinear sloshing and 6-DOF ship motion.     

双浮板液舱晃荡特性的数值研究

基于黏性流理论,采用动网格技术和6自由度模型,以及动量源方法,建立了双浮板液舱晃荡的数值模型。分别采用3种不同空间步长的网格离散计算区域,进行了网格收敛性验证。通过光滑液舱晃荡的模型试验和解析得到的爬高最大值,验证了数值模型的精确性。在载液率为50%,激励幅值为2 mm条件下,对双浮板液舱晃荡进行了数值计算,与光滑液舱相比,双浮板液舱晃荡的最大爬高明显减小。通过一个激励周期内双浮板液舱晃荡的波面显示发现,液舱晃荡模式由光滑液舱的驻波模式变为U管模式,晃荡模式的改变起到了明显地抑制液舱晃荡的效果。     

The effects of LNG-tank sloshing on the global motions of LNG carriers

The coupling and interactions between ship motion and inner-tank sloshing are investigated by a time-domain simulation scheme. For the time-domain simulation, the hydrodynamic coefficients and wave forces are obtained by a potential-thoery-based three-dimensional (3D) diffraction/radiation panel program in frequency domain. Then, the corresponding simulations of motions in time domain are carried out using convolution integral. The liquid sloshing in a tank is simulated in time domain by a Navier–Stokes solver. A finite difference method with SURF scheme is applied for the direct simulation of liquid sloshing. The computed sloshing force and moment are then applied as external excitations to the ship motion. The calculated ship motion is in turn inputted as the excitation for liquid sloshing, which is repeated for the ensuing time steps. For comparison, we independently developed a coupling scheme in the frequency domain using a sloshing code based on the linear potential theory. The hydrodynamic coefficients of the inner tanks are also obtained by a 3D panel program. The developed schemes are applied to a barge-type FPSO hull equipped with two partially filled tanks. The time-domain simulation results show similar trend when compared with MARIN's experimental results. The most pronounced coupling effects are the shift or split of peak-motion frequencies. It is also found that the pattern of coupling effects between vessel motion and liquid sloshing appreciably changes with filling level. The independent frequency-domain coupled analysis also shows the observed phenomena.     

PSME model of parametric excitation of two-layer liquid in a tank

Internal waves driven by external excitation constitute important phenomena that are often encountered in environmental fluid mechanics. In this study, a pseudospectral σ-transformation model is used to simulate parametric excitation of stratified liquid in a two-layer rectangular tank. The σ-transformation maps the physical domain including the liquid free surface, the interface between the liquid layers, and the bed, onto a pair of fixed rectangular computational domains corresponding to the two layers. The governing equation and boundary conditions are discretised using Chebyshev collocation formulae. The numerical model is verified for two analytical sloshing problems: horizontal excitation of constant density liquid in a rectangular tank, and vertical excitation of stratified liquid in a rectangular tank. A detailed analysis is provided of liquid motions in a shallow water tank due to excitations in the horizontal and the vertical directions. Also, the effect of pycnocline on the wave motions and patterns is studied. It is found that wave regimes and patterns are considerably influenced by the pycnocline, especially when the excitation frequency is large. The present study demonstrates that a pseudospectral σ-transformation is capable to model non-linear sloshing waves in a two-layer rectangular tank.     

Study on coupling effects of ship motion and sloshing

This study considers the coupling effects of ship motion and sloshing. The linear ship motion is solved using an impulse-response-function (IRF) method, while the nonlinear sloshing flow is simulated using a finite-difference method. The IRF method requires the frequency-domain solution prior to conversion to time domain, but the computational effort is much less than that of direct time-domain approaches. The developed scheme is verified by comparing the motion RAOs between the frequency-domain solution and the solution obtained by the IRF method. Furthermore, a soft-spring concept and linear roll damping are implemented to predict more realistic motions of surge, sway, yaw, and roll. For the simulation of sloshing flow in liquid tanks, a physics-based numerical approach adopted by Kim [2001. Numerical simulation of sloshing flows with impact load. Applied Ocean Research 23, 53–62] and Kim et al. [2004. Numerical study on slosh-induced impact pressures on three-dimensional prismatic tanks. Applied Ocean Research 26, 213–226] is applied. In particular, the present method focuses on the simulation of the global motion of sloshing flow, ignoring some local phenomena. The sloshing-induced forces and moments are added to wave-excitation forces and moments, and then the corresponding body motion is obtained. The developed schemes are applied for two problems: the sway motion of a box-type barge with rectangular tanks and the roll motion of a modified S175 hull with rectangular anti-rolling tank. Motion RAOs are compared with existing results, showing fair agreement. It is found that the nonlinearity of sloshing flow is very important in coupling analysis. Due to the nonlinearity of sloshing flow, ship motion shows a strong sensitivity to wave slope.     

Sloshing waves and resonance modes of fluid in a 3D tank by a time-independent finite difference method

A 3D time-independent finite difference method is developed to solve for wave sloshing in a three-dimensional tank excited by coupled surge and sway motions. The 3D equations of fluid motion are derived in a moving coordinate system. The three-dimensional tank, with an arbitrary depth and a square base, is subjected to a range of excitation frequencies with motions that exhibit multiple degrees of freedom. For demonstration purposes the numerical scheme is validated by a benchmark study. Five types of sloshing waves were observed when the tank is excited by various excitation frequencies. A spectral analysis identified the resonant frequencies of each type of wave and the results show a strong correlation between resonant modes and the occurrence of the sloshing wave types. The method can be used to simulate fluid sloshing in a 3D tank with six-degrees of freedom.     

Experimental Study on A Sloshing Mitigation Concept Using Floating Layers of Solid Foam Elements

A sloshing mitigation concept taking advantage of floating layers of solid foam elements is proposed in the present study. Physical experiments are carried out in a liquid tank to investigate the hydrodynamic mechanism of this concept. Effects of the foam-layer thickness, excitation amplitude, and excitation frequency on the sloshing properties are analyzed in detail. It is found that the floating layers of solid foam elements do not evidently affect the fundamental natural sloshing frequency of the liquid tank evidently among the considered cases. At the resonant condition, the maximum wave height and dynamic pressure are greatly reduced as the foam-layer thickness increases. Higher-order pressure components on the tank side gradually vanish with the increase of the foam-layer thickness. Cases with different excitation amplitudes are also analyzed. The phenomenon is observed when the wave breaking in the tank can be suppressed by solid foam elements.     

Motion responses of barge carrying liquid tank

The interaction between the liquid sloshing in a rectangular tank equipped inside the barge and the barge responses has been investigated through a comprehensive experimental program. The barge was subjected to both regular and random wave excitations under beam sea condition. Three relative fill levels (h_s/l) with liquid fill depth (h_s) to length of tank (l) ratio of 0.163, 0.325 and 0.488 were considered. In addition, the barge responses of equivalent dry weight condition corresponding to each fill level were measured to understand the influence of sloshing. While the excitation wave frequency equals to first mode natural sloshing frequency, a noticeable decrease in the sway response has been observed. However, the effect of sloshing oscillation on the heave response is insignificant. A split up of roll resonance was observed for the aspect ratio of 0.163 due to the coupling effect of roll motion and sloshing.     

Experimental Studies on Sloshing Mitigation Using Dual Perforated Floating Plates in A Rectangular Tank

The performance of dual perforated floating plates in a rectangular tank is investigated based on the model tests under different external excitations for different filling rates.It is found that dual perforated floating plates in the tank can remarkably mitigate violent resonant sloshing responses compared with the clean tank,especially when the external excitation frequency is in the vicinity of the first-order resonant frequency.Next,the parametric studies based on different filling rates and external excitation amplitudes are performed for the first-order resonant frequencies.The presence of dual perforated floating plates seldom shifts the sloshing natural frequencies.Further,dual perforated floating plates change the sloshing modes from the standing-wave mode in the clean tank to the Utube mode,which can arise from the sloshing reduction to some extent.     

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.     

Experimental and numerical investigation of sloshing using different free surface capturing methods

We investigated the use of numerical methods to predict liquid sloshing phenomena in a moving tank and compared our results to model test measurements. The numerical techniques for the free surface, based on the so-called finite Volume-of-Fluid (VoF) approach, comprised an incompressible VoF method, an incompressible coupled Level-Set and Volume-of-Fluid (clsVoF) method, and a compressible VoF method. We assessed the capability of these three numerical methods to achieve suitable numerical predictions of sloshing phenomena, specifically, air pockets and bubbles on the free surface inside a test tank. To observe the described sloshing phenomena, we simulated tank motions leading to well defined single impact wave motions. We performed repeated physical tests for validation purposes. Computed velocity and pressure time histories were compared to experimental data we obtained from Particle Image Velocimetry (PIV) and pressure sensor measurement. Grid sensitivity and turbulence model studies were performed. We demonstrated that the compressible VoF method was the most suitable method to obtain accurate predictions of sloshing phenomena.     

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.     

Three-dimensional liquid sloshing in a tank with baffles

A numerical model has been developed to study three-dimensional (3D) liquid sloshing in a tank with baffles. The numerical model solves the spatially averaged Navier-Stokes equations, which are constructed on a non-inertial reference frame having six degree-of-freedom (DOF) of motions. The large-eddy-simulation (LES) approach is employed to model turbulence by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. The baffles in the tank are modeled by the concept of virtual boundary force (VBF) method. The numerical model is first validated against the available analytical solution and experimental data for two-dimensional (2D) liquid sloshing in a tank without baffles. The 2D liquid sloshing in tanks with baffles is then investigated. The numerical results are compared with other results from available literatures. Good agreement is obtained. Finally, the model is used to study 3D liquid sloshing in a tank with vertical baffles. The effect of the baffle is investigated and discussed.     

液舱三维晃荡运动二阶共振理论解及特性分析

从无旋运动的理论出发,并利用微扰法,推导了液舱三维晃荡运动二阶共振问题的理论解。考虑纵荡和横荡运动情况,对液舱三维晃荡二阶共振问题进行了分析。当两个晃荡方向的和频(即其外部激发频率的和)或差频(即其外部激发频率的差值)等于液舱固有频率时,二阶共振发生;当某一晃荡方向(横荡或纵荡)外部激发频率与另一晃荡方向(纵荡或横荡方向)液舱某一固有频率的和或差值等于液舱另一固有频率时,二阶共振也会发生。进一步研究了各个二阶共振激发频率下水深变化对晃荡振幅的影响。结果表明,对于两个晃荡方向外部激发频率的和频和单一晃荡方向(纵荡或横荡)某一个激发频率与另一晃荡方向(横荡或纵荡)某一个属于奇模的固有频率的和频所引发的共振情况,水深变化对共振振幅大小的影响比较大;而对于相应差频所引发的共振情况,水深变化对共振振幅大小的影响比较小。     

A nonlinear model of parametric rolling stabilization by anti-roll tanks

The present paper describes a mathematical model in which the fluid motion inside a U-tank is nonlinearly coupled to the heave, roll and pitch motions of the ship. The main purpose of the investigation is centred on the control of roll motion in the case of parametric resonance in longitudinal waves. A transom stern small vessel, known to be quite prone to parametric amplification, is employed in the study. Four tank designs are employed in order to study the influence of tank mass, tank natural frequency and tank internal damping on the control of parametric rolling at different head seas conditions. Additionally, the influence of the vertical position of the tank is also investigated. The main results are presented in the form of limits of stability, with encounter frequency and wave amplitudes as parameters. Distinct dynamical characteristics are discussed and conclusions are drawn on the relevant parameters for the efficient control of the roll amplifications in head seas.     

基于HUST-Ship的船舶耐波性数值试验水池参数化建模方法研究

基于自研的HUST-Ship黏性流CFD求解器,提出了一种船舶耐波性数值试验水池参数化建模方法和网格生成技术。在自研的CFD软件平台中,根据波浪参数(规则波的波长、波高和不规则波的特征周期、有效波高)自动生成船舶耐波性数值试验水池。以结构化重叠网格技术为基础,分别对船舶与水池进行网格划分并进行网格的组装、重叠,实现船舶在耐波性数值试验水池中的数值“试验”,对船舶耐波性进行数值预报。根据波浪种类分为规则波数值试验水池与不规则波数值试验水池,可对各类波浪条件实现精确可靠的数值造波,模拟船舶在各种波浪环境中运动响应、耐波性事件与非线性运动特性。采用参数化建模方法生成的耐波性数值试验水池能够提高建模和网格生成效率,通过对国际标模进行数值预报并与船模试验结果对比,验证了参数化建模方法生成的耐波性数值试验水池对船舶在波浪中运动性能的计算精度。     

三维完全非线性波浪水槽的数值模拟

用有限元求解拉普拉斯方程,建立了三维完全非线性数值波浪水槽.跟踪流体自由表面的方法为满足完全非线性自由表面条件的半拉格朗日法,对离散单元采用20节点的六面体二次等参数单元.并把数值计算结果与水面初始升高产生箱体内流体运动解析解和二阶斯托克斯波理论解进行了对比,结果表明该模型是稳定的、守恒的,能精确模拟非线性波浪的产生和传播.     

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.     

Design of passive anti-roll tanks for roll stabilization in the nonlinear range

The best way of reducing roll motion is by increasing roll damping. Bilge keels are the most common devices for increasing roll damping. If more control is required, anti-roll tanks and fins are used. Tanks have the advantage of being able to function when the ship is not underway. Our objective is to develop design procedures for passive tanks for roll reduction in rough seas. This paper focuses on the design of passive U-tube tanks. The tank-liquid equation of motion is integrated simultaneously with the six-degree-of-freedom (6DOF) equations of the ship motion. The coupled set of equations is solved by using the Large Amplitude Motion Program ‘LAMP’, which is a three-dimensional time-domain simulation of the motion of ships in waves. The unstabilized and stabilized roll motions of a S60-70 ship with forward speed and beam waves have been analyzed. For high-amplitude waves, the unstabilized roll angle exhibits typical nonlinear phenomena: a shift in the resonance frequency, multi-valued responses, and jumps. The performance of a S60-70 ship with a passive tank is investigated in various sea states with different encounter wave directions. It is found that passive anti-roll tanks tuned in the linear or nonlinear ranges are very effective in reducing the roll motion in the nonlinear range. The effect of the tank damping, frequency, and mass on the tank performance is studied. Also, it is found that passive anti-roll tanks are very effective in reducing the roll motion for ships having a pitch frequency that is nearly twice the roll frequency in sea states 5 and 6.     

港口内靠码头系泊船运动的计算

本文以英国南海岸处Ｓｈｏｒｅｈａｍ港内系泊船为例，研究了港口内靠码头系泊船运动的数值计算问题。该船的实船实验和模型实验已经完成，为验证本文计算结果提供了依据。运动方程在时域内求解。在计算船体脉冲响应函数时，引入了船体阻尼系数在频率很低时的渐近表达式和一迭代算法。分析了港内共振波浪和其对船体的作用力。讨论了靠码头系泊船运动的非线性特征。计算结果与实验结果符合很好。