Intercomparison of 3D pore-scale flow and solute transport simulation methods
Intercomparison of 3D pore-scale flow and solute transport simulation methods
复制标题
DOI：
10.1016/j.advwatres.2015.09.015
复制
发表时间：
2016-09
影响因子：
4.7
通讯作者：
Xiaofan Yang;Y. Mehmani;W. Perkins;A. Pasquali;M. Schönherr;Kyungjoo Kim;M. Perego;M. Parks;Nathaniel Trask;M. Balhoff;M. Richmond;M. Geier;M. Krafczyk;Li-Shi Luo;A. Tartakovsky;Timothy Scheibe
中科院分区：
文献类型：
--
作者：
Xiaofan Yang;Y. Mehmani;W. Perkins;A. Pasquali;M. Schönherr;Kyungjoo Kim;M. Perego;M. Parks;Nathaniel Trask;M. Balhoff;M. Richmond;M. Geier;M. Krafczyk;Li-Shi Luo;A. Tartakovsky;Timothy Scheibe
多种数值方法已被开发出来模拟多孔介质的流体流动和溶质在孔隙尺度上的运移。这些方法包括：1)明确模拟孔隙空间三维几何形状的方法；2)将孔隙空间概念化为一组拓扑一致的程式化孔隙体和孔喉的方法。在之前的工作中，我们使用计算流体动力学（CFD）代码，采用标准有限体积法（FVM），对磁共振测速（MRV）的孔隙尺度速度测量验证了第一种模型。在这里，我们扩展了该验证，包括基于晶格玻尔兹曼方法（LBM）和光滑粒子流体动力学（SPH）的第一种附加模型，以及第二种模型，即孔隙网络模型（PNM）。本研究中使用的PNM方法最近得到了改进，并证明可以在二维实验中准确地模拟溶质输运。虽然PNM方法在计算上比直接数值模拟方法要求低得多，但以PNM的方式概念化复杂的三维孔隙几何形状对溶质运移的影响尚未完全确定。我们应用了所有四种方法（基于fvm的CFD、LBM、SPH和PNM）来模拟孔隙尺度的速度分布和（对于有能力的代码）非反应性溶质输运，并比较了模型结果。从宏观变量（如渗透率、溶质突破曲线）和微观变量（如局部速度和浓度）两方面进行比较。一般来说，各种方法之间达成了良好的一致，但是根据模型上下文观察到一些差异。由于代码的不同功能，相互比较的工作是具有挑战性的，并且激发了一些代码增强，以允许跨全套方法对流和传输模拟进行一致的比较。该研究为多种孔隙尺度模拟方法的可信度提供了支撑，并激励了孔隙尺度模拟方法的进一步发展和应用。
Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.