Smoothed-Particle Hydrodynamics Simulation of Ship Motion and Tank Sloshing under the Effect of Regular Waves

Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to ship motion and liquid tank sloshing has been simulated by the Smoothed-Particle Hydrodynamics(SPH)method.Firstly,the sloshing flow in a rectangular tank was simulated and the related loads were analyzed to verify and validate the accuracy of the present SPH solver.Then,a three-dimensional simplified LNG carrier model,including two prismatic liquid tanks and a wave tank,was introduced.Different conditions were examined corresponding to different wave lengths,wave heights,wave heading angles,and tank loading rates.Finally,the effects of liquid tank loading rate on LNG ship motions and sloshing loading were analyzed,thereby showing that the SPH method can effectively provide useful indications for the design of liquid cargo ships.

A modified smoothed particle hydrodynamics method considering residual stress for simulating failure and its application in layered rock mass

Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.

An Innovative Coupled Common-Node Discrete Element Method-Smoothed Particle Hydrodynamics Model Developed with LS-DYNA and Its Applications

In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.

Application of Smoothed Particle Hydrodynamics(SPH)for the Simulation of Flow-Like Landslides on 3D Terrains

Flow-type landslide is one type of landslide that generally exhibits characteristics of high flow velocities,long jump distances,and poor predictability.Simulation of its propagation process can provide solutions for risk assessment and mitigation design.The smoothed particle hydrodynamics(SPH)method has been successfully applied to the simulation of two-dimensional(2D)and three-dimensional(3D)flow-like landslides.However,the influence of boundary resistance on the whole process of landslide failure is rarely discussed.In this study,a boundary condition considering friction is proposed and integrated into the SPH method,and its accuracy is verified.Moreover,the Navier-Stokes equation combined with the non-Newtonian fluid rheologymodel was utilized to solve the dynamic behavior of the flow-like landslide.To verify its performance,the Shuicheng landslide event,which occurred in Guizhou,China,was taken as a case study.In the 2D simulation,a sensitivity analysis was conducted,and the results showed that the shearing strength parameters have more influence on the computation accuracy than the coefficient of viscosity.Afterwards,the dynamic characteristics of the landslide,such as the velocity and the impact area,were analyzed in the 3D simulation.The simulation results are in good agreement with the field investigations.The simulation results demonstrate that the SPH method performs well in reproducing the landslide process,and facilitates the analysis of landslide characteristics as well as the affected areas,which provides a scientific basis for conducting the risk assessment and disaster mitigation design.

Improved Model for Soil as a Two-Phase Mixture Based on Smoothed Particle Hydrodynamics (SPH)

It is desired to resolve soil contamination with reduced costs. “Insoluble treatment” is a soil improvement method for heavy metal containing soil, which uses soil mixers to mix soil and soil improvement liquid agents. To reduce the costs of this method, soil mixers have to be optimized. However, it is not achieved due to the lack of theoretical knowledge on mixing solid with liquid. Therefore, a numerical model to simulate the dynamic behavior of solid and liquid is on the development in this study using Smoothed Particle Hydrodynamics (SPH) method. To validate the numerical model, several experiments were carried out and numerically reproduced. The comparisons of the results showed that the numerical model replicated a liquid flow with an error rate of 2.1% and a seepage flow with an error rate up to 26.1%. Especially, the water distribution in the soil pores was highly improved with absolute gaps in volumetric water content up to 4.4% in the porosity range of 10% - 90%. For the water absorption into dry sand, the simulation result became more realistic by concerning soil suction.

基于SPH算法的热仿真系统开发

随着电子设备日益微型化和集成化,热仿真已成为其设计中的关键因素。电子封装模块的热仿真通常使用传统的有限元法FEM(finite element method),存在计算效率和精度之间的矛盾,在处理大变形问题和网格畸变方面也容易造成计算不收敛,从而导致结果错误。针对该问题,提出一种基于光滑粒子动力学SPH(smoothed particle hydrodynamics)算法的电子封装模块热仿真系统。该算法基于无网格拉格朗日数值方法,通过将热仿真对象离散为1组粒子的方式求解热传导方程,从而准确地预测电子封装模块的传热与散热,无需生成并处理大量的微小网格,不用担心网格失真等问题。SPH相对于FEM,仿真精度误差保持在1%~2%,仿真效率可提升近30倍,适合用于复杂和动态系统的模拟仿真。

