RITMARE SP1_WP2_AZ2_UO15_D04: Multi-resolution SPH 2D/3D

The present technical report is dedicated to the descrition of an algorithm for coupling SPH with an external solution. The external solution can be either another SPH solution (possibly with different discretization) or a different numerical solver or an analytical solution. The interaction between the SPH solver and the external solution is achieved through an interface region. The interface region is defined as a fixed portion of the computational domain that provides a boundary condition for the SPH solver. A ghost fluid, composed by fully lagrangian particles (i.e. ghost particles) covering the interface region, is used to impose the boundary condition. The ghost particle evolution, including its position, is integrated in time according to the field of the external solution. The physical quantities of the ghost particles needed in the integration scheme are obtained through an MLS interpolation on the field of the external solution. When a ghost particle crosses the boundary of the interface region, entering in the SPH domain, it evolves according to the SPH governing equation. The spatial distribution of the ghost particles can become largely non-uniform due to the forcing by the external solution. Thus, a packing algorithm is applied on the ghost particles in the interface region, to guarantee a particle distribution suitable for SPH operators. Since the ghost particles can exit from the interface region, a seeding algorithm is needed to introduce new ghost-particles. The algorithm is tested on several benchmarks and with the external solutions given by other SPH solvers with different discretizations and by analytical solutions. The technique is deeply investigated in terms of accuracy, efficiency and possible applications. Finally a coupled simulation involving a finite volume solver is presented.

Il presente deliverable consiste nello sviluppo di un codice di calcolo basato sul modello Smoothed Particle Hydrodynamics (SPH). Il solutore dovrà essere in grado di gestire diverse risoluzioni spaziali all'interno del dominio fluido di interesse. Essendo il modello SPH basato su una tecnica puramente lagrangiana l'utilizzo di particelle di diversa grandezza nello stesso dominio di calcolo richiede la programmazione e validazione di algoritmi dedicati che verranno sviluppati all'interno di RITMARE.