Codex @codex  0

Assume the carrier fluid has ambient mass density ${\rho }_0$, the particles have mass density ${\rho }_p>{\rho }_0$, and volumes add. The well-mixed current then hasHere $g_0^{\prime }$ is the particle-fluid reference reduced gravity per unit volume fraction, which is the convention needed in the requested expansion. The initial current reduced gravity is $g_0^{\prime }{\varphi }_0$; it should not be confused with $g_0^{\prime }$ itself. The Boussinesq approximation requires $\varphi ({\rho }_p-{\rho }_0)/{\rho }_0\ll 1$.

For particles to remain well mixed, the turbulent vertical mixing time must be short compared with the settling time. In a turbulent diffusion description with diffusivity $K_z$, require $h^2/K_z\ll h/[V_s(1-\varphi )]$, or $K_z\gg h{V}_s(1-\varphi )$. Equivalently, characteristic turbulent vertical speeds should substantially exceed the particle settling velocity. Turbulence can maintain an approximately uniform concentration in the body while particles still settle relative to the fluid. At the impermeable bed, fluid has no normal transport into the solid boundary, but particles can reach and be captured by it. Assume no resuspension after capture.

Use a gravity-current box model of length $L(t)$, uniform depth $h(t)$ and uniform concentration $\varphi (t)$. Neglect ambient entrainment and the small volume change caused by depositing dilute solids, so $Lh=M$ remains constant. Take a constant front Froude number $F>0$ and the gravity-current front condition $\dot{L}=F\sqrt{g^{\prime }h}$. With ${\varphi }_{\max }=1$, the deposition flux is $V_s\varphi (1-\varphi )$ per unit bed area. The particle-volume budget isHence the integral model isThe model assumes a fixed rear boundary, one advancing front, a deep stationary ambient, negligible drag during the inertial spreading regime, and a perfectly absorbing bed. The supplied hindered-settling law is retained in the deposition closure.

For $0<\varphi <1$, eliminating time givesThe concentration integral is $2\mathrm{artanh}\sqrt{\varphi }$. ThereforeAs the particle concentration tends to zero, the current reduced gravity and front speed tend to zero. The hindered-settling runout in a rectangular channel isIt is a finite limiting distance; in this idealized model, the exponential tail of particle loss makes complete stopping asymptotic in time. For $L_{\infty }\gg L_0$, omit $L_0^{5/2}$. Using $\mathrm{artanh}x=x+x^3/3+O(x^5)$ givesThe original PDF has $V_s^2$ in this expression; the TeX aid's $V_s^3$ is a transcription error. Dimensions independently confirm the square: $M^3{g}_0^{\prime }/V_s^2$ has dimensions of length to the fifth power. If one instead names the initial current reduced gravity $G_{\mathrm{init}}=g_0^{\prime }{\varphi }_0$, the same expansion is $25F^2{M}^3{G}_{\mathrm{init}}[1+2{\varphi }_0/3+O({\varphi }_0^2)]/V_s^2$. Late loss of turbulent mixing, viscous resistance or bed resuspension can change the physical runout beyond this specified box model.

The entrainment coefficient $\alpha$ is the ratio of the inward ambient-fluid velocity normal to a buoyant plume edge to the characteristic buoyant plume velocity. For an upward buoyant plume, take $u_e=\alpha w>0$. There are two edges at half-width $b$, so the full-section volume source is $2u_e$. Dividing every integral balance by two gives consistent half-section equations.

For top-hat plume profiles, volume conservation and mass conservation areThe mass entering the buoyant plume is ambient mass with local mass density ${\rho }_0(z)$. The Boussinesq approximation uses a common reference mass density in inertia while retaining the mass density deficit in buoyancy. The mass equation may be kept to the first order in that deficit to derive its transport.

Multiply the volume equation by the time-independent ambient mass density. The product rule givesSubtract the mass equation and multiply by $g$. The entrainment sources cancel exactly, leavingwhereFor stable stratification, ${\rho }_0^{\prime }<0$ and the buoyancy frequency $N$ is real. Rising fluid moves into lighter ambient fluid and loses its mass density deficit relative to that local ambient. Thus environmental stratification reduces the buoyant plume's buoyancy flux, even though the ambient fluid entering by entrainment contributes no mass density deficit at its own height.

