Modified dynamic boundary conditions (mDBC) for general-purpose smoothed particle hydrodynamics (SPH): application to tank sloshing, dam break and fish pass problems

Dynamic boundary conditions (DBC) for solid surfaces are standard in the weakly compressible smoothed particle hydrodynamics (SPH) code DualSPHysics. A stationary solid is simply represented by fixed particles with pressure from the equation of state. Boundaries are easy to set up and computations are relatively stable and efficient, providing robust numerical simulation for complex geometries. However, a small unphysical gap between the fluid and solid boundaries can form, decreasing the accuracy of pressures measured on the boundary. A method is presented where the density of solid particles is obtained from ghost positions within the fluid domain by linear extrapolation. With this approach, the gap between fluid and boundary is reduced and pressures in still water converge to hydrostatic, including the case of a bed with a sharp corner. The violent free-surface cases of a sloshing tank and dam break impact on an obstacle show pressures measured directly on solid surfaces in close agreement with experiments. The complex 3-D flow in a fish pass, with baffles to divert the flow, is simulated showing close agreement with measured water levels with weirs open and gates closed, but less close with gates open and weirs closed. This indicates the method is suitable for rapidly varying free-surface flows, but development for complex turbulent flows is necessary. The code with the modified dynamic boundary condition (mDBC) is available in DualSPHysics to run on CPUs or GPUs