scholarly journals Revisiting the Homogenized Lattice Boltzmann Method with Applications on Particulate Flows

Computation ◽  
2021 ◽  
Vol 9 (2) ◽  
pp. 11
Author(s):  
Robin Trunk ◽  
Timo Weckerle ◽  
Nicolas Hafen ◽  
Gudrun Thäter ◽  
Hermann Nirschl ◽  
...  

The simulation of surface resolved particles is a valuable tool to gain more insights in the behaviour of particulate flows in engineering processes. In this work the homogenized lattice Boltzmann method as one approach for such direct numerical simulations is revisited and validated for different scenarios. Those include a 3D case of a settling sphere for various Reynolds numbers. On the basis of this dynamic case, different algorithms for the calculation of the momentum exchange between fluid and particle are evaluated along with different forcing schemes. The result is an updated version of the method, which is in good agreement with the benchmark values based on simulations and experiments. The method is then applied for the investigation of the tubular pinch effect discovered by Segré and Silberberg and the simulation of hindered settling. For the latter, the computational domain is equipped with periodic boundaries for both fluid and particles. The results are compared to the model by Richardson and Zaki and are found to be in good agreement. As no explicit contact treatment is applied, this leads to the assumption of sufficient momentum transfer between particles via the surrounding fluid. The implementations are based on the open-source C++ lattice Boltzmann library OpenLB.

2014 ◽  
Vol 670-671 ◽  
pp. 747-750
Author(s):  
Zhi Jun Gong ◽  
Jiao Yang ◽  
Wen Fei Wu

For indepth study on flow characteristics for fluid bypass obstacles in micro-channel, the Lattice Boltzmann Method (LBM) was used to simulate fluid flow over two circular cylinders in side-by-side arrangement of a micro-channel. The velocity distribution and recirculation zone length under different Reynolds numbers (Re = 0~100) and different spacing ratio (H/D= 0~2.0) were obtained. The results show that the pattern of flow and the size of recirculation zone in the micro-channel depend on the combined effect of Re and H/D.


Author(s):  
Sonam Tanwar

This chapter develops a meshless formulation of lattice Boltzmann method for simulation of fluid flows within complex and irregular geometries. The meshless feature of proposed technique will improve the accuracy of standard lattice Boltzmann method within complicated fluid domains. Discretization of such domains itself may introduce significant numerical errors into the solution. Specifically, in phase transition or moving boundary problems, discretization of the domain is a time-consuming and complex process. In these problems, at each time step, the computational domain may change its shape and need to be re-meshed accordingly for the purpose of accuracy and stability of the solution. The author proposes to combine lattice Boltzmann method with a Galerkin meshfree technique popularly known as element-free Galerkin method in this chapter to remove the difficulties associated with traditional grid-based methods.


Author(s):  
Ru Yang ◽  
Chin-Sheng Wang

A Lattice Boltzmann method is employed to investigate the flow characteristics and the heat transfer phenomenon between two parallel plates separated by a micro-gap. A nine-velocity model and an internal energy distribution model are used to obtain the mass, momentum and temperature distributions. It is shown that for small Knudsen numbers (Kn), the current results are in good agreement with those obtained from the traditional Navier-Stokes equation with non-slip boundary conditions. As the value of Kn is increased, it is found that the non-slip condition may no longer be valid at the wall boundary and that the flow behavior changes to one of slip-flow. In slip flow regime, the present results is still in good agreement with slip-flow solution by Navier Stokes equations. The non-linear nature of the pressure and friction distribution for micro-channel flow is gieven. Finally, the current investigation presents a prediction of the temperature distribution for micro-channel flow under the imposed conditions of an isothermal boundary.


2005 ◽  
Vol 19 (28n29) ◽  
pp. 1551-1554 ◽  
Author(s):  
XIAOKE KU ◽  
JIANZHONG LIN

Flows over two tandem cylinders are simulated numerically based on the lattice Boltzmann method. The pressure distribution on the cylinders for varying distance between the two cylinders at different Reynolds numbers is depicted. The results show that the minimum pressure on the front cylinder does not occur at the stagnation point because of the existence of the back cylinder. The distance between the point with minimum pressure and the stagnation point becomes large with increasing Re number. The minimum pressure on the back cylinder varies with the distance between the two cylinders. The effective distance of interaction between two cylinders is less than 4d with d being the diameter of the cylinder.


2006 ◽  
Vol 22 (1) ◽  
pp. 35-42 ◽  
Author(s):  
J.-S. Wu ◽  
Y.-L. Shao

AbstractThe flows past a square cylinder in a channel are simulated using the multi-relaxation-time (MRT) model in the parallel lattice Boltzmann BGK method (LBGK). Reynolds numbers of the flow are in the range of 100 ∼ 1,850 with blockage ratio, 1/6, of cylinder height to channel height, in which the single-relaxation-time (SRT) scheme is not able to converge at higher Reynolds numbers. Computed results are compared with those obtained using the SRT scheme where it can converge. In addition, computed Strouhal numbers compare reasonably well with the numerical results of Davis (1984).


2008 ◽  
Vol 19 (11) ◽  
pp. 1683-1694 ◽  
Author(s):  
YONG RAO ◽  
YUSHAN NI ◽  
CHAOFENG LIU

