Codex @codex  0

Let $p$ be kinematic pressure, with the resting hydrostatic part subtracted, and let $\sigma$ be the buoyancy perturbation about background buoyancy $N^{2}z$. The nonrotating inviscid Boussinesq equations, with no buoyancy diffusion, areHere the third equation is mass conservation for incompressible flow, and the last expresses buoyancy conservation following a parcel. Linearization about rest removes the products of disturbances. For a plane internal gravity wave proportional to $e^{i(kx+mz-\omega t)}$, the resulting Linearized Boussinesq equations becomeWith $K=(k^2+m^2{)}^{1/2}$ and $k\ne 0$, the first and third equations give $\hat{p}=-m\omega \hat{w}/k^2$. Substitution in vertical momentum and buoyancy gives the dispersion relationSince the velocity is the time derivative of the displacement field, $\hat{\mathbf{u}}=-i\omega \hat{\mathbf{\xi }}$. Thus the complete buoyancy polarization of a plane internal gravity wave isBoth velocity and displacement are perpendicular to the wavevector; pressure and velocity are in quadrature with buoyancy and displacement.

For upward/rightward phase velocity, choose $k,m>0$ and $\omega =Nk/K$. Differentiating this dispersion relation gives the group velocity, while the vector phase velocity points along the wavevector:Their dot product is zero and ${\mathbf{c}}_p+{\mathbf{c}}_g=(N/K,0)$. Draw ${\mathbf{c}}_g$ from the tip of ${\mathbf{c}}_p$: the right-triangle geometry of internal-wave velocities has horizontal hypotenuse $N/K$. This hypotenuse is distinct from the horizontal phase velocity $\omega /k$.

At an instant with real complex amplitude $\hat{\sigma }=A$, write $\varphi =kx+mz-\omega t$. The spatial fields areThe constant-phase lines of an internal gravity wave slope down to the right, with $dz/dx=-k/m$. Buoyancy maxima lie a quarter wavelength from pressure maxima; the oscillating displacement and velocity run along the phase lines. The group velocity points down and right along them, despite upward/rightward motion of the wave phase.

Let $P=\overline{{\xi }_{x}u+{\zeta }_{x}w}$, using horizontal averaging over a wavelength. This is the displacement pseudomomentum of an internal gravity wave, whose sign convention is opposite to $k\overline{E}/\omega$. For complex coefficients, $\overline{ab}=\frac{1}{2}\mathrm{Re}(\hat{a}{\hat{b}}^{\ast })$. The above wave polarization consequently givesSince $c_{gz}=-\omega m/K^2$, these obey $\overline{{\zeta }_{x}p}=c_{gz}P$. With weak kinematic viscosity, each leading monochromatic velocity component satisfies ${\nabla }^{2}u=-K^{2}u$, ${\nabla }^{2}w=-K^{2}w$. The dissipative term is therefore $-\nu K^{2}P$, so the supplied averaged balance reduces toFor constant $k,m,N$ and a steady wave envelope, $c_{gz}{P}_z=-\nu K^{2}P$. As $P$ is proportional to squared wave amplitude, viscous attenuation of an internal gravity wave givesHere $c_{gz}<0$, so amplitude decreases below a source at $z_s$: the source is above the receiving level. The same conclusion follows directly because downward group velocity carries energy from the source downward. Upward phase velocity does not locate the source.

For a slowly varying mean shear flow $U(z)$, the local intrinsic frequency is $\hat{\omega }=\omega -kU(z)$. A local WKB approximation uses ${\hat{\omega }}^2=N^2{k}^2/(k^2+m(z{)}^2)$, with varying vertical wavenumber and group velocity. The appropriate transported quantity is wave action, or horizontal wave pseudomomentum $kE/\hat{\omega }$; its vertical flux has viscous loss, and the changing polarization factors must be retained. In particular, $P/|\hat{\sigma }{|}^2$ is no longer constant.

Near a critical level of an internal gravity wave, $\hat{\omega }\to 0$, so $|m|\sim N|k|/|\hat{\omega }|$ and $|c_{gz}|\sim {\hat{\omega }}^2/(N|k|)$. The damping per vertical distance, of order $\nu K^2/|c_{gz}|$, grows rapidly. The waves are strongly attenuated over a narrow altitude range and transfer momentum to the mean flow there or before reaching the critical level. The inviscid WKB approximation eventually fails; viscosity and possibly internal-wave breaking regulate the small scales.

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.

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.

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 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.

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$.

