TY - JOUR
T1 - Characterization of transition from Darcy to non-Darcy flow with 3D pore-level simulations
AU - Tachibana, Ikkoh
AU - Moriguchi, Shuji
AU - Takase, Shinsuke
AU - Terada, Kenjiro
AU - Aoki, Takayuki
AU - Kamiya, Kohji
AU - Kodaka, Takeshi
N1 - Funding Information:
This research was mainly supported by JSPS KAKENHI Grant No. JP16J01584 , and partly supported by KAKENHI Grant-in-Aid for Scientific Research (S) 26220002 from Japan Society for the Promotion of Science (JSPS).
Publisher Copyright:
© 2017
PY - 2017/10
Y1 - 2017/10
N2 - We characterize seepage flow in porous media via a series of three-dimensional (3D) direct numerical simulations (DNSs) of Navier-Stokes flows in representative volume elements (RVEs) at the pore level. The immersed boundary method with a fixed staggered grid is employed for calculations based on the finite difference method in the pore domain formed with spatially arranged rigid particles in RVEs. The numerical results of the DNSs with different particle sizes and different seepage flow velocities are volume-averaged over the RVEs to evaluate the permeability coefficients of the seepage flows and are ordered according to the Reynolds numbers for porous media. After the numerical scheme is validated with a simple RVE model, the transition from Darcy to non-Darcy flow, for which experimental results have reported in the literature, is demonstrated in terms of the calculated permeability coefficients and examined based on the pore-level flow characteristics to clarify the mechanism of the permeability change, i.e. the macroscopic behaviour of RVEs related to the Reynolds number.
AB - We characterize seepage flow in porous media via a series of three-dimensional (3D) direct numerical simulations (DNSs) of Navier-Stokes flows in representative volume elements (RVEs) at the pore level. The immersed boundary method with a fixed staggered grid is employed for calculations based on the finite difference method in the pore domain formed with spatially arranged rigid particles in RVEs. The numerical results of the DNSs with different particle sizes and different seepage flow velocities are volume-averaged over the RVEs to evaluate the permeability coefficients of the seepage flows and are ordered according to the Reynolds numbers for porous media. After the numerical scheme is validated with a simple RVE model, the transition from Darcy to non-Darcy flow, for which experimental results have reported in the literature, is demonstrated in terms of the calculated permeability coefficients and examined based on the pore-level flow characteristics to clarify the mechanism of the permeability change, i.e. the macroscopic behaviour of RVEs related to the Reynolds number.
KW - Immersed boundary method
KW - Micro-scale
KW - Non-Darcy flow IGC: E07
KW - Permeability
UR - http://www.scopus.com/inward/record.url?scp=85030641756&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85030641756&partnerID=8YFLogxK
U2 - 10.1016/j.sandf.2017.08.003
DO - 10.1016/j.sandf.2017.08.003
M3 - Article
AN - SCOPUS:85030641756
VL - 57
SP - 707
EP - 719
JO - Soils and Foundations
JF - Soils and Foundations
SN - 0038-0806
IS - 5
ER -