基于Smoothed Particle Hydrodynamics方法的实时流体模拟

提出基于物理的、实时的技术模拟水体,根据Smoothed Particle Hydrodynamics方法求解流体动力学的Navier-Stokes方程,并运用Marching Cubes体绘制算法重建流体表面。实验表明基于粒子的方法能模拟流体所特有的飞溅、漩涡等现象,Marching Cubes在表面重建的高效性使模拟达到实时、交互的应用。

A novel approach to determine residual stress field during FSW of AZ91 Mg alloy using combined smoothed particle hydrodynamics/neuro-fuzzy computations and ultrasonic testing

The faults in welding design and process every so often yield defective parts during friction stir welding(FSW).The development of numerical approaches including the finite element method(FEM)provides a way to draw a process paradigm before any physical implementation.It is not practical to simulate all possible designs to identify the optimal FSW practice due to the inefficiency associated with concurrent modeling of material flow and heat dissipation throughout the FSW.This study intends to develop a computational workflow based on the mesh-free FEM framework named smoothed particle hydrodynamics(SPH)which was integrated with adaptive neuro-fiizzy inference system(ANFIS)to evaluate the residual stress in the FSW process.An integrated SPH and ANFIS methodology was established and the well-trained ANIS was then used to predict how the FSW process depends on its parameters.To verify the SPH calculation,an itemized FSW case was performed on AZ91 Mg alloy and the induced residual stress was measured by ultrasonic testing.The suggested methodology can efficiently predict the residual stress distribution throughout friction stir welding of AZ91 alloy.

Numerical Investigation of Penetration in Ceramic/Aluminum Targets Using Smoothed Particle Hydrodynamics Method and Presenting a Modified Analytical Model

Radius of ceramic cone can largely contribute into final solution of analytic models of penetration into ceramic/metal targets.In the present research,a modified model based on radius of ceramic cone was presented for ceramic/aluminum targets.In order to investigate and evaluate accuracy of the presented analytic model,obtained results were compared against the results of the Florence’s analytic model and also against numerical modeling results.The phenomenon of impact onto ceramic/aluminum composites were modeled using smoothed particle hydrodynamics(SPH)implemented utilizing ABAQUS Software.Results indicated that,with increasing initial velocity and ceramic thickness and decreasing support layer thickness,the radius of ceramic cone decreases;this ends up increasing residual velocity of the projectile and penetration time and extending the area across which the pressure is distributed.These findings indicate enhanced levels of target energy absorption and the required energy for bending and tensioning the target.As such,it can be observed that,at the same thickness and areal density,the ceramic target has its efficiency enhanced with increasing ceramic thickness and decreasing the support layer thickness.Finally,the results revealed that the associated data with SPH confirm the modified analytic model at higher accuracy than the Florence’s analytic model.

Quasi-static simulation of droplet morphologies using a smoothed particle hydrodynamics multiphase model

Numerical simulation of the morphology of a droplet deposited on a solid surface requires an efficient description of the three-phase contact line. In this study, a simple method of implementing the contact angle is proposed, combined with a robust smoothed particle hydrodynamics multiphase algorithm (Zhang 2015). The first step of the method is the creation of the virtual liquid-gas interface across the solid surface by means of dummy particles, thus the calculated surface tension near the triple point serves to automatically modulate the dynarnic contact line towards the equilibrium state. We simulate the evolution process of initially square liquid lumps on fiat and curved surfaces. The predictions of droplet profiles are in good agreement with the analytical solutions provided that the macroscopic contact angle is accurately implemented. Compared to the normal correction method, the present method is straightforward without the need to manually alter the normal vectors. This study presents a robust algorithm capable of capturing the physics of the static welling. It may hold great potentials in bio-inspired superhydrophobic surfaces, oil displacement, microfluidics, ore floatation, etc.

Parametric study on smoothed particle hydrodynamics for accurate determination of drag coefficient for a circular cylinder

