Past exam of the mathematics course of the University of Cambridge / 2024 / iii / Paper 345 / 1

Let $W_f$ be the vertical velocity of the carrier fluid. Each species has its isolated Stokes settling velocity relative to that fluid, so $W_i=W_f+{\hat{W}}_i$. The mixture has no imposed vertical volumetric flow rate. Using the particle volume fractions, its zero-volume-flux condition isSubstitution of $W_i=W_f+{\hat{W}}_i$ gives the common carrier-fluid backflowConsequently the hindered settling velocities areThe same $W_f$ appears for both populations precisely because every particle is assumed to feel the same volume-averaged backflow.

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.

Setwhere $s>0$. Applying first-order perturbation of a simple eigenvalue to the flux Jacobian, or expanding its quadratic formula directly, separates the characteristic that changes total concentration from the characteristic that changes composition:At $ϵ=0$, both species move with the common hindered velocity $\hat{W}(1-s)$. Summing their conservation laws givesso the total concentration is a nonlinear kinematic wave:Taking the ratio $c={\varphi }_1/s$ instead gives the composition wave in a bidisperse suspensionand henceThese are the requested leading-order ordinary differential equations along the two characteristic families.

Let the fixed current volume per unit channel width beA suitable high-Reynolds number deep-ambient gravity-current front condition iswhere the order-one Froude number $\mathrm{Fr}$ records the selected front closure. In the dilute limit the mixture's reduced gravity iswith ambient and carrier-fluid density ${\rho }_a$.

The well-mixed particle volume of species $i$ is $V_0{\varphi }_i$. Its deposition rate through the base of length $L$ is $L{\hat{W}}_i{\varphi }_i$, where ${\hat{W}}_i<0$. The resulting gravity-current box model is thereforeIt conserves fluid volume while suspended particle volume, and therefore the driving buoyancy, decreases by deposition.

When ${\hat{W}}_1={\hat{W}}_2=\hat{W}$, both concentrations acquire the same decay factor. The total reduced gravity consequently obeysEliminating time with the front equation givesIntegration from the initial state yieldswhereThe particle-laden gravity current reaches its runout length of a gravity current when $g^{\prime }$ vanishes. Since $\hat{W}<0$,Thus $C=5\mathrm{Fr}$ for the stated front condition; the common normalization $\mathrm{Fr}=1$ gives $C=5$.

Because $\hat{W}<0$, species 1 now settles more slowly and species 2 more rapidly than in the equal-speed case. The supplied inequality says that species 2 provides the larger initial particle-density contribution to the reduced gravity; with exact Boussinesq density contrasts, the corresponding condition is ${\tilde{\varphi }}_1({\rho }_1-{\rho }_a)<{\tilde{\varphi }}_2({\rho }_2-{\rho }_a)$. The dominant buoyancy contribution is therefore removed earlier, so $g^{\prime }(t)$ and the gravity-current front condition fall below their equal-settling values at first order in $ϵ$. The runout length decreases. The slower loss of the weaker species-1 contribution only partly compensates for this effect.