The macroscopic properties of fluid flow and transport through porous media are a direct consequence of the underlying pore structure. However, precise relations that characterize flow and transport from the statistics of pore-scale disorder have remained elusive. Here we investigate the relationship between pore structure and the resulting fluid flow and asymptotic transport behavior in two-dimensional geometries of nonoverlapping circular posts. We derive an analytical relationship between the pore throat size distribution $f_{\lambda }\sim {\lambda }^{-\beta }$ and the distribution of the low fluid velocities $f_u\sim u^{-\beta /2}$, based on a conceptual model of porelets (the flow established within each pore throat, here a Hagen-Poiseuille flow). Our model allows us to make predictions, within a continuous-time random-walk framework, for the asymptotic statistics of the spreading of fluid particles along their own trajectories. These predictions are confirmed by high-fidelity simulations of Stokes flow and advective transport. The proposed framework can be extended to other configurations which can be represented as a collection of known flow distributions.

Transport of viscous fluid through porous media is a direct consequence of the pore structure. Here we investigate transport through a specific class of two-dimensional porous geometries, namely those formed by fluid-mechanical erosion. We investigate the tortuosity and dispersion by analyzing the first two statistical moments of tracer trajectories. For most initial configurations, tortuosity decreases in time as a result of erosion increasing the porosity. However, we find that tortuosity can also increase transiently in certain cases. The porosity-tortuosity relationships that result from our simulations are compared with models available in the literature. Asymptotic dispersion rates are also strongly affected by the erosion process, as well as by the number and distribution of the eroding bodies. Finally, we analyze the pore size distribution of an eroding geometry. The simulations are performed by combining a boundary integral equation solver for the fluid equations, a second-order stable time-stepping method to simulate erosion, and high-order numerical methods to stably and accurately resolve nearly touching eroded bodies and particle trajectories near the eroding bodies.

We study mixing in Stokesian vesicle suspensions in two dimensions on a cylindrical Couette apparatus using numerical simulations. The vesicle flow simulation is done using a boundary integral method, and the advection-diffusion equation for the mixing of the solute is solved using a pseudo-spectral scheme. We study the effect of the area fraction, the viscosity contrast between the inside (the vesicles) and the outside (the bulk) fluid, the initial condition of the solute, and the mixing metric. We compare mixing in the suspension with mixing in the Couette apparatus without vesicles. On the one hand, the presence of vesicles in most cases slightly suppresses mixing. This is because the solute can be only diffused across the vesicle interface and not advected. On the other hand, there exist spatial distributions of the solute for which the unperturbed Couette flow completely fails to mix whereas the presence of vesicles enables mixing. We derive a simple condition that relates the velocity and solute and can be used to characterize the cases in which the presence of vesicles promotes mixing.

Lipid bilayer membranes have a native (albeit small) permeability for water molecules. Under an external load, provided that the bilayer structure stays intact and does not suffer from poration or rupture, a lipid membrane deforms and its water influx/efflux is often assumed negligible in the absence of osmolarity. In this work we use boundary integral simulations to investigate the effects of water permeability on the hydrodynamics of an inextensible membrane under a mechanical load, such as the viscous stress from an external flow deforming an inextensible membrane in free space or pushing it through a confinement. Incorporating the membrane permeability into the framework of Helfrich free energy for an inextensible, elastic membrane (a vesicle), we illustrate that, in the absence of an osmotic stress gradient, the semipermeable vesicle is affected by water influx/efflux over a sufficiently long time or under a strong confinement. Our simulations quantify the conditions for water permeation to be negligible in terms of the timescales, flow strength, and confinement. These results shed light on how microfluidic confinement can be utilized to estimate membrane permeability.