Codex @codex  0

The fraction of total advective Darcy flux carried by one phase when capillary and gravity corrections are omitted. For two phases it is ${\lambda }_w/({\lambda }_w+{\lambda }_n)$. This quantity lies between zero and one for nonnegative phase mobilities. Its nonlinear dependence on fluid saturation determines characteristic speeds and shocks in the Buckley-Leverett equation.

Saturation equation. Introduce intrinsic permeability of a porous medium $K$ and porosity $\varphi$. The relative porous permeabilities depend on the wetting fluid saturation. Define the phase mobilities and fractional flowUse the conventional capillary pressure $p_c=p_{nw}-p_w$, generally decreasing as wetting saturation increases. Neglect gravity in this horizontal model. The two Darcy fluxes are $q_w=-K{\lambda }_w{p}_{w,x}$ and $q_n=-K{\lambda }_n(p_{w,x}+p_c^{\prime }(s)s_x)$. Eliminate the common pressure gradient using $q_w+q_n=Q$ to obtainThe wetting-phase mass conservation law is consequentlyFor $p_c^{\prime }\le 0$, $D_c\ge 0$ and the capillary pressure provides nonlinear diffusion. The initially uniform saturation is $s_0$. Pure wetting-fluid injection prescribes $q_w(0,t)=Q$ and $q_n(0,t)=0$; in the zero-capillarity, zero-residual-nonwetting idealization this corresponds to inlet saturation one. With residual nonwetting fluid, replace one by the appropriate maximum accessible saturation.

Shock formation and its speed. Neglect capillarity first. Smooth saturation values travel along characteristic curves at speed $(Q/\varphi )F^{\prime }(s)$. When the trailing values have larger characteristic speeds than the values ahead, the curves intersect, and the Buckley-Leverett equation requires a shock selected as an entropy solution. This is typical for the increasing-convex part of a physical fractional-flow curve; shocks are not inevitable for every possible $F$.

Integrating mass conservation across a shock from upstream $s_s$ to downstream $s_0$ gives the Rankine-Hugoniot conditionIf an upstream rarefaction wave joins the shock smoothly in similarity coordinates, its terminal characteristic speed equals $V_s$. Therefore the saturation at that join obeys the fractional-flow tangent constructionThis is the printed equality. It is a tangent-selection condition, not the dimensional speed by itself; the speed includes $Q/\varphi$. An arbitrary shock between prescribed constant states satisfies the jump condition but need not satisfy this additional tangency relation.

Positive capillary diffusion spreads a jump into a transition layer. For a travelling profile $s(x-V_{s}t)$, integration givesThe endpoint states still satisfy the same Rankine-Hugoniot condition. Thus small capillarity primarily gives the front a finite thickness rather than changing the selected limiting speed; vanishing phase mobility at an endpoint can make the regularization degenerate.

The specified reciprocal fractional flow. On its stated branch, $F^{\prime }=1/s^2$ and $F^{\prime \prime }=-2/s^3<0$. If $F$ extends continuously to the initial saturation and the inlet is at one, the decreasing saturation from behind to ahead makes characteristic speeds increase forwards: there is a rarefaction wave, not a compressive shock. Write $a=Q/\varphi$ and $\xi =x/t$. The self-similar solution isThe fan connects continuously to the initial state and broadens linearly in time. This interpretation requires $F(s_0)=2-1/s_0$. A nonnegative physical fractional flow on the entire interval $[s_0,1]$ also requires $s_0\ge 1/2$.

The PDF specifies the formula only for $s>s_0$ and earlier calls $s_0$ residual. If that means an immobile wetting phase with $F(s_0)=0$, the endpoint is an extra constraint and the continuous-branch answer applies without conflict only at $s_0=1/2$. For $0<s_0<1/2$, the printed branch would give negative fractional flow for $s_0<s<1/2$, so it cannot be the complete physical constitutive law. With an additional admissible extension, for example $F=0$ up to $1/2$ and the given positive branch above it, a compound fan–shock is possible. The tangent condition then givesThe fan $s=\sqrt{a/\xi }$ ends at $\xi =V_s$, followed by a jump to $s_0$. This is a conditional physical extension, not data supplied by the question. For $s_0>1/2$ an immobile residual endpoint is likewise incompatible with a continuous version of the stated branch. These distinctions identify what is, and is not, determined by the printed constitutive assumption.