Codex @codex  0

The matrix $C$ is antisymmetric, so the scalar Majorana spinor bilinear is not identically zero. For example, choose two independent odd Grassmann variables ${\lambda }_{\alpha },{\lambda }_{\beta }$ at a pair of indices with $C_{\alpha \beta }\ne 0$. Their contribution is $2C_{\alpha \beta }{\lambda }_{\alpha }{\lambda }_{\beta }$, which is nonzero as a Grassmann polynomial. This is why a Majorana mass term is possible.

A classical sequential probability ratio test addresses a fixed-start choice between two simple statistical hypotheses. It sums log-likelihood ratio increments until an upper boundary accepts the alternative or a lower boundary accepts the null; it then terminates that test. A CUSUM instead monitors indefinitely for a change at an unknown time. Its reset at zero discards sustained evidence from a previously satisfactory period, so a late change need not first undo all earlier negative scores. An upper boundary $h$ signals when the accumulated change evidence is sufficiently large.

A fast initial response CUSUM takes $0<X_0<h$, often $X_0=h/2$, in place of a zero start. Until its first reset or signal, it acts like a fixed-start sequential probability ratio test with a lower boundary at $-X_0$ and an upper boundary at $h-X_0$ for the subsequent cumulative increments. Upon reaching the lower boundary it resets and continues as an ordinary CUSUM. The head start gives faster detection if a change is already present when monitoring begins or restarts. It is useful after an intervention, a shutdown or a previous alarm when immediate renewed surveillance is desired. Its initial false-alarm behavior differs from the zero-start scheme, so the head start and signalling boundary must be calibrated together; the preceding numerical choice is a design convention, not a universal optimum.

Let $p_1$ be the alternative death probability. The odds ratio specifiesAssume $0<p_0<1$ and $k>0$. For observed batch count $y_t$, the common combinatorial factor in the two binomial likelihoods cancels, giving the CUSUM incrementThereforeFor $k>1$ large death counts add positive evidence of deterioration; for $k<1$ unusually small death counts add evidence for improvement. The score has nonpositive expectation under the null and nonnegative expectation under the specified alternative, as follows from the nonnegativity of Kullback-Leibler divergence. Independent, nonoverlapping batches avoid counting the same patient's outcome repeatedly.

The same pricing relation gives $x^{\top }{A}_0\ge 0$ whenever the terminal payment is nonnegative. Therefore negative initial cost also contradicts a strictly positive state-price density, even if the terminal payment is identically zero.

This second form is needed for directions invisible in the terminal payments. Let $L$ be the linear span of the essential values of $A_1$, equivalentlyIf $A_0\notin L$, let $v$ be its nonzero orthogonal projection on $L^{\perp }$ and set $x=-v$. Then $x^{\top }{A}_1=0$ almost surely while $x^{\top }{A}_0=-\Vert v{\Vert }^2<0$. Hence the second form exists in this case. Part (b) handles all remaining cases.

It is helpful to solve using the nonsingular total flux $J$ and volume flux $Q$ before recovering $M$. The relation $J=M/(1-Q^2/(R^{2}M))$ is quadratic:The plus sign selects the small-area pure-source branch. Its discriminant vanishes when $J=4Q^2/R^2$, givingThereforeThe specific upward plume and ambient momentum-flux contributions are equal, not opposite. The ambient's transported momentum has the opposite sign from its transported volume, so their product remains positive.

To show that the pure branch actually reaches this fold in finite height and that it is a genuine pole, introduceThen $s=q^2/m$ andThe pure-source condition is $j\sim (5/8{)}^{2/5}{q}^{4/5}$. Before the fold, $s<1/2$ and $m>2q^2$. Thus for every fixed positive $q_1$, the derivative $dj/dq$ is bounded above by a constant times $q^{-2}$ on $q\ge q_1$. If the branch continued for arbitrarily large $q$, $j$ would remain bounded, contradicting $j\ge 4q^2$. Hence it reaches $j=4q^2$ at a finite $q_o>0$. Its height is finite because $d\zeta /dq=(1-s)/(2\sqrt{m})$ is finite away from the source and behaves like a constant times $q^{-2/5}$ at the source, an integrable singularity.

