On the analysis of sloshing of water in rectangular containers: Numerical study and experimental validation

In this work the analysis of sloshing of water in rectangular open tanks has been extensively carried out. Two mathematical models are employed, respectively the Reynolds Averaged Navier Stokes Equations (RANSE) and the Shallow Water Equations (SWE). The RANSE are solved using a modified form of the well established MAC method (SIMAC) able to treat both the free surface motion and the viscous stresses over the rigid walls accurately. The Shallow Water Equations are solved by means of a simple and powerful algorithm (CE-SE) able to deal with large impacting waves over the tank walls.Successively, in order to validate the mentioned algorithms and for a better understanding of the sloshing phenomenon, experimental tests have been carried out using a 0.5 m breadth rectangular tank in periodic roll motion.It has been shown that RANSE provide more accurate solutions than SWE for small or moderate amplitudes of excitation. In particular in this paper it is proved that the shallow water approximation can be efficiently adopted within liquid depth to tank breadth RATIO = 0.15, when examining the sloshing problem. By increasing the water level inside the tank, results by SWE show large qualitative and quantitative disagreement with experiments. Nevertheless, in the case of large amplitude excitation, when sprays and large breaking waves are expected, SWE provide a fairly good estimate of the sloshing induced waves.Finally a simple baffle configuration inside the tank has been considered. By the analysis of numerical results, it has been observed that the presence of a vertical baffle at the middle of the tank dramatically changes the sloshing response compared to the unbaffled configuration. It produces a jump-like effect, resulting in a weak magnification of the dynamic loads on the vertical walls out of resonance, and a strong reduction of the dynamic loads in the resonance condition.     

MPS方法数值模拟液舱晃荡问题

基于无网格粒子法MPS方法(moving particle semi-implicit method)研究了液舱晃荡问题。针对二维矩形液舱晃荡问题进行了数值验证,结果表明MPS方法能够很好地计算晃荡产生的拍击压力。同时将MPS方法应用到带隔板的液舱晃荡问题计算中,分析了二维和三维带隔板液舱晃荡问题。计算结果表明:隔板的存在很大程度地限制了流体的水平运动,隔板附近出现了自由面的翻卷、破碎和融合现象,MPS方法能够很好地模拟这些流动现象。计算得到的波高与实验测得的波高吻合较好,表明MPS方法模拟带隔板的晃荡问题具有一定的可靠性。     

Effects of perforated baffle on reducing sloshing in rectangular tank: Experimental and numerical study

A liquid sloshing experimental rig driven by a wave-maker is designed and built to study liquid sloshing problems in a rectangular liquid tank with perforated baffle. A series of experiments are conducted in this experimental rig to estimate the free surface fluctuation and pressure distribution by changing external excitation frequency of the shaking table. An in-house CFD code is also used in this study to simulate the liquid sloshing in three-dimensional （3D） rectangular tank with perforated baffle. Good agreements of free surface elevation and pressure between the numerical results and the experimental data are obtained and presented. Spectral analysis of the time history of free surface elevation is conducted by using the fast Fourier transformation.     

Boundary element analysis of liquid sloshing characteristics in axisymmetric tanks with various porous baffles

