On the stationary macroscopic inertial effects for one phase flow in ordered and disordered porous media

Abstract

We report on the controversial dependence of the inertial correction to Darcy’s law upon the filtration velocity (or Reynolds number) for one-phase Newtonian incompressible flow in model porous media. Our analysis is performed on the basis of an upscaled form of the Navier-Stokes equation requiring the solution of both the micro-scale flow and the associated closure problem. It is carried out with a special focus on the different regimes of inertia (weak and strong inertia) and the crossover between these regimes versus flow orientation and structural parameters, namely porosity and disorder. For ordered structures, it is shown that (i) the tensor involved in the expression of the correction is generally not symmetric, despite the isotropic feature of the permeability tensor. This is in accordance with the fact that the extra force due to inertia exerted on the structure is not pure drag in the general case; (ii) the Forchheimer type of correction (which strictly depends on the square of the filtration velocity) is an approximation that does not hold at all for particular orientations of the pressure gradient with respect to the axes of the structure; and (iii) the weak inertia regime always exists as predicted by theoretical developments. When structural disorder is introduced, this work shows that (i) the quadratic dependence of the correction upon the filtration velocity is very robust over a wide range of the Reynolds number in the strong inertia regime; (ii) the Reynolds number interval corresponding to weak inertia, that is always present, is strongly reduced in comparison to ordered structures. In conjunction with its relatively small magnitude, it explains why this weak inertia regime is most of the time overlooked during experiments on natural media. In all cases, the Forchheimer correction implies that the permeability is different from the intrinsic one.International audienceWe report on the controversial dependence of the inertial correction to Darcy’s law upon the filtration velocity (or Reynolds number) for one-phase Newtonian incompressible flow in model porous media. Our analysis is performed on the basis of an upscaled form of the Navier-Stokes equation requiring the solution of both the micro-scale flow and the associated closure problem. It is carried out with a special focus on the different regimes of inertia (weak and strong inertia) and the crossover between these regimes versus flow orientation and structural parameters, namely porosity and disorder. For ordered structures, it is shown that (i) the tensor involved in the expression of the correction is generally not symmetric, despite the isotropic feature of the permeability tensor. This is in accordance with the fact that the extra force due to inertia exerted on the structure is not pure drag in the general case; (ii) the Forchheimer type of correction (which strictly depends on the square of the filtration velocity) is an approximation that does not hold at all for particular orientations of the pressure gradient with respect to the ax...

Transitional and Turbulent Flow in a Bed of Spheres as Measured with Stereoscopic Particle Image Velocimetry

“…The general tactic is to derive and investigate additional terms added to Darcy law so that high Re laminar flow through porous media can be predicted on a global level. Also, turbulent flow through porous media has been treated theoretically and a quadratic dependence of the correction term for velocity has been derived in many studies (Skjetne and Auriault 1999;Lasseux et al 2011). In Jouybari et al (2016), the model by Nakayama and Kuwahara (1999) is extended to low Re turbulent flows by using a standard numerical model for a porous medium consisting of a staggered arrangement of square cylinders.…”

Effects of porosity and inertia on the apparent permeability tensor in fibrous media

“…The directional decoupling and the diagonal property of the apparent permeability tensor have also been computationally demonstrated on a completely different REV geometry by Soulaine and Quintard (2014) . Conversely, Lasseux et al (2011) have carried out a two-dimensional study with fibers of square cross-section, finding that the off-diagonal terms are non-negligible and only about one order of magnitude smaller than the diagonal components. This result is a consequence of the non-rotationally-invariant geometry considered.…”

“…For the purposes of the present paper, we must retain that Lasseux et al (2011) have parametrized the Forchheimer correction with the Reynolds number, and have found that the inertial correction is orders of magnitude smaller than the Darcy's term, at least before the turbulent regime sets in. Moreover, Lasseux and co-workers have studied how a Forchheimer tensor, F , depends upon the direction of the macroscopic forcing term with respect to the orientation of the square cross-section of the fibers, for Re d up to 30.…”

“…Another contribution which deserves mention is that by Lasseux et al (2011) ; they have computed the permeability tensor for various Reynolds numbers, in a two-dimensional geometry with cylinders of square cross-section. Forcing the flow along the main symmetric directions of the fiber, Lasseux et al (2011) have identified different regimes:…”

“…Forcing the flow along the main symmetric directions of the fiber, Lasseux et al (2011) have identified different regimes:…”