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.     

Comparative Study of Different SPH Schemes on Simulating Violent Water Wave Impact Flows

Free surface flows are of significant interest in Computational Fluid Dynamics（CFD）. However, violent water wave impact simulation especially when free surface breaks or impacts on solid wall can be a big challenge for many CFD techniques. Smoothed Particle Hydrodynamics（SPH） has been reported as a robust and reliable method for simulating violent free surface flows. Weakly compressible SPH（WCSPH） uses an equation of state with a large sound speed, and the results of the WCSPH can induce a noisy pressure field and spurious oscillation of pressure in time history for wave impact problem simulation. As a remedy, the truly incompressible SPH（ISPH） technique was introduced, which uses a pressure Poisson equation to calculate the pressure. Although the pressure distribution in the whole field obtained by ISPH is smooth, the stability of the techniques is still an open discussion. In this paper, a new free surface identification scheme and solid boundary handling method are introduced to improve the accuracy of ISPH. This modified ISPH is used to study dam breaking flow and violent tank sloshing flows. On the comparative study of WCSPH and ISPH, the accuracy and efficiency are assessed and the results are compared with the experimental data.     

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

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

基于Level-set法和通度概念液体晃荡特性研究

液体晃荡现象在诸多工程领域具有深刻的研究意义,并引起了广泛的关注。基于Level-set方法,借助流场通度的概念,模拟了流场内具有障碍物的液体晃荡现象。选取不同形式的防晃结构分析研究对晃荡的抑制效果,得到几点关于防晃结构的设计与选择的重要结论。计算结果表明,通度系数法与Level-set方法的结合使用,能够有效地处理流场内带有障碍物的液体晃荡问题,该方法具有一定的可行性和应用前景。     

Comparison of laminar model,RANS, LES and VLES for simulation of liquid sloshing

Liquid sloshing in storage tank is a fundamental problem of great engineering importance. Sloshing motion can be laminar or turbulent. However, the necessity for inclusion of turbulence in CFD simulation of sloshing flows has not yet been established. In this paper, three roll–induced sloshing cases are studied to assess the merits and shortcomings of the laminar model and three most–commonly used turbulence models (RANS k–ε, LES and Very LES). To overcome the deficiencies in the RANS and LES, the new Very LES (VLES) model, which combines the RANS k–ε and LES, is developed in this paper. The free surface profiles are reconstructed by a coupled Level–Set and Volume–of–Fluid (CLSVOF) method. To the authors’ knowledge, the comprehensive and systematical assessment of the effect of turbulence on sloshing simulation has not been reported in the literature. The numerical results are evaluated using experimental measurements from Delorme and Souto−Iglesias. The present study indicates that the inclusion of an appropriate turbulence model has a profound influence on the simulations of violent and non–violent sloshing flows. The VLES and LES models can provide accurate predictions of free surface profiles and impact pressures, whereas the laminar flow assumption and the RANS model cannot adequately capture the energy dissipation in the sloshing simulation and lead to the inaccurate flow predictions.     

Incompressible SPH flow model for wave interactions with porous media

The paper presents an Incompressible Smoothed Particle Hydrodynamics (ISPH) method to simulate wave interactions with a porous medium. The SPH method is a mesh free particle modeling approach that is capable of tracking the large deformation of free surfaces in an easy and accurate manner. The ISPH method employs a strict incompressible hydrodynamic formulation to solve the fluid pressure and the numerical solution is obtained by using a two-step semi-implicit scheme. The ISPH flow model solves the unsteady 2D Navier–Stokes (NS) equations for the flows outside the porous media and the NS type model equations for the flows inside the porous media. The presence of porous media is considered by including additional friction forces into the equations. The developed ISPH model is first validated by the solitary and regular waves damping over a porous bed and the solitary wave interacting with a submerged porous breakwater. The convergence of the method and the sensitivity of relevant model parameters are discussed. Then the model is applied to the breaking wave interacting with a breakwater covered with a layer of porous materials. The computational results demonstrate that the ISPH flow model could provide a promising simulation tool in coastal hydrodynamic applications.     

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.     

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.     

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.     

Numerical Analysis of the Behavior of A New Aeronautical Alloy (Ti555-03) Under the Effect of A High-Speed Water Jet