For the numerator, write $m=j(1-s)$. ThenAt the fold $j=4q_o^2$, so $q_o^3\ge 5/256$. On the other hand the dimensionless numerator in the $m$ equation at $s=1/2$ isIt is strictly negative because $5/256>1/(64\sqrt{2})$. The denominator $1-2s$ approaches zero from above. Consequently $dM/dz\to -\infty$ at finite $z_o$, although the fluxes and velocities themselves remain finite. This proves the confined-plume momentum fold rather than merely noticing a possible zero denominator.

The dimensionless pure-source branch uniquely determines the remaining amplitudes. They can be specified without an elementary antiderivative by integrating the displayed scalar equation for $j(q)$ until $j=4q^2$, and settingA direct numerical integration gives $q_o\simeq 0.281610$, ${\zeta }_o\simeq 0.400049$, so the dimensional results for this specified entrainment convention areA convention using $\alpha w$ instead of $\alpha (w-v)$ would produce different numerical amplitudes; the relative-speed convention was stated in part (a).

Physically $z_o$ is a turnover or breakdown height for a slender steady plume surrounded by uniform counterflow. The opposing streams then occupy equal areas and have equal speeds; the single-valued small-area branch cannot continue smoothly. Above it one expects lateral spreading, stronger recirculation, mixing and unsteady buoyant structures rather than a continuation of the same slender top-hat column. A detailed flow or upper-layer evolution requires restoring pressure and stress dynamics and, eventually, the modified ambient stratification. The model does not predict an infinite physical velocity, nor does it require that all buoyant fluid forever remain below $z_o$.

Normalize all fluxes by $\pi$, rather than including that factor in only some of them. The specific plume volume flux, momentum flux and buoyancy flux areThe ambient transport is $-Q$ and its specific momentum flux isPut $s=b^2/R^2=Q^2/(R^{2}M)$, with $0<s<1$. The total specific momentum flux isThe first differential equation follows immediately from relative-speed fluid entrainment:For the second, differentiate $J=M/(1-s)$ and use $s^{\prime }=2s{Q}^{\prime }/Q-s{M}^{\prime }/M$. Since $J^{\prime }=b^2{g}^{\prime }=B_{s}Q/M$, this yieldsThis is the requested two-equation system, with $B=B_s$ and zero total volume flux as the additional conservation conditions. Near the pure source, $Q,M\to 0$ and $s\to 0$. Eliminating height in that limit givesThe zero integration constant selects a pure plume; a nonzero imposed source momentum would select a different branch and generally a different turning height. The leading radius is $b\sim 6\alpha z/5$, which recovers the unconfined pure-plume limit near the source.

Use the Boussinesq approximation and a top-hat plume model. Let $v(z)<0$ be the uniform ambient return velocity outside the plume. Negligible source volume flux and an impermeable cylindrical wall require zero net transport through every horizontal section:The entrainment speed must be specified relative to that moving ambient. Taking it to be $\alpha (w-v)$ givesThe surrounding unmodified fluid has zero buoyancy relative to the initial fluid. Neglecting thermal or compositional loss therefore givesLet $B_s$ denote the conserved specific source buoyancy flux: $B_s=B_o/\pi$ if $B_o$ is the actual integrated source flux, or $B_s=B_o$ if the source strength has already been normalized by $\pi$.

To obtain the usual integral-model singularity one must also state a pressure closure. Here take the whole-section mean pressure to have only its reference ambient hydrostatic pressure, and neglect wall force and extra vertical turbulent-stress fluxes in the integral budget. Internal momentum transfers between plume and ambient cancel in that whole-section balance. The remaining upward force is the plume buoyancy, givingThe two momentum fluxes add, even though $v$ is negative: vertical momentum transported downward also contributes $v^2$ to the flux through an upward-oriented horizontal section. These three flux balances, with the algebraic return-velocity constraint, govern $b,w,g^{\prime }$ in the confined plume with compensating return flow.

