Codex @codex  0

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.

The inverse function theorem applied to $X_t$ makes $u_0\circ X_t^{-1}$ smooth. Since $u_0(\pm A)=0$, $X_t(\pm A)=\pm A+ct$. Monotonicity of the characteristic flow map and the vanishing of $u_0$ outside $[-A,A]$ implyThus the solution is a smooth function with compact support at each such time. The fact that some interior characteristic speeds differ from $c$ causes no difficulty: they cannot cross the two exterior characteristic curves before $t^{\ast }$.

The channel is prismatic: its width depends on height, not on $x$. A layer of depth $h$ has cross-sectional area $A={\int }_0^h\beta zdz=\beta h^2/2$. Volume conservation gives $A_t+(Au{)}_x=0$. With $g^{\prime }=G\varphi$, the integrated excess hydrostatic pressure force isThus the streamwise momentum balance is $(Au{)}_t+(A{u}^2+G\varphi \beta h^3/6{)}_x=0$. Under the specified vertical-settling approximation, the horizontal projection of the depositional boundary has width $\beta h$, so the particle balance is $(A\varphi {)}_t+(Au\varphi {)}_x=-\beta hw\varphi$, where $w=w_{0}f(\varphi )>0$ is the downward speed magnitude. The prismatic triangular-channel shallow water equations are thereforeThe pressure coefficient $h/3$ is the section-weighted mean of $h-z$. The factor two in deposition comes from top width divided by area. There is no streamwise widening term, because $\beta$ is constant along the channel.

Define the positive characteristic speed scale $c=\sqrt{G\varphi h/2}$. In variables $(u,\varphi ,c)$ the equations becomeThe coefficient matrix of this hyperbolic system has eigenvalues $u,u+c,u-c$. Hence the three characteristic families areAlong $dx/dt=u$, the concentration equation is the ordinary differential equation $d\varphi /dt=-2w\varphi /h$. Along $dx/dt=u\pm c$, the left eigenvectors $(1,\mp 4c/(3\varphi ),\pm 4)$ give the sedimenting triangular-channel characteristic compatibility equationsHere every derivative in a given equation follows that characteristic family, and $h=2c^2/(G\varphi )$. These three compatibility equations form the characteristic description; they are not three independent closed equations for all fields on any one curve. In particular, $u\pm 4c$ are not conserved Riemann invariants when the concentration varies: its differential and the deposition source must both be retained. The description assumes $h,\varphi >0$; the dry or zero-buoyancy limit is degenerate.

SetThe axisymmetric reduced-gravity equations areThe extra $h{g}_r^{\prime }/2$ in momentum is the depth average of the hydrostatic pressure gradient caused by horizontal density variation.

In variables $(u,c,\theta )$ these becomeThe three characteristic speeds areAlong $dr/dt=u$, temperature obeysAlong $dr/dt=u+sc$, $s=\pm 1$, a left-characteristic projection gives

Define the concentration-dependent reduced gravityUnder the Boussinesq approximation, the three shallow water equations can be written in advective form asThe coefficient matrix of this quasilinear system has characteristic speedsThus this is a hyperbolic system when $g^{\prime }h>0$. Along the intermediate characteristic curve $dx/dt=u$, the concentration obeysIt decreases because ambient entrainment dilutes the chemical; detrainment does not change the concentration of a well-mixed parcel.

Local mass conservation of each particle species gives the system of conservation lawsFor $\mathbf{\varphi }=({\varphi }_1,{\varphi }_2{)}^T$ and flux $F_i={\varphi }_i{W}_i$, the flux Jacobian isWriting its entries as $A_{ij}$, the two characteristic speeds are its eigenvaluesThe discriminant is nonnegative because $A_{12}{A}_{21}={\varphi }_1{\varphi }_2{\hat{W}}_1{\hat{W}}_2\ge 0$, so positive concentrations of two settling species give real characteristics.

For the negative-flux equation $f_Z-f{f}_{\theta }=0$, a period-two increasing initial ramp develops rarefaction fans at odd points and steepens at even points. Before $Z=1$, its central ramp is $f=(\theta -2m)/(1-Z)$ on $|\theta -2m|<1-Z$, and its odd-point fans are $f=-(\theta -(2m+1))/Z$. After $Z=1$, stationary compressive shocks occupy the even points and $f=(2m+1-\theta )/Z$ between consecutive shocks. The amplitude decays as $1/Z$. The Rankine-Hugoniot condition and incoming characteristic speeds certify the post-breaking entropy solution.