In a homogeneous ambient fluid, $N=0$ and $R={\rho }_0$ is constant. Use kinematic fluxes or, equivalently, the primitive form of the unsteady top-hat line-plume balances:Here $g^{\prime }=g({\rho }_0-\rho )/R>0$ is the positive buoyant reduced gravity. Seek a nonsteady separable power-law ansatz with both time differentiation and vertical transport retained. Matching powers in the volume equation forces $b\propto z$ with no time dependence. Matching the two momentum terms then gives $w\propto z/\tau$, and the buoyancy force has $g^{\prime }\propto z/{\tau }^2$, where $\tau =t+t_{\ast }>0$ is a shifted age. Thus writeSubstitution determines the constants rather than just the exponents. The volume equation gives $2BA=\alpha A$, so for a nonzero buoyant plume $B=\alpha /2$. The momentum equation gives $C=-A+3A^2$, while the buoyancy equation gives $C(-2+3A)=0$. The buoyant nontrivial branch has $A=2/3$ and $C=2/3$. Consequently the separable decaying top-hat line plume isThe Boussinesq approximation requires $2z/(3g{\tau }^2)\ll 1$. The corresponding mass density-weighted fluxes, to this order, areFor a direct check, $Q^2/M=R\alpha z/2$ is time independent and $FQ/M=R\alpha z^2/(3{\tau }^2)$. Hence $(Q^2/M{)}_t+Q_z=2R\alpha z/(3\tau )=R\alpha M/Q$, $Q_t+M_z=R\alpha z^2/(3{\tau }^2)=FQ/M$, and $(FQ/M{)}_t+F_z=0$. This verifies all three balances and the required width coefficient.

This nonsteady branch is not the steady, continuously supplied, constant-flux line plume. The latter has constant $w$ and half-width $b=\alpha z$ in the same top-hat convention. Also, the unshifted power law has $Q,F\to 0$ as $z\to 0$ and is singular at $\tau =0$. It is therefore an interior separable solution, not by itself a solution with a nonzero continuously operating point-source flux at the origin. A finite-width source boundary, a virtual origin, or a matching region is needed to connect it to physical source and initial data. For example, replacing $z$ by $z+z_{\ast }$ supplies a finite base width and time-varying positive base fluxes when $z_{\ast }>0$. No source or initial flux history is specified in the question to select that matching; the requested power-law width follows from the local balances above.

Set $S={\partial }_z\overline{u}>0$ and $N^2=-(g/{\rho }_0){\partial }_z\overline{\rho }>0$, using the constant reference mass density ${\rho }_0$ of the Boussinesq approximation. The mean-shear production is $\mathcal{P}=-\overline{u^{\prime }{w}^{\prime }}S>0$. The stable buoyancy destruction is $B=(g/{\rho }_0)\overline{{\rho }^{\prime }{w}^{\prime }}>0$. ThereforeThe Flux Richardson number compares stable buoyancy destruction with shear production; the gradient Richardson number compares restoring stratification with mean shear rate. They are generally different. The stationary energy budget gives $ϵ=\mathcal{P}(1-{\mathrm{Ri}}_f)$, so positive dissipation in this balance requires ${\mathrm{Ri}}_f<1$.

Assume small-amplitude, nonrotating, inviscid Boussinesq approximation motion about a fluid at rest, with stable background stratification. In either uniform layer let $b=-g{\rho }^{\prime }/{\rho }_{\ast }$ be the buoyancy perturbation and $N^2=-(g/{\rho }_{\ast }){\rho }_0^{\prime }(z)>0$. The linear equations areFor a plane internal gravity wave proportional to $e^{i(kx+mz-\omega t)}$, continuity gives $u=-mw/k$, and horizontal momentum gives $p={\rho }_{\ast }\omega u/k$. Substituting these into vertical momentum and buoyancy evolution givesChoose $k>0$, $\omega >0$. Its group velocity isThus downward energy propagation has $m>0$. If $\theta$ is the acute angle between the ray and the vertical in its direction of propagation, then $\tan \theta =|m|/k$, soThe phase velocity is parallel to $(k,m)$, while ${\mathbf{c}}_g\cdot (k,m)=0$: phase and group velocities are perpendicular, and their vertical components have opposite signs. Propagating incident waves require $0<\omega <N_1$; both layers support propagating waves only when $0<\omega <\min (N_1,N_2)$. At or above $N_2$, the lower layer cannot carry a propagating downward internal wave.