The pressure-neglecting closure is an integral approximation, not a claim that the return fluid has zero acceleration. If a mean pressure departure $p_{\ast }(z)$ is retained, the whole-section momentum balance instead contains $-(R^2/{\rho }_0)p_{\ast }^{\prime }(z)$, along with any wall-force or stress contribution. The entrainment hypothesis alone does not specify that function. The singularity calculated below belongs to the stated pressure-neglecting model; it is not forced by continuity alone or asserted to be an infinite velocity in the real fluid.

The filling box model requires the first outflow to spread across the enclosure and remain above unmodified fluid, while a narrow nearly vertical plume rises below it. Several mechanisms can defeat those assumptions. A tall narrow box makes the plume a substantial fraction of the width, so its compensating return velocity is no longer small and modifies both fluid entrainment and momentum flux. A fast ceiling current can descend the sidewalls and overturn the initially stable upper region. In both cases a horizontally uniform slowly descending filling-box first front is inappropriate. Strong plume displacement from the centre, imposed crossflow, large source volume, vents or substantial wall cooling likewise alter the volume and heat budgets.

A useful diagnostic is whether the ceiling outflow's inertial acceleration, of order $U_c^2/{\ell }_c$, overwhelms the restoring buoyancy of the nascent layer. If $U_c^2/(G_H{d}_c)$ is not small for its depth $d_c$, the outflow can penetrate or overturn the stratification rather than remain trapped above it. This is a diagnostic Froude number, not a universal threshold without specifying the wall-current geometry and mixing. Low Reynolds number invalidates a constant turbulent entrainment coefficient, and thermal diffusion or wall heat loss can erase the sharp first-front jump. Thus large aspect ratio, strong return circulation, overturning, ventilation or diffusion require a model beyond the standard filling-box approximation.

Let the closed box have ceiling height $H$ above the source. The printed width $2R$ alone does not specify this height, so it must appear as an additional geometrical parameter. Take $t=0$ here when the plume has first reached the ceiling and its thin outflow has spread across the box. The leading filling box model assumes a negligible source volume, insulated boundaries, fast horizontal spreading, slow evolution compared with the plume transit time, negligible diffusion, and $b\ll R$. Below the interface the ambient still has its original temperature, so the local plume flux remains $q(z)=2\alpha w_{\ast }z$ from part (c).

The ceiling outflow creates a warm upper region. Fluid entrainment draws the original lower fluid into the plume, making that upper region deepen. The filling-box first front is the lower material boundary of fluid already modified by plume discharge, separating it from the initially unmodified ambient; it need not bound a uniformly mixed upper layer. With negligible occupied plume area, volume conservation at its height $h(t)$ givesThereforeIf source switch-on is the time origin instead, replace $t$ by elapsed time after the initial transit and spreading stage; the quasi-steady model does not resolve that stage.

The buoyancy jump requires distinguishing ambient fluid at the front from the plume crossing it. Just after the first ceiling outflow, the ambient buoyancy immediately above the nascent front isOutside the narrow plume, ambient temperature is transported without diffusion. If $G_a$ is its buoyancy relative to the original fluid, its equation isThe first front follows precisely these material trajectories, so the temperature of the fluid immediately above it is conserved. The initially unmodified fluid below has $G_a=0$. Thus the ambient jump isBy contrast, the contemporaneous plume immediately below the front has $G_p(h,t)=w_{\ast }^2/h(t)$, which grows as the front descends. That plume value is not the ambient jump across the material first front. Equating the two would inadvertently impose a uniformly mixed upper layer and discard the stable density stratification of the standard filling box model. Diffusion or extensive interfacial mixing can smooth and alter the ideal jump.

Let ${\mathcal{B}}_0>0$ be the source buoyancy flux per unit length. The conversion from the heat-source strength depends on its units: if $F$ is a kinematic temperature-volume flux, ${\mathcal{B}}_0=g\beta F$; if it is a physical heat power per unit length, ${\mathcal{B}}_0=g\beta F/({\rho }_0{c}_p)$, where $c_p$ is specific heat capacity at constant pressure. This explicit convention avoids treating heat power as a buoyancy flux without conversion.

