Abstract and Applied Analysis

Volume 2013 (2013), Article ID 131025, 7 pages

http://dx.doi.org/10.1155/2013/131025

## Lattice Boltzmann Simulation of Collision between 2D Circular Particles Suspension in Couette Flow

^{1}Institute of Fluid Engineering, Zhejiang University, Hangzhou 310027, China^{2}College of Metrology and Technology Engineering, China Jiliang University, Hangzhou 310018, China

Received 23 October 2013; Accepted 5 November 2013

Academic Editor: Jianzhong Lin

Copyright © 2013 Li-Zhong Huang and De-Ming Nie. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Collision between 2D circular particles suspension in Couette flow is simulated by using multiple-relaxation-time based lattice Boltzmann and direct forcing/fictitious domain method in this paper. The patterns of particle collisions are simulated and analyzed in detail by changing the velocity of top and bottom walls in the Couette flow. It can be seen from the simulation results that, while the velocity is large enough, the number of collisions between particles will change little as this velocity varies.

#### 1. Introduction

Multiphase flow is a very important branch of fluid mechanics, while fluid-solid two-phase flow is the main part of such flow. And fluid-solid two-phase flow is very common in nature and industry, such as raindrop formation, material science, chemical industry, aerosol deposition, fluidized beds, and injection molding machine [1, 2]. Particle collisions, which will influence the performance of product, happen occasionally in those processes. So it is very meaningful to do some research to understand particle collisions and then take control to improve the performance of product.

Even with the help of the precise and advanced particle image velocity (PIV) instrument, it is very difficult to observe the phenomenon of particle collisions in detail. Comparing the experimental measurement, numerical simulation has great advantages to investigate particle collisions in fluid-solid two-phase flow, especially for direct numerical simulation (DNS) methods. The lattice Boltzmann method (LBM), one of the best DNS methods with several remarkable advantages, was first proposed by Ladd [3] and then improved by Aidun et al. [4] to simulate particles suspended in a viscous fluid. Feng and Michaelides [5] and Tian et al. [6] united immersed boundary method (IBM) into LBM to deal with fluid-solid interface problem. In the past several decades the LBM method was proved robust and efficient for particulate flows, especially in the case of large number of particles [7–12]. Nie and Lin [13] developed a single-relaxation-time (SRT) based lattice Boltzmann-direct forcing/fictitious domain (SRT LB-DF/FD) method to simulate particle suspensions and then improved to multiple-relaxation-time (MRT LB-DF/FD).

Several papers investigated particle collisions in three-dimensional homogeneous isotropic turbulence [14–18] but not by using DNS method. In this paper, DNS method is adopted to simulate particle collisions. Because the MRT model has better ability of computing pressure and more time saving than the SRT model, the MRT LB DF/FD method is utilized. Firstly, this method is introduced in detail in Section 2. Secondly, simulation problem is described in Section 3. And finally, several simulation results and conclusions are presented in Section 4.

#### 2. Numerical Method

##### 2.1. The MRT LB DF/FD Method

The lattice Boltzmann method based on the multiple-relaxation-time (MRT) collision model is adopted in this paper [13]. The discrete equations can be written as follows: where is the diagonal collision matrix, , and is the transform matrix,

The moment space spanned by and the velocity space spanned by are related by a linear mapping, ; that is,. The D2Q9 model is used in 2D simulation, and the discrete velocities are listed as the follows:
where is the lattice speed, among which is the lattice spacing, and * is *the time step. For D2Q9 model, the corresponding nine moments are given by [19]
where is the mean density and , are the conserved moments. Moreover, the other nonconserved moments are listed as follows:

According to Guo et al. [20], the forcing term is defined as follows: where is the acceleration due to the external force and are weights related to the lattice model which are chosen as ; , ; , . The speed of sound is equal to . By applying the Taylor expansion techniques and the Chapman-Enskog analysis, (1) leads to the hydrodynamic equations that are shown as follows:

The shear viscosity and the bulk viscosity can be defined as follows [21]:

##### 2.2. Collisions Model

A collision model is needed to avoid particles overlapping. Then the short-range repulsive force developed by Wan and Turek [22] is utilized in this paper because it is easily carried out in the code. For particle-particle collisions, the short-range repulsive force is calculated by where and are the mass center coordinates of the th and th particle, and are their radius, is the distance between their mass centers which equals to , is the force range which is usually set to be one or two lattice spacing, and and are two small positive stiffness parameters for particle-particle collisions.

