Numerical analysis of submarine landslides using a smoothed particle hydrodynamics depth integral model

Submarine landslides can cause severe damage to marine engineering structures. Their sliding velocity and runout distance are two major parameters for quantifying and analyzing the risk of submarine landslides.Currently, commercial calculation programs such as BING have limitations in simulating underwater soil movements. All of these processes can be consistently simulated through a smoothed particle hydrodynamics（SPH） depth integrated model. The basis of the model is a control equation that was developed to take into account the effects of soil consolidation and erosion. In this work, the frictional rheological mode has been used to perform a simulation study of submarine landslides. Time-history curves of the sliding body＇s velocity, height,and length under various conditions of water depth, slope gradient, contact friction coefficient, and erosion rate are compared; the maximum sliding distance and velocity are calculated; and patterns of variation are discussed.The findings of this study can provide a reference for disaster warnings and pipeline route selection.

Numerical Simulation of Water Mitigation Effects on Shock Wave with SPH Method

The water mitigation effect on the propagation of shock wave was investigated numerically. The traditional smoothed particle hydrodynamics (SPH) method was modified based on Riemann solution. The comparison of numerical results with the analytical solution indicated that the modified SPH method has more advantages than the traditional SPH method. Using the modified SPH algorithm, a series of one-dimensional planar wave propagation problems were investigated, focusing on the influence of the air-gap between the high-pressure air and water and the thickness of water. The numerical results showed that water mitigation effect is significant. Up to 60% shock wave pressure reduction could be achieved with the existence of water, and the shape of shock wave was also changed greatly. It is seemly that the small air-gap between the high-pressure air and water has more influence on water mitigation effect.

Multi-resolution technique integrated with smoothed particle element method (SPEM) for modeling fluid-structure interaction problems with free surfaces

Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element method(SPEM) to simulate FSIs. In this method, both the fluid and solid regions are initially modeled using a smoothed finite element method(S-FEM) in a Lagrangian frame, whereas the fluid regions undergoing large deformations are adaptively converted into particles and modeled with an improved smoothed particle hydrodynamics(SPH) method. This approach greatly improves computational accuracy and efficiency because of the advantages of the S-FEM in efficiently treating solid/fluid regions showing small deformations and the SPH method in effectively modeling moving interfaces. In this work, we further enhance the efficiency of the SPEM while effectively capturing local fluid information by introducing a multi-resolution technique to the SPEM and developing an effective approach to treat multi-resolution element-particle interfaces. Various numerical examples demonstrate that the multiresolution SPEM can significantly reduce the computational cost relative to the original version with a constant resolution.Moreover, the novel approach is effective in modeling various incompressible flow problems involving FSIs.

Numerical modeling of sediment transport based on unsteady and steady flows by incompressible smoothed particle hydrodynamics method

The purpose of the present paper is to introduce a simple two-part multi-phase model for the sediment transport problems based on the incompressible smoothed particle hydrodynamics（ISPH） method. The proposed model simulates the movement of sediment particles in two parts. The sediment particles are classified into three categories, including the motionless particles, moving particles behave like a rigid body, and moving particles with a pseudo fluid behavior. The criterion for the classification of sediment particles is the Bingham rheological model. Verification of the present model is performed by simulation of the dam break waves on movable beds with different conditions and the bed scouring under steady flow condition. Comparison of the present model results, the experimental data and available numerical results show that it has good ability to simulate flow pattern and sediment transport.

Numerical simulation of submarine landslide tsunamis using particle based methods

This paper presents the simulation of tsunamis due to rigid and deformable landslides with consideration of submerged conditions by using particle methods. The smoothed particle hydrodynamics（SPH）, as a particle based method, is for solving problems of fast moving boundaries in the field of continuum mechanics. Other particle based methods, like the discrete element method（DEM）, are suitable for modeling the displacement and the collision related to the rigid landslides. In the present work, we use the SPH and the DEM to simulate tsunamis generated by rigid and deformable landslides with consideration of submerged conditions. The viscous free-surface flows are solved by a weakly compressible SPH and the displacement and the rotation of the rigid body slides are calculated using a multi-sphere DEM allowing for modeling solids of arbitrarily complex shapes. The fluid-solid interactions are simulated by coupling the SPH and the DEM. A rheology model combining the Papanastasiou and the Herschel-Bulkley models is applied to represent the viscoplastic behavior of the non-Newtonian flow in the submarine deformable landslide cases. Submarine landslide tsunamis due to rigid and deformable landslides are both simulated as typical landslide cases in this investigation. Our simulated results and the previous experimental results in the literatures are in good agreement, which shows that the proposed particle based methods are capable of modeling the submarine landslide tsunamis.

Numerical Simulation of Interaction of a Current with a Circular Cylinder near a Rigid Bed

The numerical modeling of a 2D flow around a horizontal cylinder near a rigid bed with the gap ratios G/D = 0, 0.2, 0.3 at Reynolds numbers 840, 4500, 9000, and 9500 is investigated by using weakly compressible smoothed particle hydrodynamics. The velocity field and the separation angles from the present simulations are compared with those obtained from the experimental measurements and are in a good agreement. The results show that the maximum value of shear stress on the bed increases as the cylinder closes the bed and suddenly decreases when the cylinder contact the wall.

