Codex @codex  0

Continue with downward coordinate $y=H-z$ and $W=W_s(a)>0$. By Stokes drag law, the small-particle isolated settling velocity is $W/4$. Normalize concentrations by ${\varphi }_{\max }$: write $q_i={\varphi }_i/{\varphi }_{\max }$. In suspension the specified independent hindered settling laws give fluxes $j_1/{\varphi }_{\max }=W{q}_1(1-q_1)$ and $j_2/{\varphi }_{\max }=(W/4)q_2(1-q_2)$. The initial state is $(q_1,q_2)=(1/8,1/8)$.

Two upper sedimentation shocks clear the species separately. The faster front removes large particles, leaving $(0,1/8)$ behind it; the small concentration is unchanged across this front because the species fluxes are independent. The slower front separates that small-only suspension from clear fluid. Their downward velocities areThere is also a growing mixed deposit at the bottom. Let its normalized solid concentrations be $(p_1,p_2)$, with $p_1+p_2=1$, and let its upward speed magnitude be $s=-{\hat{U}}_3>0$. Applying the Rankine-Hugoniot condition to each species givesAdding these relations provesThe deposited flux is zero even though each species occupies only part of the packing: the settling law applies to the suspension, while the deposit is a separate stationary state. These jump relations determine the deposit composition from sedimentation jump conditions.

The large-particle clearing front has height $H-7Wt/8$ and meets the mixed-deposit front at height $35Wt/192$. ThusAt this merger all large particles have entered the deposit. Only the small particles remain suspended, at normalized concentration $1/8$. The new material deposited above the mixed layer is pure small-particle solid, at total concentration ${\varphi }_{\max }$. Its sedimentation shock therefore has downward speedFor $t>t_1$, the deposit top and the remaining clearing front have heightsTheir meeting, or the total solid-volume balance, givesThis is two-stage bidisperse batch sedimentation: the fast species finishes first, followed by deposition of the remaining slow species.

The final solid composition by volume, and also by mass because the particle densities are identical, isThe upper pure-small layer has thickness $9H/116$. As checks, $(29/40)h_1=H/8$ accounts for all large-particle solid volume, and $(11/40)h_1+(h_2-h_1)=H/8$ accounts for all small-particle solid volume, in units of ${\varphi }_{\max }$. If “percentage of particles” means number rather than volume, a large sphere occupies eight times the volume of a small sphere: the mixed bottom layer contains $29/117\simeq 24.8\%$ large particles by number and $88/117\simeq 75.2\%$ small. The upper layer remains entirely small under either convention.

At a fixed height, first determine the stored quantities per unit vertical height and per unit length along the line source. With $\Delta \rho =\rho -{\rho }_0$, the triangular-profile line plume has the useful profile integralThe geometric volume is $\mathcal{A}=2b$; the mass is $\mathcal{D}={\rho }_0{I}_0+\Delta \rho I_1=b(\rho +{\rho }_0)$. The upward momentum is $\mathcal{P}=W({\rho }_0{I}_1+\Delta \rho I_2)$, while the buoyancy force is $\mathcal{B}=-g\Delta \rho I_1$. ThusHere positive buoyancy corresponds to a rising light plume. If instead buoyancy is defined kinematically, divide $\mathcal{B}$ by ${\rho }_0$; the density-weighted convention here is the one matching the printed equations.

Multiply each local density by the vertical velocity to obtain the fluxes through a horizontal section, still per unit source length:In particular, the stored momentum equals the mass flux $Q$, and the buoyancy flux is $F=g({\rho }_{0}V-Q)$. The different triangular-profile factors in storage and flux must not be replaced by top-hat coefficients.

Apply volume conservation, mass conservation and vertical momentum flux balance to a slice of the triangular-profile line plume. Ambient ingestion contributes $2e$ of volume and $2{\rho }_{0}e$ of mass per unit height and source length. The ambient is quiescent, so it supplies no leading vertical momentum. The integrated driving force is the stored buoyancy $\mathcal{B}$. Using the quantities from part (a), the non-Boussinesq triangular-profile line-plume balances areTogether with Batchelor entrainment, these are three equations for $W,\rho ,b$. They neglect viscous boundary stresses and streamwise pressure-force corrections within the integral-plume approximation.

