Search

This paper presents a divergence cleaning formulation for the velocity in the weakly compressible smoothed particle hydrodynamics (SPH) scheme. The proposed hyperbolic/parabolic divergence cleaning, ensures that the velocity divergence, $div(\mathbf{u})$, is minimised throughout the simulation. The divergence equation is coupled with the momentum conservation equation through a scalar field $\psi$. A parabolic term is added to the time-evolving divergence equation, resulting in a hyperbolic/parabolic form, dissipating acoustic waves with a speed of sound proportional to the local Mach number in order to maximise dissipation of the velocity divergence, preventing unwanted diffusion of the pressure field. The $div(\mathbf{u})$-SPH algorithm is implemented in the open-source weakly compressible SPH solver DualSPHysics. The new formulation is validated against a range of challenging 2-D test cases including the Taylor-Green vortices, patch impact test, jet impinging on a surface, and wave impact in a sloshing tank. The results show that the new formulation reduces the divergence in the velocity field by at least one order of magnitude which prevents spurious numerical noise and the formation of unphysical voids. The temporal evolution of the impact pressures shows that the $div(\mathbf{u})$-SPH formulation virtually eliminates unwanted acoustic pressure oscillations. Investigation of particle resolution confirms that the new $div(\mathbf{u})$-SPH formulation does not reduce the spatial convergence rate.

Inspired by transformation optics and photonic crystals, this paper presents a computational investigation into the interaction between water surface waves and array waveguides of cylinders with multiple previously unexplored lattice geometries, including, for the first time, quasiperiodic geometries. Extending beyond conventional square and hexagonal periodic arrays, transformation optics has opened up entirely new opportunities to investigate water wave propagation through arrays based on quasiperiodic lattices, and quasiperiodically arranged vacancy defects. Using the linear potential flow open-source code Capytaine, missing element and $\tau$-scaled Fibonacci square lattices, the Penrose lattice, hexagonal $H_{00}$ lattice and Amman-Beenker lattice are investigated. The existence of band gaps for all arrays is observed. An hexagonal lattice with vacancy defects transmits the least energy. Bragg diffraction consistent with rotational symmetry is observed from all arrays. Waves will distort significantly to achieve resonance with arrays, supporting transformation-based waveguides. The possible uses include adaptation to more versatile waveguides with applications such as offshore renewable energy and coastal defence., Comment: 23 pages, 6 figures

For turbulent bubbly flows, multi-phase simulations resolving both the liquid and bubbles are prohibitively expensive in the context of different natural phenomena. One example is breaking waves, where bubbles strongly influence wave impact loads, acoustic emissions, and atmospheric-ocean transfer, but detailed simulations in all but the simplest settings are infeasible. An alternative approach is to resolve only large scales, and model small scale bubbles adopting sub-resolution closures. Here we introduce a large eddy simulation (LES) Smoothed Particle Hydrodynamics (SPH) scheme for simulations of bubbly flows. The continuous liquid phase is resolved with a semi-implicit isothermally compressible SPH framework. This is coupled with a discrete Lagrangian bubble model. Bubbles and liquid interact via exchanges of volume and momentum, through turbulent closures, bubble breakup and entrainment, and free-surface interaction models. By representing bubbles as individual particles, they can be tracked over their lifetimes, allowing closure models for sub-resolution fluctuations, bubble deformation, breakup and free-surface interaction in integral form, accounting for the finite timescales over which these events occur. We investigate two flows: bubble plumes, and breaking waves, and find close quantitative agreement with published experimental and numerical data. In particular, for plunging breaking waves, our framework accurately predicts the Hinze scale, bubble size distribution, and growth rate of the entrained bubble population. This is the first coupling of an SPH framework with a discrete bubble model, with potential for cost effective simulations of wave-structure interactions and more accurate predictions of wave impact loads., Comment: Submitted to JFM June 2022, this revision submitted May 2023

