A shallow water SPH model with PML boundaries

Abstract

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.