Improved Incompressible Smoothed Particle Hydrodynamics method for simulating flow around bluff bodies
Improved Incompressible Smoothed Particle Hydrodynamics method for simulating flow around bluff bodies
复制标题
DOI：
10.1016/j.cma.2010.12.002
复制
发表时间：
2011-02
影响因子：
7.2
通讯作者：
Mostafa Safdari Shadloo;A. Zainali;S. Sadek;M. Yildiz
中科院分区：
文献类型：
--
作者：
Mostafa Safdari Shadloo;A. Zainali;S. Sadek;M. Yildiz
在这篇文章中，我们提出了一个翼型和方形障碍物的流动数值解使用不可压缩光滑粒子流体动力学（ISPH）方法与改进的固体边界处理方法，称为多边界切线（MBT）方法。结果表明，MBT边界处理技术对于处理复杂形状的边界是非常有效的。此外，我们已经提出了使用的排斥分量的Lennard-Jones势（LJP）在平流方程中，以修复粒子断裂发生在SPH方法由于SPH粒子的趋势，遵循流线轨迹。这种方法被称为人工粒子位移法。数值结果表明，由MBT方法、人工粒子位移和修正SPH离散格式组成的改进的ISPH方法能够获得非常稳定和鲁棒的SPH模拟。对攻角在0° ~ 15°范围内的方形障碍物和NACA翼型在较高雷诺数的层流流场中进行了数值模拟。结果表明，改进的ISPH方法能够自然地捕捉到钝体流动的复杂物理现象，如流动分离、尾缘尾迹形成和旋涡脱落等。SPH的结果进行了验证与网格相关的有限元法（FEM），并观察到良好的协议之间的结果。
In this article, we present numerical solutions for flow over an airfoil and a square obstacle using Incompressible Smoothed Particle Hydrodynamics (ISPH) method with an improved solid boundary treatment approach, referred to as the Multiple Boundary Tangents (MBT) method. It was shown that the MBT boundary treatment technique is very effective for tackling boundaries of complex shapes. Also, we have proposed the usage of the repulsive component of the Lennard-Jones Potential (LJP) in the advection equation to repair particle fractures occurring in the SPH method due to the tendency of SPH particles to follow the stream line trajectory. This approach is named as the artificial particle displacement method. Numerical results suggest that the improved ISPH method which is consisting of the MBT method, artificial particle displacement and the corrective SPH discretization scheme enables one to obtain very stable and robust SPH simulations. The square obstacle and NACA airfoil geometry with the angle of attacks between 0° and 15° were simulated in a laminar flow field with relatively high Reynolds numbers. We illustrated that the improved ISPH method is able to capture the complex physics of bluff-body flows naturally such as the flow separation, wake formation at the trailing edge, and the vortex shedding. The SPH results are validated with a mesh-dependent Finite Element Method (FEM) and excellent agreements among the results were observed.