A correction scheme for two-way coupled point-particle simulations on anisotropic grids

The accuracy of Lagrangian point-particle models for simulation of particle-laden flows may degrade when the particle and fluid momentum equations are two-way coupled. The exchange of force between the fluid and particle changes the velocity of the fluid at the location of the particle, thereby modifying the slip velocity and producing an erroneous prediction of coupling forces between fluid and particle. In this article, we propose a correction scheme to reduce this error and predict the undisturbed fluid velocity accurately. Conceptually, in this method, the computation cell is treated as a solid object immersed in the fluid that is subjected to the two-way coupling force and dragged at a velocity that is identical to the disturbance created by the particle. The proposed scheme is generic as it can be applied to unstructured grids with arbitrary geometry, particles that have different size and density, and arbitrary interpolation scheme. In its crudest form for isotropic grids, the present correction scheme reduces to dividing the Stokes drag by 1−0.75Λ for Λ≤1, where Λ is the ratio of the particle diameter to the grid size. The accuracy of the proposed scheme is evaluated by comparing the computed settling velocity of an individual and pair of particles under gravity on anisotropic rectilinear grids against analytical solutions. This comparison shows up to two orders of magnitude reduction in error in cases where the particle is up to 5 times larger than the grid that may have an aspect ratio of over 10. Furthermore, a comparison against a particle-resolved simulation of decaying isotropic turbulence demonstrates the excellent accuracy of the proposed scheme.