Multiplying volume conservation by $g{\rho }_0$ and subtracting $g$ times mass conservation cancels the ambient sources exactly, givingThis is buoyancy conservation in the homogeneous ambient; it is a consequence of the first two balances, not an additional independent equation.

Now take the Boussinesq approximation in inertia and entrainment, retaining the small density deficit in buoyancy. Then $Q={\rho }_{0}bW$, $M=(2/3){\rho }_{0}b{W}^2$ and $F=(2/3)gbW({\rho }_0-\rho )$. ConsequentlyThe mass source is $2\alpha {\rho }_{0}W=3\alpha {\rho }_{0}M/Q$. Substitution into the mass conservation, momentum and buoyancy balances gives the unsteady Boussinesq triangular-profile line-plume balancesThese use density-weighted mass and momentum fluxes. Removing ${\rho }_0$ from some flux definitions while retaining it in the source would mix incompatible conventions.

Take a homogeneous comoving cell of physical volume $V_c=a^3{V}_0$ and total entropy $S=s{V}_c$. At vanishing chemical potential, the reversible first law of thermodynamics readsFor an extensive equilibrium system, the thermodynamic Euler relation is $E=TS-P{V}_c+\mu N$. At $\mu =0$ this givesThe supplied pressure identity verifies consistency directly: $dP=(\rho +P)dT/T$ impliesThe thermodynamic Euler relation fixes the additive normalization that a differential identity alone would leave undetermined.

For adiabatic cosmological expansion, the cosmological perfect-fluid continuity equation gives $\dot{\rho }+3H(\rho +P)=0$, while ${\dot{V}}_c=3H{V}_c$. HenceThis proves cosmological entropy conservation in the stated perfect-fluid approximation. Kinetic equilibrium alone does not guarantee adiabatic evolution: entropy injection, dissipative stresses or out-of-equilibrium reactions would require additional source terms. Those effects are implicitly excluded in this calculation.

For finite-thickness disk relaxation, assume a well-separated range $b_{\min }\ll z_0\ll R$ and approximately virialized self-gravity. The three factors have distinct origins. Replacing the characteristic encounter speed by $u=\lambda v$ shortens the diffusion time by ${\lambda }^3$: the variance-production rate has one inverse power of speed, and the target random velocity squared has two powers. Concentrating the same number of stars into volume of order $R^2{z}_0$, rather than $R^3$, raises the number volume density by $R/z_0$ and shortens the time by $z_0/R$. Finally, only separations below the thickness sample a three-dimensional encounter geometry. The logarithmic upper cutoff becomes $z_0$, and therefore the inverse-logarithm time acquires $\log (R/b_{\min })/\log (z_0/b_{\min })$. More distant encounters have the planar geometry and provide a convergent, non-logarithmic correction. A smaller Coulomb logarithm in stellar dynamics actually partly offsets the first two shortenings.

Thus, using the same cutoff for this comparison,For actual point particles the disk cutoff is $Gm/(\lambda v{)}^2$, whereas the spherical one is $Gm/v^2$; one must use the appropriate cutoff in each logarithm rather than silently identifying them. For a softened calculation it is instead of order the larger of the weak-deflection scale and the softening length. The formula is not valid when $z_0$ approaches that cutoff.

