Variable density and viscosity, miscible displacements in horizontal Hele-Shaw cells. Part 1. Linear stability analysis

A computational investigation of variable density and viscosity, miscible displacements in horizontal Hele-Shaw cells is presented. As a first step, two-dimensional base states are obtained by means of simulations of the Stokes equations, which are nonlinear due to the dependence of the viscosity on the local concentration. Here, the vertical position of the displacement front is seen to reach a quasisteady equilibrium value, reflecting a balance between viscous and gravitational forces. These base states allow for two instability modes: first, there is the familiar tip instability driven by the unfavourable viscosity contrast of the displacement, which is modulated by the presence of density variations in the gravitational field; second, a gravitational instability occurs at the unstably stratified horizontal interface along the side of the finger. Both of these instability modes are investigated by means of a linear stability analysis. The gravitational mode along the side of the finger is characterized by a wavelength of about one half to one full gap width. It becomes more unstable as the gravity parameter increases, even though the interface is shifted closer to the wall. The growth rate is largest far behind the finger tip, where the interface is both thicker, and located closer to the wall, than near the finger tip. The competing influences of interface thickness and wall proximity are clarified by means of a parametric stability analysis. The tip instability mode represents a gravity-modulated version of the neutrally buoyant mode. The analysis shows that in the presence of density stratification its growth rate increases, while the dominant wavelength decreases. This overall destabilizing effect of gravity is due to the additional terms appearing in the stability equations, which outweigh the stabilizing effects of gravity onto the base state.