| Numerical models for the simulation of longitudinal and transverse Stokes flow in cylindrical periodic porous media are presented. The models, which are based on a finite-volume formulation in primitive variables, utilize digital image representations of the geometries to simulate, making them particularly well-suited for the rapid automated analysis of creeping flow in porous media with complex morphologies. Complete details of the model formulations are given, including extensive treatment of the pressure boundary conditions at the solid-liquid interface needed to guarantee convergence with all possible geometries. The convergence behavior of both models is tested, and the models are shown to be second-order accurate.; The models are used to simulate flow over the whole range of volume fractions of liquid in several regular geometries. The longitudinal model is used to simulate flow in square arrays of circular and square ducts, and both models are used to simulate flow in square and hexagonal arrays of circular cylinders and square arrays of square cylinders rotated by varying amounts. For each of the geometries, accurate empirical expressions for the Darcy permeability as a function of volume fraction solid are presented. Where applicable, model predictions of permeability are compared to existing analytical results.; Subsequently, the models are used to simulate Stokes flow in random domains over a wide range of fractions liquid. The sequential random packing algorithm is used to generate 1000 random packings of circular cylinders at each of 14 fractions of liquid, and longitudinal and transverse flow simulations are performed for each geometry. Histograms and summary statistics are computed for the permeability for each fraction liquid, and empirical expressions for mean permeability as a function of fraction liquid are given. The autocorrelation structure of the geometry and of the fluid velocity is analyzed, and an analysis of the scaling of longitudinal permeability variance is presented. In transverse flow at high packing densities, it is found that lightning-like patterns emerge in the fluid velocity. It is also found that the details of flows in such geometries are strongly sensitive to the placement of individual solid obstacles. |
