New studying of Lattice Boltzmann method for two-phase driven in porous media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｗｏ＿ｐｈａｓｅｄｒｉｖｅｎｉｎｐｏｒｏｕｓｍｅｄｉａｉｓａｎｉｍｐｏｒｔａｎｔｓｕｂｊｅｃｔｏｆｗｉｄｅｉｎｔｅｒｄｉｓｃｉｐｌｉｎａｒｙｃｏｎｃｅｒｎ ,ｓｕｃｈａｓｐｅｔｒｏｌｅｕｍｉｎｄｕｓｔｒｙ ,ｓｕｂｔｅｒｒａｎｅａｎ ,ｃｈｅｍｉｃａｌｐｒｏｊｅｃｔｅｔｃ..Ｉｔ’ｓｖｅｒｙｄｉｆｆｉｃｕｌｔｉｆｗｅａｄｏｐｔｎｕｍｅｒｉｃａｌｉｎｔｅｇｒａｌｍｅｔｈｏｄｆｏｒｔｈｅｃａｓｅｏｆｍｕｌｔｉｐｈａｓｅａｎｄｍｕｌｔｉｃｏｍｐｏｓｉｔｉｏｎ ,ｔｈｅｃｒｕｘｏｆ…     

Intercomparison of boundary schemes in Lattice Boltzmann method for flow simulation in porous media

The Lattice Boltzmann method has been widely adopted to simulate flow in porous media. The choice of appropriate boundary schemes is essential to achieve simulation accuracy; however, the criteria for the most suitable boundary treatment in the simulation of flow in porous media flow remain unresolved. Here, three types of the most commonly used boundary conditions are tested: interpolation bounce back (IBB), partial saturated method (PSM), and immersed boundary method (IBM). The dimensionless drag of face-centered cubic (FCC) sphere array and the dimensionless permeability of a random closely packed (RCP) sphere array are calculated and compared at different viscosities and resolutions. In the FCC sphere array case where spheres are not contacted, the IBB and PSM exhibit the same accuracy and both are of the second-order convergence rate. The IBM is less accurate and is of the first-order convergence rate. In the RCP sphere array case where the spheres are contacted, the IBB shows finer results and a second-order convergence rate. PSM underestimates the dimensionless permeability and increases resolution only slightly improved the results of PSM. The IBM overestimates the dimensionless permeability. These results indicate that among the three methods, the IBB is the most accurate. The PSM has the same accuracy as the IBB when sediments are not contacted; however, it loses its accuracy in the simulation of flow in closely packed porous media. This work could serve as a benchmark for further research in choosing the most appropriate method in the simulation of flow in porous media.     

Lattice Boltzmann method for the fractional sub‐diffusion equation

A lattice Boltzmann model for the fractional sub‐diffusion equation is presented. By using the Chapman–Enskog expansion and the multiscale time expansion, several higher‐order moments of equilibrium distribution functions and a series of partial differential equations in different time scales are obtained. Furthermore, the modified partial differential equation of the fractional sub‐diffusion equation with the second‐order truncation error is obtained. In the numerical simulations, comparisons between numerical results of the lattice Boltzmann models and exact solutions are given. The numerical results agree well with the classical ones. Copyright © 2015 John Wiley & Sons, Ltd.     

Lattice Boltzmann modeling of pore‐scale fluid flow through idealized porous media

In order to capture the hydro‐mechanical impacts on the solid skeleton imposed by the fluid flowing through porous media at the pore‐scale, the flow in the pore space has to be modeled at a resolution finer than the pores, and the no‐slip condition needs to be enforced at the grain–fluid interface. In this paper, the lattice Boltzmann method (LBM), a mesoscopic Navier–Stokes solver, is shown to be an appropriate pore‐scale fluid flow model. The accuracy and lattice sensitivity of LBM as a fluid dynamics solver is demonstrated in the Poiseuille channel flow problem (2‐D) and duct flow problem (3‐D). Well‐studied problems of fluid creeping through idealized 2‐D and 3‐D porous media (J. Fluid Mech. 1959; 5 (2):317–328, J. Fluid Mech. 1982; 115 :13–26, Int. J. Multiphase Flow 1982; 8 (4):343–360, Phys. Fluids A 1989; 1 (1):38–46, Int. J. Numer. Anal. Meth. Geomech. 1999; 23 :881–904, Int. J. Numer. Anal. Meth. Geomech. 2010; DOI: 10.1002/nag.898, Int. J. Multiphase Flow 1982; 8 (3):193–206) are then simulated using LBM to measure the friction coefficient for various pore throats. The simulation results agree well with the data reported in the literature. The lattice sensitivity of the frictional coefficient is also investigated. Copyright © 2010 John Wiley & Sons, Ltd.     

