scholarly journals Application Research of an Efficient and Stable Boundary Processing Method for the SPH Method

Water ◽  
2019 ◽  
Vol 11 (5) ◽  
pp. 1110 ◽  
Author(s):  
Xing Huang ◽  
Wu Chen ◽  
Zhe Hu ◽  
Xing Zheng ◽  
Shanqin Jin ◽  
...  

The boundary truncation of the kernel function affects the numerical accuracy and calculation stability of the smooth particle hydrodynamics (SPH) method and has been one of the key research fields for this method. In this paper, an efficient and stable boundary processing method for the SPH method was introduced by adopting an improved boundary interpolation method (i.e., the improved Shepard method) which needs only the sum of direct accumulation for fixed-boundary particles to improve the numerical stability and computational efficiency of the fixed ghost particle method. The improvement effect of the method was demonstrated by comparing it with different interpolation methods using the cases of still water, a wave generated by dam-breaking, and a solitary wave attacking problem with fixed walls and a moveable wall. The results showed that the new boundary processing method for SPH can help remarkably improve the efficiency of calculation and reduce the oscillations of pressure when simulating various flows.

2011 ◽  
Vol 94-96 ◽  
pp. 1638-1641 ◽  
Author(s):  
Gui Ming Rong ◽  
Hiroyuki Kisu

In the smoothed particle hydrodynamics (SPH) method, the particle inconsistency problem significantly influences the calculation accuracy. In the present study, we investigate primarily the influence of the particle inconsistency on the first derivative of field functions and discuss the behavior of several methods of addressing this problem. In addition, we propose a new approach by which to compensate for this problem, especially for functions having a non-zero second derivative, that is less computational demanding, as compared to the finite particle method (FPM). A series of numerical studies have been carried out to verify the performance of the new approach.


2018 ◽  
Vol 15 (03) ◽  
pp. 1850011 ◽  
Author(s):  
Xiufeng Yang ◽  
Song-Charng Kong

The purpose of this paper is to present and illustrate a smoothed particle hydrodynamics (SPH) method to study the process of a drop impacting on a dry solid surface. SPH is a Lagrangian mesh-free particle method that offers advantages in modeling the evolution of the liquid surface during drop impact. A new surface tension model is used. The artificial viscosity is also used, which is demonstrated to be, approximately, a linear function of the dynamic viscosity of the liquid. The SPH method is used to simulate different liquid drops impacting on dry surfaces. The numerical results agree with experimental data obtained from the literature. The influence of various parameters on the drop impact, including impact velocity, diameter, viscosity, surface tension, and density of the drop, is also studied. The results show that the dimensionless spreading diameter of the drop increases if the impact velocity, diameter, or density increases, while the increase in viscosity and surface tension decreases the spreading diameter. The results indicate that the drop impact depends more strongly on the viscosity and impact velocity than on the diameter, surface tension, and density of the drop. In addition to the impact of a spherical drop, the impact of an ellipsoidal drop on a dry surface is also studied. The results show that the aspect ratio of the drop has a significant influence on the outcome of drop impact.


2006 ◽  
Vol 306-308 ◽  
pp. 595-600
Author(s):  
Fei Xu ◽  
Masanori Kikuchi

Smoothed Particle Hydrodynamics (SPH) is a relatively new technique for simulating the dynamic response of solids, especially for high velocity impact and fracture problem. However, closer examination of SPH reveals some severe problems. The major difficulties are: (1) tensile instability; (2) zero-energy mode; (3) boundary deficiency; (4) less accuracy. One solution to these major difficulties with SPH is to improve the consistency of the kernel function. Based on the Reproducing Kernel Particle Method (RKPM), the concept of the proposed simplified linear consistency is introduced. The most attractive feature of the simplified linear consistency is the ease and cheapness of doing 3D calculation. One contribution of this paper is to show clearly the accuracy of solution gradually improved by increasing the order of the consistency. Simple 3D impacting models are established with different geometries and higher accurate results are obtained by using higher consistency kernel functions. Other features as numerical convergence, computational efficiency, etc. and some considerations of the simplified linear consistency kernel function are also discussed.


