## Abstract

Time-dependent flows of a Newtonian fluid through periodic arrays of spheres were simulated using the lattice-Boltzmann scheme. By applying a constant body force per unit mass to the fluid, a steady background fluid flow through the array of stationary spheres was first established. Subsequently, a small-amplitude perturbation to the body force, which varied periodically in time, was added and the long-time behaviour of the unsteady flow fields and the forces on the particles were determined. From the simulations, the pressure and friction (shear) forces acting on the particles were determined for a range of conditions. Results on simple cubic lattices are presented. Computations spanned a range of particle volume fractions (0.1 < φ < 0.4) background flow Reynolds numbers (0.25 ≤ Re_{p} ≤ 60, where Re_{p}= 2au_{∫}/ν) and oscillatory flow Reynolds numbers (0.9 ≤ Re_{ω} ≤ 420 with Re_{ω}= 2a^{2}ω/ν). Here u∫ is the superficial velocity of the fluid through the bed, a is the particle radius, ν is the kinematic viscosity of the fluid, and ω is the oscillation frequency. In the limit of Re_{ω} → 0 the quasi-steady-state drag force was obtained. At low Re_{p} this force approached the steady-state drag force, while its increase with Re_{p} was stronger than the steady-state drag force, similar to that for isolated spheres given by Mei et al. (J. Fluid Mech., vol. 233, 1991, p. 613). Th e unsteady force was decomposed into pressure and friction components. The phase angles of these components in the limit Re_{ω} → ∞ indicate that the virtual mass force contributes to the unsteady pressure force while the history force contributes to the friction force. The remainder of the unsteady friction and pressure forces is attributed to unsteady drag force. The a pparent virtual mass coefficient was found to vary from ∼0.5, which is the well-known limit for isolated spheres in inviscid flows, to ∼ 1.0 at low Re_{ω}. This change is clearly a consequence of viscous effects. The Re_{ω} at which the transition between these limits occurs increases with φ. The history force exhibits a strong decay towards lower values of Re_{ω} in accordance with the results of Mei et al. (1991) for isolated spheres; however, the Re_{ω} value at which this decay sets in increases appreciably with φ. This φ-dependence is associated with the limited separation between the particles available for the Stokes boundary layer. It was found that the unsteady drag coefficient β varies with Re_{ω}. At low Re_{p}, the drag coefficient initially decreases with increasing Re_{ω}, passes through a minimum and then increases strongly. With increasing Re_{ω} the relative contribution of pressure and friction forces to the unsteady drag force changes.

Original language | English (US) |
---|---|

Pages (from-to) | 257-287 |

Number of pages | 31 |

Journal | Journal of Fluid Mechanics |

Volume | 552 |

DOIs | |

State | Published - Apr 2006 |

## All Science Journal Classification (ASJC) codes

- Condensed Matter Physics
- Mechanics of Materials
- Mechanical Engineering
- Applied Mathematics