In this paper, we present a numerical simulation of a water jet impacting a new aeronautical material Ti555-03 plate.The Computational Fluid Dynamics(CFD) behavior of the jet is investigated using a FV(Finite Volume) method.The Fluid–Structure Interaction(FSI) is studied using a coupled SPH(Smoothed Particle Hydrodynamics)-FE(Finite Element) method. The jets hit the metal sheet with an initial velocity 500 m/s. Two configurations which differ from each other by the position(angle of inclination) of the plate relatively to the axis of revolution of the jet inlet are investigated in this study. The objective of this study is to predict the impact of the fluid produced at high pressure and high speed especially at the first moment of impact. Numerical simulations are carried out under ABAQUS. We have shown in this study that the inclination of the titanium alloy plate by 45° stimulates the phenomenon of recirculation of water. This affects the velocity profile, turbulence and boundary layers in the impact zone. The stagnation zone and the phenomenon of water recirculation are strongly influenced by the slope of the plate which gives a pressure gradient and displacement very important between the two configurations. Fluctuations of physical variables(displacement and pressure) prove the need for a noise and vibratory study. These predictions will subsequently be used for the modeling of the problem of an orthogonal cut in a high-speed machining process assisted by high-pressure water jet.     

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.     

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.     

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.     

Numerical simulation of sloshing with large deforming free surface by MPS-LES method

Moving particle semi-implicit (MPS) method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation (LES) by coupling with a sub-particle-scale (SPS) turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is "neighbor particle". In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no "neighbor particles" in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.     

Green water overtopping analyzed with a SPH model

Wave overtopping on the decks of offshore platforms and ships can cause severe damage due to the high forces generated by the water. This phenomenon is analyzed within the framework of the Smoothed Particle Hydrodynamics (SPH) method. The presence of a fixed horizontal deck above the mean water level modifies strongly the wave kinematics. In particular, the flow in the wave crest is split into two, showing a different behavior above and below the deck. Numerical results generated by the SPH method are compared to laboratory experiments. The formation of a jet in the rear of the deck after overtopping is observed under extreme conditions.     

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.     

容器固有频率对液体晃荡的影响

针对矩形容器内液体晃荡问题,采用了时域高阶边界元方法建立自由水面满足完全非线性边界条件的数学模型。求解中采用混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过将计算得到的波面结果与实验数据、解析解和已发表结果对比,吻合良好,验证了本方法的准确性。进而采用谱分析方法分析了波面时间历程,得到容器各阶固有频率对液体晃荡的影响。研究发现,基频对液体晃荡的影响最大,且非线性越强,更高阶容器固有频率的影响越大。     

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.     

Simulation of wave overtopping by an incompressible SPH model

The paper presents an incompressible Smoothed Particle Hydrodynamics (SPH) model to investigate the wave overtopping of coastal structures. The SPH method is a grid-less Lagrangian approach which is capable of tracking the large deformations of the free surface with good accuracy. The incompressible algorithm of the model is implemented by enforcing the constant particle density in the pressure projection. The SPH model is employed to reproduce a transient wave overtopping over a fixed horizontal deck and the regular/irregular waves overtopping of a sloping seawall. The computations are validated against the experimental and numerical data and a good agreement is observed. The SPH modelling is shown to provide a promising tool to predict the overtopping characteristics of different waves. The present model is expected to be of practical purpose if further improvement in the spatial resolution and CPU time can be adequately made.     

Energy dissipation in sloshing waves in a rolling rectangular tank — III. Results and applications

The hydrodynamic aspects of the motion of a viscous fluid having a free surface in a rolling tank have been investigated. In a sequence of three papers, an analytical technique together with a numerical solution method will be presented to describe the sloshing phenomenon accurately and efficiently. This first paper introduces a linear theory of viscous liquid sloshing and formulates a boundary value problem subject to appropriate conditions. Viscosity is included in the problem formulation and its effects are properly accounted for. The second paper will describe a solution of the problem in function space by the truncation of infinite series. Boundary conditions are satisfied through the use of Fourier series expansions. However, the no-slip condition at the side walls can also be treated in a least-squares sense.Among the results that will be reported in the third paper are the effects of viscosity on liquid sloshing phenomenon and the dependence of viscous dissipation on the Reynolds and Froude numbers. Furthermore, the influence of the tank aspect ratio on viscous dissipation has been explored. These results demonstrate some unknown features of the functional relationships that exist between the dissipated energy and the Reynolds and Froude numbers. Similarly, the dependence of the dissipated energy on the aspect ratio has been analytically studied. The results obtained agree with the physical laws for the range of parameters investigated.