This paper aims to investigate the effects of the porous baffles on the suppression of sloshing for the tanks with axisymmetric geometries under lateral excitation. Based on the assumptions of inviscid, irrotational, incompressible liquid and small amplitude sloshing, an axisymmetric boundary element method (BEM) for 3D Laplace equation is derived by using the Green's theorem together with the weighted residual method. And a zoning method is employed to model fluid domain in the tanks with complex porous baffles. Meanwhile, the porous baffles are treated motioning together with the tanks, and the velocity across the porous baffle is assumed to be linearly proportional to the pressure gradient between each side of the porous baffle. And the mechanism of suppressing the sloshing response is mainly the energy dissipation of the fluid passing through the porous baffle. Moreover, the linear free surface boundary conditions are also used to solve the governing equations. Compared with other numerical methods, the most prominent advantage of the BEM in solving axisymmetric potential problem is that only the boundaries of half the cross-section instead of the entire problem domain should be discretized, which can cut down large amount of memory and time costs. The present method is verified by comparing the numerical results with the existing literatures, and excellent agreements are obtained. Meanwhile, the proposed models are applied to investigate the effects of the porous baffles on sloshing response in circular cylindrical, annular cylindrical and conical tanks. The effects of the porous baffle length, porous-effect parameter, installation angle and baffle height on the sloshing force, natural frequency and surface elevation are studied. Additionally, some typical sloshing pressure distributions, velocity potential contours and velocity fields are plotted. The results show that swirls at the tips of the baffles can be observed in many cases, and the top-mounted porous baffle makes more significant suppression effects on sloshing response than that of bottom-mounted porous baffle, while increasing the number of ring porous baffles can achieve better restraint effects on sloshing response. And increasing the baffle length of the horizontal wall-mounted ring porous baffle can significantly decrease the sloshing frequencies, as well as the first non-dimensional natural frequency decreases with decrease in porous-effect parameter of the coaxial porous baffle. In addition, remarkable effects on sloshing can be obtained when reasonable designed by selecting the optimal porous-effect parameter, installation angle and baffle height. And this paper can be a useful guide for the seismic design and analysis of many actual liquid storage tanks (such as the Advanced Passive PWR, large water cooling tower, etc.).     

Investigation of the Effects of Baffles on the Shallow Water Sloshing in A Rectangular Tank Using A 2D Turbulent ISPH Method

Liquid sloshing is a common phenomenon in the liquid tanks transportation. Liquid waves lead to fluctuating forces on the tank wall. Uncontrolled fluctuations lead to large forces and momentums. Baffles can control these fluctuations. A numerical method, which has been widely used to model this phenomenon, is Smoothed Particle Hydrodynamics(SPH). The Lagrangian nature of this method makes it suitable for simulating free surface flows. In the present study, an accurate Incompressible Smoothed Particle Hydrodynamics(ISPH) method is developed and improved using the kernel gradient correction tensors, particle shifting algorithms, k–ε turbulence model, and free surface particle detectors. Comparisons with the experimental data approve the ability of the present algorithm for simulating shallow water sloshing. The main aim of this study is to investigate the effects of the vertical baffle on the damping of liquid sloshing. Results show that baffles number has a major role in sloshing fluctuation damping.     

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.     

Wave impacts on structures with rectangular geometries: Part 2 decks,baffles and seawalls with impermeable or porous surfaces

This paper considers wave impacts on baffles, on baffles or decks adjacent to a vertical wall, and on porous seawalls and/or sea beds. For seawalls and vertical baffles, impacts can occur in steep waves, whilst a deck can be struck from below by a rising wave crest either in open sea or in a tank with standing waves (sloshing). A simple analytical model for the pressure impulse, P, due to a wave of idealized geometry and dynamics is developed and applied to the following geometries with impermeable surfaces:

• •horizontal wave impact onto a vertical wall with a deck at the waterline,
• •vertical wave impact under a deck in the same configuration (equivalent to vertical water impact of a horizontal plate),
• •horizontal wave impact onto a surface-piercing vertical baffle in open sea,
• •as for 3. but with the baffle in front of a wall,
• •as for 4. but with a deck extending from the vertical wall to the baffle,
• •bottom-mounted baffle in front of a wall with impact occurring on the wall.

We also consider cases that complement part 1 of this paper to include the effect on impacts on a seawall with a porous sea bed and/or sea wall with/without a berm. Finally we reconsider case 3) above but with a porous baffle.The method uses eigenfunction expansions in each of the rectangular regions that satisfy some of the impermeable or porous surface conditions, and a simplified free-surface condition. Their unknown coefficients are determined from the impact boundary condition, impermeable or porous boundary conditions and by matching the solutions, in any two neighbouring rectangles, along their common boundary. Although the fluid motion is treated rather crudely, the method yields the pressure impulse throughout the entire region. Impulses, I, and moment impulses, M, on all or parts of the structure are also presented.     

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.     

弹性液舱内液体晃荡实验研究

