Paéz-García Catherine Teresa, Valdés-Parada Francisco J, Lasseux Didier
Departamento de I.P.H., Universidad Autónoma Metropolitana-Iztapalapa. Av. San Rafael Atlixco 186, Col. Vicentina, CDMX, 09340, Mexico.
CNRS, UMR 5295, Univ. Bordeaux, Esplanade des Arts et Métiers, 33405 Talence, Cedex, France.
Phys Rev E. 2017 Feb;95(2-1):023101. doi: 10.1103/PhysRevE.95.023101. Epub 2017 Feb 1.
Modeling flow in porous media is usually focused on the governing equations for mass and momentum transport, which yield the velocity and pressure at the pore or Darcy scales. However, in many applications, it is important to determine the work (or power) needed to induce flow in porous media, and this can be achieved when the mechanical energy equation is taken into account. At the macroscopic scale, this equation may be postulated to be the result of the inner product of Darcy's law and the seepage velocity. However, near the porous medium boundaries, this postulate seems questionable due to the spatial variations of the effective properties (velocity, permeability, porosity, etc.). In this work we derive the macroscopic mechanical energy equation using the method of volume averaging for the simple case of incompressible single-phase flow in porous media. Our analysis shows that the result of averaging the pore-scale version of the mechanical energy equation at the Darcy scale is not, in general, the expected product of Darcy's law and the seepage velocity. As a matter of fact, this result is only applicable in the bulk region of the porous medium and, in the derivation of this result, the properties of the permeability tensor are determinant. Furthermore, near the porous medium boundaries, a more novel version of the mechanical energy equation is obtained, which incorporates additional terms that take into account the rapid variations of structural properties taking place in this particular portion of the system. This analysis can be applied to multiphase and compressible flows in porous media and in many other multiscale systems.
多孔介质中的流动建模通常聚焦于质量和动量传输的控制方程,这些方程可得出孔隙尺度或达西尺度下的速度和压力。然而,在许多应用中,确定在多孔介质中引发流动所需的功(或功率)很重要,而考虑机械能方程时就能实现这一点。在宏观尺度上,该方程可假定为达西定律与渗流速度内积的结果。然而,在多孔介质边界附近,由于有效属性(速度、渗透率、孔隙率等)的空间变化,这一假定似乎存在问题。在这项工作中,我们针对多孔介质中不可压缩单相流的简单情况,采用体积平均法推导宏观机械能方程。我们的分析表明,在达西尺度下对机械能方程的孔隙尺度形式进行平均的结果,一般而言并非达西定律与渗流速度的预期乘积。事实上,该结果仅适用于多孔介质的主体区域,且在推导此结果时,渗透率张量的属性起决定性作用。此外,在多孔介质边界附近,会得到一个更新颖的机械能方程版本,其中包含考虑系统这一特定部分结构属性快速变化的附加项。该分析可应用于多孔介质中的多相流和可压缩流以及许多其他多尺度系统。