Simple absorbing layer conditions for shallow wave simulations with Smoothed Particle Hydrodynamics

We study and implement a simple method, based on the Perfectly Matched Layer approach, to treat non reflecting boundary conditions with the Smoothed Particles Hydrodynamics numerical algorithm. The method is based on the concept of physical damping operating on a fictitious layer added to the computational domain. The method works for both 1D and 2D cases, but here we illustrate it in the case of 1D and 2D time dependent shallow waves propagating in a finite domain

PANORMUS-SPH. A new Smoothed Particle Hydrodynamics solver for incompressible flows

A new Smoothed Particle Hydrodynamics (SPH) solver is presented, fully integrated within the PANORMUS package [7], originally developed as a Finite Volume Method (FVM) solver. The proposed model employs the fully Incompressible SPH approach, where a Fractional Step Method is used to make the numerical solution march in time. The main novelty of the proposed model is the use of a general and highly flexible procedure to account for different boundary conditions, based on the discretization of the boundary surfaces with a set of triangles and the introduction of mirror particles with suitable hydrodynamic properties. Both laminar and turbulent flows can be solved (the latter using the ε turbulence closure) and considerable flexibility in the solver algorithms is guaranteed, achieved through the use of the available Graphical User Interface. The integration of the FVM and SPH solvers within one code easily allows to develop an hybrid approach, using a simple and efficient procedure described in the paper. Model performance is tested for a series of benchmark cases for which analytical, numerical and/or experimental comparison results are available, demonstrating the ability of the solver to provide reliable solutions of incompressible flows

A shallow water SPH model with PML boundaries

We focus on the study and implementation of Smoothed Particle Hydrodynamics (SPH) numerical code to deal with non-reflecting boundary conditions, starting from the Perfect Matched Layer (PML) approach. Basically, the method exploits the concept of a physical damping which acts on a fictitious layer added to the edges of computational domain. In this paper, we develop the study of time dependent shallow waves propagating on a finite 2D-XY plane domain and their behavior in the presence of circular and, more generic, rectangular boundary absorbing layers. In particular, an analysis of variation of the layer׳s thickness versus the absorbing efficiency is conducted. In our model, the magnitude of absorbtion of a specific layer in which two types of damping functions (linear and hyperbolic) are activated is compared with the one produced by the antithetical cases of total reflecting and open boundaries. The results obtained indicate the good applicability of PML approach to SPH numerical scheme showing high absorption values with reasonable thickness of the absorbing layers