But for particle-wall collisions, the corresponding short-range repulsive force is determined as follows: where is the coordinate of the nearest imaginary particle located on the boundary nearby the th particle, is the distance between the imaginary particle and the real one which equals to , and and are the other two small positive stiffness parameters for particle-wall collisions. In this paper, , , , and are all set to be 1e-7.

#### 3. Problem Description

The main purpose of this paper is to study collisions between 2D circular particles suspension in Couette flow. The schematic diagram of geometry is shown in Figure 1 which is constructed by four solid walls, cm, and mesh element is 640000. The left and right walls are fixed, while the top wall moves in the right direction with a velocity of and the bottom wall moves in the left direction at the same velocity. 128 circular particles (8 rows, 16 columns) are placed in the middle of this domain, and the distance of center of two neighborhood particles in horizontal and vertical direction is two times of the diameter of particle. In order to ignore the influence of gravity, the density of particles is set to be same as the suspension fluid, g/cm^{3}. The kinematic viscosity of the fluid is set to be 0.01 cm^{2} s^{−1}. For each particle, the diameter is 0.2 cm and the arrangement of Lagrangian points is shown in Figure 2.

#### 4. Simulation Results and Conclusions

In this paper, when , the circular particles are fixed to generate a fully developed Couette flow which can be seen from Figures 3(a), 4(a), and 5(a). And when , those particles begin to move in the suspension fluid. At first, is set to be very small which equals to 0.1 m/s and the instantaneous vorticity contours can be seen from Figures 3(a)–3(i). Obviously, the vorticity is not strong enough which makes the number of collisions very big.

Instantaneous vorticity contours for m/s at several interval times are shown in Figures 4(a)–4(i). It can be seen from those figures that the particles are moving by the influence of the two main symmetry vortices. Even when , the vortex contour is rotational symmetry from the center of the simulation domain.

But when m/s, Figures 5(a)–5(i) show that the circular particles are agglomerated near the four sidewalls by the influence of the two main vortices. And when , the particles are no longer rotational symmetry from the center of the simulation domain, distributed at completely random. The vortex patterns are totally different when varies from 0.1 m/s to 0.3 m/s; however, the number of particle collisions varies little when is large enough which can be seen from Figure 6.

#### Acknowledgment

This work is supported by the National Natural Science Foundation of China with Grant 11132008.

#### References

