2019/06/06 by Georgios Fourtakas, José M. Domínguez, Jose M. Dominguez +2 · 237 citations
Engineering · Mathematics · #Boundary value problem #Classical mechanics #Computational science #Computer science #Curvature #Flow (mathematics) #Fluid Dynamics Simulations and Interactions #Fluid Dynamics and Heat Transfer #Geometry #Graphics processing unit #Hagen–Poiseuille equation #Lattice Boltzmann Simulation Studies #Mathematical analysis #Mathematics #Mechanics #Parallel computing #Physics #Robustness (evolution) #Smoothed-particle hydrodynamics #Stencil
paper · pdf · doi:10.1016/j.compfluid.2019.06.009
published in Computers & Fluids 190, 346-361 (Elsevier BV)
openalex publication_date 2019/06/06 · crossref created 2019/06/06 · crossref issued 2019/08/01 · crossref published 2019/08/01 · crossref published-print 2019/08/01 · crossref deposited 2025/09/23 · openalex created_date 2025/10/10 · crossref indexed 2026/08/03 · openalex updated_date 2026/08/04
This paper presents the development of a new boundary treatment for free-surface hydrodynamics using the smoothed particle hydrodynamics (SPH) method accelerated with a graphics processing unit (GPU). The new solid boundary formulation uses a local uniform stencil (LUST) of fictitious particles that surround and move with each fluid particle and are only activated when they are located inside a boundary. This addresses the issues currently affecting boundary conditions in SPH, namely the accuracy, robustness and applicability while being amenable to easy parallelization such as on a GPU. In 3-D, the methodology uses triangles to represent the geometry with a ray tracing procedure to identify when the LUST particles are activated. A new correction is proposed to the popular density diffusion term treatment to correct for pressure errors at the boundary. The methodology is applicable to complex arbitrary geometries without the need of special treatments for corners and curvature is presented. The paper presents the results from 2-D and 3-D Poiseuille flows showing convergence rates typical for weakly compressible SPH. Still water in a complex 3-D geometry with a pyramid demonstrates the robustness of the technique with excellent agreement for the pressure distributions. The method is finally applied to the SPHERIC benchmark of a dry-bed dam-break impacting an obstacle showing satisfactory agreement and convergence for a violent flow.