Lattice Boltzmann methods in porous media simulations: From laminar to turbulent flow
Lattice Boltzmann methods in porous media simulations: From laminar to turbulent flow
复制标题
DOI：
10.1016/j.compfluid.2016.10.007
复制
发表时间：
2016-11-25
影响因子：
2.8
通讯作者：
Helmig, Rainer
中科院分区：
文献类型：
--
作者：
Fattahi, Ehsan;Waluga, Christian;Helmig, Rainer
格子Boltzmann方法已成为一种流行的工具，用于确定在很宽的雷诺数和固体体积分数的多孔介质中的阻力和渗透率的相关性。为了获得准确和相关的结果，不仅要实现非常有效的代码，而且要选择最合适的仿真设置。此外，它是必要的，以准确地评估的边界条件和碰撞模型，是有效的从斯托克斯制度的惯性和湍流制度。在本文中，我们比较各种无滑移边界计划和碰撞算子，以评估其效率和精度。而不是假设一个恒定的体积力驱动的流量，一个周期性的压降边界条件来模仿压力驱动的流量通过简单的球包在一个周期性的域。我们首先考虑斯托克斯流场中具有不同碰撞算子的各种边界条件的收敛率。此外，我们选择不同的边界条件，代表第一至第三阶计划在弯曲的边界，以评估其收敛速度数值惯性和湍流。我们发现，多次反射边界条件产生二阶惯性流，而它收敛于三阶斯托克斯制度。考虑到计算量和精度要求，我们选择中心线性插值反弹格式结合双弛豫碰撞模型。该组合的特征在于提供粘度无关的结果和二阶空间收敛。该方法被应用于执行模拟触摸球排列在一个简单的立方体阵列。全模板和简化模板网格模型，即，对D-3（Q（27））和D-3（Q（19））分别进行了比较，并给出了雷诺数在0.001到2，477范围内的阻力和摩擦系数结果。使用这两种不同的网格模型计算的阻力在本研究中考虑的最高雷诺数下的相对差异低于3%。使用的评估结果，我们证明了模型和软件在两个大规模的计算，第一个通过非结构化填料的球形颗粒的流量，第二个可渗透床的湍流的灵活性。(C)2016爱思唯尔有限公司版权所有
The lattice Boltzmann method has become a popular tool for determining the correlations for drag force and permeability in porous media for a wide range of Reynolds numbers and solid volume fractions. In order to achieve accurate and relevant results, it is important not only to implement very efficient code, but also to choose the most appropriate simulation setup. Moreover, it is essential to accurately evaluate the boundary conditions and collision models that are effective from the Stokes regime to the inertial and turbulent flow regimes. In this paper, we compare various no-slip boundary schemes and collision operators to assess their efficiency and accuracy. Instead of assuming a constant volume force driving the flow, a periodic pressure drop boundary condition is employed to mimic the pressure-driven flow through the simple sphere pack in a periodic domain. We first consider the convergence rates of various boundary conditions with different collision operators in the Stokes flow regime. Additionally, we choose different boundary conditions that are representatives of first- to third-order schemes at curved boundaries in order to evaluate their convergence rates numerically for both inertial and turbulent flow. We find that the multi-reflection boundary condition yields second order for inertial flow while it converges with third order in the Stokes regime. Taking into account the both computational cost and accuracy requirements, we choose the central linear interpolation bounce-back scheme in combination with the two-relaxation-time collision model. This combination is characterized by providing viscosity independent results and second order spatial convergence. This method is applied to perform simulations of touching spheres arranged in a simple cubic array. Full- and reduced-stencil lattice models, i.e., the D-3(Q(27) and D-3(Q(19), respectively, are compared and the drag force and friction factor results are presented for Reynolds numbers in the range of 0.001 to 2,477. The drag forces computed using these two different lattice models have a relative difference below 3% for the highest Reynolds number considered in this study. Using the evaluation results, we demonstrate the flexibility of the models and software in two large scale computations, first a flow through an unstructured packing of spherical particles, and second for the turbulent flow over a permeable bed. (C) 2016 Elsevier Ltd. All rights reserved.