The flow around two square cylinders arranged side by side has been investigated through lattice Boltzmann method under different Reynolds numbers and various space ratios (s = d/D, d is the separation distance between two cylinders, D is the characteristic length) from 1.0, 1.1 to 2.7, including 18 space ratios. It is found that the flip-flop regime occurs at small space ratios and the synchronized regime occurs at large space ratios. Wide and narrow wakes at small spacing are formed and intermittently change behind the cylinders, and the biased flow in the gap is bistable. The frequency of vortex shedding is different in two wakes. The upper frequency is smaller than the lower frequency for small space ratios (s < 1.4), and the time-averaged drag and lift coefficients of cylinders are also different. When the space ratios increase, two distinct vortex streets occur behind the cylinders, and the frequency of vortex shedding is almost equal in two wakes. Also the difference of time-averaged drag and lift coefficients of the cylinders decreases with the increase in space ratios; in this case the flow shows synchronized regime. The transition between flip-flop and synchronized regimes occurs at s = 1.5. When s < 1.5, the flow shows flip-flop regime; otherwise, it shows synchronized regime. When s = 2.0 and 2.5, the curves for the time-averaged drag and lift coefficient with different Reynolds numbers are smooth. When s = 1.5 and 1.8, the curves are also smooth under Re ≤ 140, but that will be fluctuant under Re > 140 because of the nonlinear interaction between the wakes, and the instability of flow becomes stronger with the increase in Reynolds numbers. On the other hand, the vortex shedding type from the cylinder occurs in-phase when s < 2.5 and s = 2.5 for Re < 190, whereas that occurs anti-phase when s = 2.5 for Re ≥190. In addition, the pressure varies a little on the left surfaces and greatly on the right surfaces of both cylinders with the increase in Reynolds number at s = 2.5.


Particuology ◽  
2015 ◽  
Vol 21 ◽  
pp. 90-98 ◽  
Author(s):  
Marzie Babaie Rabiee ◽  
Shahram Talebi ◽  
Omid Abouali ◽  
Ehsan Izadpanah

2017 ◽  
Vol 14 (01) ◽  
pp. 1750002 ◽  
Author(s):  
Yi-Kun Wei ◽  
Xu-Qu Hu

Two-dimensional simulations of channel flow past an array of cylinders are carried out at high Reynolds numbers. Considering the thickness fluctuating effect on the equation of motion, a modified lattice Boltzmann method (LBM) is proposed. Special attention is paid to investigate the thickness fluctuations and vortex shedding mechanisms between 11 cylinders. Results for the velocity and vorticity differences are provided, as well as for the energy density and enstrophy spectra. The numerical results coincide very well with some published experimental data that was obtained by turbulent soap films. The spectra extracted from the velocity and vorticity fields are displayed from simulations, along with the thickness fluctuation spectrum H(k). Our results show that the statistics of thickness fluctuations resemble closely those of a passive scalar in turbulent flows.


2019 ◽  
Vol 142 (1) ◽  
Author(s):  
Majid Hassan Khan ◽  
Atul Sharma ◽  
Amit Agrawal

Abstract This article reports flow behavior around a suspended cube obtained using three-dimensional (3D) lattice Boltzmann method (LBM)-based simulations. The Reynolds number (Re) range covered is from 84 to 770. Four different flow regimes are noted based on the flow structure in this range of Re: steady axisymmetric (84 ≤ Re ≤ 200), steady nonaxisymmetric (215 ≤ Re ≤ 250), unsteady nonaxisymmetric in one plane and axisymmetric in the other plane (276 ≤ Re ≤ 300), and unsteady nonaxisymmetric in streamwise orthogonal planes (339 ≤ Re ≤ 770). Recirculation length and drag coefficient follow inverse trend in the steady flow regime. The unsteady flow regime shows hairpin vortices for Re ≤ 300 and then it becomes structureless. The nature of force coefficients has been examined at various Reynolds numbers. Temporal behavior of force coefficients is presented along with phase dependence of side force coefficients. The drag coefficient decreases with increase in Reynolds number in the steady flow regime and the side force coefficients are in phase. Drag coefficients are compared with established correlations for flow around a cube and a sphere. The side force coefficients are perfectly correlated at Re = 215 and they are anticorrelated at Re = 250. At higher Reynolds numbers, side force coefficients are highly uncorrelated. This work adds to the existing understanding of flow around a cube reported earlier at low and moderate Re and extends it further to unsteady regime at higher Re.


2019 ◽  
Vol 11 (03) ◽  
pp. 1950028 ◽  
Author(s):  
N. M. Sangtani Lakhwani ◽  
F. C. G. A. Nicolleau ◽  
W. Brevis

Lattice Boltzmann Method (LBM) simulations for turbulent flows over fractal and non-fractal obstacles are presented. The wake hydrodynamics are compared and discussed in terms of flow relaxation, Strouhal numbers and wake length for different Reynolds numbers. Three obstacle topologies are studied, Solid (SS), Porous Regular (PR) and Porous Fractal (FR). In particular, we observe that the oscillation present in the case of the solid square can be annihilated or only pushed downstream depending on the topology of the porous obstacle. The LBM is implemented over a range of four Reynolds numbers from 12,352 to 49,410. The suitability of LBM for these high Reynolds number cases is studied. Its results are compared to available experimental data and published literature. Compelling agreements between all three tested obstacles show a significant validation of LBM as a tool to investigate high Reynolds number flows in complex geometries. This is particularly important as the LBM method is much less time consuming than a classical Navier–Stokes equation-based computing method and high Reynolds numbers need to be achieved with enough details (i.e., resolution) to predict for example canopy flows.


Sign in / Sign up

Export Citation Format

Share Document