For uniform ambient $N^2=0$, so $\mathcal{B}={\mathcal{B}}_0$. The ideal pure line-source limit has negligible initial volume and momentum transport. The flux equations admit $w=w_{\ast }$ constant, $q=2\alpha w_{\ast }z$ and $b=\alpha z$. Substituting into the momentum equation gives $2\alpha w_{\ast }^2=2bG$, and the conserved buoyancy flux then fixes $2\alpha w_{\ast }^3={\mathcal{B}}_0$. ThusThis verifies all three differential equations, not just their dimensional exponents. The plume grows in width while maintaining constant upward velocity; its temperature anomaly decreases like $z^{-1}$. The zero-height singularity is the ideal point-source limit. A finite source or nonzero momentum input introduces a plume virtual origin or a forced near-source region; the heat strength alone cannot uniquely determine that finite-source region.

In the top-hat plume model, vertical velocity $w(z)$ and reduced gravity $G(z)$ are uniform over $-b<x<b$. The ambient has zero vertical velocity, negligible ambient turbulence and a prescribed $T_0(z)$. Thenand the equations becomeEquivalently, the ordinary differential equations for the three unknowns areThe last equation displays dilution by fluid entrainment as well as loss of relative buoyancy while rising through a stably stratified ambient. The model assumes a steady slender two-dimensional turbulent plume, large Reynolds number, the Boussinesq approximation, constant $\beta$ and $\alpha$, negligible heat loss, negligible imposed crossflow, and negligible pressure departure from ambient hydrostatic pressure. The factors of two count both entraining edges; $b$ is a half-width. A finite source also requires its volume, momentum and buoyancy fluxes as boundary data.

Use $G=g\beta (T-T_0(z))$ for the plume reduced gravity, and take $z$ upwards. For a symmetric line plume define fluxes per unit source length,The integrals span the plume; its edges can move horizontally with $z$. The fluid entrainment rate at one moving edge is the inward normal transport relative to that edge. For an edge $x=b(z)$ it is $E=w{b}^{\prime }-u$; symmetry supplies equal transport through the other edge. The Batchelor entrainment hypothesis closes its mean value as $E=\alpha w_c$, where $w_c$ is a representative vertical velocity and $\alpha$ is the entrainment coefficient. The integrated volume conservation equation is therefore $q^{\prime }=2\alpha w_c$.

Using volume conservation, the vertical acceleration equation has conservative form ${\partial }_x(uw)+{\partial }_z(w^2)=G$. Incoming ambient has zero vertical momentum; the integrated momentum conservation equation is consequently $m^{\prime }=\int Gdx$. For the temperature anomaly, $T_0$ depends on height andAmbient fluid enters with $G=0$, giving the buoyancy flux balance. The integral model isHere $N$ is the buoyancy frequency when $N^2>0$; the same algebra permits unstable temperature gradients with $N^2<0$.

These are mean integral balances. The advective equations printed in the question omit explicit turbulent stresses and heat fluxes; a turbulent plume interpretation requires the usual entrainment closure for that unresolved transverse transport. Assuming vanishing edge stresses and heat fluxes, and negligible additional axial turbulent-flux contributions, gives the displayed bulk model. It would be inconsistent to treat a discontinuous top-hat profile as a smooth exact pointwise solution of the inviscid temperature-advection equation.

With $W=0$ and $V>0$, $Q=Q_0$ and $\varphi =e^{-kx}$, where $k=V/Q_0$. No carrier momentum is withdrawn; hence the integrated momentum flux is exactly constant:On an initially small-Froude number branch, the pressure term dominates. At leading orderA slightly more accurate large-distance amplitude is $\sqrt{2K/g_0^{\prime }}$, so $h\sim \sqrt{2K/g_0^{\prime }}{e}^{kx/2}$. The Froude number satisfies ${\mathrm{Fr}}^2\sim {\mathrm{Fr}}_0^2{e}^{-kx/2}$ and remains small. The growing depth can eventually violate confinement or the shallow water assumption; this does not make the asymptotic momentum calculation inconsistent within its intended regime.

