Graphics processing unit-accelerated smoothed particle hydrodynamics—Finite difference method and the application for the flow around a cylinder with forced motions

Sijie Wang, Hu Jun, Can Huang, Yong Yu

Physics of Fluids · 2021 · 22 citations · 56 references

Concepts

TL;DR

The finite difference method offers high accuracy and efficiency on grids, while smoothed particle hydrodynamics provides meshless handling of free surfaces and moving boundaries. The study develops a hybrid FDM–SPH coupling to leverage the strengths of both methods. The hybrid method is implemented on GPUs using a unified device architecture, with particle shifting to preserve particle distribution, and applied to low‑Reynolds‑number cylinder flows (static, rotating, oscillating, and translating) to study velocity effects. The GPU‑accelerated hybrid accurately reproduces benchmark results and shows that a translating cylinder’s velocity accelerates separation, speeds vortex shedding, and increases forces, demonstrating its potential as an effective flow‑control strategy.

Abstract

The finite difference method (FDM) is a widely used grid-based method with high computational accuracy and efficiency. Smoothed particle hydrodynamics (SPH) is a meshless particle method, which can easily deal with the fluid flow with the free surface and moving boundary conditions. Considering the different advantages of the two methods, a hybrid approach is developed by coupling FDM with SPH. To improve the computational efficiency, compute unified device architecture is used to develop a graphics processing unit parallel computing model for the SPH-FDM coupling method. At low Reynolds number, some flow cases (static, rotation, and oscillation) of a cylinder are simulated. The particle shifting technology is used to maintain the uniform distribution of the particles when the cylindrical boundary moves. The results obtained by the coupling method are in good agreement with reference results, which indicates the acceptable accuracy and stability of the coupling method. Then, the flow past a translating cylinder with a constant velocity is simulated, and the influence of the translating velocity on the flow field is discussed. The results indicate that the translational velocity can accelerate the flow separation near the boundary layer of the object, speed up the vortex shedding process, and enhance the forces acting on the object. Therefore, the translating velocity can be used as an effective flow control method.

References

56