Numerical simulations of shock wave reflection phenomena in non-stationary flows using regularized smoothed particle hydrodynamics (RSPH)

In this paper we wish to demonstrate to what extent the numerical method regularized smoothed particle hydrodynamics (RSPH) is capable of modelling shocks and shock reflection patterns in a satisfactory manner. The use of SPH based methods to model shock wave problems has been relatively sparse, both due to historical reasons, as the method was originally developed for studies of astrophysical gas dynamics, but also due to the fact that boundary treatment in Lagrangian methods may be a difficult task. The boundary conditions have therefore been given special attention in this paper. Results presented for one quasi-stationary and three non-stationary flow tests reveal a high degree of similarity, when compared to published numerical and experimental data. The difference is found to be below 5, in the case where experimental data was found tabulated. The transition from regular reflection (RR) to Mach reflection (MR) and the opposite transition from MR to RR are studied. The results are found to be in close agreement with the results obtained from various empirical and semi-empirical formulas published in the literature. A convergence test shows a convergence rate slightly steeper than linear, comparable to what is found for other numerical methods when shocks are involved.     

An enhanced treatment of boundary conditions in implicit smoothed particle hydrodynamics简

In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme （Implicit SPH）. Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares （MLS） interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.     

基于SPH方法的表面张力修正算法及其应用

针对传统SPEI方法中基于CSF模型的表面张力算法,在计算边界、尖角等粒子缺失部位的曲率时存在偏差较大,且粒子秩序较差,对大变形问题表面张力计算精度较低的问题,在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到了表面张力修正方程组.应用本文方法模拟了水溶液中初始...     

Numerical simulation of soil–water interaction using smoothed particle hydrodynamics (SPH) method

An application of smoothed particle hydrodynamics (SPH) to simulation of soil–water interaction is presented. In this calculation, water is modeled as a viscous fluid with week compressibility and soil is modeled as an elastic–perfectly plastic material. The Mohr–Coulomb failure criterion is applied to describe the stress states of soil in the plastic flow regime. Dry soil is modeled by one-phase flow while saturated soil is modeled by separate water and soil phases. Interaction between soil and water is taken into account by means of pore water pressure and seepage force. Simulation tests of soil excavation by a water jet are calculated as a challenging example to verify the broad applicability of the SPH method. The excavations are carried out in two different soil models, one is dry soil and the other is fully saturated soil. Numerical results obtained in this paper have shown that the gross discontinuities of soil failure can be simulated without any difficulties. This supports the feasibility and attractiveness of this a new approach in geomechanics applications. Advantages of the method are robustness, conceptual simplicity and relative ease of incorporating new physics.     

Répartition des contraintes au sein d'une plaque d'athérome coronarienne chez l'homme obtenue à partir d'un modèle 3D aux éléments finisModelling of the stress distribution in an artherosclerotic plaque in man using a Finite Element Analysis

To quantify the stress distribution in an artherosclerotic plaque, it is usual to perform a finite element analysis (FEA) of the cross section having the highest stenosis severity. The purpose of this study is to compare the stress distributions obtained from the real 3D-plaque geometry with those obtained from the 2D FEA. Our calculations shown that for a complex and heterogeneous plaque, the 2D FEA of several cross sections is needed in order to be in agreement with the 3D FEA results, and to localize the region in which plaque rupture may occur. To cite this article: F. Treyve et al., C. R. Mecanique 331 (2003).