Analysis of Fluid Flow in Porous Media Using the Lattice Boltzmann Method: Inertial Flow Regime

Author(s):  
Huizhe Zhao ◽  
Aydin Nabovati ◽  
Cristina H. Amon

In this work, we use the lattice Boltzmann method to study inertial flow in three-dimensional random fibrous porous materials. In order to validate the methodology, inertial flow in two-dimensional hexagonal arrangements of circular cylinders is simulated, and the results are compared against those previously reported in the literature. The three-dimensional fibrous porous materials are then constructed by randomly placing straight cylindrical fibers inside the computational domain. Inertial effects are studied systematically for a wide range of pore Reynolds numbers in materials with porosities between 0.60 and 0.95. A previously proposed semi-empirical relation is modified to represent the inertial effects in three-dimensional fibrous materials. Three distinct regimes of constant, quadratic, and linear relations between the inverse of the permeability and pore Reynolds number are observed for both two- and three-dimensional simulations. The critical Reynolds number, beyond which the inertial effects are strong and this relation is linear, is shown to be smaller in three-dimensional simulations, when compared to the critical Reynolds number in two-dimensional simulations.

Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three-dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio of in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations fail to capture completely the total expansion effect on the flow, which couples both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When the two-dimensional geometry was consistent with the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations was compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, helping to sustain laminar flow symmetry to higher Reynolds numbers in comparison to nominally two-dimensional double-expansion geometries. Lastly, and most importantly, when the logarithm of the critical Reynolds number from all these studies was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


1999 ◽  
Vol 385 ◽  
pp. 41-62 ◽  
Author(s):  
DEWEI QI

A lattice-Boltzmann method has been developed to simulate suspensions of both spherical and non-spherical particles in finite-Reynolds-number flows. The results for sedimentation of a single elliptical particle are shown to be in excellent agreement with the results of Huang, Hu & Joseph (1998) who used a finite-element method. Sedimentation of two-dimensional circular and rectangular particles in a two-dimensional channel and three-dimensional spherical particles in a tube with square cross-section is simulated. Computational results are consistent with experimentally observed phenomena, such as drafting, kissing and tumbling.


Author(s):  
Masoud A. Al Rmah ◽  
Abdulmajeed A. Mohamad

A laminar two-dimensional thermal-flow generated by multiple confined jets impinging on an isothermal plate is investigated numerically using the lattice Boltzmann method. The impinging plate is kept at a constant high temperature while a cold air is flowing through the jets. The effect of different parameters (Reynolds number (Re), and the ratio between the jets height (H) and jets width (W))on the hydrodynamics and thermal characteristics of the flow field is discussed. The Reynolds number is ranging from 50 to 400, and the (H/W) ratio is varying between 1 and 3W.


2013 ◽  
Vol 444-445 ◽  
pp. 395-399
Author(s):  
Di Bo Dong ◽  
Sheng Jun Shi ◽  
Zhen Xiu Hou ◽  
Wei Shan Chen

A lattice Boltzmann method (LBM) with single-relaxation time and on-site boundary condition is used for the simulation of viscous flow in a three-dimensional (3D) lid-driven cavity. Firstly, this algorithm is validated by compared with the benchmark experiments for a standard cavity, and then the results of a cubic cavity with different inflow angles are presented. Steady results presented are for the inflow angle of and, and the Reynolds number is selected as 500. It is found that for viscous flow under moderate Reynolds number, there exists a primary vortex near the center and a secondly vortex at the lower right corner on each slice when, namely in a standard 3D lid-driven cavity, which cant be found when. So it can be thought that the flow pattern in a 3D lid-driven cavity depends not only on the Reynolds number but also the inflow angle.


2003 ◽  
Vol 14 (07) ◽  
pp. 925-944 ◽  
Author(s):  
C. SHU ◽  
X. D. NIU ◽  
Y. T. CHEW

