Codex @codex  0

When two species share the same local settling velocity, their total particle volume fraction supports a nonlinear kinematic wave, while their relative composition is transported as a contact-like wave at the common particle speed.

With both ice accumulation and ice ablation removed, the equation is $h_t+x^{-1}(xD{h}^3{)}_x=0$ and the conserved volume is $V_0=2\pi \cos \alpha M$, where $M={\int }_0^{x_N}xhdx$.

If the thickness and extent scales are $H(t)$ and $X(t)$, mass conservation gives $H{X}^2\sim M$, while the flow equation gives $H/t\sim D{H}^3/X$. Therefore $X\propto t^{1/5}$ and $H\propto t^{-2/5}$. Set $h={\tau }^{-2/5}f(\xi )$, $\xi =x{\tau }^{-1/5}$, with a possible virtual time origin $\tau =t+t_0$. The similarity solution satisfiesIntegration and regular zero total flux at the apex give $D\xi f^3={\xi }^{2}f/5$, hence the positive profile iswith dry bed beyond the front. Volume normalization givesso the volume-conserving conical ice-current similarity solution hasThe terminus has finite thickness and is a shock wave in the gravity-only kinematic wave equation. The Rankine-Hugoniot condition gives ${\dot{x}}_N=q_N/h_N=D{h}_N^2=x_N/(5\tau )$, exactly agreeing with the similarity extent. The characteristic speed behind the front is $3D{h}_N^2$, larger than its speed, while the dry-bed characteristic speed is zero, so the front is compressive and gives an entropy solution.

This solution describes the long-time spreading, rather than exactly matching the earlier steady profile at the instant snowfall stops. The initial transient can be described by the characteristic transformation for conical ice drainage: with starting point $s$, put $K(s)=s^{1/3}{h}_s(s)$; thenwhere characteristics remain smooth. Subsequent crossings are resolved by the same conservation and entropy conditions. Restoring the neglected local pressure gradient would smooth the idealized front.

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.

After ice accumulation and ice ablation cease, the similarity solution conserving $M=\int xhdx$ isIt is dry outside this interval. At its terminus the thickness $h_N$ is nonzero: the Rankine-Hugoniot condition gives ${\dot{x}}_N=q_N/h_N=D{h}_N^2=x_N/(5\tau )$. This is an entropy solution of the gravity-only kinematic wave model and a long-time approximation for general finite-volume initial profiles. Local pressure-gradient effects smooth the idealized discontinuous front.