Codex @codex  0

Stable stratification is a constraint on motion, not an isolated parcel spring. A light-over-heavy density profile has $N^2=-(g/{\rho }_0){\bar{\rho }}_z>0$. An adiabatically displaced parcel has buoyancy anomaly $b=-N^2{\xi }_z$ relative to its surroundings. If one neglects the pressure response and movement of surrounding fluid, the vertical equation would be ${\ddot{\xi }}_z=-N^2{\xi }_z$. A finite spherical blob cannot satisfy that isolated-oscillator picture: incompressible flow requires surrounding fluid motion, and the pressure perturbation communicates its acceleration to the rest of the fluid.

This can be shown directly with the linear Boussinesq equations. For a uniform background, take a Fourier wavevector $({\mathbf{k}}_h,m)$. Pressure removes the component of the buoyancy force parallel to this wavevector. Taking the divergence of momentum to determine pressure and projecting the vertical component givessoThis pressure projection in a displaced stratified blob produces frequencies less than $N$ except for wavevectors with $m=0$. A localized spherical perturbation contains a range of directions and therefore a range of frequencies. It deforms and emits internal gravity waves, rather than remaining one sphere oscillating coherently at $N$. The added mass of a rigid sphere is a useful inertia analogy, but its numerical correction does not define an exact oscillation frequency for a deformable material blob in stratified fluid.

Internal waves and energy transport. The restoring force supports anisotropic internal gravity waves: their frequency depends on the direction, not just the length, of the wavevector. For a two-dimensional positive-frequency branch, $\omega =Nk/\sqrt{k^2+m^2}$ with $k>0$. Differentiation givesThe vertical group and phase velocities have opposite signs. Wave energy therefore travels along a direction different from the direction of phase propagation, a point essential in selecting radiation conditions for wave-generating obstacles.

For constant $N$, the linear perturbation energy density per reference mass isDot momentum with velocity and multiply buoyancy evolution by $b/N^2$; the $bw$ terms cancel, leaving ${\partial }_t\mathcal{E}+\nabla \cdot (\Pi \mathbf{u})=0$. The second term is the available potential energy stored by displacement against the stable density profile. In an ideal fluid, that energy is exchanged with kinetic energy or transported away, not irreversibly dissipated by buoyancy diffusion.

Steady flows, layers and topography. Stable stratification penalizes vertical displacement and favors flows with small vertical motion when the stratified Froude number $U/(NH)$ is small. The natural vertical displacement scale associated with kinetic energy $U^2$ is of order $U/N$; obstacles requiring much larger displacements can cause blocking and diversion instead of simple overpassing. In the ideal model fluid retains its buoyancy label, so motion tends to follow density surfaces. This inhibits vertical mixing but does not forbid horizontal rearrangement or turbulence powered by a shear flow.

A uniform current crossing topography generates stationary internal waves when its intrinsic frequency lies in the internal-wave range. For horizontal wavenumber $k$, stationarity gives $m^2=N^2/U^2-k^2$. Height-dependent wind and stratification replace this by the Scorer parameter relation derived in Question 1. A decrease of that parameter can produce a turning level, evanescence and vertical trapping, explaining lee-wave trains confined near a lower boundary. The exact channel superpositions in Question 1 also show that nonlinearity does not automatically destroy every wave solution.

Stability of the background density profile does not imply stability of every flow built on it. Sufficiently large waves can overturn density surfaces, and shear can supply energy to disturbances. Ideal advection can create fine structure, but irreversible homogenization of buoyancy in a real fluid additionally involves diffusion or dissipation. A finite-amplitude ideal solution is consequently distinct from a guaranteed stable or physically realizable state.

How rotation changes the picture. With the vertical Coriolis parameter $f$, the corresponding uniform-background inertia–gravity dispersion becomesWhen $N>|f|$, the propagating branch lies between $|f|$ and $N$. Rotation also allows a slow balanced sector: horizontal pressure and Coriolis forces balance, vertical pressure and buoyancy balance, and the circulation can be described through potential-vorticity inversion. Small Rossby number and small stratified Froude number lead to the quasi-geostrophic equations rather than to the fast internal-wave dynamics. Vertical shear is then tied to horizontal buoyancy gradients through thermal-wind balance.

