TY - JOUR

T1 - An efficient parallel simulation of interacting inertial particles in homogeneous isotropic turbulence

AU - Onishi, Ryo

AU - Takahashi, Keiko

AU - Vassilicos, J. C.

PY - 2013/6/1

Y1 - 2013/6/1

N2 - This study has conducted parallel simulations of interacting inertial particles in statistically-steady isotropic turbulence using a newly-developed efficient parallel simulation code. Flow is computed with a fourth-order finite-difference method and particles are tracked with the Lagrangian method. A binary-based superposition method has been developed and implemented in the code in order to investigate the hydrodynamic interaction among many particles. The code adopts an MPI library for a distributed-memory parallelization and is designed to minimize the MPI communication, which leads to a high parallel performance. The code has been run to obtain collision statistics of a monodisperse system with St=0.4 particles, where St is the Stokes number representing the particle relaxation time relative to the Kolmogorov time. The attained Taylor-microscale based Reynolds number R λ ranges from 54.9 to 527. The largest simulation computed the flow on 20003 grids and 10003 (one billion) particles. Numerical results have shown that the collision kernel increases for R λ<100 then decreases as R λ increases. This Reynolds dependency is attributed to that of the radial distribution function at contact, which measures the contribution of particle clustering to the collision kernel. The results have also shown that the hydrodynamic interaction for St=0.4 particles decreases both the radial relative velocity and radial distribution function at contact, leading the collision efficiency less than unity. The collision efficiency increases from 0.65 to 0.75 as R λ increases for R λ<200 and then saturates.

AB - This study has conducted parallel simulations of interacting inertial particles in statistically-steady isotropic turbulence using a newly-developed efficient parallel simulation code. Flow is computed with a fourth-order finite-difference method and particles are tracked with the Lagrangian method. A binary-based superposition method has been developed and implemented in the code in order to investigate the hydrodynamic interaction among many particles. The code adopts an MPI library for a distributed-memory parallelization and is designed to minimize the MPI communication, which leads to a high parallel performance. The code has been run to obtain collision statistics of a monodisperse system with St=0.4 particles, where St is the Stokes number representing the particle relaxation time relative to the Kolmogorov time. The attained Taylor-microscale based Reynolds number R λ ranges from 54.9 to 527. The largest simulation computed the flow on 20003 grids and 10003 (one billion) particles. Numerical results have shown that the collision kernel increases for R λ<100 then decreases as R λ increases. This Reynolds dependency is attributed to that of the radial distribution function at contact, which measures the contribution of particle clustering to the collision kernel. The results have also shown that the hydrodynamic interaction for St=0.4 particles decreases both the radial relative velocity and radial distribution function at contact, leading the collision efficiency less than unity. The collision efficiency increases from 0.65 to 0.75 as R λ increases for R λ<200 and then saturates.

KW - Homogeneous isotropic turbulence

KW - Hydrodynamic interaction

KW - Parallel computing

KW - Particle collision

UR - http://www.scopus.com/inward/record.url?scp=84875795579&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84875795579&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2013.02.027

DO - 10.1016/j.jcp.2013.02.027

M3 - Article

AN - SCOPUS:84875795579

VL - 242

SP - 809

EP - 827

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -