Improved coupling of time integration and hydrodynamic interaction in particle suspensions using the lattice Boltzmann and discrete element methods. (1st April 2018)

• Record Type: Journal Article
• Authors: Wang, Duo; Leonardi, Christopher R.; Aminossadati, Saiied M.
• Is Part Of: Computers & mathematics with applications. Volume 75:issue 7(2018)
• Page Start: 2593
• Page End: 2606
• Publication Date: 2018-04-01
• Subjects: Lattice Boltzmann method -- Discrete element method -- Partially saturated boundary -- Two-relaxation-time -- Fluid–particle interactions; Electronic data processing -- Periodicals; Mathematics -- Data processing -- Periodicals; 510.28541
• DOI: 10.1016/j.camwa.2018.01.002
• Languages: English
• ISSNs: 0898-1221

Abstract

Abstract: This paper introduces improvements to the simulation of particle suspensions using the lattice Boltzmann method (LBM) and the discrete element method (DEM). First, the benefit of using a two-relaxation-time (TRT) collision operator, instead of the popular Bhatnagar–Gross–Krook (BGK) collision operator, is demonstrated. Second, a modified solid weighting function for the partially saturated method (PSM) for fluid–solid interaction is defined and tested. Results are presented for a range of flow configurations, including sphere packs, duct flows, and settling spheres, with good accuracy and convergence observed. Past research has shown that the drag, and consequently permeability, predictions of the LBM exhibit viscosity-dependence when used with certain boundary conditions such as bounce-back or interpolated bounce-back, and this is most pronounced when the BGK collision operator is employed. The improvements presented here result in a range of computational viscosities, and therefore relaxation parameters, within which drag and permeability predictions remain invariant. This allows for greater flexibility in using the relaxation parameter to adjust the LBM timestep, which can subsequently improve synchronisation with the time integration of the DEM. This has significant implications for the simulation of large-scale suspension phenomena, where the limits of computational hardware persistently constrain the resolution of the LBM lattice.