Practical Considerations for Selection of Representative Elementary Volumes for Fluid Permeability in Fibrous Porous Media

We investigate the lower bound of the area of a square-shaped representative elementary volume (REV) for the permeability tensor for transverse Stokes flow through randomly packed, parallel, and monodisperse cylinders. The investigation is significant to flow models using small calculation regions for fibrous porous media, such as modeling defect formation during directional solidification in the mushy zone of dendritic alloys. Using 90 ensembles of 1,000 domains, where each ensemble comprises domains with the same number and size of cylinders, we develop correlations between the permeability tensor invariants and macroscopic features of the domain. We find that for ensembles of domains with fewer than 200 cylinders, the eigenvectors of the permeability tensors exhibit preferential alignment with the domain axes, demonstrating that the estimated permeability is significantly affected by the periodic boundary conditions for these cases. Our results also suggest that the anisotropy of the permeability tensor may not be insignificant even for large sampling volumes. These results provide a practical lower bound for the calculation volumes used in permeability simulations in fibrous porous media, and also suggest that modelers should consider using an anisotropic tensor for small calculation volumes if phenomena such as channeling are important.