Numerical modeling of dam-break flood through intricate city layouts including underground spaces using GPU-based SPH method

This paper applies the meshfree Smoothed Particle Hydrodynamics （SPH） method with Graphical Processing Unit （GPU） parallel computing technique to investigate the highly complex 3-D dam-break flow in urban areas including underground spaces. Taking the advantage of GPUs parallel computing techniques, simulations involving more than 107 particles can be achieved. We use a virtual geometric plane boundary to handle the outermost solid wall in order to save considerable video card memory for the GPU computing. To evaluate the accuracy of the new GPU-based SPH model, qualitative and quantitative comparison to a real flooding experiment is performed and the results of a numerical model based on Shallow Water Equations （SWEs） is given with good accuracy. With the new GPU-based SPH model, the effects of the building layouts and underground spaces on the propagation of dambreak flood through an intricate city layout are examined.

TRACKING METHODS FOR FREE SURFACE AND SIMULATION OF A LIQUID DROPLET IMPACTING ON A SOLID SURFACE BASED ON SPH

With some popular tracking methods for free surface, simulations of several typical examples are carried out under various flow field conditions. The results show that the Smoothed Particle Hydrodynamics （SPH） method is very suitable in simulating the flow problems with a free surface. A viscous liquid droplet with an initial velocity impacting on a solid surface is simulated based on the SPH method, and the surface tension is considered by searching the free surface particles, the initial impact effect is considered by using the artificial viscosity method, boundary virtual particles and image virtual particles are introduced to deal with the boundary problem, and the boundary defect can be identified quite well. The comparisons of simulated results and experimental photographs show that the SPH method can not only exactly simulate the spreading process and the rebound process of a liquid droplet impacting on a solid surface but also accurately track the free surface particles, simulate the free-surface flow and determine the shape of the free surface due to its particle nature.

An Energy Stable SPH Method for Incompressible Fluid Flow

In this paper,a novel unconditionally energy stable Smoothed Particle Hydrodynamics(SPH)method is proposed and implemented for incompressible fluid flows.In this method,we apply operator splitting to break the momentum equation into equations involving the non-pressure term and pressure term separately.The idea behind the splitting is to simplify the calculation while still maintaining energy stability,and the resulted algorithm,a type of improved pressure correction scheme,is both efficient and energy stable.We show in detail that energy stability is preserved at each full-time step,ensuring unconditionally numerical stability.Numerical examples are presented and compared to the analytical solutions,suggesting that the proposed method has better accuracy and stability.Moreover,we observe that if we are interested in steady-state solutions only,our method has good performance as it can achieve the steady-state solutions rapidly and accurately.

A numerical model for bird strike on sidewall structure of an aircraft nose

In order to examine the potential of using the coupled smooth particles hydrodynamic （SPH） and finite element （FE） method to predict the dynamic responses of aircraft structures in bird strike events, bird-strike tests on the sidewall structure of an aircraft nose are carried out and numerically simulated. The bird is modeled with SPH and described by the Murnaghan equation of state, while the structure is modeled with finite elements. A coupled SPH-FE method is developed to simulate the bird-strike tests and a numerical model is established using a commercial software PAM-CRASH. The bird model shows no signs of instability and correctly modeled the break-up of the bird into particles. Finally the dynamic response such as strains in the skin is simulated and compared with test results, and the simulated deformation and fracture process of the sidewall structure is compared with images recorded by a high speed camera. Good agreement between the simulation results and test data indicates that the coupled SPH-FE method can provide a very powerful tool in predicting the dynamic responses of aircraft structures in events of bird strike.

FLUID-STRUCTURE INTERACTION OF HYDRODYNAMIC DAMPER DURING THE RUSH INTO THE WATER CHANNEL

The hydrodynamic damper is a device to decrease the motion of armament carrier by use of the water resistance. When hydrodynamic damper rushes into the water channel with high velocity, it is a complicated flow phenomenon with fluid-structure interaction, free surface and moving interface. Numerical simulation using the Smoothed Particle Hydrodynamics （SPH） method coupled with the Finite Element （FE） method was successfully conducted to predict the dynamic characteristics of hydrodynamic damper. The water resistance, the pressure in the interface and the stress of structure were investigated, and the relationship among the peak of water resistance, initial velocity and actual draught was also discussed. The empirical formula was put forward to predict the water resistance. And it is found that the resistance coefficient is commonly in the range of 0.3 ≤ C ≤ 0.5, when the initial velocity is larger than 50 m/s. It can be seen that the SPH method coupled with the FE method has many obvious advantages over other numerical methods for this complicated flow problem with fluid-structure interaction.

An ISPH model for flow-like landslides and interaction with structures

An incompressible smoothed particle hydrodynamics（ISPH） model has been developed to investigate the flow-like landslide phenomena The landslide mass is idealized as rigid and perfectly-plastic material with a constant density. Unlike the widely-used explicit smoothed particle hydrodynamics（SPH） models for landslides, the Chorin＇s projection method is used herein to implicitly solve the normal stress via a pressure Poisson equation, leading to a realistic distribution of instantaneous stress fields free of spurious fluctuations. The capability of the model is demonstrated through three case studies, including the idealized granular flow, landslide interaction with a rigid barrier, and a historical cut-slope landslide in Hong Kong.