Dynamic multiphase flow triggers chaotic mixing in porous media

Solute mixing plays a pivotal role in a broad spectrum of chemical and biological processes across natural and engineered porous media. However, current understanding of mixing dynamics remains largely constrained to steady flows in fully or partially water-saturated environments. Multiphase flow systems are generally unsteady, with moving fluid interfaces and flow paths that change in time. Despite the widespread occurrence of dynamic multiphase flows, their impacts on solute mixing are largely unknown. Here, we use experiments and numerical simulations to investigate the effect of dynamic two-phase flow on the stretching and folding of fluid elements, a fundamental mechanism driving solute mixing and reactions in porous media. We find that dynamic two-phase flows induce chaotic mixing, characterized by exponential stretching of fluid elements, leading to strongly enhanced mixing compared to steady single-phase flows. By extensive numerical multiphase flow simulations, we establish dynamic steady states where we reliably measure the mean fluid stretching rate as a function of flow rate. We show that stretching is maximized at an optimum flow rate which balances fluid shear deformation against the frequency of flow reorientation by the intermittent motion of the fluid interface. The findings are rationalized by a mechanistic model linking basic multiphase flow characteristics to the stretching rate, opening perspectives to understand and control mixing and reactions in a wide range of multiphase flow systems.