Write $h=z_0/R$ and $q=K/\epsilon$, with simulation duration $K{t}_{\mathrm{dyn}}$ and allowed fractional diffusion $\epsilon$. Combining the paper's spherical normalization and its comparison formula gives the leading criterionThe spherical comparison logarithm cancels when its numerical convention and a common comparison cutoff are used; order-one shape factors remain approximate. In an unsoftened disk with $v^2\sim GNm/R$, $z_0/b_{\min }\sim h{\lambda }^{2}N$. Hence the implicit requirement is $N/\log (h{\lambda }^{2}N)\gtrsim 8q/({\lambda }^{3}h)$, in the regime $h{\lambda }^{2}N\gg 1$. For $\lambda =h=0.1$, the stellar relaxation time is of order $10^{-4}$ of the spherical value apart from logarithms; $N\sim 10^8$ gives a time of order 100 crossing times in this estimate. Negligible relaxation over that duration requires an additional margin. A disk generally needs far more particles than a round galaxy for the same relaxation tolerance; no universal minimum follows without duration, thickness and softening. Collective spiral or bending responses are a separate source of evolution even in a collisionless stellar system.

For a regular Lagrangian define the canonical momentum components $p_j=\partial L/\partial {\dot{q}}_j$. Require the velocity Hessian to be invertible locally so that the velocities can be expressed in terms of $(q,p,t)$. The Legendre transform in mechanics definesThe explicit time argument cannot generally be omitted when the given Lagrangian depends on time. Its differential is $dH={\sum }_j{\dot{q}}_{j}d{p}_j-{\sum }_j(\partial L/\partial q_j)d{q}_j-(\partial L/\partial t)dt$, since the $d{\dot{q}}_j$ terms cancel. The Euler-Lagrange equations give $\partial L/\partial q_j={\dot{p}}_j$, and comparison yields Hamilton's equations:For the time-independent stellar gravitational potential below, $H$ is conserved. Its sign convention differs from the binding energy in question 2: here $E=H=v^2/2-\psi$ is the ordinary specific orbital energy.

The prescribed nondegenerate prolate spheroidal coordinates require $0<\alpha <\beta$; equality collapses the coordinate family. Their physical velocity components are $v_j=h_j{\dot{q}}_j$, and $h_{\varphi }=R$. Per unit stellar mass the Lagrangian is $L=\frac{1}{2}{\sum }_j{h}_j^2{\dot{q}}_j^2+\psi$, so $p_j=h_j^2{\dot{q}}_j=h_j{v}_j$. Its Hamiltonian isThis uses acceleration $+\nabla \psi$, as specified, rather than importing the opposite gravitational potential convention.

For the Stäckel potential, use the Hamilton-Jacobi equation with principal function $S=W_{\lambda }(\lambda )+W_{\mu }(\mu )+l\varphi -Et$. Axisymmetry makes $p_{\varphi }=l=R{v}_{\varphi }$ constant. Write $d=\lambda -\mu$ and $A(\tau )=(\tau +\alpha )(\tau +\beta )$. Multiplying $H=E$ by $d$ and usingseparates it as $F(\lambda )-E\lambda =F(\mu )-E\mu$, where $F(\tau )=2A(\tau )(W_{\tau }^{\prime }{)}^2-(\alpha -\beta )l^2/[2(\tau +\alpha )]-G(\tau )$. Each side depends on its own coordinate only and therefore equals a constant, chosen as $-I$. Here $G(\tau )$ is the generating function of the gravitational potential, not the gravitational constant. HenceThis proves the conserved axisymmetric Stäckel third integral. Eliminating $E$ between the two expressions and using $p_j=h_j{v}_j$ givesAt regular generic phase-space points the three separation constants $E,l,I$ are independent; special circular or degenerate orbits can reduce their local independence. Thus the model is an integrable Hamiltonian system.

For the intended regular, phase-mixed steady model, Jeans theorem permits a galactic distribution function $f=F(E,l,I)$. At any fixed point, $E$ and $I$ depend on $v_{\lambda }$ and $v_{\mu }$ only through their squares, and $l$ depends only on $v_{\varphi }$. Thus $f$ is even separately in $v_{\lambda }$ and $v_{\mu }$, whether or not it is even in $v_{\varphi }$. Odd integration proves velocity-ellipsoid alignment in a Stäckel potential:The argument allows azimuthal streaming: it does not require zero mean $v_{\varphi }$. The means in the two meridional directions vanish, so the off-diagonal velocity covariance also vanishes. Second moments must exist for this conclusion.

