Eulerian–Lagrangian and Eulerian–Eulerian approaches for the simulation of particle-laden free surface flows using the lattice Boltzmann method

Heidler, Václav; Bublík, Ondřej; Pecka, Aleš; Vimmr, Jan

doi:10.1016/j.cam.2021.113672

Cited by 2 publications

(1 citation statement)

References 37 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The lattice Boltzmann method (LBM) [1][2][3] applies for fluid modeling within a wide range of engineering, biological and physical problems with complex static and moving surfaces, such as particle-laden ones [4][5][6][7][8][9], suspensions of soft particles [10], red blood cells [11] and pulsatile [12,13] flows, porous flow in materials [14], synthetic structures [15][16][17][18][19][20], or natural rocks [21][22][23][24]. These problems are essentially described by the Stokes and finite Reynolds number regimes, and characterized by a coarse grid resolution over a narrow fluid path.…”

Section: Introductionmentioning

confidence: 99%

Unified directional parabolic-accurate lattice Boltzmann boundary schemes for grid-rotated narrow gaps and curved walls in creeping and inertial fluid flows

Ginzburg

Silva²,

Marson³

et al. 2023

Phys. Rev. E

View full text Add to dashboard Cite

The goal of this work is to advance the characteristics of existing lattice Boltzmann Dirichlet velocity boundary schemes in terms of the accuracy, locality, stability, and mass conservation for arbitrarily grid-inclined straight walls, curved surfaces, and narrow fluid gaps, for both creeping and inertial flow regimes. We reach this objective with two infinite-member boundary classes: (1) the single-node "Linear Plus" (LI + ) and (2) the two-node "Extended Multireflection" (EMR). The LI + unifies all directional rules relying on the linear combinations of up to three pre-or postcollision populations, including their "ghost-node" interpolations and adjustable nonequilibrium approximations. On this basis, we propose three groups of LI + nonequilibrium local corrections:(1) the LI + 1 is parametrized, meaning that its steady-state solution is physically consistent: the momentum accuracy is viscosity-independent in Stokes flow, and it is fixed by the Reynolds number (Re) in inertial flow;(2) the LI + 3 is parametrized, exact for arbitrary grid-rotated Poiseuille force-driven Stokes flow and thus most accurate in porous flow; and (3) the LI + 4 is parametrized, exact for pressure and inertial term gradients, and hence advantageous in very narrow porous gaps and at higher Reynolds range. The directional, two-relaxation-time collision operator plays a crucial role for all these features, but also for efficiency and robustness of the boundary schemes due to a proposed nonequilibrium linear stability criterion which reliably delineates their suitable coefficients and relaxation space. Our methodology allows one to improve any directional rule for Stokes or Navier-Stokes accuracy, but their parametrization is not guaranteed. In this context, the parametrized two-node EMR class enlarges the single-node schemes to match exactness in a grid-rotated linear Couette flow modeled with an equilibrium distribution designed for the Navier-Stokes equation (NSE). However, exactness of a grid-rotated Poiseuille NSE flow requires us to perform (1) the modification of the standard NSE term for exact bulk solvability and (2) the EMR extension towards the third neighbor node. A unique relaxation and equilibrium exact configuration for grid-rotated Poiseuille NSE flow allows us to classify the Galilean invariance characteristics of the boundary schemes without any bulk interference; in turn, its truncated solution suggests how, when increasing the Reynolds number, to avoid a deterioration of the mass-leakage rate and momentum accuracy due to a specific Reynolds scaling of the kinetic relaxation collision rate. The optimal schemes and strategies for creeping and inertial regimes are then singled out through a series of numerical tests, such as grid-rotated channels and rotated Couette flow with wall-normal injection, cylindrical porous array, and Couette flow between concentric cylinders, also comparing them against circular-shape fitted FEM solutions.

show abstract