Permeability of microscale fibrous porous media using the lattice Boltzmann method

The permeabilities of microscale fibrous porous media were calculated using the multiple-relaxation-time (MRT) lattice Boltzmann method (LBM). Two models of the microscale fibrous porous media were constructed based on overlapping fibers (simple cubic, body-centered cubic). Arranging the fibers in skew positions yielded two additional models comprising non-overlapping fibers (skewed simple cubic, skewed body-centered cubic). As the fiber diameter increased, the fibers acted as granular inclusions. The effects of the overlapping fibers on the media permeability were investigated. The overlapping fibers yielded permeability values that were a factor of 2.5 larger than those obtained from non-overlapping fibers, but the effects of the fiber arrangement were negligible. Two correlations were obtained for the overlapping and non-overlapping fiber models, respectively. The effects of the rarefaction and slip flow are also discussed. As the Knudsen number increased, the dimensionless permeability increased; however, the increase differed depending on the fiber arrangement. In the slip flow regime, the fiber arrangement inside the porous media became an important factor.     

微通道人工泡状流的Lattice Boltzmann数值模拟

采用高频电控热激发汽泡的方式构造微通道人工泡状流,可以有效抑制微通道沸腾流动的不稳定性和强化传热。本文基于Lattice Boltzmann大密度比多相流复合模型,数值研究了通道内人工泡状流的流动和传热,通过比较分析不同发泡频率的泡状流,量化分析了汽泡运动和增长对微通道流动与传热的相互影响。一方面着重分析了汽泡运动对微通道运动边界层以及汽泡相变增长对热边界层的影响,另一方面也研究了边界层对汽泡动力行为的影响,所得结论对研究抑制微通道沸腾流动不稳定性和强化传热有参考意义。     

Lattice Boltzmann simulation of MHD natural convection in a cavity with porous media and sinusoidal temperature distribution

The lattice Boltzmann method (LBM) is used to simulate the effect of magnetic field on the natural convection in a porous cavity. The sidewalls of the cavity are heated sinusoidally with a phase derivation, whereas the top and bottom walls are thermally insulated. Numerical simulation is performed, and the effects of the pertinent parameters, e.g., the Hartmann number, the porosity, the Darcy number, and the phase deviation, on the fluid flow and heat transfer are investigated. The results show that the heat transfer is affected by the temperature distribution on the sidewalls clearly. When the Hartmann number is 0, the maximum average Nusselt number is obtained at the phase deviation 90°. Moreover, the heat transfer enhances when the Darcy number and porosity increase, while decreases when the Hartman number increases.     

格子Boltzmann方法中基于内插值的曲线边界处理新方法

对格子Boltzmann方法提出了一种新的曲面边界条件处理方法。在笛卡尔坐标系中,这种处理方法是现有的格子Boltzmann方法有关边界条件处理与浸入式边界条件的混合,它采用内插值方法计算靠近物理边界的网格点速度,使其保证最低精度为二阶,然后利用格子Boltzmann方法中的边界条件技术得到相应的分布函数。由理论推导和数值计算表明,本文提出的方法比其他方法更稳定且具有二阶精度。     

Three-dimensional flow driven pore-crack networks in porous composites: Boltzmann Lattice method and hybrid hypersingular integrals

