We investigate the average drag, lift, and torque on static assemblies of capsule-like particles of aspect ratio 4. The performed simulations are from Stokes flow to high Reynolds numbers (0.1 ≤ Re ≤ 1,000) at different solids volume fraction (0.1 ≤ ɛ_s ≤ 0.5). Individual particle forces as a function of the incident angle ϕ with respect to the average flow are scattered. However, the average particle force as a function of ϕ is found to be independent of mutual particle orientations for all but the highest volume fractions. On average, a sine-squared scaling of drag and sine-cosine scaling of lift holds for static multiparticle systems of elongated particles. For a packed bed, our findings can be utilized to compute the pressure drop with knowledge of the particle-orientation distribution, and the average particle drag at ϕ = 0° and 90°. We propose closures for average forces to be used in Euler–Lagrange simulations of particles of aspect ratio 4.

In this paper, we present a number of key numerical methods that can be used to study elongated particles in fluid flows, with a specific emphasis on fluidised beds. Fluidised beds are frequently used for the production of biofuels, bioenergy, and other products from biomass particles, which often have an approximate elongated shape. This raises numerous issues in a numerical approach such as particle-particle contact detection and the accurate description of the various hydrodynamic forces, such as drag, lift, and torque, that elongated particles experience when moving in a fluid flow. The modelling is further complicated by a separation of length scales where industrial flow structures that can extend for many metres evolve subject to solid-solid and solid-fluid interactions at the millimetre scale. As a result, it is impossible to simulate both length scales using the same numerical approach, and a multiscale approach is necessary. First, we outline the direct numerical simulation (DNS) approach that may be employed to estimate hydrodynamic force closures for elongated particles in a fluid flow. We then describe the key aspects of a CFD-DEM approach, which can be used to simulate laboratory scale fluidisation processes, that must be addressed to study elongated particles. Finally, we briefly consider how current industrial-scale models, which concretely assume particle sphericity, could be adapted for the simulation of large collections of elongated particles subject to fluidisation.

Various curved no-slip boundary conditions available in literature improve the accuracy of lattice Boltzmann simulations compared to the traditional staircase approximation of curved geometries. Usually, the required unknown distribution functions emerging from the solid nodes are computed based on the known distribution functions using interpolation or extrapolation schemes. On using such curved boundary schemes, there will be mass loss or gain at each time step during the simulations, especially apparent at high Reynolds numbers, which is called mass leakage. Such an issue becomes severe in periodic flows, where the mass leakage accumulation would affect the computed flow fields over time. In this paper, we examine mass leakage of the most well-known curved boundary treatments for high-Reynolds-number flows. Apart from the existing schemes, we also test different forced mass conservation schemes and a constant density scheme. The capability of each scheme is investigated and, finally, recommendations for choosing a proper boundary condition scheme are given for stable and accurate simulations.

Accurate direct numerical simulations are performed to determine the drag, lift and torque coefficients of non-spherical particles. The numerical simulations are performed using the lattice Boltzmann method with multi-relaxation time. The motivation for this work is the need for accurate drag, lift and torque correlations for high Re regimes, which are encountered in Euler-Lagrangian simulations of fluidization and pneumatic conveying of larger non-spherical particles. The simulations are performed in the Reynolds number range 0.1 ≤ Re ≤ 2000 for different incident angles ϕ. Different tests are performed to analyse the influence of grid resolution and confinement effects for different Re. The measured drag, lift and torque coefficients are utilized to derive accurate correlations for specific non-spherical particle shapes, which can be used in unresolved simulations. The functional forms for the correlations are chosen to agree with the expected physics at Stokes flow as well as the observed leveling off of the drag coefficient at high Re flows. Therefore the fits can be extended to regimes outside the Re regimes simulated. We observe sine-squared scaling of the drag coefficient for the particles tested even at Re=2000 with C_D,ϕ=C_D,ϕ=0^∘ +(C_D,ϕ=90^∘ −C_D,ϕ=0^∘ )sin²ϕ. Furthermore, we also observe that the lift coefficient approximately scales as C_L,ϕ=(C_D,ϕ=90^∘ −C_D,ϕ=0^∘ )sinϕcosϕ for the elongated particles. The current work would greatly improve the accuracy of Euler-Lagrangian simulations of larger non-spherical particles considering the existing literature is mainly limited to steady flow regimes and lower Re.