Calculation of the permeability in porous media using the lattice Boltzmann method
Calculation of the permeability in porous media using the lattice Boltzmann method
复制标题
DOI：
10.1016/j.ijheatfluidflow.2016.05.010
复制
发表时间：
2016-12-01
影响因子：
2.6
通讯作者：
Thevenin, Dominique
中科院分区：
文献类型：
--
作者：
Eshghinejadfard, Amir;Daroczy, Laszlo;Thevenin, Dominique
本文采用格子玻尔兹曼法（LBM）来模拟多孔介质中的三维层流并计算相关的渗透率。这项研究采用了使用消息传递接口技术的内部并行代码。研究了三种不同的流动配置：首先，通过手动指定面心立方体（FCC）中的固体单元；然后，在体心立方体（BCC）中做同样的事情；最后，从一组实验性 2D 计算机断层扫描图像中读取固体单元以获得真实的 3D 几何形状。在所有模拟中，雷诺数都保持远低于 1。发现当前的 LBM 模拟可以对渗透率值进行良好的估计。还研究了所用力方案和单次或多次弛豫时间（SRT、MRT）的影响。尽管每种力方案（Guo-SRT、Guo-MRT 和 Shan-Chen-SRT）可能在某些区域显示出更好的结果，但 SRT 模型对弛豫时间的强烈依赖性表明，力方案、弛豫时间和域分辨率的正确选择是所需精度和计算成本之间的折衷。首先，正如预期的那样，更高的分辨率会带来越来越准确的结果，但需要更多的计算成本和时间。其次，与 SRT 模型相比，MRT 模型表现出较低的粘度依赖性，但速度稍慢。此外，结果对较粗域的弛豫时间值更敏感。此外，较低的弛豫时间需要更多的迭代次数才能达到稳定状态。还观察到，FCC 和 BCC 结构的渗透率在大孔隙率下非常接近。最后一个具有真实几何形状的测试案例表明，LB 模拟可以提供 3D 多孔介质中的流体动力学特性，这对于依赖经典纳维-斯托克斯方程的求解器来说是一个非常具有挑战性的问题。 (C) 2016 Elsevier Inc. 保留所有权利。
In this paper, the lattice Boltzmann method (LBM) is used to simulate three-dimensional laminar flows in porous media and to calculate the associated permeability. An in-house, parallelized code using the message passing interface technique is employed for the study. Three different flow configurations are studied: first, by manually specifying solid cells in a face-centered cube (FCC); then, doing the same in a body-centered cube (BCC); and finally by reading the solid cells for a real 3D geometry from a set of experimental 2D computed tomography images. In all simulations, the Reynolds number is kept well below 1. It was found that the current LBM simulations yield good estimates for the permeability value. The impact of the employed force scheme and single- or multiple-relaxation time (SRT, MRT) was also studied. Although each force scheme (Guo-SRT, Guo-MRT and Shan-Chen-SRT) may show better results in some regions, the strong dependency of SRT models on relaxation time suggests that the proper choice of the force scheme, relaxation time and domain resolution is a compromise between the required accuracy and computational cost. First, higher resolutions lead as expected to increasingly accurate results but requires more computational cost and time. Second, the MRT model shows a lower viscosity dependence in comparison with SRT models but is somewhat slower. Also, the results are more sensitive to the relaxation time value for coarser domains. Furthermore, lower relaxation times necessitate a higher number of iterations to reach the steady state. It was also observed that permeability of both FCC and BCC structures are very close at large porosities. The last test case with the real geometry demonstrates that LB simulations can deliver the hydrodynamic properties in a 3D porous media, a very challenging issue for solvers when relying on classical Navier-Stokes equations. (C) 2016 Elsevier Inc. All rights reserved.