This study introduces a hybrid hypersingular integral equation-Lattice Boltzmann method (HHIE-LBM) for analyzing extended 3D flow driven pore-crack networks problem in various porosity composites. First, the extended hybrid electronic-ionic, thermal, magnetic, electric and force coupled fields’ pressure and velocity boundary conditions for HHIE-LBM model are established, and the closed form solutions of extended distribution functions are given. Second, an extended 3D flow driven pore-crack networks problem in various porosity composites is translated into a coupled of HHIE-LBM equations, and the pore-crack networks propagation parameters are analyzed. Third, the extended dynamic stress intensity factors (EDSIFs) are calculated by using the parallel numerical technology and the visualization results are presented. Last, the relationship between the EDSIFs and the differential porosity is discussed, and several rules have been found, which can be utilized to understand the extended fluid flow mechanism in various porosity composites and analyze the extended fluid flow varying mechanism on coseismal slip.     

Lattice Boltzmann method for polymer kinetic theory

The purpose of this work is to extend the applicability of the lattice Boltzmann method (LBM) to the field of polymer kinetic theory or more generally suspensions that could be described in the Fokker–Planck formalism. This method has been, in a first time, used for gas kinetic theory, where the resolution space corresponds to the physical space coordinate. In a second time is has been generalized to be applied to fluid flow involving different behaviours: turbulence, porous media, multiphase flow, etc. However this powerful, parallel, and efficient algorithm has not been applied for solving Fokker–Planck equations widely used to describe suspension kinetic theory. In this scale, molecular models involve a high computational costs because of the multidimensionality of the fully coupled micro–macro complex flow. The originality of this work consists to apply the lattice Boltzmann technique for solving Fokker–Planck equation based on a discretization of the configuration space where the resolution coordinates correspond to the microscopic configuration space (and not the physical coordinates). The result of this work emphasizes the optimality of the used technique that, in addition to its parallel ability, gathers the simplicity of the stochastic simulation and the robustness of the traditional fixed mesh support (such as the finite element method). Accuracy and convergence of the LBM will be compared to the stochastic and the finite element techniques for homogeneous shear flow.     

A Lattice Boltzmann study of non-newtonian flow in digitally reconstructed porous domains

In the present study, the Lattice Boltzmann Method (LBM) is applied to simulate the flow of non-Newtonian shear-thinning fluids in three-dimensional digitally reconstructed porous domains. The non-Newtonian behavior is embedded in the LBM through a dynamical change of the local relaxation time. The relaxation time is related to the local shear rate in such a way that the power law rheology is recovered. The proposed LBM is applied to the study of power-law fluids in ordered sphere packings and stochastically reconstructed porous domains. A linear relation is found between the logarithm of the average velocity and the logarithm of the body force with a curve slope approximately equal to the inverse power-law index. The validity of the LBM for the flow of shear thinning fluids in porous media is also tested by comparing the average velocity with the well known semi-empirical Christopher–Middleman correlation. Good agreement is observed between the numerical results and the Christopher–Middleman correlation, indicating that the LBM combined with digital reconstruction constitutes a powerful tool for the study of non-Newtonian flow in porous media.     

Lattice Boltzmann simulations for flow and heat/mass transfer problems in a three‐dimensional porous structure

The lattice Boltzmann method (LBM) for a binary miscible fluid mixture is applied to problems of transport phenomena in a three‐dimensional porous structure. Boundary conditions for the particle distribution function of a diffusing component are described in detail. Flow characteristics and concentration profiles of diffusing species at a pore scale in the structure are obtained at various Reynolds numbers. At high Reynolds numbers, the concentration profiles are highly affected by the flow convection and become completely different from those at low Reynolds numbers. The Sherwood numbers are calculated and compared in good agreement with available experimental data. The results indicate that the present method is useful for the investigation of transport phenomena in porous structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

Lattice Boltzmann method for melting/solidification problems

The present work uses the Lattice Boltzmann method for solving solid/liquid phase change problems. The computed results demonstrate a good agreement with the existing benchmark solution for natural convection and with the experimental solution for solid/liquid interface interacting with the flow field. To cite this article: E. Semma et al., C. R. Mecanique 335 (2007).     

The zonal linearization method for multidimensional mass transfer problems in porous media