On an initially large-Froude number branch, horizontal momentum flux dominates instead. The leading depth is constant, $\boxed{h\sim h_0}$, while the decaying reduced gravity makes the Froude number still larger. To resolve the small change, put $\delta ={\mathrm{Fr}}_0^{-2}\ll 1$ and $y=h/h_0$:soThe depth decreases slightly and tends exactly to $h_{\infty }=Q_0^2/K=h_0/(1+\delta /2)$. Its Froude number grows like ${\mathrm{Fr}}_0{e}^{kx/2}$. Both approximations agree with the exact depth equation $h^{\prime }=kh/[2(1-{\mathrm{Fr}}^2)]$: the subcritical flow deepens, while the supercritical flow thins.

At $V=0$, the particle budget gives $\varphi =1$, so $g^{\prime }=g_0^{\prime }$ is constant. The primitive momentum conservation equation integrates toIntroduce $q=Q/Q_0=1-Wx/Q_0$, $y=h/h_0$ and $\epsilon ={\mathrm{Fr}}_0^2=Q_0^2/(g_0^{\prime }{h}_0^3)\ll 1$. The exact relation on the subcritical flow branch isExpanding uniformly for $0\le q\le 1$ givesThe depth increases slightly because horizontal kinetic energy decreases as fluid is removed. In fact the depth equation has $h^{\prime }\ge 0$ on this branch, andThus the small-Froude number approximation becomes better downstream, including as the flux vanishes. At leading order $u=(Q_0-Wx)/h_0$.

Let $W\ge 0$ be the downward fluid suction speed and $V\ge 0$ the downward particle settling velocity relative to that fluid. Set $Q=uh$ and $g^{\prime }=g_0^{\prime }\varphi$. The carrier volume conservation and particle-volume balances per unit channel width areThe particle sink includes suction as well as relative settling: the particles follow the well-mixed withdrawn fluid, in addition to settling through it. Subtracting $\varphi Q^{\prime }$ yields $Q{\varphi }^{\prime }=-V\varphi$. Thus for $W>0$, as long as $Q>0$,For $W=0$ the continuous limit is $Q=Q_0$ and $\varphi =e^{-Vx/Q_0}$. In particular suction alone does not change the well-mixed particle fraction: $V=0$ implies $\varphi =1$.

Withdrawn fluid carries its local horizontal momentum, so the momentum flux and hydrostatic pressure force obeyExpand this and use $Q^{\prime }=-W$. The suction momentum terms cancel in the primitive acceleration equation, leavingThese are the governing equations for a particle-laden current with suction. Substituting $u=Q/h$ and the concentration equation gives the convenient depth equationThe steady carrier flux is exhausted at the formal length $\boxed{x_{\ast }=Q_0/W}$ when $W>0$. Positive $Q$ cannot continue beyond it. Near that limit a chosen depth branch, mixing or shallow-layer assumptions may fail, so this is the mass-budget extent, not a guarantee that every arbitrary inlet state has a regular shallow continuation to it. When $W=0$, settling makes the driving reduced gravity exponentially small but does not remove carrier fluid; this ideal steady model has no finite mass-budget runout.

Use a gravity-current box model: approximate the current by uniform depth $h(t)$ and length $L(t)$, with constant reduced gravity, no deposition and no fluid entrainment. The triangular cross-sectional area is $h^2/2$, so the fixed released volume isThe Benjamin deep-ambient front condition closes the integral model:For a finite initial length $L_0$, direct integration givesThe ideal point-release limit takes $L_0=0$ and has $L\propto t^{4/5}$, $h\propto t^{-2/5}$. Its divergent initial depth is outside the shallow water approximation, so a real finite reservoir regularizes the early stage. This is the finite-volume triangular-channel current integral law; it does not assert a spatially uniform exact solution of the local equations.

Put $c_0=\sqrt{g^{\prime }{h}_0/2}$. The undisturbed reservoir has $u=0$ and Riemann invariants $u\pm 4c=\pm 4c_0$. The right-going characteristic curves supply the invariant $u+4c=4c_0$ throughout the release. The spreading disturbance is the other family, whose speed is $u-c$. In its rarefaction wave let $\xi =x/t=u-c$. ThenAt the reservoir edge $u=0,c=c_0$, so the back of the rarefaction wave travels at $-c_0$. At the dry edge $c=0$ the Saint-Venant dry-front condition givesThus the fan occupies $-c_{0}t\le x\le 4c_{0}t$, while $x<-c_{0}t$ remains undisturbed. Its full length grows at $5c_0$; there is no finite-depth uniform shelf between the fan and the dry front. The optional interior depth is $h=h_0[(4c_0-x/t)/(5c_0){]}^2$.