Simulations of two-dimensional(2D) flow past a circular cylinder with the smoothed particle hydrodynamics(SPH) method were conducted in order to accurately determine the drag coefficient. The fluid was modeled as a viscous liquid with weak compressibility. Boundary conditions,such as a no-slip solid wall, inflow and outflow, and periodic boundaries, were employed to resemble the physical problem. A sensitivity analysis, which has been rarely addressed in previous studies, was conducted on several SPH parameters. Hence, the effects of distinct parameters, such as the kernel choices and the domain dimensions, were investigated with the goal of obtaining highly accurate results. A range of Reynolds numbers(1-500) was simulated, and the results were compared with existing experimental data. It was observed that the domain dimensions and the resolution of SPH particles, in comparison to the obstacle size, affected the obtained drag coefficient significantly. Other parameters, such as the background pressure, influenced the transient condition, but did not influence the steady state at which the drag coefficient was determined.

颗粒流运动SPH方法及滑坡破碎效应研究

山体滑坡是全球范围内多发的一种地质灾害。由于滑坡具有突发性和破坏性强的特点,建立有效的数值分析模型将有助于制定针对性的防治策略。本文针对滑坡运动过程中表现出的颗粒流特性,基于μ(I)模型提出了针对浅层滑坡的动态摩擦系数表达式,并构建了对应的光滑粒子流体动力学(SPH)方法求解框架。考虑到颗粒破碎对滑坡运动性的显著影响,结合基于破坏势能的颗粒破碎法则,建立μ(I)模型中基底摩擦力与颗粒分布之间的关系。通过两个经典的三维斜面模型试验,验证了μ(I)模型在滑坡运动分析中的应用价值,并进行了颗粒破碎效应参数敏感性分析,为后续滑坡灾害的防治工作提供参考。

Modified smoothed particle hydrodynamics approach for modelling dynamic contact angle hysteresis

Dynamic wetting plays an important role in the physics of multiphase flow, and has a significant influence on many industrial and geotechnical applications. In this work, a modified smoothed particle hydrodynamics (SPH) model is employed to simulate surface tension, contact angle and dynamic wetting effects at meso-scale. The wetting and dewetting phenomena are simulated in a capillary tube, where the liquid particles are raised or withdrawn by a shifting substrate. The SPH model is modified by introducing a newly developed viscous force formulation at the liquid-solid interface to reproduce the rate-dependent behaviour of the moving contact line. Dynamic contact angle simulations with the interfacial viscous force are conducted to verify the effectiveness and accuracy of this new formulation. In addition, the influence of interfacial viscous forces with different magnitude on the contact angle dynamics is examined by empirical power-law correlations;the derived constants suggest that the dynamic contact angle changes monotonically with the interfacial viscous force. The simulation results are consistent with experimental observations and theoretical predictions, implying that the interfacial viscous force can be associated with the slip length of flow and the microscopic surface roughness. This work demonstrates that the modified SPH model can successfully account for the rate-dependent effects of a moving contact line, and can be used for realistic multiphase flow simulations under dynamic conditions.

Coupling of discrete-element method and smoothed particle hydrodynamics for liquid-solid flows

Particle based methods can be used for both the simulations of solid and fluid phases in multiphase medium, such as the discrete-element method for solid phase and the smoothed particle hydrodynamics for fluid phase. This paper presents a computational method combining these two methods for solid-liquid medium. The two phases are coupled by using an improved model from a reported Lagrangian-Eulerian method. The technique is verified by simulating liquid-solid flows in a two-dimensional lid-driven cavity.

Numerical Simulation of Bubble Formation at a Single Orifice in Gas-fluidized Beds with Smoothed Particle Hydrodynamics and Finite Volume Coupled Method

A coupled method describing gas–solid two-phase flow has been proposed to numerically study the bubble formation at a single orifice in gas-fluidized beds.Solid particles are traced with smoothed particle hydrodynamics,whereas gas phase is discretized by finite volume method.Drag force,gas pressure gradient,and volume fraction are used to couple the two methods.The effect of injection velocities,particle sizes,and particle densities on bubble growth is analyzed using the coupled method.The simulation results,obtained for two-dimensional geometries,include the shape and diameter size of a bubble as a function of time;such results are compared with experimental data,previous numerical results,and other approximate model predictions reported in the literature.Moreover,the flow profiles of gas and particle phases and the temperature distribution by the heat transfer model around the forming bubble are also discussed.All results show that the coupled method efficiently describes of the bubble formation in fluidized beds.The proposed method is applicable for solving gas–solid two-phase flow in fluidization.