通过物理模型实验,对水平简谐激励下弹性矩形液舱内液体晃荡问题进行了研究。实验中给出了不同液深和不同激励振幅下弹性液舱内液体最低阶固有频率。采用影像采集与分析系统获得液舱内自由液面的形状和高度,通过压力采集系统得到舱壁上的压力分布。实验中分析了不同液深、不同激励频率下弹性和刚性液舱内自由液面高度和晃荡压力的变化特性。比较了弹性和刚性液舱内不同位置处自由液面高度随外界激励频率的变化规律,以及舱壁上不同位置处的压力峰值随外界激励频率的变化规律。结果表明:非共振情况下,弹性和刚性实验结果及理论值三者的波高和压力较为接近。共振情况下,弹性和刚性波高基本相等,在近舱壁处二者明显大于理论值,而在远离舱壁处二者与理论值存在一定差别;弹性压力较刚性压力略小,但二者与理论值差别较大。     

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.     

New model to determine forces at on-bottom slender pipelines

The present paper proposes a numerical model to determine horizontal and vertical components of the hydrodynamic forces on a slender submarine pipeline lying at the sea bed and exposed to non-linear waves plus a current. The new model is an extension of the Wake II type model, originally proposed for sinusoidal waves (Soedigdo et al., 1999) and for combined sinusoidal waves and currents (Sabag et al., 2000), to the case of periodic or random waves, even with a superimposed current. The Wake II type model takes into account the wake effects on the kinematic field and the time variation of drag and lift hydrodynamic coefficients. The proposed extension is based on an evolutional analysis carried out for each half period of the free stream horizontal velocity at the pipeline. An analytical expression of the wake velocity is developed starting from the Navier–Stokes and the boundary layer equations. The time variation of the drag and lift hydrodynamic coefficients is obtained using a Gaussian integration of the start-up function. A reduced scale laboratory investigation in a large wave flume has been conducted in order to calibrate the empirical parameters involved in the proposed model. Different wave and current conditions have been considered and measurements of free stream horizontal velocities and dynamic pressures on a bottom-mounted pipeline have been conducted. The comparison between experimental and numerical hydrodynamic forces shows the accuracy of the new model in evaluating the time variation of peaks and phase shifts of the horizontal and vertical wave and current induced forces.     

Finite element analysis of two-dimensional nonlinear sloshing problems in random excitations

A two-dimensional nonlinear random sloshing problem is analyzed by the fully nonlinear wave velocity potential theory based on the finite element method. A rectangular container filled with liquid subjected to specified horizontal random oscillations is studied. Both wave elevation and hydrodynamic force are obtained. The spectra of random waves and forces have also been investigated, and the effects of the peak frequencies and spectral width of the specified spectrum used for the generation of the random oscillations are discussed. It is found that the energy mainly concentrates at the natural frequencies of the container and is dominant at the ith order natural frequency when the peak frequency is close to the ith order natural frequency. Some results are compared between the fully nonlinear solutions, the linear solutions and the linear plus second-order solutions.     

A modal approach to second-order analysis of sloshing using boundary element method

This paper aims at developing a modal approach for the non-linear analysis of sloshing in an arbitrary-shape tank under both horizontal and vertical excitations. For this purpose, the perturbation technique is employed and the potential flow is adopted as the liquid sloshing model. The first- and second-order kinematic and dynamic boundary conditions of the liquid-free surface are used along with a boundary element model which is formulated in terms of the velocity potential of the liquid-free surface. The boundary element model is used to determine the natural mode shapes of sloshing and their corresponding frequencies. Using the modal analysis technique, a non-linear model is presented for the calculation of the first- and second-order potential which can be used to obtain a reduced-order model for the sloshing dynamics. The results of the presented model are verified with the analytical solution for the second-order analysis of sloshing in a rectangular tank and very good results were obtained. Also, the second-order sloshing in some other example tanks with complex bed shapes is studied. The second-order resonance conditions of liquid sloshing in the example tanks are investigated and some conclusions are drawn.     

Experimental Study on Sloshing Characteristics in the Elastic Tank Based on Morlet Wavelet Transform

