Saffman–Taylor instability

The Saffman–Taylor instability occurs when a more mobile fluid displaces a less mobile fluid in a porous medium or Hele–Shaw cell. Viscosity contrast amplifies interface corrugations, while surface tension suppresses sufficiently short wavelengths.

For two fluids with dynamic viscosities ${\mu }_1<{\mu }_2$, uniform permeability of a porous medium $k$, porosity $\varphi$, and imposed Darcy velocity $U$, the Saffman–Taylor instability of a planar interface has growth rate $\sigma =\alpha U({\mu }_2-{\mu }_1)/[\varphi ({\mu }_2+{\mu }_1)]$ without capillarity. With constant interfacial tension $\gamma$ and the stabilizing pressure-jump convention, a perturbation of transverse wavenumber $\alpha >0$ instead hasThe fastest-growing mode satisfies ${\alpha }_m^2=({\mu }_2-{\mu }_1)U/(3k\gamma )$.

For a sharp thermal front moving at fraction $\Gamma$ of the Darcy velocity $U$, a cold fluid of viscosity $\mu$ displacing its hot counterpart of viscosity $b\mu$ has the displayed growth rate in the local, unbounded planar approximation without thermal conduction. The harmonic pressure perturbations decay on either side. Pressure continuity and continuity of perturbed Darcy flux give the velocity response $\delta U=U|a|(b-1)\eta /(b+1)$; the thermal front's kinematic law is ${\eta }_t=\Gamma \delta U$. The prefactor is the thermal-front speed, not the material pore velocity. Increasing viscosity with heating therefore makes the trailing thermal interface unstable even when the leading displacement of the original fluid is stable.

Let fluid 1 lie above fluid 2, with mass densities ${\rho }_1<{\rho }_2$ and dynamic viscosities ${\mu }_1<{\mu }_2$. Define $\Delta \rho ={\rho }_2-{\rho }_1>0$, $\Delta \mu ={\mu }_2-{\mu }_1>0$, porosity $\varphi$, permeability of a porous medium $k$, and downward interface speed $W$. With Darcy flux $\varphi W$, no capillarity, and two unbounded layers, Darcy's law givesFor a Fourier mode of interface displacement, the two perturbation pressures decay exponentially into the respective layers. The kinematic boundary condition makes their amplitudes proportional to $\varphi {\mu }_1\sigma /(k|\alpha |)$ and $-\varphi {\mu }_2\sigma /(k|\alpha |)$. Linearized pressure continuity adds the base-gradient difference $\varphi W\Delta \mu /k-g\Delta \rho$. Eliminating the amplitudes proves the formula. The critical interface speed is $kg\Delta \rho /(\varphi \Delta \mu )$; the critical Darcy velocity is $kg\Delta \rho /\Delta \mu$. Without surface tension this model has unbounded short-wave growth rate on its unstable side.

Suppose the apparent capillary pressure is $p_1-p_2=\gamma (1+\beta x)({\kappa }_0+\nabla \cdot \mathbf{n})$, with normal towards the displaced fluid and positive apparent tension at the flat front $x=X$. Linearizing the Darcy law interface problem gives the local growth rateThe pressure increase experienced by a forward protrusion is stabilizing when ${\kappa }_0\beta >0$. It suppresses every positive wavenumber when $U\le k\gamma {\kappa }_0\beta /({\mu }_2-{\mu }_1)$. For a fixed spatial gradient, $X=Ut/\varphi$ makes the rate time dependent; constant exponential growth is a frozen-position approximation.