Impact of structured heterogeneities on reactive two-phase porous flow

Two-phase flow through heterogeneous media leads to scale-free distributions of irregularly shaped pockets of one fluid trapped within the other. Although reactions within these fluids are often modeled at the homogeneous continuum scale, there exists no current framework for upscaling from the pore scale that accounts for the complex and scale-free geometry of the bubbles. In this paper, we apply a linear-kinetics reaction-diffusion model to characterize the steady-state chemical environment inside the irregular pockets. Using a combination of theory and invasion-percolation simulations, we derive scaling laws describing the distribution of diffusion times within bubbles. We show that chemical concentrations within the bubbles are determined by the Laplace transform of the entire distribution of diffusion times from each location. This serves as a means to compute average concentrations of reactant within a bubble of unique geometry and size. Furthermore, the overall system size imposes upper bounds on the distribution of bubble sizes, thereby imposing a system-size dependence on the statistics and average concentrations. These conclusions have profound implications for continuum models of porous reactive flow, where kinetic and equilibrium parameters are often chosen from laboratory measurements made at centimeter scales.