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.     

Approximation of near-resonant wave motion using a damped harmonic oscillator model

In this work, transient resonant motions excited in linearised scattering interactions are approximated using a simple damped harmonic oscillator model. The scattering interactions considered involve the diffraction of an incident wave-train with a dominant regular time-harmonic component by structures which enclose a portion of the free surface. Provided a single resonant mode only is excited significantly, the fluid oscillation in the vicinity of the structure will primarily be composed of the resonant mode and incident wave mode contributions. The forced damped harmonic oscillator equation is used to predict the fluid motion and, in particular, the elevation of the internal free surface. The predictions are compared to the results from a numerical time-domain solver based on the linearised water-wave equations. It is shown that, given a good estimate of the location of the resonance in the complex frequency plane and a priori knowledge of the form of the incident wave, the model can successfully predict the time-dependent behaviour of the motion. Both two- and three-dimensional scattering problems are considered with a variety of scatterers in each case.     

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.     

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.     

Roll performance of a small fishing vessel with live fish tank

The concept of live fish tanks in trawlers is to use the catch in a better condition and to reduce marine pollution. It also reduces the infrastructure meant to freeze the catch to preserve it for longer period. But the presence of additional free surface in the vessel challenges the stability of the vessel. This is besides the sloshing effect due to the moving liquid mass in the tank. Roll motions are initiated due to various factors related to the hull characteristics of the vessel, loading and operating conditions and its interaction with the environment. Location of fish tank, its orientation, arrangement of baffles inside the tank to reduce the free surface affects and careful design of tank opening are to be given priority during the design, manufacturing and tank testing. The results obtained from tank test of model are compared with that of analytical method. The non-linear roll performance become further complicated due to the free surface and sloshing effects of the mass in the live fish tank. Wave makers are used for generating waves under laboratory conditions compatible with the scaled down model of the trawler model. The tests are conducted in the towing tank of Pusan National University.     

Nonlinear modeling of liquid sloshing in a moving rectangular tank

A nonlinear liquid sloshing inside a partially filled rectangular tank has been investigated. The fluid is assumed to be homogeneous, isotropic, viscous, Newtonian and exhibit only limited compressibility. The tank is forced to move harmonically along a vertical curve with rolling motion to simulate the actual tank excitation. The volume of fluid technique is used to track the free surface. The model solves the complete Navier–Stokes equations in primitive variables by use of the finite difference approximations. At each time step, a donor–acceptor method is used to transport the volume of fluid function and hence the locations of the free surface. In order to assess the accuracy of the method used, computations are verified through convergence tests and compared with the theoretical solutions and experimental results.     

Analysis of second-order resonance in wave interactions with floating bodies through a finite-element method

A time-domain method is employed to analyse the resonant oscillations of the liquid confined within the two floating bodies. The velocity potentials at each time step are obtained through a finite-element method (FEM) with quadratic shape functions. The matrix equation of the FEM is solved through an iteration. The radiation condition is satisfied through a combination of the damping zone method and the Sommerfeld–Orlanski equation. A detailed analysis is made for two rectangular floating cylinders undergoing forced oscillation. The first-order potential reveals the resonant behaviour of the wave motion at certain frequencies ω_i, which is similar to sloshing in a tank. More interestingly, the second-order theory further reveals that when the oscillation frequency is at ω_i/2 or half of the resonant frequency, no first-order resonance is observed as expected, but the second-order resonant motion becomes evident, which does not seem to have been extensively investigated so far. Detailed results for two rectangular cylinders are provided to show some insights into the resonant effect due to the interaction between the bodies. The first- and second-order resonant phenomena have been observed and the result has shown that the second-order components have significant influence on the wave and force in some cases, especially at the second-order resonance.     

Numerical simulation of non-linear regular and focused waves in an infinite water-depth

Inviscid three-dimensional free surface wave motions are simulated using a novel quadratic higher order boundary element model (HOBEM) based on potential theory for irrotational, incompressible fluid flow in an infinite water-depth. The free surface boundary conditions are fully non-linear. Based on the use of images, a channel Green function is developed and applied to the present model so that two lateral surfaces of an infinite-depth wave tank can be excluded from the calculation domain. In order to generate incident waves and dissipate outgoing waves, a non-reflective wave generator, composed of a series of vertically aligned point sources in the computational domain, is used in conjunction with upstream and downstream damping layers. Numerical experiments are carried out, with linear and fully non-linear, regular and focused waves. It can be seen from the results that the present approach is effective in generating a specified wave profile in an infinite water-depth without reflection at the open boundaries, and fully non-linear numerical simulations compare well with theoretical solutions. The present numerical technique is aimed at efficient modelling of the non-linear wave interactions with ocean structures in deep water.     

Hydroelastic analysis of floating structures with liquid tanks and comparison with experimental tests

In this study, we develop a numerical method for a 3D linear hydroelastic analysis of floating structures with liquid tanks subjected to surface regular water waves and compare the numerical results with experimental tests. Considering direct couplings among structural motion, sloshing, and water waves, a mathematical formulation and a numerical method extended from a recent work [1] are developed. The finite element method is employed for the floating structure and internal fluid in tanks, and the boundary element method is used for the external fluid. The resulting formulation completely incorporates all the interaction terms including hydrostatic stiffness and the irregular frequency effect is removed by introducing the extended boundary integral equations. Through various numerical tests, we verify the proposed numerical method. We also performed 3D hydroelastic experimental tests of a floating production unit (FPU) model in an ocean basin. The measured dynamic motions are compared with the numerical results obtained using the proposed method.     