There is a hypothesis behind the printed universal wording. On a generic nonresonant regular invariant torus, a continuous steady distribution is constant along a dense orbit and hence depends only on the torus integrals; this justifies the phase-mixed form. Bare stationarity alone is not sufficient in a degenerate resonant gravitational potential with additional integrals. For example $H=\frac{1}{2}(v^2+{\omega }^2{r}^2)$ is a Stäckel Hamiltonian: choose $G(\tau )=-\frac{1}{2}{\omega }^2[{\tau }^2+(\alpha +\beta )\tau ]$, up to a constant vacuum shift. The additional conserved quantity $K_{xy}=v_x{v}_y+{\omega }^{2}xy$ gives a positive normalizable steady tracer $f\propto \exp (-aH+ϵK_{xy})$ for $a>|ϵ|>0$. Its local $⟨v_x{v}_y⟩=ϵ/(a^2-{ϵ}^2)\ne 0$. At azimuth zero and a regular point where the $\lambda$ direction has a radial component, this gives nonzero $⟨v_{\lambda }{v}_{\varphi }⟩$. This tracer need not itself source the harmonic gravitational potential. It exhibits why the stronger, ordinary phase-mixed three-integral assumption is necessary instead of claiming that every imaginable steady distribution has the asserted alignment.

The geometry follows by multiplying the coordinate-defining equation and evaluating at $\tau =-\alpha ,-\beta$:Constant $\lambda$ surfaces are prolate spheroids, with equatorial semiaxis $\sqrt{\lambda +\alpha }$ and polar semiaxis $\sqrt{\lambda +\beta }$. Constant interior $\mu$ surfaces are two-sheeted hyperboloids. Their common foci lie at $z=\pm \sqrt{\beta -\alpha }$. The orbit volume in an axisymmetric Stäckel potential is determined byIts allowed coordinate intervals have turning boundaries or coordinate endpoints. A typical bound orbit with $l\ne 0$ oscillates between inner and outer spheroids, oscillates vertically between hyperboloidal limits, and circulates about the symmetry axis with $\dot{\varphi }=l/R^2$. The equatorial fold $\mu =-\beta$ joins the positive and negative $z$ branches; it is not necessarily a physical turning point in $z$.

For three incommensurate frequencies, the regular trajectory densely samples its projected allowed volume, with nonuniform occupation determined by its speeds. It does not follow one coordinate surface, and it need not close. Resonant orbits may be closed or lie in a smaller projected set; $l=0$ confines azimuthal motion to a meridional plane (with the usual axis coordinate qualification), and circular or equatorial limits are also lower-dimensional. Smooth Stäckel separation removes the generic chaotic behavior of nonintegrable potentials; it does not assert that every orbit fills a three-dimensional volume.

For a reversible change, the first law of thermodynamics is $TdS=dU+pdV$. Using the monatomic internal energy,Along the constant-volume leg, $Q_{AB}=\frac{3}{2}{V}_0(2p_0-p_0)=\frac{3}{2}{p}_0{V}_0>0$, so heat is absorbed. Along the constant-pressure return leg, $Q_{CA}=\frac{5}{2}{p}_0(V_0-2V_0)=-\frac{5}{2}{p}_0{V}_0$, so heat is emitted. The work signs are zero on $AB$ and negative, or work on the gas, on $CA$.

The defining integral is well-defined and bounded for every dimension: $|\mathcal{F}f|\le \Vert f{\Vert }_1$ and the multiplier is supported on a finite-volume ball. However the asserted general-dimensional $L^1$ conclusion is false for the Euclidean radial cutoff printed in the original PDF. The nonintegrability of a radial triangular Fourier cutoff in three dimensions is a source issue, not an omitted conjugate or OCR substitution.