In this balanced sector, a planetary or topographic PV gradient supports Rossby waves, whose restoring mechanism is conservation of PV under displacement across that gradient. Stable stratification couples horizontal levels through the ${\partial }_z[(f_0^2/N^2){\psi }_z]$ term, and the mountain-wave calculation in Question 3 illustrates its vertical propagation filter. Thus stratification can support rapid buoyancy-restored waves, constrain slowly evolving balanced circulation, and shape interactions with boundaries; which description applies depends on the velocity, length, time and rotation scales.

Once $C$ is removed, the remaining connected segments form the chain $K\to L\to N\to M$, of total length $3\ell$. For $RY\gg 1$, each hydraulic resistance draws a small flow, so to leading order the far end at $M$ is reflecting, with zero segment flow there. Successive wave transfer matrices of an arterial segment giveThe incident-plus-reflected inlet relation remains $q_K=2Y(2P_A-p_K)$. Therefore, defining $D=3e^{3i\beta }+e^{-3i\beta }$,The modulus of $D$ equals $\gamma =|3e^{4i\beta }+e^{-2i\beta }|$, because the expressions differ by a unit-modulus factor. Comparison with the normal flows gives the amplitude ratiosFor small phases the cosines are positive and the absolute values can be removed. The original PDF prints an additional $\cos (\beta /2)$ in the middle numerator. That factor does not follow from the displayed equal-length network: the transfer relation from $M$ to $K$ above fixes $p_M=8P_A/D$, with no such factor. The corrected middle ratio is shown here.

Since ${\gamma }^2=10+6\cos 6\beta$, one has $\gamma =4[1-27{\beta }^2/8+O({\beta }^4)]$. Thus the three corrected ratios, in the order $N,M,L$, areAll three exceed one for sufficiently small nonzero $\beta$; they equal one at zero frequency. This is an oscillatory arterial wave reflection effect, not an increase in steady conductance. Removing a branch changes the standing wave pattern and the total inlet pressure associated with a fixed incident wave. Reflected waves can enhance junction pressures and hence peripheral flow amplitudes. Neither total inlet pressure nor total inlet flow was held fixed, and no extra energy source is implied.

Take positive volumetric flow rate in the direction of increasing $x$. For an arterial pressure wave, the characteristic admittance of an arterial wave givesThe backward wave contributes the opposite flow sign. Write $p_j$ for each junction's complex pressure amplitude, suppressing the common $e^{i\omega t}$ factor, and put $r=RY$, $C_{\beta }=\cos \beta$, $S_{\beta }=\sin \beta$. The topology is the two paths $K\to L\to N$ and $K\to M\to N$; their intermediate and common distal junctions each have a hydraulic resistance. The coordinate convention in the source puts the upstream end at $-\ell$ and the downstream end at zero in each segment. Consequently the wave transfer matrix of an arterial segment isEquivalently $q_u=Y(C_{\beta }{p}_u-p_d)/(i{S}_{\beta })$ and $q_d=Y(p_u-C_{\beta }{p}_d)/(i{S}_{\beta })$, with their continuous limits understood when $S_{\beta }=0$.

Symmetry gives $p_L=p_M=p_s$. Distal mass conservation and the resistive relation $Q_N=p_N/R$ giveAt $L$, the incoming flow supplies both its bed and segment $D$, soChoose the incident amplitude $P_A$ at $K$. If $P_{AR}$ is the reflected amplitude in vessel $A$, pressure continuity equation gives $p_K=P_A+P_{AR}$, while inlet flow is $2Y(P_A-P_{AR})=2Y(2P_A-p_K)$. Equating this to the two outgoing segment flows yieldsEliminating $p_s,p_K$ from these three relations gives the wave transmission through a resistively loaded arterial loop resultFor large $r$ at fixed phase, $p_N=2P_A{e}^{-2i\beta }+O(r^{-1})$ and $p_M=C_{\beta }{p}_N+O(r^{-1})$. HenceThese leading terms must not be used as nonzero denominators near a phase at which the normal intermediate flow vanishes.