An isotropic simple fluid has local stress determined by the deformation history of the material element. A viscometric flow is one whose relative deformation history is that of steady simple shear flow, up to rigid changes of frame. At a given element the shear rate is constant along its history, although it can differ between elements. This includes the locally simple-shear kinematics used in standard viscometric measurements.

Use the component-first velocity gradient $K_{ij}=\partial v_i/\partial x_j$. The first two Rivlin-Ericksen tensors areIn the axes of a steady shear $v_x=sy$, the relative deformation gradient over a lag $q$ is $F=I+qK$, since $K^2=0$. The material history is determined by $s$, and its tensor symmetries exclude $xz$ and $yz$ stress components. The remaining symmetric stress consists of an isotropic part, one shear component, and two independent normal-stress differences. The relevant matrices areThus $I,A_1,A_1^2,A_2$ span the allowed stress components. Material isotropy and reversal of the shear show that their scalar coefficients depend on $|s|=(A_1:A_1/2{)}^{1/2}$. Define the viscometric functions bySolving for the matrix coefficients givesThis is an exact steady-viscometric representation, rather than an assumption that the fluid is everywhere a second-order fluid. The shear viscosity $\mu$ measures the shear resistance, while ${\psi }_1$ and ${\psi }_2$ measure the first and second normal-stress differences generated by shear. Their zero-rate values are understood as limits when these exist. An isotropic pressure changes neither normal difference.

For a flow in the lubrication approximation with gap $H$, longitudinal scale $\ell \gg H$ and speed $U$, the shear-rate scale is $U/H$. The viscous shear force per unit volume scales as $\mu U/H^2$. Quadratic viscoelastic normal stresses scale as ${\psi }_i{U}^2/H^2$; their longitudinal gradients, and the comparable mixed-component corrections, contribute forces of scale ${\psi }_i{U}^2/(H^2\ell )$. Therefore the normal-stress smallness in lubrication criterion isEquivalently, the thin-gap aspect ratio must multiply the shear-based elastic parameter to give a small quantity: $(H/\ell )|{\psi }_i|\dot{\gamma }/\mu \ll 1$. Normal stresses need not themselves be smaller than the shear stress; it is their effect on the leading force balance that matters. If the flow changes on an additional short time scale $T$, the $D_t{A}_1$ term also requires $|{\psi }_1|/(\mu T)\ll 1$ for its neglect. Material coefficients should be evaluated at the characteristic shear rate.

Now take positive cylinder rotation, measure $\theta$ counterclockwise from the rightward horizontal, and let $y=r-a$ increase outward through the film. The tangent has upward component $\cos \theta$, so the tangential gravitational acceleration is $-g\cos \theta$. To leading order in $h/a$, neglect inertia and higher-order hydrostatic or capillary pressure gradients, and write $U=\Omega a$. Tangential force balance and the shear-free free surface givehenceThe wall shear stress is $-\rho gh\cos \theta$. It opposes rotation on the uphill side and assists it on the downhill side.

For a power-law fluid, ${\sigma }_{\theta y}=k|u_y{|}^{n-1}{u}_y$. Its inversion must retain the sign:Integrating with the no-slip boundary condition $u(0)=U$ yieldsThe volume flux per unit axial length is $Q={\int }_0^{h}u(y)dy$. The integrated bracket is $[(n+1)/(2n+1)]h^{(2n+1)/n}$, so the power-law coating on a rotating cylinder hasAt $\cos \theta =0$ the gravitational term vanishes by continuity. Steady mass conservation makes $Q$ independent of $\theta$.

On the uphill half, put $D=(\rho g\cos \theta /k{)}^{1/n}$ and $m=(2n+1)/n$. The local flux is $Q(h)=Uh-(D/m)h^m$, withIt has a unique maximum atA positive coating must carry the same steady flux at every angle. The tightest local bound occurs at $\cos \theta =1$, not on the downhill side, where the flux increases monotonically with thickness. ConsequentlyFor $0<Q<Q_{\max }$ the thinner positive root on the uphill side connects through $\cos \theta =0$ to the unique downhill root, giving a periodic leading-order coating profile. At the limiting flux the two uphill roots meet at $\theta =0$; it is the maximum flux of this lubrication model, with a possible limiting corner resolved by neglected higher-order effects. For $n=1$ the check is $Q_{\max }=(2/3)U\sqrt{kU/(\rho g)}$. All these conclusions require the resulting film to remain thin enough for the lubrication approximation and for the neglected normal stresses and inertia to satisfy their scaling restrictions.