Water entry of a finite width wedge near a floating body

The problem of a two-dimensional finite-width wedge entering water near a freely floating body is considered through the velocity potential theory for the incompressible liquid with the fully nonlinear boundary conditions on the free surface. The problem is solved by using the boundary element method in the time domain. The numerical process is divided into two phases based on whether the interaction between the wedge and floating body is significant. In the first phase, when the single wedge enters water at initial stage, only a small part near its tip is in the fluid, the problem is studied in a stretched coordinate system and the presence of the floating body has no major effect. In the second phase, the disturbance by water entry of the wedge has reached the floating body, and both are considered together in the physical system. The auxiliary function method is adopted to decouple the nonlinear mutual dependence between the motions of the wedge and floating body, both in three degrees of freedom, and the fluid flow, as well as the interaction effects between them. Case studies are undertaken for a wedge entering water in forced or free fall motion, vertically or obliquely. Results are provided for the accelerations, velocities, pressure distribution and free surface deformation, and the interaction effects are discussed.     

A finite element model of interaction between viscous free surface waves and submerged cylinders

This paper describes the simulation of the flow of a viscous incompressible Newtonian liquid with a free surface. The Navier–Stokes equations are formulated using a streamline upwind Petrov–Galerkin scheme, and solved on a Q-tree-based finite element mesh that adapts to the moving free surface of the liquid. Special attention is given to fitting the mesh correctly to the free surface and solid wall boundaries. Fully non-linear free surface boundary conditions are implemented. Test cases include sloshing free surface motions in a rectangular tank and progressive waves over submerged cylinders.     

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.     

Two phase modal analysis of nonlinear sloshing in a rectangular container

Sloshing, or liquid free surface oscillation, in containers has many important applications in a variety of engineering fields. The modal method can be used to solve linear sloshing problems and is the most efficient reduced order method that has been used during the previous decade. In the present article, the modal method is used to solve a nonlinear sloshing problem. The method is based on a potential flow solution that implements a two-phase analysis on sloshing in a rectangular container. According to this method, the solution to the mass conservation equation, with a nonpenetration condition at the tank walls, results in velocity potential expansion; this is similar to the mode shapes used in modal method. The kinematic and dynamic boundary conditions create a set of two-space-dimensional differential equations with respect to time. The numerical solution of this set of differential equations, in the time domain, predicts the time response of interfacial oscillations. Modal method solutions for the time response of container sloshing due to lateral harmonic oscillations show a good agreement with experimental and numerical results reported in the literature.     

On damping models in free and forced rolling motion

In this paper two different models for the damping moment to introduce in the rolling equation of the ship are proposed. They contain two terms, respectively linear-quadratic and linear-cubic in the angular velocity, and furthermore they foresee a non-linear term representing the dependence of the damping from the heeling angle. These models constitute a generalization of all the models up to now used in the naval literature.With the Bogoliubov-Krilov asymptotic method approximate relations, describing the decay curve of the free oscillations and the maximum roll amplitude in synchronism condition, are obtained. The analysis shows that the choice of the more realistic damping model cannot be based on the simple verification of a good fitting of the free oscillation decay curves. It is necessary to examine also the behaviour of the forced oscillations in synchronism.Finally, a plan of experiments which allows the determination of separate values for the different non-linear damping coefficients is proposed.     

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

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

Coupled analysis of nonlinear sloshing and ship motions

A coupled numerical model considering nonlinear sloshing flows and the linear ship motions has been developed based on a boundary element method. Hydrodynamic performances of a tank containing internal fluid under regular wave excitations in sway are investigated by the present time-domain simulation model and comparative model tests. The numerical model features well the hydrodynamic performance of a tank and its internal sloshing flows obtained from the experiments. In particular, the numerical simulations of the strong nonlinear sloshing flows at the natural frequency have been validated. The influence of the excitation wave height and wave frequency on ship motions and internal sloshing has been investigated. The magnitude of the internal sloshing increases nonlinearly as the wave excitation increases. It is observed that the asymmetry of the internal sloshing relative to still water surface becomes more pronounced at higher wave excitation. The internal sloshing-induced wave elevation is found to be amplitude-modulated. The frequency of the amplitude modulation envelope is determined by the difference between the incident wave frequency and the natural frequency of the internal sloshing. Furthermore, the coupling mechanism between ship motions and internal sloshing is discussed.     

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.     

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.     

Fluid-induced inertia and damping in vibrating offshore structures

The problem addressed in this paper is that of the fluid-structure interaction in offshore structures for which effects due to fluid viscosity are negligible. A method for evaluating the dynamic response of structures is described in which the generalized coordinates for the analysis are associated with the ‘dry’ modes of the structure. All hydrodynamic actions, including those induced by motions of the structure, are retained on the right-hand side of the equations of motion. This formulation requires the determination of generalized added mass and damping matrices associated with motions of the structure in its first few ‘dry’ modes. By these means, free surface and three dimensional fluid flow effects are retained in the analysis. Examples of the resulting frequency dependent matrices, computed using the boundary integral method, are presented for some common structural forms, and the validity of this approach is demonstrated through some experiments on a deforming column structure. Finally, results for the wave-excited dynamic response of a typical offshore structure obtained using (i) conventional strip theory ‘wet’ modes, and (ii) ‘dry’ modes, are compared.