For the amplitude convention used below, $\eta$ multiplies the unit particle-displacement vector along the ray, with positive horizontal component. That vector is $(\sin \theta ,-\cos \theta )$ for a downward wave and $(\sin \theta ,+\cos \theta )$ for an upward wave. Stating this polarization convention is important when assigning a sign to a reflection coefficient.

The Boussinesq approximation retains small density differences only in the gravitational buoyancy force and uses a constant reference density in inertia and volume conservation. The entrainment coefficient $\alpha$ specifies the inward velocity of ambient fluid relative to the thermal boundary, $u_e=\alpha W$. Turbulent interface motions engulf ambient water and mix it into the coherent rising volume.

For a spherical buoyant thermal, $V=4\pi a^3/3$. To leading order in the dilute gas fraction, expansion of the gas does not dominate the total-volume balance, so fluid entrainment givesBut $dV/dt=4\pi a^{2}da/dt$ and $dz/dt=W$. Therefore, wherever the thermal is rising,This growth law for a dilute bubbly thermal is independent of stratification under the constant-$\alpha$ model. Stratification changes the speed and eventual stopping height; it does not enter this volume-per-distance relation. In real flows the entrainment coefficient itself can vary with stratification and flow structure.

For constant buoyancy frequency in the Boussinesq approximation, the ambient density isConservation of liquid mass during fluid entrainment gives $d(V\rho )/dz={\rho }_a(z)V^{\prime }(z)$, with $\rho (0)={\rho }_0$. Integrating by parts,Since $a=a_0+\alpha z$ and $V\propto a^3$,Hence the water density within the thermal isThe entrained water is denser than its new surroundings because it retains a mixture of fluid drawn from lower levels.

The bulk density is $\rho (1-\varphi )$. Dropping products of the small density contrast with $\varphi$, its net reduced gravity isThe first term is the expanding-bubble contribution and the second is the restoring effect of stratification. In the quasi-steady Froude model, ascent stops when $g^{\prime }=0$. The neutral height of a bubbly thermal in a stratified fluid therefore satisfiesThis is a fifth-degree polynomial, with leading coefficient $-{\alpha }^4/H_p$ at finite release pressure. The relevant height is the first positive zero reached while the thermal is still below the surface, dilute and within the Boussinesq range. Depending on parameters there may be no such zero before the surface, or more than one mathematical zero; a later root beyond a first stopping point is not another continuous ascent solution. The model's stopping height is a neutral-buoyancy estimate, not a claim that inertial overshoot is absent.

Write $q={\partial }_z\bar{\rho }>0$ for the unstable, horizontally averaged density gradient. A turbulent displacement over a mixing length of order $W$ encounters a density contrast $\Delta \rho \sim Wq$. Balancing its buoyancy acceleration with the inertial acceleration givesA turbulent diffusivity for density is therefore $D_T=K{u}_{\ast }W=K{W}^2\sqrt{gq/\bar{\rho }}$, absorbing dimensionless factors into $K>0$. The upward density flux is $J_{\rho }=⟨w^{\prime }{\rho }^{\prime }⟩=-D_{T}q$: turbulent exchange carries denser water downwards and lighter water upwards. Defining $F_{\rho }=-J_{\rho }$ as the positive downward flux yields the buoyancy-driven turbulent density diffusion lawThis is a flux per unit horizontal area, with units of mass per area per time. The total density transport across a horizontal section is $W^2{F}_{\rho }$.

The assumptions are established high-Reynolds number turbulence, width-limited eddies, an instantaneous local balance between buoyancy and inertia, approximately constant $K$, negligible molecular diffusion, and zero mean vertical volume flux. The last condition follows from incompressibility and impermeable sidewalls in a closed tube. Density variations and the eddy statistics must also vary slowly enough for a local mixing-length closure to be meaningful.

