Log In Sign Up

Efficient and Accurate Adaptive Resolution for Weakly-Compressible SPH

by   Abhinav Muta, et al.

In this paper we propose an accurate, and computationally efficient method for incorporating adaptive spatial resolution into weakly-compressible Smoothed Particle Hydrodynamics (SPH) schemes. Particles are adaptively split and merged in an accurate manner while ensuring that the number of particles is not large for a given resolution. Critically, the method ensures that the number of neighbors of each particle is optimal, leading to an efficient algorithm. A set of background particles is used to specify either geometry-based spatial resolution or solution-based adaptive resolution. This allows us to simulate problems using particles having length variations of the order of 1:250 with much fewer particles than currently reported with other techniques. The method is designed to automatically adapt when any solid bodies move. The algorithms employed are fully parallel. We consider a suite of benchmark problems to demonstrate the accuracy of the approach. We then consider the classic problem of the flow past a circular cylinder at a range of Reynolds numbers and show that the proposed method produces accurate results with a significantly reduced number of particles. We provide an open source implementation and a fully reproducible manuscript.


page 23

page 27

page 28

page 35


An improved particle packing algorithm for complex geometries in SPH

Accurate simulation of fluid flow in the presence of complex geometries ...

A fast Chebyshev method for the Bingham closure with application to active nematic suspensions

Continuum kinetic theories provide an important tool for the analysis an...

New results on particle filters with adaptive number of particles

In this paper, we present new results on particle filters with adaptive ...

Adaptive symplectic model order reduction of parametric particle-based Vlasov-Poisson equation

High-resolution simulations of particle-based kinetic plasma models typi...

Gathering by Repulsion

We consider a repulsion actuator located in an n-sided convex environmen...

Generalized smoothed particle hydrodynamics with overset methods in total Lagrangian formulations

This study proposes a generalized coordinates based smoothed particle hy...

LEoPart: a particle library for FEniCS

This paper introduces LEoPart, an add-on for the open source finite elem...