共查询到20条相似文献，搜索用时 625 毫秒
1.
The lattice Boltzmann method is a popular tool for pore-scale simulation of flow. This is likely due to the ease of including complex geometries such as porous media and representing multiphase and multifluid flows. Many advancements, including multiple relaxation times, increased isotropy, and others have improved the accuracy and physical fidelity of the method. Additionally, the lattice Bolzmann method is computationally very efficient, thanks to the explicit nature of the algorithm and relatively large amount of local work. The combination of many algorithmic options and efficiency means that a software framework enabling the usage and comparison of these advancements on computers from laptops to large clusters has much to offer. In this paper, we introduce Taxila LBM, an open-source software framework for lattice Boltzmann simulations. We discuss the design of the framework and lay out the features available, including both methods in the literature and a few new enhancements which generalize methods to complex geometries. We discuss the trade-off of accuracy and performance in various methods, noting how the Taxila LBM makes it easy to perform these comparisons for real problems. And finally, we demonstrate a few common applications in pore-scale simulation, including the characterization of permeability of a Berea sandstone and analysis of multifluid flow in heterogenous micromodels. 相似文献
2.
文章的目的是对格子玻尔兹曼方法进行系统的介绍，格子玻尔兹曼方法(Lattice Boltzmann Method)的出现直接来源于20世纪60年代的元胞自动机(Cellular Automata)思想，而这一方法用于解决流动现象时，又可以追溯到19世纪的分子运动论，求解的是Boltzmann提出的玻尔兹曼输运方程，因此将这一方法称为格子玻尔兹曼方法，之前也被称为格子气自动机(Lattice Gas Automaton)。该方法多用于研究复杂现象，如材料晶体凝聚时的生长过程、城市土地利用的演化等方面。在20世纪70年代由Hardy、Pomeau和Pazzis建立了第一个用于研究流体运动的格子气自动机，此后，这一方法被广泛用来模拟各种流动问题，诸如二相流、孔隙介质中的渗流等，并根据这一方法开发了相应的商业软件PowerFlow。同时，格子玻尔兹曼方法由于其在微观水平描述运动的特点，成为研究湍流的一个很好的数值计算工具，特别是用其进行直接数值模拟(DNS)计算，成为继传统的差分法、有限体积法和谱方法之后的又一有力的手段。而作为大气运动的一个主要现象的大气湍流，比普通湍流更加复杂，在这里着重介绍了大气湍流的特点和应用格子玻尔兹曼方法模拟湍流的发展过程。 相似文献
3.
The multiscale transport mechanism of methane in unconventional reservoirs is dominated by slip and transition flows resulting from the ultra-low permeability of micro/nano-scale pores, which requires consideration of the microscale and rarefaction effects. Traditional continuum-based computational fluid dynamics (CFD) becomes problematic when modeling micro-gaseous flow in these multiscale pore networks because of its disadvantages in the treatment of cases with a complicated boundary. As an alternative, the lattice Boltzmann method (LBM), a special discrete form of the Boltzmann equation, has been widely applied to model the multi-scale and multi-mechanism flows in unconventional reservoirs, considering its mesoscopic nature and advantages in simulating gas flows in complex porous media. Consequently, numerous LBM models and slip boundary schemes have been proposed and reported in the literature. This study investigates the predominately reported LBM models and kinetic boundary schemes. The results of these LBM models systematically compare to existing experimental results, analytical solutions of Navier-Stokes, solutions of the Boltzmann equation, direct simulation of Monte Carlo (DSMC) and information-preservation DSMC (IP_DSMC) results, as well as the numerical results of the linearized Boltzmann equation by the discrete velocity method (DVM). The results point out the challenges and limitations of existing multiple-relaxation-times LBM models in predicting micro-gaseous flow in unconventional reservoirs. 相似文献
4.
Two‐phase fluid distributions in fractured porous media were studied using a single‐component multiphase (SCMP) lattice Boltzmann method (LBM), which was selected among three commonly used numerical approaches through a comparison against the available results of micro X‐ray computed tomography. The influence of the initial configuration and the periodic boundary conditions in the SCMP LBM for the fluid distribution analysis were investigated as well. It was revealed that regular porous media are sensitive to the initial distribution, whereas irregular porous media are insensitive. Moreover, to eliminate the influence of boundaries, the model's buffer size of an SCMP LBM simulation was suggested to be taken as approximately 12.5 times the average particle size. Then, the two‐phase fluid distribution of a porous medium was numerically studied using the SCMP LBM. Both detailed distribution patterns and macroscopic morphology parameters were reasonably well captured. Finally, the two‐phase fluid distributions in a fractured porous media were investigated. The influence of the degree of saturation, fracture length, and fracture width on the fluid distributions and migration was explored. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
5.
A. M. Tartakovsky N. Trask K. Pan B. Jones W. Pan J. R. Williams 《Computational Geosciences》2016,20(4):807-834
Smoothed particle hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and advection-diffusion-reaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-scale multiphase flows and reactive transport in porous and fractured media. 相似文献
6.
This paper presents a computational method able to effectively model both the simultaneous processes typically observed in backward erosion piping, ie, the pipe tip propagation and the conduit cross-section enlargement. The numerical method is based on the novel formulation of a problem of localized erosion along a line propagating in a multidimensional porous medium. In this line, a conduit with evolving transverse size is embedded, which conveys a multiphase flow. The two systems, porous medium and pipe, are bridged by exchange terms of multiphase fluid mass and by a shared fluid pressure field. On the contrary, different fields are considered to describe flows, which are assumed as Darcian in the porous medium and turbulent in the conduit. These two flows drive pipe propagation and enlargement, respectively, as modeled by means of proper erosion kinetic laws. The corresponding numerical formulation is based on the combination between one- and multidimensional finite elements, to model the erosion conduit and the porous medium, respectively. Several simulations are proposed to demonstrate the ability of the proposed approach in reproducing available experimental data of real-scale tests on levees. Our results point out the crucial role played by the combined influence of pipe propagation and enlargement, as well as of three-dimensional (3D) effects. We also assess the mesh independence of the proposed numerical solution, particularly as concerns the calculated pipe propagation history. 相似文献
7.
Robert Eymard Cindy Guichard Raphaele Herbin Roland Masson 《Computational Geosciences》2012,16(4):987-1005
This paper concerns the discretisation on general 3D meshes of multiphase compositional Darcy flows in heterogeneous anisotropic porous media. Extending Coats’ formulation [15] to an arbitrary number of phases, the model accounts for the coupling of the mass balance of each component with the pore volume conservation and the thermodynamical equilibrium and dynamically manages phase appearance and disappearance. The spatial discretisation of the multiphase compositional Darcy flows is based on a generalisation of the Vertex Approximate Gradient scheme, already introduced for single-phase diffusive problems in [24]. It leads to an unconditionally coercive scheme for arbitrary meshes and permeability tensors. The stencil of this vertex-centred scheme typically comprises 27 points on topologically Cartesian meshes, and the number of unknowns on tetrahedral meshes is considerably reduced, compared with the usual cell-centred approaches. The efficiency of our approach is exhibited on several examples, including the nearwell injection of miscible CO2 in a saline aquifer taking into account the vaporisation of H2O in the gas phase as well as the precipitation of salt. 相似文献
8.
When fluid flows in porous media under subsurface conditions, significant deformation can occur. Such deformation is dependent on structural and phase characteristics. In this paper, we investigate the effect of multiphase flow on the deformation of porous media at the pore scale by implementing a strongly coupled partitioned solver discretized with finite volume (FV) technique. Specifically, the role of capillary forces on grain deformation in porous media is investigated. The fluid and solid subdomains are meshed using unstructured independent grids. The model is applied for solving multiphase coupled equations and is capable of capturing pore scale physics during primary drainage by solving the Navier-Stokes equation and advecting fluid indicator function using volume of fluid (VOF) while the fluid is interacting with a nonlinear elastic solid matrix. The convergence of the coupled solver is accelerated by Aitken underrelaxation. We also reproduce geomechanical stress conditions, at the pore scale, by applying uniaxial stress on the solid while simultaneously solving the multiphase fluid-solid interaction problem to investigate the effect of external stress on fluid occupancy, velocity-field distribution, and relative permeability. We observe that the solid matrix exhibits elasto-capillary behavior during the drainage sequence. Relative permeability endpoints are shifted on the basis of the external stress exerted. 相似文献
9.
10.
We describe a new approach for simulation of multiphase flows through heterogeneous porous media, such as oil reservoirs.
The method, which is based on the wavelet transformation of the spatial distribution of the single-phase permeabilities, incorporates
in the upscaled computational grid all the relevant data on the permeability, porosity, and other important properties of
a porous medium at all the length scales. The upscaling method generates a nonuniform computational grid which preserves the resolved structure
of the geological model in the near-well zones as well as in the high-permeability sectors and upscales the rest of the geological
model. As such, the method is a multiscale one that preserves all the important information across all the relevant length
scales. Using a robust front-detection method which eliminates the numerical dispersion by a high-order total variation diminishing
method (suitable for the type of nonuniform upscaled grid that we generate), we obtain highly accurate results with a greatly
reduced computational cost. The speed-up in the computations is up to over three orders of magnitude, depending on the degree
of heterogeneity of the model. To demonstrate the accuracy and efficiency of our methods, five distinct models (including
one with fractures) of heterogeneous porous media are considered, and two-phase flows in the models are studied, with and
without the capillary pressure. 相似文献
11.
Raphaël di Chiara Roupert Gerhard Schäfer Philippe Ackerer Michel Quintard Guy Chavent 《Comptes Rendus Geoscience》2010,342(11):855-863
Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand. 相似文献
12.
The intrinsic permeability of microcracks in porous solids: Analytical models and Lattice Boltzmann simulations 
下载免费PDF全文
下载免费PDF全文 The paper presents a synthesis of analytical modeling and computational simulations of the intrinsic permeability of microcracks, embedded in porous materials taking into account the interaction of the fluid flow in the microcrack with the surrounding porous material. In the first part of the paper, using the DARCY , STOKES , BRINKMAN , and the BEAVERS–JOSEPH approximations, we derive the intrinsic permeability of a plain non‐rough microcrack in terms of the microcrack geometry and the permeability of the porous material surrounding the microcrack. In the second part of the paper, the intrinsic permeability of a microcrack is determined by means of computational simulations using the framework of the lattice Boltzmann method with partial bounceback conditions. The comparison of predictions from the analytical model and the numerical simulations show an excellent agreement. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
13.
格子Boltzmann方法地震波场模拟 总被引：3，自引：0，他引：3
格子Ｂｏｌｔｚｍａｎｎ方法是细胞自动机在某些学科中的具体化和应用。它根据微观运动过程的某些基本特征建立简化的、时间和空间完全离散的动力学模型，这种模型的平行行为符合宏观的微分方程。 相似文献
14.
路基冻胀问题严重影响寒区高速铁路的安全服役,而成冰相变过程是解释冻胀机制的关键。基于介观尺度的格子Boltzmann方法,将修正的孔隙水冻结温度算法与焓法固液相变格子Boltzmann模型相结合,模拟了悬浮液滴冻结和冻土孔隙水成冰两个过程,分别揭示了液态水在自由状态和孔隙束缚状态下冰水相变的细观机制。计算结果表明：土体孔隙中冰晶由中心向外生长的过程与悬浮在空气中的液滴冻结过程截然不同,并且孔隙水越接近颗粒表面,其冻结温度越低。相同粒径颗粒按照不同排列方式得到的冻结特征曲线（soil freezing characteristic curves,简称SFCC）具有明显差异;不同粒径的SFCC随着颗粒增大残余水含量逐渐变少,形态更加陡峭。通过与文献试验结果对比,验证了格子Boltzmann方法的有效性,表明该方法能够为研究多孔介质水气迁移与相变过程提供介观尺度的新手段。 相似文献
15.
Janez Perko 《Computational Geosciences》2018,22(6):1423-1432
In this paper, we describe a single-relaxation-time (SRT) lattice Boltzmann formulation, which can be effectively applied to anisotropic advection-dispersion equations (AADE). The formulation can be applied to space and time variable anisotropic hydrodynamic dispersion tensor. The approach utilizes diffusion velocity lattice Boltzmann formulation which in the case of AADE can represent anisotropic diagonal and off-diagonal elements of the dispersion matrix by the coupling of advective and diffusive fluxes in equilibrium function. With this approach, AADE can be applied to the SRT lattice Boltzmann formulation using the same equilibrium function and without any changes to collision step nor in the application of boundary conditions. The approach shows good stability even for highly anisotropic dispersion tensor and is tested on selected illustrative examples which demonstrate the accuracy and applicability of the proposed method. 相似文献
16.
用格子波耳兹曼方法模拟双重孔隙介质中的流体迁移 总被引：3，自引：0，他引：3
作者在本文中介绍了基于格子波耳兹曼模型的双重孔隙介质中流体运移的数值模拟计算方法。我们从格子波耳兹曼碰撞模型出发，利用格子波耳兹曼方程、Chapman-Enskog展开，以及多尺度技术，得到了描述双重孔隙介质中流体迁移的二维扩散方程。利用格子气自动机方法计算该扩散方程，实现了对双重孔隙介质中流体运移过程的数值模拟仿真。数值实验表明，我们所使用的方法正确、有效。 相似文献
17.
为了研究土体的细观渗流特性,假设土体是完全饱和且在渗流过程中水分的流动始终处于层流状态。考虑宏观统计参数（孔隙率、渗透率及有效黏滞系数等）的影响,基于表征体元（REV）尺度的格子博尔兹曼（Boltzmann）方法,建立了压力作用下土体细观渗流的数值模型。采用D2Q9模型考虑水分流动的离散速度分布,宏观边界条件为左右侧面为不透水边界 ,上下边界设置不同的密度来控制压力边界,在微观边界条件上采用非平衡态外推格式。编制相应的计算程序,将计算区域内的多孔介质材料设置成流体（孔隙率 1.0,渗透率 ）,验证了经典的Poiseuille流。此外,结合算例分别讨论了土体在压力作用下孔隙率、渗透率及渗透压力等影响因素与渗流速度的相互关系,研究表明该数值方法与Darcy定律得到的计算结果较为吻合。因此,基于REV尺度的格子Boltzmann方法可以有效地模拟土体的渗流机制,为进一步研究土体渗流特性提供了一种新的研究手段。 相似文献
18.
This paper describes the application of a three-dimensional lattice Boltzmann method (LBM) to Newtonian and non-Newtonian (Bingham fluid in this work) flows with free surfaces. A mass tracking algorithm was incorporated to capture the free surface, whereas Papanastasiou’s modified model was used for Bingham fluids. The lattice Boltzmann method was first validated using two benchmarks: Newtonian flow through a square cross-section tube and Bingham flow through a circular cross-section tube. Afterward, the dam-break problem for the Newtonian fluid and the slump test for Bingham fluid were simulated to validate the free-surface-capturing algorithm. The numerical results were in good agreement with analytical results, as well as other simulations, thereby proving the validity and correctness of the current method. The proposed method is a promising substitute for time-consuming and costly physical experiments to solve problems encountered in geotechnical and geological engineering, such as the surge and debris flow induced by a landslide or earthquake. 相似文献
19.
The lattice Boltzmann (LB) method is an efficient technique for simulating fluid flow through individual pores of complex porous media. The ease with which the LB method handles complex boundary conditions, combined with the algorithm’s inherent parallelism, makes it an elegant approach to solving flow problems at the sub-continuum scale. However, the realities of current computational resources can limit the size and resolution of these simulations. A major research focus is developing methodologies for upscaling microscale techniques for use in macroscale problems of engineering interest. In this paper, we propose a hybrid, multiscale framework for simulating diffusion through porous media. We use the finite element (FE) method to solve the continuum boundary-value problem at the macroscale. Each finite element is treated as a sub-cell and assigned permeabilities calculated from subcontinuum simulations using the LB method. This framework allows us to efficiently find a macroscale solution while still maintaining information about microscale heterogeneities. As input to these simulations, we use synchrotron-computed 3D microtomographic images of a sandstone, with sample resolution of 3.34 μm. We discuss the predictive ability of these simulations, as well as implementation issues. We also quantify the lower limit of the continuum (Darcy) scale, as well as identify the optimal representative elementary volume for the hybrid LB–FE simulations. 相似文献
20.
We consider the impact of using time-lapse seismic data in addition to production data for permeability estimation in a porous
medium with multiphase fluid flows, such as a petroleum reservoir under water-assisted production. Since modeling seismic
wave propagation in addition to modeling fluid flows in the reservoir is quite involved, it is assumed that the time-lapse
seismic data have already been inverted into fluid saturation differences (pseudoseismic data). Because an inversion process
often leads to considerable error growth, we will consider pseudoseismic data with large uncertainties. The impact of pseudoseismic
data is assessed through permeability estimation with and without such data and through application of some uncertainty measures
for the estimated parameters. A multiscale algorithm is used for the parameter estimations, so that potential differences
in attainable permeability resolution will be easily revealed. The numerical examples clearly indicate that the permeability
estimation problem is stabilized at a higher level of resolution when pseudoseismic data are applied in addition to production
data, even if the pseudoseismic data have large associated uncertainties. Use of the parameter uncertainty measures confirm
these results. 相似文献