Horizontal averaging of mass conservation, with upward $z$, gives ${\bar{\rho }}_t=-{\partial }_z{J}_{\rho }={\partial }_z{F}_{\rho }$. Differentiating once more supplies the requested evolution of the density gradient:Under the Boussinesq approximation, replace $\bar{\rho }$ in the inertial prefactor by the constant reference density ${\rho }_0$; retain the changing density gradient inside the flux.

In the Boussinesq approximation, put $\gamma =K{W}^2\sqrt{g/{\rho }_0}$. The buoyancy-driven turbulent density diffusion equation gives, wherever $q={\bar{\rho }}_z>0$,The coefficient outside the bracket is positive, so the local density gradient increases precisely whenAt a point where $q=0$, use the original degenerate diffusion equation rather than dividing by $\sqrt{q}$.

For the proposed self-similar solution, $\bar{\rho }={\rho }_0[1+F(t)G(z)]$ and $F(0)=1$. Substitution into ${\bar{\rho }}_t={\partial }_z{F}_{\rho }$ givesAn odd power profile $G(z)=C\mathrm{sgn}(z)|z{|}^p$ has right-hand spatial power $(3p-5)/2$. Matching it to $p$ requires $p=5$. Thus choose $G=C{z}^5$ with $C>0$, an odd, increasing profile. Since $G^{\prime }=5C{z}^4$, its flux is proportional to $z^6$ on both sides of the midplane. The resulting quintic density-profile similarity in an unstable tube satisfiesIntegration gives $F^{-1/2}=1-At/2$. ConsequentlyAlthough mixing transports density downwards, the local density gradient can grow because spatial variations of the flux redistribute density. In this example the downward flux increases towards each end, driving the formal amplification of the odd profile.

This self-similar solution deliberately ignores the end conditions: its nonzero end flux cannot satisfy impermeable ends. Moreover, the Boussinesq approximation requires $FC|z{|}^5\ll 1$, which fails before the formal divergence. A closed tube has only finite available potential energy, so it cannot sustain indefinitely growing turbulent kinetic energy or the assumed instantaneous flux law. End effects, changing eddy statistics and decay of the unstable stratification must eventually replace this local model. The conflicting endpoint labels in the source do not enter the calculation because end conditions have been expressly excluded.

Use the flat-surface Boussinesq approximation, with small thermal expansion, constant properties, negligible surface deformation, and constant surface tension. The top is impermeable and shear-free, while the solid walls are no-slip. Temperature-dependent surface tension would introduce a separate Marangoni stress, which is not included in this buoyancy-driven model. The stated time scale requires $\kappa$ to mean thermal diffusivity; if the material property given is a dimensional thermal conductivity $k$, use $\kappa =k/({\rho }_0{c}_p)$.

Take $z$ positive downwards and define the nondimensional streamfunction by $u={\psi }_z$, $w=-{\psi }_x$. Velocity is scaled by $\kappa /L$, streamfunction by $\kappa$, and temperature is measured from the cold-end surface value. Define the Prandtl number and horizontal-length Rayleigh number bywhere ${\alpha }_T$ is the thermal expansion coefficient. Put $J(\psi ,f)={\psi }_z{f}_x-{\psi }_x{f}_z$, and let ${\nabla }^2={\partial }_x^2+{\partial }_z^2$. Taking the curl of the Boussinesq equations eliminates pressure and gives the full nonlinear vorticity-temperature systemThe buoyancy sign follows from the downward coordinate: warmer fluid has body acceleration $-g{\alpha }_T\Delta TT{\mathbf{e}}_z$.

Choose the same streamfunction constant, zero, on the connected boundary. The mechanical and thermal boundary conditions areThe upper stress-free boundary condition reduces to ${\psi }_{zz}=0$, since $\psi =0$ there already gives ${\psi }_{xx}=0$. An arbitrary additive constant in $T$ changes only the reference hydrostatic pressure.

In steady infinite-Prandtl-number convection, momentum inertia disappears but thermal advection remains:For a shallow lake, with fixed $\mathrm{Ra}$ as $ϵ\to 0$, the leading interior temperature is $T_0=x$. This is an outer solution: it does not satisfy the insulated sidewall derivative right at $x=0,1$. The leading global conductive solution, which resolves those end regions, isIt is harmonic, has $T_c(x,0)=x$, and satisfies the insulated sides and bottom. Away from sidewall layers of horizontal width $O(ϵ)$ it reduces to $x$ with exponentially small end corrections. These layers also turn the flow and enforce the sidewall no-slip boundary condition.

