DualSPHysics model

Beach, Dike

velocity

pressure

density

position

Navier-Stokes
Lagrangian derivative
discrete summation

Fluid mass balance (continuity)

Fluid momentum and energy balance (Navier-Stokes)

Equation of state (Tait’s, weakly compressible)

Viscosity closure (artificial viscosity, laminar)

Turbulence closure (Sub-Particle Scale)

Interpolation procedure (kernel function)

Free surface (intrinsic, no need of tracking)

Dynamic boundary conditions (Crespo et al., 2007)

Floating objects (Monaghan et al., 2003)

Periodic open boundaries (Gomez-Gesteira et al., 2012)

Coupling with Discrete Element Method (Canelas et al., 2016)

Multi-phase: Two-phase liquid-sediment implementation in DualSPHysics (Fourtakas and Rogers, 2016)

Measured data:

Steering signal (wave generator)

Solid properties (in case of floating objects)

For validation:

Water surface elevation

Wave forces

Profiles of velocity

Model generated data:

Binary format by default (but also ascii, csv, text, vtk)

Single and double precision results

No time-consuming mesh generation necessary

No higher order time schemes available (2nd order has bugs)

Mature but continuously under development, unsuitable for external users (poorly documented)

Boundary conditions are still an open problem with the SPH community

Convergence

Time-consuming but GPU helps a lot!!

Parameters

Initial particle interspace

Kernel length

Coefficient of sound

Viscosity

delta-SPH

Procedure

1-2: usually known (depending on the problem)

3-5: start from recommended values, modify if necessary

Procedure

Comparison of time series

Validation test cases available

Analytical solutions

Wave flume experiments

Wave tank experiments

OWC experiments

Sloshing experiments

Debris flows