### Abstract

We present results of direct numerical simulations of travelling waves in dense assemblies of monodisperse spherical particles fluidized by a liquid. The cases we study have been derived from the experimental work of others. In these simulations, the flow of interstitial fluid is solved by the lattice- Boltzmann method (LBM) and the particles move under the influence of gravity, hydrodynamic forces stemming from the LBM, subgrid-scale lubrication forces and hard-sphere collisions. We first show that the propagating inhomogeneous structures seen in the simulations are in agreement with those observed experimentally. We then use the detailed information contained in the simulation results to assess aspects of two-fluid model closures, namely, fluid-particle drag, and the various contributions to the effective stresses. We show that the rates of compaction and dilation of the particle phase in the travelling waves are comparable to the rate at which the microstructure relaxes, and that there is a pronounced effect of the rate of compaction on the average collisional normal stress. Although this effect can be expressed as an effective bulk viscosity term, this approach would require the use of a path-dependent bulk viscosity. We also find that the effective fluid-particle drag coefficient can be described well with the often-used closure motivated by the experiments of Richardson & Zaki (Trans. Inst. Chem. Engng vol. 32, 1954, p. 35). In this respect, the effect of the system size for determining the drag requires specific care. The shear viscosity of the particle phase manifests small, but clearly noticeable dependence on the rate of compaction/dilation of the particle phase. Our observations point to the need for higher-order closures that recognize the slow evolution of the microstructure in these flows and account for the effects of non-equilibrium microstructure on the stresses.

Original language | English |
---|---|

Pages (from-to) | 303-336 |

Number of pages | 34 |

Journal | Journal of Fluid Mechanics |

Volume | 587 |

DOIs | |

Publication status | Published - 25 Sep 2007 |

### Keywords

- DISCRETIZED BOLTZMANN-EQUATION
- HIGH PARTICLE INERTIA
- PARTICULATE SUSPENSIONS
- REYNOLDS-NUMBER
- MONODISPERSE SPHERES
- SOLID FLOWS
- FORCE
- PRESSURE
- ARRAYS
- MOTION

### Cite this

*Journal of Fluid Mechanics*,

*587*, 303-336. https://doi.org/10.1017/S0022112007007094

**Direct numerical simulations of dense suspensions : wave instabilities in liquid-fluidized beds.** / Derksen, J. J.; Sundaresan, S.

Research output: Contribution to journal › Article

*Journal of Fluid Mechanics*, vol. 587, pp. 303-336. https://doi.org/10.1017/S0022112007007094

}

TY - JOUR

T1 - Direct numerical simulations of dense suspensions

T2 - wave instabilities in liquid-fluidized beds

AU - Derksen, J. J.

AU - Sundaresan, S.

PY - 2007/9/25

Y1 - 2007/9/25

N2 - We present results of direct numerical simulations of travelling waves in dense assemblies of monodisperse spherical particles fluidized by a liquid. The cases we study have been derived from the experimental work of others. In these simulations, the flow of interstitial fluid is solved by the lattice- Boltzmann method (LBM) and the particles move under the influence of gravity, hydrodynamic forces stemming from the LBM, subgrid-scale lubrication forces and hard-sphere collisions. We first show that the propagating inhomogeneous structures seen in the simulations are in agreement with those observed experimentally. We then use the detailed information contained in the simulation results to assess aspects of two-fluid model closures, namely, fluid-particle drag, and the various contributions to the effective stresses. We show that the rates of compaction and dilation of the particle phase in the travelling waves are comparable to the rate at which the microstructure relaxes, and that there is a pronounced effect of the rate of compaction on the average collisional normal stress. Although this effect can be expressed as an effective bulk viscosity term, this approach would require the use of a path-dependent bulk viscosity. We also find that the effective fluid-particle drag coefficient can be described well with the often-used closure motivated by the experiments of Richardson & Zaki (Trans. Inst. Chem. Engng vol. 32, 1954, p. 35). In this respect, the effect of the system size for determining the drag requires specific care. The shear viscosity of the particle phase manifests small, but clearly noticeable dependence on the rate of compaction/dilation of the particle phase. Our observations point to the need for higher-order closures that recognize the slow evolution of the microstructure in these flows and account for the effects of non-equilibrium microstructure on the stresses.

AB - We present results of direct numerical simulations of travelling waves in dense assemblies of monodisperse spherical particles fluidized by a liquid. The cases we study have been derived from the experimental work of others. In these simulations, the flow of interstitial fluid is solved by the lattice- Boltzmann method (LBM) and the particles move under the influence of gravity, hydrodynamic forces stemming from the LBM, subgrid-scale lubrication forces and hard-sphere collisions. We first show that the propagating inhomogeneous structures seen in the simulations are in agreement with those observed experimentally. We then use the detailed information contained in the simulation results to assess aspects of two-fluid model closures, namely, fluid-particle drag, and the various contributions to the effective stresses. We show that the rates of compaction and dilation of the particle phase in the travelling waves are comparable to the rate at which the microstructure relaxes, and that there is a pronounced effect of the rate of compaction on the average collisional normal stress. Although this effect can be expressed as an effective bulk viscosity term, this approach would require the use of a path-dependent bulk viscosity. We also find that the effective fluid-particle drag coefficient can be described well with the often-used closure motivated by the experiments of Richardson & Zaki (Trans. Inst. Chem. Engng vol. 32, 1954, p. 35). In this respect, the effect of the system size for determining the drag requires specific care. The shear viscosity of the particle phase manifests small, but clearly noticeable dependence on the rate of compaction/dilation of the particle phase. Our observations point to the need for higher-order closures that recognize the slow evolution of the microstructure in these flows and account for the effects of non-equilibrium microstructure on the stresses.

KW - DISCRETIZED BOLTZMANN-EQUATION

KW - HIGH PARTICLE INERTIA

KW - PARTICULATE SUSPENSIONS

KW - REYNOLDS-NUMBER

KW - MONODISPERSE SPHERES

KW - SOLID FLOWS

KW - FORCE

KW - PRESSURE

KW - ARRAYS

KW - MOTION

U2 - 10.1017/S0022112007007094

DO - 10.1017/S0022112007007094

M3 - Article

VL - 587

SP - 303

EP - 336

JO - Journal of Fluid Mechanics

JF - Journal of Fluid Mechanics

SN - 0022-1120

ER -