Research Papers: Multiphase Flows

Flow Modulation by Finite-Size Neutrally Buoyant Particles in a Turbulent Channel Flow

[+] Author and Article Information
Lian-Ping Wang

Department of Mechanical Engineering,
University of Delaware,
Newark, DE 19716-3140
e-mail: lwang@udel.edu

Cheng Peng

Department of Mechanical Engineering,
University of Delaware,
Newark, DE 19716-3140
e-mail: cpengxpp@udel.edu

Zhaoli Guo

National Laboratory of Coal Combustion,
Huazhong University of Science and Technology,
Wuhan 430074, China
e-mail: zlguo@hust.edu.cn

Zhaosheng Yu

Department of Mechanics,
Zhejiang University,
Hangzhou 310027, China
e-mail: yuzhaosheng@zju.edu.cn

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received January 11, 2015; final manuscript received August 22, 2015; published online December 8, 2015. Assoc. Editor: E. E. Michaelides.

J. Fluids Eng 138(4), 041306 (Dec 08, 2015) (15 pages) Paper No: FE-15-1016; doi: 10.1115/1.4031691 History: Received January 11, 2015; Revised August 22, 2015

A fully mesoscopic, multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) is developed to perform particle-resolved direct numerical simulation (DNS) of wall-bounded turbulent particle-laden flows. The fluid–solid particle interfaces are treated as sharp interfaces with no-slip and no-penetration conditions. The force and torque acting on a solid particle are computed by a local Galilean-invariant momentum exchange method. The first objective of the paper is to demonstrate that the approach yields accurate results for both single-phase and particle-laden turbulent channel flows, by comparing the LBM results to the published benchmark results and a full-macroscopic finite-difference direct-forcing (FDDF) approach. The second objective is to study turbulence modulations by finite-size solid particles in a turbulent channel flow and to demonstrate the effects of particle size. Neutrally buoyant particles with diameters 10% and 5% the channel width and a volume fraction of about 7% are considered. We found that the mean flow speed was reduced due to the presence of the solid particles, but the local phase-averaged flow dissipation was increased. The effects of finite particle size are reflected in the level and location of flow modulation, as well as in the volume fraction distribution and particle slip velocity near the wall.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Corsini, A. , Rispoli, F. , Sheard, A. G. , and Venturini, P. , 2013, “ Numerical Simulation of Coal Fly-Ash Erosion in an Induced Draft Fan,” ASME J. Fluids Eng., 135(1), p. 081303. [CrossRef]
Patro, P. , and Dash, S. K. , 2014, “ Computations of Particle-Laden Turbulent Jet Flows Based on Eulerian Model,” ASME J. Fluids Eng., 136(1), p. 011301. [CrossRef]
Balachandar, S. , and Eaton, J. K. , 2010, “ Turbulent Dispersed Multiphase Flow,” Annu. Rev. Fluid Mech., 42(1), pp. 111–133. [CrossRef]
Maxey, M. R. , and Riley, J. J. , 1983, “ Equation of Motion for a Small Rigid Sphere in a Nonuniform Flow,” Phys. Fluids, 26(4), pp. 883–889. [CrossRef]
Wang, L.-P. , Rosa, B. , Gao, H. , He, G. W. , and Jin, G. D. , 2009, “ Turbulent Collision of Inertial Particles: Point-Particle Based, Hybrid Simulations and Beyond,” Int. J. Multiphase Flow, 35(9), pp. 854–867. [CrossRef]
Squires, K. D. , and Eaton, J. K. , 1991, “ Preferential Concentration of Particles by Turbulence,” Phys. Fluids A, 3(5), pp. 1169–1179. [CrossRef]
Wang, L.-P. , and Maxey, M. R. , 1993, “ Settling Velocity and Concentration Distribution of Heavy Particles in Homogeneous Isotropic Turbulence,” J. Fluid Mech., 256, pp. 27–68. [CrossRef]
Squires, K. D. , and Eaton, J. D. , 1990, “ Particle Response and Turbulence Modification in Isotropic Turbulence,” Phys. Fluids A, 2(7), pp. 1191–1203. [CrossRef]
Elghobashi, S. , and Truesdell, G. , 1993, “ On the 2-Way Interaction Between Homogeneous Turbulence and Dispersed Solid Particles. I: Turbulence Modification,” Phys. Fluids A, 5(7), pp. 1790–1801. [CrossRef]
Sundaram, S. , and Collins, L. R. , 1997, “ Collision Statistics in an Isotropic Particle-Laden Turbulent Suspension. Part 1—Direct Numerical Simulations,” J. Fluid Mech., 335, pp. 75–109. [CrossRef]
Zhou, Y. , Wexler, A. S. , and Wang, L.-P. , 2001, “ Modeling Turbulent Collision of Bidisperse Inertial Particles,” J. Fluid Mech., 433, pp. 77–104. [CrossRef]
Ayala, O. , Grabowski, W. W. , and Wang, L.-P. , 2007, “ A Hybrid Approach for Simulating Turbulent Collisions of Hydrodynamically-Interacting Particles,” J. Comput. Phys., 225(1), pp. 51–73. [CrossRef]
Jin, G. D. , He, G. W. , and Wang, L.-P. , 2010, “ Large Eddy Simulation of Collisional Statistics of Inertial Particles in Isotropic Turbulence,” Phys. Fluids., 22(5), p. 055106. [CrossRef]
Burton, T. M. , and Eaton, J. K. , 2005, “ Fully Resolved Simulations of Particle–Turbulence Interaction,” J. Fluid Mech., 545, pp. 67–111. [CrossRef]
Uhlmann, M. , 2005, “ An Immersed Boundary Method With Direct Forcing for the Simulation of Particulate Flows,” J. Comput. Phys. 209(2), pp. 448–476. [CrossRef]
Uhlmann, M. , 2008, “ Interface-Resolved Direct Numerical Simulation of Vertical Particulate Channel Flow in the Turbulent Regime,” Phys. Fluids, 20(5), p. 053305. [CrossRef]
Shao, X. , Wu, T. , and Yu, Z. , 2012, “ Fully Resolved Numerical Simulation of Particle-Laden Turbulent Flow in a Horizontal Channel at a Low Reynolds Number,” J. Fluid Mech., 693, pp. 319–344. [CrossRef]
Yang, J. M. , and Stern, F. , 2014, “ A Sharp Interface Direct Forcing Immersed Boundary Approach for Fully Resolved Simulations of Particulate Flows,” ASME J. Fluids Eng., 136(4), p. 040904. [CrossRef]
Zhang, Z. , and Prosperetti, A. , 2005, “ A Second-Order Method for Three-Dimensional Particle Simulation,” J. Comput. Phys., 210(1), pp. 292–324. [CrossRef]
Yeo, K. , Dong, S. , Climent, E. , and Maxey, M. R. , 2010, “ Modulation of Homogeneous Turbulence Seeded With Finite Size Bubbles or Particles,” Int. J. Multiphase Flow, 36(3), pp. 221–233. [CrossRef]
Homann, H. , and Bec, J. , 2010, “ Finite-Size Effects in the Dynamics of Neutrally Buoyant Particles in Turbulent Flow,” J. Fluid Mech., 651, pp. 81–91. [CrossRef]
Gao, H. , Li, H. , and Wang, L.-P. , 2013, “ Lattice Boltzmann Simulation of Turbulent Flow Laden With Finite-Size Particles,” Comput. Math. Appl., 65(2), pp. 194–210. [CrossRef]
Wang, L.-P. , Ayala, O. , Gao, H. , Andersen, C. , and Mathews, K. , 2014, “ Study of Forced Turbulence and Its Modulation by Finite-Size Solid Particles Using the Lattice Boltzmann Approach,” Comput. Math. Appl., 67(2), pp. 363–380. [CrossRef]
Gillissen, J. , 2013, “ Turbulent Drag Reduction Using Fluid Spheres,” J. Fluid Mech., 716, pp. 83–95. [CrossRef]
Ten Cate, A. , Derksen, J. J. , Portela, L. M. , and van den Akker, H. E. A. , 2004, “ Fully Resolved Simulations of Colliding Monodisperse Spheres in Forced Isotropic Turbulence,” J. Fluid Mech., 519, pp. 233–271. [CrossRef]
Ernst, M. , Dietzel, M. , and Sommerfeld, M. , 2013, “ A Lattice Boltzmann Method for Simulating Transport and Agglomeration of Resolved Particles,” Acta Mech., 224(10), pp. 2425–2449. [CrossRef]
Pan, Y. , and Banerjee, S. , 1997, “ Numerical Investigation of the Effects of Large Particles on Wall Turbulence,” Phys. Fluids, 9(12), pp. 3786–3807. [CrossRef]
Vowinckel, B. , Kempe, T. , and Froehlich, J. , 2014, “ Fluid–Particle Interaction in Turbulent Open Channel Flow With Fully-Resolved Mobile Beds,” Adv. Water Resour., 72, pp. 32–44. [CrossRef]
Rashidi, M. , Hetstoni, G. , and Banerjee, S. , 1990, “ Particle Turbulence Interaction in a Boundary Layer,” Int. J. Multiphase Flow, 16(6), pp. 935–949. [CrossRef]
Kaftori, D. , Hetsroni, G. , and Banerjee, S. , 1998, “ The Effect of Particles on Wall Turbulence,” Int. J. Multiphase Flow, 24(3), pp. 359–386. [CrossRef]
Rogers, C. B. , and Eaton, J. K. , 1991, “ The Effect of Small Particles on Fluid Turbulence in a Flat-Plate, Turbulent Boundary Layer in Air,” Phys. Fluids, 3(5), pp. 928–937. [CrossRef]
Zisselmar, R. , and Molerus, O. , 1979, “ Investigation of Solid–Liquid Pipe Flow With Regard to Turbulence Modification,” Chem. Eng. J., 18(3), pp. 233–239. [CrossRef]
Belt, R. J. , Daalmans, A. C. L. M. , and Portela, M. , 2012, “ Experimental Study of Particle-Driven Secondary Flow in Turbulent Pipe Flows,” J. Fluid Mech., 709, pp. 1–36. [CrossRef]
Gore, R. A. , and Crowe, C. T. , 1989, “ Effect of Particle Size on Modulating Turbulent Intensity,” Int. J. Multiphase Flow, 15(2), pp. 279–285. [CrossRef]
Hetsroni, G ., 1989, “ Particle Turbulence Interaction,” Int. J. Multiphase Flow, 15(5), pp. 735–746. [CrossRef]
Matas, J. P. , Morris, J. F. , and Guazzelli, E. , 2003, “ Transition to Turbulence in Particulate Pipe Flow,” Phys. Rev. Lett., 90(1), p. 014501. [CrossRef] [PubMed]
Yu, Z. , and Shao, X. , 2007, “ A Direct-Forcing Fictitious Domain Method for Particulate Flows,” J. Comput. Phys., 227(1), pp. 292–314. [CrossRef]
Chen, S. , and Doolen, G. , 1998, “ Lattice Boltzmann Method for Fluid Flows,” Annu. Rev. Fluid Mech., 30(1), pp. 329–364. [CrossRef]
Qian, Y. H. , d’Humières, D. , and Lallemand, P. , 1992, “ Lattice BGK Models for Navier–Stokes Equation,” Europhys. Lett., 17(6BIS), pp. 479–484. [CrossRef]
Peng, Y. , Liao, W. , Luo, L.-S. , and Wang, L.-P. , 2010, “ Comparison of the Lattice Boltzmann and Pseudo-Spectral Methods for Decaying Turbulence: Low-Order Statistics,” Comput. Fluids, 39(4), pp. 568–591. [CrossRef]
Lammers, P. , Beronov, K. N. , Volkert, R. , Brenner, G. , and Durst, F. , 2006, “ Lattice BGK Direct Numerical Simulation of Fully Developed Turbulence in Incompressible Plane Channel Flow,” Comput. Fluids, 35(10), pp. 1137–1153. [CrossRef]
Peng, C. , Teng, Y. , Hwang, B. , Guo, Z. , and Wang, L.-P. , 2015, “ Implementation Issues and Benchmarking of Lattice Boltzmann Method for Moving Particle Simulations in a Viscous Flow,” Comput. Math. Appl. (in press).
d’Humières, D. , Ginzburg, I. , Krafczyk, M. , Lallemand, P. , and Luo, L.-S. , 2002, “ Multiple-Relaxation-Time Lattice Boltzmann Models in Three-Dimensions,” Philos. Trans. R. Soc. London A, 360(1792), pp. 437–451. [CrossRef]
Guo, Z. , Zheng, C. , and Shi, B. , 2002, “ Discrete Lattice Effects on the Forcing Term in the Lattice Boltzmann Method,” Phys. Rev. E., 65(4), p. 046308. [CrossRef]
Lu, J. , Han, H. , Shi, B. , and Guo, Z. , 2012, “ Immersed Boundary Lattice Boltzmann Model Based on Multiple Relaxation Times,” Phys. Rev. E., 85(1), p. 016711. [CrossRef]
He, X. , and Luo, L.-S. , 1997, “ Theory of the Lattice Boltzmann Method: From the Boltzmann Equation to the Lattice Boltzmann Equation,” Phys. Rev. E., 56(6), pp. 6811–6817. [CrossRef]
Shan, X. W. , Yuan, X. F. , and Chen, H. D. , 2006, “ Kinetic Theory Representation of Hydrodynamics: A Way Beyond the Navier–Stokes Equation,” J. Fluid Mech., 550, pp. 413–441. [CrossRef]
Lallemand, P. , and Luo, L.-S. , 2003, “ Lattice Boltzmann Method for Moving Boundaries,” J. Comput. Phys., 184(2), pp. 406–421. [CrossRef]
Caiazzo, A. , 2008, “ Analysis of Lattice Boltzmann Nodes Initialisation in Moving Boundary Problems,” Prog. Comput. Fluid Dyn., 8(1–4), pp. 3–10. [CrossRef]
Wen, B. , Zhang, C. , Tu, Y. , Wang, C. , and Fang, H. , 2014, “ Galilean Invariant Fluid–Solid Interfacial Dynamics in Lattice Boltzmann Simulations,” J. Comput. Phys., 266, pp. 161–170. [CrossRef]
Nguyen, N. Q. , and Ladd, A. J. C. , 2002, “ Lubrication Corrections for Lattice-Boltzmann Simulations of Particle Suspensions,” Phys. Rev. E, 66(4), p. 046708. [CrossRef]
Dance, S. L. , and Maxey, M. R. , 2003, “ Incorporation of Lubrication Effects into the Force-Coupling Method for Particulate Two-Phase Flow,” J. Comput. Phys., 189(1), pp. 212–238. [CrossRef]
Glowinski, R. , Pan, T.-W. , Hesla, T. I. , and Joseph, D. D. , 1999, “ A Distributed Lagrange Multiplier/Fictitious Domain Method for Particulate Flows,” Int. J. Multiphase Flow, 25(5), pp. 755–794. [CrossRef]
Wang, L.-P. , Peng, C. , Guo, Z. , and Yu, Z. , 2015, “ Lattice Boltzmann Simulation of Particle-Laden Turbulent Channel Flow,” Comput. Fluids (in press).
Wu, T. , Shao, X. , and Yu, Z. , 2011, “ Fully Resolved Numerical Simulation of Turbulent Pipe Flows Laden With Large Neutrally-Buoyant Particles,” J. Hydrodyn. Ser. B, 23(1), pp. 21–25. [CrossRef]
Yu, Z. , Wu, T. , Shao, X. , and Lin, J. , 2013, “ Numerical Studies of the Effects of Large Neutrally Buoyant Particles on the Flow Instability and Transition to Turbulence in Pipe Flow,” Phys. Fluids, 25(4), p. 043305. [CrossRef]
Kim, J. , Moin, P. , and Moser, R. , 1987, “ Turbulence Statistics in Fully-Developed Channel Flow at Low Reynolds-Number,” J. Fluid Mech., 177, pp. 133–166. [CrossRef]
Moser, R. , Kim, J. , and Mansour, N. N. , 1999, “ Direct Numerical Simulation of Turbulent Channel Flow up to Re-Tau = 590,” Phys. Fluids., 11(4), pp. 943–945. [CrossRef]
Picano, F. , Breugem, W. P. , and Brandt, L. , 2015, “ Turbulent Channel Flow of Dense Suspensions of Neutrally Buoyant Spheres,” J. Fluid Mech., 764, pp. 463–487. [CrossRef]
Li, Y. M. , McLaughlin, J. B. , Kontomaris, K. , and Portela, L. , 2001, “ Numerical Simulation of Particle-Laden Turbulent Channel Flow,” Phys. Fluids, 13(10), pp. 2957–2967. [CrossRef]
Segré, G. , and Silberberg, A. , 1961, “ Radial Particle Displacements in Poiseuille Flow of Suspensions,” Nature, 189(476), pp. 209–210. [CrossRef]
Segré, G. , and Silberberg, A. , 1962, “ Behavior of Macroscopic Rigid Spheres in Poiseuille Flow,” J. Fluid Mech., 14(1), pp. 136–157. [CrossRef]
Shao, X. , Yu, Z. , and Sun, B. , 2008, “ Inertial Migration of Spherical Particles in Circular Poiseuille Flow at Moderately High Reynolds Numbers,” Phys. Fluids, 20(10), p. 103307. [CrossRef]


Grahic Jump Location
Fig. 1

Sketches of (a) the coordinate system used for the channel flow simulation and (b) the two-dimensional domain decomposition for message passing interface parallel implementation

Grahic Jump Location
Fig. 4

Comparison of the RMS velocity profiles of the simulated single-phase turbulent channel flows

Grahic Jump Location
Fig. 7

The relative change in the mean velocity of the flow due to the presence of particles: (a) linear–linear plot and (b) log–linear plot

Grahic Jump Location
Fig. 5

The time evolution of mean flow speed (averaged over y) in the LBM simulations. The two vertical dash lines indicate the stationary stage (32.2<t*<56.1) used to obtain average statistics for LBM–PLL. The two blue vertical lines indicate the stationary stage (31.9<t*<53.9) used to obtain average statistics for LBM–PLS. The three horizontal lines mark 15.74, 15.02, and 14.82, respectively.

Grahic Jump Location
Fig. 2

Turbulent Reynolds shear stress profiles in half of the channel. All quantities are normalized by u*2. The thin diagonal line denotes the total shear stress, so the difference between this straight line and the data represents the viscous shear stress due to the mean flow.

Grahic Jump Location
Fig. 3

Comparison of the mean flow velocity profiles of the simulated single-phase turbulent channel flows

Grahic Jump Location
Fig. 8

