The δplus-SPH model: Simple procedures for a further improvement of the SPH scheme

• Published in: Computer methods in applied mechanics and engineering Vol. 315; pp. 25 - 49
• Main Authors: Sun, P.N., Colagrossi, A., Marrone, S., Zhang, A.M.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.03.2017; Elsevier BV
• Subjects: [formula omitted]-SPH; Compressibility; Continuity (mathematics); Continuity equation; Dissipation; Fragmentation; Free-surface flows; High frequencies; Irregularities; Mathematical models; Particle shifting; Smoothed particle hydrodynamics; Stability; Tensile instability; Weakly-compressibility; Particle shifting; Smoothed particle hydrodynamics; Tensile instability; Free-surface flows; δ-SPH; Weakly-compressibility
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2016.10.028

The present work is dedicated to the improvement of the δ-SPH scheme. This is an enhanced weakly-compressible SPH model widely used in recent years thanks to its benefits to the standard SPH scheme, to its low CPU costs and to its ease of implementation. Nonetheless, the δ-SPH still presents some drawbacks as other SPH models. For example, in some critical conditions it does not prevent the tensile instability and the consequent numerical fragmentation. Furthermore, even if the use of a diffusive term in the SPH continuity equation is able to reduce numerical high frequencies on the pressure field, the velocity gradients are generally noisy because of the irregularities of the particle spatial configurations, which, in specific flow conditions, can induce also extra numerical-dissipation. For these reasons a particle shifting technique is used to improve the model and a special treatment has been developed for particles that are close to the free-surface region. The introduction of the particle-shifting procedure is generalized in the context of multi-resolutions for which a novel algorithm is formulated to handle the particle re-positioning in the different resolution levels. The proposed algorithms can be straightforwardly implemented in an SPH model without requiring cumbersome code modifications. The δ+-SPH is validated on seven different benchmarks giving a wide panorama on the improvements of this new SPH model. •The present work is dedicated to the improvement of the δ-SPH scheme.•A particle shifting technique is introduced in the context of weakly-compressible models.•A special treatment for particles that are close to the free-surface region is developed.•The particle-shifting procedure is generalized for multi-resolutions approaches.•The δ+-SPH is validated on seven different benchmarks.