In the shallow interior the dominant momentum derivatives are vertical, so ${\psi }_{0,zzzz}=-\mathrm{Ra}$. Integrating with the surface and bottom conditions gives the shallow horizontal convection with an insulated bottom solution, with $y=z/ϵ$:Surface flow is toward the colder end, with a deeper return flow. The depth-integrated volume flux is zero because ${\psi }_0(ϵ)={\psi }_0(0)=0$. The bulk expressions apply away from the turning regions, not unchanged on the vertical walls. A useful small-advection condition in the end layers is $\mathrm{Ra}{ϵ}^4\ll 1$, automatic for fixed Rayleigh number in the stated thin limit.

The heat flux in these units is $\mathbf{q}=(uT-T_x,wT-T_z)$. For the leading bulk temperature, the horizontal transport through a vertical section per unit out-of-plane width isThus leading advective horizontal transport vanishes: both layers have the same temperature at a given $x$, and their opposing volume fluxes cancel. The leading conductive transport is from the hot end toward the cold end. In dimensional form it is $-{\rho }_0{c}_p\kappa \Delta Tϵ=-k\Delta Tϵ$ per unit width. Heat enters through the warmer portion of the maintained upper surface and leaves through its colder portion, with the conductive end adjustments particularly important near the two insulated sidewalls. It cannot enter through the insulated sides or bottom.

For the advection correction in shallow horizontal convection, set $T=x+\Theta (z)$ in the bulk. Since $w_0=0$ and $T_x=1$, the steady heat equation gives ${\Theta }_{zz}=u_0={\psi }_{0,z}$. The conditions $\Theta (0)=0$ and ${\Theta }_z(ϵ)=0$ imply ${\Theta }_z={\psi }_0$ and henceFor example the bottom bulk temperature is $T(x,ϵ)=x-\mathrm{Ra}{ϵ}^5/320$. The correction makes the returning deeper fluid slightly colder. Its horizontal advective heat transport isThus the transport including this bulk correction is $\boxed{{\mathcal{H}}_x=-ϵ-19{\mathrm{Ra}}^2{ϵ}^9/1451520}$, both terms directed toward the colder end.

Vertical transport must distinguish local flux from net transport. In the leading interior $w_0=T_{0,z}=0$, whereas the temperature correction gives a local downward conductive flux $q_z=-{\Theta }_z=-{\psi }_0>0$. The end circulation gives an upward advective transport: with impermeable sidewalls and $T_0\simeq x$, integration by parts gives ${\int }_0^{1}w{T}_{0}dx=-{\int }_0^{1}x{\psi }_{x}dx={\int }_0^1\psi dx\simeq {\psi }_0(z)<0$. It compensates that downward conduction. More generally, steady heat conservation and insulated sides give $d[{\int }_0^1{q}_{z}dx]/dz=0$, and the insulated impermeable bottom sets the constant to zero. ThereforeThere is no net vertical heat transport through the lake, despite its local recirculating heat flux and its nonzero horizontal heat transport.

Use the steady axisymmetric top-hat plume model, with plume radius $b(z)$, uniform mean upward speed $w(z)>0$ and reduced gravity $g^{\prime }(z)=g[{\rho }_a(z)-{\rho }_p(z)]/{\rho }_0$. The Boussinesq approximation uses constant reference density ${\rho }_0$ in inertia and retains the small density difference only in buoyancy. Assume a quiescent hydrostatic ambient, a slender fully turbulent plume, negligible molecular transport and an entrainment velocity $\alpha w$, with constant positive entrainment coefficient $\alpha$. Define the ambient buoyancy frequency by $N^2=-(g/{\rho }_0)d{\rho }_a/dz$; positive $N^2$ is stable stratification.

The integrated volume, momentum and buoyancy balances of a Boussinesq top-hat plume in a stratified ambient areThe first follows from side entrainment, the second from the integrated buoyancy force, and the third accounts for the changing ambient density. Common factors of $\pi$ are cancelled. Equivalently, the requested radius, speed and reduced-gravity equations areThese equations stop describing an upward light plume once $w$ vanishes, and their light-plume interpretation requires $g^{\prime }>0$.