The fan's minus-family characteristic curves are rays $x=\xi t$. For the plus family in the fan, $dx/dt=u+c=(8c_0+3x/t)/5$; integration gives $x=4c_{0}t+C{t}^{3/5}$, with constants selected at entry from the reservoir. These curved paths transport the constant plus Riemann invariant.

Constant density in each layer and negligible ambient acceleration far from the front giveThis primitive momentum conservation equation remains valid with slowly varying $B$. In conservative form it iswhere the right side is the longitudinal component of the hydrostatic pressure force on the sloping sidewalls. Omitting that term while retaining variable $B$ would give the wrong acceleration.

Together with volume conservation, the principal coefficient matrix in $(h,u)$ isIts real eigenvalues and characteristic curves areThese are downstream and upstream gravity waves relative to the moving fluid. The factor $1/2$ is geometric: the hydraulic depth is $A/b(h)=h/2$. Both waves move downstream in a sufficiently fast current, whereas one can carry information upstream in a slower current.

For $B=1$, differentiate $c^2=g^{\prime }h/2$ and combine the two equations. The Riemann invariants areFor variable $B$ they instead obeywhich also checks that the invariant claim depends on a prismatic channel.

Write $g^{\prime }=g({\rho }_f-{\rho }_a)/{\rho }_f>0$. The Saint-Venant dry-front condition describes a liquid advancing over a dry bed, or into an ambient whose inertia is negligible, as when ${\rho }_a/{\rho }_f\ll 1$. The edge has $h_f=0$; its velocity is obtained by continuing the appropriate shallow water Riemann invariant to zero depth. There is no independent universal finite-depth head Froude number. In the triangular geometry a release of a resting layer gives $u_f=2\sqrt{2g^{\prime }{h}_0}$, whereas a rectangular release gives $2\sqrt{g^{\prime }{h}_0}$.

The Benjamin deep-ambient front condition includes the inertia of an ambient that must turn around a finite-depth gravity current head. In the ideal deep-ambient Boussinesq approximation, ${\rho }_a/{\rho }_f$ is close to one andHere $h_f$ is the current depth just behind its head, $u_f$ its speed relative to the far ambient, and $F$ the ideal head Froude number. The ambient is much deeper than $h_f$; a finite upper-layer depth changes the closure. The lower density must exceed the ambient density for a bottom current. A non-Boussinesq head requires an appropriate density-ratio correction rather than using $F=\sqrt{2}$ at every ratio. The Saint-Venant and Benjamin models therefore describe different front physics, even though both use the same long-wave equations in the interior.

Let $\ell$ and $\mathcal{T}$ denote longitudinal length and time scales, and let $U$ be a typical streamwise velocity. The shallow water approximation requires $h/\ell \ll 1$, slow changes in the channel geometry $h|B^{\prime }/B|\ll 1$, and small longitudinal slopes. The vertical velocity is then of order $Uh/\ell$; its acceleration must be small compared with the restoring acceleration, for example $Uh/(\mathcal{T}\ell )+U^{2}h/{\ell }^2\ll g^{\prime }$ in a reduced-gravity layer. The leading pressure is hydrostatic pressure, and cross-sectional horizontal velocity is approximated by $u(x,t)$. High Reynolds number permits negligible interior viscous terms but does not itself guarantee the geometric or hydrostatic assumptions. We also require $B>0$ and $h>0$ in the wetted interior. The dimensionless $B$ need not be small: it measures the cross-section's opening angle, while the smallness condition concerns streamwise variation.

The cross-sectional area and volume transport areFor an impermeable bed and no fluid entrainment, volume conservation givesor, equivalently,The geometrical source term expresses the depth change needed when a moving layer encounters a wider or narrower triangular section.