The stability of immiscible viscous fingering in Hele-Shaw cells with spatially varying permeability

In this paper, we investigate the stability of immiscible viscous fingering in Hele-Shaw cells with spatially varying permeability, across a range of capillary numbers. We utilise a coupled boundary element - radial basis function (BE –RBF) numerical method that adapts and moves with the growing interface, providing an efficient, high accuracy scheme to track the interfacial displacement of immiscible fluids. By comparing the interfacial evolution and growth rate in varying permeability cells to that in uniform cells, we can assess the relative stability of the perturbations as a consequence of the variable permeability.Numerical experiments in Hele-Shaw cells with gradually varying permeability highlight 3 aperture effects that control the interfacial stability: (1) Gradients in the capillary pressure (2) Local changes in fluid mobility (3) Variation in the viscous pressure gradient. In low capillary number regimes, we find that aperture effect 1 and 2 dominate, which (relatively) stabilise interfacial perturbations in converging geometries and destabilise perturbations in diverging geometries. In high capillary number regimes, aperture effect 3 dominates meaning the relative stability transitions; the interface is destabilised in converging cells and stabilised in diverging cells. We find an upper bound critical capillary number Cagt at which the relative stability transitions in our gradually varying cell as 100