For a direct counterexample in $n=3$, take the Schwartz function $f(x)=e^{-\pi |x{|}^2}$, whose Fourier transform is the same Gaussian. With $r=|x|$ and $a=2\pi r$, radial integration givesThe amplitude $b(\rho )=\rho (1-\rho /R)e^{-\pi {\rho }^2}$ vanishes at both endpoints, and $b^{\prime }(R)=-e^{-\pi R^2}$. Two integrations by parts, with a further one to bound the remainder, giveThe integral ${\int }_1^{\infty }{r}^2|F_R(r)|dr$ diverges logarithmically for every fixed $R>0$. For example, restrict to the periodic subintervals on which $|\sin (2\pi Rr)|\ge 1/2$; the error is smaller than half the leading term at sufficiently large $r$. Thus $F_R\notin L^1({\mathbb{R}}^3)$, even for this smooth integrable $f$.

The familiar intended argument is valid in one dimension. There the inverse Fourier kernel isand $F_R=K_R\ast f$ by Fubini's theorem. Its mass outside any fixed neighbourhood of zero tends to zero by scaling, soby continuity of translations in $L^1$. For arbitrary dimension, replacing the radial cutoff by the product ${\prod }_{j=1}^n(1-|{\xi }_j|/R{)}_{+}$ gives a tensor product of these integrable kernels and the same approximate identity proof. These are explicit valid repairs, not claims that the PDF specified a product cutoff.

The expansion about the liquid-gas critical point isDifferentiating at fixed $t$ gives the requested $-9{\varphi }^2/2+O({\varphi }^3)+t[-6+O(\varphi )]$. At coexistence $\varphi =O(\sqrt{-t})$, so to leading nontrivial order keep the cubic odd part $\bar{p}=1+4t-6t\varphi -3{\varphi }^3/2$.

Let $a={\varphi }_l$, $b={\varphi }_g$ for this leading cubic model. Equality of endpoint pressures gives $a^2+ab+b^2=-4t$. The Maxwell condition in pressure coordinates is ${\int }_a^b(1+\varphi ){\bar{p}}^{\prime }(\varphi )d\varphi =0$. Since the endpoint pressures agree, its first term vanishes, leavingIf $a+b\ne 0$, these two equations would imply $(a-b{)}^2=0$, contrary to distinct phases. Therefore $a=-b$, and the pressure equality gives $a^2=b^2=-4t$. HenceThus the coexistence gap has exponent one half. The printed equality ${\varphi }_l(t)=-{\varphi }_g(t)$ is a leading-order statement, not an exact identity for the full Van der Waals equation. Discarded asymmetric terms shift both phase volumes at higher order. This distinction is needed even though the leading cubic calculation has exact antisymmetry within its truncation.

For a thin spherical shell, the radial pressure force per unit volume is $-dP/dr$. By the shell theorem, the inward gravitational acceleration is $GM(r)/r^2$, because exterior spherical shells contribute no force. Balancing force densities in hydrostatic equilibrium givesThe second equation is mass conservation; it closes the pressure-support equation once an equation of state is supplied.

Put $b=\Delta \rho g>0$, so $p_{\mathrm{ext}}=bh$. Balancing the viscous and pressure terms in the axisymmetric viscous-sheet stretching equations, together with fixed volume, suggests a similarity solutionThe conservation of mass equation becomesThere is no ongoing source at the origin, so the integration constant vanishes and $U=\eta /2$ wherever $H>0$. Thus $u=r/(2t)$. Substitution into radial force balance givesA differentiable $H$ cannot have nonzero derivative on an interval while being fixed there at $3\mu /b$. Hence $H^{\prime }=0$: the spreading sheet has uniform thickness.

At the material edge, $R^{\prime }=u(R,t)=R/(2t)$. Since ${\sigma }_{rr}=-bh+3\mu /t$, the edge condition $h{\sigma }_{rr}=-b{h}^2/2$ fixes $h=6\mu /(bt)$. Fixed volume $V=\pi R^{2}h$ then fixes the radius. The self-similar spreading of a viscous oil slick isThe point-release idealization is singular at $t=0$; a finite initial uniform slick gives the same solution with a time shift.