View source for Physics:Saffman–Taylor instability

[[File:Viscous fingers in a TiO2 sol-gel thin film formed from a Saffman-Taylor instability. 200x magnification reflected light.tif|thumb|An example of viscous fingering in a Hele-Shaw cell.]]
The '''Saffman–Taylor instability''', also known as '''viscous fingering''', is the formation of patterns in a morphologically unstable interface between two fluids in a porous medium, described mathematically by [[Biography:Philip Saffman|Philip Saffman]] and [[Biography:G. I. Taylor|G. I. Taylor]] in a paper of 1958.<ref name=":0">{{Cite journal|last1=Saffman|first1=Philip Geoffrey|last2=Taylor|first2=Geoffrey Ingram|date=1958-06-24|title=The penetration of a fluid into a porous medium or Hele-Shaw cell containing a more viscous liquid|url=https://royalsocietypublishing.org/doi/10.1098/rspa.1958.0085|journal=Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences|volume=245|issue=1242|pages=312–329|doi=10.1098/rspa.1958.0085|bibcode=1958RSPSA.245..312S|s2cid=95750900}}</ref><ref>{{Cite journal|last=Homsy|first=G M|date=1987-01-01|title=Viscous Fingering in Porous Media|url=https://www.annualreviews.org/doi/abs/10.1146/annurev.fl.19.010187.001415|journal=Annual Review of Fluid Mechanics|language=en|volume=19|issue=1|pages=271–311|doi=10.1146/annurev.fl.19.010187.001415|bibcode=1987AnRFM..19..271H|issn=0066-4189}}</ref> This situation is most often encountered during drainage processes through media such as soils.<ref>{{cite journal | last1=Li| first1=S|display-authors=0| title= Dynamics of Viscous Entrapped Saturated Zones in Partially Wetted Porous Media | journal= Transport in Porous Media | year=2018 | volume=125| issue= 2| pages= 193–210| doi= 10.1007/s11242-018-1113-3 | arxiv=1802.07387| s2cid=53323967}}</ref>  It occurs when a less viscous fluid is injected, displacing a more viscous fluid; in the inverse situation, with the more viscous displacing the other, the interface is stable and no instability is seen. Essentially the same effect occurs driven by gravity (without injection) if the interface is horizontal and separates two fluids of different densities, the heavier one being above the other: this is known as the Rayleigh-Taylor instability. In the rectangular configuration the system evolves until a single finger (the Saffman–Taylor finger) forms, whilst in the radial configuration the pattern grows forming fingers by successive tip-splitting.<ref>{{Cite journal|last1=Lajeunesse|first1=E.|last2=Couder|first2=Y.|date=2000-09-01|title=On the tip-splitting instability of viscous fingers|url=https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/abs/on-the-tipsplitting-instability-of-viscous-fingers/C174B7B4E8F0D20C6426A00C7ECEF846|journal=Journal of Fluid Mechanics|language=en|volume=419|issue=1|pages=125–149|doi=10.1017/S0022112000001324|bibcode=2000JFM...419..125L|s2cid=121812154|issn=1469-7645}}</ref>

Most experimental research on viscous fingering has been performed on Hele-Shaw cells, which consist of two closely spaced, parallel sheets of glass containing a viscous fluid. The two most common set-ups are the channel configuration, in which the less viscous fluid is injected at one end of the channel, and the radial configuration, in which the less viscous fluid is injected at the centre of the cell.  Instabilities analogous to viscous fingering can also be self-generated in biological systems.<ref name="Mather">{{Cite journal | doi = 10.1103/PhysRevLett.104.208101| title = Streaming Instability in Growing Cell Populations| journal = Physical Review Letters| volume = 104| issue = 20| pages = 208101| year = 2010| last1 = Mather | first1 = W. | last2 = Mondragón-Palomino | first2 = O. | last3 = Danino | first3 = T. | last4 = Hasty | first4 = J. | last5 = Tsimring | first5 = L. S. |bibcode = 2010PhRvL.104t8101M | pmid=20867071 | pmc=2947335}}</ref>

== Derivation for a planar interface ==
The simplest case of the instability arises at a planar interface within a porous medium or Hele-Shaw cell, and was treated by Saffman and Taylor<ref name=":0" /> but also earlier by other authors.<ref>{{Cite journal|last=Hill|first=S.|date=1952|title=Channeling in packed columns|url=https://dx.doi.org/10.1016/0009-2509%2852%2987017-4|journal=Chemical Engineering Science|language=en|volume=1|issue=6|pages=247–253|doi=10.1016/0009-2509(52)87017-4|issn=0009-2509}}</ref> A fluid of viscosity <math>\mu_1</math> is driven in the <math>x</math>-direction into another fluid of viscosity <math>\mu_2</math> at some velocity <math>V</math>. Denoting the permeability of the porous medium as a constant, isotropic, <math>\Pi</math>, [[Physics:Darcy's law|Darcy's law]] gives the unperturbed pressure fields in the two fluids <math>i=1,\,2</math> to be<math display="block">p^{(0)}_i = p_{\text{int}} - \frac{V \mu_i}{\Pi}x,</math>where <math>p_{\text{int}}</math> is the pressure at the planar interface, working in a frame where this interface is instantaneously given by <math>x=0</math>. Perturbing this interface to <math>x=\eta_0 \exp{\left(\mathrm{i}ky + \sigma t\right)}</math> (decomposing into [[Physics:Normal mode|normal mode]]s in the <math>x-y</math> plane, and taking <math>\left|\eta_0\right| \ll 1</math>), the pressure fields become<math display="block">p_i = p^{(0)}_i + \tilde{p}_i\left(x\right)\exp{\left(\mathrm{i}ky+\sigma t\right)}.</math>As a consequence of the incompressibility of the flow and Darcy's law, the pressure fields must be [[Harmonic function|harmonic]], which, coupled with the requirement that the perturbation decay as <math>x \to \pm \infty</math>, fixes <math>\tilde{p}_1 = \tilde{p}_1 e^{kx}</math> and <math>\tilde{p}_2 = \tilde{p}_2 e^{-kx}</math>, with the constants <math>\tilde{p}</math> to be determined by continuity of pressure. Upon [[Linearization|linearization]], the kinematic boundary condition at the interface (that fluid velocity in the <math>x</math> direction must match the velocity of the fluid interface), coupled with Darcy's law, gives<math display="block">-\frac{\Pi}{\mu_i}\left.\frac{\partial \tilde{p}_i}{\partial x}\right|_{x=0} = \sigma \eta_0,</math>and thus that <math>\tilde{p}_1 = -\frac{\sigma\eta_0\mu_1}{\Pi k}</math> and <math>\tilde{p}_2 = \frac{\sigma\eta_0\mu_2}{\Pi k}</math>. Matching the pressure fields at the interface gives<math display="block">-V \mu_1 - \frac{\sigma \mu_1}{k} = -V \mu_2 + \frac{\sigma \mu_2}{k},</math>and so <math>\sigma = kV\left(\mu_2 - \mu_1\right)/\left(\mu_1 + \mu_2\right)</math>, leading to growth of the perturbation when <math>\mu_2 > \mu_1</math> - i.e. when the injected fluid is less viscous than the ambient fluid. There are problems with this basic case: namely that the most unstable mode has infinite wavenumber <math>k</math> and grows at an infinitely fast rate, which can be rectified by the introduction of [[Physics:Surface tension|surface tension]]<ref>{{Cite journal|last1=Chuoke|first1=R. L.|last2=van Meurs|first2=P.|last3=van der Poel|first3=C.|date=1959-12-01|title=The Instability of Slow, Immiscible, Viscous Liquid-Liquid Displacements in Permeable Media|url=https://www.onepetro.org/journal-paper/SPE-1141-G|journal=Transactions of the AIME|language=en|volume=216|issue=1|pages=188–194|doi=10.2118/1141-G|issn=0081-1696}}</ref> (which provides a jump condition in pressures across the fluid interface through the [[Young–Laplace equation]]), which has the effect of modifying the growth rate to

<math display="block">\sigma = \frac{kV\left(\mu_2 - \mu_1\right) - \gamma H_f k^3}{\mu_1 + \mu_2},</math>

with surface tension <math>\gamma</math> and <math>H_f</math> the [[Mean curvature#Mean curvature in fluid mechanics|mean curvature]]. This suppresses small-wavelength (high-wavenumber) disturbances, and we would expect to see instabilities with wavenumber <math>k</math> close to the value of <math>k</math> which results in the maximal value of <math>\sigma</math>; in this case with surface tension, there is a unique maximal value.

== In radial geometry ==
The Saffman–Taylor instability is usually seen in an axisymmetric context as opposed to the simple planar case derived above.<ref>{{Cite journal|last=Wilson|first=S. D. R|date=1975-06-01|title=A note on the measurement of dynamic contact angles|url=https://dx.doi.org/10.1016/0021-9797%2875%2990151-4|journal=Journal of Colloid and Interface Science|language=en|volume=51|issue=3|pages=532–534|doi=10.1016/0021-9797(75)90151-4|bibcode=1975JCIS...51..532W|issn=0021-9797}}</ref><ref>{{Cite journal|last=Paterson|first=Lincoln|date=1981-12-01|title=Radial fingering in a Hele Shaw cell|url=https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/abs/radial-fingering-in-a-hele-shaw-cell/17685AC1EA92429CFC7B4E94CD824B04|journal=Journal of Fluid Mechanics|language=en|volume=113|pages=513–529|doi=10.1017/S0022112081003613|bibcode=1981JFM...113..513P|s2cid=122222282 |issn=1469-7645}}</ref> The mechanisms for the instability remain the same in this case, and the selection of the most unstable wavenumber in this case corresponds to a given number of fingers (an integer).

==See also==
*[[Physics:Kelvin–Helmholtz instability|Kelvin–Helmholtz instability]]
*[[Physics:Darcy's law|Darcy's law]]

==References==
{{Reflist}}<br />

{{DEFAULTSORT:Saffman-Taylor instability}}
[[Category:Fluid dynamics]]

{{Sourceattribution|Saffman–Taylor instability}}