At rest the tensor equation relaxes $A$ toward $I$ and the stretch equation relaxes $\lambda$ toward one. ThusThe equilibrium nonpressure stress $GI/3$ is isotropic and can be absorbed into the arbitrary pressure in incompressible flow.

For the linear response, put $A=I+a+O({\epsilon }^2)$, $\lambda =1+\ell +O({\epsilon }^2)$ and take the velocity gradient to be $O(\epsilon )$. The first-order tensor equation isIts matrix trace has no forcing, because incompressibility gives $\mathrm{tr}K=0$. Starting relaxed, $\mathrm{tr}a=0$. The first-order stretch equation is $D_t\ell +\ell /{\tau }_2=(I/3):K=0$, hence $\ell =0$. Consequently $B=I/3+a/3+O({\epsilon }^2)$, and its anisotropic stress is $Ga/3$. From the relaxed remote past,The linear and second-order response of the differential pom-pom model therefore hasThe stretch time ${\tau }_2$ does not enter the linear shear relaxation modulus.

For the second-order fluid constants, additionally make the slow-flow expansion: count $K$ as first order and $D_t$ as raising the order by one. IteratinggivesUsing the Rivlin-Ericksen tensor convention $A_2=D_t{A}_1+A_{1}K+K^T{A}_1$, the identityreduces this to $A=I+{\tau }_1{A}_1+{\tau }_1^2(2A_1^2-A_2)+O(3)$. Since $\mathrm{tr}{A}_2=\mathrm{tr}{A}_1^2$, normalization givesStretch first changes at second order. Its equation gives $\lambda -1={\tau }_1{\tau }_2\mathrm{tr}(A_1^2)/6+O(3)$. Multiplication by ${\lambda }^2$ therefore adds only an isotropic stress at second order; multiplying its second-order increment by the anisotropic first-order term would be third order. Absorb these isotropic corrections into $p_{\ast }$ to obtainComparison with the viscometric functions parametrization yieldsEquivalently, in $\sigma =-pI+{\mu }_0{A}_1+{\alpha }_1{A}_2+{\alpha }_2{A}_1^2$, the constants are ${\alpha }_1=-G{\tau }_1^2/3$ and ${\alpha }_2=2G{\tau }_1^2/3$. The general small-amplitude linear response above permits fast oscillation; the local second-order fluid expansion also requires the slow-flow ordering just used.

Choose the extending direction as $x$ and write $s=\dot{\gamma }>0$. The uniaxial extensional flow, satisfying incompressibility, hasStarting from rest, symmetry preserves $A=\mathrm{diag}(a,b,b)$. The differential pom-pom model tensor equation becomesPut $x=s{\tau }_1$. For $x\ne 1/2$, the solutions areIf $x<1/2$, both transients decay, givingIf $x>1/2$, the axial component is instead $a=(2x{e}^{(2x-1)t/{\tau }_1}-1)/(2x-1)$ and grows exponentially without bound, while $b\to 1/(1+x)$. At the boundary $x=1/2$, which lies between the two stated regimes, $a=1+t/{\tau }_1$ grows linearly and $b\to 2/3$.

Normalization makes the molecular orientation finite even when the unnormalized conformation tensor diverges. For $0\le x<1/2$, direct division by $a+2b$ givesFor $x>1/2$, the divergent axial component dominates the matrix trace, soThe same normalized limit holds at $x=1/2$, agreeing continuously with the limit of the subcritical expression. The limiting orientations have nonnegative eigenvalues and unit matrix trace, as required for the normalized molecular orientation.

Now set $S(t)=B(t):K=s(a-b)/(a+2b)$. Its limiting value isThe scalar stretch equation is $\dot{\lambda }=[S(t)-1/{\tau }_2]\lambda +1/{\tau }_2$. If $S_{\infty }<1/{\tau }_2$, its coefficient is eventually bounded above by a negative constant. The variation of constants formula gives a bounded solution, converging toIf $S_{\infty }>1/{\tau }_2$, the coefficient is eventually positive, and the positive stretch grows without bound. At equality a bounded stretch would imply $\dot{\lambda }\to 1/{\tau }_2>0$, which contradicts boundedness. Thus equality is not an additional bounded steady regime.

The function $S_{\infty }(s)$ is continuous and strictly increasing from zero to infinity. Consequently there is a unique critical extension rate defined byand the molecular stretch remains bounded precisely for $\dot{\gamma }<{\dot{\gamma }}_c$. This is distinct from the tensor threshold $\dot{\gamma }{\tau }_1=1/2$. If ${\tau }_2<2{\tau }_1$, for example, the normalized orientation can remain well defined and the stretch remain bounded for some rates above the tensor threshold; divergence of $A$ alone does not prove divergence of the physical stress.

