A simulation approach including under-resolved scales for multi-component fluid flows in multi-scale porous structures

TL;DR

This paper introduces a simulation method for multi-component fluid flows in multi-scale porous structures that accounts for under-resolved regions using physics-based modeling, enabling accurate flow predictions without high computational costs.

Contribution

The study presents a novel simulation approach that incorporates physics-based models for under-resolved scales in multi-scale porous media, improving accuracy and efficiency.

Findings

Accurately captures imbibition patterns and entry pressure.

Produces quantitatively consistent permeability results.

Effective in both single- and multi-component flow simulations.

Abstract

In this study, we develop computational models and methodology for accurate multi-component-flow simulation in under-resolved multi-scale porous structures. It is generally impractical to fully resolve the flow in porous structures with large length-scale difference due to tremendously high computational expense. The flow contributions from under-resolved scales need to be accounted for with proper physics modeling as well as simulation processes. Using pre-computed physical properties such as the absolute permeability, K0, the capillary-pressure-saturation curve, and the relative permeability, Kr, in typically resolved porous structures, local fluid force is conjectured and applied to simulation in the under-resolved regions that are represented by porous media. By doing so, accurate simulation of flow in multi-scale porous structures becomes feasible. In order to check the accuracy and robustness of this method, a set of benchmark test cases are performed for both single-component and multi-component flows in artificially constructed multi-scale porous structures, and simulation results are compared with analytic solutions and/or results with much finer resolution resolving the porous structures. Quantitatively consistent results are obtained with proper input of K0, capillary pressure, and Kr in all tested cases. Specifically, imbibition patterns, entry pressure, residual component patterns, and the absolute and relative permeability are accurately captured with this approach.