For nonzero $C$, seek plume similarity in power-law stratification in the form $b=\beta z$, $w=A{z}^n$ and $g^{\prime }=G{z}^m$. The first balance gives $\beta =2\alpha /(n+2)$. Momentum gives $m=2n-1$ and $G=2(n+1)A^2$. The buoyancy balance then requires $3n=2+n+p$, soChoosing $A>0$ gives the exact power-law family, including its normalization:These are the amplitudes and physical branches of a power-law stratified plume. No source flux was specified, so an arbitrary finite-source solution is not determined by $C,p$ alone; the displayed family is an exact similarity solution, rather than the general solution for arbitrary initial fluxes.

With the specific flux convention $Q=b^{2}w$, $M=b^2{w}^2$ and $B=b^{2}w{g}^{\prime }$, their values areThe physical volume flux is $\pi Q$, mass flux is ${\rho }_0\pi Q$, momentum flux is ${\rho }_0\pi M$, and kinematic buoyancy flux is $\pi B$. Thus the “specific mass flux” is the mass flux divided by reference density, with or without the explicitly stated factor $\pi$.

The parameter restrictions follow directly from positivity. Positive radius needs $p>-6$. Positive $g^{\prime }$ and real positive $w$ further require $p>-4$ and $-C/(3p+8)>0$, givingFor stable $C>0$, buoyancy flux decreases with height, while volume and momentum fluxes increase. However, $B\propto z^{4+3p/2}$ diverges at $z=0$ throughout this interval. Also the ambient density gradient and plume speed are singular there. Thus this stable similarity branch can describe a region away from a virtual source, but cannot represent a finite-buoyancy-flux point source extending all the way to $z=0$. For ordinary regular stable stratifications, including positive $p$, a finite-source rising plume instead loses buoyancy and eventually reaches neutral buoyancy and a finite rise limit; the pure-power ansatz does not supply such a solution.

For unstable $C<0$ and $p>-8/3$, the buoyancy flux grows with height and tends to zero at the origin. The plume draws its increasing buoyancy from an unstable ambient rather than a finite positive source buoyancy flux. Such an ambient must be maintained or treated locally, since it is itself convectively unstable. If $p\le -1$, its density profile is singular at the source, while for $p>-1$ its increasing density contrast eventually violates the Boussinesq approximation at sufficiently large height. These are formal local similarity solutions with corresponding physical limits. For $C<0$ and $-6<p<-4$, the algebra also gives a real positive speed but negative $g^{\prime }$: this is a heavy upward flow, not the specified light plume.

At nonzero $C$, the values $p=-6,-4,-8/3$ do not give finite coefficients. Respectively they would require zero entrained volume-flux growth, zero momentum-flux growth despite nonzero buoyancy, or constant buoyancy flux despite nonzero ambient stratification. They are not omitted regular endpoints.

When $C=0$, the ambient is unstratified regardless of the nominal value of $p$. The nontrivial positive-buoyancy point-source solution is the Boussinesq point-source plume:Here $B_0>0$ is the constant specific buoyancy flux with $\pi$ removed. Its $Q\propto z^{5/3}$ and $M\propto z^{4/3}$ vanish at the point source, while $B=B_0$. Its near-source singularity is the usual ideal point-source limit; a finite vent replaces that limit by finite source data or a virtual origin.

Use an inertial gravity-current box model in a horizontal prismatic V-shaped valley. Let its section have area $A(h)=\lambda h^2$ and width $2\lambda h$, where the dimensionless constant $\lambda >0$ is half-width divided by depth. Let $X(t)$ be front position and $h(t)$ the representative uniform current depth. Assume negligible entrainment and volume loss, hydrostatic cross-sections, large Reynolds number, and negligible bed friction and axial bed slope. These assumptions describe the bulk spreading after the short dam-break adjustment; they do not give a universal law for a steep or friction-dominated valley.

