View source for Buckley-Leverett equation

An equation for fluid flow in a porous medium, proposed in [[#References|[a1]]], that models the displacement of oil by water in a one-dimensional porous medium (e.g., sand). If <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b1109901.png" /> denotes the saturation of water, i.e., the fraction of pore volume occupied by water, and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b1109902.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b1109903.png" /> are the mobilities of water and oil, respectively, i.e., the ratio of relative permeability and viscosity, then the Buckley–Leverett equation 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/b/b110/b110990/b1109904.png" /></td> </tr></table>

<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/b/b110/b110990/b1109905.png" /></td> </tr></table>

Under standard conditions, the flux function <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b1109906.png" /> has an  "S"  shape with one inflection point. The equation is a quasi-linear hyperbolic partial differential equation (cf. [[Quasi-linear hyperbolic equations and systems|Quasi-linear hyperbolic equations and systems]] and [[Hyperbolic partial differential equation|Hyperbolic partial differential equation]]). In the case of displacement of oil by water, the initial data are such that they give rise to the Buckley–Leverett profile, consisting of a shock (cf. [[Shock waves, mathematical theory of|Shock waves, mathematical theory of]]) followed by a rarefaction wave. The solution method was simplified by H.J. Welge [[#References|[a2]]], who showed that one could get the shock speed by constructing a particular tangent to the flux function. This technique is equivalent to the Rankine–Hugoniot condition.

The Buckley–Leverett equation is derived under the assumptions of incompressible fluids and immiscible flow. Capillary pressure and gravitation are neglected, and the permeable medium is assumed to obey Darcy's law, stating that the velocity of each phase is proportional to the gradient of its pressure. Let <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b1109907.png" /> denote the porosity of the medium, and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b1109908.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b1109909.png" /> denotes the phase. Conservation of mass of each phase gives that

<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/b/b110/b110990/b11099010.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b11099011.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b11099012.png" /> denote the velocity and density of each phase, respectively. Darcy's law 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/b/b110/b110990/b11099013.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b11099014.png" /> is the (absolute) permeability and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b11099015.png" /> the pressure of phase <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/b/b110/b110990/b11099016.png" />. The pressures are identical in the absence of capillary forces. In one space dimension one can eliminate the pressure, and subsequently scale the time variable, making the total velocity constant, thereby obtaining the Buckley–Leverett equation.

====References====
<table><TR><TD valign="top">[a1]</TD> <TD valign="top">  S.E. Buckley,   M.C. Leverett,   "Mechanism of fluid displacement in sands"  ''Trans. AIME'' , '''146'''  (1942)  pp. 187–196</TD></TR><TR><TD valign="top">[a2]</TD> <TD valign="top">  H.J. Welge,   "A simplified method for computing oil recovery by gas or water drive"  ''Trans. AIME'' , '''195'''  (1952)  pp. 97–108</TD></TR></table>