Hydroelastic effect of sloshing is studied through an experimental investigation. Different excitation frequencies are considered with low-fill-depth and large amplitude. Morlet wavelet transform is introduced to analyze the free surface elevations and sloshing pressures. It focuses on variations and distributions of the wavelet energy in elastic tanks. The evolutions of theoretical and experimental wavelet spectra are discussed and the corresponding Fourier spectrums are compared. Afterwards, average values of the wavelet spectra are extracted to do a quantitative study at various points. From the wavelet analysis, sloshing energies are mainly distributed around the external excitation frequency and expanded to high frequencies under violent condition. In resonance, experimental wavelet energy of the elevation in elastic tanks is obviously less than that in the rigid one; for sloshing pressures, the elastic wavelet energy close to the rigid one and conspicuous impulse is observed. It recommends engineers to concern the primary natural frequency and impulsive peak pressures.     

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.     

直墙前矩形月池水波振荡研究分析

月池内流体存在活塞和晃荡两类振荡现象。基于线性势流理论,推导了波浪斜向入射下,直墙前矩形月池辐射和绕射问题的解析解。通过分离变量法和特征函数展开法求解了速度势函数,根据边界条件来确定速度势函数中的未知系数,由速度势函数计算斜向波与矩形月池相互作用的水动力系数和波浪激励力,对它们的变化规律进行了分析讨论,研究了底部开口大小、波浪入射角度对矩形月池水动力特性的影响以及直墙远近对波浪力的影响。结果表明,月池底部开口大小对流体水平作用的影响较小,而对流体垂直作用的影响较大;波浪入射角度的变化对矩形月池横荡和横摇运动时的水动力特性有一定的影响;在一定条件下,直墙的存在会使得月池在水平方向所受到的波浪力比开敞水域中的要大。     

The air–water sloshing problem: Fundamental analysis and parametric studies on excitation and fill levels

The problem of liquid sloshing has gained recent attention with the proliferation of liquefied natural gas (LNG) carriers transporting liquids in partially filled tanks. Impact pressures caused by sloshing depend on the tank fill level, period and amplitude of oscillation of the tank. In this paper, we first present the rudiments of a linear potential theory for sloshing motions in a two-dimensional rectangular tank, due to small amplitude sway motions. Although this topic is fundamental, we clarify inconsistencies in the published literature and texts.Numerical investigations were carried out on the sloshing motions in a two-dimensional tank in the sway excitation. The fluid domain was modeled using a finite volume approximation, and the air–water interface was tracked using a volume-of-fluid (VOF) technique. Computational results for free surface elevation and impact pressure are found to be in good agreement with theory and published data. The fill levels were varied from 10% to 95%, and the excitation time periods were varied from 0.8 to 2.8 s for a constant sway amplitude of 0.25 m (peak–peak) at 1:30 scale. The results of the parametric study are compared with theoretical predictions and suggestions are made on incorporating sloshing effects in standard seakeeping analysis for LNG carriers.     

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.     

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.     

Complete two-dimensional analysis of sea-wave-induced fully non-linear sloshing fluid in a rigid floating tank

A time-independent finite-difference method and a fifth-order Runge–Kutta–Felhberg scheme were used to analyze the dynamic responses of sea-wave-induced fully non-linear sloshing fluid in a floating tank. The interaction effect between the fully non-linear sloshing fluid and the floating tank associated with coupled surge, heave and pitch motions of the tank are analyzed for the first time in the present pilot study. For the analysis of fluid motion in the tank, the coordinate system is moving (translating and rotating) with tank motion. The time-dependent water surface of the sloshing fluid is transformed to a horizontal plane and the flow field is mapped on to a rectangular region. The Euler equations as well as the fully non-linear kinematic free surface condition were used in the analysis of the sloshing fluid. The strip theory for linearized harmonic sea-wave loading was adopted to evaluate the regular encounter wave force. In addition, the dynamic coefficients used in the dynamic equations of tank motion were also derived based on strip theory and a harmonic motion of the tank. The characteristics of free and forced tank motions with and without the sloshing effect are studied. By the damping effect, the response of free oscillation will damp out and that of forced oscillation will approach a steady state. Without sea-wave action, the contribution of the sloshing load would enlarge the angular response of tank motion as well as the rise of free surface and the sloshing effect will delay the damping effect on angular displacement. On the contrary, under sea-wave action, the sloshing effect will decrease the dynamic response of tank motion and rise of free surface. The interaction, sloshing and coupling effects are found to be significant and should be considered in the analysis and design of floating tanks.