Efficient 3D DNS of gas-solid flows on Fermi GPGPU

Three-dimensional (3D) gas-solid Direct Numerical Simulation (DNS) requires huge computational resources which imposes a great challenge to both current hardware and software conditions. In this article, an efficient implementation of 3D gas-solid DNS with Lattice Boltzmann Method and Discrete Element Method is developed on a Fermi GPGPU. An Immersed Moving Boundary approach is utilized to impose the no-slip condition at particle-fluid interfaces. Optimization strategies such as changing the sequence of collision and propagation in grid evolution and making multiple kernel executing concurrently are discussed in detail. This algorithm is demonstrated to be competitive both in terms of accuracy and performance. Approximately 131 Millions of Lattices Update Per Second have been achieved, indicating that this GPGPU implementation is very suitable for 3D gas-solid DNS. (C) 2012 Elsevier Ltd. All rights reserved