基于SPH方法的船体海洋飞沫生成形态模拟

为了研究海洋飞沫的生成形态,采用光滑粒子法(SPH)的粒子捕捉技术,分析4种常见的船首特征对海洋飞沫轨迹的影响。研究发现,船首的设计在飞沫的生成和运动特性上起着决定性作用。球鼻艏和破冰艏通过镇压波浪减少飞沫的覆盖范围,而飞剪艏和前倾艏则利用其独特的流线型设计有效地对飞沫进行引流,但却导致飞沫的生成高度较高,覆盖范围较广。此外,研究还揭示了不同类型的飞沫(如崩波破碎、卷波和激散波破碎)在船首附近的工况和轨迹模式。研究结果可为船舶设计优化及防除冰技术提供重要理论支持。

Hydrodynamic Coefficients for a 3-D Uniform Flexible Barge UsingWeakly Compressible Smoothed Particle Hydrodynamics

The numerical modelling of the interactions between water waves and floating structures is significant for different areas of the marine sector, especially seakeeping and prediction of wave-induced loads. Seakeeping analysis involving severe flow fluctuations is still quite challenging even for the conventional RANS method. Particle method has been viewed as alternative for such analysis especially those involving deformable boundary, wave breaking and fluid fragmentation around hull shapes. In this paper, the weakly compressible smoothed particle hydrodynamics(WCSPH), a fully Lagrangian particle method, is applied to simulate the symmetric radiation problem for a stationary barge treated as a flexible body. This is carried out by imposing prescribed forced simple harmonic oscillations in heave, pitch and the two-and three-node distortion modes. The resultant,radiation force predictions, namely added mass and fluid damping coefficients, are compared with results from 3-D potential flow boundary element method and 3-D RANS CFD predictions, in order to verify the adopted modelling techniques for WCSPH.WCSPH were found to be in agreement with most results and could predict the fluid actions equally well in most cases.

Application of Smooth Particle Hydrodynamics Method for Modelling Blood Flow with Thrombus Formation

Thrombosis plays a crucial role in atherosclerosis or in haemostasis when a blood vessel is injured.This article focuses on using a meshless particle-based Lagrangian numerical technique,the smoothed particles hydrodynamic(SPH)method,to study the flow behaviour of blood and to explore the flow parameters that induce formation of a thrombus in a blood vessel.Due to its simplicity and effectiveness,the SPH method is employed here to simulate the process of thrombogenesis and to study the effect of various blood flow parameters.In the present SPH simulation,blood is modelled by two sets of particles that have the characteristics of plasma and of platelet,respectively.To simulate coagulation of platelets which leads to a thrombus,the so-called adhesion and aggregation mechanisms of the platelets during this process are modelled by an inter-particle force model.The transport of platelets in the flowing blood,platelet adhesion and aggregation processes are coupled with viscous blood flow for various low Reynolds number scenarios.The numerical results are compared with the experimental observations and a good agreement is found between the simulated and experimental results.

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.

Improved kernel gradient free-smoothed particle hydrodynamics and its applications to heat transfer problems

Kernel gradient free-smoothed particle hydrodynamics （KGF-SPH） is a modified smoothed particle hydrodynamics （SPH） method which has higher precision than the conventional SPH. However, the Laplacian in KGF-SPH is approximated by the two-pass model which increases computational cost. A new kind of discretization scheme for the Laplacian is proposed in this paper, then a method with higher precision and better stability, called Improved KGF-SPH, is developed by modifying KGF-SPH with this new Laplacian model. One-dimensional （1D） and two-dimensional （2D） heat conduction problems are used to test the precision and stability of the Improved KGF-SPH. The numerical results demonstrate that the Improved KGF-SPH is more accurate than SPH, and stabler than KGF-SPH. Natural convection in a closed square cavity at different Rayleigh numbers are modeled by the Improved KGF-SPH with shifting particle position, and the Improved KGF-SPH results are presented in comparison with those of SPH and finite volume method （FVM）. The numerical results demonstrate that the Improved KGF-SPH is a more accurate method to study and model the heat transfer problems.