Hyperbolic regions in flows through three-dimensional pore structures.

Finite time Lyapunov exponents are used to determine expanding, contracting, and hyperbolic regions in computational simulations of laminar steady-state fluid flows within realistic three dimensional pore structures embedded within an impermeable matrix. These regions correspond approximately to pores where flow converges (contraction) or diverges (expansion), and to throats between pores where the flow mixes (hyperbolic). The regions are sparse and disjoint from one another, occupying only a small percentage of the pore space. Nonetheless, nearly every percolating fluid particle trajectory passes through several hyperbolic regions indicating that the effects of in-pore mixing are distributed throughout an entire pore structure. Furthermore, the observed range of fluid dynamics evidences two scales of heterogeneity within each of these flow fields. There is a larger scale that affects dispersion of fluid particle trajectories across the connected network of pores and a relatively small scale of nonuniform distributions of velocities within an individual pore.

Heterogeneities of flow in stochastically generated porous media.

Heterogeneous flows are observed to result from variations in the geometry and topology of pore structures within stochastically generated three dimensional porous media. A stochastic procedure generates media comprising complex networks of connected pores. Inside each pore space, the Navier-Stokes equations are numerically integrated until steady state velocity and pressure fields are attained. The intricate pore structures exert spatially variable resistance on the fluid, and resulting velocity fields have a wide range of magnitudes and directions. Spatially nonuniform fluid fluxes are observed, resulting in principal pathways of flow through the media. In some realizations, up to 25% of the flux occurs in 5% of the pore space depending on porosity. The degree of heterogeneity in the flow is quantified over a range of porosities by tracking particle trajectories and calculating their attributes including tortuosity, length, and first passage time. A representative elementary volume is first computed so the dependence of particle based attributes on the size of the domain through which they are followed is minimal. High correlations between the dimensionless quantities of porosity and tortuosity are calculated and a logarithmic relationship is proposed. As the porosity of a medium increases the flow field becomes more uniform.