Steady mass conservation makes $\alpha u=Q$. For uniform area, $u=U=Q/{\alpha }_0$ is constant and the pressure-area law of a collapsible tube gives $p_x=p_e^{\prime }(x)$. The acceleration vanishes, so momentum balance requiresThe external pressure must therefore decrease linearly along a positive-flow tube. Its additive constant is free; the constant transmural pressure is $\tilde{P}({\alpha }_0)$.

The first equation is mass conservation for an incompressible flow whose area may change: local area storage balances the divergence of the volume flux $\alpha u$. The two left-hand momentum terms are local acceleration and convective acceleration, while $-p_x$ is the driving pressure gradient and $-R(\alpha )u$ models viscous drag per unit mass. The pressure-area law of a collapsible tube expresses transverse wall balance using internal minus prescribed external pressure.

For example, choose a reference area ${\mathcal{A}}_{\ast }$, length $L$ and velocity $U_{\ast }$. If dimensional variables are $\mathcal{A},s,t_d,v,P$, takeThe dimensionless external pressure uses the same scale. A dimensional friction rate $r_f(\mathcal{A})$ becomes $R(\alpha )=L{r}_f({\mathcal{A}}_{\ast }\alpha )/U_{\ast }$, and the dimensional wall pressure law is divided by $\rho U_{\ast }^2$. These choices make both inertia and the pressure-gradient coefficient unity; other equivalent reference scales are possible. On the mechanically stable wall branch, increasing area requires increased transmural pressure, whereas an enlarged tube reduces resistance. Thus $\boxed{{\tilde{P}}^{\prime }(\alpha )>0,R^{\prime }(\alpha )<0}$.

Complete the square in the dispersion relation:If the imaginary part of the square root has magnitude $b$, thenOne root grows exactly when $b>R_0/2$. For $k\ne 0$, squaring the equivalent positive inequality gives $n^2{U}^2>c_0^2$. Thus positive flow is unstable whenever $\boxed{U>c_0/n}$. This proves the entire stated range, not just a small neighborhood of threshold. Zero wavenumber has only a neutral uniform-area change and a damped velocity mode.

Near threshold write $nU=c_0(1+\delta )$ and follow the neutral branch as $\nu =k{c}_0+\delta {\nu }_1+O({\delta }^2)$. Substitution yieldsTaking its imaginary part gives the growth rateArea increases reduce friction, causing a velocity increase that feeds the area disturbance through mass conservation. The pressure-driven base flow supplies the energy; the positive friction coefficient does not by itself prevent instability of that maintained flow.

In the travelling coordinate $X=x-ct$, mass conservation becomes $[\alpha (u-c){]}^{\prime }=0$. Thus the conserved relative flux is $Q$ andFor a periodic peristaltic pumping solution, the neglected-inertia momentum equation isTake ${\alpha }_0$ to be the spatial mean area, so $⟨{\alpha }_1⟩=⟨{\alpha }_2⟩=0$. This normalization, and zero net pressure drop over a period, are implicit in the specified periodic forcing. The constant leading state has no pressure gradient, giving $Q_0=-{\alpha }_{0}c$. Put $B={\tilde{P}}_0^{\prime }$ and $d=R_{0}c/{\alpha }_0$. At first order,Period averaging gives $Q_1=0$, and the zero-mean solution isFor the second-order mean balance, expand the velocity and friction:Averaging the exact pressure derivative gives zero, so $⟨R_0{u}_2+R_1{u}_1⟩=0$. Hence the weak peristaltic pumping of a collapsible tube result isThe mean laboratory volumetric flow rate is $c⟨\alpha ⟩+Q={ϵ}^2{Q}_2+O({ϵ}^3)$, forward with the wave when $c>0$. The negative leading relative flux is therefore consistent with positive mean pumping.