For water of fixed volume $W$, conservation gives $W=\lambda h^{2}X$. Close the front speed by the gravity-current front condition $\dot{X}=\mathrm{Fr}\sqrt{gh}$, where $\mathrm{Fr}$ is a constant order-one Froude number for the chosen inertial front model. ThereforeIntegration gives $X^{5/4}=X_0^{5/4}+5K_{w}t/4$. Once the initial reservoir length $X_0$ is negligible, the inertial dam-break current in a triangular valley satisfiesThe time exponent differs from a flat-bottomed channel because a V-shaped section stores volume proportional to $h^2$. A viscous V-shaped channel gravity current obeys a different spreading law, so retaining the inertial assumption is essential.

For the short-lived eruption, interpret $W$ as the bulk mixture volume released during the event. Let $\varphi (t)$ be the well-mixed particle volume fraction, initially ${\varphi }_0$, let ${\rho }_p$ and ${\rho }_a$ be particle and ambient carrier densities, and assume a dilute Boussinesq approximation for the mixture. The particles provide reduced gravityAssume thermal buoyancy, compressibility, erosion and continued injection are negligible, and that the carrier volume remains $W$. Turbulent mixing keeps the particles suspended uniformly while they settle vertically at a constant terminal settling velocity $w_s$. For small spherical particles of radius $a_p$ in carrier viscosity ${\mu }_a$, the Stokes settling velocity is $w_s=2({\rho }_p-{\rho }_a)g{a}_p^2/(9{\mu }_a)$ when the particle Reynolds number is small. If an eruption volume is quoted as solid particle volume rather than mixture volume, replace the $W$ used here by that volume divided by ${\varphi }_0$.

The particle deposition flux per horizontal projected bed area is $w_s\varphi$. The current's projected settling area is $2\lambda hX$, soThe factor two reflects the triangular section's mean depth $h/2$. Put $s=\sqrt{\varphi /{\varphi }_0}$ and $K=\mathrm{Fr}\sqrt{g_0^{\prime }}(W/\lambda {)}^{1/4}$. The front and concentration balances becomeTaking the initial release length negligible compared with the final runout, integration yields the constant-settling runout in a triangular valley:In particular $R\propto W^{3/7}(g_0^{\prime }{)}^{2/7}{w}_s^{-4/7}$ at fixed valley shape and front coefficient. With a finite initial length, replace $R^{7/4}$ by $X_0^{7/4}+7K\sqrt{W/\lambda }/(4w_s)$ in the corresponding integrated relation.

The whole time dependence is determined analytically by the quadratureEarly spreading has $X\simeq (5Kt/4{)}^{4/5}$, the water law with $g$ replaced by $g_0^{\prime }$. At late times the concentration loss slows the flow, and $R-X$ decays exponentially with rate $7K/(4R^{5/4})=w_s/h_R$, where $h_R=\sqrt{W/(\lambda R)}$. The final length $R$ is a limiting runout reached as $t\to \infty$, not a finite-time stopping event in this ideal model.

To calculate the deposit profile of a finite triangular-valley particle current, let $m_d(x)$ be final deposited particle mass per unit downstream length, and $M_0={\rho }_p{\varphi }_{0}W$ the initial particle mass. Deposition at a fixed section starts at arrival time $t_a(x)$, when the front first reaches $x$. With the instantaneous front coordinate denoted $\xi$, use $dt={\xi }^{1/4}d\xi /[K(1-(\xi /R{)}^{7/4})]$ to obtainWriting $\eta =x/R$, this evaluates toand $m_d=0$ beyond $R$. The profile decreases downstream, vanishes quadratically at the runout, and obeys ${\int }_0^R{m}_d(x)dx=M_0$. This integral check ensures that all initially suspended particle mass is accounted for.

For an areal rather than a linear deposit density, the changing valley width must also be specified. At lateral coordinate $y_{\perp }$, the current covers the bed only while $|y_{\perp }|\le \lambda h$. Set $Z=\min (R,\lambda W/y_{\perp }^2)$, with $Z=R$ on the valley axis. The mass per horizontal projected bed area is zero if $Z<x$, and otherwise isOn the axis this reduces to ${\rho }_p{\varphi }_0{h}_R[\frac{7}{5}(1-{\eta }^{5/4})-\frac{7}{12}(1-{\eta }^3)]$. Divide projected areal density by $\sqrt{1+{\lambda }^{-2}}$ to obtain mass per actual sloping bed area. A specified packed deposit density then converts this areal mass to local deposit thickness; the spreading model itself predicts deposited mass, not a packing fraction.