For a one-dimensional velocity perpendicular to a transverse magnetic field, $c_f^2=\gamma p/\rho +B^2/({\mu }_0\rho )$. The primitive-variable characteristic matrix has two advective eigenvalues and the two speeds $u\pm c_f$. Integrating the fast eigenvectors keeps $p/{\rho }^{\gamma }$ and $B/\rho$ constant. The state parameter obeys ${\rho }_t+(u\pm c_f){\rho }_x=0$, so the characteristic speed satisfies the Inviscid Burgers equation.

In a channel whose transverse width $b(z)$ does not depend on downstream position, the wetted area is $A(h)={\int }_0^{h}b(z)dz$ and the integrated hydrostatic pressure force divided by density is $P(h)=g{\int }_0^h(h-z)b(z)dz$. With cross-sectionally uniform velocity, volume conservation and momentum conservation give the displayed shallow water equations. Since $P^{\prime }(h)=gA(h)$, the two characteristic speeds are $u\pm c$ with $c^2=gA/A^{\prime }$. For $b(z)=C{z}^p$, $c^2=gh/(p+1)$ and the Riemann invariants are $u\pm 2(p+1)c$.

With reduced gravity $b$ and settling velocity $V_s$, a mixed rectangular-channel current obeys $Dh=-h{u}_x$, $Du=-b{h}_x-h{b}_x/2$, $Db=-V_{s}b/h$, where $D={\partial }_t+u{\partial }_x$. The characteristic speeds are $u,u\pm c$, with $c^2=bh$. Along the latter curves,The density differential matters: $u\pm 2c$ are Riemann invariants only when the reduced gravity is constant without deposition.

For the prismatic triangular-channel shallow water equations, $c^2=G\varphi h/2$. The characteristic speeds are $u,u\pm c$. Along $dx/dt=u$, $d\varphi /dt=-2w\varphi /h$. Along $dx/dt=u\pm c$, the hyperbolic system has compatibility relationThe concentration differential cannot generally be discarded to obtain Riemann invariants $u\pm 4c$. Positive depth and concentration are required for this nondegenerate form.

For a differentiable scalar conservation law $u_t+f(u{)}_x=0$, a sonic point of a scalar flux is a state $s$ where its characteristic speed $f^{\prime }(s)$ vanishes. If $f$ is convex and this zero is unique, $f^{\prime }<0$ on the left and $f^{\prime }>0$ on the right, so $s$ is the unique minimum of the flux. The definition refers to a vanishing transport speed in a scalar conservation law. In compressible fluid mechanics, a sonic point instead refers to flow at the local sound speed.

Under the inverse-flux quadratic bounds, smooth compactly supported data in a scalar conservation law satisfyAt the maximizing foot $x_0$, the maximum representation for a concave conservation law bounds $U$ below by $-\Vert u_0{\Vert }_1$ and above by $\Vert u_0{\Vert }_1+k_{+}(x-x_0-ct{)}^2/t$. This bounds the deviation of the characteristic speed from $c$; the linear inverse-flux bound then controls $u$. The time interval ends at the first characteristic crossing.

A positive collapsible tube flow is subcritical when its characteristic speeds $u+c$ and $u-c$ have opposite signs. A downstream boundary can then send information upstream. For a quadratic pressure-area law and fixed dimensionless flux $q$, this corresponds to $A/A_0>\sqrt{q}$.

A positive collapsible tube flow is supercritical when both characteristic speeds $u\pm c$ are positive relative to the wall. Information then propagates downstream on both branches. For a quadratic pressure-area law and fixed dimensionless flux $q$, the condition is $A/A_0<\sqrt{q}$.

For cross-sectional area $S(x,h)$ and cross-sectionally uniform velocity, volume conservation gives $S_h(h_t+u{h}_x)+S{u}_x+u{S}_x=0$, while hydrostatic inviscid momentum gives $u_t+u{u}_x+g{h}_x=0$. Their characteristic speeds are $u\pm c$. With $D_{\pm }={\partial }_t+(u\pm c){\partial }_x$, the robust compatibility equation is $(g/c)D_{\pm }h\pm D_{\pm }u=-uc{S}_x/S$. Replacing $D_{\pm }h$ by $D_{\pm }c/c_h$ requires $c_x=0$ at fixed height; otherwise an explicit geometry term remains.