A second-order accurate immersed boundary method for fully resolved simulations of particle-laden flows

详细信息	查看全文 | 推荐本文 |

• 作者：Wim-Paul Breugem ; ^{w.p.breugem@tudelft.nl}
• 关键词：Immersed boundary method ; Particle-laden flow ; Finite-volume method ; Grid locking ; Multidirect forcing scheme ; Retraction distance ; Order of grid convergence
• 刊名：Journal of Computational Physics
• 出版年：2012
• 期刊代码：136_00219991
• 类别：cp
• 出版时间：1 May, 2012
• 卷：231
• 期：13
• 页码：4469-4498
• 文件大小：4529 K

摘要

An immersed boundary method (IBM) with second-order spatial accuracy is presented for fully resolved simulations of incompressible viscous flows laden with rigid particles. The method is based on the computationally efficient direct-forcing method of Uhlmann [M. Uhlmann, An immersed boundary method with direct forcing for simulation of particulate flows, J. Comput. Phys. 209 (2005) 448-476] that is embedded in a finite-volume/pressure-correction method. The IBM consists of two grids: a fixed uniform Eulerian grid for the fluid phase and a uniform Lagrangian grid attached to and moving with the particles. A regularized delta function is used to communicate between the two grids and proved to be effective in suppressing grid locking. Without significant loss of efficiency, the original method is improved by: (1) a better approximation of the no-slip/no-penetration (ns/np) condition on the surface of the particles by a multidirect forcing scheme, (2) a correction for the excess in the effective particle diameter by a slight retraction of the Lagrangian grid from the surface towards the interior of the particles with a fraction of the Eulerian grid spacing, and (3) an enhancement of the numerical stability for particle-fluid mass density ratios near unity by a direct account of the inertia of the fluid contained within the particles. The new IBM contains two new parameters: the number of iterations N_s of the multidirect forcing scheme and the retraction distance r_d. The effect of N_s and r_d on the accuracy is investigated for five different flows. The results show that r_d has a strong influence on the effective particle diameter and little influence on the error in the ns/np condition, while exactly the opposite holds for N_s. A novel finding of this study is the demonstration that r_d has a strong influence on the order of grid convergence. It is found that for spheres the choice of r_d = 0.3螖x yields second-order accuracy compared to first-order accuracy of the original method that corresponds to r_d = 0. Finally, N_s = 2 appears optimal for reducing the error in the ns/np condition and maintaining the computational efficiency of the method.