A number of methods have been proposed in recent years for calculating the combined flows of immiscible and miscible liquids in strata to systems of boreholes. We propose a method which can naturally be called the zonal linearization method [1]. It is more compact than the usual finite-difference method and has high accuracy, in particular, in the neighborhood of a borehole, since it is closely similar to the method of characteristics. The method can be applied to both continuous and discontinuous flows and in principle makes it possible to investigate the formation and breakdown of discontinuities. As distinct from the method of characteristics, it is well suited to programming and implementation on a computer, and it also makes it possible to obtain an approximate analytic solution of the problem in many cases and to estimate the accuracy of the solution. The method is based on the zonal linearization of the equation for mass conservation in the total flow between chosen surfaces or contour lines (lines of equal saturation or concentration). Determination of the dynamics of the contour surfaces leads to a Cauchy problem for a system of integrodifferential equations involving partial derivatives. The zonal linearization method is a development of the scheme described in [2–4], and the method of solving the Cauchy problem is a generalization of the methods described in [4–13]. The essence of the method and its convergence are illustrated by means of two-dimensional problems in two-phase filtration.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 66–80, July–August, 1973.     

A ghost fluid Lattice Boltzmann method for complex geometries

A ghost fluid Lattice Boltzmann method (GF‐LBM) is developed in this study to represent complex boundaries in Lattice Boltzmann simulations of fluid flows. Velocity and density values at the ghost points are extrapolated from the fluid interior and domain boundary via obtaining image points along the boundary normal inside the fluid domain. A general bilinear interpolation algorithm is used to obtain values at image points which are then extrapolated to ghost nodes thus satisfying hydrodynamic boundary conditions. The method ensures no‐penetration and no‐slip conditions at the boundaries. Equilibrium distribution functions at the ghost points are computed using the extrapolated values of the hydrodynamic variables, while non‐equilibrium distribution functions are extrapolated from the interior nodes. The method developed is general, and is capable of prescribing Dirichlet as well as Neumann boundary conditions for pressure and velocity. Consistency and second‐order accuracy of the method are established by running three test problems including cylindrical Couette flow, flow between eccentric rotating cylinders and flow over a cylinder in a confined channel. Copyright © 2011 John Wiley & Sons, Ltd.     

格子玻尔兹曼法多块网格的交界结构优化研究

基于格子玻尔兹曼方法LBM(Lattice Boltzmann Method)对多块网格方法(Multi-Block)的粗细网格交界结构进行了研究,提出了一种新的优化处理方案。解决了原有网格交界结构存在的三个问题,即两套插值运算造成的程序结构复杂的问题,存储前几个时间步的节点流场数据以备插值运算造成内存浪费的问题和基于时间插值结果进行空间插值计算造成插值误差积累的问题。用一次多点二维空间插值的方式,将原方法的空间和时间双插值,简并成一次空间插值。通过对经典的非定常的圆柱绕流算例和定常的标准顶盖方腔驱动流算例的仿真模拟,验证了交界面处质量、动量及应力的连续性以及网格交界面数据过渡的流畅度,最终验证了改进方法的正确性。数值模拟结果表明,改进后多块算法可实现局部网格细化,进一步推动LBM方法在实际工程问题中的应用。     

Impact of local diffusion on macroscopic dispersion in three-dimensional porous media

While macroscopic longitudinal and transverse dispersion in three-dimensional porous media has been simulated previously mostly under purely advective conditions, the impact of diffusion on macroscopic dispersion in 3D remains an open question. Furthermore, both in 2D and 3D, recurring difficulties have been encountered due to computer limitation or analytical approximation. In this work, we use the Lagrangian velocity covariance function and the temporal derivative of second-order moments to study the influence of diffusion on dispersion in highly heterogeneous 2D and 3D porous media. The first approach characterizes the correlation between the values of Eulerian velocity components sampled by particles undergoing diffusion at two times. The second approach allows the estimation of dispersion coefficients and the analysis of their behaviours as functions of diffusion. These two approaches allowed us to reach new results. The influence of diffusion on dispersion seems to be globally similar between highly heterogeneous 2D and 3D porous media. Diffusion induces a decrease in the dispersion in the direction parallel to the flow direction and an increase in the dispersion in the direction perpendicular to the flow direction. However, the amplification of these two effects with the permeability variance is clearly different between 2D and 3D. For the direction parallel to the flow direction, the amplification is more important in 3D than in 2D. It is reversed in the direction perpendicular to the flow direction.