Fast free-surface detection and level-set function definition in SPH solvers

Marrone, S.; Colagrossi, A.; Touzé, David Le; Graziani, G.

doi:10.1016/j.jcp.2010.01.019

Cited by 230 publications

(137 citation statements)

References 7 publications

Supporting

Mentioning

123

Contrasting

Unclassified

Order By: Relevance

“…This technique was tested on the triangular wedge case (Section 4.3). Other techniques exist to track the free-surface (see [28] for example), but the present one proved sufficient to ensure the impermeability of the walls while solving properly the pressure Poisson equation.…”

Section: New Isph Algorithmmentioning

confidence: 96%

Unified semi-analytical wall boundary conditions applied to 2-D incompressible SPH

Leroy

Violeau

Ferrand³

et al. 2014

Journal of Computational Physics

108

View full text Add to dashboard Cite

This work aims at improving the 2-D incompressible SPH model (ISPH) by adapting it to the unified semi-analytical wall boundary conditions proposed by Ferrand et al. [10]. The ISPH algorithm considered is as proposed by Lind et al. [25], based on the projection method with a divergence-free velocity field and using a stabilising procedure based on particle shifting. However, we consider an extension of this model to Reynolds-Averaged Navier-Stokes equations based on the k − ǫ turbulent closure model, as done in [10]. The discrete SPH operators are modified by the new description of the wall boundary conditions. In particular, a boundary term appears in the Laplacian operator, which makes it possible to accurately impose a von Neumann pressure wall boundary condition that corresponds to impermeability. The shifting and free-surface detection algorithms have also been adapted to the new boundary conditions. Moreover, a new way to compute the wall renormalisation factor in the frame of the unified semi-analytical boundary conditions is proposed in order to decrease the computational time. We present several verifications to the present approach, including a lid-driven cavity, a water column collapsing on a wedge and a periodic schematic fish-pass. Our results are compared to Finite Volumes methods, using Volume of Fluids in the case of free-surface flows. We briefly investigate the convergence of the method and prove its ability to model complex free-surface and turbulent flows. The results are generally improved when compared to a weakly compressible SPH model with the same boundary conditions, especially * Corresponding author. tel: +33 (0)6 67 88 92 13Email addresses: agnes.leroy@edf.fr (A. Leroy), damien.violeau@edf.fr (D. Violeau),

show abstract

Section: New Isph Algorithmmentioning

confidence: 96%

Unified semi-analytical wall boundary conditions applied to 2-D incompressible SPH

Leroy

Violeau

Ferrand³

et al. 2014

Journal of Computational Physics

108

View full text Add to dashboard Cite

show abstract

“…Actually, in truly incompressible formulations the pressure solution is obtained through an implicit method involving linear system solution, and requires one to impose the dynamic free-surface condition at the system boundary. In a meshless context the latter means to first detect the free surface, which is not an easy operation; see, e.g., [10]. Conversely, in weakly compressible Lagrangian formulations a fully explicit method is used and free-surface conditions are supposed to be implicitly verified.…”

Section: Governing Equations a Field Equationsmentioning

confidence: 99%

Theoretical analysis and numerical verification of the consistency of viscous smoothed-particle-hydrodynamics formulations in simulating free-surface flows

et al. 2011

View full text Add to dashboard Cite

The theoretical formulation of the smoothed particle hydrodynamics (SPH) method deserves great care because of some inconsistencies oceurring when considering free-surface inviscid flows. Actually, in SPH formulations one usually assumes that (i) surface integral terms on the boundary of the interpolation kernel support are neglected, (ii) free-surface conditions are implicitly verified. These assumptions are studied in detail in the present work for free-surface Newtonian viscous flow. The consistency of classical viscous weakly compressible SPH formulations is investigated. In particular, the principie of virtual work is used to study the verification of the free-surface boundary conditions in a weak sense. The latter can be related to the global energy dissipation induced by the viscous term formulations and their consistency. Numerical verification of this theoretical analysis is provided on three free-surface test cases including a standing wave, with the three viscous term formulations investigated.

show abstract

“…Thus this method is easy to implement and has a cheap computational expense since it is using the existing neighboring lists of particles and the same kernel function for interpolation. The freesurface tracking technique is similar to the algorithm proposed by Marrone et al (2010) and is used at each time step of numerical simulation without a significant increase of the computational time.…”

Section: Free-surface Particles Detectionmentioning

confidence: 99%

SPH Modeling of Mean Velocity Circulation in a Rip Current System

Farahani

Dalrymple

Hérault

et al. 2012

Int. Conf. Coastal. Eng.

View full text Add to dashboard Cite

A Lagrangian numerical model called Smoothed ParticleHydrodynamics is used to analyze rip current system generated by a single bar and a rip channel. The pattern of the wave-induced circulation cell over the bar, the oppositely-rotating circulation cell on-shore and a strong seaward-directed current in the rip channel is modeled numerically. The mean horizontal variations of rip current system as well as three-dimensional circulations are studied. The results in three-dimensional space reveal the wave-current interaction and flow patterns in different parts of rip channel, bar, and the trough located near shore. For comparison to experimental data, Eulerian nodes are introduced to the numerical model and SPH interpolation over neighboring Lagrangian particles is implemented to find fluid parameters at those specific nodes. This methodology leads to a better understanding of depth-integrated flows and a more accurate comparison of numerical results with experimental results. Model predictions are compared to laboratory measurements of Drønen et al. (2002) and show good agreement, including mean velocity profiles, mean surface elevation and three-dimensional velocity components.

show abstract

Fast free-surface detection and level-set function definition in SPH solvers

Cited by 230 publications

References 7 publications

Unified semi-analytical wall boundary conditions applied to 2-D incompressible SPH

Unified semi-analytical wall boundary conditions applied to 2-D incompressible SPH

Theoretical analysis and numerical verification of the consistency of viscous smoothed-particle-hydrodynamics formulations in simulating free-surface flows

SPH Modeling of Mean Velocity Circulation in a Rip Current System

Contact Info

Product

Resources

About