Numerical simulation and analysis of mixing enhancement due to chaotic advection using an adaptive approach for approximating the dilution index

A velocity field characterized by chaotic advection induces stretching and folding processes that increase the solute-solvent interface available for diffusion. Identifying chaotic flow fields with optimized mixing enhancement is relevant for applications like groundwater remediation, microfluidics, and many more. This work uses the dilution index to quantify the temporal increase in mixing of a solute within its solvent. We introduce a new approach to select a suitable grid size for each time step in the dilution index approximation, motivated by the theory of representative elementary volumes. It preserves the central feature of the dilution index, which is monotonically increasing in time and hence leads to reliable results. Our analysis highlights the importance of a suitable choice for the grid size in the dilution index approximation. We use this approach to demonstrate the mixing enhancement for two chaotic injection-extraction systems that exhibit chaotic structures: a source-sink dipole and a rotated potential mixing. By analyzing the chaotic flow fields, we identify Kolmogorov--Arnold--Moser (KAM) islands, non-mixing regions that limit the chaotic area in the domain and, thereby, the mixing enhancement. Using our new approach, we assess the choice of design parameters of the injection-extraction systems to effectively engineer chaotic mixing. We demonstrate the important role of diffusion in filling the KAM islands and reaching complete mixing in the systems.