Author(s):  
Zhi-Gang Feng ◽  
Miguel Cortina Ponton

Smoothed Particle Hydrodynamics (SPH) is a meshless Lagrangian method in which the domain is represented by particles. Each particle is assigned properties such as mass, pressure, density, temperature, and velocity. These properties are then evaluated at the particle positions using a smoothing kernel that integrates over the values of the surrounding particles. In the present study, the SPH method is first used to study the flow between infinite plates at different temperatures; velocity and temperature profiles are compared to the analytical solution and good agreements are found. Then, the SPH is further extended to study the heat transfer of an isothermal cylinder in a periodic flow. The velocity profile is found to be in good agreement with results provided in the literature, and the temperature field obtained seems reasonable. Finally, the natural convection of a heated cylinder in a square enclosure is studied. For this, the discretized version of the heat equation is used in conjunction with the Boussinesq approximation. The results between the SPH with and without a tensile instability are compared.


Water ◽  
2019 ◽  
Vol 11 (11) ◽  
pp. 2314 ◽  
Author(s):  
Shu Wang ◽  
Anping Shu ◽  
Matteo Rubinato ◽  
Mengyao Wang ◽  
Jiping Qin

Non-homogeneous viscous debris flows are characterized by high density, impact force and destructiveness, and the complexity of the materials they are made of. This has always made these flows challenging to simulate numerically, and to reproduce experimentally debris flow processes. In this study, the formation-movement process of non-homogeneous debris flow under three different soil configurations was simulated numerically by modifying the formulation of collision, friction, and yield stresses for the existing Smoothed Particle Hydrodynamics (SPH) method. The results obtained by applying this modification to the SPH model clearly demonstrated that the configuration where fine and coarse particles are fully mixed, with no specific layering, produces more fluctuations and instability of the debris flow. The kinetic and potential energies of the fluctuating particles calculated for each scenario have been shown to be affected by the water content by focusing on small local areas. Therefore, this study provides a better understanding and new insights regarding intermittent debris flows, and explains the impact of the water content on their formation and movement processes.


Author(s):  
Amin Rahmat ◽  
Hossein Nasiri ◽  
Marjan Goodarzi ◽  
Ehsan Heidaryan

Purpose This paper aims to introduce a numerical investigation of aquatic locomotion using the smoothed particle hydrodynamics (SPH) method. Design/methodology/approach To model this problem, a simple improved SPH algorithm is presented that can handle complex geometries using updatable dummy particles. The computational code is validated by solving the flow over a two-dimensional cylinder and comparing its drag coefficient for two different Reynolds numbers with those in the literature. Findings Additionally, the drag coefficient and vortices created behind the aquatic swimmer are quantitatively and qualitatively compared with available credential data. Afterward, the flow over an aquatic swimmer is simulated for a wide range of Reynolds and Strouhal numbers, as well as for the amplitude envelope. Moreover, comprehensive discussions on drag coefficient and vorticity patterns behind the aquatic are made. Originality/value It is found that by increasing both Reynolds and Strouhal numbers separately, the anguilliform motion approaches the self-propulsion condition; however, the vortices show different pattern with these increments.


2014 ◽  
Author(s):  
Zhen Chen ◽  
Li Zou ◽  
Zhi Zong