The two-dimensional form of the Taylor series expansion- and least square-based lattice Boltzmann method (TLLBM) was recently presented by Shu et al.8 TLLBM is based on the standard lattice Boltzmann method (LBM), Taylor series expansion and the least square optimization. The final formulation is an algebraic form and essentially has no limitation on the mesh structure and lattice model. In this paper, TLLBM is extended to the three-dimensional case. The resultant form keeps the same features as the two-dimensional one. The present form is validated by its application to simulate the three-dimensional lid-driven cavity flow at Re=100, 400 and 1000. Very good agreement was achieved between the present results and those of Navier–Stokes solvers.


2021 ◽  
Vol 8 (1) ◽  
Author(s):  
Seika Tanji ◽  
Masaru Inatsu ◽  
Tsubasa Okaze

AbstractWe developed a snowdrift model to evaluate the snowdrift height around snow fences, which are often installed along roads in snowy, windy locations. The model consisted of the conventional computational fluid dynamics solver that used the lattice Boltzmann method and a module for calculating the snow particles’ motion and accumulation. The calculation domain was a half channel with a flat free-slip boundary on the top and a non-slip boundary on the bottom, and an inflow with artificially generated turbulence from one side to the outlet side was imposed. In addition to the reference experiment with no fence, experiments were set up with a two-dimensional and a three-dimensional fence normal to the dominant wind direction in the channel center. The estimated wind flow over the two-dimensional fence was characterized by a swirling eddy in the cross section, whereas the wind flow in the three-dimensional fence experiment was horizontally diffluent with a dipole vortex pair on the leeward side of the fence. Almost all the snowdrift formed on the windward side of the two-dimensional and three-dimensional fences, although the snowdrift also formed along the split streaks on the leeward side of the three-dimensional fence. Our results suggested that the fence should be as long as possible to avoid snowdrifts on roads.


2021 ◽  
Author(s):  
Seika Tanji ◽  
Masaru Inatsu ◽  
Tsubasa Okaze

Abstract This study developed a new snowdrift model to evaluate the snowdrift height around a snow fence, often installed along a road in a snowy and windy environment. The model consisted of the conventional computational fluid dynamics (CFD) solver by the Lattice Boltzmann method (LBM) and a module for snow particles’ motion and accumulation. The calculation domain was a half channel with a flat free-slip boundary on the top and a non-slip boundary on the bottom, imposing an inflow with artificially generated turbulence from one side to the other outlet side. Besides the reference experiment with no fence, the experiment was set up with a two-dimensional and a three-dimensional fences normal to the dominant wind direction in the channel center. The estimated wind flow over the two-dimensional fence was characterized by a swirling eddy in the cross-section, whereas the wind flow in the three-dimensional fence experiment was horizontally diffluent with a dipole vortex pair in the leeward of the fence. As a result, almost all of snowdrift was formed in the windward of the two-dimensional and three-dimensional fences, but it was also formed as the split streak in the leeward of the three-dimensional fence. The result suggested that the fence should be as long as possible to avoid the snowdrift on roads.


Author(s):  
Pitambar Randive ◽  
Amaresh Dalal ◽  
Partha P. Mukherjee

The displacement of a three-dimensional immiscible blob subject to oscillatory acoustic excitation in a channel is studied with the Lattice Boltzmann method. The effects of amplitude of the force, viscosity and frequency on blob dynamics are investigated. The trend for variation of mean displacement of blob and frequency response is in agreement to that of the previous two-dimensional studies reported in literature. The response of the blob with pinned contact line shows underdamped behavior. It is also found that increasing the amplitude of the force increases the mean displacement and frequency response.


2005 ◽  
Vol 128 (4) ◽  
pp. 671-679 ◽  
Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations failed to capture completely the total expansion effect on the flow, which is influenced by both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When two-dimensional simulations were performed using the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations were compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, thus sustaining laminar flow symmetry to higher Reynolds numbers. Last, and most important, when the logarithm of the critical Reynolds number was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


Sign in / Sign up

Export Citation Format

Share Document