Simple Iterative Incompressible Smoothed Particle Hydrodynamics

08/06/2019
by   Abhinav Muta, et al.
0

In this paper a simple, robust, and general purpose approach to implement the Incompressible Smoothed Particle Hydrodynamics (ISPH) method is proposed. The new approach is well suited for implementation on CPUs and GPUs. The method is matrix-free and uses an iterative formulation to setup and solve the pressure Poisson equation. A novel approach is used to ensure homogeneous particle distributions and improved boundary conditions. The new formulation enables the use of solid wall boundary conditions from the weakly-compressible SPH schemes. The new method is fast and runs on GPUs without the need for complex integration with sparse linear solvers. Several benchmark problems that illustrate the robustness, performance, and wide range of applicability of the new scheme are demonstrated. An open source implementation is provided and the manuscript is fully reproducible.

READ FULL TEXT

page 28

page 31

research
08/23/2022

How to train your solver: Verification of boundary conditions for smoothed particle hydrodynamics

The weakly compressible smoothed particle hydrodynamics (WCSPH) method h...
research
09/10/2019

PySPH: a Python-based framework for smoothed particle hydrodynamics

PySPH is an open-source, Python-based, framework for particle methods in...
research
01/28/2020

Automatic weak imposition of free slip boundary conditions via Nitsche's method: application to nonlinear problems in geodynamics

Imposition of free slip boundary conditions in science and engineering s...
research
01/26/2020

A Variational Staggered Particle Framework for Incompressible Free-Surface Flows

Smoothed particle hydrodynamics (SPH) has been extensively studied in co...
research
11/30/2022

Corner Cases of the Generalized Tau Method

Polynomial spectral methods provide fast, accurate, and flexible solvers...
research
05/17/2023

An unambiguous and robust formulation for Wannier localization

We provide a new variational definition for the spread of an orbital und...

Please sign up or login with your details

Forgot password? Click here to reset