Fluid flow through 3D-printed particle beds: a new technique for understanding, validating, and improving predictability of permeability from empirical equations

Abstract

The application of 3D technology for fabrication of artificial porous media samples improves porous media flow studies. The geometrical characteristics of a porous media pore channel: the channel shape, size, porosity, specific surface, expansion ratio, contraction ratio, and the tortuous pathway of the channel can be controlled through advanced additive manufacturing techniques (3D printing), computed tomography imagery (CT imaging) and image analysis methods. These 3D technologies have here been applied to construct and analyze four homogeneous porous media samples with predefined geometrical properties that are otherwise impossible to construct with conventional methods. Uncertainties regarding the geometrical properties are minimized because the 3D-printed porous media samples can be evaluated with CT imaging after fabrication. It is this combination of 3D technology that improves the data acquisition and data interpretation and contributes to new insight into the phenomenon of fluid flow through porous media. The effects of the individual geometrical properties on the fluid flow are then accounted for in permeability experiments in a Hassler flow cell. The results of the experimental work are used to test the theoretical foundation of the Kozeny–Carman equation and the extended version known as the Ergun equation. These equations are developed from analogies to the Hagen–Poiseuille flow equation. Based on the results from the laboratory experiments in this study, an analytical equation based on the analytical Navier–Stokes equations is presented as an alternative to the Hagen–Poiseuille analogy for porous media channels with non-uniform channel geometries. The agreement between experiment and the new equation reveals that the dissipating losses of mechanical energy in porous media flows are not a result of frictional shear alone. The mechanical losses are also a result of pressure dissipation that arise due to the non-uniformity of the channel geometry, which induced spatial variations to the strain rate field and induce acceleration of the velocity field in the flow through the porous medium. It is this acceleration that causes a divergence from linear flow conditions as the Stokes flow criterion (Re ≪ 1) is breached and causes the convective acceleration term to affect the flow behavior. The suggested modifications of theory and the presented experiments prove that the effects of surface roughness (1) do not alter the flow behavior in the Darcy flow regime or (2) in the Forchheimer flow regime. This implies that the flow is still laminar for the Forchheimer flow velocities tested.