DualSPHysics is a weakly compressible smoothed particle hydrodynamics (SPH) Navier-Stokes solver initially conceived to deal with coastal engineering problems, especially those related to wave impact with coastal structures. Since the first release back in 2011, DualSPHysics has shown to be robust and accurate for simulating extreme wave events along with a continuous improvement in efficiency thanks to the exploitation of hardware such as graphics processing units (GPUs) for scientific computing or the coupling with wave propagating models such as SWASH and OceanWave3D. Numerous additional functionalities have also been included in the DualSPHysics package over the last few years which allow the simulation of fluid-driven objects. The use of the discrete element method (DEM) has allowed the solver to simulate the interaction among different bodies (sliding rocks, for example), which provides a unique tool to analyse debris flows. In addition, the recent coupling with other solvers like Project Chrono or MoorDyn has been a milestone in the development of the solver. Project Chrono allows the simulation of articulated structures with joints, hinges, sliders and springs and MoorDyn allows simulating moored structures. Both functionalities make DualSPHysics one of the meshless model world leaders in the simulation of offshore energy harvesting devices. Lately, the present state of maturity of the solver goes beyond single phase simulations, allowing multi-phase simulations with gas-liquid and a combination of Newtonian and non-Newtonian models expanding further the capabilities and range of applications for the DualSPHysics solver. These advances and functionalities make DualSPHysics a state-of-the-art meshless solver with emphasis on free-surface flow modelling., Comment: Comp. Part. Mech. (2021)

The semi-analytical wall boundary conditions present a mathematically rigorous framework to prescribe the influence of solid walls in SPH for fluid flows. In this paper they are investigated with respect to the skew-adjoint property which implies exact energy conservation. It will be shown that this property holds only in the limit of the continuous SPH approximation, whereas in the discrete SPH formulation it is only approximately true, leading to numerical noise. This noise, interpreted as form of "turbulence", is treated using an additional volume diffusion term in the continuity equation which we show is equivalent to an approximate Riemann solver. Subsequently two extensions to the boundary conditions are presented. The first dealing with a variable driving force when imposing a volume flux in a periodic flow and the second showing a generalization of the wall boundary condition to Robin type and arbitrary-order interpolation. Two modifications for free-surface flows are presented for the volume diffusion term as well as the wall boundary condition. In order to validate the theoretical constructs numerical experiments are performed showing that the present volume flux term yields results with an error 5 orders of magnitude smaller then previous methods while the Robin boundary conditions are imposed correctly with an error depending on the order of the approximation. Furthermore, the proposed modifications for free-surface flows improve the behaviour at the intersection of free surface and wall as well as prevent free-surface detachment when using the volume diffusion term. Finally, this paper is concluded by a simulation of a dam break over a wedge demonstrating the improvements proposed in this paper., Comment: 30 pages, 12 figures, accepted for publication in Computer Physics Communications

Starting from the single graphics processing unit (GPU) version of the Smoothed Particle Hydrodynamics (SPH) code DualSPHysics, a multi-GPU SPH program is developed for free-surface flows. The approach is based on a spatial decomposition technique, whereby different portions (sub-domains) of the physical system under study are assigned to different GPUs. Communication between devices is achieved with the use of Message Passing Interface (MPI) application programming interface (API) routines. The use of the sorting algorithm radix sort for inter-GPU particle migration and sub-domain halo building (which enables interaction between SPH particles of different subdomains) is described in detail. With the resulting scheme it is possible, on the one hand, to carry out simulations that could also be performed on a single GPU, but they can now be performed even faster than on one of these devices alone. On the other hand, accelerated simulations can be performed with up to 32 million particles on the current architecture, which is beyond the limitations of a single GPU due to memory constraints. A study of weak and strong scaling behaviour, speedups and efficiency of the resulting programis presented including an investigation to elucidate the computational bottlenecks. Last, possibilities for reduction of the effects of overhead on computational efficiency in future versions of our scheme are discussed., Comment: Accepted by the Journal of Parallel and Distributed Computing July 2012. Available online 20 August 2012. http://dx.doi.org/10.1016/j.jpdc.2012.07.010