View source for Viscous fingering

An interface between two fluids presents a resistance towards an increase of its surface area. This interface has a surface tension which opposes such an increase. This is the property giving a liquid drop its spherical shape: A sphere has a minimal surface for a given volume of fluid. T. Young and P.S. Laplace rationalized this scenario by postulating the Young–Laplace law.

In viscous fingering, where a less viscous fluid, such as air, pushes a more viscous fluid, such as an oil, in thin rectangular cells (known as Hele-Shaw cells), surface tension plays an important role in determining the shape and properties of the fingers. The basic experimental situation to obtain viscous fingers is as follows (here, the focus is on the linear channel geometry, but such experiments can also be obtained in circular cells). One usually uses two thick glass plates and a thin plastic sheet to obtain a linear Hele-Shaw cell [[#References|[a8]]]. The plastic sheet serves as a spacer between the glass plates once they are brought into contact with each other. A long rectangular canal can be cut in the middle of the plastic sheet, which is then sandwiched between the glass plates to obtain a thin linear rectangular channel. The thickness of this canal is fixed by the thickness of the plastic sheet. Through holes drilled in the upper glass plates one then injects the viscous fluid to fill up the cell. Subsequently, air is injected at constant pressure or constant flux to push the oil out of this canal (see Fig.a2 for a scheme of the cell). The interface between the air and the oil is then unstable: the flat interface takes a curved shape. This instability is well known and from the pioneering work of P.G. Saffman and G.I. Taylor in the 1950s [[#References|[a7]]], it is known that the finger of air has a width which decreases with the velocity of the finger. The most mysterious aspect observed by Saffman and Taylor is that the limiting width of the fingers obtained at high velocities is very close to half the channel width. This property was not explained until the 1980s, although Saffman and Taylor did calculate the shape of the finger they obtained in their experiments.

<img style="border:1px solid;" src="https://www.encyclopediaofmath.org/legacyimages/common_img/v130070a.gif" />

Figure: v130070a

A Saffman–Taylor finger obtained with air injected in a viscous fluid. This finger has a width very close to half the channel width. The gray region delimiting the finger contains the fluid, while the white region contains air. The black dots represent the shape obtained from the Saffman–Taylor solution

<img style="border:1px solid;" src="https://www.encyclopediaofmath.org/legacyimages/common_img/v130070b.gif" />

Figure: v130070b

Scheme of the linear Hele-Shaw cell: two thick glass plates with two tubes attached to inject the less viscous fluid and to recover the fluid being pushed. The gray part is the plastic sheet of thickness <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300701.png" /> into which a linear channel of width <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300702.png" /> is cut

See below for the basic equations governing this problem and some hints as to how it was resolved using techniques from complex analysis and perturbation theory in the 1980s. See [[#References|[a11]]], [[#References|[a6]]], [[#References|[a3]]] for accessible reviews of the experiment and phenomena.

The use of a Hele-Shaw cell permits one to use a simplified version of the equations governing the flow of the fluids. Basically, the [[Navier–Stokes equations|Navier–Stokes equations]] are simplified to keep just the pressure gradient term and the dissipative term. This approximation is known as Darcy's law and is obeyed by flows in porous media, for example. In a sense, the linear Hele-Shaw cell is a simple porous medium, as was shown by H.J.S. Hele-Shaw [[#References|[a8]]]. This equation is written as follows:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300703.png" /></td> </tr></table>

Here, the gradient is in the plane of the cell, denoted by <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300704.png" /> (the finger propagates along the <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300705.png" /> direction with <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300706.png" /> being the transverse direction), <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300707.png" /> is the pressure, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300708.png" /> is the fluid shear viscosity, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v1300709.png" /> is the velocity of the fluid in the <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007010.png" />-plane far from the fingertip and averaged over the vertical gap of thickness <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007011.png" /> (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007012.png" /> being the spacing between the plates fixed by the thickness of the plastic sheet). One can rewrite this equation in terms of a velocity potential <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007013.png" /> such that <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007014.png" />; due to the incompressibility of the fluid this potential obeys the [[Laplace equation|Laplace equation]] <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007015.png" />. Suppose the following boundary conditions on the interface apply:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007016.png" /></td> </tr></table>

Here, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007017.png" /> is the normal to the interface delimiting the finger, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007018.png" /> is the unit vector in the direction of propagation of the finger and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007019.png" /> is the velocity of the tip. The finger is supposed to advance at constant velocity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007020.png" /> and keeps a constant width <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007021.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007022.png" /> is the width of the channel and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007023.png" /> is the relative width of the finger [[#References|[a7]]]. The potential at the interface is

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007024.png" /></td> </tr></table>

Here, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007025.png" /> is the surface tension of the interface and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007026.png" /> is the curvature of the interface. The additional term in the potential expresses the effect of surface tension due to the Young–Laplace law. <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007027.png" /> is the potential inside the less viscous fluid (which can be taken as inviscid without loss of generality). It is clear that far behind the fingertip the potential is zero while far ahead of the fingertip the viscous fluid is moving with velocity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007028.png" />. The problem is now to determine the relation between the coordinates of the interface <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007029.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007030.png" />. In 1980, J.W. McLean and P.G. Saffman [[#References|[a5]]] suggested a way of rewriting the equations in complex representation which came in handy for later work. One writes the coordinates in the complex plane as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007031.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007032.png" /> (where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007033.png" /> is the [[Stream function|stream function]]). The transformation <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007034.png" /> is a [[Conformal mapping|conformal mapping]] from the flow region delimited by the finger and the walls of the channel into an infinite strip of unit width in the potential plane. The conformal mapping <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007035.png" /> maps the flow region into the upper <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007036.png" />-plane. This transforms the interface between the tails and the tip of the finger into the real segment <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007037.png" />. The boundary <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007038.png" /> corresponds to the tail of the finger, for which <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007039.png" />, and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007040.png" /> corresponds to the tip, for which <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007041.png" />. The problem now reduces to finding <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007042.png" /> in terms of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007043.png" />. However, it turns out to be simpler to study the quantity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007044.png" />; it is the complex velocity <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007045.png" /> and may be written as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007046.png" />, since the velocity is tangent to the interface (with <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007047.png" /> the angle between <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007048.png" /> and the tangent to the interface). The modulus <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007049.png" /> at the interface can be written as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007050.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007051.png" /> measures arclength along the interface starting from the tip. By invoking the analytic properties of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007052.png" />, the following relation can be established between <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007053.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007054.png" />:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007055.png" /></td> </tr></table>

Besides this relation there is a second one, dictated by the Young–Laplace law:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007056.png" /></td> </tr></table>

where

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007057.png" /></td> </tr></table>

These two equations are to be solved with the following boundary conditions:

at the tail of the finger: <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007058.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007059.png" />;

at the tip of the finger: <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007060.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007061.png" />. If a solution of these equations is found, the value of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007062.png" /> (which is the quantity that needs to be determined) as well as its dependence on physical parameters, such as the velocity of the finger, the viscosity of the fluid and especially the surface tension, can be found from the relation:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007063.png" /></td> </tr></table>

At this point the equations are complete. Of course, the difficulty is in solving these equations. For the case of zero surface tension (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007064.png" />), one finds a continuum of solutions for the finger shape, parametrized by the value of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007065.png" /> which remains unknown. The family of solutions first found by Saffman and Taylor is written as:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007066.png" /></td> </tr></table>

with <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007067.png" />. In the coordinates <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007068.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007069.png" />, this solution for the interface shape reads:

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007070.png" /></td> </tr></table>

The problem now is to find a selection mechanism for the allowed values of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007071.png" />. This problem was not resolved until 1986, when three different groups [[#References|[a9]]], [[#References|[a10]]], [[#References|[a1]]], [[#References|[a2]]] studied the small-<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007072.png" /> limit of these equations using techniques from perturbation theory and showed that the continuum of solutions breaks into a discrete set of solutions with the value of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007073.png" /> decreasing to its well-known value of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007074.png" /> as <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007075.png" /> approaches <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/v/v130/v130070/v13007076.png" />. This calculation is rather involved; see [[#References|[a9]]], [[#References|[a10]]], [[#References|[a1]]], [[#References|[a2]]].

The use of complex fluids such as, e.g., gels, clays, polymeric solutions, mixtures of water and surfactants can give rise to new morphologies of the interfaces in Hele-Shaw cells [[#References|[a6]]].

In the above, the viscosity of the fluid as well as the surface tension of the interface were considered to be constant. For complicated fluids these physical quantities may depend on the velocity of the fluid and therefore become variable. This complicates the mathematical analysis of the instability; however, some progress has been made [[#References|[a4]]].

====References====
<table><TR><TD valign="top">[a1]</TD> <TD valign="top">  B.I. Shraiman,   ''Phys. Rev. Lett.'' , '''56'''  (1986)  pp. 2028</TD></TR><TR><TD valign="top">[a2]</TD> <TD valign="top">  D.C. Hong,   J. Langer,   ''Phys. Rev. Lett.'' , '''56'''  (1986)  pp. 2032</TD></TR><TR><TD valign="top">[a3]</TD> <TD valign="top">  D. Bensimon,   L.P. Kadanoff,   S. Liang,   B.I. Shraiman,   C. Tang,   ''Rev. Mod. Phys.'' , '''58'''  (1986)  pp. 977</TD></TR><TR><TD valign="top">[a4]</TD> <TD valign="top">  D. Bonn,   H. Kellay,   J. Meunier,   ''Philos. Mag. B'' , '''78'''  (1998)  pp. 131–142</TD></TR><TR><TD valign="top">[a5]</TD> <TD valign="top">  J.W. McLean,   P.G. Saffman,   ''J. Fluid Mech.'' , '''102'''  (1981)  pp. 455</TD></TR><TR><TD valign="top">[a6]</TD> <TD valign="top">  K.V. McCloud,   J.V. Maher,   ''Phys. Rept.'' , '''260'''  (1995)  pp. 139</TD></TR><TR><TD valign="top">[a7]</TD> <TD valign="top">  P.G. Saffman,   G.I. Taylor,   ''Proc. Royal Soc. London A'' , '''245'''  (1958)  pp. 312</TD></TR><TR><TD valign="top">[a8]</TD> <TD valign="top">  H.J.S. Hele-Shaw,   ''Nature'' , '''58'''  (1898)  pp. 34</TD></TR><TR><TD valign="top">[a9]</TD> <TD valign="top">  R. Combescot,   T. Dombre,   V. Hakim,   Y. Pomeau,   A. Pumir,   ''Phys. Rev. Lett.'' , '''56'''  (1986)  pp. 2036</TD></TR><TR><TD valign="top">[a10]</TD> <TD valign="top">  R. Combescot,   T. Dombre,   V. Hakim,   Y. Pomeau,   A. Pumir,   ''Phys. Rev. A'' , '''37'''  (1988)  pp. 1270</TD></TR><TR><TD valign="top">[a11]</TD> <TD valign="top">  P. Pelcé,   A. Libchaber,   "Dynamics of curved fronts" , ''Perspectives in Physics'' , Acad. Press  (1997)</TD></TR></table>