An immersed boundary/multi-relaxation time lattice Boltzmann method on adaptive octree grids for the particle-resolved simulation of particle-laden flows

• Published in: Journal of computational physics Vol. 471; p. 111669
• Main Authors: Cheng, Zihao, Wachs, Anthony
• Format: Journal Article
• Language: English
• Published: Elsevier Inc 15.12.2022
• Subjects: Adaptive mesh refinement; Fluid-solid interaction; Immersed boundary method; Lattice Boltzmann method; Parallel computing; Particle-laden flow; Fluid-solid interaction; Parallel computing; Adaptive mesh refinement; Lattice Boltzmann method; Immersed boundary method; Particle-laden flow
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2022.111669

We present an immersed boundary/multi-relaxation time lattice Boltzmann method for the particle-resolved simulation of particle-laden flows. The no-slip boundary condition is handled by an explicit feedback immersed boundary method for fixed and moving particles with arbitrary shape. Two special treatments are applied to improve the simulation stability and accuracy: (i) a smoothed discrete delta function (Yang et al. (2009) [47]) is adopted to suppress the spurious force oscillations and (ii) a multi-relaxation time collision operator is utilized to fix the viscosity-dependent numerical slip error of the immersed boundaries. To further reduce the computational effort, we extend the method from fixed and uniform Cartesian grids to adaptive quadtree/octree grids by implementing it in the open-source software Basilisk. A Lax-Wendroff streaming operator is adopted in our method to retain second-order accuracy in both space and time on non-uniform grids, as well as to maintain a unified time scale over the entire computational domain. Consequently, one collision-streaming operation only is required per time step at all grid levels. We test our proposed method on a set of validation cases corresponding to assorted incompressible flows with fixed and moving rigid particles in both 2D and 3D. The accuracy and robustness of our method are demonstrated via thorough comparison with the analytical, experimental and numerical data provided in the literature. •Development of an IB/LBM to simulate 2D and 3D flows laden with fixed/moving particles.•Extension of the method to adaptive mesh refinement on quad/octree grids.•Implementation in the open-source software Basilisk.•Satisfactory parallel computing performance for large-scale computations.•Comprehensive validations to verify the robustness and accuracy of the method.