TY - JOUR

T1 - Multiscale formulation of pore-scale compressible Darcy-Stokes flow

AU - Guo, Bo

AU - Mehmani, Yashar

AU - Tchelepi, Hamdi A.

N1 - Funding Information:
This work is supported in part by TOTAL through the Stanford TOTAL enhanced modeling of source rock (STEMS) project. We thank the Center for Computational Earth & Environmental Science (CEES) at Stanford University for access to computational resources.
Publisher Copyright:
© 2019 Elsevier Inc.

PY - 2019/11/15

Y1 - 2019/11/15

N2 - Direct numerical simulation (DNS) of fluid dynamics in digital images of porous materials is challenging due to the cut-off length issue where interstitial voids below the resolution of the imaging instrument cannot be resolved. Such subresolution microporosity can be critical for flow and transport because they could provide important flow pathways. A micro-continuum framework can be used to address this problem, which applies to the entire domain a single momentum equation, i.e., Darcy-Brinkman-Stokes (DBS) equation, that recovers Stokes equation in the resolved void space (i.e., macropores) and Darcy equation in the microporous regions. However, the DBS-based micro-continuum framework is computationally demanding. Here, we develop an efficient multiscale method for the compressible Darcy-Stokes flow arising from the micro-continuum approach. The method decomposes the domain into subdomains that either belong to the macropores or the microporous regions, on which Stokes or Darcy problems are solved locally, only once, to build basis functions. The nonlinearity from compressible flow is accounted for in a local correction problem on each subdomain. A global interface problem is solved to couple the local bases and correction functions to obtain an approximate global multiscale solution, which is in excellent agreement with the reference single-scale solution. The multiscale solution can be improved through an iterative strategy that guarantees convergence to the single-scale solution. The method is computationally efficient and well-suited for parallelization to simulate fluid dynamics in large high-resolution digital images of porous materials.

AB - Direct numerical simulation (DNS) of fluid dynamics in digital images of porous materials is challenging due to the cut-off length issue where interstitial voids below the resolution of the imaging instrument cannot be resolved. Such subresolution microporosity can be critical for flow and transport because they could provide important flow pathways. A micro-continuum framework can be used to address this problem, which applies to the entire domain a single momentum equation, i.e., Darcy-Brinkman-Stokes (DBS) equation, that recovers Stokes equation in the resolved void space (i.e., macropores) and Darcy equation in the microporous regions. However, the DBS-based micro-continuum framework is computationally demanding. Here, we develop an efficient multiscale method for the compressible Darcy-Stokes flow arising from the micro-continuum approach. The method decomposes the domain into subdomains that either belong to the macropores or the microporous regions, on which Stokes or Darcy problems are solved locally, only once, to build basis functions. The nonlinearity from compressible flow is accounted for in a local correction problem on each subdomain. A global interface problem is solved to couple the local bases and correction functions to obtain an approximate global multiscale solution, which is in excellent agreement with the reference single-scale solution. The multiscale solution can be improved through an iterative strategy that guarantees convergence to the single-scale solution. The method is computationally efficient and well-suited for parallelization to simulate fluid dynamics in large high-resolution digital images of porous materials.

UR - http://www.scopus.com/inward/record.url?scp=85070213174&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85070213174&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2019.07.047

DO - 10.1016/j.jcp.2019.07.047

M3 - Article

AN - SCOPUS:85070213174

VL - 397

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

M1 - 108849

ER -