An efficient Discrete Element Lattice Boltzmann model for simulation of particle-fluid, particle-particle interactions

Abstract

In this study, an efficient Discrete Element Lattice Boltzmann Model (DE-LBM) is introduced to simulate mechanical behaviours of multiphase systems involving particle-fluid and particle-particle interactions. The LBM is based on the Multiple Relaxation Time (MRT-LBM) formalism for the fluid phase and the Discrete Element Method for particle motions. A novel algorithm is developed for detecting the particle contact base on particle overlapping areas computed directly from the grid-based LBM data. This contact algorithm achieves the same accuracy in determining the particle contact as provided by the Hertz contact model but is far more efficient computationally. The DE-LBM coupling approach is also modified to unify the different schemes developed previously. A modified Verlet List method for updating the solid occupation fraction is proposed to further speed up the simulation. The new model is validated by a series of simulations including the single particle settling and well-known 'Drafting, Kissing and Tumbling' (DKT) phenomenon found in suspensions. The settling of a large number (2500) of particles in a still fluid is also simulated with predicted concentration profiles matching well the analytic solution. These applications demonstrate the potential of the present DE-LBM model as a powerful numerical tool for simulating multiphase particulate systems encountered in many engineering and science * Corresponding author

Study on Flow and Heat Transfer Characteristics of Porous Media in Engine Particulate Filters Based on Lattice Boltzmann Method

“…Obviously, e i = −e i+4 . Depending on the probability of movement in directions 1, 2, 3, and 4, the solid particle may remain at the original node or move to the nearest adjacent node X * p = X P + θ 1 e 1 + θ 2 e 2 + θ 3 e 3 + θ 4 e 4 (18) where θ i is a Boolean variable, which equal to 1 with probability p i . As shown in Figure 6, p 1 > 0, p 2 > 0, p 3 = 0, p 4 = 0, so the particle may remain stationary with probability…”

“…For the lattice Boltzmann method of pure particle flow, Eshghinejadfar et al [17] used LBM model to consider the interaction between fluid and particles, and used thermal lattice calculation without considering the interaction between particles. Zhang et al [18] used the lattice Boltzmann method as well as discrete element method to describe the interaction between particulate matter and fluid, particulate matter and particulate matter, in the meantime, they have carried out two-dimensional flow calculation without considering the effect of heat transfer.…”

Intercomparison of boundary schemes in Lattice Boltzmann method for flow simulation in porous media

“…The method of Zhang et al is adopted for the cases of RCP spheres array, where one lattice cell may be occupied by multiple spheres. ϵ is summed up by each sphere's solid volume fraction in this lattice cell and the maximum value is 1 44 …”

A simple collision algorithm for arbitrarily shaped objects in particle‐resolved flow simulation using an immersed boundary method

“…Therefore, several contact detection algorithms have been proposed to reduce the computational cost . Zhang et al proposed an efficient collision algorithm for spheroidal particles for the DEM‐CFD (LBM) flow solver. They reduced the computational cost for collision detection by using a lattice grid for the lattice Boltzmann simulation to detect the contact of the object.…”