In this paper, the impact pressures of two different base forms are comparatively studied using Smoothed Particle Hydrodynamics (SPH) method. It is summarized from previous works that the improved weakly compressible SPH model shows better performances than incompressible SPH model in numerical simulations of free surface flows accompany with large deformations and strong discontinuities. Such advantages are observed in numerical accuracy, stability and efficiency. The weakly compressible SPH model used in this paper is equipped with some new correction algorithms, among which include the density reinitialization algorithm and a new coupled dynamic Solid Boundary Treatment (SBT) on solid boundaries. The new boundary treatment combines the advantages of both the repulsive boundary treatment and the dynamic boundary treatment, intending to obtain more stable and accurate numerical results. A benchmark test of dam breaking is conducted to prove the reliability of the numerical model used in this paper. Two representative cases, among which one has one cavity and the other one has three cavities, are numerically investigated and compared to support the conclusion that the base form with cavities generally experience lower local and overall impact pressures than the base form of flat plate. It is found that with the application of cavities on the bottom, the peak values of the boundary pressure near central bottom significantly decrease, leading to smaller force load and better structural stability. The mechanisms of such phenomenon might be the pressure absorption effect conducted by the cavities.


2020 ◽  
Vol 105 (4) ◽  
pp. 1119-1147
Author(s):  
G. Chaussonnet ◽  
T. Dauch ◽  
M. Keller ◽  
M. Okraschevski ◽  
C. Ates ◽  
...  

AbstractThis paper illustrates recent progresses in the development of the smoothed particle hydrodynamics (SPH) method to simulate and post-process liquid spray generation. The simulation of a generic annular airblast atomizer is presented, in which a liquid sheet is fragmented by two concentric counter swirling air streams. The accent is put on how the SPH method can bridge the gap between the CAD geometry of a nozzle and its characterization, in terms of spray characteristics and dynamics. In addition, the Lagrangian nature of the SPH method allows to extract additional data to give further insight in the spraying process. First, the sequential breakup events can be tracked from one large liquid blob to very fine stable droplets. This is herein called the tree of fragmentation. From this tree of fragmentation, abstract quantities can be drawn such as the breakup activity and the fragmentation spectrum. Second, the Lagrangian coherent structures in the turbulent flow can be determined easily with the finite-time Lyapunov exponent (FTLE). The extraction of the FTLE is particularly feasible in the SPH framework. Finally, it is pointed out that there is no universal and ultimate non-dimensional number that can characterize airblast primary breakup. Depending on the field of interest, a non-dimensional number (e.g. Weber number) might be more appropriate than another one (e.g. momentum flux ratio) to characterize the regime, and vice versa.


2018 ◽  
Vol 10 (4) ◽  
pp. 159-169
Author(s):  
Zhu Yue ◽  
Jiang Shengyao ◽  
Yang Xingtuan ◽  
Duan Riqiang

The moving particle semi-implicit method is a meshless particle method for incompressible fluid and has proven useful in a wide variety of engineering applications of free-surface flows. Despite its wide applicability, the moving particle semi-implicit method has the defects of spurious unphysical pressure oscillation. Three various divergence approximation formulas, including basic divergence approximation formula, difference divergence approximation formula, and symmetric divergence approximation formula are proposed in this paper. The proposed three divergence approximation formulas are then applied for discretization of source term in pressure Poisson equation. Two numerical tests, including hydrostatic pressure problem and dam-breaking problem, are carried out to assess the performance of different formulas in enhancing and stabilizing the pressure calculation. The results demonstrate that the pressure calculated by basic divergence approximation formula and difference divergence approximation formula fluctuates severely. However, application of symmetric divergence approximation formula can result in a more accurate and stabilized pressure.


2022 ◽  
Author(s):  
Binghui Cui ◽  
Liaojun Zhang

Abstract Flow-type landslide is one type of landslide that generally exhibits characteristics of high flow velocities, long jump distances, and poor predictability. Simulation of it facilitates propagation analysis and provides 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 algorithm considering the friction is proposed, and integrated into the boundary condition of the SPH method, and its accuracy is verified. Moreover, the Navier-Stokes equation combined with the non-Newtonian fluid rheology model 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 in comparison with 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.


Sign in / Sign up

Export Citation Format

Share Document