- J. Lin, X. Shi, and Z. Yu, “The motion of fibers in an evolving mixing layer,”
*International Journal of Multiphase Flow*, vol. 29, no. 8, pp. 1355–1372, 2003. View at Publisher · View at Google Scholar · View at Scopus - L. Jianzhong, Z. Weifeng, and Y. Zhaosheng, “Numerical research on the orientation distribution of fibers immersed in laminar and turbulent pipe flows,”
*Journal of Aerosol Science*, vol. 35, no. 1, pp. 63–82, 2004. View at Publisher · View at Google Scholar · View at Scopus - A. J. C. Ladd, “Numerical simulations of particulate suspensions via a discretized Boltzmann equation. I. Theoretical foundation,”
*Journal of Fluid Mechanics*, vol. 271, pp. 285–309, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. K. Aidun, Y. Lu, and E.-J. Ding, “Direct analysis of particulate suspensions with inertia using the discrete Boltzmann equation,”
*Journal of Fluid Mechanics*, vol. 373, pp. 287–311, 1998. View at Google Scholar · View at Scopus - Z.-G. Feng and E. E. Michaelides, “The immersed boundary-lattice Boltzmann method for solving fluid-particles interaction problems,”
*Journal of Computational Physics*, vol. 195, no. 2, pp. 602–628, 2004. View at Publisher · View at Google Scholar · View at Scopus - F.-B. Tian, H. Luo, L. Zhu, J. C. Liao, and X.-Y. Lu, “An efficient immersed boundary-lattice Boltzmann method for the hydrodynamic interaction of elastic filaments,”
*Journal of Computational Physics*, vol. 230, no. 19, pp. 7266–7283, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Jianzhong, S. Xing, and Y. Zhenjiang, “Effects of the aspect ratio on the sedimentation of a fiber in Newtonian fluids,”
*Journal of Aerosol Science*, vol. 34, no. 7, pp. 909–921, 2003. View at Publisher · View at Google Scholar · View at Scopus - X. Shi, J. Lin, and Z. Yu, “Discontinuous Galerkin spectral element lattice Boltzmann method on triangular element,”
*International Journal for Numerical Methods in Fluids*, vol. 42, no. 11, pp. 1249–1261, 2003. View at Publisher · View at Google Scholar · View at Scopus - J.-Z. Lin, Y.-L. Wang, and J. A. Olsen, “Sedimentation of rigid cylindrical particles with mechanical contacts,”
*Chinese Physics Letters*, vol. 22, no. 3, pp. 628–631, 2005. View at Publisher · View at Google Scholar · View at Scopus - L. Jianzhong and K. Xiaoke, “Fiber orientation distributions in a suspension flow through a parallel plate channel containing a cylinder,”
*Journal of Composite Materials*, vol. 43, no. 12, pp. 1373–1390, 2009. View at Publisher · View at Google Scholar · View at Scopus - X. K. Ku and J. Z. Lin, “Inertial effects on the rotational motion of a fibre in simple shear flow between two bounding walls,”
*Physica Scripta*, vol. 80, no. 2, Article ID 025801, 2009. View at Publisher · View at Google Scholar · View at Scopus - X. Ku and J. Lin, “Numerical simulation of the flows over two tandem cylinders by lattice Boltzmann method,”
*Modern Physics Letters B*, vol. 19, no. 28-29, pp. 1551–1554, 2005. View at Publisher · View at Google Scholar · View at Scopus - D. Nie and J. Lin, “A LB-DF/FD method for particle suspensions,”
*Communications in Computational Physics*, vol. 7, no. 3, pp. 544–563, 2010. View at Publisher · View at Google Scholar · View at Scopus - L. I. Zaichik, O. Simonin, and V. M. Alipchenkov, “Collision rates of bidisperse inertial particles in isotropic turbulence,”
*Physics of Fluids*, vol. 18, no. 3, Article ID 035110, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. I. Zaichik, O. Simonin, and V. M. Alipchenkov, “Two statistical models for predicting collision rates of inertial particles in homogeneous isotropic turbulence,”
*Physics of Fluids*, vol. 15, no. 10, pp. 2995–3005, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L.-P. Wang, A. S. Wexler, and Y. Zhou, “On the collision rate of small particles in isotropic turbulence. I. Zero-inertia case,”
*Physics of Fluids*, vol. 10, no. 1, pp. 266–276, 1998. View at Google Scholar · View at Scopus - Y. Zhou, A. S. Wexler, and L.-P. Wang, “On the collision rate of small particles in isotropic turbulence. II. Finite inertia case,”
*Physics of Fluids*, vol. 10, no. 5, pp. 1206–1216, 1998. View at Google Scholar · View at Scopus - C. N. Franklin, P. A. Vaillancourt, M. K. Yau, and P. Bartello, “Collision rates of cloud droplets in turbulent flow,”
*Journal of the Atmospheric Sciences*, vol. 62, no. 7, pp. 2451–2466, 2005. View at Publisher · View at Google Scholar · View at Scopus - L.-S. Luo, W. Liao, X. Chen, Y. Peng, and W. Zhang, “Numerics of the lattice Boltzmann method: effects of collision models on the lattice Boltzmann simulations,”
*Physical Review E*, vol. 83, no. 5, Article ID 056710, 2011. View at Publisher · View at Google Scholar · View at Scopus - Z. Guo, C. Zheng, and B. Shi, “Discrete lattice effects on the forcing term in the lattice Boltzmann method,”
*Physical Review E*, vol. 65, no. 4, Article ID 046308, 6 pages, 2002. View at Publisher · View at Google Scholar · View at Scopus - D. d'Humières, I. Ginzburg, M. Krafczyk, P. Lallemand, and L.-S. Luo, “Multiple-relaxation-time lattice Boltzmann models in three dimensions,”
*The Royal Society of London*, vol. 360, no. 1792, pp. 437–451, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Wan and S. Turek, “An efficient multigrid-FEM method for the simulation of solid-liquid two phase flows,”
*Journal of Computational and Applied Mathematics*, vol. 203, no. 2, pp. 561–580, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet