Research Papers: Multiphase Flows

Study of Local Turbulence Profiles Relative to the Particle Surface in Particle-Laden Turbulent Flows

[+] Author and Article Information
Lian-Ping Wang

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

Oscar G. C. Ardila

School of Mechanical Engineering,
Universidad del Valle,
Cali 12345, Colombia
e-mail: oscar.gerardo.castro@correounivalle.edu.co

Orlando Ayala

Department of Engineering Technology,
Old Dominion University,
Norfolk, VA 23529
e-mail: oayala@odu.edu

Hui Gao

United Technologies Research Center,
East Hartford, CT 06118
e-mail: gaoh@utrc.utc.com

Cheng Peng

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

1Corresponding author.

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

J. Fluids Eng 138(4), 041307 (Dec 08, 2015) (14 pages) Paper No: FE-15-1081; doi: 10.1115/1.4031692 History: Received February 01, 2015; Revised August 22, 2015

As particle-resolved simulations (PRSs) of turbulent flows laden with finite-size solid particles become feasible, methods are needed to analyze the simulated flows in order to convert the simulation data to a form useful for model development. In this paper, the focus is on turbulence statistics at the moving fluid–solid interfaces. An averaged governing equation is developed to quantify the radial transport of turbulent kinetic energy when viewed in a frame moving with a solid particle. Using an interface-resolved flow field solved by the lattice Boltzmann method (LBM), we computed each term in the transport equation for a forced, particle-laden, homogeneous isotropic turbulence. The results illustrate the distributions and relative importance of volumetric source and sink terms, as well as pressure work, viscous stress work, and turbulence transport. In a decaying particle-laden flow, the dissipation rate and kinetic energy profiles are found to be self-similar.

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


Balachandar, S. , and Eaton, J. K. , 2010, “ Turbulent Dispersed Multiphase Flow,” Annu. Rev. Fluid Mech. 42(1), pp. 111–133. [CrossRef]
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]
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]
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]
Crowe, C. T. , 2000, “ On Models for Turbulence Modulation in Fluid-Particle Flows,” Int. J. Multiphase Flow, 26(5), pp. 719–727. [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]
Wang, L.-P. , Peng, C. , Guo, Z. , and Yu, Z. , 2015, “ Flow Modulation by Finite-Size Neutrally Buoyant Particles in a Turbulent Channel Flow,” ASME J. Fluids Eng. (in press).
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).
Lucci, F. , Ferrante, A. , and Elghobashi, S. , 2010, “ Modulation of Isotropic Turbulence by Particles of Taylor Length-Scale Size,” J. Fluid Mech., 650, pp. 5–55. [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]
d'Humières, D. , Ginzburg, I. , Krafczyk, M. , Lallemand, P. , and Luo, L.-S. , 2002, “ Multiple-Relaxation-Time Lattice Boltzmann Models in Three-Dimensions,” Phil. Trans. R. Soc. London A, 360(1792), pp. 437–451. [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).
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]
Lallemand, P. , and Luo, L.-S. , 2003, “ Lattice Boltzmann Method for Moving Boundaries,” J. Comp. Phys., 184(2), pp. 406–421. [CrossRef]
Chen, L. , Yu, Y. , Lu, J. , and Hou, G. , 2014, “ A Comparative Study of Lattice Boltzmann Methods Using Bounce-Back Schemes and Immersed Boundary Ones for Flow Acoustic Problems,” Int. J. Numer. Meth. Fluids, 74(6), pp. 439–467. [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]
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]
Crowe, C. T. , Schwarzkopf, J. D. , Sommerfeld, M. , and Tsuji, Y. , 2012, Multiphase Flows With Droplets and Particles, 2nd ed., CRC Press, Boca Raton.
Fox, R. O. , 2012, “ Large-Eddy-Simulation Tools for Multiphase Flows,” Annu. Rev. Fluid Mech., 44(1), pp. 47–76. [CrossRef]
Mei, R. , Luo, L.-S. , Lallemand, P. , and d'Humières, D. , 2006, “ Consistent Initial Conditions for Lattice Boltzmann Simulations,” Comput. Fluids, 35(8–9), pp. 855–862. [CrossRef]
Yu, H. , Luo, L.-S. , and Girimaji, S. S. , 2006, “ LES of Turbulent Square Jet Flow Using an MRT Lattice Boltzmann Model,” Comput. Fluids, 35(8–9), pp. 957–965. [CrossRef]


Grahic Jump Location
Fig. 1

Time evolution of high-order statistics of turbulence for both particle-laden and unladen flow: (a) skewness and (b) flatness

Grahic Jump Location
Fig. 2

Zoom-in view of vorticity contour and particle location on a plane-cut of z = 255.5 in the 5123 simulation: (c) 0.47Te,0 and (d) 0.79Te,0. Note that the presence of particles is associated with high vorticity values (represented by the colors toward the red end in the online version (or the dark end in black and white)), indicating relatively large dissipation near particle surfaces. The corresponding vorticity contours for single-phase flow at the two times are shown in (a) and (b), respectively.

Grahic Jump Location
Fig. 3

Time evolution of (a) Taylor microscale Reynolds number and (b) normalized dissipation rate for both particle-laden and unladen flow

Grahic Jump Location
Fig. 4

Normalized (a) dissipation rate and (b) turbulent kinetic energy as a function of distance from the center of the solid particle. The horizontal line marks the level of 1.0.

Grahic Jump Location
Fig. 5

(a) Contour of the normalized local strain rate defined as ε/〈ε〉fluid and particle positions in an x–y plane at z = 128.5, from a PRS of three-dimensional, particle-laden, forced homogeneous isotropic turbulence. This figure is taken from Ref.[25]. Only a quarter of the plane is shown. (b) Profiles conditioned on the particle surface. Each is normalized by its respective field mean.

Grahic Jump Location
Fig. 6

Profiles of volumetric source and sink terms conditioned on the particle surface. All quantities are normalized by the fluid-phase average dissipation rate.

Grahic Jump Location
Fig. 7

Profiles of turbulent transport terms conditioned on the particle surface. All quantities are normalized by the fluid-phase average dissipation rate.

Grahic Jump Location
Fig. 8

Profile of viscous stress work conditioned on the particle surface. It is normalized by the fluid-phase average dissipation rate.

Grahic Jump Location
Fig. 9

Profiles of the pressure work and other terms in the balance equation (Eq. (22)). All quantities are normalized by the fluid-phase average dissipation rate.

Grahic Jump Location
Fig. 10

(a) The PDF of the bin node numbers and (b) a replot of Fig. 9 when weighted by the PDF in (a)

Grahic Jump Location
Fig. 11

Profile of bin-averaged 〈eijeij〉 as a function of distance from the particle center, obtained from two grid resolutions

Grahic Jump Location
Fig. 12

The normalized (a) dissipation rate and (b) turbulent kinetic energy as a function of the distance from the center of the nearest solid particle. The horizontal line marks the level of 1.0.

Grahic Jump Location
Fig. 13

Profiles of bin-averaged (a) dissipation rate and (b) turbulent kinetic energy as a function of distance from the center of a particle. Each is normalized by the field mean of the particle-laden flow at the same time instant.

Grahic Jump Location
Fig. 14

Profiles conditioned on the particle surface for the Stokes disturbance flow



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