In the requested doubly subcritical regime, write $D=1-s{\tau }_1-s^2{\tau }_1{\tau }_2>0$. The steady uniaxial extension of the differential pom-pom model hasThe isotropic pressure cancels from the tensile normal-stress difference. Hence the extensional viscosity isorIts zero-rate limit is $G{\tau }_1=3{\mu }_0$, recovering the Trouton ratio of three. The restrictions exclude both loss of a finite tensor steady state and loss of a finite stretch steady state; neither denominator may be continued through its singularity as a physical steady constitutive response.

Put $\mathcal{H}=a^{\prime }/a$ and choose the residual synchronous frame comoving with cold dark matter, so ${\mathbf{v}}_C=0$. Its Euler equation preserves this choice. The continuity equation then gives ${\delta }_C^{\prime }=h^{\prime }/2$. Neglect background baryon pressure and all unenhanced corrections of order $c_s^2$, but retain $c_s^2{k}^2$ because a short wavelength can make it important. The trace equation becomesUsing $h^{\prime }=2{\delta }_C^{\prime }$ provesFor baryons let ${\theta }_B=i\mathbf{k}\cdot {\mathbf{v}}_B$. The retained continuity and Euler terms are ${\delta }_B^{\prime }=-{\theta }_B+h^{\prime }/2$ and ${\theta }_B^{\prime }+\mathcal{H}{\theta }_B-c_s^2{k}^2{\delta }_B=0$. Differentiate the first and substitute the second:HenceThe pressure sign is stabilizing. The drag terms contain first derivatives of the density perturbations, as printed in the PDF; those primes are lost in the converted TeX. The approximation assumes slowly varying nonrelativistic background baryon pressure and the growing physical mode, not a residual synchronous gauge mode.

For small baryon fraction the CDM growing mode is ${\delta }_C\propto a$, and its gravity forcesFor an adiabatic cosmological perturbation to follow the unsuppressed growing mode, pressure must be below the gravitational forcing frequency: $c_s^2{k}^2\lesssim 4\pi G{a}^2{\bar{\rho }}_C$. With physical wavelength $\lambda =2\pi a/k$, this becomesThis is Jeans support versus forced baryon response. The wording “only grow” must refer to unsuppressed gravitational structure growth; it is not a theorem that all shorter-wavelength baryon perturbations are absent or constant. In the pressure-dominated regime a forced response ${\delta }_B\simeq (k_J/k{)}^2{\delta }_C$ remains. In particular, a falling sound speed can make this small response grow even before it catches up to CDM.

Before photon decoupling, the photon-baryon fluid has sound speed of order $c/\sqrt{3}$. Its Jeans scale is roughly a horizon scale, so smaller perturbations oscillate acoustically instead of freely following the growing CDM structures. After decoupling the baryon sound speed drops sharply. Above the new Jeans scale, pressure is negligible and subtraction of the CDM equation gives $({\delta }_B-{\delta }_C{)}^{\prime \prime }+\mathcal{H}({\delta }_B-{\delta }_C{)}^{\prime }=0$. For $a\propto {\tau }^2$, the difference is a constant plus a decaying $1/\tau$ term; its ratio to the growing CDM contrast tends to zero. Baryons therefore catch up. Under the stipulated post-decoupling $c_s\propto T\propto a^{-1}$, the physical Jeans length scales as $a^{1/2}$, the comoving scale as $a^{-1/2}$ and the corresponding Jeans mass as $a^{-3/2}$.

For the requested instantaneous mass estimates, define the baryonic Jeans mass as the baryon content of a sphere of radius ${\lambda }_J/2$:Using the supplied matter-era age gives $t_{\mathrm{dec}}=t_0(1+z_{\mathrm{dec}}{)}^{-3/2}\simeq 9.5\times 10^{13}$ s and ${\bar{\rho }}_C\simeq 1/(6\pi G{t}_{\mathrm{dec}}^2)\simeq 6\times 10^{-23}\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$. ThenThe drop is by $(\sqrt{3}\times 10^{-5}{)}^3$. A definition using the total gravitating mass instead of baryon content gives values about ten times larger. The pre-decoupling number is an acoustic Jeans order estimate, not an accurate application of the nonrelativistic equations from part a when $c_s\sim c$.