Inertial Sensitivity of Porous Microstructures

• Published in: Transport in porous media 2018-11, Vol.125 (2), p.211-238
• Main Authors: Pauthenet, Martin, Davit, Yohan, Quintard, Michel, Bottaro, Alessandro
• Format: Article
• Language: English
• Subjects: Civil Engineering; Classical and Continuum Physics; Computational fluid dynamics; Computer simulation; Deviation; Drag; Earth and Environmental Science; Earth Sciences; Flow simulation; Fluid flow; Fluid mechanics; Geotechnical Engineering & Applied Earth Sciences; Hydrogeology; Hydrology/Water Resources; Industrial Chemistry/Chemical Engineering; Inertia; Inertial sensing devices; Mechanics; Microstructure; Parameter sensitivity; Physics; Porous media; Porous media flow; Reynolds number; Sensitivity; Simulation; Turbulent flow
• ISSN: 0169-3913; 1573-1634
• DOI: 10.1007/s11242-018-1115-1

Fluid flows through porous media are subject to different regimes, ranging from linear creeping flows to unsteady, chaotic turbulence. These different flow regimes at the pore scale have repercussions at larger scales, with the macroscale drag force experienced by a fluid moving through the medium becoming a nonlinear function of the average velocity beyond the creeping flow regime. Accurate prediction of the transition between different flow regimes is an important challenge with repercussions onto many engineering applications. Here, we are interested in the first deviation from Darcy’s law, when inertia effects become sizeable. Our goal is to define a Reynolds number, R e C , so that the inertial deviation occurs when R e C ∼ 1 for any microstructure. The difficulty in doing so is to reduce the multiple length scales characterizing the geometry of the porous structure to a single length scale, ℓ . We analyze the problem using the method of volume averaging and identify a length scale in the form ℓ = C λ K λ ϵ β , with C λ a parameter that indicates the sensitivity of the microstructure to inertia. The main advantage of this definition is that an explicit formula for C λ is given; C λ is computed from a creeping flow simulation in the porous medium; and R e C can be used to predict the transition to a non-Darcian regime more accurately than by using Reynolds numbers based on alternative length scales. The theory is validated numerically with data from flow simulations for a variety of microstructures.