Efficient and Accurate Adaptive Resolution for Weakly-Compressible SPH
arXiv:2107.01276 · doi:10.1016/j.cma.2022.115019
Abstract
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. 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, where the resolution is a function of distance to a solid body, or solution-based adaptive resolution, where the resolution is a function of the computed solution. 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.
arXiv comment51 pages, 40 figures, 4 tables
References in corpus (11)
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- Smoothed Particle Hydrodynamics in Astrophysics
- A smoothed particle hydrodynamics method for evaporating multiphase flows
- Adaptive resolution for multiphase smoothed particle hydrodynamics
- PySPH: a Python-based framework for smoothed particle hydrodynamics
- Entropically Damped Artificial Compressibility for SPH
- A Consistent Spatially Adaptive Smoothed Particle Hydrodynamics Method for Fluid-Structure Interactions
- Algorithms for uniform particle initialization in domains with complex boundaries
- automan: a simple, Python-based, automation framework for numerical computing
- How to train your solver: A method of manufactured solutions for weakly-compressible SPH
- A corrected transport-velocity formulation for fluid and structural mechanics with SPH
Cited by in corpus (5)
- Techniques for second order convergent weakly-compressible smoothed particle hydrodynamics schemes without boundaries
- How to train your solver: Verification of boundary conditions for smoothed particle hydrodynamics
- Robust Solid Boundary Treatment for Compressible Smoothed Particle Hydrodynamics
- Adaptive Compressible Smoothed Particle Hydrodynamics
- Rapid Variable Resolution Particle Initialization for Complex Geometries