The mean velocity profiles from the LBM simulations and fitting coefficients in the inertial sublayer by the logarithmic law. The two vertical lines mark the region where linear regression is performed to obtain a fit of the form U+=lny+/k+A. The two vertical lines mark the region (30<y+<130) used to perform the logarithmic fitting.

Grahic Jump Location
Fig. 16

Comparison of the phase-partitioned mean velocity profiles: (a) Log-linear plot and (b) linear-linear plot

Grahic Jump Location
Fig. 6

The mean velocity profiles from both the LBM and FDDF simulations with a domain size of 4H×2H×2H : (a) log–linear plot and (b) linear–linear plot

Grahic Jump Location
Fig. 9

The RMS velocity profiles in the flows simulated by LBM: (a) linear–linear plot and (b) log–linear plot

Grahic Jump Location
Fig. 10

The RMS velocity profiles in the flows simulated by the FDDF approach: (a) linear–linear plot and (b) log–linear plot

Grahic Jump Location
Fig. 11

The relative changes of RMS velocity fluctuations due to the presence of solid particles in the LBM simulations, relative to the single-phase flow simulated by the same method using a same domain size: (a) linear–linear plot and (b) log–linear plot

Grahic Jump Location
Fig. 12

The relative changes of RMS velocity fluctuations due to the presence of solid particles in the FDDF simulations, relative to the single-phase flow simulated by the same method using a same domain size: (a) linear–linear plot and (b) log–linear plot

Grahic Jump Location
Fig. 13

The average local volume fraction Ψp of the particulate phase as a function of y+: (a) linear–linear plot and (b) log–linear plot. The profile from Picano et al. [59] at ϕV=0.1 and ap/H=1/18=0.056 is also shown for comparison. The vertical lines mark the location of y = ap.

Grahic Jump Location
Fig. 15

The average forces normalized by [ρf(u∗)2dp2], acting on the particle along the transverse direction as a function of (a) y+ and (b) y/ap

Grahic Jump Location
Fig. 17

Visualization of instantaneous spanwise vorticity ωz, normalized by u*2/ν, on an xy slice: (a) LBM–SP, (b) LBM–PLL, and (c) LBM–PLS

Grahic Jump Location
Fig. 14

Replot of Fig. 13 as a function of y/ap: (a) log–linear plot and (b) linear–linear plot. The profiles from Uhlmann [16] and Picano et al. [59] are also shown for comparison.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In